Co-occurrence of cohesin complex and Ras signaling mutations during progression from myelodysplastic syndromes to secondary acute myeloid leukemia.

Myelodysplastic syndromes (MDS) are hematological disorders at high risk of progression to secondary acute myeloid leukemia (sAML). However, the mutational dynamics and clonal evolution underlying disease progression are poorly understood at present. To elucidate the mutational dynamics of pathways and genes occurring during the evolution to sAML, next generation sequencing was performed on 84 serially paired samples of MDS patients who developed sAML (discovery cohort) and 14 paired samples from MDS patients who did not progress to sAML during follow-up (control cohort). Results were validated in an independent series of 388 MDS patients (validation cohort). We used an integrative analysis to identify how mutations, alone or in combination, contribute to leukemic transformation. The study showed that MDS progression to sAML is characterized by greater genomic instability and the presence of several types of mutational dynamics, highlighting increasing (STAG2) and newly-acquired (NRAS and FLT3) mutations. Moreover, we observed cooperation between genes involved in the cohesin and Ras pathways in 15-20% of MDS patients who evolved to sAML, as well as a high proportion of newly acquired or increasing mutations in the chromatin-modifier genes in MDS patients receiving a disease-modifying therapy before their progression to sAML.

Clinical, biological, and prognostic implications of SF3B1 co-occurrence mutations in very low/low- and intermediate-risk MDS patients.

SF3B1 is a highly mutated gene in myelodysplastic syndrome (MDS) patients, related to a specific subtype and parameters of good prognosis in MDS without excess blasts. More than 40% of MDS patients carry at least two myeloid-related gene mutations but little is known about the impact of concurrent mutations on the outcome of MDS patients. In applying next-generation sequencing (NGS) with a 117 myeloid gene custom panel, we analyzed the co-occurrence of SF3B1 with other mutations to reveal their clinical, biological, and prognostic implications in very low/low- and intermediate-risk MDS patients. Mutations in addition to those of SF3B1 were present in 80.4% of patients (median of 2 additional mutations/patient, range 0-5). The most frequently mutated genes were as follows: TET2 (39.2%), DNMT3A (25.5%), SRSF2 (10.8%), CDH23 (5.9%), and ASXL1, CUX1, and KMT2D (4.9% each). The presence of at least two mutations concomitant with that of SF3B1 had an adverse impact on survival compared with those with the SF3B1 mutation and fewer than two additional mutations (median of 54 vs. 87 months, respectively: p = 0.007). The co-occurrence of SF3B1 mutations with specific genes is also linked to a dismal prognosis: SRSF2 mutations were associated with shorter overall survival (OS) than SRSF2wt (median, 27 vs. 75 months, respectively; p = 0.001), concomitant IDH2 mutations (median OS, 11 [mut] vs. 75 [wt] months; p = 0.001), BCOR mutations (median OS, 11 [mut] vs. 71 [wt] months; p = 0.036), and NUP98 and STAG2 mutations (median OS, 27 and 11 vs. 71 months, respectively; p = 0.008 and p = 0.002). Mutations in CHIP genes (TET2, DNMT3A) did not significantly affect the clinical features or outcome. Our results suggest that a more comprehensive NGS study in low-risk MDS SF3B1mut patients is essential for a better prognostic evaluation.

Spanish Guidelines for the use of targeted deep sequencing in myelodysplastic syndromes and chronic myelomonocytic leukaemia.

The landscape of medical sequencing has rapidly changed with the evolution of next generation sequencing (NGS). These technologies have contributed to the molecular characterization of the myelodysplastic syndromes (MDS) and chronic myelomonocytic leukaemia (CMML), through the identification of recurrent gene mutations, which are present in >80% of patients. These mutations contribute to a better classification and risk stratification of the patients. Currently, clinical laboratories include NGS genomic analyses in their routine clinical practice, in an effort to personalize the diagnosis, prognosis and treatment of MDS and CMML. NGS technologies have reduced the cost of large-scale sequencing, but there are additional challenges involving the clinical validation of these technologies, as continuous advances are constantly being made. In this context, it is of major importance to standardize the generation, analysis, clinical interpretation and reporting of NGS data. To that end, the Spanish MDS Group (GESMD) has expanded the present set of guidelines, aiming to establish common quality standards for the adequate implementation of NGS and clinical interpretation of the results, hoping that this effort will ultimately contribute to the benefit of patients with myeloid malignancies.

