Somatic nuclear mitochondrial DNA insertions are prevalent in the human brain and accumulate over time in fibroblasts.

The transfer of mitochondrial DNA into the nuclear genomes of eukaryotes (Numts) has been linked to lifespan in nonhuman species and recently demonstrated to occur in rare instances from one human generation to the next. Here, we investigated numtogenesis dynamics in humans in 2 ways. First, we quantified Numts in 1,187 postmortem brain and blood samples from different individuals. Compared to circulating immune cells (n = 389), postmitotic brain tissue (n = 798) contained more Numts, consistent with their potential somatic accumulation. Within brain samples, we observed a 5.5-fold enrichment of somatic Numt insertions in the dorsolateral prefrontal cortex (DLPFC) compared to cerebellum samples, suggesting that brain Numts arose spontaneously during development or across the lifespan. Moreover, an increase in the number of brain Numts was linked to earlier mortality. The brains of individuals with no cognitive impairment (NCI) who died at younger ages carried approximately 2 more Numts per decade of life lost than those who lived longer. Second, we tested the dynamic transfer of Numts using a repeated-measures whole-genome sequencing design in a human fibroblast model that recapitulates several molecular hallmarks of aging. These longitudinal experiments revealed a gradual accumulation of 1 Numt every ~13 days. Numtogenesis was independent of large-scale genomic instability and unlikely driven by cell clonality. Targeted pharmacological perturbations including chronic glucocorticoid signaling or impairing mitochondrial oxidative phosphorylation (OxPhos) only modestly increased the rate of numtogenesis, whereas patient-derived SURF1-mutant cells exhibiting mtDNA instability accumulated Numts 4.7-fold faster than healthy donors. Combined, our data document spontaneous numtogenesis in human cells and demonstrate an association between brain cortical somatic Numts and human lifespan. These findings open the possibility that mito-nuclear horizontal gene transfer among human postmitotic tissues produces functionally relevant human Numts over timescales shorter than previously assumed.

Cortical proteins may provide motor resilience in older adults.

Motor resilience proteins may be a high value therapeutic target that offset the negative effects of pathologies on motor function. This study sought to identify cortical proteins associated with motor decline unexplained by brain pathologies that provide motor resilience. We studied 1226 older decedents with annual motor testing, postmortem brain pathologies and quantified 226 proteotypic peptides in prefrontal cortex. Twenty peptides remained associated with motor decline in models controlling for ten brain pathologies (FDR < 0.05). Higher levels of nine peptides and lower levels of eleven peptides were related to slower decline. A higher motor resilience protein score based on averaging the levels of all 20 peptides was related to slower motor decline, less severe parkinsonism and lower odds of mobility disability before death. Cortical proteins may provide motor resilience. Targeting these proteins in further drug discovery may yield novel interventions to maintain motor function in old age.

BIN1 protein isoforms are differentially expressed in astrocytes, neurons, and microglia: neuronal and astrocyte BIN1 are implicated in tau pathology.

BACKGROUND: Identified as an Alzheimer's disease (AD) susceptibility gene by genome wide-association studies, BIN1 has 10 isoforms that are expressed in the Central Nervous System (CNS). The distribution of these isoforms in different cell types, as well as their role in AD pathology still remains unclear. METHODS: Utilizing antibodies targeting specific BIN1 epitopes in human post-mortem tissue and analyzing mRNA expression data from purified microglia, we identified three isoforms expressed in neurons and astrocytes (isoforms 1, 2 and 3) and four isoforms expressed in microglia (isoforms 6, 9, 10 and 12). The abundance of selected peptides, which correspond to groups of BIN1 protein isoforms, was measured in dorsolateral prefrontal cortex, and their relation to neuropathological features of AD was assessed. RESULTS: Peptides contained in exon 7 of BIN1's N-BAR domain were found to be significantly associated with AD-related traits and, particularly, tau tangles. Decreased expression of BIN1 isoforms containing exon 7 is associated with greater accumulation of tangles and subsequent cognitive decline, with astrocytic rather than neuronal BIN1 being the more likely culprit. These effects are independent of the BIN1 AD risk variant. CONCLUSIONS: Exploring the molecular mechanisms of specific BIN1 isoforms expressed by astrocytes may open new avenues for modulating the accumulation of Tau pathology in AD.

