Prediction of outcomes for high-count monoclonal B lymphocytosis using an epigenetic and immunogenetic signature.

ABSTRACT: Monoclonal B-cell lymphocytosis (MBL) progresses to chronic lymphocytic leukemia (CLL) requiring therapy at 1% to 5% per year. Improved prediction of progression would greatly benefit individuals with MBL. Patients with CLL separate into 3 distinct epigenetic subtypes (epitypes) with high prognostic significance, and recently the intermediate epitype has been shown to be enriched for high-risk immunoglobulin lambda variable (IGLV) 3-21 rearrangements, impacting outcomes for these patients. Here, we employed this combined strategy to generate the epigenetic and light chain immunoglobulin (ELCLV3-21) signature to classify 219 individuals with MBL. The ELCLV3-21 high-risk signature distinguished MBL individuals with a high probability of progression (39.9% and 71.1% at 5 and 10 years, respectively). ELCLV3-21 improved the accuracy of predicting time to therapy for individuals with MBL compared with other established prognostic indicators, including the CLL international prognostic index (c-statistic, 0.767 vs 0.668, respectively). Comparing ELCLV3-21 risk groups in MBL vs a cohort of 226 patients with CLL revealed ELCLV3-21 high-risk individuals with MBL had significantly shorter time to therapy (P = .003) and reduced overall survival (P = .03) compared with ELCLV3-21 low-risk individuals with CLL. These results highlight the power of the ELCLV3-21 approach to identify individuals with a higher likelihood of adverse clinical outcome and may provide a more accurate approach to classify individuals with small B-cell clones.

DNA methylation epitypes highlight underlying developmental and disease pathways in acute myeloid leukemia.

Acute myeloid leukemia (AML) is a molecularly complex disease characterized by heterogeneous tumor genetic profiles and involving numerous pathogenic mechanisms and pathways. Integration of molecular data types across multiple patient cohorts may advance current genetic approaches for improved subclassification and understanding of the biology of the disease. Here, we analyzed genome-wide DNA methylation in 649 AML patients using Illumina arrays and identified a configuration of 13 subtypes (termed "epitypes") using unbiased clustering. Integration of genetic data revealed that most epitypes were associated with a certain recurrent mutation (or combination) in a majority of patients, yet other epitypes were largely independent. Epitypes showed developmental blockage at discrete stages of myeloid differentiation, revealing epitypes that retain arrested hematopoietic stem-cell-like phenotypes. Detailed analyses of DNA methylation patterns identified unique patterns of aberrant hyper- and hypomethylation among epitypes, with variable involvement of transcription factors influencing promoter, enhancer, and repressed regions. Patients in epitypes with stem-cell-like methylation features showed inferior overall survival along with up-regulated stem cell gene expression signatures. We further identified a DNA methylation signature involving STAT motifs associated with FLT3-ITD mutations. Finally, DNA methylation signatures were stable at relapse for the large majority of patients, and rare epitype switching accompanied loss of the dominant epitype mutations and reversion to stem-cell-like methylation patterns. These results show that DNA methylation-based classification integrates important molecular features of AML to reveal the diverse pathogenic and biological aspects of the disease.

Developmental subtypes assessed by DNA methylation-iPLEX forecast the natural history of chronic lymphocytic leukemia.

