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

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-416) contains supplementary material, which is available to authorized users.     

Promoter methylation analysis of WNT/β-catenin pathway regulators and its association with expression of DNMT1 enzyme in colorectal cancer

Background

Aberrant DNA methylation as the most important reason making epigenetic silencing of genes is a main mechanism of gene inactivation in patients with colorectal cancer. In this study, we decided to identify promoter methylation status of ten genes encoding WNT negative regulators, and measure the expression of DNMT1 enzyme in colorectal cancer samples.

Results

Aberrant methylation of APC gene was statistically significant associated with age over 50 (p = 0.017), DDK3 with male (p < 0.0001), SFRP4, WIF1, and WNT5a with increasing tumor stage (p = 0.004, p = 0.029, and p = 0.004), SFRP4 and WIF1 with tumor differentiation (p = 0.009 and p = 0.031) and SFRP2 and SFRP5 with histological type (p = 0.001 and p = 0.025). The increasing number of methylated genes correlated with the expression levels of the DNMT1 mRNA.

Conclusions

The rate of gene promoter methylation of WNT pathway regulators is high in colorectal cancer cells. Hyper-methylation is associated with increased expression of the DNMT1 enzyme.     

Epigenome-wide association study reveals longitudinally stable DNA methylation differences in CD4+ T cells from children with IgE-mediated food allergy

Food allergy is mediated by a combination of genetic and environmental risk factors, potentially mediated by epigenetic mechanisms. CD4+ T-cells are key drivers of the allergic response, and may therefore harbor epigenetic variation in association with the disease phenotype. Here we retrospectively examined genome-wide DNA methylation profiles (~450 000 CpGs) from CD4+ T-cells on a birth cohort of 12 children with IgE-mediated food allergy diagnosed at 12-months, and 12 non-allergic controls. DNA samples were available at two time points, birth and 12-months. Case:control comparisons of CD4+ methylation profiles identified 179 differentially methylated probes (DMP) at 12-months and 136 DMP at birth (FDR-adjusted P value < 0.05, delta β > 0.1). Approximately 30% of DMPs were coincident with previously annotated SNPs. A total of 96 allergy-associated non-SNP DMPs were present at birth when individuals were initially disease-free, potentially implicating these loci in the causal pathway. Pathway analysis of differentially methylated genes identified several MAP kinase signaling molecules. Mass spectrometry was used to validate 15 CpG sites at 3 candidate genes. Combined analysis of differential methylation with gene expression profiles revealed gene expression differences at some but not all allergy associated differentially methylated genes. Thus, dysregulation of DNA methylation at MAPK signaling-associated genes during early CD4+ T-cell development may contribute to suboptimal T-lymphocyte responses in early childhood associated with the development of food allergy.     

Genome-wide DNA methylation variability in adolescent monozygotic twins followed since birth

DNA methylation patterns are characterized by highly conserved developmental programs, but allow for divergent gene expression resulting from stochastic epigenetic drift or divergent environments. Genome-wide methylation studies in monozygotic (MZ) twins are providing insight into the extent of epigenetic variation that occurs, irrespective of genotype. However, little is known about the variability of DNA methylation patterns in adolescence, a period involving significant and rapid physical, emotional, social, and neurodevelopmental change. Here, we assessed genome-wide DNA methylation using the 450 K Illumina BeadChip in a sample of 37 MZ twin pairs followed longitudinally since birth to investigate: 1) the extent of variation in DNA methylation in identical genetic backgrounds in adolescence and; 2) whether these variations are randomly distributed or enriched in particular functional pathways. We also assessed stability of DNA methylation over 3–6 months to distinguish stable trait-like and variable state-like genes. A pathway analysis found high within-pair variability in genes associated with development, cellular mechanisms, tissue and cell morphology, and various disorders. Test-retest analyses performed in a sub-sample of 8 twin pairs demonstrated enrichment in gene pathways involved in organismal development, cellular growth and proliferation, cell signaling, and particular disorders. The variability found in functional gene pathways may plausibly underlie phenotypic differences in this adolescent MZ twin sample. Furthermore, we assessed stability of methylation over 3–6 months and found that some genes were stable while others were unstable, suggesting that the methylome remains dynamic in adolescence and that dynamic sites tend to be organized in certain gene pathways.     

Methylome sequencing for fibrolamellar hepatocellular carcinoma depicts distinctive features

