Skip to main content
Download PDF

Short telomere length and its correlation with gene mutations in myelodysplastic syndrome

Journal of Hematology & Oncology volume 9, Article number: 62 (2016) Cite this article

Abstract

Background

Telomere erosion can lead to genomic instability and cancer progression. It has been suggested that the shortest telomere, not the average telomere length (TL), is critical for cell viability. Some studies have shown shorter TL in myelodysplastic syndrome (MDS) patients but the critically short telomeres, the variability of TL within individual patient has not been evaluated. Thus, we aimed to investigate the TL of MDS patients and assessed the association of TL with recurrent genetic mutations in MDS.

Methods

We measured the TL of bone marrow nucleated cells for diagnostic samples at a single-cell level by quantitative fluorescence in situ hybridization (Q-FISH) for 58 MDS patients and analyzed the minimum, median, average, standard deviation, average of the 0th to 10th percentile TL within a patient, and the proportion of cells with TL that is shorter than the lowest 10th percentile of the normal control (NC). The correlations of TL to clinical parameters, cytogenetic results, and genetic mutations were assessed.

Results

MDS patients showed eroded telomeres and narrow distribution compared to the NC (P < 0.001, P = 0.018, respectively). Patients with mutation showed significantly lesser cells with short TL, below the lowest 10th percentile of the NC (P = 0.017), but no differences in TL were found according to mutations/cytogenetic abnormalities except for CSF3R mutation. However, those patients with a high percentage (≥80 %) of cells with short TL showed poorer overall survival (P = 0.021), and this was an independent prognostic factor, along with TP53, U2AF1 mutation, and high BM blast count (P = 0.044, 0.001, 0.004, 0.012, respectively).

Conclusions

The shortest TL, which determines the fate of the cell, was significantly shorter, and higher burden of cells with short TL were found in MDS, which correlated with poor survival, suggesting the need to measure TL in single cells by Q-FISH.

Background

Telomeres are non-coding, repetitive sequences of DNA at the ends of the chromosomes of eukaryotic cells which become shorter as cells divide, and when telomere attrition reaches its limit, cell proliferation arrest, senescence, and apoptosis can occur [1, 2]. Shortened telomeres can cause end-to-end fusion of chromosomes, which results in genomic instability and contributes to carcinogenesis [3]. Thus, telomeres and their maintenance have been studied in many hematologic malignancies [4–6].

Myelodysplastic syndromes (MDS) are clonal hematopoietic stem cell diseases characterized by cytopenia, dysplasia, and ineffective hematopoiesis, with an increased risk of transformation to acute myeloid leukemia (AML) [7, 8]. Due to its predisposition to AML, genomic instability has been studied as the key to the pathogenesis of MDS, and thus, many studies have assessed the telomere components of MDS [9–11]. MDS patients have shown a shorter average telomere length (TL) compared to normal control (NC), but the association of TL to other clinical parameters in patients with MDS has been inconsistent [9, 11, 12].

Many somatic mutations have been identified in patients with MDS, involving various pathways including epigenetic regulation, RNA splicing machinery, and transcription factors [13–15], which help elucidate the pathophysiology of the disease. A recent study showed that telomere dysfunction causes aberrant RNA splicing by repressing the splicing gene expression [16]. However, the correlation of TL with various genetic mutations other than telomerase complex genes has not been thoroughly explored [17, 18].

Moreover, most previous studies have used Southern blot analysis, quantitative PCR, or flow fluorescence in situ hybridization (FISH) for the analysis of TL [9, 11, 19, 20]. However, these studies only assessed the mean or median TL of the patient whereas it is believed that the cell with the shortest telomere is actually the critical factor leading telomere dysfunction [3, 21]. Performing quantitative-FISH (Q-FISH) for telomeres allows the measurement of TL at a single-cell level with the measurement of the distribution of TL of within individuals as well as the shortest TL [22]. Telomeres of individual chromosomes are stained with peptic nucleic acid probes for telomeres, and the fluorescence intensity of the telomeric signal is measured with fluorescence-based automatic microscope [23]. The fluorescence intensity is normalized by the use of an additional centromeric probe, and the telomere to centromere (T/C) fluorescence intensity ratio is considered to be the surrogate of TL since the fluorescence intensity of the telomeric signal is known to be proportional to the TL [24]. Inclusion of the centromeric probe also allows the selection of single cells in which two centromeric signals are observed and the normalization of the signal for hybridization [23]. Performing Q-FISH in interphase nuclei further allows measurement of each value in numerous nucleated cells. Therefore, we measured the TL of bone marrow (BM) nucleated cells of patients with MDS at diagnosis by Q-FISH, evaluated the correlation of the TL with genetic mutations and clinical parameters, and explored the usefulness of TL parameters as a prognostic marker.

Methods

Patients

This study examined 58 patients diagnosed with MDS (37 men and 21 women) at Seoul National University Hospital between January 2004 and December 2014. This study was approved by the institutional review board of Seoul National University Hospital (1604-082-754). The samples included in this study were from diagnostic samples of MDS. The characteristics of the MDS patients are shown in Table 1. BM aspirates from 23 individuals who underwent marrow biopsy and showed no evidence of hematologic malignancy were included as the NC. The bone marrow of the NC showed no evidence of bone marrow involvement and showed normal karyotype.

Table 1 Characteristics of MDS patients (n = 58)
Full size table

Karyotyping—Giemsa banding