Alterations in global DNA methylation patterns are a major hallmark of cancer and represent attractive biomarkers for personalized risk stratification. Chronic lymphocytic leukemia (CLL) risk stratification studies typically focus on time to first treatment (TTFT), time to progression (TTP) after treatment, and overall survival (OS). Whereas TTFT risk stratification remains similar over time, TTP and OS have changed dramatically with the introduction of targeted therapies, such as the Bruton tyrosine kinase inhibitor ibrutinib. We have shown that genome-wide DNA methylation patterns in CLL are strongly associated with phenotypic differentiation and patient outcomes. Here, we developed a novel assay, termed methylation-iPLEX (Me-iPLEX), for high-throughput quantification of targeted panels of single cytosine guanine dinucleotides from multiple independent loci. Me-iPLEX was used to classify CLL samples into 1 of 3 known epigenetic subtypes (epitypes). We examined the impact of epitype in 1286 CLL patients from 4 independent cohorts representing a comprehensive view of CLL disease course and therapies. We found that epitype significantly predicted TTFT and OS among newly diagnosed CLL patients. Additionally, epitype predicted TTP and OS with 2 common CLL therapies: chemoimmunotherapy and ibrutinib. Epitype retained significance after stratifying by biologically related biomarkers, immunoglobulin heavy chain mutational status, and ZAP70 expression, as well as other common prognostic markers. Furthermore, among several biological traits enriched between epitypes, we found highly biased immunogenetic features, including IGLV3-21 usage in the poorly characterized intermediate-programmed CLL epitype. In summary, Me-iPLEX is an elegant method to assess epigenetic signatures, including robust classification of CLL epitypes that independently stratify patient risk at diagnosis and time of treatment.

Lenalidomide-mediated enhanced translation of C/EBPα-p30 protein up-regulates expression of the antileukemic microRNA-181a in acute myeloid leukemia.

Recently, we showed that increased miR-181a expression was associated with improved outcomes in cytogenetically normal acute myeloid leukemia (CN-AML). Interestingly, miR-181a expression was increased in CN-AML patients harboring CEBPA mutations, which are usually biallelic and associate with better prognosis. CEBPA encodes the C/EBPα transcription factor. We demonstrate here that the presence of N-terminal CEBPA mutations and miR-181a expression are linked. Indeed, the truncated C/EBPα-p30 isoform, which is produced from the N-terminal mutant CEBPA gene or from the differential translation of wild-type CEBPA mRNA and is commonly believed to have no transactivation activity, binds to the miR-181a-1 promoter and up-regulates the microRNA expression. Furthermore, we show that lenalidomide, a drug approved for myelodysplastic syndromes and multiple myeloma, enhances translation of the C/EBPα-p30 isoform, resulting in higher miR-181a levels. In xenograft mouse models, ectopic miR-181a expression inhibits tumor growth. Similarly, lenalidomide exhibits antitumorigenic activity paralleled by increased miR-181a expression. This regulatory pathway may explain an increased sensitivity to apoptosis-inducing chemotherapy in subsets of AML patients. Altogether, our data provide a potential explanation for the improved clinical outcomes observed in CEBPA-mutated CN-AML patients, and suggest that lenalidomide treatment enhancing the C/EBPα-p30 protein levels and in turn miR-181a may sensitize AML blasts to chemotherapy.

inv(16)/t(16;16) acute myeloid leukemia with non-type A CBFB-MYH11 fusions associate with distinct clinical and genetic features and lack KIT mutations.

The inv(16)(p13q22)/t(16;16)(p13;q22) in acute myeloid leukemia results in multiple CBFB-MYH11 fusion transcripts, with type A being most frequent. The biologic and prognostic implications of different fusions are unclear. We analyzed CBFB-MYH11 fusion types in 208 inv(16)/t(16;16) patients with de novo disease, and compared clinical and cytogenetic features and the KIT mutation status between type A (n = 182; 87%) and non-type A (n = 26; 13%) patients. At diagnosis, non-type A patients had lower white blood counts (P = .007), and more often trisomies of chromosomes 8 (P = .01) and 21 (P < .001) and less often trisomy 22 (P = .02). No patient with non-type A fusion carried a KIT mutation, whereas 27% of type A patients did (P = .002). Among the latter, KIT mutations conferred adverse prognosis; clinical outcomes of non-type A and type A patients with wild-type KIT were similar. We also derived a fusion-type-associated global gene-expression profile. Gene Ontology analysis of the differentially expressed genes revealed-among others-an enrichment of up-regulated genes involved in activation of caspase activity, cell differentiation and cell cycle control in non-type A patients. We conclude that non-type A fusions associate with distinct clinical and genetic features, including lack of KIT mutations, and a unique gene-expression profile.