Genome-wide transcriptomics leads to the identification of deregulated genes after deferasirox therapy in low-risk MDS patients.

The iron chelator deferasirox is widely used in patients with iron overload. Patients with low-grade myelodysplastic syndromes (MDS) get transfusion dependency and need to be treated with deferasirox to avoid iron overload. Moreover, in some patients an increase in both erythroid and platelets have been observed after deferasirox therapy. However, the mechanisms involved in these clinical findings are poorly understood. The aim of this work was to analyze, in patients treated with deferasirox, the changes in the gene-expression profile after receiving the treatment. A total of 15 patients with the diagnosis of low-grade MDS were studied. Microarrays were carried out in RNA from peripheral blood before and after 14 weeks of deferasirox therapy. Changes in 1457 genes and 54 miRNAs were observed: deferasirox induced the downregulation of genes related to the Nf kB pathway leading of an overall inactivation of this pathway. In addition, the iron chelator also downregulated gamma interferon. Altogether these changes could be related to the improvement of erythroid response observed in these patients after therapy. Moreover, the inhibition of NFE2L2/NRF2, which was predicted in silico, could be playing a critical role in the reduction of reactive oxygen species (ROS). Of note, miR-125b, overexpressed after deferasirox treatment, could be involved in the reduced inflammation and increased hematopoiesis observed in the patients after treatment. In summary this study shows, for the first time, the mechanisms that could be governing deferasirox impact in vivo.

Chronic graft-versus-host disease could ameliorate the impact of adverse somatic mutations in patients with myelodysplastic syndromes and hematopoietic stem cell transplantation.

Somatic mutations in patients with myelodysplastic syndromes (MDS) undergoing allogeneic hematopoietic stem cell transplantation (HSTC) are associated with adverse outcome, but the role of chronic graft-versus-host disease (cGVHD) in this subset of patients remains unknown. We analyzed bone marrow samples from 115 patients with MDS collected prior to HSCT using next-generation sequencing. Seventy-one patients (61%) had at least one mutated gene. We found that patients with a higher number of mutated genes (more than 2) had a worse outcome (2 years overall survival [OS] 54.8% vs. 31.1%, p = 0.035). The only two significant variables in the multivariate analysis for OS were TET2 mutations (p = 0.046) and the development of cGVHD, considered as a time-dependent variable (p < 0.001), correlated with a worse and a better outcome, respectively. TP53 mutations also demonstrated impact on the cumulative incidence of relapse (CIR) (1 year CIR 47.1% vs. 9.8%, p = 0.006) and were related with complex karyotype (p = 0.003). cGVHD improved the outcome even among patients with more than 2 mutated genes (1-year OS 88.9% at 1 year vs. 31.3%, p = 0.02) and patients with TP53 mutations (1-year CIR 20% vs. 42.9%, p = 0.553). These results confirm that cGVHD could ameliorate the adverse impact of somatic mutations in patients with MDS with HSCT.

Mutations in TP53 and JAK2 are independent prognostic biomarkers in B-cell precursor acute lymphoblastic leukaemia.