A standard protocol was used to perform Giemsa banding with heparinized whole blood marrow samples [25]. A minimum of 20 metaphase cells were karyotyped according to the International System for Human Cytogenetic Nomenclature (ISCN) [26, 27].

Interphase fluorescence in situ hybridization for MDS

Interphase FISH tests (del(5q)/−5, del(7q)/−7, trisomy 8, del(20q), trisomy 1/1q+) were performed. The probes used were Vysis LSI EGR1 (5q31)/D5S23, D5S721 Dual color probe kit, Vysis LSI D7S522 (7q31)/CEP7 FISH probe, Vysis CEP8(Dz82) Spectrum Green Probe, Vysis D20S108 FISH probe, and Vysis LSI 1p36/LSI1q25 and LSI 19q13/19p13 Dual color probe (Abbott Molecular, Downers Grove, IL, USA). These five types of FISH items were selected based on the previous report of recurrent cytogenetic changes in Korean patients with MDS [28]. FISH was performed with mononuclear cells of BM aspirates at the time of BM examination as described in a previous study [28].

Telomere quantitative fluorescence in situ hybridization

Cryopreserved BM nucleated cells in fixative after FISH preparation were used for telomere analysis. Q-FISH was performed using a Cy3-labelled telomere PNA (peptide nucleic acid) FISH kit (DakoCytomation Denmark A/S, Glostrup, Denmark) and a FITC-labeled PNA probe for the centromere of chromosome 2 (kindly provided by DakoCytomation). One microliter of the probe for the centromere of chromosome 2 was added to 10 μL of the telomere probe. Telomere and centromere Q-FISH hybridizations were performed according to the manufacturer’s instructions. Interphase Q-FISH images were captured with a Zeiss Axioplan 2 imaging microscope (Carl Zeiss MicroImaging GmbH, Munchen, Germany) equipped with ISIS software (MetaSystems GmbH, Altlussheim, Germany) (Fig. 1). For TL measurements, the ISIS-Telomere module (MetaSystems) was used as described previously [23]. The software calculates a T/C fluorescence intensity ratio, which is a measure of TL, for individual cells. At least 25 interphase nuclei were scanned for each patient to measure the TL. The cell with the shortest TL is represented as the minimum TL and the lower 25th percentile value of each individual is expressed as the “Q1” TL. The median, the average and the distribution of TL, standard distribution (SD), and the average of the 0th to 10th percentile TL for each individual was calculated from the total cells of which the TL was measured. Furthermore, the percentage of cells with TL shorter than the 10th percentile of the normal control was assessed and stated as “short” TL.

Fig. 1
figure 1

Interphase telomere fluorescence in situ hybridization (FISH) image. FISH results obtained with pan-telomere (cyanine (Cy)3-labeled) and centromere 2 specific (fluorescein isothiocyanate (FITC)-labeled) peptic nucleic acid probe in interphase nuclei. Counter stain was performed with 4′,6-diamidino-2-phenylindole (DAPI) to identify the nuclei

Full size image

Targeted sequencing

Eighty-seven genes were selected for targeted sequencing (Additional file 1: Table S1). Genes that were previously reported as mutational targets in MDS (n = 49) [13, 14, 29] or those mutated in other myeloid and lymphoid hematologic malignancies (n = 11 and n = 27) were included [13, 14, 29–35]. The BM aspirates collected in a sodium ethylenediaminetetraacetic acid tube were treated with ammonium chloride for RBC lysis, and the cell pellets were kept in RNAlater solution (Qiagen, Valencia, CA, USA) at −80 °C until genetic analysis. The genomic DNA (gDNA) was extracted from the buffy coat of BM aspirate using QIAamp DNA Blood Mini Kit (Qiagen), according to the manufacturer’s instructions. The gDNA shearing, standard library production, and hybridization were performed by Celemics Inc. (Seoul, Korea). The quality of the gDNA was assessed by the Agilent 2200 TapeStation System (Agilent, Santa Clara, CA, USA). A total target length of the 259-kb region using the paired-end 150-bp rapid-run sequencing mode was performed on an Illumina Hiseq 2500 platform (Illumina, San Diego, CA, USA). The sequencing data calling process is described in Additional file 1.

Statistical analysis

Categorization of MDS patients was performed according to the 2008 World Health Organization (WHO) classification of tumors of hematopoietic and lymphoid tissues [8], the International Prognostic Scoring System (IPSS), and the revised IPSS et al. [36, 37]. The correlation of TL with various clinical parameters, BM findings, cytogenetic abnormalities, and mutations found in multi-gene panels was assessed.

The gene mutations were evaluated separately and in subgroups according to the function of genes. Ten gene subgroups are splicing machinery, DNA methylation, chromatin modification, transcription factor, receptor/kinases, RAS pathway, cell signaling, DNA repair/cell cycle, cohesion, and miscellaneous. The genes in each subgroup are listed in the Additional file 1: Table S1. The χ2 test and Fisher’s exact test were used to compare categorical variables. Pairwise correlations between TLs and gene mutations, cytogenetic results, and categorical clinical parameters were evaluated using Kendall tau-b (T b) correlation. For linear parameters, Pearson’s r correlation was used. The comparisons of TLs between patients with and without a certain mutation were made with the Mann-Whitney method for genes that showed mutation in more than 5.0 % of the patients.