Genome-wide methylation profiling in decitabine-treated patients with acute myeloid leukemia.

The outcome of older (≥ 60 years) acute myeloid leukemia (AML) patients is poor, and novel treatments are needed. In a phase 2 trial for older AML patients, low-dose (20 mg/m(2) per day for 10 days) decitabine, a DNA hypomethylating azanucleoside, produced 47% complete response rate with an excellent toxicity profile. To assess the genome-wide activity of decitabine, we profiled pretreatment and post treatment (day 25/course 1) methylomes of marrow samples from patients (n = 16) participating in the trial using deep-sequencing analysis of methylated DNA captured by methyl-binding protein (MBD2). Decitabine significantly reduced global methylation compared with pretreatment baseline (P = .001). Percent marrow blasts did not correlate with global methylation levels, suggesting that hypomethylation was related to the activity of decitabine rather than to a mere decrease in leukemia burden. Hypomethylation occurred predominantly in CpG islands and CpG island-associated regions (P ranged from .03 to .04) A significant concentration (P < .001) of the hypomehtylated CpG islands was found in chromosome subtelomeric regions, suggesting a differential activity of decitabine in distinct chromosome regions. Hypermethylation occurred much less frequently than hypomethylation and was associated with low CpG content regions. Decitabine-related methylation changes were concordant with those previously reported in distinct genes. In summary, our study supports the feasibility of methylome analyses as a pharmacodynamic endpoint for hypomethylating therapies.

miR-3151 interplays with its host gene BAALC and independently affects outcome of patients with cytogenetically normal acute myeloid leukemia.

High BAALC expression levels are associated with poor outcome in cytogenetically normal acute myeloid leukemia (CN-AML) patients. Recently, miR-3151 was discovered in intron 1 of BAALC. To evaluate the prognostic significance of miR-3151 expression levels and to gain insight into the biologic and prognostic interplay between miR-3151 and its host, miR-3151 and BAALC expression were measured in pretreatment blood of 179 CN-AML patients. Gene-expression profiling and miRNA-expression profiling were performed using microarrays. High miR-3151 expression was associated with shorter disease-free and overall survival, whereas high BAALC expression predicted failure of complete remission and shorter overall survival. Patients exhibiting high expression of both miR-3151 and BAALC had worse outcome than patients expressing low levels of either gene or both genes. In gene-expression profiling, high miR-3151 expressers showed down-regulation of genes involved in transcriptional regulation, posttranslational modification, and cancer pathways. Two genes, FBXL20 and USP40, were validated as direct miR-3151 targets. The results of the present study show that high expression of miR-3151 is an independent prognosticator for poor outcome in CN-AML and affects different outcome end points than its host gene, BAALC. The combination of both markers identified a patient subset with the poorest outcome. This interplay between an intronic miR and its host may have important biologic implications.

Global assessment of promoter methylation in a mouse model of cancer identifies ID4 as a putative tumor-suppressor gene in human leukemia.

DNA methylation is associated with malignant transformation, but limitations imposed by genetic variability, tumor heterogeneity, availability of paired normal tissues and methodologies for global assessment of DNA methylation have limited progress in understanding the extent of epigenetic events in the initiation and progression of human cancer and in identifying genes that undergo methylation during cancer. We developed a mouse model of T/natural killer acute lymphoblastic leukemia that is always preceded by polyclonal lymphocyte expansion to determine how aberrant promoter DNA methylation and consequent gene silencing might be contributing to leukemic transformation. We used restriction landmark genomic scanning with this mouse model of preleukemia reproducibly progressing to leukemia to show that specific genomic methylation is associated with only the leukemic phase and is not random. We also identified Idb4 as a putative tumor-suppressor gene that is methylated in most mouse and human leukemias but in only a minority of other human cancers.