BACKGROUND: In B-cell precursor acute lymphoblastic leukaemia (B-ALL), the identification of additional genetic alterations associated with poor prognosis is still of importance. We determined the frequency and prognostic impact of somatic mutations in children and adult cases with B-ALL treated with Spanish PETHEMA and SEHOP protocols. METHODS: Mutational status of hotspot regions of TP53, JAK2, PAX5, LEF1, CRLF2 and IL7R genes was determined by next-generation deep sequencing in 340 B-ALL patients (211 children and 129 adults). The associations between mutation status and clinicopathological features at the time of diagnosis, treatment outcome and survival were assessed. Univariate and multivariate survival analyses were performed to identify independent prognostic factors associated with overall survival (OS), event-free survival (EFS) and relapse rate (RR). RESULTS: A mutation rate of 12.4% was identified. The frequency of adult mutations was higher (20.2% vs 7.6%, P=0.001). TP53 was the most frequently mutated gene (4.1%), followed by JAK2 (3.8%), CRLF2 (2.9%), PAX5 (2.4%), LEF1 (0.6%) and IL7R (0.3%). All mutations were observed in B-ALL without ETV6-RUNX1 (P=0.047) or BCR-ABL1 fusions (P<0.0001). In children, TP53mut was associated with lower OS (5-year OS: 50% vs 86%, P=0.002) and EFS rates (5-year EFS: 50% vs 78.3%, P=0.009) and higher RR (5-year RR: 33.3% vs 18.6% P=0.037), and was independently associated with higher RR (hazard ratio (HR)=4.5; P=0.04). In adults, TP53mut was associated with a lower OS (5-year OS: 0% vs 43.3%, P=0.019) and a higher RR (5-year RR: 100% vs 61.4%, P=0.029), whereas JAK2mut was associated with a lower EFS (5-year EFS: 0% vs 30.6%, P=0.035) and a higher RR (5-year RR: 100% vs 60.4%, P=0.002). TP53mut was an independent risk factor for shorter OS (HR=2.3; P=0.035) and, together with JAK2mut, also were independent markers of poor prognosis for RR (TP53mut: HR=5.9; P=0.027 and JAK2mut: HR=5.6; P=0.036). CONCLUSIONS: TP53mut and JAK2mut are potential biomarkers associated with poor prognosis in B-ALL patients.

Identification of expression patterns in the progression of disease stages by integration of transcriptomic data.

BACKGROUND: In the study of complex diseases using genome-wide expression data from clinical samples, a difficult case is the identification and mapping of the gene signatures associated to the stages that occur in the progression of a disease. The stages usually correspond to different subtypes or classes of the disease, and the difficulty to identify them often comes from patient heterogeneity and sample variability that can hide the biomedical relevant changes that characterize each stage, making standard differential analysis inadequate or inefficient. RESULTS: We propose a methodology to study diseases or disease stages ordered in a sequential manner (e.g. from early stages with good prognosis to more acute or serious stages associated to poor prognosis). The methodology is applied to diseases that have been studied obtaining genome-wide expression profiling of cohorts of patients at different stages. The approach allows searching for consistent expression patterns along the progression of the disease through two major steps: (i) identifying genes with increasing or decreasing trends in the progression of the disease; (ii) clustering the increasing/decreasing gene expression patterns using an unsupervised approach to reveal whether there are consistent patterns and find genes altered at specific disease stages. The first step is carried out using Gamma rank correlation to identify genes whose expression correlates with a categorical variable that represents the stages of the disease. The second step is done using a Self Organizing Map (SOM) to cluster the genes according to their progressive profiles and identify specific patterns. Both steps are done after normalization of the genomic data to allow the integration of multiple independent datasets. In order to validate the results and evaluate their consistency and biological relevance, the methodology is applied to datasets of three different diseases: myelodysplastic syndrome, colorectal cancer and Alzheimer's disease. A software script written in R, named genediseasePatterns, is provided to allow the use and application of the methodology. CONCLUSION: The method presented allows the analysis of the progression of complex and heterogeneous diseases that can be divided in pathological stages. It identifies gene groups whose expression patterns change along the advance of the disease, and it can be applied to different types of genomic data studying cohorts of patients in different states.

Single nucleotide polymorphism array karyotyping: a diagnostic and prognostic tool in myelodysplastic syndromes with unsuccessful conventional cytogenetic testing.