With the goal of studying epigenetic alterations in fibrolamellar hepatocellular carcinoma (FLC) and establish an associated DNA methylation signature, we analyzed LINE-1 methylation in a cohort of FLC and performed next-generation sequencing of DNA methylation in a training set of pure-FLCs and non-cirrhotic hepatocellular carcinomas (nc-HCC). DNA methylation was correlated with gene expression. Furthermore, we established and validated an epigenetic signature differentiating pure-FLC from other HCCs. LINE-1 methylation correlated with shorter recurrence-free survival and overall survival in resected pure-FLC patients. Unsupervised clustering using CG sites located in islands distinguished pure-FLC from nc-HCC. Major DNA methylation changes occurred outside promoters, mainly in gene bodies and intergenic regions located in the vicinity of liver developmental genes (i.e., SMARCA4 and RXRA). Partially methylated domains were more prone to DNA methylation changes. Furthermore, we identified several putative tumor suppressor genes (e.g., DLEU7) and oncogenes (e.g., DUSP4). While ∼70% of identified gene promoters gaining methylation were marked by bivalent histone marks (H3K4me3/H3K27me3) in embryonic stem cells, ∼70% of those losing methylation were marked by H3K4me3. Finally, we established a pure FLC DNA methylation signature and validated it in an independent dataset. Our analysis reveals a distinct epigenetic signature of pure FLC as compared to nc-HCC, with DNA methylation changes occurring in the vicinity of liver developmental genes. These data suggest new options for targeting FLC based on cancer epigenome aberrations.     

Genome-wide methylation and expression differences in HPV(+) and HPV(?) squamous cell carcinoma cell lines are consistent with divergent mechanisms of carcinogenesis

Oncogenic human papillomaviruses (HPV) are associated with nearly all cervical cancers and are increasingly important in the etiology of oropharyngeal tumors. HPV-associated head and neck squamous cell carcinomas (HNSCC) have distinct risk profiles and appreciate a prognostic advantage compared to HPV-negative HNSCC. Promoter hypermethylation is widely recognized as a mechanism in the progression of HNSCC, but the extent to which this mechanism is consistent between HPV(+) and HPV(−) tumors is unknown. To investigate the epigenetic regulation of gene expression in HPV-induced and carcinogen-induced cancers, we examined genome-wide DNA methylation and gene expression in HPV(+) and HPV(−) SCC cell lines. We used two platforms: the Illumina Infinium Methylation BeadArray and tiling arrays, and confirmed illustrative examples with pyrosequencing and quantitative PCR. These analyses indicate that HPV(+) cell lines have higher DNA methylation in genic and LINE-1 regions than HPV(−) cell lines. Differentially methylated loci between HPV(+) and HPV(−) cell lines significantly correlated with HPV-typed HNSCC primary tumor DNA methylation levels. Novel findings include higher promoter methylation of polycomb repressive complex 2 target genes in HPV(+) cells compared to HPV(−) cells and increased expression of DNMT3A in HPV(+) cells. Additionally, CDKN2A and KRT8 were identified as interaction hubs among genes with higher methylation and lower expression in HPV(−) cells. Conversely, RUNX2, IRS-1 and CCNA1 were major hubs with higher methylation and lower expression in HPV(+) cells. Distinct HPV(+) and HPV(−) epigenetic profiles should provide clues to novel targets for development of individualized therapeutic strategies.Key words: epigenetics, human papillomavirus, HNSCC, DNA methylation, squamous cell carcinoma, gene expression, microarrays, illumina infinium humanmethylation27 beadarray     

DNA methylation-associated colonic mucosal immune and defense responses in treatment-na?ve pediatric ulcerative colitis

Inflammatory bowel diseases (IBD) are emerging globally, indicating that environmental factors may be important in their pathogenesis. Colonic mucosal epigenetic changes, such as DNA methylation, can occur in response to the environment and have been implicated in IBD pathology. However, mucosal DNA methylation has not been examined in treatment-naïve patients. We studied DNA methylation in untreated, left sided colonic biopsy specimens using the Infinium HumanMethylation450 BeadChip array. We analyzed 22 control (C) patients, 15 untreated Crohn’s disease (CD) patients, and 9 untreated ulcerative colitis (UC) patients from two cohorts. Samples obtained at the time of clinical remission from two of the treatment-naïve UC patients were also included into the analysis. UC-specific gene expression was interrogated in a subset of adjacent samples (5 C and 5 UC) using the Affymetrix GeneChip PrimeView Human Gene Expression Arrays. Only treatment-naïve UC separated from control. One-hundred-and-twenty genes with significant expression change in UC (> 2-fold, P &lt; 0.05) were associated with differentially methylated regions (DMRs). Epigenetically associated gene expression changes (including gene expression changes in the IFITM1, ITGB2, S100A9, SLPI, SAA1, and STAT3 genes) were linked to colonic mucosal immune and defense responses. These findings underscore the relationship between epigenetic changes and inflammation in pediatric treatment-naïve UC and may have potential etiologic, diagnostic, and therapeutic relevance for IBD.     

