Digital gene expression profiling of primary acute lymphoblastic leukemia cells

Integrated systems biology analysis of acute lymphoblastic leukemia: unveiling molecular signatures and drug repurposing opportunities

Acute lymphoblastic leukemia (ALL) is a hematological malignancy characterized by aberrant proliferation and accumulation of lymphoid precursor cells within the bone marrow. The tyrosine kinase inhibitor (TKI), imatinib mesylate, has played a significant role in the treatment of Philadelphia chromosome-positive ALL (Ph + ALL). However, the achievement of durable and sustained therapeutic success remains a challenge due to the development of TKI resistance during the clinical course.The primary objective of this investigation is to propose a novel and efficacious treatment approach through drug repositioning, targeting ALL and its Ph + subtype by identifying and addressing differentially expressed genes (DEGs). This study involves a comprehensive analysis of transcriptome datasets pertaining to ALL and Ph + ALL in order to identify DEGs associated with the progression of these diseases to identify possible repurposable drugs that target identified hub proteins.The outcomes of this research have unveiled 698 disease-related DEGs for ALL and 100 for Ph + ALL. Furthermore, a subset of drugs, specifically glipizide for Ph + ALL, and maytansine and isoprenaline for ALL, have been identified as potential candidates for therapeutic intervention. Subsequently, cytotoxicity assessments were performed to confirm the in vitro cytotoxic effects of these selected drugs on both ALL and Ph + ALL cell lines.In conclusion, this study offers a promising avenue for the management of ALL and Ph + ALL through drug repurposed drugs. Further investigations are necessary to elucidate the mechanisms underlying cell death, and clinical trials are recommended to validate the promising results obtained through drug repositioning strategies.

BACH2-mediated CD28 and CD40LG axes contribute to pathogenesis and progression of T-cell lymphoblastic leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive subtype of ALL characterized by its high heterogeneity and unfavorable clinical features. Despite improved insights in genetic and epigenetic landscapes of T-ALL, the molecular mechanisms that drive malignant T-cell development remain unclear. BTB and CNC homology 2 (BACH2) is a lymphoid-specific transcription repressor recognized as a tumor suppressor in B-cell malignancies, but little is known about its function and regulatory network in T-ALL. Here we found extremely low levels of BACH2 in T-ALL clinical samples and cell lines compared to normal T cells. Overexpression of BACH2 in T-ALL cells not only induced cell growth retardation but also inhibited cancer progression and infiltration in xenografts. Further RNA sequencing (RNA-seq) analysis revealed significant alterations in regulation of defense and immune responses in T-ALL cells upon BACH2 overexpression. Strikingly, CD28 and CD40LG, two essential stimulatory molecules on T cells, were for the first time identified as novel downstream targets repressed by BACH2 in T-ALL cells. Interestingly, both CD28 and CD40LG were indispensable for T-ALL survival, since largely or completely silencing CD28 and CD40LG led to rapid cell death, whereas partial knockdown of them resulted in cell-cycle arrest and enhanced apoptosis. More importantly, BACH2-mediated CD28 and CD40LG signals contributed to cell migration and dissemination of T-ALL cells to the bone marrow, thus adding a new layer to the BACH2-mediated tumor immunoregulation in T-cell malignancies.

Aberrant LncRNA Expression in Leukemia

Leukemia is a common malignant cancer of the hematopoietic system, whose pathogenesis has not been fully elucidated. Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nucleotides without protein-coding function. Recent studies report their role in cellular processes such as the regulation of gene expression, as well as in the carcinogenesis, occurrence, development, and prognosis of various tumors. Evidence indicating relationships between a variety of lncRNAs and leukemia pathophysiology has increased dramatically in the previous decade, with specific lncRNAs expected to serve as diagnostic biomarkers, novel therapeutic targets, and predictors of clinical outcomes. Furthermore, these lncRNAs might offer insight into disease pathogenesis and novel treatment options. This review summarizes progress in studies on the role(s) of lncRNAs in leukemia.