AML1-ETO requires enhanced C/D box snoRNA/RNP formation to induce self-renewal and leukaemia.

Leukaemogenesis requires enhanced self-renewal, which is induced by oncogenes. The underlying molecular mechanisms remain incompletely understood. Here, we identified C/D box snoRNAs and rRNA 2'-O-methylation as critical determinants of leukaemic stem cell activity. Leukaemogenesis by AML1-ETO required expression of the groucho-related amino-terminal enhancer of split (AES). AES functioned by inducing snoRNA/RNP formation via interaction with the RNA helicase DDX21. Similarly, global loss of C/D box snoRNAs with concomitant loss of rRNA 2'-O-methylation resulted in decreased leukaemia self-renewal potential. Genomic deletion of either C/D box snoRNA SNORD14D or SNORD35A suppressed clonogenic potential of leukaemia cells in vitro and delayed leukaemogenesis in vivo. We further showed that AML1-ETO9a, MYC and MLL-AF9 all enhanced snoRNA formation. Expression levels of C/D box snoRNAs in AML patients correlated closely with in vivo frequency of leukaemic stem cells. Collectively, these findings indicate that induction of C/D box snoRNA/RNP function constitutes an important pathway in leukaemogenesis.

Integrative analysis of multiple genomic variables using a hierarchical Bayesian model.

MOTIVATION: Genes showing congruent differences in several genomic variables between two biological conditions are crucial to unravel causalities behind phenotypes of interest. Detecting such genes is important in biomedical research, e.g. when identifying genes responsible for cancer development. Small sample sizes common in next-generation sequencing studies are a key challenge, and there are still only very few statistical methods to analyze more than two genomic variables in an integrative, model-based way. Here, we present a novel bioinformatics approach to detect congruent differences between two biological conditions in a larger number of different measurements such as various epigenetic marks or mRNA transcript levels. RESULTS: We propose a coefficient quantifying the degree to which genes present consistent alterations in multiple (more than two) genomic variables when comparing samples presenting a condition of interest (e.g. cancer) to a reference group. A hierarchical Bayesian model is employed to assess uncertainty on a gene level, incorporating information on functional relationships between genes. We demonstrate the approach on different data sets containing RNA-seq gene transcripton and up to four ChIP-seq histone modification measurements. Both the coefficient-based ranking and the inference based on the model lead to a plausible prioritizing of candidate genes when analyzing multiple genomic variables. AVAILABILITY AND IMPLEMENTATION: BUGS code in the Supplement. CONTACT: m.schaefer@uni-duesseldorf.de. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Loss of the histone methyltransferase EZH2 induces resistance to multiple drugs in acute myeloid leukemia.

In acute myeloid leukemia (AML), therapy resistance frequently occurs, leading to high mortality among patients. However, the mechanisms that render leukemic cells drug resistant remain largely undefined. Here, we identified loss of the histone methyltransferase EZH2 and subsequent reduction of histone H3K27 trimethylation as a novel pathway of acquired resistance to tyrosine kinase inhibitors (TKIs) and cytotoxic drugs in AML. Low EZH2 protein levels correlated with poor prognosis in AML patients. Suppression of EZH2 protein expression induced chemoresistance of AML cell lines and primary cells in vitro and in vivo. Low EZH2 levels resulted in derepression of HOX genes, and knockdown of HOXB7 and HOXA9 in the resistant cells was sufficient to improve sensitivity to TKIs and cytotoxic drugs. The endogenous loss of EZH2 expression in resistant cells and primary blasts from a subset of relapsed AML patients resulted from enhanced CDK1-dependent phosphorylation of EZH2 at Thr487. This interaction was stabilized by heat shock protein 90 (HSP90) and followed by proteasomal degradation of EZH2 in drug-resistant cells. Accordingly, inhibitors of HSP90, CDK1 and the proteasome prevented EZH2 degradation, decreased HOX gene expression and restored drug sensitivity. Finally, patients with reduced EZH2 levels at progression to standard therapy responded to the combination of bortezomib and cytarabine, concomitant with the re-establishment of EZH2 expression and blast clearance. These data suggest restoration of EZH2 protein as a viable approach to overcome treatment resistance in this AML patient population.

