Large-scale benchmarking of circRNA detection tools reveals large differences in sensitivity but not in precision.

The detection of circular RNA molecules (circRNAs) is typically based on short-read RNA sequencing data processed using computational tools. Numerous such tools have been developed, but a systematic comparison with orthogonal validation is missing. Here, we set up a circRNA detection tool benchmarking study, in which 16 tools detected more than 315,000 unique circRNAs in three deeply sequenced human cell types. Next, 1,516 predicted circRNAs were validated using three orthogonal methods. Generally, tool-specific precision is high and similar (median of 98.8%, 96.3% and 95.5% for qPCR, RNase R and amplicon sequencing, respectively) whereas the sensitivity and number of predicted circRNAs (ranging from 1,372 to 58,032) are the most significant differentiators. Of note, precision values are lower when evaluating low-abundance circRNAs. We also show that the tools can be used complementarily to increase detection sensitivity. Finally, we offer recommendations for future circRNA detection and validation.

circRIP: an accurate tool for identifying circRNA-RBP interactions.

Circular ribonucleic acids (RNAs) (circRNAs) are formed by covalently linking the downstream splice donor and the upstream splice acceptor. One of the most important functions of circRNAs is mainly exerted through binding RNA-binding proteins (RBPs). However, there is no efficient algorithm for identifying genome-wide circRNA-RBP interactions. Here, we developed a unique algorithm, circRIP, for identifying circRNA-RBP interactions from RNA immunoprecipitation sequencing (RIP-Seq) data. A simulation test demonstrated the sensitivity and specificity of circRIP. By applying circRIP, we identified 95 IGF2BP3-binding circRNAs based on the IGF2BP3 RIP-Seq dataset. We further identified 2823 and 1333 circRNAs binding to >100 RBPs in K562 and HepG2 cell lines, respectively, based on enhanced cross-linking immunoprecipitation (eCLIP) data, demonstrating the significance to survey the potential interactions between circRNAs and RBPs. In this study, we provide an accurate and sensitive tool, circRIP (https://github.com/bioinfolabwhu/circRIP), to systematically identify RBP and circRNA interactions from RIP-Seq and eCLIP data, which can significantly benefit the research community for the functional exploration of circRNAs.

CSCD2: an integrated interactional database of cancer-specific circular RNAs.

The significant function of circRNAs in cancer was recognized in recent work, so a well-organized resource is required for characterizing the interactions between circRNAs and other functional molecules (such as microRNA and RNA-binding protein) in cancer. We previously developed cancer-specific circRNA database (CSCD), a comprehensive database for cancer-specific circRNAs, which is widely used in circRNA research. Here, we updated CSCD to CSCD2 (http://geneyun.net/CSCD2 or http://gb.whu.edu.cn/CSCD2), which includes significantly more cancer-specific circRNAs identified from a large number of human cancer and normal tissues/cell lines. CSCD2 contains >1000 samples (825 tissues and 288 cell lines) and identifies a large number of circRNAs: 1 013 461 cancer-specific circRNAs, 1 533 704 circRNAs from only normal samples and 354 422 circRNAs from both cancer and normal samples. In addition, CSCD2 predicts potential miRNA-circRNA and RBP-circRNA interactions using binding motifs from >200 RBPs and 2000 microRNAs. Furthermore, the potential full-length and open reading frame sequence of these circRNAs were also predicted. Collectively, CSCD2 provides a significantly enhanced resource for exploring the function and regulation of circRNAs in cancer.

ADEIP: an integrated platform of age-dependent expression and immune profiles across human tissues.

Gene expression and immune status in human tissues are changed with aging. There is a need to develop a comprehensive platform to explore the dynamics of age-related gene expression and immune profiles across tissues in genome-wide studies. Here, we collected RNA-Seq datasets from GTEx project, containing 16 704 samples from 30 major tissues in six age groups ranging from 20 to 79 years old. Dynamic gene expression along with aging were depicted and gene set enrichment analysis was performed among those age groups. Genes from 34 known immune function categories and immune cell compositions were investigated and compared among different age groups. Finally, we integrated all the results and developed a platform named ADEIP (http://gb.whu.edu.cn/ADEIP or http://geneyun.net/ADEIP), integrating the age-dependent gene expression and immune profiles across tissues. To demonstrate the usage of ADEIP, we applied two datasets: severe acute respiratory syndrome coronavirus 2 and human mesenchymal stem cells-assoicated genes. We also included the expression and immune dynamics of these genes in the platform. Collectively, ADEIP is a powerful platform for studying age-related immune regulation in organogenesis and other infectious or genetic diseases.

Numerical and experimental investigations on the propagation property of a mid-infrared 7 × 1 multimode fiber combiner.

Mid-infrared fiber combiners have great potential in power and spectral combination. However, studies on mid-infrared transmission optical field distributions using these combiners are limited. In this study, we designed and fabricated a 7 × 1 multimode fiber combiner based on sulfur-based glass fibers and observed approximately 80% per-port transmission efficiency at 4.778 µm wavelength. We investigated the propagation properties of the prepared combiners and explored the effects of transmission wavelength, output fiber length, and fusion deviation on the transmitted optical field and beam quality factor M2. Additionally, we assessed the effect of coupling on the excitation mode and spectral combination of the mid-infrared fiber combiner for multiple light sources. Our results provide an in-depth understanding of the propagation properties of the mid-infrared multimode fiber combiners, which may find applications in high-beam-quality laser devices.

Silencing of circIgf1r plays a protective role in neuronal injury via regulating astrocyte polarization during epilepsy.

Epilepsy is a common brain disorder, repeated seizures of epilepsy may lead to a series of brain pathological changes such as neuronal or glial damage. However, whether circular RNAs are involved in neuronal injury during epilepsy is not fully understood. Here, we screened circIgf1r in the status epilepticus model through circRNA sequencing, and found that it was upregulated after the status epilepticus model through QPCR analysis. Astrocytes polarizing toward neurotoxic A1 phenotype and neurons loss were observed after status epilepticus. Through injecting circIgf1r siRNA into the lateral ventricle, it was found that knocking down circIgf1r in vivo would induce the polarization of astrocytes to phenotype A2 and reduce neuronal loss. The results in vitro further confirmed that inhibiting the expression of circIgf1r in astrocytes could protect neurons by converting reactive astrocytes from A1 to the protective A2. In addition, knocking down circIgf1r in astrocytes could functionally promote astrocyte autophagy and relieve the destruction of 4-AP-induced autophagy flux. In terms of mechanism, circIgf1r promoted the polarization of astrocytes to phenotype A1 by inhibiting autophagy. Taken together, our results reveal circIgf1r may serve as a potential target for the prevention and treatment of neuron damage after epilepsy.

CircRNA inhibits DNA damage repair by interacting with host gene.

BACKGROUND: Deregulated circular RNAs (circRNAs) are associated with the development of cancer and therapy resistance. However, functional research of circRNAs mostly focus on potential miRNA or protein binding and more potential regulation of circRNA on host gene DNA in cancers are yet to be inspected. METHOD: We performed total RNA sequencing on clinical breast cancer samples and identified the expression patterns of circRNAs and corresponding host genes in patient blood, tumor and adjacent normal tissues. qPCR, northern blot and in situ hybridization were used to validate the dysregulation of circRNA circSMARCA5. A series of procedures including R-loop dot-blotting, DNA-RNA immunoprecipitation and mass spectrum, etc. were conducted to explore the regulation of circSMARCA5 on the transcription of exon 15 of SMARCA5. Moreover, immunofluorescence and in vivo experiments were executed to investigate the overexpression of circSMARCA5 with drug sensitivities. RESULTS: We found that circRNAs has average higher expression over its host linear genes in peripheral blood. Compared to adjacent normal tissues, circSMARCA5 is decreased in breast cancer tissues, contrary to host gene SMARCA5. The enforced expression of circSMARCA5 induced drug sensitivity of breast cancer cell lines in vitro and in vivo. Furthermore, we demonstrated that circSMARCA5 can bind to its parent gene locus, forming an R-loop, which results in transcriptional pausing at exon 15 of SMARCA5. CircSMARCA5 expression resulted in the downregulation of SMARCA5 and the production of a truncated nonfunctional protein, and the overexpression of circSMARCA5 was sufficient to improve sensitivity to cytotoxic drugs. CONCLUSION: Our results revealed a new regulatory mechanism for circRNA on its host gene and provided evidence that circSMARCA5 may serve as a therapeutic target for drug-resistant breast cancer patients.

CircView: a visualization and exploration tool for circular RNAs.

Circular RNAs (circRNAs) are novel rising stars of noncoding RNAs, which are highly abundant and evolutionarily conserved across species. Number of publications related to circRNAs increased sharply in recent years, representing emerging focuses in the field. Therefore, tools, pipelines and databases have been developed to identify and store circRNAs. However, there is no existing tool to visualize and explore circRNAs. Therefore, we introduce CircView, a user-friendly visualization tool for circRNAs detected from existing tools. CircView enables users to visualize circRNAs and to quantify number of samples with detected circRNAs. CircView allows users to explore circRNAs detected by unique or multiple tools. Furthermore, CircView allows users to view the regulatory elements, such as microRNA response elements and RNA-binding protein binding sites. CircView is a unique tool to visualize and explore circRNAs, which helps users to better understand potential functions of circRNAs and design the functional experiments.

Genome-wide characterization of lncRNAs in acute myeloid leukemia.

Long noncoding RNAs (lncRNAs) are a large family of noncoding RNAs that play a critical role in various normal bioprocesses as well as tumorigenesis. However, the expression patterns and biological functions of lncRNAs in acute leukemia have not been well studied. Here, we performed transcriptome-wide lncRNA expression profiling of acute myeloid leukemia (AML) patient samples, along with non-leukemia control hematopoietic samples. We found that lncRNAs were differentially expressed in AML samples relative to control samples. Notably, we identified that lncRNAs upregulated in AML (relative to the control samples) are associated with a lower degree of DNA methylation and a higher ratio of being bound by transcription factors such as SP1, STAT4, ATF-2 and ELK-1 compared with those downregulated in AML. Moreover, an enrichment of H3K4me3 and a depletion of H3K27me3 were observed in upregulated lncRNAs in AML. Expression patterns of three types of lncRNAs (antisense, enhancer and intergenic lncRNAs) have previously been characterized. Of the identified lncRNAs, we found that high expression level lncRNA LOC285758 is associated with the poor prognosis in AML patients. Furthermore, we found that LOC285758 regulates proliferation of AML cell lines by enhancing the expression of HDAC2, a key factor in carcinogenesis. Collectively, our study depicts a landscape of important lncRNAs in AML and provides novel potential therapeutic targets and prognostic markers for AML treatment.

CSCD: a database for cancer-specific circular RNAs.

Circular RNA (circRNA) is a large group of RNA family extensively existed in cells and tissues. High-throughput sequencing provides a way to view circRNAs across different samples, especially in various diseases. However, there is still no comprehensive database for exploring the cancer-specific circRNAs. We collected 228 total RNA or polyA(-) RNA-seq samples from both cancer and normal cell lines, and identified 272 152 cancer-specific circRNAs. A total of 950 962 circRNAs were identified in normal samples only, and 170 909 circRNAs were identified in both tumor and normal samples, which could be further used as non-tumor background. We constructed a cancer-specific circRNA database (CSCD, http://gb.whu.edu.cn/CSCD). To understand the functional effects of circRNAs, we predicted the microRNA response element sites and RNA binding protein sites for each circRNA. We further predicted potential open reading frames to highlight translatable circRNAs. To understand the association between the linear splicing and the back-splicing, we also predicted the splicing events in linear transcripts of each circRNA. As the first comprehensive cancer-specific circRNA database, we believe CSCD could significantly contribute to the research for the function and regulation of cancer-associated circRNAs.

Loss of miR-146b-5p promotes T cell acute lymphoblastic leukemia migration and invasion via the IL-17A pathway.

Metastatic disease remains the primary cause of death for individuals with T cell acute lymphoblastic leukemia (T-ALL). microRNAs (miRNAs) play important roles in the pathogenesis of T-ALL by inhibiting gene expression at posttranscriptional levels. The goal of the current project is to identify any significant miRNAs in T-ALL metastasis. We observed miR-146b-5p to be downregulated in T-ALL patients and cell lines, and bioinformatics analysis implicated miR-146b-5p in the hematopoietic system. miR-146b-5p inhibited the migration and invasion in T-ALL cells. Interleukin-17A (IL-17A) was predicted to be a target of miR-146b-5p; this was confirmed by luciferase assays. Interestingly, T-ALL patients and cell lines secreted IL-17A and expressed the IL-17A receptor (IL-17RA). IL-17A/IL-17RA interactions promoted strong T-ALL cell migration and invasion responses. Gene set enrichment analysis (GSEA) and quantitative polymerase chain reaction (qPCR) analysis indicated that matrix metallopeptidase-9 (MMP9), was a potential downstream effector of IL-17A activation, and nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) signaling was also implicated in this process. Moreover, IL-17A activation promoted T-ALL cell metastasis to the liver in IL17A -/- mouse models. These results indicate that reduced miR-146b-5p expression in T-ALL may lead to the upregulation of IL-17A, which then promotes T-ALL cell migration and invasion by upregulating MMP9 via NF-κB signaling.

The prognostic value of matrix metalloproteinase-7 and matrix metalloproteinase-15 in acute myeloid leukemia.

Matrix metalloproteinases (MMPs), a family of zinc-dependent endopeptidases, are involved in a variety of physiological and pathological processes. We analyzed 11 data sets from Gene Expression Omnibus Database and found that MMP7 and MMP15 were highly expressed in multiple carcinomas. GSE13204 showed that MMP7 and MMP15 were overexpressed in acute myeloid leukemia (AML) patients. The Cancer Genome Atlas data set exhibited that high expression of MMP7 or MMP15 in bone marrow (BM) of AML patients predicted poor overall survival. The χ 2 test results indicated that high expression level of MMP7 and MMP15 were correlated with high-risk stratification and high BM blast cell percentage in AML patients. To confirm these findings, we performed reverse-transcription quantitative polymerase chain reaction (RT-qPCR) and found that MMP7 and MMP15 were highly expressed in three AML cell lines. Further study showed that MMP7 and MMP15 were highly expressed both in BM and peripheral blood in collected AML samples compared with healthy individuals. Additionally, long noncoding RNA (lncRNA) microarray of BM samples of AML patients revealed that multiple lncRNAs were correlated with MMP7 and MMP15, suggesting that lncRNAs might be involved in the pathogenesis of AML via modulating MMPs. In conclusion, our study uncovers the potential roles of MMP7 and MMP15 in the prognosis of AML.

Genome-wide identification of cancer-specific alternative splicing in circRNA.

Circular RNA (circRNA) is a group of RNA families generated by RNA circularization, which was discovered ubiquitously across different cancers. However, the internal structure of circRNA is difficult to determine due to alternative splicing that occurs in its exons and introns. Furthermore, cancer-specific alternative splicing of circRNA is less likely to be identified. Here, we proposed a de novo algorithm, CircSplice, that could identify internal alternative splicing in circRNA and compare differential circRNA splicing events between different conditions ( http://gb.whu.edu.cn/CircSplice or https://github.com/GeneFeng/CircSplice ). By applying CircSplice in clear cell renal cell carcinoma and bladder cancer, we detected 4498 and 2977 circRNA alternative splicing (circ-AS) events in the two datasets respectively and confirmed the expression of circ-AS events by RT-PCR. We further inspected the distributions and patterns of circ-AS in cancer and adjacent normal tissues. To further understand the potential functions of cancer-specific circ-AS, we classified those events into tumor suppressors and oncogenes and performed pathway enrichment analysis. This study is the first comprehensive view of cancer-specific circRNA alternative splicing, which could contribute significantly to regulation and functional research of circRNAs in cancers.

Comprehensive characterization of tissue-specific circular RNAs in the human and mouse genomes.

Circular RNA (circRNA) is a group of RNA family generated by RNA circularization, which was discovered ubiquitously across different species and tissues. However, there is no global view of tissue specificity for circRNAs to date. Here we performed the comprehensive analysis to characterize the features of human and mouse tissue-specific (TS) circRNAs. We identified in total 302 853 TS circRNAs in the human and mouse genome, and showed that the brain has the highest abundance of TS circRNAs. We further confirmed the existence of circRNAs by reverse transcription polymerase chain reaction (RT-PCR). We also characterized the genomic location and conservation of these TS circRNAs and showed that the majority of TS circRNAs are generated from exonic regions. To further understand the potential functions of TS circRNAs, we identified microRNAs and RNA binding protein, which might bind to TS circRNAs. This process suggested their involvement in development and organ differentiation. Finally, we constructed an integrated database TSCD (Tissue-Specific CircRNA Database: http://gb.whu.edu.cn/TSCD) to deposit the features of TS circRNAs. This study is the first comprehensive view of TS circRNAs in human and mouse, which shed light on circRNA functions in organ development and disorders.

IGF1/MAPK/ERK signaling pathway-mediated programming alterations of adrenal cortex cell proliferation by prenatal caffeine exposure in male offspring rats.

Our previous study proposed a glucocorticoid-insulin-like growth factor 1 (GC-IGF1) axis programming mechanism for prenatal caffeine exposure (PCE)-induced adrenal developmental dysfunction. Here, we focused on PCE-induced cell proliferation changes of the adrenal cortex in male offspring rats before and after birth and clarified the intrauterine programming mechanism. On gestational day (GD) 20, the PCE group had an elevated serum corticosterone level reduced fetal bodyweight, maximum adrenal sectional area, and elevated adrenal corticosterone and aldosterone contents. However, in postnatal week (PW) 6, the serum corticosterone level was decreased, and the bodyweight, with catch-up growth, adrenal cortex maximum cross-sectional area and aldosterone content were relatively increased, while the adrenal corticosterone content was lower. On GD20, the expression of adrenal IGF1, IGF1R and proliferating cell nuclear antigen (PCNA) were decreased, while the expression of these factors at PW6 were increased in the PCE group. Fetal adrenal gene chip analysis suggested that the mitogen-activated protein kinase/extracellular regulated protein kinase (MAPK/ERK) signal pathway was suppressed in the PCE group. Moreover, in the rat primary adrenal cells, corticosterone (rather than caffeine) was shown to significantly inhibit cell proliferation, IGF1 and PCNA expression, and ERK phosphorylation, which could be reversed by exogenous IGF1. Meanwhile, the effects of exogenous IGF1 were reversed by the ERK pathway inhibitor (PD184161). In conclusion, PCE could induce programming alterations in adrenal cortical cell proliferation before and after birth in male offspring rats. The underlying mechanism is associated with the inhibition of fetal adrenal IGF1-related MAPK/ERK signaling pathway caused by high glucocorticoid levels.

TET1 plays an essential oncogenic role in MLL-rearranged leukemia.

The ten-eleven translocation 1 (TET1) gene is the founding member of the TET family of enzymes (TET1/2/3) that convert 5-methylcytosine to 5-hydroxymethylcytosine. Although TET1 was first identified as a fusion partner of the mixed lineage leukemia (MLL) gene in acute myeloid leukemia carrying t(10,11), its definitive role in leukemia is unclear. In contrast to the frequent down-regulation (or loss-of-function mutations) and critical tumor-suppressor roles of the three TET genes observed in various types of cancers, here we show that TET1 is a direct target of MLL-fusion proteins and is significantly up-regulated in MLL-rearranged leukemia, leading to a global increase of 5-hydroxymethylcytosine level. Furthermore, our both in vitro and in vivo functional studies demonstrate that Tet1 plays an indispensable oncogenic role in the development of MLL-rearranged leukemia, through coordination with MLL-fusion proteins in regulating their critical cotargets, including homeobox A9 (Hoxa9)/myeloid ecotropic viral integration 1 (Meis1)/pre-B-cell leukemia homeobox 3 (Pbx3) genes. Collectively, our data delineate an MLL-fusion/Tet1/Hoxa9/Meis1/Pbx3 signaling axis in MLL-rearranged leukemia and highlight TET1 as a potential therapeutic target in treating this presently therapy-resistant disease.

miR-9 is an essential oncogenic microRNA specifically overexpressed in mixed lineage leukemia-rearranged leukemia.

MicroRNAs (miRNAs), small noncoding RNAs that regulate target gene mRNAs, are known to contribute to pathogenesis of cancers. Acute myeloid leukemia (AML) is a group of heterogeneous hematopoietic malignancies with various chromosomal and/or molecular abnormalities. AML with chromosomal translocations involving the mixed lineage leukemia (MLL) gene are usually associated with poor survival. In the present study, through a large-scale, genomewide miRNA expression assay, we show that microRNA-9 (miR-9) is the most specifically up-regulated miRNA in MLL-rearranged AML compared with both normal control and non-MLL-rearranged AML. We demonstrate that miR-9 is a direct target of MLL fusion proteins and can be significantly up-regulated in expression by the latter in human and mouse hematopoietic stem/progenitor cells. Depletion of endogenous miR-9 expression by an appropriate antagomiR can significantly inhibit cell growth/viability and promote apoptosis in human MLL-rearranged AML cells, and the opposite is true when expression of miR-9 is forced. Blocking endogenous miR-9 function by anti-miRNA sponge can significantly inhibit, whereas forced expression of miR-9 can significantly promote, MLL fusion-induced immortalization/transformation of normal mouse bone marrow progenitor cells in vitro. Furthermore, forced expression of miR-9 can significantly promote MLL fusion-mediated leukemogenesis in vivo. In addition, a group of putative target genes of miR-9 exhibited a significant inverse correlation of expression with miR-9 in a series of leukemia sample sets, suggesting that they are potential targets of miR-9 in MLL-rearranged AML. Collectively, our data demonstrate that miR-9 is a critical oncomiR in MLL-rearranged AML and can serve as a potential therapeutic target to treat this dismal disease.

Pravastatin reverses obesity-induced dysfunction of induced pluripotent stem cell-derived endothelial cells via a nitric oxide-dependent mechanism.

AIMS: High-fat diet-induced obesity (DIO) is a major contributor to type II diabetes and micro- and macro-vascular complications leading to peripheral vascular disease (PVD). Metabolic abnormalities of induced pluripotent stem cell-derived endothelial cells (iPSC-ECs) from obese individuals could potentially limit their therapeutic efficacy for PVD. The aim of this study was to compare the function of iPSC-ECs from normal and DIO mice using comprehensive in vitro and in vivo assays. METHODS AND RESULTS: Six-week-old C57Bl/6 mice were fed with a normal or high-fat diet. At 24 weeks, iPSCs were generated from tail tip fibroblasts and differentiated into iPSC-ECs using a directed monolayer approach. In vitro functional analysis revealed that iPSC-ECs from DIO mice had significantly decreased capacity to form capillary-like networks, diminished migration, and lower proliferation. Microarray and ELISA confirmed elevated apoptotic, inflammatory, and oxidative stress pathways in DIO iPSC-ECs. Following hindlimb ischaemia, mice receiving intramuscular injections of DIO iPSC-ECs had significantly decreased reperfusion compared with mice injected with control healthy iPSC-ECs. Hindlimb sections revealed increased muscle atrophy and presence of inflammatory cells in mice receiving DIO iPSC-ECs. When pravastatin was co-administered to mice receiving DIO iPSC-ECs, a significant increase in reperfusion was observed; however, this beneficial effect was blunted by co-administration of the nitric oxide synthase inhibitor, N(ω)-nitro-l-arginine methyl ester. CONCLUSION: This is the first study to provide evidence that iPSC-ECs from DIO mice exhibit signs of endothelial dysfunction and have suboptimal efficacy following transplantation in a hindlimb ischaemia model. These findings may have important implications for future treatment of PVD using iPSC-ECs in the obese population.

PBX3 is an important cofactor of HOXA9 in leukemogenesis.

Although PBX proteins are known to increase DNA-binding/transcriptional activity of HOX proteins through their direct binding, the functional importance of their interaction in leukemogenesis is unclear.We recently reported that overexpression of a 4-homeobox-gene signature (ie, PBX3/HOXA7/HOXA9/HOXA11) is an independent predictor of poor survival in patients with cytogenetically abnormal acute myeloid leukemia (CA-AML). Here we show that it is PBX3, but not PBX1 or PBX2, that is consistently coexpressed with HOXA9 in various subtypes of CA-AML, particularly MLL-rearranged AML, and thus appears as a potential pathologic cofactor of HOXA9 in CA-AML. We then show that depletion of endogenous Pbx3 expression by shRNA significantly inhibits MLL-fusion-mediated cell transformation, and coexpressed PBX3 exhibits a significantly synergistic effect with HOXA9 in promoting cell transformation in vitro and leukemogenesis in vivo. Furthermore, as a proof of concept, we show that a small peptide, namely HXR9, which was developed to specifically disrupt the interactions between HOX and PBX proteins, can selectively kill leukemic cells with overexpression of HOXA/PBX3 genes. Collectively, our data suggest that PBX3 is a critical cofactor of HOXA9 in leukemogenesis, and targeting their interaction is a feasible strategy to treat presently therapy resistant CA-AML (eg, MLL-rearranged leukemia) in which HOXA/PBX3 genes are overexpressed.

MiR-495 is a tumor-suppressor microRNA down-regulated in MLL-rearranged leukemia.

Acute myeloid leukemia (AML) is a heterogeneous group of hematopoietic malignancies with variable response to treatment. AMLs bearing MLL (mixed lineage leukemia) rearrangements are associated with intermediate or poor survival. MicroRNAs (miRNAs), a class of small noncoding RNAs, have been postulated to be important gene expression regulators virtually in all biological processes, including leukemogenesis. Through a large-scale, genome-wide miRNA expression profiling assay of 85 human AML and 15 normal control samples, we show that among 48 miRNAs that are significantly differentially expressed between MLL- and non-MLL-rearranged AML samples, only one (miR-495) is expressed at a lower level in MLL-rearranged AML than in non-MLL-rearranged AML; meanwhile, miR-495 is also significantly down-regulated in MLL-rearranged AML samples compared with normal control samples. Through in vitro colony-forming/replating assays and in vivo bone marrow transplantation studies, we show that forced expression of miR-495 significantly inhibits MLL-fusion-mediated cell transformation in vitro and leukemogenesis in vivo. In human leukemic cells carrying MLL rearrangements, ectopic expression of miR-495 greatly inhibits cell viability and increases cell apoptosis. Furthermore, our studies demonstrate that PBX3 and MEIS1 are two direct target genes of miR-495, and forced expression of either of them can reverse the effects of miR-495 overexpression on inhibiting cell viability and promoting apoptosis of human MLL-rearranged leukemic cells. Thus, our data indicate that miR-495 likely functions as a tumor suppressor in AML with MLL rearrangements by targeting essential leukemia-related genes.