The overall survival (OS) was calculated from the date of diagnosis to the date of death from any cause and compared by Kaplan-Meier method (log-rank test). Univariable and multivariable Cox analyses were performed to assess variables as prognostic factors for survival. Backward LR selection was performed for multivariable Cox analysis. Harrell’s C-index was calculated for the proposed multivariable Cox model and the revised IPSS. Statistical analyses were performed using the SPSS version 22.0 (IBM Inc., Chicago, IL, USA) and R statistical program (http://www.r-project.org). P values <0.05 were considered to be statistically significant.

Results

Telomere lengths and clinical profiles

We compared the telomere lengths (TLs) of the MDS patients (n = 58) with those of the NC (n = 23), according to several parameters (Table 2). The clinical characteristics of the MDS patients included in this study are shown in Table 1. There was no significant difference in age and sex (P = 0.175, 0.796, respectively) between the MDS patients and NC. However, the number of cells assessed in the MDS patients was smaller (P = 0.004) than that assessed in the NC due to insufficient nuclei in the patient samples. All of the TL parameters were significantly lower in the MDS patients than the NC (P < 0.001 for all the TL except SD, P = 0.018 for SD). The distribution of the telomere length of each individual was narrower in the MDS patients than the NC (P = 0.018).

Table 2 Telomere length parameters of MDS patients and normal control
Full size table

We also calculated the lowest 10th percentile TL value of the NC group and found out that 52.1 (±24.9)% of MDS patients’ cells were under this length (Table 2). These results show that not only the mean TL of the MDS patients is shorter than that of the NC but also a large proportion of MDS patients’ cells are shorter than the lowest 10th percentile TL of the NC.

TL was compared for WHO categories and IPSS risk groups. The TL parameters were not different among the clinical subgroups (the average of the 0th–10th percentile TL and the median TL shown in Fig. 2, the other TL parameters are in the Additional file 1: Table S2). Additional comparison of TL was made according to the complete blood count and BM findings. Patients with low hemoglobin (<8 g/dL) showed a higher percentage of cells with TL less than the 10th percentile of normal control (P = 0.020). However, the other TL parameters were not significantly different in those patients with hemoglobin values less than vs. equal to/greater than 8 g/dL. There were no significant differences in the TL according to absolute neutrophil count, platelet count, BM blast count, fibrosis, or cellularity (Additional file 1: Table S3).

Fig. 2
figure 2

Telomere lengths according to the World Health Organization (WHO) categories, International Prognostic Scoring System (IPSS), and revised IPSS subgroups. T/C telomere to centromere

Full size image

Multi-gene panel results

Among the 87 genes, we found mutations of 43 genes in more than one patient. Forty-five patients (77.6 %) had showed mutation in one or more of the genes. Only the gene mutations found in more than 5 % of the samples are shown in Fig. 3. Grouping the 87 genes by gene function, genes related to splicing were most frequently mutated (34.5 %), following the genes for transcription factors (31.0 %), chromatin modification (25.9 %), DNA methylation (19.0 %), receptor/kinases (15.5 %), and RAS pathways (8.6 %). The associations of gene mutations and/or cytogenetic results are shown in Fig. 4. Kendall tau-b (T b) correlation coefficients greater than 0.40 were present between SRSF2 and SETBP1 mutation (T b = 0.432, P = 0.001), IDH2 and SRSF2 mutation (T b = 0.432, P = 0.001), TP53 mutation and del(5q)/−5 (T b = 0.520, P < 0.001), TP53 mutation and del(7q)/−7 (T b = 0.564, P < 0.001), TP53 mutation and abnormalities of 17p (T b = 0.583, P < 0.001), del(5q)/−5 and abnormality of 17p (T b = 0.431, P = 0.001). Significant correlations were also present between DNMT3A vs. IDH2 mutation, DNMT3A vs. SETBP1 mutation (T b = 0.390, P = 0.003, respectively), U2AF1 mutation vs. trisomy 8 (T b = 0.379, P = 0.005), del(7q)/−7 vs. abnormal 17p (T b = 0.357, P = 0.008), SF3B1 vs. TET2 mutation (T b = 0.343, P = 0.010), DNMT3A vs. TET2 mutation (T b = 0.263, P = 0.047), IDH2 vs. RUNX1 mutation, RUNX1 vs. SETBP1 mutation, IDH2 vs. SETBP1 mutation (T b = 0.297, P = 0.025, respectively). Comparing the association between gene categories and cytogenetic abnormalities; DNA methylation vs. cell signaling-related genes (T b = 0.391, P = 0.003), DNA methylation vs. DNA repair/cell cycle-related genes (T b = 0.389, P = 0.003), RAS pathway vs. splicing-related genes (T b = 0.294, P = 0.026), splicing-related genes vs. trisomy 8 (T b = 0.338, P = 0.011), transcription factor genes vs. del(5q)/−5 (T b = 0.276, P = 0.039), transcription factor genes vs. del(7q)/−7 (T b = 0.399, P = 0.003), transcription factor genes vs. abnormality of 17p (T b = 0.362, P = 0.007) were shown to have significant correlations (Additional file 1: Figure S1).

Fig. 3
figure 3

The percentage of samples that was found with specific gene mutations in genes with mutations in more than 5 % of the patients

Full size image
Fig. 4
figure 4

Correlations between gene mutations and cytogenetic changes. Significant correlations are marked with bold border

Full size image

Telomere length and genetic mutations/cytogenetic abnormalities

We analyzed whether somatic mutations had a relationship with TL. Table 3 shows the TL of patients according to whether or not they have a certain somatic mutation. Patients without gene mutations had shorter TL than those with mutation. All of the TL parameters showed significant differences (minimum, Q1, median, average, average of 0–10 percentile, SD, percentage of cells under the 10th percentile of the NC) between those patients with or without gene mutations (P = 0.043, 0.041, 0.025, 0.014, 0.030, 0.007, 0.017, respectively). Looking at specific mutations, the difference in TL was only present for those with or without CSF3R mutation. Patients without CSF3R mutation (n = 55) had a significantly higher percentage of cells whose length was below the lowest 10th percentile of the NC than those with CSF3R mutation (n = 3) (P = 0.037). In addition, comparisons were made for the TL according to gene groups. However, there were no significant differences in TL parameters by 10 different gene groups (Additional file 1: Table S4). Since TL varies according to age in the normal population, the age was compared between the mutational status of the genes or gene groups. Significant differences were only found for TET2 mutational status and the presence/absence of mutation in DNA methylation genes category. Patients with TET2 mutation and mutations of DNA methylation were older than those without mutation (P = 0.021, 0.008, respectively). The TL was not significantly different between patients with or without certain cytogenetic abnormalities by karyotyping or FISH (Table 3). In summary, we found that the presence of any gene mutation, but not the cytogenetic results, was related to TL.

Table 3 Telomere lengths according to gene mutations or cytogenetic changes
Full size table

Survival analysis

The OS was compared among different clinical subgroups (Fig. 5). There were significant differences in survival among different IPSS and the revised IPSS risk groups (P = 0.010, 0.004, respectively). However, no difference in OS was observed among the WHO categories (P = 0.462) by Kaplan-Meier survival analysis.

Fig. 5
figure 5

Kaplan-Meier survival curves according to WHO subtypes, IPSS subgroups, telomere lengths, and TP53 and U2AF1 mutation status

Full size image

The OS between those with or without certain gene mutations was compared in genes with mutations that were found in more than 5 % of the patients. We found that patients with TP53 mutation and U2AF1 mutation had a significantly shorter OS compared to those without mutation (median survival 32.6 vs. 6.0 months, P < 0.001 for TP53 mutation, 32.6 vs. 15.5, P = 0.035 for U2AF1 mutation). We also compared the OS of patients with ≥80 vs. <80 % of their cells with TL below the lowest 10th percentile TL of normal controls, and the former was significantly shorter (12.9 vs. 32.6 months, P = 0.017). The OS was significantly different in patients with ≥80 vs. <80 % of their cells with TL below the lowest 10th percentile TL of the NC even within WHO subtypes and IPSS subgroups (P = 0.038, 0.022, respectively). The presence of other mutations was not related to the OS.

Patients with complex karyotype (P = 0.015), del(7q)/−7 by karyotype (P < 0.001), abnormalities of 17p by karyotype (P = 0.001), abnormalities on FISH for chromosome 5 and 7 (P = 0.014, P < 0.001, respectively), high BM blast count (≥10 %) (P = 0.014), and those patients with ≥80 % of TL under the 10th percentile of the NC (P = 0.021) showed significantly lower survival than those without abnormalities. Univariable and multivariable Cox analyses for overall survival were performed with the clinical parameters, mutational profiles and TL parameters (Table 4). Excluding those parameters with significant correlation and many missing values, multivariable Cox analysis showed that the TP53 mutation, U2AF1 mutation, high BM blast count (≥10 %), and high percentage (≥80 %) of cells with TL less than the 10th percentile of the NC were independent prognostic factors for survival in patients with MDS (Table 4). The Harrell’s C-index was 0.649 (0.579–0.719) for the revised IPSS only and 0.717 (0.654–0.779) for the model including TL, TP53, and U2AF1 mutation and blast count showing similar prediction power for survival in our model including TL.

Table 4 Univariable and multivariable Cox analysis of the overall survival of 58 MDS patients
Full size table

Discussion

The findings of shortened TL in many MDS patients have led to the studies of telomere dynamics [9, 11, 38, 39] and its relationship with MDS [16]. The shortest TL in an individual, which is considered to be critical in cellular survival and telomere dysfunction [21, 40] was measured by Q-FISH in a single cell level in patients with MDS. We have found out that in addition to the mean TL, the minimum was also shorter in the MDS patients, and the patients showed a narrower range of TL than the NC. In average, 52.2 % of each patient had a TL of lower than the 10th percentile of the NC, suggesting a high burden of cells with critically short TL in the patients with MDS.

Few studies have analyzed the association of mutation status; only that of the telomerase complex-related genes and TL in MDS patients [17, 18]. Some patients with telomere attrition had TERC or TERT mutations, but only a small number of patients were included in each study [17]. Therefore, we evaluated TL in relation to mutation of myeloid neoplasia-related genes and genes commonly found in hematologic malignancies. Patients without gene mutations had shorter TL than those with mutation (regardless of the gene), but comparing the TL between those and those without each specific gene mutation, the difference was only present for CSF3R gene. In the present study, the population with TL below the 10th percentile TL of the NC was markedly low in patients with CSF3R mutation. It is well known that mutation of CSF3R is accompanied in congenital neutropenia. The congenital neutropenia patients that did not develop leukemia had an incidence of CSF3R mutation in 34 %, compared to 78 % for patients who showed transformation into acute leukemia [41]. Mutation of CSF3R could confer the proliferative activity to the hematopoietic stem cell, which possibly leads to a transformation to AML. This may be reflected in the TL of patients with CSF3R, having lesser portion of cells with short telomeres and one of our patients showed AML transformation shortly after diagnosis. The mutation site of CSF3R in MDS patients was different from the known site of mutation in congenital neutropenia but two of the patients showed CSF3R mutation reported in myeloid neoplasms such as chronic neutrophilic leukemia or chronic myelomonocytic leukemia [42].

The hypothesis of this study was that TL erosion might be more significant in certain groups of patients with gene mutations that are known be correlated with telomere dysfunction such as SRSF2 [16]. However, as TL is known to be affected by many other factors including age, environment, gender, stress, oxidative stress, and emotion [20, 43, 44], we were unable to show a correlation of TL in relation to a single gene mutation except for CSF3R. Contrary to our expectation, we found out that regardless of the gene mutation type, those patients with more than one gene mutation had higher TL than those without mutation. This finding somewhat contrasts with the recently reported result of aplastic anemia patients. Dumitriu et al. reported that telomere attrition and somatic mutations precede monosomy 7 and that the accumulation of short telomere in each chromosome may have a link to the development of aneuploidy in severe aplastic anemia (AA) [45]. This explanation can be applied to MDS evolved from previous AA. However, we infer that telomere attrition in AA exert a different influence on hematopoietic cells in case of MDS. AA has a decreased hematopoietic stem cell pool, but MDS has a dysfunctional hematopoietic stem cell pool, which has already acquired oncogenic mutations. Although the patients with gene mutations had longer TL than those without mutation in this study, the TL was still shorter than those of the NC. Comparisons within MDS groups may be less useful than the comparison of TL within an individual since telomere lengths are affected by many factors [20, 43, 44] including therapy, which was shown in a study that evaluated TL changes with treatment [46]. Another explanation may be that the mutations found in these patients were in the genes related to splicing machinery, DNA methylation, chromatin modification, or cohesion complexes and they are known to be mutated in MDS patients as driver/oncogenic mutations [14]. Therefore, these oncogenic mutations may have given the MDS cells the characteristics of “cancer” cells, which maintain telomere lengths [47], thus showing longer TL in patients with gene mutation. On the other hand, another explanation may be that these patients without gene mutations may have incorporated a mechanism that further shortens the length of the telomeres in this group. Mutations in other genes not included in this study such as TERC or TERT may give clue to these patients [17, 18]. Since no other study compared the TL of MDS patients and the presence of genetic mutation commonly found in hematologic malignancies, further studies would help elucidate the relationship between the telomere length and somatic mutations.

In addition, as TL is also regulated by the telomerase activity (TA) [1], we believe that the TA results would give additional information about the telomere dynamics. However, due to the sample limitation, TA could not be measured. TA was measured in MDS in few other studies, and MDS showed normal to low levels of telomerase activity [9, 48] despite the short TL. These results suggest the disruption of telomere maintenance in patients with MDS and suggest the maintenance of TL may be to some extent independent of TA incorporating alternative lengthening of telomere mechanisms, as is found in 10–15 % of the tumors [47].

The evaluation of TL with disease progression may be more informative and some reported telomere attrition during progression to acute myeloid leukemia [38, 49] but were shown in a small number of patients, thus requiring further study. In our study, there was no correlation of TL with various clinical subgroups and cytogenetic abnormalities, which is in line with previous studies showing inconsistent association between clinical factors and TL [9, 12]. Our finding may partly be due to the small number of patients included in our study, which is a major limitation of this study. Since cells with critically shortened telomeres may show repeated fusion bridge cycles, which leads to accumulation of chromosomal abnormalities, the relationship between chromosomal aberrations [50] and TL were also assessed. However, no significant change was found in relation to those with or without chromosomal aberrations. In addition to TP53 mutation, U2AF1 mutation was also suggested to be associated with a poor risk [51, 52] and was incorporated in the prognostic model along with the blast count and having high percentage of cells (≥80 %) with TL lower than the 10th percentile of the NC. Although the shortest telomere, the minimum TL, which is thought to be critical in the maintenance of cell survival, did not show prognostic significance in our study, we have found that having a high percentage of cells with short TL compared to NC may be of prognostic significance. The short TL was defined by having shorter than the lower 10th percentile value of the NC and this value may already suggest a portion of cells with critically short telomeres. Despite the limitation that Q-FISH measures the relative TL, the comparison made with the NC allowed us to confirm the shortened TL in MDS. Measuring TL in individual cells by Q-FISH also enabled us to identify the portion of cells with the short telomeres, which is not obtainable when measured by Southern blot or QPCR, which can only show the average TL. Moreover, our model including the TL parameter was comparable in the prediction of survival by Harrell’s C-index suggesting a probability of incorporation of TL and somatic mutation as a useful marker for prognostication.

Conclusions

In summary, we evaluated the telomere lengths of MDS patients at a single-cell level by quantitative FISH assessing the shortest TL in each individual along with the average and distribution of TLs. Patients with myelodysplastic syndrome showed short telomere lengths compared to the control with a significant proportion of cells having shorter telomeres than the 10th percentile of the normal control. A proportion of cell population with shorter TL below the 10th percentile of the normal control was significantly larger. We conclude that MDS is characterized by concentrated hematopoietic cells with short telomere, a kind of senescent cells. There were no specific somatic mutation or cytogenetic aberration correlating with short TL, but patients without somatic mutation harbored short telomere than those with somatic mutation. Interestingly, CSF3R mutation correlated with longer telomere length, and proportion of short telomere length was significantly small. TL parameter showed independent prognostic significance and incorporating somatic mutation and TL in addition to known prognostic markers may give additional prediction power for prognosis.

Abbreviations

AML, acute myeloid leukemia; BM, bone marrow; FISH, fluorescence in situ hybridization; IPSS, International Prognostic Scoring System; NC, normal control; MDS, myelodysplastic syndrome; OS, overall survival; SD, standard distribution; T/C, telomere to centromere; TL, telomere length; WHO, World Health Organization

References

  1. Greider CW. Telomere length regulation. Ann Rev Biochem. 1996;65(1):337–65.

    Article  CAS  PubMed  Google Scholar 

  2. McClintock B. The behavior in successive nuclear divisions of a chromosome broken at meiosis. Proc Natl Acad Sci U S A. 1939;25(8):405–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Capper R, Britt-Compton B, Tankimanova M, Rowson J, Letsolo B, Man S, et al. The nature of telomere fusion and a definition of the critical telomere length in human cells. Genes Dev. 2007;21(19):2495–508.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Bechter OE, Eisterer W, Pall G, Hilbe W, Kühr T, Thaler J. Telomere length and telomerase activity predict survival in patients with B cell chronic lymphocytic leukemia. Cancer Res. 1998;58(21):4918–22.

    CAS  PubMed  Google Scholar 

  5. Wu KD, Orme LM, Shaughnessy J, Jacobson J, Barlogie B, Moore MA. Telomerase and telomere length in multiple myeloma: correlations with disease heterogeneity, cytogenetic status, and overall survival. Blood. 2003;101(12):4982–9.

    Article  CAS  PubMed  Google Scholar 

  6. Aalbers AM, Calado RT, Young NS, Zwaan CM, Wu C, Kajigaya S, et al. Telomere length and telomerase complex mutations in pediatric acute myeloid leukemia. Leukemia. 2013;27(8):1786–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Hofmann WK, Koeffler HP. Myelodysplastic syndrome. Annu Rev Med. 2005;56:1–16.

    Article  CAS  PubMed  Google Scholar 

  8. Vardiman J, Arber D, Brunning R, Larson R, Matutes E, Baumann I. WHO classification of tumours of haematopoietic and lymphoid tissues. Lyon: International Agency for Research on Cancer; 2008.

    Google Scholar 

  9. Ohyashiki JH, Iwama H, Yahata N, Ando K, Hayashi S, Shay JW, et al. Telomere stability is frequently impaired in high-risk groups of patients with myelodysplastic syndromes. Clin Cancer Res. 1999;5(5):1155–60.

    CAS  PubMed  Google Scholar 

  10. Rollison DE, Epling-Burnette P, Park JY, Lee JH, Park H, Jonathan K, et al. Telomere length in myelodysplastic syndromes. Leuk Lymphoma. 2011;52(8):1528–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Boultwood J, Fidler C, Kusec R, Rack K, Elliott P, Atoyebi O, et al. Telomere length in myelodysplastic syndromes. Am J Hematol. 1997;56(4):266–71.

    Article  CAS  PubMed  Google Scholar 

  12. Lange K, Holm L, Vang Nielsen K, Hahn A, Hofmann W, Kreipe H, et al. Telomere shortening and chromosomal instability in myelodysplastic syndromes. Genes Chromosomes Cancer. 2010;49(3):260–9.

    CAS  PubMed  Google Scholar 

  13. Haferlach T, Nagata Y, Grossmann V, Okuno Y, Bacher U, Nagae G, et al. Landscape of genetic lesions in 944 patients with myelodysplastic syndromes. Leukemia. 2014;28(2):241–7.

    Article  CAS  PubMed  Google Scholar 

  14. Papaemmanuil E, Gerstung M, Malcovati L, Tauro S, Gundem G, Van Loo P, et al. Clinical and biological implications of driver mutations in myelodysplastic syndromes. Blood. 2013;122(22):3616–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Graubert TA, Shen D, Ding L, Okeyo-Owuor T, Lunn CL, Shao J, et al. Recurrent mutations in the U2AF1 splicing factor in myelodysplastic syndromes. Nat Genet. 2012;44(1):53–7.

    Article  CAS  Google Scholar 

  16. Colla S, Ong DST, Ogoti Y, Marchesini M, Mistry NA, Clise-Dwyer K, et al. Telomere dysfunction drives aberrant hematopoietic differentiation and myelodysplastic syndrome. Cancer Cell. 2015;27(5):644–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Yamaguchi H. Mutations of the human telomerase RNA gene (TERC) in aplastic anemia and myelodysplastic syndrome. Blood. 2003;102(3):916–8.

    Article  CAS  PubMed  Google Scholar 

  18. Kirwan M, Vulliamy T, Marrone A, Walne AJ, Beswick R, Hillmen P, et al. Defining the pathogenic role of telomerase mutations in myelodysplastic syndrome and acute myeloid leukemia. Hum Mut. 2009;30(11):1567–73.

    Article  CAS  PubMed  Google Scholar 

  19. Rigolin GM, Porta MD, Bugli AM, Castagnari B, Mauro E, Bragotti LZ, et al. Flow cytometric detection of accelerated telomere shortening in myelodysplastic syndromes: correlations with aetiological and clinical–biological findings*. Eur J Haematol. 2004;73(5):351–8.

    Article  PubMed  Google Scholar 

  20. Sanders JL, Newman AB. Telomere length in epidemiology: a biomarker of aging, age-related disease, both, or neither? Epidemiol Rev. 2013;35(1):112–31.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Hemann MT, Strong MA, Hao LY, Greider CW. The shortest telomere, not average telomere length, is critical for cell viability and chromosome stability. Cell. 2001;107(1):67–77.

    Article  CAS  PubMed  Google Scholar 

  22. Perner S, Brüderlein S, Hasel C, Waibel I, Holdenried A, Ciloglu N, et al. Quantifying telomere lengths of human individual chromosome arms by centromere-calibrated fluorescence in situ hybridization and digital imaging. Am J Pathol. 2003;163(5):1751–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Narath R, Lörch T, Greulich‐Bode KM, Boukamp P, Ambros PF. Automatic telomere length measurements in interphase nuclei by IQ‐FISH. Cytometry A. 2005;68(2):113–20.

    Article  PubMed  Google Scholar 

  24. Lansdorp PM, Verwoerd NP, Van De Rijke FM, Dragowska V, Little MT, Dirks RW, et al. Heterogeneity in telomere length of human chromosomes. Hum Mol Genet. 1996;5(5):685–91.

    Article  CAS  PubMed  Google Scholar 

  25. Czepulkowski B, Bhatt B, Rooney D. Basic techniques for the preparation and analysis of chromosomes from bone marrow and leukaemic blood. Human cytogenetics: malignancy and acquired abnormalities. 3rd ed. Oxford: Oxford University Press; 2001. p. 1–26.

    Google Scholar 

  26. Shaffer L, Slovak M, Campbell L. ISCN 2009. An international system for human cytogenetics nomenclature. 2009.

    Google Scholar 

  27. Shaffer LG, McGowan-Jordan, J, Schmid, M. ISCN 2013: an international system for human cytogenetic nomenclature (2013). Basel: Karger Medical and Scientific Publishers; 2013.

  28. Lee DS, Kim SH, Seo EJ, Park CJ, Chi HS, Ko EK, et al. Predominance of trisomy 1q in myelodysplastic syndromes in Korea: is there an ethnic difference? A 3-year multi-center study. Cancer Genet Cytogenet. 2002;132(2):97–101.

    Article  CAS  PubMed  Google Scholar 

  29. Walter M, Shen D, Shao J, Ding L, White B, Kandoth C, et al. Clonal diversity of recurrently mutated genes in myelodysplastic syndromes. Leukemia. 2013;27(6):1275–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Chapman MA, Lawrence MS, Keats JJ, Cibulskis K, Sougnez C, Schinzel AC, et al. Initial genome sequencing and analysis of multiple myeloma. Nature. 2011;471(7339):467–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Maxson JE, Gotlib J, Pollyea DA, Fleischman AG, Agarwa A, Eide CA, et al. Oncogenic CSF3R mutations in chronic neutrophilic leukemia and atypical CML. N Engl J Med. 2013;368(19):1781–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Wang L, Lawrence MS, Wan Y, Stojanov P, Sougnez C, Stevenson K, et al. SF3B1 and other novel cancer genes in chronic lymphocytic leukemia. N Engl J Med. 2011;365(26):2497–506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Landau DA, Carter SL, Stojanov P, McKenna A, Stevenson K, Lawrence MS, et al. Evolution and impact of subclonal mutations in chronic lymphocytic leukemia. Cell. 2013;152(4):714–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Doménech E, Gómez-López G, Gzlez-Peña D, López M, Herreros B, Menezes J, et al. New mutations in chronic lymphocytic leukemia identified by target enrichment and deep sequencing. Plos One. 2012;7(6):e38158.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Tenedini E, Bernardis I, Artusi V, Artuso L, Roncaglia E, Guglielmelli P, et al. Targeted cancer exome sequencing reveals recurrent mutations in myeloproliferative neoplasm. Leukemia. 2014;28:1052–9.

    Article  CAS  PubMed  Google Scholar 

  36. Greenberg PL, Tuechler H, Schanz J, Sanz G, Garcia-Manero G, Solé F, et al. Revised international prognostic scoring system for myelodysplastic syndromes. Blood. 2012;120(12):2454–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Greenberg P, Cox C, LeBeau MM, Fenaux P, Morel P, Sanz G, et al. International scoring system for evaluating prognosis in myelodysplastic syndromes. Blood. 1997;89(6):2079–88.

    CAS  PubMed  Google Scholar 

  38. Ohyashik K, Iwama H, Yahata N, Tauchi T, Kawakubo K, Shimamoto T, et al. Telomere dynamics in myelodysplastic syndromes and acute leukemic transformation. Leuk Lymphoma. 2001;42(3):291–9.

    Article  Google Scholar 

  39. Sieglová Z, Žilovcová S, Čermák J, Řı́hová H, Březinová D, Dvořáková R, et al. Dynamics of telomere erosion and its association with genome instability in myelodysplastic syndromes (MDS) and acute myelogenous leukemia arising from MDS: a marker of disease prognosis? Leuk Res. 2004;28(10):1013–21.

    Article  PubMed  Google Scholar 

  40. Abdallah P, Luciano P, Runge KW, Lisby M, Géli V, Gilson E, et al. A two-step model for senescence triggered by a single critically short telomere. Nat Cell Biol. 2009;11(8):988–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Zeidler C, Germeshausen M, Klein C, Welte K. Clinical implications of ELA2, HAX1- and G-CSF-receptor (CSF3R) mutations in severe congenital neutropenia. Brit J Haematol. 2009;144:459–67.

    Article  CAS  Google Scholar 

  42. Germeshausen M, Ballmaier M, Welte K. Incidence of CSF3R mutations in severe congenital neutropenia and relevance for leukemogenesis: results of a long-term survey. Blood. 2007;109(1):93–9.

    Article  CAS  PubMed  Google Scholar 

  43. Gardner M, Bann D, Wiley L, Cooper R, Hardy R, Nitsch D, et al. Gender and telomere length: systematic review and meta-analysis. Exp Gerontol. 2014;51:15–27.

    Article  CAS  PubMed  Google Scholar 

  44. Puterman E, Lin J, Blackburn E, O’Donovan A, Adler N, Epe E. The power of exercise: buffering the effect of chronic stress on telomere length. PLoS One. 2010;5(5):e10837.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Dumitriu B, Feng X, Townsley DM, Ueda Y, Yoshizato T, Calado RT, et al. Telomere attrition and candidate gene mutations preceding monosomy 7 in aplastic anemia. Blood. 2015;125(4):706–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Beier F, Masouleh BK, Buesche G, Ferreira MSV, Schneider RK, Ziegler P, et al. Telomere dynamics in patients with del (5q) MDS before and under treatment with lenalidomide. Leuk Res. 2015;39(11):1292–8.

    Article  CAS  Google Scholar 

  47. Jafri MA, Ansari SA, Alqahtani MH, Shay JW. Roles of telomeres and telomerase in cancer, and advances in telomerase-targeted therapies. Genome Med. 2016;8(1):69.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Kim M, Oh B, Kim TY, Yoon SS, Kim SY, Hwang SM, Lee DS. Elevated telomerase activity in essential thrombocythemia with extreme thrombocytosis. Clin Biochem. 2014;47(6):389–92.

    Article  CAS  PubMed  Google Scholar 

  49. Göhring G, Lange K, Hofmann W, Nielsen K, Hellström-Lindberg E, Roy L, et al. Telomere shortening, clonal evolution and disease progression in myelodysplastic syndrome patients with 5q deletion treated with lenalidomide. Leukemia. 2012;26(2):356–8.

    Article  PubMed  Google Scholar 

  50. Artandi SE, Chang S, Lee S-L, Alson S, Gottlieb GJ, Chin L, et al. Telomere dysfunction promotes non-reciprocal translocations and epithelial cancers in mice. Nature. 2000;406(6796):641–5.

    Article  CAS  PubMed  Google Scholar 

  51. Wu SJ, Tang JL, Lin CT, Kuo YY, Li LY, Tseng MH, et al. Clinical implications of U2AF1 mutation in patients with myelodysplastic syndrome and its stability during disease progression. Am J Hematol. 2013;88(11):E277–82.

    Article  CAS  PubMed  Google Scholar 

  52. Bejar R, Stevenson KE, Caughey RA, Abdel-Wahab O, Steensma DP, Galili N, et al. Validation of a prognostic model and the impact of mutations in patients with lower-risk myelodysplastic syndromes. J Clin Oncol. 2012;30:3376–82.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank the Division of Statistics in Medical Research Collaborating Center at Seoul National University Bundang Hospital for the statistical assistance.

Funding

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (NRF-2014R1A2A1A10052286) and the Ministry of Education (NRF-2015R1D1A1A02062419) of South Korea. The funding sources played no role in the design; in the collection, analysis, and interpretation of the data; in the writing of the manuscript; and in the decision to submit the manuscript for publication.

Availability of data and materials

Although all our data is de-identified, we opt not to share the data and materials in public due to further study on this subject. However, we will share the data in request by other researchers if necessary. All of the methods including the software programs or reagents used in this study are on market, which are accessible by other researchers.

Authors’ contributions

SMH, SYK, KOI, SNP, JAK, HSP, SMK, and KK performed the research. JAK, SMH, SYK, and DSL designed the research study. SNP, KOI, KK, SMK, and SMH analyzed the data. SMH, SYK, and DSL wrote the paper. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

This study was approved by the institutional review board of Seoul National University Hospital (1604-082-754), and the participants gave informed consent.

Author information

Authors and Affiliations

  1. Department of Laboratory Medicine, Seoul National University Bundang Hospital, Seongnam, Republic of Korea

    Sang Mee Hwang

  2. Department of Laboratory Medicine, Seoul National University College of Medicine, Seoul, Republic of Korea

    Sang Mee Hwang, Jung Ah Kim, Hee-Sue Park & Dong Soon Lee

  3. Department of Laboratory Medicine, Chungnam National University College of Medicine, Daejeon, Republic of Korea

    Seon Young Kim

  4. Cancer Research Institute, Seoul National University College of Medicine, Seoul, Republic of Korea

    Si Nae Park, Kyongok Im, Kwantae Kim, Sung-Min Kim & Dong Soon Lee

Authors
  1. Sang Mee Hwang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Seon Young Kim
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jung Ah Kim
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Hee-Sue Park
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Si Nae Park
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Kyongok Im
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Kwantae Kim
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Sung-Min Kim
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Dong Soon Lee
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Dong Soon Lee.

Additional file

Additional file 1:

Supplementary method. Figure S1. Correlation of gene categories and cytogenetic abnormalities. Table S1. Genes included in the target sequencing and its gene categories. Table S2. Telomere lengths according to WHO categories or IPSS subgroups. Table S3. Telomere lengths according to laboratory results. Table S4. Telomere lengths according to mutations in gene groups. (DOC 803 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hwang, S.M., Kim, S.Y., Kim, J.A. et al. Short telomere length and its correlation with gene mutations in myelodysplastic syndrome. J Hematol Oncol 9, 62 (2016). https://doi.org/10.1186/s13045-016-0287-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13045-016-0287-9

Keywords

Journal of Hematology & Oncology

ISSN: 1756-8722

Contact us