Genomic landscape of liposarcoma.

Liposarcoma (LPS) is the most common type of soft tissue sarcoma accounting for 20% of all adult sarcomas. Due to absence of clinically effective treatment options in inoperable situations and resistance to chemotherapeutics, a critical need exists to identify novel therapeutic targets. We analyzed LPS genomic landscape using SNP arrays, whole exome sequencing and targeted exome sequencing to uncover the genomic information for development of specific anti-cancer targets. SNP array analysis indicated known amplified genes (MDM2, CDK4, HMGA2) and important novel genes (UAP1, MIR557, LAMA4, CPM, IGF2, ERBB3, IGF1R). Carboxypeptidase M (CPM), recurrently amplified gene in well-differentiated/de-differentiated LPS was noted as a putative oncogene involved in the EGFR pathway. Notable deletions were found at chromosome 1p (RUNX3, ARID1A), chromosome 11q (ATM, CHEK1) and chromosome 13q14.2 (MIR15A, MIR16-1). Significantly and recurrently mutated genes (false discovery rate < 0.05) included PLEC (27%), MXRA5 (21%), FAT3 (24%), NF1 (20%), MDC1 (10%), TP53 (7%) and CHEK2 (6%). Further, in vitro and in vivo functional studies provided evidence for the tumor suppressor role for Neurofibromin 1 (NF1) gene in different subtypes of LPS. Pathway analysis of recurrent mutations demonstrated signaling through MAPK, JAK-STAT, Wnt, ErbB, axon guidance, apoptosis, DNA damage repair and cell cycle pathways were involved in liposarcomagenesis. Interestingly, we also found mutational and copy number heterogeneity within a primary LPS tumor signifying the importance of multi-region sequencing for cancer-genome guided therapy. In summary, these findings provide insight into the genomic complexity of LPS and highlight potential druggable pathways for targeted therapeutic approach.

Deep Sequencing in Conjunction with Expression and Functional Analyses Reveals Activation of FGFR1 in Ewing Sarcoma.

PURPOSE: A low mutation rate seems to be a general feature of pediatric cancers, in particular in oncofusion gene-driven tumors. Genetically, Ewing sarcoma is defined by balanced chromosomal EWS/ETS translocations, which give rise to oncogenic chimeric proteins (EWS-ETS). Other contributing somatic mutations involved in disease development have only been observed at low frequency. EXPERIMENTAL DESIGN: Tumor samples of 116 Ewing sarcoma patients were analyzed here. Whole-genome sequencing was performed on two patients with normal, primary, and relapsed tissue. Whole-exome sequencing was performed on 50 Ewing sarcoma and 22 matched normal tissues. A discovery dataset of 14 of these tumor/normal pairs identified 232 somatic mutations. Recurrent nonsynonymous mutations were validated in the 36 remaining exomes. Transcriptome analysis was performed in a subset of 14 of 50 Ewing sarcomas and DNA copy number gain and expression of FGFR1 in 63 of 116 Ewing sarcomas. RESULTS: Relapsed tumors consistently showed a 2- to 3-fold increased number of mutations. We identified several recurrently mutated genes at low frequency (ANKRD30A, CCDC19, KIAA0319, KIAA1522, LAMB4, SLFN11, STAG2, TP53, UNC80, ZNF98). An oncogenic fibroblast growth factor receptor 1 (FGFR1) mutation (N546K) was detected, and the FGFR1 locus frequently showed copy number gain (31.7%) in primary tumors. Furthermore, high-level FGFR1 expression was noted as a characteristic feature of Ewing sarcoma. RNA interference of FGFR1 expression in Ewing sarcoma lines blocked proliferation and completely suppressed xenograft tumor growth. FGFR1 tyrosine kinase inhibitor (TKI) therapy in a patient with Ewing sarcoma relapse significantly reduced 18-FDG-PET activity. CONCLUSIONS: FGFR1 may constitute a promising target for novel therapeutic approaches in Ewing sarcoma.