Clinicopathological features and prognostic value of SOX11 in childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia is marked by aberrant transcriptional features that alter cell differentiation, self-renewal, and proliferative features. We sought to identify the transcription factors exhibiting altered and subtype-specific expression patterns in B-ALL and report here that SOX11, a developmental and neuronal transcription factor, is aberrantly expressed in the ETV6-RUNX1 and TCF3-PBX1 subtypes of acute B-cell leukemias. We show that a high expression of SOX11 leads to alterations of gene expression that are typically associated with cell adhesion, migration, and differentiation. A high expression is associated with DNA hypomethylation at the SOX11 locus and a favorable outcome. The results indicate that SOX11 expression marks a group of patients with good outcomes and thereby prompts further study of its use as a biomarker.

Long Non-Coding RNA and Acute Leukemia

Acute leukemia (AL) is the main type of cancer in children worldwide. Mortality by this disease is high in developing countries and its etiology remains unanswered. Evidences showing the role of the long non-coding RNAs (lncRNAs) in the pathophysiology of hematological malignancies have increased drastically in the last decade. In addition to the contribution of these lncRNAs in leukemogenesis, recent studies have suggested that lncRNAs could be used as biomarkers in the diagnosis, prognosis, and therapeutic response in leukemia patients. The focus of this review is to describe the functional classification, biogenesis, and the role of lncRNAs in leukemogenesis, to summarize the evidence about the lncRNAs which are playing a role in AL, and how these genes could be useful as potential therapeutic targets.

Differential genome-wide profiling of alternative polyadenylation sites in nasopharyngeal carcinoma by high-throughput sequencing

Background

Alternative polyadenylation (APA) is a widespread phenomenon in the posttranscriptional regulation of gene expression that generates mRNAs with alternative 3′-untranslated regions (3’UTRs). APA contributes to the pathogenesis of various diseases, including cancer. However, the potential role of APA in the development of nasopharyngeal carcinoma (NPC) remains largely unknown.

Methods

A strategy of sequencing APA sites (SAPAS) based on second-generation sequencing technology was carried out to explore the global patterns of APA sites and identify genes with tandem 3’UTRs in samples from 6 NPC and 6 normal nasopharyngeal epithelial tissue (NNET). Sequencing results were then validated using quantitative RT-PCR in a larger cohort of 16 NPC and 16 NNET samples.

Results

The sequencing data showed that the use of tandem APA sites was prevalent in NPC, and numerous genes with APA-switching events were discovered. In total, we identified 195 genes with significant differences in the tandem 3’UTR length between NPC and NNET; including 119 genes switching to distal poly (A) sites and 76 genes switching to proximal poly (A) sites. Several gene ontology (GO) terms were enriched in the list of genes with switched APA sites, including regulation of cell migration, macromolecule catabolic process, protein catabolic process, proteolysis, small conjugating protein ligase activity, and ubiquitin-protein ligase activity.

Conclusions

APA site-switching events are prevalent in NPC. APA-mediated regulation of gene expression may play an important role in the development of NPC, and more detailed studies targeting genes with APA-switching events may contribute to the development of novel future therapeutic strategies for NPC.

Electronic supplementary material

The online version of this article (10.1186/s12929-018-0477-6) contains supplementary material, which is available to authorized users.

Epigenetics in pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common malignancy in children. ALL arises from the malignant transformation of progenitor B- and T-cells in the bone marrow into leukemic cells, but the mechanisms underlying this transformation are not well understood. Recent technical advances and decreasing costs of methods for high-throughput DNA sequencing and SNP genotyping have stimulated systematic studies of the epigenetic changes in leukemic cells from pediatric ALL patients. The results emerging from these studies are increasing our understanding of the epigenetic component of leukemogenesis and have demonstrated the potential of DNA methylation as a biomarker for lineage and subtype classification, prognostication, and disease progression in ALL. In this review, we provide a concise examination of the epigenetic studies in ALL, with a focus on DNA methylation and mutations perturbing genes involved in chromatin modification, and discuss the future role of epigenetic analyses in research and clinical management of ALL.

Deregulation of the non-coding genome in leukemia