Cytogenetic aberrations identified by metaphase cytogenetics (MC) have diagnostic, prognostic, and therapeutic implications in myelodysplastic syndromes (MDS). However, in some MDS patients MC study is unsuccesful. Single nucleotide polymorphism array (SNP-A) based karyotyping could be helpful in these cases. We performed SNP-A in 62 samples from bone marrow or peripheral blood of primary MDS with an unsuccessful MC study. SNP-A analysis enabled the detection of aberrations in 31 (50%) patients. We used the copy number alteration information to apply the International Prognostic Scoring System (IPSS) and we observed differences in survival between the low/intermediate-1 and intermediate-2/high risk patients. We also saw differences in survival between very low/low/intermediate and the high/very high patients when we applied the revised IPSS (IPSS-R). In conclusion, SNP-A can be used successfully in PB samples and the identification of CNA by SNP-A improve the diagnostic and prognostic evaluation of this group of MDS patients.

Outcomes and effect of somatic mutations after erythropoiesis stimulating agents in patients with lower-risk myelodysplastic syndromes.

Background: Erythropoiesis stimulating agents (ESAs) are the first-line therapy in patients with lower-risk myelodysplastic syndromes (LR-MDS). Some predictive factors for ESAs response have been identified. Type and number of somatic mutations have been associated with prognosis and response to therapies in MDS patients. Objectives: The objective was to evaluate the outcomes after ESAs in patients with LR-MDS and to address the potential predictive value of somatic mutations in ESAs-treated patients. Design: Multi-center retrospective study of a cohort of 722 patients with LR-MDS included in the SPRESAS (Spanish Registry of Erythropoietic Stimulating Agents Study) study. Retrospective analysis of 65 patients with next generation sequencing (NGS) data from diagnosis. Methods: ESAs' efficacy and safety were evaluated in patients receiving ESAs and best supportive care (BSC). To assess the potential prognostic value of somatic mutations in erythroid response (ER) rate and outcome, NGS was performed in responders and non-responders. Results: ER rate for ESAs-treated patients was 65%. Serum erythropoietin (EPO) level <200 U/l was the only variable significantly associated with a higher ER rate (odds ratio, 2.45; p = 0.036). Median overall survival (OS) in patients treated with ESAs was 6.7 versus 3.1 years in patients receiving BSC (p < 0.001). From 65 patients with NGS data, 57 (87.7%) have at least one mutation. We observed a trend to a higher frequency of ER among patients with a lower number of mutated genes (40.4% in <3 mutated genes versus 22.2% in ⩾3; p = 0.170). The presence of ⩾3 mutated genes was also significantly associated with worse OS (hazard ratio, 2.8; p = 0.015), even in responders. A higher cumulative incidence of acute myeloid leukemia progression at 5 years was also observed in patients with ⩾3 mutated genes versus <3 (33.3% and 10.7%, respectively; p < 0.001). Conclusion: This large study confirms the beneficial effect of ESAs and the adverse effect of somatic mutations in patients with LR-MDS.

Rearrangements involving 11q23.3/KMT2A in adult AML: mutational landscape and prognostic implications - a HARMONY study.

Balanced rearrangements involving the KMT2A gene (KMT2Ar) are recurrent genetic abnormalities in acute myeloid leukemia (AML), but there is lack of consensus regarding the prognostic impact of different fusion partners. Moreover, prognostic implications of gene mutations co-occurring with KMT2Ar are not established. From the HARMONY AML database 205 KMT2Ar adult patients were selected, 185 of whom had mutational information by a panel-based next-generation sequencing analysis. Overall survival (OS) was similar across the different translocations, including t(9;11)(p21.3;q23.3)/KMT2A::MLLT3 (p = 0.756). However, independent prognostic factors for OS in intensively treated patients were age >60 years (HR 2.1, p = 0.001), secondary AML (HR 2.2, p = 0.043), DNMT3A-mut (HR 2.1, p = 0.047) and KRAS-mut (HR 2.0, p = 0.005). In the subset of patients with de novo AML < 60 years, KRAS and TP53 were the prognostically most relevant mutated genes, as patients with a mutation of any of those two genes had a lower complete remission rate (50% vs 86%, p < 0.001) and inferior OS (median 7 vs 30 months, p < 0.001). Allogeneic hematopoietic stem cell transplantation in first complete remission was able to improve OS (p = 0.003). Our study highlights the importance of the mutational patterns in adult KMT2Ar AML and provides new insights into more accurate prognostic stratification of these patients.