BRCC3 mutations in myeloid neoplasms.

Next generation sequencing technologies have provided insights into the molecular heterogeneity of various myeloid neoplasms, revealing previously unknown somatic genetic events. In our cohort of 1444 cases analyzed by next generation sequencing, somatic mutations in the gene BRCA1-BRCA2-containing complex 3 (BRCC3) were identified in 28 cases (1.9%). BRCC3 is a member of the JAMM/MPN+ family of zinc metalloproteases capable of cleaving Lys-63 linked polyubiquitin chains, and is implicated in DNA repair. The mutations were located throughout its coding region. The average variant allelic frequency of BRCC3 mutations was 30.1%, and by a serial sample analysis at two different time points a BRCC3 mutation was already identified in the initial stage of a myelodysplastic syndrome. BRCC3 mutations commonly occurred in nonsense (n=12), frameshift (n=4), and splice site (n=5) configurations. Due to the marginal male dominance (odds ratio; 2.00, 0.84-4.73) of BRCC3 mutations, the majority of mutations (n=23; 82%) were hemizygous. Phenotypically, BRCC3 mutations were frequently observed in myelodysplastic syndromes and myelodysplastic/myeloproliferative neoplasms and associated with -Y abnormality (odds ratio; 3.70, 1.25-11.0). Clinically, BRCC3 mutations were also related to higher age (P=0.01), although prognosis was not affected. Knockdown of Brcc3 gene expression in murine bone marrow lineage negative, Sca1 positive, c-kit positive cells resulted in 2-fold more colony formation and modest differentiation defect. Thus, BRCC3 likely plays a role as tumor-associated gene in myelodysplastic syndromes and myelodysplastic/myeloproliferative neoplasms.

Inherited and Somatic Defects in DDX41 in Myeloid Neoplasms.

Most cases of adult myeloid neoplasms are routinely assumed to be sporadic. Here, we describe an adult familial acute myeloid leukemia (AML) syndrome caused by germline mutations in the DEAD/H-box helicase gene DDX41. DDX41 was also found to be affected by somatic mutations in sporadic cases of myeloid neoplasms as well as in a biallelic fashion in 50% of patients with germline DDX41 mutations. Moreover, corresponding deletions on 5q35.3 present in 6% of cases led to haploinsufficient DDX41 expression. DDX41 lesions caused altered pre-mRNA splicing and RNA processing. DDX41 is exemplary of other RNA helicase genes also affected by somatic mutations, suggesting that they constitute a family of tumor suppressor genes.

Epigenetic dysregulation of KCa 3.1 channels induces poor prognosis in lung cancer.

Epigenomic changes are an important feature of malignant tumors. How tumor aggressiveness is affected by DNA methylation of specific loci is largely unexplored. In genome-wide DNA methylation analyses, we identified the KCa 3.1 channel gene (KCNN4) promoter to be hypomethylated in an aggressive non-small-cell lung carcinoma (NSCLC) cell line and in patient samples. Accordingly, KCa 3.1 expression was increased in more aggressive NSCLC cells. Both findings were strong predictors for poor prognosis in lung adenocarcinoma. Increased KCa 3.1 expression was associated with aggressive features of NSCLC cells. Proliferation and migration of pro-metastatic NSCLC cells depended on KCa 3.1 activity. Mechanistically, elevated KCa 3.1 expression hyperpolarized the membrane potential, thereby augmenting the driving force for Ca(2+) influx. KCa 3.1 blockade strongly reduced the growth of xenografted NSCLC cells in mice as measured by positron emission tomography-computed tomography. Thus, loss of DNA methylation of the KCNN4 promoter and increased KCa 3.1 channel expression and function are mechanistically linked to poor survival of NSCLC patients.