Methodological advances that allow deeper characterization of non-coding elements in the genome have started to reveal the full spectrum of deregulation in cancer. We generated an inducible cell model to track transcriptional changes after induction of a well-known leukemia-inducing fusion gene, ETV6-RUNX1. Our data revealed widespread transcriptional alterations outside coding elements in the genome. This adds to the growing list of various alterations in the non-coding genome in cancer and pinpoints their role in diseased cellular state.

A L1-regularized feature selection method for local dimension reduction on microarray data

Dimension reduction is a crucial technique in machine learning and data mining, which is widely used in areas of medicine, bioinformatics and genetics. In this paper, we propose a two-stage local dimension reduction approach for classification on microarray data. In first stage, a new L1-regularized feature selection method is defined to remove irrelevant and redundant features and to select the important features (biomarkers). In the next stage, PLS-based feature extraction is implemented on the selected features to extract synthesis features that best reflect discriminating characteristics for classification. The suitability of the proposal is demonstrated in an empirical study done with ten widely used microarray datasets, and the results show its effectiveness and competitiveness compared with four state-of-the-art methods. The experimental results on St Jude dataset shows that our method can be effectively applied to microarray data analysis for subtype prediction and the discovery of gene coexpression.

Sense-antisense gene-pairs in breast cancer and associated pathological pathways

More than 30% of human protein-coding genes form hereditary complex genome architectures composed of sense-antisense (SA) gene pairs (SAGPs) transcribing their RNAs from both strands of a given locus. Such architectures represent important novel components of genome complexity contributing to gene expression deregulation in cancer cells. Therefore, the architectures might be involved in cancer pathways and, in turn, be used for novel drug targets discovery. However, the global roles of SAGPs in cancer pathways has not been studied. Here we investigated SAGPs associated with breast cancer (BC)-related pathways using systems biology, prognostic survival and experimental methods. Gene expression analysis identified 73 BC-relevant SAGPs that are highly correlated in BC. Survival modelling and metadata analysis of the 1161 BC patients allowed us to develop a novel patient prognostic grouping method selecting the 12 survival-significant SAGPs. The qRT-PCR-validated 12-SAGP prognostic signature reproducibly stratified BC patients into low- and high-risk prognostic subgroups. The 1381 SAGP-defined differentially expressed genes common across three studied cohorts were identified. The functional enrichment analysis of these genes revealed the GABPA gene network, including BC-relevant SAGPs, specific gene sets involved in cell cycle, spliceosomal and proteasomal pathways. The co-regulatory function of GABPA in BC cells was supported using siRNA knockdown studies. Thus, we demonstrated SAGPs as the synergistically functional genome architectures interconnected with cancer-related pathways and associated with BC patient clinical outcomes. Taken together, SAGPs represent an important component of genome complexity which can be used to identify novel aspects of coordinated pathological gene networks in cancers.

Digital gene expression profiling analysis of childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most commonly diagnosed malignancy in children. It is a heterogeneous disease, and is determined by multiple gene alterations and chromosomal rearrangements. To improve current understanding of the underlying molecular mechanisms of ALL, the present study profiled genome‑wide digital gene expression (DGE) in a population of children with ALL in China. Using second‑generation sequencing technology, the profiling revealed that 2,825 genes were upregulated and 1,952 were downregulated in the ALL group. Based on the DGE profiling data, the present study further investigated seven genes (WT1, RPS26, MSX1, CD70, HOXC4, HOXA5 and HOXC6) using reverse transcription‑quantitative polymerase chain reaction analysis. Gene Ontology analysis suggested that the differentially expressed genes were predominantly involved in immune cell differentiation, metabolic processes and programmed cell death. The results of the present study provided novel insights into the gene expression patterns in children with ALL.

Transcriptome profiling and pathway analysis of hepatotoxicity induced by tris (2-ethylhexyl) trimellitate (TOTM) in mice