A novel approach to the discovery of survival biomarkers in glioblastoma using a joint analysis of DNA methylation and gene expression

Glioblastoma multiforme (GBM) is the most aggressive of all brain tumors, with a median survival of less than 1.5 years. Recently, epigenetic alterations were found to play key roles in both glioma genesis and clinical outcome, demonstrating the need to integrate genetic and epigenetic data in predictive models. To enhance current models through discovery of novel predictive biomarkers, we employed a genome-wide, agnostic strategy to specifically capture both methylation-directed changes in gene expression and alternative associations of DNA methylation with disease survival in glioma. Human GBM-associated DNA methylation, gene expression, IDH1 mutation status, and survival data were obtained from The Cancer Genome Atlas. DNA methylation loci and expression probes were paired by gene, and their subsequent association with survival was determined by applying an accelerated failure time model to previously published alternative and expression-based association equations. Significant associations were seen in 27 unique methylation/expression pairs with expression-based, alternative, and combinatorial associations observed (10, 13, and 4 pairs, respectively). The majority of the predictive DNA methylation loci were located within CpG islands, and all but three of the locus pairs were negatively correlated with survival. This finding suggests that for most loci, methylation/expression pairs are inversely related, consistent with methylation-associated gene regulatory action. Our results indicate that changes in DNA methylation are associated with altered survival outcome through both coordinated changes in gene expression and alternative mechanisms. Furthermore, our approach offers an alternative method of biomarker discovery using a priori gene pairing and precise targeting to identify novel sites for locus-specific therapeutic intervention.     

De novo assembly and delivery to mouse cells of a 101 kb functional human gene

Design and large-scale synthesis of DNA has been applied to the functional study of viral and microbial genomes. New and expanded technology development is required to unlock the transformative potential of such bottom-up approaches to the study of larger mammalian genomes. Two major challenges include assembling and delivering long DNA sequences. Here, we describe a workflow for de novo DNA assembly and delivery that enables functional evaluation of mammalian genes on the length scale of 100 kilobase pairs (kb). The DNA assembly step is supported by an integrated robotic workcell. We demonstrate assembly of the 101 kb human HPRT1 gene in yeast from 3 kb building blocks, precision delivery of the resulting construct to mouse embryonic stem cells, and subsequent expression of the human protein from its full-length human gene in mouse cells. This workflow provides a framework for mammalian genome writing. We envision utility in producing designer variants of human genes linked to disease and their delivery and functional analysis in cell culture or animal models.     

ChIPing the cistrome of PXR in mouse liver

The pregnane X receptor (PXR) is a key regulator of xenobiotic metabolism and disposition in liver. However, little is known about the PXR DNA-binding signatures in vivo, or how PXR regulates novel direct targets on a genome-wide scale. Therefore, we generated a roadmap of hepatic PXR bindings in the entire mouse genome [chromatin immunoprecipitation (ChIP)-Seq]. The most frequent PXR DNA-binding motif is the AGTTCA-like direct repeat with a 4bp spacer [direct repeat (DR)-4)]. Surprisingly, there are also high motif occurrences with spacers of a periodicity of 5 bp, forming a novel DR-(5n + 4) pattern for PXR binding. PXR-binding overlaps with the epigenetic mark for gene activation (histone-H3K4-di-methylation), but not with epigenetic marks for gene suppression (DNA methylation or histone-H3K27-tri-methylation) (ChIP-on-chip). After administering a PXR agonist, changes in mRNA of most PXR-direct target genes correlate with increased PXR binding. Specifically, increased PXR binding triggers the trans-activation of critical drug-metabolizing enzymes and transporters. The mRNA induction of these genes is absent in PXR-null mice. The current work provides the first in vivo evidence of PXR DNA-binding signatures in the mouse genome, paving the path for predicting and further understanding the multifaceted roles of PXR in liver.