A new workflow for whole-genome sequencing of single human cells.

Unbiased amplification of the whole-genome amplification (WGA) of single cells is crucial to study cancer evolution and genetic heterogeneity, but is challenging due to the high complexity of the human genome. Here, we present a new workflow combining an efficient adapter-linker PCR-based WGA method with second-generation sequencing. This approach allows comparison of single cells at base pair resolution. Amplification recovered up to 74% of the human genome. Copy-number variants and loss of heterozygosity detected in single cell genomes showed concordance of up to 99% to pooled genomic DNA. Allele frequencies of mutations could be determined accurately due to an allele dropout rate of only 2%, clearly demonstrating the low bias of our PCR-based WGA approach. Sequencing with paired-end reads allowed genome-wide analysis of structural variants. By direct comparison to other WGA methods, we further endorse its suitability to analyze genetic heterogeneity.

DNA methyltransferase inhibition reverses epigenetically embedded phenotypes in lung cancer preferentially affecting polycomb target genes.

PURPOSE: Cancer cell phenotypes are partially determined by epigenetic specifications, such as DNA methylation. Metastasis development is a late event in cancerogenesis and might be associated with epigenetic alterations. EXPERIMENTAL DESIGN: An in vivo selection approach was used to generate highly aggressive non-small cell lung cancer (NSCLC) cell lines (A549 and HTB56) followed by genome-wide DNA methylation analysis. Furthermore, the therapeutic effects of the epigenetic agent azacytidine on DNA methylation patterns and the in vivo phenotypes were explored. RESULTS: Widespread changes of DNA methylation were observed during development of highly aggressive cell lines. Up to 2.5% of the CpG-rich region was differentially methylated as identified by reduced representation bisulfite sequencing compared with the less aggressive parental cell lines. DNA methyltransferase inhibition by azacytidine reversed the prometastatic phenotype; this was highly associated with the preferential loss of DNA methylation at sites that were hypermethylated during the in vivo selection. Of note, polycomb (PRC2) binding sites were particularly affected by DNA methylation changes after azacytidine exposure that persisted over time. CONCLUSIONS: We could show that metastatic capability of NSCLC is closely associated with DNA methylome alterations. Because inhibition of DNA methyltransferase reversed metastasis-prone phenotype, epigenetic modulation seems to be a potential therapeutic approach to prevent metastasis formation.

Endoplasmic reticulum protein GliPR1 regulates G protein signaling and the cell cycle and is overexpressed in AML.

Glioma pathogenesis­related protein 1 (GliPR1) is a pleiotropic protein involved in cell proliferation, tumor growth and apoptosis. The aim of the present study was to further characterize GliPR1 in regard to its subcellular localization and its overall effect on cellular gene expression. Knockdown of GliPR1 and Affymetrix microarray mRNA expression analysis revealed 262 GliPR1­dependent differentially expressed genes, of which 40 were induced and 222 were suppressed. Differentially expressed genes were overrepresented in five Gene Ontology categories: G protein signaling pathways, regulation of cyclin­dependent protein kinase activity, ER to Golgi vesicle-mediated transport, axon guidance and dephosphorylation. GliPR1-EGFP fusion protein co­localized with the endoplasmic reticulum (ER) or with cytoplasmic vesicles as demonstrated by confocal microscopy. GliPR1 expression was found to be significantly increased in acute myeloid leukemia (AML) bone marrow samples, while markedly reduced in acute lymphoblastic leukemia, unchanged in myelodysplastic syndrome and slightly decreased in chronic lymphocytic leukemia as well as in chronic myelocytic leukemia (CML) when compared to normal samples. GliPR1 was localized and involved in the ER secretory protein pathway. GliPR1 affects G protein signaling and cell cycle regulation. Based on the observed overexpression in AML samples, GliPR1 should be further explored as a potential target for AML.