Tris (2-ethylhexyl) trimellitate (TOTM) is commonly used as an alternative plasticizer for medical devices. But very little information was available on its biological effects. In this study, we investigated toxicity effects of TOTM on hepatic differential gene expression analyzed by using high-throughput sequencing analysis for over-represented functions and phenotypically anchored to complementary histopathologic, and biochemical data in the liver of mice. Among 1668 candidate genes, 694 genes were up-regulated and 974 genes were down-regulated after TOTM exposure. Using Gene Ontology analysis, TOTM affected three processes: the cell cycle, metabolic process and oxidative activity. Furthermore, 11 key genes involved in the above processes were validated by real time PCR. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that these genes were involved in the cell cycle pathway, lipid metabolism and oxidative process. It revealed the transcriptome gene expression response to TOTM exposure in mouse, and these data could contribute to provide a clearer understanding of the molecular mechanisms of TOTM-induced hepatotoxicity in human.

DNA methylation-based subtype prediction for pediatric acute lymphoblastic leukemia

BACKGROUND

We present a method that utilizes DNA methylation profiling for prediction of the cytogenetic subtypes of acute lymphoblastic leukemia (ALL) cells from pediatric ALL patients. The primary aim of our study was to improve risk stratification of ALL patients into treatment groups using DNA methylation as a complement to current diagnostic methods. A secondary aim was to gain insight into the functional role of DNA methylation in ALL.

RESULTS

We used the methylation status of ~450,000 CpG sites in 546 well-characterized patients with T-ALL or seven recurrent B-cell precursor ALL subtypes to design and validate sensitive and accurate DNA methylation classifiers. After repeated cross-validation, a final classifier was derived that consisted of only 246 CpG sites. The mean sensitivity and specificity of the classifier across the known subtypes was 0.90 and 0.99, respectively. We then used DNA methylation classification to screen for subtype membership of 210 patients with undefined karyotype (normal or no result) or non-recurrent cytogenetic aberrations ('other' subtype). Nearly half (n = 106) of the patients lacking cytogenetic subgrouping displayed highly similar methylation profiles as the patients in the known recurrent groups. We verified the subtype of 20% of the newly classified patients by examination of diagnostic karyotypes, array-based copy number analysis, and detection of fusion genes by quantitative polymerase chain reaction (PCR) and RNA-sequencing (RNA-seq). Using RNA-seq data from ALL patients where cytogenetic subtype and DNA methylation classification did not agree, we discovered several novel fusion genes involving ETV6, RUNX1, and PAX5.

CONCLUSIONS

Our findings indicate that DNA methylation profiling contributes to the clarification of the heterogeneity in cytogenetically undefined ALL patient groups and could be implemented as a complementary method for diagnosis of ALL. The results of our study provide clues to the origin and development of leukemic transformation. The methylation status of the CpG sites constituting the classifiers also highlight relevant biological characteristics in otherwise unclassified ALL patients.

Cell and molecular determinants of in vivo efficacy of the BH3 mimetic ABT-263 against pediatric acute lymphoblastic leukemia xenografts

PURPOSE

Predictive biomarkers are required to identify patients who may benefit from the use of BH3 mimetics such as ABT-263. This study investigated the efficacy of ABT-263 against a panel of patient-derived pediatric acute lymphoblastic leukemia (ALL) xenografts and utilized cell and molecular approaches to identify biomarkers that predict in vivo ABT-263 sensitivity.

EXPERIMENTAL DESIGN

The in vivo efficacy of ABT-263 was tested against a panel of 31 patient-derived ALL xenografts composed of MLL-, BCP-, and T-ALL subtypes. Basal gene expression profiles of ALL xenografts were analyzed and confirmed by quantitative RT-PCR, protein expression and BH3 profiling. An in vitro coculture assay with immortalized human mesenchymal cells was utilized to build a predictive model of in vivo ABT-263 sensitivity.

RESULTS

ABT-263 demonstrated impressive activity against pediatric ALL xenografts, with 19 of 31 achieving objective responses. Among BCL2 family members, in vivo ABT-263 sensitivity correlated best with low MCL1 mRNA expression levels. BH3 profiling revealed that resistance to ABT-263 correlated with mitochondrial priming by NOXA peptide, suggesting a functional role for MCL1 protein. Using an in vitro coculture assay, a predictive model of in vivo ABT-263 sensitivity was built. Testing this model against 11 xenografts predicted in vivo ABT-263 responses with high sensitivity (50%) and specificity (100%).

CONCLUSION