Low expression of MN1 associates with better treatment response in older patients with de novo cytogenetically normal acute myeloid leukemia.

Low MN1 expression bestows favorable prognosis in younger adults with cytogenetically normal acute myeloid leukemia (CN-AML), but its prognostic significance in older patients is unknown. We analyzed pretherapy MN1 expression in 140 older (≥ 60 years) de novo CN-AML patients treated on cytarabine/daunorubicin-based protocols. Low MN1 expressers had higher complete remission (CR) rates (P = .001), and longer overall survival (P = .03) and event-free survival (EFS; P = .004). In multivariable models, low MN1 expression was associated with better CR rates and EFS. The impact of MN1 expression on overall survival and EFS was predominantly in patients 70 years of age or older, with low MN1 expressers with mutated NPM1 having the best outcome. The impact of MN1 expression was also observed in the Intermediate-I, but not the Favorable group of the European LeukemiaNet classification, where low MN1 expressers had CR rates and EFS similar to those of Favorable group patients. MN1 expresser-status-associated gene- and microRNA-expression signatures revealed underexpression of drug resistance and adverse outcome predictors, and overexpression of HOX genes and HOX-gene-embedded microRNAs in low MN1 expressers. We conclude that low MN1 expression confers better prognosis in older CN-AML patients and may refine the European LeukemiaNet classification. Biologic features associated with MN1 expression may help identify new treatment targets.

ASXL1 mutations identify a high-risk subgroup of older patients with primary cytogenetically normal AML within the ELN Favorable genetic category.

The associations of mutations in the enhancer of trithorax and polycomb family gene ASXL1 with pretreatment patient characteristics, outcomes, and gene-/microRNA-expression profiles in primary cytogenetically normal acute myeloid leukemia (CN-AML) are unknown. We analyzed 423 adult patients for ASXL1 mutations, other prognostic gene mutations, and gene-/microRNA-expression profiles. ASXL1 mutations were 5 times more common in older (≥ 60 years) patients (16.2%) than those younger than 60 years (3.2%; P < .001). Among older patients, ASXL1 mutations associated with wild-type NPM1 (P < .001), absence of FLT3-internal tandem duplications (P = .002), mutated CEBPA (P = .01), and with inferior complete remission (CR) rate (P = .04), disease-free survival (DFS; P = .03), overall survival (OS; P = .006), and event-free survival (EFS; P = .002). Within the European LeukemiaNet (ELN) genetic categories of older CN-AML, ASXL1 mutations associated with inferior CR rate (P = .02), OS (P < .001), and EFS (P < .001) among ELN Favorable, but not among ELN Intermediate-I patients. Multivariable analyses confirmed associations of ASXL1 mutations with unfavorable CR rate (P = .03), DFS (P < .001), OS (P < .001), and EFS (P < .001) among ELN Favorable patients. We identified an ASXL1 mutation-associated gene-expression signature, but no microRNA-expression signature. This first study of ASXL1 mutations in primary CN-AML demonstrates that ASXL1-mutated older patients, particularly within the ELN Favorable group, have unfavorable outcomes and may be candidates for experimental treatment approaches.

FLT3 internal tandem duplication associates with adverse outcome and gene- and microRNA-expression signatures in patients 60 years of age or older with primary cytogenetically normal acute myeloid leukemia: a Cancer and Leukemia Group B study.