Detection of significantly differentially methylated regions in targeted bisulfite sequencing data.

MOTIVATION: Bisulfite sequencing is currently the gold standard to obtain genome-wide DNA methylation profiles in eukaryotes. In contrast to the rapid development of appropriate pre-processing and alignment software, methods for analyzing the resulting methylation profiles are relatively limited so far. For instance, an appropriate pipeline to detect DNA methylation differences between cancer and control samples is still required. RESULTS: We propose an algorithm that detects significantly differentially methylated regions in data obtained by targeted bisulfite sequencing approaches, such as reduced representation bisulfite sequencing. In a first step, this approach tests all target regions for methylation differences by taking spatial dependence into account. A false discovery rate procedure controls the expected proportion of incorrectly rejected regions. In a second step, the significant target regions are trimmed to the actually differentially methylated regions. This hierarchical procedure detects differentially methylated regions with increased power compared with existing methods. AVAILABILITY: R/Bioconductor package BiSeq. SUPPLEMENTARY INFORMATION: Supplementary Data are available at Bioinformatics online.

Robustness of amplicon deep sequencing underlines its utility in clinical applications.

We investigated the robustness of amplicon deep sequencing to study its utility in routine clinical applications offering patient-specific individualized assays for molecular disease characterization and monitoring. Amplicons were designed targeting RUNX1, CEBPA, CBL, NRAS, KRAS, DNMT3A, EZH2, and TP53 using different PCR amplification strategies and Roche GS FLX Titanium and Illumina MiSeq sequencing platforms. Thirty-three patients with leukemia were selected as an exemplary cohort representing heterogeneous cancer specimens. Both standard two-primer amplification and four-primer microfluidics PCRs yielded highly linear characteristics in detecting molecular alterations in series of dilution experiments. By fitting a linear mixed-effects model to the logarithmized data, a slope ß of -1.000 (95% CI, ±0.046) was obtained for two-primer assays and of -0.998 (95% CI, ±0.105) was obtained for four-primer assays, which represented a near-perfect decrease of the mutation load. Furthermore, data are presented on technical precision, limit of detection, and occurrence of small subclones in TP53- and RUNX1-mutated patients to identify clonal disease progression and residual disease. We demonstrate that, depending on the local sequence context for each amplicon, the limit of detection of the assay cannot be lower than a range of 0.25% to 3.5%. In conclusion, amplicon deep sequencing enabled the assessment of mutations in a highly robust manner and across a broad range of relative frequencies of mutations. This assay detects residual disease or identifies mutations with predictive relevance to direct treatment strategies.

Genomic instability may originate from imatinib-refractory chronic myeloid leukemia stem cells.

Genomic instability is a hallmark of chronic myeloid leukemia in chronic phase (CML-CP) resulting in BCR-ABL1 mutations encoding resistance to tyrosine kinase inhibitors (TKIs) and/or additional chromosomal aberrations leading to disease relapse and/or malignant progression. TKI-naive and TKI-treated leukemia stem cells (LSCs) and leukemia progenitor cells (LPCs) accumulate high levels of reactive oxygen species (ROS) and oxidative DNA damage. To determine the role of TKI-refractory LSCs in genomic instability, we used a murine model of CML-CP where ROS-induced oxidative DNA damage was elevated in LSCs, including quiescent LSCs, but not in LPCs. ROS-induced oxidative DNA damage in LSCs caused clinically relevant genomic instability in CML-CP-like mice, such as TKI-resistant BCR-ABL1 mutations (E255K, T315I, H396P), deletions in Ikzf1 and Trp53, and additions in Zfp423 and Idh1. Despite inhibition of BCR-ABL1 kinase, imatinib did not downregulate ROS and oxidative DNA damage in TKI-refractory LSCs to the levels detected in normal cells, and CML-CP-like mice treated with imatinib continued to accumulate clinically relevant genetic aberrations. Inhibition of class I p21-activated protein kinases by IPA3 downregulated ROS in TKI-naive and TKI-treated LSCs. Altogether, we postulate that genomic instability may originate in the most primitive TKI-refractory LSCs in TKI-naive and TKI-treated patients.