These results highlight the in vivo efficacy of ABT-263 against a broad range of pediatric ALL subtypes and shows that a combination of in vitro functional assays can be used to predict its in vivo efficacy.

Molecular mechanisms of reduced nerve toxicity by titanium dioxide nanoparticles in the phoxim-exposed brain of Bombyx mori

Bombyx mori (B. mori), silkworm, is one of the most important economic insects in the world, while phoxim, an organophosphorus (OP) pesticide, impact its economic benefits seriously. Phoxim exposure can damage the brain, fatbody, midgut and haemolymph of B. mori. However the metabolism of proteins and carbohydrates in phoxim-exposed B. mori can be improved by Titanium dioxide nanoparticles (TiO2 NPs). In this study, we explored whether TiO2 NPs treatment can reduce the phoxim-induced brain damage of the 5th larval instar of B. mori. We observed that TiO2 NPs pretreatments significantly reduced the mortality of phoxim-exposed larva and relieved severe brain damage and oxidative stress under phoxim exposure in the brain. The treatments also relieved the phoxim-induced increases in the contents of acetylcholine (Ach), glutamate (Glu) and nitric oxide (NO) and the phoxim-induced decreases in the contents of norepinephrine (NE), Dopamine (DA), and 5-hydroxytryptamine (5-HT), and reduced the inhibition of acetylcholinesterase (AChE), Na+/K+-ATPase, Ca2+-ATPase, and Ca2+/Mg2+-ATPase activities and the activation of total nitric oxide synthase (TNOS) in the brain. Furthermore, digital gene expression profile (DGE) analysis and real time quantitative PCR (qRT-PCR) assay revealed that TiO2 NPs pretreatment inhibited the up-regulated expression of ace1, cytochrome c, caspase-9, caspase-3, Bm109 and down-regulated expression of BmIap caused by phoxim; these genes are involved in nerve conduction, oxidative stress and apoptosis. TiO2 NPs pretreatment also inhibited the down-regulated expression of H+ transporting ATP synthase and vacuolar ATP synthase under phoxim exposure, which are involved in ion transport and energy metabolism. These results indicate that TiO2 NPs pretreatment reduced the phoxim-induced nerve toxicity in the brain of B. mori.

RNA sequencing: from tag-based profiling to resolving complete transcript structure

Technological advances in the sequencing field support in-depth characterization of the transcriptome. Here, we review genome-wide RNA sequencing methods used to investigate specific aspects of gene expression and its regulation, from transcription to RNA processing and translation. We discuss tag-based methods for studying transcription, alternative initiation and polyadenylation events, shotgun methods for detection of alternative splicing, full-length RNA sequencing for the determination of complete transcript structures, and targeted methods for studying the process of transcription and translation. With the ensemble of technologies available, it is now possible to obtain a comprehensive view on transcriptome complexity and the regulation of transcript diversity.

Epigenetic deregulation in pediatric acute lymphoblastic leukemia

Similar to most cancers, genome-wide DNA methylation profiles are commonly altered in pediatric acute lymphoblastic leukemia (ALL); however, recent observations highlight that a large portion of malignancy-associated DNA methylation alterations are not accompanied by related gene expression changes. By analyzing and integrating the methylome and transcriptome profiles of pediatric B-cell ALL cases and primary tissue controls, we report 325 genes hypermethylated and downregulated and 45 genes hypomethylated and upregulated in pediatric B-cell ALL, irrespective of subtype. Repressed cation channel subunits and cAMP signaling activators and transducers are overrepresented, potentially indicating a reduced cellular potential to receive and propagate apoptotic signals. Furthermore, we report specific DNA methylation alterations with concurrent gene expression changes within individual ALL subtypes. The ETV6-RUNX1 translocation was associated with downregulation of ASNS and upregulation of the EPO-receptor, while Hyperdiploid patients (> 50 chr) displayed upregulation of B-cell lymphoma (BCL) members and repression of PTPRG and FHIT. In combination, these data indicate genetically distinct B-cell ALL subtypes contain cooperative epimutations and genome-wide epigenetic deregulation is common across all B-cell ALL subtypes.