The clinical impact of FLT3-internal tandem duplications (ITDs), an adverse prognostic marker in adults aged < 60 years with primary cytogenetically normal acute myeloid leukemia (CN-AML), requires further investigation in older patients. In CN-AML patients aged ≥ 60 years treated on Cancer and Leukemia Group B frontline trials, we found that FLT3-ITD remained associated with shorter disease-free survival (P < .001; hazard ratio = 2.10) and overall survival (P < .001; hazard ratio = 1.97) in multivariable analyses. This impact on shorter disease-free survival and overall survival was in patients aged 60-69 (P < .001, each) rather than in those aged ≥ 70 years. An FLT3-ITD-associated gene-expression signature revealed overexpression of FLT3, homeobox genes (MEIS1, PBX3, HOXB3), and immunotherapeutic tar-gets (WT1, CD33) and underexpression of leukemia-associated (MLLT3, TAL1) and erythropoiesis-associated (GATA3, EPOR, ANK1, HEMGN) genes. An FLT3-ITD-associated microRNA-expression signature included overexpressed miR-155 and underexpressed miR-144 and miR-451. FLT3-ITD identifies older CN-AML patients with molecular high risk and is associated with gene- and microRNA-expression signatures that provide biologic insights for novel therapeutic approaches.

BAALC and ERG expression levels are associated with outcome and distinct gene and microRNA expression profiles in older patients with de novo cytogenetically normal acute myeloid leukemia: a Cancer and Leukemia Group B study.

BAALC and ERG expression levels are prognostic markers in younger (< 60 years) cytogenetically normal acute myeloid leukemia (CN-AML) adults; their prognostic impact in older (≥ 60 years) patients requires further investigation. We evaluated pretreatment expression of BAALC and ERG in 158 de novo patients treated on cytarabine/daunorubicin-based protocols. The patients were also characterized for other established molecular prognosticators. Low BAALC and ERG expression levels were associated with better outcome in univariable and multivariable analyses. Expression levels of both BAALC and ERG were the only factors significantly associated with overall survival upon multivariable analysis. To gain biological insights, we derived gene expression signatures associated with BAALC and ERG expression in older CN-AML patients. Furthermore, we derived the first microRNA expression signatures associated with the expression of these 2 genes. In low BAALC expressers, genes associated with undifferentiated hematopoietic precursors and unfavorable outcome predictors were down-regulated, whereas HOX genes and HOX-gene-embedded microRNAs were up-regulated. Low ERG expressers presented with down-regulation of genes involved in the DNA-methylation machinery, and up-regulation of miR-148a, which targets DNMT3B. We conclude that in older CN-AML patients, low BAALC and ERG expression associates with better outcome and distinct gene and microRNA expression signatures that could aid in identifying new targets and novel therapeutic strategies for older patients.

Mutations of the Wilms tumor 1 gene (WT1) in older patients with primary cytogenetically normal acute myeloid leukemia: a Cancer and Leukemia Group B study.

We previously reported the adverse prognostic impact of Wilms tumor 1 gene (WT1) mutations in younger adult cytogenetically normal acute myeloid leukemia (CN-AML). Here, we investigated 243 older (> or = 60 years) primary CN-AML patients. WT1 mutated (WT1mut) patients (7%) had FLT3-ITD more frequently (P < .001), lower hemoglobin (P = .01), higher white blood cell count (P = .03) and percentage blood blasts (P = .03), and a shorter overall survival (P = .08) than WT1 wild-type (WT1wt) patients. Comparing older and younger WT1mut patients, they had similar pretreatment characteristics and outcome. By contrast, among WT1wt CN-AML, younger patients had a significantly better outcome. A WT1 mutation-associated gene-expression signature, reported here for the first time, included CD96, a leukemia stem cell-specific marker, and genes involved in gene regulation (eg, MLL, PML, and SNRPN) and in proliferative and metabolic processes (eg, INSR, IRS2, and PRKAA1), supporting the role of mutated WT1 in deregulating multiple homeostatic processes. Our results indicate that WT1mut CN-AML represents a distinct entity with poor treatment response across age groups. This study has been registered at www.clinicaltrials.gov as #NCT00900224.

Identification and Targeting of the Developmental Blockade in Extranodal Natural Killer/T-cell Lymphoma.