Prognostic impact of the number of methylated genes in myelodysplastic syndromes and acute myeloid leukemias treated with azacytidine.

The prognostic impact of the aberrant hypermethylation in response to azacytidine (AZA) remains to be determined. Therefore, we have analyzed the influence of the methylation status prior to AZA treatment on the overall survival and clinical response of myeloid malignancies. DNA methylation status of 24 tumor suppressor genes was analyzed by methylation-specific multiplex ligation-dependent probe amplification in 63 patients with myelodysplastic syndromes and acute myeloid leukemia treated with azacytidine. Most patients (73 %) showed methylation of at least one gene, but only 12 % of patients displayed ≥3 methylated genes. The multivariate analysis demonstrated that the presence of a high number (≥2) of methylated genes (P = 0.022), a high WBC count (P = 0.033), or anemia (P = 0.029) were independent prognostic factors associated with shorter overall survival. The aberrant methylation status did not correlate with the response to AZA, although four of the five patients with ≥3 methylated genes did not respond. By contrast, favorable cytogenetics independently influenced the clinical response to AZA as 64.7 % of patients with good-risk cytogenetic abnormalities responded (P = 0.03). Aberrant methylation status influences the survival of patients treated with AZA, being shorter in those patients with a high number of methylated genes.

Imatinib therapy of chronic myeloid leukemia restores the expression levels of key genes for DNA damage and cell-cycle progression.

BACKGROUND: Chronic myeloid leukemia (CML) is a malignant clonal disorder of the hematopoietic system caused by the expression of the BCR/ABL fusion oncogene. It is well known that CML cells are genetically unstable. However, the mechanisms by which these cells acquire genetic alterations are poorly understood. Imatinib mesylate is the standard therapy for newly diagnosed CML patients. Imatinib mesylate targets the oncogenic kinase activity of BCR-ABL. OBJECTIVE: To study the gene expression profile of bone marrow hematopoietic cells in the same patients with CML before and 1 month after imatinib therapy. METHODS: Samples from patients with CML were analyzed using Affymetrix GeneChip Expression Arrays. RESULTS: A total of 594 differentially expressed genes, most of which (393 genes) were downregulated, as a result of imatinib therapy were observed. CONCLUSION: The blockade of oncoprotein Bcr-Abl by imatinib could cause a decrease in the expression of key DNA repair genes and substantially modify the expression profile of the bone marrow cells in the first days of therapy.

Multicenter Next-Generation Sequencing Studies between Theory and Practice: Harmonization of Data Analysis Using Real-World Myelodysplastic Syndrome Data.

In the age of personalized medicine, genetic testing by means of targeted sequencing has taken a key role. However, when comparing different sets of targeted sequencing data, these are often characterized by a considerable lack of harmonization. Laboratories follow their own best practices, analyzing their own target regions. The question on how to best integrate data from different sites remains unanswered. Studying the example of myelodysplastic syndrome (MDS), we analyzed 11 targeted sequencing sets, collected from six different centers (n = 831). An intersecting target region of 43,076 bp (30 genes) was identified; whereas, the original target regions covered up to 499,097 bp (117 genes). Considering a region of interest in the context of MDS, a target region of 55,969 bp (31 genes) was identified. For each gene, coverage and sequencing data quality was evaluated, calculating a sequencing score. Analyses revealed huge differences between different data sets as well as between different genes. Analysis of the relation between sequencing score and mutation frequency in MDS revealed that most genes with high frequency in MDS could be sequenced without expecting low coverage or quality. Still, no gene appeared consistently unproblematic for all data sets. To allow for comparable results in a multicenter setting analyzing MDS, we propose to use a predefined target region of interest and to perform centralized data analysis using harmonized criteria.

