Stability of gene expression and epigenetic profiles highlights the utility of patient-derived paediatric acute lymphoblastic leukaemia xenografts for investigating molecular mechanisms of drug resistance
- Nicholas C Wong†1, 2, 4,
- Vivek A Bhadri†3,
- Jovana Maksimovic1,
- Mandy Parkinson-Bates1,
- Jane Ng1,
- Jeff M Craig1, 4,
- Richard Saffery†1, 4 and
- Richard B Lock†3Email author
© Wong et al.; licensee BioMed Central Ltd. 2014
Received: 29 November 2013
Accepted: 20 May 2014
Published: 1 June 2014
Patient-derived tumour xenografts are an attractive model for preclinical testing of anti-cancer drugs. Insights into tumour biology and biomarkers predictive of responses to chemotherapeutic drugs can also be gained from investigating xenograft models. As a first step towards examining the equivalence of epigenetic profiles between xenografts and primary tumours in paediatric leukaemia, we performed genome-scale DNA methylation and gene expression profiling on a panel of 10 paediatric B-cell precursor acute lymphoblastic leukaemia (BCP-ALL) tumours that were stratified by prednisolone response.
We found high correlations in DNA methylation and gene expression profiles between matching primary and xenograft tumour samples with Pearson’s correlation coefficients ranging between 0.85 and 0.98. In order to demonstrate the potential utility of epigenetic analyses in BCP-ALL xenografts, we identified DNA methylation biomarkers that correlated with prednisolone responsiveness of the original tumour samples. Differential methylation of CAPS2, ARHGAP21, ARX and HOXB6 were confirmed by locus specific analysis. We identified 20 genes showing an inverse relationship between DNA methylation and gene expression in association with prednisolone response. Pathway analysis of these genes implicated apoptosis, cell signalling and cell structure networks in prednisolone responsiveness.
The findings of this study confirm the stability of epigenetic and gene expression profiles of paediatric BCP-ALL propagated in mouse xenograft models. Further, our preliminary investigation of prednisolone sensitivity highlights the utility of mouse xenograft models for preclinical development of novel drug regimens with parallel investigation of underlying gene expression and epigenetic responses associated with novel drug responses.
Despite progress in the treatment of several cancers over recent decades, the lack of clinically relevant tumour models for individual subtypes of human cancer has proven to be a major impediment in the development of effective anti-cancer therapies . Approaches that facilitate development of novel rational therapies targeting specific tumours (or specific features of tumours) remain an urgent priority. Traditional models of human cancer involving the analysis of immortalised cell lines have given way in recent years to more clinically relevant studies in models that mirror the features of primary tumours . The two main approaches have been the generation of primary tumour-derived cell lines, and the generation of mouse models, either via transgenic approaches or through the engraftment of primary human tumour into immune-compromised mouse models . Mouse models have been used extensively in this regard, for preclinical testing of drug efficacy and toxicity prior to establishing clinical trials. A broad panel of xenografts with known treatment responsiveness, and well-defined molecular profiles, would provide an excellent adjunct to these models .
Mouse xenograft models of haematological malignancies, established by the transplantation of donor cells into non-obese diabetic/severe combined immunodeficient (NOD/SCID) or NOD/SCID/IL-2 receptor gamma chain−/− (NSG) mice, are recognised as one of the most clinically relevant systems for investigating leukaemia biology and testing new treatments [5–12]. This is due to the faithful recapitulation of many aspects of the human disease, including kinetics of engraftment in the bone marrow (BM), with subsequent infiltration of the spleen, peripheral blood and other organs [10, 13, 14]. For these reasons, patient-derived xenografts (PDXs) are considered superior to in vitro immortalised cancer cell lines that show many differences to primary tumours, including gene expression, drug responsiveness and epigenetic profiles , which is most likely due to the selective processes associated with long term culturing. PDXs have become increasingly popular as evidence mounts that they accurately recapitulate many of the features of patient tumours, such as tumour microenvironment, differentiation state and morphology, architecture and in some instances molecular signatures of the original patient tumour (reviewed in [1, 2]).
To establish the relevance of PDX models to primary tumours, high density molecular profiling of gene expression and epigenetic markers should be performed. This was recently demonstrated for gene expression both between two tissue types, bone marrow and spleen and between independently engrafted mice for T-ALL .
As a first step towards examining the equivalence of epigenetic profiles between primary tumour and xenograft, we carried out parallel DNA methylation and gene expression profiling on a panel of childhood B-cell precursor acute lymphoblastic leukaemia (BCP-ALL) selected by their clinical responses to prednisolone. This panel consisted of five individuals who had a good response to prednisolone (PGR) and five who had a poor response (PPR). By comparing DNA methylation and gene expression profiles between primary and derived, single-passaged xenograft lines, we report the stability of both gene expression and DNA methylation in the xenograft, further highlighting their potential for exploring gene expression and epigenetic changes associated with responses to established and novel drugs.
Patient samples, characteristics and xenograft model generation
Patient demographics of xenografts used in this study
Day 8 blasts
(Diagnostic patient sample)
Genomic DNA and total RNA extraction
Genomic DNA was extracted from the primary bone marrow biopsies used for xenografting and from cells harvested from the spleens of engrafted animals for each xenograft using standard phenol/chloroform extraction and isopropanol precipitation. Total RNA was extracted using TriZol Reagent (Life Technologies, Carlsbad, USA) according to manufacturer’s instructions. Quality and yield were measured using a Nanodrop spectrophotometer.
Sodium bisulphite conversion of genomic DNA
Genomic DNA was converted for DNA methylation analysis using the MethylEasy Xceed Kit (Human Genetic Signatures, Sydney, Australia) according to manufacturer’s instructions. Converted DNA was used for downstream Illumina Infinium DNA methylation BeadArray analysis and SEQUENOM EpitTYPER validation.
Genome-scale DNA methylation analysis
Converted genomic DNA was processed and analysed for Illumina Infinium HumanMethylation27 BeadArray (Illumina, San Diego, USA) according to manufacturer’s instructions (ServiceXS, Leiden, The Netherlands). This BeadArray platform interrogates 27,578 CpG sites across the human genome. The arrays were scanned using an Illumina BeadArray Reader and subsequently processed using the Illumina GenomeStudio V.1 software package. The Bioconductor Lumi package was used for downstream data processing and normalisation . Briefly, DNA probe methylation data were quality checked and then colour balance adjusted, background corrected and scaled based on the mean of all probes, using the methylation simple scaling normalization (SSN) implemented within the Lumi package. CpG sites with at least one sample having a detection p-value > 0.01 were excluded from subsequent analyses, leaving 27,341 CpG sites. Differential methylation analysis was performed using the LIMMA package from Bioconductor . Significantly differentially methylated probes were selected based on a Benjamini-Hochberg adjusted p-value < 0.05. The methylation microarray data have been deposited into Gene Expression Ominibus (http://www.ncbi.nlm.nih.gov/geo/) with the identifier GSE57581.
Gene expression Illumina array analysis
Total RNA was extracted from the primary and xenograft tumours and amplified using the Illumina TotalPrep RNA amplification kit (Ambion, Austin, USA). The amplified total RNA was analysed using Illumina WG-6_V3 chips (Illumina, San Diego, USA) according to manufacturer’s instructions. The sample probe profiles with no normalisation or background correction were exported from BeadStudio (version 3.0.14, Illumina), and the data were pre-processed using quantile normalisation. Probes with detection p-value greater than 0.01 on all arrays were deemed as non-expressed probes and filtered out. Differential gene expression was determined using LIMMA with the positive False Discovery Rate (FDR) correction for multiple testing (Benjamini-Hochberg adjusted p-value < 0.05). The gene expression microarray data have been deposited into Gene Expression Ominibus (http://www.ncbi.nlm.nih.gov/geo/) with the identifier GSE57491.
SEQUENOM MassArray EpiTYPER analysis
Primers (detailed in Additional file 1: Table S1) were designed to generate PCR amplicons from bisulphite converted genomic DNA suitable for SEQUENOM EpiTYPER chemistry as per the manufacturer’s protocol. Samples were analysed using MALDI-TOF mass spectrometry, DNA methylation information was collected using EpiTYPER Viewer Software (v 1.0.5). Non-analysable and poor quality CpG sites were removed from downstream analysis as previously described .
Xenograft models of BCP-ALL are an accurate reflection of DNA methylation and gene expression status of the corresponding primary tumour
To identify differential DNA methylation between primary tumour and matching xenograft samples we applied a linear model with empirical Bayes estimation and found 1564 probes to be differentially methylated between matching primary tumour and xenograft sample after correction for multiple testing (adjusted p-value < 0.05, Additional file 2: Table S2). The majority of these probes demonstrated a small change in DNA methylation with the average difference across individuals ranging from 0.4 to 8.6% (Additional file 3: Figure S1A).
We also looked for differential gene expression between matching primary and xenograft cell lines again applying a linear model with empirical Bayes estimation on the genome-scale gene expression microarray results. We found 3441 probes from 3208 genes to be differentially expressed between primary and xenograft lines (adjusted p-value < 0.05, Additional file 4: Table S3). However, as we observed with DNA methylation, the differences in expression of these probes between primary and xenograft were minimal with an average fold difference in expression between primary and xenograft tumours of 1.12 (Additional file 3: Figure S1B).
Using DAVID (http://david.abcc.ncifcrf.gov/), the differentially methylated and differentially expressed genes between primary and matching xenograft lines were found to be mainly involved in haematological and cell signalling processes that could be accounted for given the cellular origins of the primary (bone marrow) and xenograft (spleen) samples.
Given the relatively small number of differentially methylated probes (6%) and differentially expressed probes (17%), and the minimal absolute differences in DNA methylation and expression (Additional file 3: Figure S1A and S1B), our results indicate that xenograft models largely recapitulate the DNA methylation and gene expression profile of the corresponding primary tumour. This highlights the potential utility of xenograft cell lines for modelling primary disease.
Molecular biomarkers associated with prednisolone response
Differential probes associated with prednisolone response
DNA methylation probes
GLI pathogenesis-related 1 like 2
phosphorylase kinase, gamma 1 (muscle)
interferon regulatory factor 7
apoptosis-inducing factor, mitochondrion-associated, 2
leukocyte immunoglobulin-like receptor, subfamily B (with TM and ITIM domains), member 5
sorting nexin 7
BCL2/adenovirus E1B 19kDa interacting protein 2
DEAH (Asp-Glu-Ala-His) box polypeptide 15
nuclear factor I/X (CCAAT-binding transcription factor)
glycoprotein 2 (zymogen granule membrane)
roundabout homolog 4, magic roundabout (Drosophila)
Rho GTPase activating protein 21
gap junction protein, gamma 2, 47kDa
programmed cell death 1
leucine rich repeat containing 32
solute carrier family 22, member 17
cysteine-rich secretory protein LCCL domain containing 1
aquaporin 3 (Gill blood group)
vesicle-associated membrane protein 8 (endobrevin)
C-type lectin domain family 4, member F
mastermind-like domain containing 1
ORM1-like 3 (S. cerevisiae)
G protein-coupled receptor, family C, group 6, member A
aristaless related homeobox
coiled-coil domain containing 19
Gene expression probes
PRKC, apoptosis, WT1, regulator
Microtubule associated serine/threonine kinase family member 4
Quiescin Q6 sulfhydryl oxidase 2
Neuroblastoma breakpoint family, member 1
Lipopolysaccharide-induced TNF factor
Myosin VIIA and Rab interacting protein
Protein kinase, cAMP-dependent, regulatory, type II, alpha
Zinc fingers and homeoboxes 2
Collagen, type V, alpha 1
Plasminogen activator, urokinase
Sema domain, immunoglobulin domain (Ig), short basic domain, secreted, (semaphorin) 3F
Homo sapiens rhophilin-like protein mRNA, complete cds.
Adenylate cyclase 4
ArfGAP with RhoGAP domain, ankyrin repeat and PH domain 3
1-acylglycerol-3-phosphate O-acyltransferase 2 (lysophosphatidic acid acyltransferase, beta)
SH3 domain containing ring finger 3
Adenomatosis polyposis coli down-regulated 1
Zeta-chain (TCR) associated protein kinase 70kDa
Sperm associated antigen 16
Probes both differentially methylated and expressed in association with prednisolone response
Upregulated and less methylated in PPR
H2A histone family, member Y2
(>2-fold expression, <-0.2 methylation)
H2A histone family, member Y2
PRKC, apoptosis, WT1, regulator
PRKC, apoptosis, WT1, regulator
MOCO sulphurase C-terminal domain containing 1
Nerve growth factor receptor (TNFRSF16) associated protein 1
Myristoylated alanine-rich protein kinase C substrate
Chemokine (C-C motif) receptor 7
Plastin 3 (T isoform)
Downregulated and more methylated in PPR
POU class 4 homeobox 1
(<-2-fold expression, >0.2 methylation)
Cytochrome b5 reductase 2
Transmembrane emp24 protein transport domain containing 6
Collapsin response mediator protein 1
Collapsin response mediator protein 1
Iroquois homeobox 3
Leucine zipper, down-regulated in cancer 1
Validation of DNA methylation biomarkers associated with prednisolone response
It is becoming clear that the complexity of genetic, epigenetic, and subsequent gene expression disruption associated with human cancer is immense. As such, many mouse models of tumourigenesis are limited in their capacity to faithfully mimic human disease. In light of this, patient derived tumour tissue xenograft models are increasingly recognised as offering the most robust approach for testing tumour responses to various chemotherapeutic regimens, evaluating the efficacy of novel therapeutic agents, analysing the process of tumour progression at the cellular and molecular level and the identification of new therapeutic targets . However, as with most mouse xenograft models, the stability of molecular profiles (gene expression and epigenetic) that regulate all aspects of tumour function remains to be determined. Confirmation of this stability is crucial in order identify molecular responses to treatment within the xenograft that could be extrapolated back to patients.
Here, we have determined the stability of genome wide DNA methylation and gene expression profiles between primary tumour cells and matching xenograft tumour cells from a small number of paediatric ALL cases with differential response to prednisolone. A high correlation in both DNA methylation and gene expression profiles was observed in all cases, confirming the stability of these molecular features of primary tumours in the mouse system. Differences in DNA methylation and gene expression between primary and xenograft samples were negligible in magnitude (Additional file 3: Figure S1) and comprised of a small fraction of probes for each array platform. The differentially methylated genes include MYOD1, GPR6 and SLC27A6 (Table 1). Many genes associated with minor expression differences were part of the globin gene family and genes involved in oxygen transport and include HBB, AHSP, HBD, HBA2 (Table 2). This is likely to have arisen by the differences in cellular composition as the primary tumour samples contained a milieu of haematopoietic cells, including human erythrocytes that were absent in the xenograft samples that comprised of mononuclear cells derived from the murine spleen. Given the high degree of correlation and clustering of matching primary and xenograft samples after unsupervised hierarchical clustering of the most varied probes for DNA methylation and gene expression, the xenografts described in this study are an accurate reflection of their corresponding primary tumours.
While a number of candidate genes whose DNA methylation and/or gene expression status were associated with prednisolone response, given the small sample numbers and inherent genetic heterogeneity of the tumours, the significance of these genes remains unclear. Using hierarchical clustering, the most significant probes for DNA methylation discriminated prednisolone response while the gene expression probes did not (Figure 2), reflecting the more variable nature of gene expression compared to DNA methylation [21, 22]. Using SEQUENOM, we were able to replicate DNA methylation changes at four genes associated with prednisolone response indicative of a potential DNA methylation biomarker. Taking methylation and expression status together, 20 genes were differentially regulated between good and poor responders to prednisolone (Table 3). While the genes were found to be part of apoptotic and cell signalling pathways, their significance remains unclear given the small numbers in each group. PAWR demonstrated significant overexpression and hypomethylation across PPRs compared to PGRs. This is a WT1 interacting protein that also functions as a transcriptional repressor with pro-apoptotic functions and tumour resistance . While the down regulation of PAWR confers poor prognosis in a range of solid tumours [24, 25], its role in haematological malignancy is less clear, with expression detectable in a range of leukaemias . Our results warrant further investigation of PAWR to determine a potential role in prednisolone response and responses to other novel drug regimens in an expanded xenograft cohort.
Another gene with potential interest is POU4F1, which appears to be differentially regulated according to prednisolone response (Table 3). However in our analysis, statistical significance was not achieved with the modest sample size of our panel. POU4F1 has been shown to have a role in regulating the expression of B-cell markers in t(8;21) positive acute myeloid leukaemia [27–29]. Its role in B-cell ALL response to prednisolone remains unclear and could be a potential gene target for further characterisation in an expanded B-cell ALL xenograft panel.
While our study did not identify statistically significant genes associated with prednisolone response, we present here a first pass analysis using low-resolution microarray platforms to interrogate DNA methylation and gene expression across our model system. We demonstrate that our B-cell ALL xenograft panel recapitulates the DNA methylation and gene expression profiles of the primary tumour and will facilitate future genome-wide interrogation of gene expression and DNA methylation using next generation sequencing methodology.
Patient-derived tumour xenograft models offer superior utility as preclinical models over cell line systems with their ability to recapitulate the milieu and microenvironment of the primary tumour. However, the extent of gene expression and epigenetic stability within the xenograft has remained unclear at least in the haematological setting. We have demonstrated that the gene expression and DNA methylation profiles of cells taken from the spleens of engrafted mice are highly correlated to the original primary tumour. Given the similarity to the primary tumour, our study confirms the opportunity to investigate gene expression and DNA methylation biomarkers in response to novel treatment strategies.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its additional files. All microarray data presented in this paper have been deposited into Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) with the identifiers GSE57581 and GSE57491.
This work was supported by an NHMRC Project Grant to JC, RS and RBL, and an NHMRC Research Fellowship to RBL. NW has been supported by the Leukaemia Foundation of Australia, My Room and the Children’s Cancer Centre Foundation. VB was supported by a scholarship from the Leukaemia Foundation of Australia. The Murdoch Childrens Research Institute is supported by the Victorian Operational Infrastructure Grant. We thank Dr Benjamin Ong for assistance with the SEQUENOM Facility. Children’s Cancer Institute Australia for Medical Research is affiliated with the University of New South Wales and the Sydney Children’s Hospitals Network.
- Williams SA, Anderson WC, Santaguida MT, Dylla SJ: Patient-derived xenografts, the cancer stem cell paradigm, and cancer pathobiology in the 21st century. Lab Invest. 2013, 93: 970-982. 10.1038/labinvest.2013.92.PubMedView ArticleGoogle Scholar
- Tentler JJ, Tan A-C, Weekes CD, Jimeno A, Leong S, Pitts TM, Arcaroli JJ, Messersmith WA, Eckhardt SG: Patient-derived tumour xenografts as models for oncology drug development. Nat Rev Clin Oncol. 2012, 9: 338-350. 10.1038/nrclinonc.2012.61.PubMed CentralPubMedView ArticleGoogle Scholar
- Vandamme T: Use of rodents as models of human diseases. J Pharm Bioallied Sci. 2014, 6: 2-10.4103/0975-7406.124301.PubMed CentralPubMedView ArticleGoogle Scholar
- Bachmann PS, Lock RB: In vivo models of childhood leukemia for preclinical drug testing. Curr Drug Targets. 2007, 8: 773-783. 10.2174/138945007780830809.PubMedView ArticleGoogle Scholar
- Anderson K, Lutz C, Van Delft FW, Bateman CM, Guo Y, Colman SM, Kempski H, Moorman AV, Titley I, Swansbury J, Kearney L, Enver T, Greaves M: Genetic variegation of clonal architecture and propagating cells in leukaemia. Nature. 2010, 469: 356-361.PubMedView ArticleGoogle Scholar
- Carol H, Boehm I, Reynolds CP, Kang MH, Maris JM, Morton CL, Gorlick R, Kolb EA, Keir ST, Wu J, Wozniak AE, Yang Y, Manfredi M, Ecsedy J, Wang J, Neale G, Houghton PJ, Smith MA, Lock RB: Efficacy and pharmacokinetic/pharmacodynamic evaluation of the Aurora kinase A inhibitor MLN8237 against preclinical models of pediatric cancer. Cancer Chemother Pharmacol. 2011, 68: 1291-1304. 10.1007/s00280-011-1618-8.PubMed CentralPubMedView ArticleGoogle Scholar
- Carol H, Szymanska B, Evans K, Boehm I, Houghton PJ, Smith MA, Lock RB: The anti-CD19 antibody-drug conjugate SAR3419 prevents hematolymphoid relapse postinduction therapy in preclinical models of pediatric acute lymphoblastic leukemia. Clin Cancer Res. 2013, 19: 1795-1805. 10.1158/1078-0432.CCR-12-3613.PubMed CentralPubMedView ArticleGoogle Scholar
- Clappier E, Gerby B, Sigaux F, Delord M, Touzri F, Hernandez L, Ballerini P, Baruchel A, Pflumio F, Soulier J: Clonal selection in xenografted human T cell acute lymphoblastic leukemia recapitulates gain of malignancy at relapse. J Exp Med. 2011, 208: 653-661. 10.1084/jem.20110105.PubMed CentralPubMedView ArticleGoogle Scholar
- Liem NLM, Papa RA, Milross CG, Schmid MA, Tajbakhsh M, Choi S, Ramirez CD, Rice AM, Haber M, Norris MD, MacKenzie KL, Lock RB: Characterization of childhood acute lymphoblastic leukemia xenograft models for the preclinical evaluation of new therapies. Blood. 2004, 103: 3905-3914. 10.1182/blood-2003-08-2911.PubMedView ArticleGoogle Scholar
- Lock RB, Liem N, Farnsworth ML, Milross CG, Xue C, Tajbakhsh M, Haber M, Norris MD, Marshall GM, Rice AM: The nonobese diabetic/severe combined immunodeficient (NOD/SCID) mouse model of childhood acute lymphoblastic leukemia reveals intrinsic differences in biologic characteristics at diagnosis and relapse. Blood. 2002, 99: 4100-4108. 10.1182/blood.V99.11.4100.PubMedView ArticleGoogle Scholar
- Maude SL, Tasian SK, Vincent T, Hall JW, Sheen C, Roberts KG, Seif AE, Barrett DM, Chen IM, Collins JR, Mullighan CG, Hunger SP, Harvey RC, Willman CL, Fridman JS, Loh ML, Grupp SA, Teachey DT: Targeting JAK1/2 and mTOR in murine xenograft models of Ph-like acute lymphoblastic leukemia. Blood. 2012, 120: 3510-3518. 10.1182/blood-2012-03-415448.PubMed CentralPubMedView ArticleGoogle Scholar
- Notta F, Mullighan CG, Wang JCY, Poeppl A, Doulatov S, Phillips LA, Ma J, Minden MD, Downing JR, Dick JE: Evolution of human BCR-ABL1 lymphoblastic leukaemia-initiating cells. Nature. 2011, 469: 362-367. 10.1038/nature09733.PubMedView ArticleGoogle Scholar
- Kamel-Reid S, Letarte M, Doedens M, Greaves A, Murdoch B, Grunberger T, Lapidot T, Thorner P, Freedman MH, Phillips RA: Bone marrow from children in relapse with pre-B acute lymphoblastic leukemia proliferates and disseminates rapidly in scid mice. Blood. 1991, 78: 2973-2981.PubMedGoogle Scholar
- Nijmeijer BA, Mollevanger P, Van Zelderen-Bhola SL, Kluin-Nelemans HC, Willemze R, Falkenburg JH: Monitoring of engraftment and progression of acute lymphoblastic leukemia in individual NOD/SCID mice. Exp Hematol. 2001, 29: 322-329. 10.1016/S0301-472X(00)00669-X.PubMedView ArticleGoogle Scholar
- Saferali A, Grundberg E, Berlivet S, Beauchemin H, Morcos L, Polychronakos C, Pastinen T, Graham J, McNeney B, Naumova AK: Cell culture-induced aberrant methylation of the imprinted IG DMR in human lymphoblastoid cell lines. Epigenetics. 2010, 5: 50-60. 10.4161/epi.5.1.10436.PubMedView ArticleGoogle Scholar
- Samuels AL, Peeva VK, Papa RA, Firth MJ, Francis RW, Beesley AH, Lock RB, Kees UR: Validation of a mouse xenograft model system for gene expression analysis of human acute lymphoblastic leukaemia. BMC Genomics. 2010, 11: 256-10.1186/1471-2164-11-256.PubMed CentralPubMedView ArticleGoogle Scholar
- Bhadri VA, Cowley MJ, Kaplan W, Trahair TN, Lock RB: Evaluation of the NOD/SCID xenograft model for glucocorticoid-regulated gene expression in childhood B-cell precursor acute lymphoblastic leukemia. BMC Genomics. 2011, 12: 565-10.1186/1471-2164-12-565.PubMed CentralPubMedView ArticleGoogle Scholar
- Du P, Kibbe WA, Lin SM: lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008, 24: 1547-1548. 10.1093/bioinformatics/btn224.PubMedView ArticleGoogle Scholar
- Wettenhall JM, Smyth GK: limmaGUI: A graphical user interface for linear modeling of microarray data. Bioinformatics. 2004, 20: 3705-3706. 10.1093/bioinformatics/bth449.PubMedView ArticleGoogle Scholar
- Ollikainen M, Smith KR, Joo EJ-H, Ng H-K, Andronikos R, Novakovic B, Abdul Aziz NK, Carlin JB, Morley R, Saffery R, Craig JM: DNA methylation analysis of multiple tissues from newborn twins reveals both genetic and intrauterine components to variation in the human neonatal epigenome. Hum Mol Genet. 2010, 19: 4176-4188. 10.1093/hmg/ddq336.PubMedView ArticleGoogle Scholar
- Gutierrez-Arcelus M, Lappalainen T, Montgomery SB, Buil A, Ongen H, Yurovsky A, Bryois J, Giger T, Romano L, Planchon A, Falconnet E, Bielser D, Gagnebin M, Padioleau I, Borel C, Letourneau A, Makrythanasis P, Guipponi M, Gehrig C, Antonarakis SE, Dermitzakis ET: Passive and active DNA methylation and the interplay with genetic variation in gene regulation. Elife. 2013, 2: e00523-e00523. 10.7554/eLife.00523.PubMed CentralPubMedGoogle Scholar
- Pai AA, Bell JT, Marioni JC, Pritchard JK, Gilad Y: A Genome-wide study of DNA methylation patterns and gene expression levels in multiple human and chimpanzee tissues. PLoS Genet. 2011, 7: e1001316-10.1371/journal.pgen.1001316.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhao Y, Rangnekar VM: Apoptosis and tumor resistance conferred by Par-4. Cancer Biol Ther. 2008, 7: 1867-1874. 10.4161/cbt.7.12.6945.PubMed CentralPubMedView ArticleGoogle Scholar
- Nagai MA, Gerhard R, Salaorni S, Fregnani JHTG, Nonogaki S, Netto MM, Soares FA: Down-regulation of the candidate tumor suppressor gene PAR-4 is associated with poor prognosis in breast cancer. Int J Oncol. 2010, 37: 41-49.PubMedView ArticleGoogle Scholar
- Moreno-Bueno G, Fernandez-Marcos PJ, Collado M, Tendero MJ, Rodriguez-Pinilla SM, Garcia-Cao I, Hardisson D, Diaz-Meco MT, Moscat J, Serrano M, Palacios J: Inactivation of the candidate tumor suppressor par-4 in endometrial cancer. Cancer Res. 2007, 67: 1927-1934. 10.1158/0008-5472.CAN-06-2687.PubMedView ArticleGoogle Scholar
- Boehrer S, Chow KU, Puccetti E, Ruthardt M, Godzisard S, Krapohl A, Schneider B, Hoelzer D, Mitrou PS, Rangnekar VM, Weidmann E: Deregulated expression of prostate apoptosis response gene-4 in less differentiated lymphocytes and inverse expressional patterns of par-4 and bcl-2 in acute lymphocytic leukemia. Hematol J. 2001, 2: 103-107. 10.1038/sj.thj.6200089.PubMedView ArticleGoogle Scholar
- Fortier JM, Payton JE, Cahan P, Ley TJ, Walter MJ, Graubert TA: POU4F1 is associated with t(8;21) acute myeloid leukemia and contributes directly to its unique transcriptional signature. Leukemia. 2010, 24: 950-957. 10.1038/leu.2010.61.PubMed CentralPubMedView ArticleGoogle Scholar
- Dunne J, Gascoyne DM, Lister TA, Brady HJM, Heidenreich O, Young BD: AML1/ETO proteins control POU4F1/BRN3A expression and function in t(8;21) acute myeloid leukemia. Cancer Res. 2010, 70: 3985-3995. 10.1158/0008-5472.CAN-09-3604.PubMed CentralPubMedView ArticleGoogle Scholar
- Dunne J, Mannari D, Farzaneh T, Gessner A, Van Delft FW, Heidenreich O, Young BD, Gascoyne DM: AML1/ETO and POU4F1 synergy drives B-lymphoid gene expression typical of t(8;21) acute myeloid leukemia. Leukemia. 2012, 26: 1131-1135. 10.1038/leu.2011.316.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.