Extranodal natural killer/T-cell lymphoma (ENKTL) is an aggressive, rare lymphoma of natural killer (NK) cell origin with poor clinical outcomes. Here we used phenotypic and molecular profiling, including epigenetic analyses, to investigate how ENKTL ontogeny relates to normal NK-cell development. We demonstrate that neoplastic NK cells are stably, but reversibly, arrested at earlier stages of NK-cell maturation. Genes downregulated in the most epigenetic immature tumors were associated with polycomb silencing along with genomic gain and overexpression of EZH2. ENKTL cells exhibited genome-wide DNA hypermethylation. Tumor-specific DNA methylation gains were associated with polycomb-marked regions, involving extensive gene silencing and loss of transcription factor binding. To investigate therapeutic targeting, we treated novel patient-derived xenograft (PDX) models of ENKTL with the DNA hypomethylating agent, 5-azacytidine. Treatment led to reexpression of NK-cell developmental genes, phenotypic NK-cell differentiation, and prolongation of survival. These studies lay the foundation for epigenetic-directed therapy in ENKTL. SIGNIFICANCE: Through epigenetic and transcriptomic analyses of ENKTL, a rare, aggressive malignancy, along with normal NK-cell developmental intermediates, we identified that extreme DNA hypermethylation targets genes required for NK-cell development. Disrupting this epigenetic blockade in novel PDX models led to ENKTL differentiation and improved survival. This article is highlighted in the In This Issue feature, p. 85.

Clinical outcome and gene- and microRNA-expression profiling according to the Wilms tumor 1 (WT1) single nucleotide polymorphism rs16754 in adult de novo cytogenetically normal acute myeloid leukemia: a Cancer and Leukemia Group B study.

BACKGROUND: The alleles of the Wilms tumor 1 (WT1) polymorphism rs16754 harbor adenine (A) or guanine (G). Recently, rs16754 has been reported to affect the outcome of patients with cytogenetically normal acute myeloid leukemia. To validate this finding, we investigated pretreatment features and outcome associated with rs16754 in a large cohort of patients with cytogenetically normal acute myeloid leukemia. DESIGN AND METHODS: Four-hundred and thirty-three intensively treated and molecularly characterized cytogenetically normal patients with de novo acute myeloid leukemia (18-83 years old) were analyzed for rs16754. To gain biological insights, we studied the gene- and microRNA-expression profiles for associations with rs16754. RESULTS: Three-hundred and nine (71%) patients were homozygous for A (WT1(AA)), 112 (26%) were heterozygous (WT1(AG)) and 12 (3%) were homozygous for G (WT1(GG)). For comparison with previous studies, we grouped WT1(AG) and WT1(GG) patients and compared them with WT1(AA) patients divided into younger (<60 years) and older (≥60 years) adults. We found no independent prognostic impact of WT1(AA). However, WT1(GG) patients, who were less often Caucasian than WT1(AG) (P=0.001) or WT1(AA) (P=0.008) patients, and had TET2 mutations more often than WT1(AG) (P=0.02) patients, had, among patients with FLT3-internal tandem duplication and/or NPM1 wild-type, better disease-free (P=0.02) and overall survival (P=0.04) than WT1(AA) and WT1(AG) patients combined. Unsupervised and supervised analyses of the gene- and microRNA-expression profiles suggested that there were no distinct expression patterns associated with any rs16754 genotype. CONCLUSIONS: We did not observe the previously reported adverse impact of WT1(AA) but found favorable outcomes associated with the homozygous WT1(GG). Considering its low frequency, confirmatory studies are necessary. The biological significance of rs16754 remains questionable as no distinct expression profiles were associated with the genotypes.

DNA methylation changes in a human cell model of breast cancer progression.