Characterization of ZC3H15 as a potential TRAF-2-interacting protein implicated in the NFκB pathway and overexpressed in AML.

The gene product of the zinc finger CCCH-type containing 15 (ZC3H15) gene, an immediate early erythropoietin response gene (synonymous: LEREPO4), was further characterized. ZC3H15 was expressed ubiquitously in all human tissues tested by northern blotting and showed mainly a diffuse cytoplasmic distribution by immune fluorescence microscopy and western blotting of subcellular protein fractions. The expression of ZC3H15 was downregulated effectively in HeLa cells to ≤13% of the control by transfection of specific small interfering RNA (siRNA). Subsequent Affymetrix microarray analysis revealed 202 differentially expressed genes including 114 induced (≥3-fold) genes and 88 suppressed (≤0.3-fold) genes. The gene ontology (GO) categories containing an over-representation of differentially expressed genes comprised cell growth, transcription, cell adhesion, regulation of NF-κB, regulation of MAPK, cell cycle arrest and immune response. ZC3H15 interacted with the signaling adapter protein tumor necrosis factor receptor associated factor 2 (TRAF-2) as shown by co-immunoprecipitation. ZC3H15 expression was found to be significantly increased in acute myeloid leukemia (AML) samples compared to MDS, CML, ALL and normal bone marrow samples using the Leukemia Gene Atlas (LGA) database. Based on these data, it is hypothesized that ZC3H15 may interact with TRAF-2 functionally within the NF-κB pathway, and may be explored as a potential target in AML.

Cancer gene prioritization by integrative analysis of mRNA expression and DNA copy number data: a comparative review.

A variety of genome-wide profiling techniques are available to investigate complementary aspects of genome structure and function. Integrative analysis of heterogeneous data sources can reveal higher level interactions that cannot be detected based on individual observations. A standard integration task in cancer studies is to identify altered genomic regions that induce changes in the expression of the associated genes based on joint analysis of genome-wide gene expression and copy number profiling measurements. In this review, we highlight common approaches to genomic data integration and provide a transparent benchmarking procedure to quantitatively compare method performances in cancer gene prioritization. Algorithms, data sets and benchmarking results are available at http://intcomp.r-forge.r-project.org.

DNA methylation changes are a late event in acute promyelocytic leukemia and coincide with loss of transcription factor binding.

The origin of aberrant DNA methylation in cancer remains largely unknown. In the present study, we elucidated the DNA methylome in primary acute promyelocytic leukemia (APL) and the role of promyelocytic leukemia-retinoic acid receptor α (PML-RARα) in establishing these patterns. Cells from APL patients showed increased genome-wide DNA methylation with higher variability than healthy CD34(+) cells, promyelocytes, and remission BM cells. A core set of differentially methylated regions in APL was identified. Age at diagnosis, Sanz score, and Flt3-mutation status characterized methylation subtypes. Transcription factor-binding sites (eg, the c-myc-binding sites) were associated with low methylation. However, SUZ12- and REST-binding sites identified in embryonic stem cells were preferentially DNA hypermethylated in APL cells. Unexpectedly, PML-RARα-binding sites were also protected from aberrant DNA methylation in APL cells. Consistent with this, myeloid cells from preleukemic PML-RARα knock-in mice did not show altered DNA methylation and the expression of PML-RARα in hematopoietic progenitor cells prevented differentiation without affecting DNA methylation. Treatment of APL blasts with all-trans retinoic acid also did not result in immediate DNA methylation changes. The results of the present study suggest that aberrant DNA methylation is associated with leukemia phenotype but is not required for PML-RARα-mediated initiation of leukemogenesis.