Genome-wide signatures of differential DNA methylation in pediatric acute lymphoblastic leukemia

BACKGROUND

Although aberrant DNA methylation has been observed previously in acute lymphoblastic leukemia (ALL), the patterns of differential methylation have not been comprehensively determined in all subtypes of ALL on a genome-wide scale. The relationship between DNA methylation, cytogenetic background, drug resistance and relapse in ALL is poorly understood.

RESULTS

We surveyed the DNA methylation levels of 435,941 CpG sites in samples from 764 children at diagnosis of ALL and from 27 children at relapse. This survey uncovered four characteristic methylation signatures. First, compared with control blood cells, the methylomes of ALL cells shared 9,406 predominantly hypermethylated CpG sites, independent of cytogenetic background. Second, each cytogenetic subtype of ALL displayed a unique set of hyper- and hypomethylated CpG sites. The CpG sites that constituted these two signatures differed in their functional genomic enrichment to regions with marks of active or repressed chromatin. Third, we identified subtype-specific differential methylation in promoter and enhancer regions that were strongly correlated with gene expression. Fourth, a set of 6,612 CpG sites was predominantly hypermethylated in ALL cells at relapse, compared with matched samples at diagnosis. Analysis of relapse-free survival identified CpG sites with subtype-specific differential methylation that divided the patients into different risk groups, depending on their methylation status.

CONCLUSIONS

Our results suggest an important biological role for DNA methylation in the differences between ALL subtypes and in their clinical outcome after treatment.

Integration of high-resolution methylome and transcriptome analyses to dissect epigenomic changes in childhood acute lymphoblastic leukemia

B-cell precursor acute lymphoblastic leukemia (pre-B ALL) is the most common pediatric cancer. Although the genetic determinants underlying disease onset remain unclear, epigenetic modifications including DNA methylation are suggested to contribute significantly to leukemogenesis. Using the Illumina 450K array, we assessed DNA methylation in matched tumor-normal samples of 46 childhood patients with pre-B ALL, extending single CpG-site resolution analysis of the pre-B ALL methylome beyond CpG-islands (CGI). Unsupervised hierarchical clustering of CpG-site neighborhood, gene, or microRNA (miRNA) gene-associated methylation levels separated the tumor cohort according to major pre-B ALL subtypes, and methylation in CGIs, CGI shores, and in regions around the transcription start site was found to significantly correlate with transcript expression. Focusing on samples carrying the t(12;21) ETV6-RUNX1 fusion, we identified 119 subtype-specific high-confidence marker CpG-loci. Pathway analyses linked the CpG-loci-associated genes with hematopoiesis and cancer. Further integration with whole-transcriptome data showed the effects of methylation on expression of 17 potential drivers of leukemogenesis. Independent validation of array methylation and sequencing-derived transcript expression with Sequenom Epityper technology and real-time quantitative reverse transcriptase PCR, respectively, indicates more than 80% empirical accuracy of our genome-wide findings. In summary, genome-wide DNA methylation profiling enabled us to separate pre-B ALL according to major subtypes, to map epigenetic biomarkers specific for the t(12;21) subtype, and through a combined methylome and transcriptome approach to identify downstream effects on candidate drivers of leukemogenesis.

Genome-scale technology driven advances to research into normal and malignant haematopoiesis

Haematopoiesis or blood development has long served as a model system for adult stem cell biology. Moreover, when combined, the various cancers of the blood represent one of the commonest human malignancies. Large numbers of researchers have therefore dedicated their scientific careers to studying haematopoiesis for more than a century. Throughout this period, many new technologies have first been applied towards the study of blood cells, and the research fields of normal and malignant haematopoiesis have also been some of the earliest adopters of genome-scale technologies. This has resulted in significant new insights with implications ranging from basic biological mechanisms to patient diagnosis and prognosis and also produced lessons likely to be relevant for many other areas of biomedical research. This paper discusses the current state of play for a range of genome-scale applications within haemopoiesis research, including gene expression profiling, ChIP-sequencing, genomewide association analysis, and cancer genome sequencing. A concluding outlook section explores likely future areas of progress as well as potential technological and educational bottlenecks.