Epigenetic inactivation of genes by DNA hypermethylation plays an important role in carcinogenesis. An in vitro model of human breast epithelial cell transformation was used to study epigenetic changes induced by estradiol during the neoplastic process. Different stages of tumor initiation and progression are represented in this model being MCF-10F the normal stage; trMCF cells, the transformed stage; bsMCF cells, the invasive stage and, caMCF cells, the tumor stage. Global methylation studies by restriction landmark genomic scanning (RLGS) showed an increased DNA methylation during the in the invasive and tumor stages. Expression studies showed that NRG1 (neuregulin 1), CSS3 (chondroitin sulfate synthase 3) and SNIP (SNAP-25-interacting protein) were downregulated in the invasive and tumor cells. The transformed cells showed low expression of STXBP6 (amysin) compared to the parental cells MCF-10F. The treatment of these cells with the demethylating agent 5-aza-dC alone or in combination with the histone deacetylase inhibitor trichostatin increased the expression of NRG1, STXBP6, CSS3 and SNIP confirming that DNA methylation plays an important role in the regulation of the expression of these genes. The NRG1 exon 1 has a region located between -136 and +79 (considering +1, the translational initiation site) rich in CpG sites that was analyzed by methylation specific PCR (MSP). NRG1 exon 1 showed progressive changes in the methylation pattern associated with the progression of the neoplastic process in this model; NRG1 exon 1 was unmethylated in MCF-10F and trMCF cells, becoming hypermethylated in the invasive (bsMCF) and tumor (caMCF) stages. Studies of human breast tissue samples showed that NRG1 exon 1 was partially methylated in 14 out of 17 (82.4%) invasive carcinomas although it was unmethylated in normal tissues (8 out of 10 normal breast tissue samples). Furthermore, NRG1 exon 1 was partially methylated in 9 out of 14 (64.3%) morphologically normal tissue samples adjacent to invasive carcinomas.

CpG island methylation in a mouse model of lymphoma is driven by the genetic configuration of tumor cells.

Hypermethylation of CpG islands is a common epigenetic alteration associated with cancer. Global patterns of hypermethylation are tumor-type specific and nonrandom. The biological significance and the underlying mechanisms of tumor-specific aberrant promoter methylation remain unclear, but some evidence suggests that this specificity involves differential sequence susceptibilities, the targeting of DNA methylation activity to specific promoter sequences, or the selection of rare DNA methylation events during disease progression. Using restriction landmark genomic scanning on samples derived from tissue culture and in vivo models of T cell lymphomas, we found that MYC overexpression gave rise to a specific signature of CpG island hypermethylation. This signature reflected gene transcription profiles and was detected only in advanced stages of disease. The further inactivation of the Pten, p53, and E2f2 tumor suppressors in MYC-induced lymphomas resulted in distinct and diagnostic CpG island methylation signatures. Our data suggest that tumor-specific DNA methylation in lymphomas arises as a result of the selection of rare DNA methylation events during the course of tumor development. This selection appears to be driven by the genetic configuration of tumor cells, providing experimental evidence for a causal role of DNA hypermethylation in tumor progression and an explanation for the tremendous epigenetic heterogeneity observed in the evolution of human cancers. The ability to predict genome-wide epigenetic silencing based on relatively few genetic alterations will allow for a more complete classification of tumors and understanding of tumor cell biology.

Diverse histone modifications on histone 3 lysine 9 and their relation to DNA methylation in specifying gene silencing.