Correction: Mutations in the DNA methylation pathway and number of driver mutations predict response to azacitidine in myelodysplastic syndromes.

[This corrects the article DOI: 10.18632/oncotarget.22157.].

A two-step approach for sequencing spliceosome-related genes as a complementary diagnostic assay in MDS patients with ringed sideroblasts.

Our study aimed to analyze the presence of mutations in SF3B1 and other spliceosome-related genes in myelodysplastic syndromes with ringed sideroblasts (MDS-RS) by combining conventional Sanger and next-generation sequencing (NGS) methods, and to determine the feasibility of this approach in a clinical setting. 122 bone marrow samples from MDS-RS patients were studied. Initially, exons 14 and 15 of the SF3B1 gene were analyzed by Sanger sequencing. Secondly, they were studied by NGS covering besides SF3B1, SRSF2, U2AF1 and ZRSR2 genes. An 86% of all patients showed mutations in the SF3B1 gene. Six of them, which were not identifiable by conventional sequencing in the first diagnostic step, were revealed by NGS. In addition, 19.5% of cases showed mutations in other splicing genes: SRSF2, U2AF1, and ZRSR2. Furthermore, 8.7% of patients had two mutations in SF3B1, SF3B1 and SRSF2, and SF3B1 and U2AF1, while 5.7% showed no mutations in the four spliceosome-related genes analyzed. The combined use of conventional Sanger and NGS allows the identification of mutations in spliceosome-related genes in almost all MDS patients with RS. This two-step approach is affordable and could be useful as a complementary technique in cases with an unclear diagnosis.

Mutations in the DNA methylation pathway and number of driver mutations predict response to azacitidine in myelodysplastic syndromes.

We evaluated the association of mutations in 34 candidate genes and response to azacitidine in 84 patients with myelodysplastic syndrome (MDS), with 217 somatic mutations identified by next-generation sequencing. Most patients (93%) had ≥1 mutation (mean=2.6/patient). The overall response rate to azacitidine was 42%. No clinical characteristic was associated with response to azacitidine. However, total number of mutations/patient was negatively associated with overall drug response (odds ratio [OR]: 0.56, 95% confidence interval [CI]: 0.33-0.94; p=0.028), and a positive association was found for having ≥1 mutation in a DNA methylation-related gene: TET2, DNMT3A, IDH1 and/or IDH2 (OR: 4.76, 95%CI: 1.31-17.27; p=0.017). Mutations in TP53 (hazard ratio [HR]: 3.88; 95%CI: 1.94-7.75) and EZH2 (HR: 2.50; 95%CI: 1.23-5.09) were associated with shorter overall survival. Meta-analysis of 6 studies plus present data (n=815 patients) allowed assessment of the association of drug response with mutations in 9 candidate genes: ASXL1, CBL, EZH2, SF3B1, SRSF2, TET2, DNMT3A, IDH1/2 and TP53. TET2 mutations predicted a more favorable drug response compared with 'wild-type' peers (pooled OR: 1.67, 95%CI: 1.14-2.44; p=0.01). In conclusion, mutations in the DNA methylation pathway, especially TET2 mutations, and low number of total mutations are associated with a better response to azacitidine.

Genome-Wide DNA Copy Number Analysis of Acute Lymphoblastic Leukemia Identifies New Genetic Markers Associated with Clinical Outcome.

