Email updates

Keep up to date with the latest news and content from Molecular Cancer and BioMed Central.

Open Access Research

Non random distribution of genomic features in breakpoint regions involved in chronic myeloid leukemia cases with variant t(9;22) or additional chromosomal rearrangements

Francesco Albano1*, Luisa Anelli1, Antonella Zagaria1, Nicoletta Coccaro1, Paola Casieri1, Antonella Russo Rossi1, Laura Vicari2, Vincenzo Liso1, Mariano Rocchi3 and Giorgina Specchia1

Author Affiliations

1 Hematology, University of Bari, 70124, Bari, Italy

2 Service of Medical Genetics, Cardarelli Hospital, via Cardarelli 9, 80131 Naples, Italy

3 Department of Genetics and Microbiology, University of Bari, 70126 Bari, Italy

For all author emails, please log on.

Molecular Cancer 2010, 9:120  doi:10.1186/1476-4598-9-120


The electronic version of this article is the complete one and can be found online at: http://www.molecular-cancer.com/content/9/1/120


Received:19 January 2010
Accepted:25 May 2010
Published:25 May 2010

© 2010 Albano et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The t(9;22)(q34;q11), generating the Philadelphia (Ph) chromosome, is found in more than 90% of patients with chronic myeloid leukemia (CML). As a result of the translocation, the 3' portion of the ABL1 oncogene is transposed from 9q34 to the 5' portion of the BCR gene on chromosome 22 to form the BCR/ABL1 fusion gene. At diagnosis, in 5-10% of CML patients the Ph chromosome is derived from variant translocations other than the standard t(9;22).

Results

We report a molecular cytogenetic study of 452 consecutive CML patients at diagnosis, that revealed 50 cases identifying three main subgroups: i) cases with variant chromosomal rearrangements other than the classic t(9;22)(q34;q11) (9.5%); ii) cases with cryptic insertions of ABL1 into BCR, or vice versa (1.3%); iii) cases bearing additional chromosomal rearrangements concomitant to the t(9;22) (1.1%). For each cytogenetic group, the mechanism at the basis of the rearrangement is discussed.

All breakpoints on other chromosomes involved in variant t(9;22) and in additional rearrangements have been characterized for the first time by Fluorescence In Situ Hybridization (FISH) experiments and bioinformatic analyses. This study revealed a high content of Alu repeats, genes density, GC frequency, and miRNAs in the great majority of the analyzed breakpoints.

Conclusions

Taken together with literature data about CML with variant t(9;22), our findings identified several new cytogenetic breakpoints as hotspots for recombination, demonstrating that the involvement of chromosomes other than 9 and 22 is not a random event but could depend on specific genomic features. The presence of several genes and/or miRNAs at the identified breakpoints suggests their potential involvement in the CML pathogenesis.

Background

Chronic myeloid leukemia (CML) is characterized by the constitutive expression of the 5'BCR/3'ABL1 fusion gene resulting from the t(9;22)(q34;q11); this translocation is evident in more than 90% of patients and produces the Philadelphia chromosome (Ph)[1].

In 5-10% of CML patients, the 5'BCR/3'ABL1 fusion gene arises from complex variant rearrangements which may involve one or more chromosomes in addition to 9 and 22 [2,3]. In some variant t(9;22) cases, additional material is transferred onto the Ph chromosome, resulting in a "masked" Ph whereas other CML patients show a classic Ph and an atypical der(9) chromosome as a consequence of a rearrangement between the der(9)t(9;22) and another chromosome [4,5]. Serial translocations or a single simultaneous event are alternative hypotheses proposed to justify the occurrence of these complex rearrangements [6].

In a subset of CML patients, cryptic rearrangements have been postulated to induce the chimeric gene formation, such as a nonreciprocal insertion between chromosomes 9 and 22 or two sequential translocations restoring the partner chromosomes morphology [7-11].

Microdeletions on the der(9) chromosome next to the t(9;22) breakpoint have been described in patients with classic and variant Ph translocations, and appear to be a valuable prognostic factor [12-17]. Recently, the frequency of such deletions has been investigated in the subgroup of CML patients with a masked Ph chromosome [18]. Additional genomic deletions on the third derivative chromosome have also been described in CML cases with variant translocations [19,20].

To our knowledge, an accurate breakpoints identification and bioinformatic analysis of other chromosomes involved in variant t(9;22) or in concomitant chromosomal rearrangements apart from the t(9;22) has never been performed in CML.

In this paper, a detailed molecular cytogenetic characterization of 50 (11.1%) out of 452 chronic phase (CP) CML patients was carried out to define the precise breakpoints on chromosomes other than 9 and 22. Bioinformatic analysis of breakpoint regions was performed to investigate the presence of repeated elements (Alu, LINE), GC content, Segmental Duplications (SDs), miRNAs, and known genes. Our findings, taken together with a review of literature data, allowed us to identify new cytogenetic hotspots in CML cases with variant t(9;22).

Methods

Patients

The study included 452 CP-CML patients. All of them were newly diagnosed at our hospital between 1990 and 2009. As a consequence of the long time span for sample accrual, several therapeutic regimens (hydroxyurea, interferon-α, imatinib, nilotinib, and dasatinib) were adopted.

Out of the 452 CP-CML cases, 50 showed variant t(9;22) or additional chromosomal rearrangements, 9 of which have been characterized in previous reports by our group [19-23].

Conventional cytogenetics

Conventional cytogenetic analysis of a 24-48 hour culture was performed at diagnosis of CML on bone marrow cells by standard techniques and evaluated by Giemsa-Trypsin-Giemsa (GTG) banding at about the 400-band level according to the ISCN [24]. At least 25 metaphases were analyzed for each case.