BACKGROUND: Previous studies of individual genes have shown that in a self-enforcing way, dimethylation at histone 3 lysine 9 (dimethyl-H3K9) and DNA methylation cooperate to maintain a repressive mode of inactive genes. Less clear is whether this cooperation is generalized in mammalian genomes, such as mouse genome. Here we use epigenomic tools to simultaneously interrogate chromatin modifications and DNA methylation in a mouse leukemia cell line, L1210. RESULTS: Histone modifications on H3K9 and DNA methylation in L1210 were profiled by both global CpG island array and custom mouse promoter array analysis. We used chromatin immunoprecipitation microarray (ChIP-chip) to examine acetyl-H3K9 and dimethyl-H3K9. We found that the relative level of acetyl-H3K9 at different chromatin positions has a wider range of distribution than that of dimethyl-H3K9. We then used differential methylation hybridization (DMH) and the restriction landmark genome scanning (RLGS) to analyze the DNA methylation status of the same targets investigated by ChIP-chip. The results of epigenomic profiling, which have been independently confirmed for individual loci, show an inverse relationship between DNA methylation and histone acetylation in regulating gene silencing. In contrast to the previous notion, dimethyl-H3K9 seems to be less distinct in specifying silencing for the genes tested. CONCLUSION: This study demonstrates in L1210 leukemia cells a diverse relationship between histone modifications and DNA methylation in the maintenance of gene silencing. Acetyl-H3K9 shows an inverse relationship between DNA methylation and histone acetylation in regulating gene silencing as expected. However, dimethyl-H3K9 seems to be less distinct in relation to promoter methylation. Meanwhile, a combination of epigenomic tools is of help in understanding the heterogeneity of epigenetic regulation, which may further our vision accumulated from single-gene studies.

Restriction landmark genomic scanning (RLGS) spot identification by second generation virtual RLGS in multiple genomes with multiple enzyme combinations.

BACKGROUND: Restriction landmark genomic scanning (RLGS) is one of the most successfully applied methods for the identification of aberrant CpG island hypermethylation in cancer, as well as the identification of tissue specific methylation of CpG islands. However, a limitation to the utility of this method has been the ability to assign specific genomic sequences to RLGS spots, a process commonly referred to as "RLGS spot cloning." RESULTS: We report the development of a virtual RLGS method (vRLGS) that allows for RLGS spot identification in any sequenced genome and with any enzyme combination. We report significant improvements in predicting DNA fragment migration patterns by incorporating sequence information into the migration models, and demonstrate a median Euclidian distance between actual and predicted spot migration of 0.18 centimeters for the most complex human RLGS pattern. We report the confirmed identification of 795 human and 530 mouse RLGS spots for the most commonly used enzyme combinations. We also developed a method to filter the virtual spots to reduce the number of extra spots seen on a virtual profile for both the mouse and human genomes. We demonstrate use of this filter to simplify spot cloning and to assist in the identification of spots exhibiting tissue-specific methylation. CONCLUSION: The new vRLGS system reported here is highly robust for the identification of novel RLGS spots. The migration models developed are not specific to the genome being studied or the enzyme combination being used, making this tool broadly applicable. The identification of hundreds of mouse and human RLGS spot loci confirms the strong bias of RLGS studies to focus on CpG islands and provides a valuable resource to rapidly study their methylation.

Microarray analysis of epigenetic silencing of gene expression in the KAS-6/1 multiple myeloma cell line.

The epigenetic control of gene transcription in cancer has been the theme of many recent studies and therapeutic approaches. Carcinogenesis is frequently associated with hypermethylation and consequent down-regulation of genes that prevent cancer, e.g., those that control cell proliferation and apoptosis. We used the demethylating drug zebularine to induce changes in DNA methylation, then examined patterns of gene expression using cDNA array analysis and Restriction Landmark Genomic Scanning followed by RNase protection assay and reverse transcription-PCR to confirm the results. Microarray studies revealed that many genes were epigenetically regulated by methylation. We concluded that methylation decreased the expression of, or silenced, several genes, contributing to the growth and survival of multiple myeloma cells. For example, a number of genes (BAD, BAK, BIK, and BAX) involved in apoptosis were found to be suppressed by methylation. Sequenced methylation-regulated DNA fragments identified by Restriction Landmark Genomic Scanning were found to contain CpG islands, and some corresponded to promoters of genes that were regulated by methylation. We also observed that after the removal of the demethylating drug, the addition of interleukin 6 restored CpG methylation and re-established previously silenced gene patterns, thus implicating a novel role of interleukin 6 in processes regulating epigenetic gene repression and carcinogenesis.