UNLABELLED: Identifying additional genetic alterations associated with poor prognosis in acute lymphoblastic leukemia (ALL) is still a challenge. AIMS: To characterize the presence of additional DNA copy number alterations (CNAs) in children and adults with ALL by whole-genome oligonucleotide array (aCGH) analysis, and to identify their associations with clinical features and outcome. Array-CGH was carried out in 265 newly diagnosed ALLs (142 children and 123 adults). The NimbleGen CGH 12x135K array (Roche) was used to analyze genetic gains and losses. CNAs were analyzed with GISTIC and aCGHweb software. Clinical and biological variables were analyzed. Three of the patients showed chromothripsis (cth6, cth14q and cth15q). CNAs were associated with age, phenotype, genetic subtype and overall survival (OS). In the whole cohort of children, the losses on 14q32.33 (p = 0.019) and 15q13.2 (p = 0.04) were related to shorter OS. In the group of children without good- or poor-risk cytogenetics, the gain on 1p36.11 was a prognostic marker independently associated with shorter OS. In adults, the gains on 19q13.2 (p = 0.001) and Xp21.1 (p = 0.029), and the loss of 17p (p = 0.014) were independent markers of poor prognosis with respect to OS. In summary, CNAs are frequent in ALL and are associated with clinical parameters and survival. Genome-wide DNA copy number analysis allows the identification of genetic markers that predict clinical outcome, suggesting that detection of these genetic lesions will be useful in the management of patients newly diagnosed with ALL.

Chromothripsis Is a Recurrent Genomic Abnormality in High-Risk Myelodysplastic Syndromes.

To explore novel genetic abnormalities occurring in myelodysplastic syndromes (MDS) through an integrative study combining array-based comparative genomic hybridization (aCGH) and next-generation sequencing (NGS) in a series of MDS and MDS/myeloproliferative neoplasms (MPN) patients. 301 patients diagnosed with MDS (n = 240) or MDS/MPN (n = 61) were studied at the time of diagnosis. A genome-wide analysis of DNA copy number abnormalities was performed. In addition, a mutational analysis of DNMT3A, TET2, RUNX1, TP53 and BCOR genes was performed by NGS in selected cases. 285 abnormalities were identified in 71 patients (23.6%). Three high-risk MDS cases (1.2%) displayed chromothripsis involving exclusively chromosome 13 and affecting some cancer genes: FLT3, BRCA2 and RB1. All three cases carried TP53 mutations as revealed by NGS. Moreover, in the whole series, the integrative analysis of aCGH and NGS enabled the identification of cryptic recurrent deletions in 2p23.3 (DNMT3A; n = 2.8%), 4q24 (TET2; n = 10%) 17p13 (TP53; n = 8.5%), 21q22 (RUNX1; n = 7%), and Xp11.4 (BCOR; n = 2.8%), while mutations in the non-deleted allele where found only in DNMT3A (n = 1), TET2 (n = 3), and TP53 (n = 4). These cryptic abnormalities were detected mainly in patients with normal (45%) or non-informative (15%) karyotype by conventional cytogenetics, except for those with TP53 deletion and mutation (15%), which had a complex karyotype. In addition to well-known copy number defects, the presence of chromothripsis involving chromosome 13 was a novel recurrent change in high-risk MDS patients. Array CGH analysis revealed the presence of cryptic abnormalities in genomic regions where MDS-related genes, such as TET2, DNMT3A, RUNX1 and BCOR, are located.

aCGH-MAS: analysis of aCGH by means of multiagent system.

There are currently different techniques, such as CGH arrays, to study genetic variations in patients. CGH arrays analyze gains and losses in different regions in the chromosome. Regions with gains or losses in pathologies are important for selecting relevant genes or CNVs (copy-number variations) associated with the variations detected within chromosomes. Information corresponding to mutations, genes, proteins, variations, CNVs, and diseases can be found in different databases and it would be of interest to incorporate information of different sources to extract relevant information. This work proposes a multiagent system to manage the information of aCGH arrays, with the aim of providing an intuitive and extensible system to analyze and interpret the results. The agent roles integrate statistical techniques to select relevant variations and visualization techniques for the interpretation of the final results and to extract relevant information from different sources of information by applying a CBR system.