Identification of cytogenetic hotspots

To identify new cytogenetic hotspots, an estimate of the Haploid Autosomal Length (HAL) of the bands involved in variant t(9;22) cases was performed [25,26]. We calculated the number of breaks expected (E) in any band, given the null hypothesis of a random distribution of all breaks across the genome. Reviewing large series of CML patients with variant t(9;22) we assessed the number of breaks observed (O) in each band and divided this value by the expected (E) value to determine an O/E ratio. An O/E ratio >1 identified new cytogenetic hotspots.

FISH analysis

FISH analysis was performed on bone marrow samples of all CP-CML patients at diagnosis using "home-brew" FISH probes specific for ABL1 and BCR genes, validated in previous papers [13,16,27]. Breakpoints characterization and deletions size definition were carried out with additional bacterial artificial chromosome (BAC) and Phage P1-derived artificial chromosome (PAC) probes. All clones were selected according to the University of California Santa Cruz (UCSC http://genome.ucsc.edu/index.html webcite; March 2006 release) database [28]; the mapping of each clone was first tested on normal human metaphases. Chromosome preparations were hybridized in situ with probes labeled with biotin by nick translation [29]. Briefly, 500 ng of labeled probe were used for FISH experiments; hybridization was performed at 37°C in 2× standard saline citrate (SSC), 50% (vol/vol) formamide, 10% (wt/vol) dextran sulphate, 5 μg COT1 DNA (Bethesda Research Laboratories, Gaithersburg, MD), and 3 μg sonicated salmon sperm DNA in a volume of 10 μL. Post-hybridization washing was done at 60°C 0.1× SSC. Biotin-labeled DNA was detected with Cy3-conjugated avidin. In cohybridization experiments, other probes were directly labeled with fluorescein. Chromosomes were identified by 4',6-diamidino-2-phenylindole (DAPI) staining. Digital images were obtained using a Leica DMRXA epifluorescence microscope equipped with a cooled CCD camera (Princeton Instruments, Boston, MA). Cy3 (red; New England Nuclear, NJ), fluorescein (green; NEN Life Science Products, Boston, MA), and DAPI (blue) fluorescence signals, which were detected using specific filters, were recorded separately as gray-scale images. Pseudocoloring and merging of images were performed with Adobe Photoshop software.

Bioinformatic analysis

Breakpoint regions on other chromosomes involved in variant t(9;22) and additional rearrangements were included in 250 Kb size intervals, according to the resolution limit of the BAC clones used for breakpoints definition. Each interval was checked for the presence of interspersed repeats classes (Alu and LINE repeats), SDs, GC content, and gene density. The UCSC Table Browser [28] was queried for summary analysis about the items belonging to the tracks "RepeatMasker", "Segmental Dups", "GC Percent", and "RefSeq Genes". For SDs and RefSeq gene analysis, both "Item count" and "Item Bases" values were considered, to assess their number and the bases percentage involved in SDs or coding sequences, respectively. For each genomic feature, the obtained value was normalized to the mean value for the examined chromosome. For example, in case 1, the breakpoint mapped in 1q32.1 (chr1:203,949,120-204,199,120) showed an Alu frequency of 13.47%. As the mean Alu content inside chromosome 1 is estimated to be 11.9%, the normalized value will be 1.13 (i.e. 13.47/11.90). Therefore, greater or lesser values than 1 correspond to regions with a richer or poorer content of a specific genomic feature than those observed along the entire chromosome.

In view of the known low miRNAs density in the human genome, regions spanning 2 Mb proximally and distally to breakpoints were investigated by querying the UCSC database at the track "sno/miRNA". For each chromosome the expected miRNA density within a 4 Mb interval was established according to the following formula: number of miRNA along the entire chromosome/size in bp of the chromosome × 4000000 bp. The identification of the predicted miRNAs target genes was performed by querying the miRGen database http://www.diana.pcbi.upenn.edu/cgi-bin/miRGen/v3/Targets.cgi webcite. Intersection data from the three widely used target prediction programs (miRanda, PicTar, TargetScan) were considered. The definition of target genes function as oncogenes or tumor suppressor genes (TSGs) was made according to the National Center for Biotechnology Information (NCBI, http://www.ncbi.nlm.nih.gov/gene/ webcite) database.

Results

FISH data

Cytogenetic analysis and FISH experiments with specific BAC/PAC probes for the ABL1 and BCR genes allowed us to detect 50 (11.1%) out of 452 cases, that identify 3 main subgroups of CML patients showing variant t(9;22) rearrangements, the occurrence of cryptic insertions of the ABL1 in the BCR region (or vice versa), and the presence of additional chromosomal abnormalities, respectively (Table 1).

Table 1. Cytogenetic groups identified in a large series of CML patients at diagnosis

1. Variant t(9;22) rearrangements

Forty-three (9.5%) out of 452 CML patients showed the involvement of one (90.7%) or more chromosomes (9.3%) in addition to 9 and 22. These complex variant translocations generated a classic Ph together with a masked der(9) in 36 out of 43 cases (83.7%) and a masked Ph in association with a classic der(9) chromosome in 7 patients (16.3%) (Table 1). Cases with a masked der(9) showed the presence of additional material belonging to partner chromosomes other than chromosome 22 (Fig. 1A, B). Several chromosomes were involved in these variant translocations, with a prevalence of chromosomes 4, 6, 12, and 17 (Table 1). The 5'BCR/3'ABL1 fusion gene was localized on the Ph chromosome in all these cases, whereas the 5'ABL1/3'BCR gene was retained on the der(9) only in 4 (11.1%) out of 36 patients (Table 1). In the remaining 32 (88.9%) cases, the 5'ABL1/3'BCR gene was not detected on der(9) due to deletions and/or 3'BCR transfer onto partner chromosomes (Table 1; Fig. 1A, B). Molecular cytogenetic characterization performed to verify the presence of microdeletions at the level of the rearrangements breakpoints revealed sequences loss in 18 out of the 43 (42%) cases. Among these 18 patients, 10 (55.6%) showed microdeletions of sequences belonging to the third partner chromosome, revealing a high incidence of this kind of deletions in t(9;22) variant rearrangement cases (Table 1).

thumbnailFigure 1. FISH pattern observed with specific probes for the ABL1 and BCR genes on bone marrow metaphases from the analyzed CML patients. Examples from each of the identified cytogenetic groups are shown: "masked der(9)" (A, B), "masked Ph" (C, D), "cryptic insertions" (E, F), and "concomitant rearrangements" (G, H).

Among 7 CML cases with a "masked Ph" chromosome, 3 showed the 5'BCR/3'ABL1 fusion signal on 22q11, the second breakpoint on the derivative chromosome 22 mapping inside chromosome 9 sequences distal to the ABL1 gene (Table 1; Fig. 1C). In 3 cases the fusion gene was detected on the third partner chromosome, the second chromosome 22 breakpoint being localized centromerically to the BCR gene (Table 1; Fig. 1D). In patient #42 the insertion of 5' BCR into the ABL1 gene caused the 5'BCR/3'ABL1 localization on der(9). The 5'ABL1/3'BCR gene was detected on the der(9) in 3 cases with masked Ph, was deleted in case #39 whereas in the remaining patients the 5'ABL1 gene was retained on the der(9) and the 3'BCR gene was transferred onto other derivative chromosomes (Table 1). Chromosome 9 sequences loss next to the rearrangement breakpoint was observed in case #43 and an unusual loss of a region of about 400 Kb localized telomerically to the ABL1 gene was detected in case #41 [22] (Table 1).

2. Cryptic insertions

Six (1.3%) out of the 452 CML cases showed cryptic insertions of ABL1 into BCR, or vice versa, as the cause of the 5'BCR/3'ABL1 fusion gene generation (Fig. 1E, F). Four (66.7%) of these cases are indicated as "Ph negative" (Ph-), with chromosome 22 appearing normal without the presence of additional genomic material (Table 1). Two (33.3%) out of these 6 cases were also included in the previous group as they showed variant rearrangements generating a masked Ph (Table 1). The 5'BCR/3'ABL1 gene was detected on the der(9) or on the der(22) at a ratio of 1:1 as a consequence of 5' BCR insertion in 9q34 or 3' ABL1 insertion in 22q11, respectively (Table 1).

3. Chromosomal rearrangements concomitant to the presence of 5'BCR/3'ABL1

Conventional and molecular cytogenetic analysis showed 5 (1.1%) out of 452 CML cases bearing additional chromosomal rearrangements concomitant to the generation of the 5'BCR/3'ABL1 fusion gene (Table 1; Fig. 1G, H). Cases #46 and #47 were also included in the previous patients group as they showed cryptic insertions; the remaining 3 CML cases carried a classic Ph chromosome.

Bioinformatic analysis of breakpoints on other chromosomes involved in variant t(9;22) or in concomitant chromosomal rearrangements

FISH experiments with BAC clones specific for other chromosomes involved in variant or additional chromosomal rearrangements revealed a total number of 58 breakpoints. These breakpoints were mapped within a single BAC clone or in the region between two overlapping or adjacent clones (Table 2). In cases with sequences loss, two different breakpoints were defined at the level of the deleted regions boundaries.

Table 2. Bioinformatic analysis of the analyzed breakpoints

Interestingly, the majority of breakpoints on chromosomes involved in variant or additional chromosomal rearrangements showed a high frequency of Alu repeats (Table 2; Fig. 2A). In fact, 41 out of 58 (71%) breakpoints showed an Alu content of more than one whereas the remaining 17 out of 58 (29%) had a content of less than one. Instead, the LINE content was lower than one in 44 out of 58 (76%) breakpoints (Table 2). Thirty-five out of 41 breakpoints (85%) with Alu >1 showed a LINE amount < 1 (Table 2).

thumbnailFigure 2. Histograms representing the frequency distribution of genomic features in the breakpoint regions. All values shown for Alu (A), RefSeq Genes Item count (B), GC (C), and miRNAs (D) are reported in logarithmic scale. Chromosomal bands with a null value were excluded from the analysis. The observed distribution of miRNAs within 4 Mb intervals (brown) was compared to the expected mean value for each chromosome (orange).

Most of the analyzed breakpoints map within gene-rich regions as a RefSeq Genes Item count of more than one was observed in 45 out of 58 (78%) breakpoints (Table 2; Fig. 2B). Moreover, 40 out of 58 (69%) breakpoints showed a RefSeq Genes Item bases value of more than 1 (Table 2). It is worthy of note that 34 out of 41 bp (83%) with Alu >1 showed a RefSeq Genes Item count >1 (Table 2). The number of known genes localized at breakpoints and their function as oncogenes and/or TSGs are reported in Table 3.

Table 3. Genes with known function identified in breakpoint regions

In the search for SDs, 49 out of 58 (84%) and 51 out of 58 (88%) breakpoints revealed SDs Item count and SDs Item bases of less than one, respectively (Table 2). In cases showing the presence of SDs within breakpoint regions no specific association with chromosomes 9 and 22 was detected, as the duplicated elements recognized several chromosomal regions.

Finally, a GC content >1 was detected in 43 out of 58 (74%) breakpoints (Table 2; Fig. 2C). A GC content of more than one was identified in 34 out of 41 (83%) breakpoints with Alu >1 and in 34 out of 45 (76%) with a RefSeq Genes Item count >1 (Table 2).

The search for miRNAs revealed a different density from the expected value in 32 out of 58 (55%) breakpoint regions (Fig. 2D). In detail, in 30 (94%) and 2 out of 32 (6%) breakpoints a higher or lower number of miRNA than the expected value was identified, respectively (Fig. 2D). In the remaining 26 out of 58 (45%) breakpoints no miRNA was revealed in the 4 Mb analyzed intervals. It is noteworthy that in case #49 with an additional t(14;15)(q32;q24) a miRNA cluster of 54 elements was revealed in the 14q32 breakpoint region. In this patient a microdeletion of about 450 Kb was detected on 14q32, resulting in the loss of almost the entire miRNA cluster. The list of miRNAs found at the breakpoints is reported in Table 4; in addition to the 14q32 miRNA cluster a total number of 63 known miRNA was identified, 8 (13%) of which show involvement in hematological malignancies. Moreover, querying the miRGen database (the intersection data from the miRanda, PicTar, and TargetScan programs) allowed the identification of the predicted target genes in 19 out of 63 (30%) analyzed miRNAs (see Additional File 1). Among the identified target genes, several play a role as oncogenes or TSGs (see Additional File 1). Noteworthy, some miRNAs share the same target oncogenes or TSGs; for example, PPM1D (protein phosphatase, Mg2+/Mn2+ dependent, 1D) and AKT3 (v-akt murine thymoma viral oncogene homolog 3) genes are the most frequent miRNAs targets (see Additional File 1).

Table 4. Known miRNA mapped at the breakpoint regions.

Additional file 1. Oncogenes and TSGs regulated by the analyzed miRNAs. The number of the predicted target genes for each analyzed miRNA was reported according to the miRGen database. Target genes with a role as oncogenes or TSGs were identified by querying the NCBI database.

Format: XLS Size: 23KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Identification of cytogenetic hotspots

Our study revealed 46 cytogenetic breakpoints on other chromosomes involved in variant t(9;22) rearrangements (see Additional File 2). The assessment of the O/E ratio for each breakpoint allowed us to identify 24 hotspots, 12 of which have been previously described in literature [26] (see Additional File 2). Notably, 4 out of 12 new hotspots showed a ratio >2 involving the chromosomal bands 4q12, 9p11, 11q21 and 21q22 (see Additional File 2).

Additional file 2. Variant t(9;22) Breakpoints. The chromosomal bands involved in variant t(9;22) are shown according to our study and recent large series of CML patients reported in literature. The number of total breaks observed in each band is shown, together with the observed/expected (O/E) ratio. In bold are indicated cytogenetic hotspots, O/E ratio being > 1. The symbol * represents hotspots previously reported [26].

Format: XLS Size: 81KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

To investigate the breakpoints distribution in the genome, a review of literature data about variant t(9;22) following the study by Fisher et al. was carried out [4,30-33]. In total, 60 new hotspots were identified, 18 of which have already been reported [26]. However, 10 previously published hotspots were not supported by our literature review (see Additional File 2). Among the 60 new hotspots, 27 showed an O/E ratio > 2.

Treatment response

Data on the response to treatment in the analyzed CML patients were only available for about 50% of the cases; a summary is shown in Table 5. All the cases evaluable for the response to interferon-α therapy were non responders whereas 11 out of 17 (65%) cases treated with imatinib achieved cytogenetic response. Among patients resistant to imatinib, 3 (75%) treated with dasatinib achieved CCyR.

Table 5. Response to treatment of 50 CML patients included in the study

Discussion

Literature data indicate that breakpoints on additional chromosomes involved in CML cases with variant t(9;22) are not distributed randomly in the genome but show hotspots [26]. Several genomic features such as the density of CpG islands, genes, Alu repeats, recombination events, openness of the chromatin structure and transcription activity have been correlated to the occurrence of breakpoints in variant t(9;22) cases [26,33,34].

In this study, we have performed for the first time a precise molecular cytogenetic characterization of breakpoints involved in variant t(9;22) or in additional rearrangements, in 50 CML cases. To identify genomic elements with a role in the occurrence of chromosomal translocations, bioinformatic analysis was carried out to investigate the distribution and density of several genomic features, such as Alu, LINE, GC, SDs, miRNAs, and genes at breakpoint regions. To date, according to the miRBase database http://www.mirbase.org webcite[35] the total number of known miRNAs is very low (about 720) as compared to the human genome size (3.1 × 109 bp). In this study the miRNAs density within the 4 Mb analyzed intervals resulted higher than the expected value in 32 out of 58 (55%) breakpoint regions. These findings suggest a potential role for miRNAs in the pathogenesis of CML cases with variant or additional chromosomal rearrangements. Few miRNAs located at breakpoint regions have previously been described in several hematological malignancies [36-48]. However, none of them was involved in CML. It is worth noting the presence of the miRNA cluster next to the breakpoint region in 14q32 (case #49). miRNAs in this region are organized in an imprinted domain regulated by a differentially methylated region located upstream of the miRNA cluster. It has been reported that these miRNAs act as tumor suppressor genes and that changes in their methylation status could promote tumor development [49]. Querying of miRGen and NCBI databases showed the involvement of interesting target oncogenes or TSGs implicated in a wide variety of biological processes including cell proliferation, differentiation, apoptosis, and tumorigenesis.

Increasing evidence shows a high density of interspersed repetitive elements, such as Alu and LINE, at some chromosomal translocation breakpoints, suggesting a mediator role of some recurrent rearrangements in tumors [50]. Because a much higher density of Alu repeats has been observed in the DNA sequences flanking the ABL1 and BCR genes, it has been hypothesized that Alu elements provide hotspots for non allelic homologous recombination and mediate chromosomal translocation in CML [34,50]. Our data, supported by bioinformatic evidence, suggest that the high density of Alu repeats could increase the propensity to undergo rearrangements also of other chromosomes involved in variant t(9;22). In our CML series, a high Alu density was detected in 71% of the analyzed breakpoints. Moreover, a rich content of Alu repeats was revealed also on breakpoint regions identified in chromosomal rearrangements concomitant to the t(9;22).

Literature data revealed a preferential breakpoints distribution in CML cases with variant t(9;22) within the CG-richest regions of the genome corresponding to the G-light banding karyotype [26,33]. Our data confirmed this association, as 83% of the identified cytogenetic breakpoints mapped inside G-light bands. Moreover, we report the first bioinformatic evidence of the association between GC-content and breaks in cases with variant t(9;22), as 73% of the molecular breakpoints showed a GC content >1. In addition, these data showed that CG richness was related to other genomic features such as Alu content and a greater gene density than the mean expected value.

The search for SDs revealed a low density in the majority of the analyzed breakpoints, without showing any specific association with chromosomes 9 and 22 regions, unlike what has been reported about the occurrence of the t(9;22) in CML [51].

Moreover, our study provided an outline of the frequency and molecular features of the most relevant cytogenetic groups identified in a very large series of CML patients at diagnosis. Three-way translocations were the most frequent among variant t(9;22) rearrangements, chromosomes 4, 6, 12, and 17 being common partners. However, no cytogenetic breakpoints clustering was revealed when the same partner chromosome was rearranged, except for the 3p21 band, that was involved in 3 CML cases with variant t(9;22).

As to the mechanisms involved in the formation of the variant t(9;22) rearrangements, our data indicated that the most probable mechanism, identified in cases with a "masked der(9)" chromosome, is a single event consisting of multiple simultaneous breaks and rejoins (one-step model). In fact, splitting of the 5'ABL1/3'BCR fusion gene signal was observed in the majority (27 out of 36, 75%) of analyzed cases. A two-step mechanism was hypothesized in about 11% of cases bearing a "masked der(9)" chromosome; the permanence of the 5'ABL1/3'BCR gene on the der(9) suggests that a second break occurred inside the chromosome 22 sequence telomeric to the BCR gene. On the contrary, in 71.4% of cases (#37 - #41) with a "masked Ph" chromosome a second break located proximally to BCR or distally to ABL1 was identified, suggesting the occurrence of a two-step mechanism in the majority of the CML patients included in this group.

In our study, FISH 'walking' with BAC/PAC contigs belonging to the chromosomes 9 and 22 next to the t(9;22) breakpoint regions allowed us to assess the frequency of deletions in three main cytogenetic subgroups of CML patients and the size of these microdeletions. Confirming the deletion frequency reported in literature [12], 12 out of the 36 (33%) cases with a "masked der(9)" chromosome showed chromosome 9 and/or 22 sequences loss. Moreover, in about 55% of these patients we found extensive genomic deletions on the third chromosome, in addition to deletions on der(9). Chromosome 9 sequences deletions were detected in 3 out of 6 (50%) cases with a masked Ph (#39, #41, and #43) and in 1 out of 4 (25%) Ph- cases (#44). These frequencies are higher than the value recently reported in literature [18].

The biological significance and the prognostic impact of the cytogenetic molecular heterogeneity occurring in the generation of the 5'BCR/3'ABL1 fusion gene remain to be clarified. However, the bioinformatic analysis performed in this study on a large number of breakpoints in CML cases with variant t(9;22) or additional chromosomal alterations revealed that the rearranged regions are characterized by an elevated content of miRNAs, Alu repeats, GC and known genes.

In conclusion, this genomic analysis of breakpoint regions provides clues to a better understanding of the pathogenetic mechanisms that underlie CML onset. Further analyses will be needed to demonstrate the functional meaning of these genomic features.

List of abbreviations

(CP-CML): Chronic Myeloid Leukemia; (FISH): Fluorescence In Situ Hybridization; (CML): Chronic myeloid leukemia; (Ph): Philadelphia chromosome; (SDs): Segmental Duplications; (GTG): Giemsa-Trypsin-Giemsa; (HAL): Haploid Autosomal Length; (E): Number of breaks expected; (O): Number of breaks observed; (BAC): Bacterial artificial chromosome; (PAC): Phage P1-derived artificial chromosome; (UCSC): University of California Santa Cruz; (Ph-): Ph negative.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

FA, LA, and AZ were involved in the design and execution of the experiments, wrote the manuscript and contributed to the overall experimental design. NC conducted most of the FISH experiments. PC and LV performed conventional cytogenetic analysis; ARR contributed to clinical data collection. VL, MR and GS participated in the design of the study and supervised the manuscript preparation. All authors have read and approved the final manuscript.

Acknowledgements

The authors thank Ms MVC Pragnell, B.A. for language revision of the manuscript. The financial support of Associazione Italiana contro le Leucemie (AIL)-BARI is gratefully acknowledged.

References

  1. Melo JV, Barnes DJ: Chronic myeloid leukaemia as a model of disease evolution in human cancer.

    Nat Rev Cancer 2007, 7:441-453. PubMed Abstract | Publisher Full Text OpenURL

  2. Huret JL: Complex translocations, simple variant translocations and Ph-negative cases in chronic myelogenous leukaemia.

    Hum Genet 1990, 85:565-568. PubMed Abstract | Publisher Full Text OpenURL

  3. Mitelman F, Johansson B, Mertens F: The Mitelman database of chromosome aberrations in cancer. [http://www.cgap.nci.nih.gov/Chromosomes/Mitelman] webcite

  4. Gorusu M, Benn P, Li Z, Fang M: On the genesis and prognosis of variant translocations in chronic myeloid leukemia.

    Cancer Genet Cytogenet 2007, 173:97-106. PubMed Abstract | Publisher Full Text OpenURL

  5. Virgili A, Brazma D, Reid AG, Howard-Reeves J, Valgañón M, Chanalaris A, De Melo VA, Marin D, Apperley JF, Grace C, Nacheva EP: FISH mapping of Philadelphia negative BCR/ABL1 positive CML.

    Mol Cytogenet 2008, 1:14. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  6. Fitzgerald PH, Morris CM: Complex chromosomal translocations in the Philadelphia chromosome leukemias: serial translocations or a concerted genomic rearrangement?

    Cancer Genet Cytogenet 1991, 57:143-151. PubMed Abstract | Publisher Full Text OpenURL

  7. Plas DC, Grosveld G, Hagemeijer A: Review of clinical, cytogenetic, and molecular aspects of Ph-negative CML.

    Cancer Genet Cytogenet 1991, 52:143-156. PubMed Abstract | Publisher Full Text OpenURL

  8. Chase A, Huntly BJP, Cross NCP: Cytogenetics of chronic myeloid leukemia.

    Best Pract Res Clin Haematol 2001, 14:553-571. PubMed Abstract | Publisher Full Text OpenURL

  9. Haigh S, Cuthbert G: Fluorescence in situ hybridization characterization of different cryptic BCR-ABL rearrangements in chronic myeloid leukemia.

    Cancer Genet Cytogenet 2004, 155:132-137. PubMed Abstract | Publisher Full Text OpenURL

  10. Reddy KS, Grove B: A Philadelphia-negative chronic myeloid leukemia with a BCR/ABL fusion gene on chromosome 9.

    Cancer Genet Cytogenet 1998, 107:48-50. PubMed Abstract | Publisher Full Text OpenURL

  11. Storlazzi CT, Anelli L, Surace C, Lonoce A, Zagaria A, Nanni M, Curzi P, Rocchi M: Molecular cytogenetic characterization of a complex rearrangement involving chromosomes 9 and 22 in a case of Ph negative chronic myeloid leukemia.

    Cancer Genet Cytogenet 2002, 136:141-145. PubMed Abstract | Publisher Full Text OpenURL

  12. Huntly BJ, Reid AG, Bench AJ, Campbell LJ, Telford N, Shepherd P, Szer J, Prince HM, Turner P, Grace C, Nacheva EP, Green AR: Deletions of the derivative chromosome 9 occur at the time of the Philadelphia translocation and provide a powerful and independent prognostic indicator in chronic myeloid leukemia.

    Blood 2001, 98:1732-1738. PubMed Abstract | Publisher Full Text OpenURL

  13. Storlazzi CT, Specchia G, Anelli L, Albano F, Pastore D, Zagaria A, Rocchi M, Liso V: Breakpoint characterization of der(9) deletions in CML patients.

    Genes Chromosomes Cancer 2002, 35:271-276. PubMed Abstract | Publisher Full Text OpenURL

  14. Huntly BJ, Guilhot F, Reid AG, Vassiliou G, Hennig E, Franke C, Byrne J, Brizard A, Niederwieser D, Freeman-Edward J, Cuthbert G, Bown N, Clark RE, Nacheva EP, Green AR, Deininger MW: Imatinib improves but may not fully reverse the poor prognosis of patients with CML with derivative chromosome 9 deletions.

    Blood 2003, 102:2205-2212. PubMed Abstract | Publisher Full Text OpenURL

  15. Reid AG, Huntly BJ, Grace C, Green AR, Nacheva EP: Survival implications of molecular heterogeneity in variant Philadelphia-positive chronic myeloid leukaemia.

    Br J Haematol 2003, 121:419-427. PubMed Abstract | Publisher Full Text OpenURL

  16. Specchia G, Albano F, Anelli L, Storlazzi CT, Zagaria A, Liso A, Pannunzio A, Pastore D, Mestice A, Greco G, Liso V, Rocchi M: Derivative chromosome 9 deletions in chronic myeloid leukemia are associated with loss of tumor suppressor genes.

    Leuk Lymphoma 2004, 45:689-694. PubMed Abstract | Publisher Full Text OpenURL

  17. Quintas-Cardama A, Kantarjian H, Talpaz M, O'brien S, Garcia-Manero G, Verstovsek S, Rios MB, Hayes K, Glassman A, Bekele BN, Zhou X, Cortes J: Imatinib mesylate therapy may overcome the poor prognostic significance of deletions of derivative chromosome 9 in patients with chronic myelogenous leukemia.

    Blood 2005, 105:2281-2286. PubMed Abstract | Publisher Full Text OpenURL

  18. De Melo VA, Milojkovic D, Marin D, Apperley JF, Nacheva EP, Reid AG: Deletions adjacent to BCR and ABL1 breakpoints occur in a substantial minority of chronic myeloid leukemia patients with masked Philadelphia rearrangements.

    Cancer Genet Cytogenet 2008, 182:111-115. PubMed Abstract | Publisher Full Text OpenURL

  19. Albano F, Specchia G, Anelli L, Zagaria A, Storlazzi CT, Buquicchio C, Roberti MG, Liso V, Rocchi M: Genomic deletions on other chromosomes involved in variant t(9;22) chronic myeloid leukemia cases.

    Genes Chromosomes Cancer 2003, 36:353-360. PubMed Abstract | Publisher Full Text OpenURL

  20. Anelli L, Albano F, Zagaria A, Liso A, Roberti MG, Rocchi M, Specchia G: A chronic myelocytic leukemia case bearing deletions on the three chromosomes involved in a variant t(9;22;11).

    Cancer Genet Cytogenet 2004, 148:137-140. PubMed Abstract | Publisher Full Text OpenURL

  21. Zagaria A, Anelli L, Albano F, Storlazzi CT, Liso A, Roberti MG, Buquicchio C, Liso V, Rocchi M, Specchia G: A fluorescence in situ hybridization study of complex t(9;22) in two chronic myelocytic leukemia cases with a masked Philadelphia chromosome.

    Cancer Genet Cytogenet 2004, 150:81-85. PubMed Abstract | Publisher Full Text OpenURL

  22. Zagaria A, Anelli L, Albano F, Vicari L, Schiavone EM, Annunziata M, Pane F, Liso V, Rocchi M, Specchia G: Molecular cytogenetic characterization of deletions on der(9) in chronic myelocytic leukemia.

    Cancer Genet Cytogenet 2006, 167:97-102. PubMed Abstract | Publisher Full Text OpenURL

  23. Albano F, Zagaria A, Anelli L, Pannunzio A, Manodoro F, Coccaro N, Russo Rossi A, Liso V, Rocchi M, Specchia G: The double deceit generated by an insertion mechanism in chronic myeloid leukemia with t(9;9;22).

    Ann Hematol 2008, 87:923-926. PubMed Abstract | Publisher Full Text OpenURL

  24. Mitelman F: ISCN. An International System for Human Cytogenetic Nomenclature. Basel, Switzerland: Karger; 1995.

  25. Daniel A: Structural differences in reciprocal translocations. Potential for a model of risk in Rcp.

    Hum Genet 1979, 51:171-182. PubMed Abstract | Publisher Full Text OpenURL

  26. Fisher AM, Strike P, Scott C, Moorman AV: Breakpoints of variant 9;22 translocations in chronic myeloid leukemia locate preferentially in the CG-richest regions of the genome.

    Genes Chromosomes Cancer 2005, 43:383-389. PubMed Abstract | Publisher Full Text OpenURL

  27. Albano F, Anelli L, Zagaria A, Archidiacono N, Liso V, Specchia G, Rocchi M: "Home-brew" FISH assay shows higher efficiency than BCR-ABL dual color, dual fusion probe in detecting microdeletions and complex rearrangements associated with t(9;22) in chronic myeloid leukemia.

    Cancer Genet Cytogenet 2007, 174:121-126. PubMed Abstract | Publisher Full Text OpenURL

  28. The University of California Santa Cruz database [http://genome.ucsc.edu/index.html] webcite

  29. Lichter P, Tang Chang CJ, Call K, Hermanson G, Evans GA, Housman D, Ward DC: High resolution mapping of human chromosomes 11 by in situ hybridization with cosmid clones.

    Science 1990, 247:64-69. PubMed Abstract | Publisher Full Text OpenURL

  30. El-Zimaity MM, Kantarjian H, Talpaz M, O'Brien S, Giles F, Garcia-Manero G, Verstovsek S, Thomas D, Ferrajoli A, Hayes K, Nebiyou Bekele B, Zhou X, Rios MB, Glassman AB, Cortes JE: Results of imatinib mesylate therapy in chronic myelogenous leukaemia with variant Philadelphia chromosome.

    Br J Haematol 2004, 125:187-195. PubMed Abstract | Publisher Full Text OpenURL

  31. Costa D, Carrió A, Madrigal I, Arias A, Valera A, Colomer D, Aguilar JL, Teixido M, Camós M, Cervantes F, Campo E: Studies of complex Ph translocations in cases with chronic myelogenous leukemia and one with acute lymphoblastic leukemia.

    Cancer Genet Cytogenet 2006, 166:89-93. PubMed Abstract | Publisher Full Text OpenURL

  32. Richebourg S, Eclache V, Perot C, Portnoi MF, Akker J, Terré C, Maareck O, Soenen V, Viguié F, Laï JL, Andrieux J, Corm S, Roche-Lestienne C, Fi-LMC Group: Mechanisms of genesis of variant translocation in chronic myeloid leukemia are not correlated with ABL1 or BCR deletion status or response to imatinib therapy.

    Cancer Genet Cytogenet 2008, 182:95-102. PubMed Abstract | Publisher Full Text OpenURL

  33. Bennour A, Sennana H, Laatiri MA, Elloumi M, Khelif A, Saad A: Molecular cytogenetic characterization of variant Philadelphia translocations in chronic myeloid leukemia: genesis and deletion of derivative chromosome 9.

    Cancer Genet Cytogenet 2009, 194:30-37. PubMed Abstract | Publisher Full Text OpenURL

  34. Jeffs AR, Wells E, Morris CM: Nonrandom distribution of interspersed repeat elements in the BCR and ABL1 genes and its relation to breakpoint cluster regions.

    Genes Chromosomes Cancer 2001, 32:144-154. PubMed Abstract | Publisher Full Text OpenURL

  35. The miRBase database [http://www.mirbase.org] webcite

  36. Visone R, Rassenti LZ, Veronese A, Taccioli C, Costinean S, Aguda BD, Volinia S, Ferracin M, Palatini J, Balatti V, Alder H, Negrini M, Kipps TJ, Croce CM: Karyotype-specific microRNA signature in chronic lymphocytic leukemia.

    Blood 2009, 114:3872-3879. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Zanette DL, Rivadavia F, Molfetta GA, Barbuzano FG, Proto-Siqueira R, Silva WA Jr, Falcao RP, Zago MA: miRNA expression profiles in chronic lymphocytic and acute lymphocytic leukemia.

    Braz J Med Biol Res 2007, 40:1435-1440. PubMed Abstract | Publisher Full Text OpenURL

  38. Wang Y, Li Z, He C, Wang D, Yuan X, Chen J, Jin J: MicroRNAs expression signatures are associated with lineage and survival in acute leukemias.

    Blood Cells Mol Dis 2010, 44:191-197. PubMed Abstract | Publisher Full Text OpenURL

  39. Garzon R, Volinia S, Liu CG, Fernandez-Cymering C, Palumbo T, Pichiorri F, Fabbri M, Coombes K, Alder H, Nakamura T, Flomenberg N, Marcucci G, Calin GA, Kornblau SM, Kantarjian H, Bloomfield CD, Andreeff M, Croce CM: MicroRNA signatures associated with cytogenetics and prognosis in acute myeloid.

    Blood 2008, 111:3183-3189. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Fulci V, Colombo T, Chiaretti S, Messina M, Citarella F, Tavolaro S, Guarini A, Foa' R, Macino G: Characterization of B- and T-lineage acute lymphoblastic leukemia by integrated analysis of MicroRNA and mRNA expression profiles.

    Genes Chromosomes Cancer 2009, 48:1069-1082. PubMed Abstract | Publisher Full Text OpenURL

  41. Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, Aldler H, Rattan S, Keating M, Rai K, Rassenti L, Kipps T, Negrini M, Bullrich F, Croce CM: Frequent deletions and down-regulation of micro- RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia.

    Proc Natl Acad Sci USA 2002, 99:15524-1529. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Pons A, Nomdedeu B, Navarro A, Gaya A, Gel B, Diaz T, Valera S, Rozman M, Belkaid M, Montserrat E, Monzo M: Hematopoiesis-related microRNA expression in myelodysplastic syndromes.

    Leuk Lymphoma 2009, 50:1854-1859. PubMed Abstract | Publisher Full Text OpenURL

  43. Hanlon K, Rudin CE, Harries LW: Investigating the targets of MIR-15a and MIR-16-1 in patients with chronic lymphocytic leukemia (CLL).

    PLoS One 2009, 4:e7169. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Corthals SL, Jongen-Lavrencic M, de Knegt Y, Peeters JK, Beverloo HB, Lokhorst HM, Sonneveld P: Micro-RNA-15a and micro-RNA-16 expression and chromosome 13 deletions in multiple myeloma.

    Leuk Res 2010, 34:677-681. PubMed Abstract | Publisher Full Text OpenURL

  45. Dixon-McIver A, East P, Mein CA, Cazier JB, Molloy G, Chaplin T, Andrew Lister T, Young BD, Debernardi S: Distinctive patterns of microRNA expression associated with karyotype in acute myeloid leukaemia.

    PLoS One 2008, 3:e2141. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Katrincsakova B, Takeda H, Urbankova H, Michaux L, Jarosova M, Vandenberghe P, Georges M, Charlier C, Wlodarska I: Methylation analysis of the imprinted DLK1-GTL2 domain supports the random parental origin of the IGH-involving del(14q) in B-cell malignancies.

    Epigenetics 2009, 4:469-475. PubMed Abstract | Publisher Full Text OpenURL

  47. Robbiani DF, Bunting S, Feldhahn N, Bothmer A, Camps J, Deroubaix S, McBride KM, Klein IA, Stone G, Eisenreich TR, Ried T, Nussenzweig A, Nussenzweig MC: AID produces DNA double-strand breaks in non-Ig genes and mature B cell lymphomas with reciprocal chromosome translocations.

    Mol Cell 2009, 36:631-641. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Bellon M, Lepelletier Y, Hermine O, Nicot C: Deregulation of microRNA involved in hematopoiesis and the immune response in HTLV-I adult T-cell leukemia.

    Blood 2009, 113:4914-4917. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Saito Y, Liang G, Egger G, Friedman JM, Chuang JC, Coetzee GA, Jones PA: Specific activation of microRNA-127 with downregulation of the proto-oncogene BCL6 by chromatin-modifying drugs in human cancer cells.

    Cancer Cell 2006, 9:435-443. PubMed Abstract | Publisher Full Text OpenURL

  50. Kolomietz E, Meyn MS, Pandita A, Squire JA: The role of Alu repeat clusters as mediators of recurrent chromosomal aberrations in tumors.

    Genes Chromosomes Cancer 2002, 35:97-112. PubMed Abstract | Publisher Full Text OpenURL

  51. Albano F, Anelli L, Zagaria A, Coccaro N, D'Addabbo P, Liso V, Rocchi M, Specchia G: Genomic segmental duplications on the basis of the t(9;22) rearrangement in Chronic Myeloid Leukemia.

    Oncogene 2010, 29:2509-2516. PubMed Abstract | Publisher Full Text OpenURL