Dysregulation of H/ACA ribonucleoprotein components in chronic lymphocytic leukemia

Abstract

Telomeres are protective repeats of TTAGGG sequences located at the end of human chromosomes. They are essential to maintain chromosomal integrity and genome stability. Telomerase is a ribonucleoprotein complex containing an internal RNA template (hTR) and a catalytic subunit (hTERT). The human hTR gene consists of three major domains; among them the H/ACA domain is essential for telomere biogenesis. H/ACA ribonucleoprotein (RNP) complex is composed of four evolutionary conserved proteins, including dyskerin (encoded by DKC1 gene), NOP10, NHP2 and GAR1. In this study, we have evaluated the expression profile of the H/ACA RNP complex genes: DKC1, NOP10, NHP2 and GAR1, as well as hTERT and hTR mRNA levels, in patients with chronic lymphocytic leukemia (CLL). Results were correlated with the number and type of genetic alteration detected by conventional cytogenetics and FISH (fluorescence in situ hybridization), IGHV (immunoglobulin heavy chain variable region) mutational status, telomere length (TL) and clinico-pathological characteristics of patients. Our results showed significant decreased expression of GAR1, NOP10, DKC1 and hTR, as well as increased mRNA levels of hTERT in patients compared to controls (p≤0.04). A positive correlation between the expression of GAR1-NHP2, GAR1-NOP10, and NOP10-NHP2 (p<0.0001), were observed. The analysis taking into account prognostic factors showed a significant increased expression of hTERT gene in unmutated-IGHV cases compared to mutated-CLL patients (p = 0.0185). The comparisons among FISH groups exhibited increased expression of DKC1 in cases with two or more alterations with respect to no abnormalities, trisomy 12 and del13q14, and of NHP2 and NOP10 compared to those with del13q14 (p = 0.03). The analysis according to TL showed a significant increased expression of hTERT (p = 0.0074) and DKC1 (p = 0.0036) in patients with short telomeres compared to those with long TL. No association between gene expression and clinical parameters was found. Our results suggest a role for these telomere associated genes in genomic instability and telomere dysfunction in CLL.

About the Authors:

Patricia Carolina Dos Santos

Affiliation: Laboratorio de Genética de Neoplasias Linfoides, Instituto de Medicina Experimental, CONICET-Academia Nacional de Medicina, Buenos Aires, Argentina

ORCID http://orcid.org/0000-0002-6607-7620

Julieta Panero

Affiliation: Laboratorio de Genética de Neoplasias Linfoides, Instituto de Medicina Experimental, CONICET-Academia Nacional de Medicina, Buenos Aires, Argentina

Carmen Stanganelli

Affiliation: División Patología Molecular, Instituto de Investigaciones Hematológicas-Academia Nacional de Medicina, Buenos Aires, Argentina

Virginia Palau Nagore

Affiliation: Laboratorio de Genética de Neoplasias Linfoides, Instituto de Medicina Experimental, CONICET-Academia Nacional de Medicina, Buenos Aires, Argentina

Flavia Stella

Affiliation: Laboratorio de Genética de Neoplasias Linfoides, Instituto de Medicina Experimental, CONICET-Academia Nacional de Medicina, Buenos Aires, Argentina

Raimundo Bezares

Affiliation: Servicio de Hematología, Hospital Teodoro Álvarez, Buenos Aires, Argentina

Irma Slavutsky

* E-mail: [email protected]

Affiliation: Laboratorio de Genética de Neoplasias Linfoides, Instituto de Medicina Experimental, CONICET-Academia Nacional de Medicina, Buenos Aires, Argentina

Introduction

Telomeres are protective repeats of TTAGGG sequences located at the end of chromosomes [1]. One important function of the telomere is to prevent end-to-end joining of chromosomes, which otherwise could lead to the accumulation of chromosomal breaks during mitosis. Due to the end-replication problem, telomeres become progressively shorter with repeated cell division. Telomere shortening leads to telomere dysfunction and elicitation of a DNA damage response. Moreover, it has implications in oncogenesis and cellular aging due to limited cellular proliferation below a certain critical telomere length (TL) [2].

TL is maintained by the enzymatic complex telomerase. This enzyme participates in the formation and maintenance of telomeres in normal progenitor cells [3] and is active in most of tumors and immortal cell lines, endowing cancer cells with an infinite replicative capacity [4, 5]. Telomerase consists of two core components: a catalytic subunit, human telomerase reverse transcriptase (hTERT) with reverse transcriptase activity and a RNA component, human telomerase RNA component (hTR) that serves as a template for the elongation of telomeres. As described in detail previously [6] human hTR extends for 451 nucleotides, the last 240 at the 3´end form a consensus H/ACA two hairpin structure while the 5' half folds into a pseudoknot containing the template for the reverse transcriptase. While it is not needed for telomerase activity in vitro, the H/ACA domain of hTR is required for its in vivo accumulation and stability [6, 7]. Telomerase activity requires additional ribonucleoprotein (RNP) factors that associate with hTERT and hTR to form the holoenzyme. Thus, H/ACA RNPs consist of four evolutionarily conserved proteins, Dyskerin [encoded by the gene DKC1 (Dyskerin Pseudouridine Synthase 1)], NHP2 (NHP2 ribonucleoprotein), NOP10 (NP10 ribonucleoprotein), and GAR1 (GAR1 ribonucleoprotein), and a function-specifying, noncoding H/ACA RNAs [8, 9]. All these RNPs are concentrated in nucleoli and Cajal bodies of mammalian cells, reflecting the location of H/ACA RNPs [10]. In reference to its functions, H/ACA RNPs contribute to telomerase assembly and stabilization, posttranscriptional processing of nascent ribosomal RNA and pre-mRNA splicing [6]. Dyskerin, NHP2 and NOP10 are essential for the cellular accumulation of the H/ACA-motif RNAs, while GAR1 binds only to Dyskerin and is necessary for the nucleolar localization and function of the RNP complex [11, 12]. Particularly, Dyskerin functions as a pseudouridine synthase and is involved in post-transcriptional modifications of rRNA [13]. In addition, Dyskerin is part of the telomerase complex and is required for the correct activity of the enzyme, which directly binds to and stabilizes hTR within the complex [14], playing a role in the maintenance of telomere integrity. Point mutations in DKC1 gene cause the X-linked form of dyskeratosis congenital, a disease with increased predisposition to cancer.

Chronic lymphocytic leukemia (CLL) is the most frequent leukemia in adults in Western countries; it is a heterogeneous disease with variable clinical presentation and evolution [15]. At the molecular level, two major subtypes can be distinguished, mutated (M) and unmutated (UM), characterized respectively by a high or low number of somatic hypermutations in the variable region of immunoglobulin genes and different clinical evolution [16, 17]. Cytogenetic and FISH (fluorescence in situ hybridization) studies have also proved to be important tools in the biologic characterization of this disease, allowing the identification of distinct cytogenetic risk groups and providing information for clinical outcome [18, 19]. In addition, short telomeres have been associated with genetic complexity and a high risk of genomic aberrations [20–23], and ultimately implicated in disease outcome [21, 22, 24]. Increased telomerase expression and activity were also observed in CLL patients [25, 26]. However, data concerning the participation of H/ACA RNPs in telomere dysfunction in CLL are lacking. In the present study, we investigated the expression profile of DKC1, NHP2, NOP10 and GAR1, in association with hTR and hTERT mRNA levels. Results were correlated with the number and type of genetic alterations detected by conventional cytogenetics and FISH, the mutational status of IGHV (immunoglobulin heavy chain variable region), TL and clinico-pathological characteristics of patients.

Materials and methods

Patients

This study was based on a retrospective series of 71 CLL patients at diagnosis (31 females; mean age 67 years, range: 44–88 years) derived to our Laboratory from 2006 to 2016, and from whom RNA samples were available. Patients were diagnosed according to the National Cancer Institute-Working Group CLL criteria [27]. Rai stages [28] were available for 54 cases with the following distribution: 0: 19 (35.2%); I–II: 30 (55.5%); and III–IV: 5 (9.3%). Table 1 shows clinical and biological characteristics of patients. In addition, 21 normal controls (NC) (without a medical history of cancer) (10 females; mean age 62 years, range: 40–89 years), were properly selected to this study. The study was approved by the Ethics Committee of the Institutes of the National Academy of Medicine of Buenos Aires (Approval number: T.I.N° 12042/15/X). All patients provided their written informed consent.

[Figure omitted. See PDF.]

Table 1. Clinical and biological characteristics of CLL patients.

https://doi.org/10.1371/journal.pone.0179883.t001

RNA extraction and reverse transcription (RT)

Total RNA was extracted from mononuclear cells isolated on a Ficoll-Paque Plus (GE Healthcare Bio-Sciences) density gradient of peripheral blood (PB) samples of patients, normal controls (NC) and the K-562 cell line, using the Trizol reagent. K-562 cell line was obtained from the Laboratory of Genetics in Hematology, Institute of Experimental Medicine, CONICET-National Academy of Medicine (Buenos Aires, Argentina). The complementary DNA (cDNA) synthesis was carried out by using Moloney Murine Leukemia Virus reverse transcriptase (MMLV) and random primer (Promega). The cDNA synthesis was performed in a final volume of 20 μl, containing 1μg of the total RNA, for 10 minutes at 70°C, for 60 minutes at 37°C and 10 minutes at 95°C to inactive the enzyme. cDNA was stored at -20°C until use.

Gene expression analysis by quantitative real-time RT-PCR

The analysis of four genes of the RNP complex GAR1, NHP2, NOP10 and DKC1, as well as both telomerase units, hTR and hTERT, were performed using real-time quantitative PCR (qPCR) on a Rotor GENE Q (Qiagen), based on SybrGreen methodology. Primer sequences for hTR and hTERT were previously described [29, 30]. For the remaining genes, primers were designed for this work and detailed in Table 2. All PCR runs were performed in duplicate, using 2.5 μl of RT reaction, 2X SYBR^® Select Master Mix (Applied Biosystems^®), 10 mM of each primer, in a 20 μl final volume.

[Figure omitted. See PDF.]

Table 2. Primer sequences.

https://doi.org/10.1371/journal.pone.0179883.t002

The constitutive expression gene β-Actin was used to normalize sample-to-sample differences in cDNA input, RNA quality and RT efficiency [31]. For all targets, the thermal cycling conditions were 2 minutes at 50°C, followed by 2 minutes at 95°C; continued by 45 cycles at 95°C for 20 seconds and 60°C for 1 minute, ending with a melting curve from 50°C to 99°C. All measurements included a determination of the standards and no-template as a negative control, in which water was substituted for the cDNA. Standard curves were constructed with fivefold serial dilutions of the cDNA from the K-562 cell line.

Cytogenetics and FISH analysis

Chromosome analyses were performed on PB lymphocytes, cultured for 96 h at 37°C in F-12 medium supplemented with 15% of fetal calf serum, stimulated with Pokeweed mitogen. Slides were prepared by conventional method. G-banding technique was used. For FISH analysis, slides were hybridized with SE 12, OLE13q14 D13S319, OLE11q22.3 ATM, and 17p13.1 TP53 DNA probes (LiVE-LEXEL, Buenos Aires, Argentina) according to the manufacturer’s protocol. Two hundred interphase nuclei were analyzed for each probe. The cut-offs for positive values (mean of normal control plus 3 standard deviations) determined from ten cytogenetically normal donors were as follows: 3%, 10%, 7.5%, and 5.5% for trisomy 12, monosomies of D13S319, ATM, and TP53, respectively.

Absolute telomere length analysis

High-molecular weight genomic DNA was obtained from PB mononuclear cells of CLL patients and controls, separated on a Ficoll-Paque Plus (GE Healthcare Bio-Sciences, Uppsala, Sweden) density gradient. Absolute TL measurement was carried out by real-time quantitative PCR (qPCR) in a LightCycler system (Roche Diagnostics) according to Panero et al. [32]. For each DNA sample, two consecutive reactions were performed: the first one to amplify a single copy gene, RPLP0 (ribosomal protein lateral stalk subunit P0) (12q24.2), and the second one for telomeric sequence. Briefly, both PCRs were performed in a final volume of 20 μl using 1X SYBRGreen Master Mix (Roche), 100 nmol/L of primers, and 20 ng of DNA from patients and controls. The PCR conditions were as follows: 95°C 10 min, followed by 45 cycles of 95°C 15 s, 60°C 1 min. The melting curve was performed with 1 cycle of 95°C 20 s, 50°C 15 s, and 98°C with a temperature ramp of 0.1°C/s. All samples were analyzed in duplicate. Absolute TL (kb/diploid genome) values were calculated as previously described [33]. The conversion of absolute TL values into the equivalent Southern blot TL values (kb) was done by using the regression formula: y = 0.0746x + 0.5285, where x = Log [TL (qPCR)] and y = Log [TL (TRF)] as was previously explained [243].

IGHV mutational status

The IGHV gene sequences were determined as previously described [34]. Briefly, amplification of IGHV regions by polymerase chain reaction was performed on cDNA by using VH framework region 1 consensus family specific primers (VH1-VH6) and JH primers. When amplifications of these primers were unsuccessful, an alternative set of primers that anneal to sequences in the leader region (LH1-LH6) and one antisense Cμ-primer were used. Thermal cycling conditions were 3 minutes at 93°C, followed by 33 cycles at 94°C for 30 seconds, 62°C for 30 seconds, 72°C for 30 seconds, elongation at 72°C for 7 minutes, and a final step at 4°C for 10 minutes. Polymerase chain reaction products were purified in 2% agarose gels, sequenced bidirectionally, and analyzed on an automated DNA sequence analyzer (377 ABI Prism, PE biosystem, Foster City, CA). Sequence data were analyzed by using IgBLAST (immunoglobulin BLAST) (http://www.ncbi.nlm.nih.gov/igblast) and the ImMunoGeneTics database (IMGT) (http://imgt.cines.fr). IGHV sequences with <98% homology with respect to the germ line counterpart were considered as M, whereas those with homology of 98% or higher were classified as UM [17, 18].

Statistical analysis

Statistical analysis was performed using GraphPad Prism Version 5.0 (2008). Comparisons of data from patients and controls were performed using the Mann-Whitney test meanwhile Kruskal-Wallis test was used for comparisons among subgroups. Groupwise comparison of the distribution of clinical and laboratory variables was performed with the Student t test (for quantitative variables) and the χ2 or Fisher’s exact test (for categorical variables). Correlations between gene expression and TL or clinical variables were assessed by using the Kendall’s coefficient. The cut-off point for gene expression was selected according to receiver operating characteristic (ROC) analysis. A hierarchical clustering was employed to segregate patients groups based on the qPCR expression patterns using the Heatmapper (www.heatmapper.ca). Treatment free survival (TFS), calculated from the date of diagnosis to the first CLL-specific treatment were estimated by the Kaplan-Meier method and compared with the Log-rank test. For all tests, p<0.05 was considered statistically significant.

Results

Expression profile of RNP complex genes and telomerase subunits

The expression profile of the four RNP complex genes: GAR1, NHP2, NOP10 and DKC1, as well as the two subunits of telomerase, hTR and hTERT, were evaluated in 71 CLL patients and 21 NC. The analysis showed a decreased expression of RNP complex genes in CLL patients compared to NC, with significant differences for: GAR1 (p = 0.0156), DKC1 (p = 0.0025) and NOP10 (p = 0.0048). No statistical differences were found for NHP2 mRNA expression (p = 0.182) (Fig 1A). In reference to telomerase subunits, an increase of hTERT (p = 0.0441) and a decrease of hTR (p = 0.0005) in patients compared to NC were observed (Fig 1B). In addition, a positive correlation between the expression of GAR1-NHP2, GAR1-NOP10, and NOP10-NHP2 (p<0.0001) were found, indicating a strong interaction among them (Supporting information; S1 Fig). A more integrative vision of telomere-associated gene expression profiles can be observed in Supporting information (Heat map; S2 Fig).

[Figure omitted. See PDF.]

Fig 1. Gene expression profiles in CLL patients and normal controls (NC).

(A) GAR1, NHP2, DKC1 and NOP10 genes. (B) hTR and hTERT genes.

https://doi.org/10.1371/journal.pone.0179883.g001

Correlation between expression profiles and prognostic factors in CLL

We would like to determine if the expression profile of these genes correlated with prognostic factors of relevance in CLL. Thus, we analyzed the association between transcription levels and different parameters, including the number and type of genetic alterations detected by conventional cytogenetics and FISH, the mutational status of IGHV, TL and clinico-pathological characteristics of patients.

IGHV mutational status was evaluated in 66 cases: 34 showed M-CLL, 31 were UM-CLL and one case had a stop codon. Expression analysis of telomerase genes according to this parameter revealed significant increased expression of hTERT gene in UM-CLL patients (p = 0.0185) compared to M-CLL (Fig 2).

[Figure omitted. See PDF.]

Fig 2. Increased expression of hTERT gene in unmutated (UM)-CLL patients compared to mutated (M)-CLL.

https://doi.org/10.1371/journal.pone.0179883.g002

Cytogenetic and/or FISH analysis were performed in 56 patients. Eight cases had two or more cytogenetic alterations in the karyotype, 14 showed normal karyotype and no FISH alterations, 14 had del13q14 as the only abnormality, 7 cases showed trisomy 12, and the remaining 13 had two or more FISH alterations. No differences in gene expression profiles between patients with two or more alterations by cytogenetics and FISH analysis were found, thus they were considered as a single group (2 or more CA). The analysis among groups showed significant differences for DKC1 (p = 0.0028) and hTR (p = 0.0129) and a tendency for NHP2 (p = 0.057) and NOP10 (p = 0.068). Interestingly, when risk genetic groups were compared, an increased expression of DKC1 in cases with 2 or more CA with respect to those with no abnormalities (NA) (p = 0.0120), trisomy 12 (p = 0.0362) and del13q14 (p = 0.0304) (Fig 3A), was observed. Moreover, overexpression of NHP2 and NOP10 genes in patients with 2 or more CA compared to those with del13q14 (p = 0.003 and p = 0.0375, respectively) were also found (Fig 3B and 3C).

[Figure omitted. See PDF.]

Fig 3. Analysis of gene expression profiles of telomere associated genes taking into account the number and type of genetic alterations in CLL patients.

(A) Significant differences for DKC1 in cases with two or more chromosome alterations (2 or more CA) with respect to no abnormalities (NA), trisomy 12 and del13q14. In addition, significant differences between normal controls (NC) and CLL patients with NA (p = 0.01134), trisomy 12 (p = 0.0274) and del13q14 (p = 0.0081) were observed; (B) Increased expression of NHP2 in patients with 2 or more CA compared to those with del13q14; (C) Increased expression of NOP10 in patients with 2 or more CA compared to those with del13q14. Furthermore, significant differences between NC and del13q14 (p = 0.0115) was observed.

https://doi.org/10.1371/journal.pone.0179883.g003

TL was measured in 32 patients of our cohort. Short telomeres were defined as below the tenth percentile of normal controls (6.63 Kb). Twelve cases had short TL and 20 showed long TL. This analysis showed a significant increased expression of hTERT (p = 0.0074) and DKC1 (p = 0.0036) in patients with short telomeres compared to those with long TL. In contrast, the remaining genes did not showed significant differences (Fig 4).

[Figure omitted. See PDF.]

Fig 4. Analysis of hTERT and DKC1 gene expression profiles tacking into account telomere length (TL).

Significant increased expression for hTERT and DKC1 in patients with short TL compared to those with long TL.

https://doi.org/10.1371/journal.pone.0179883.g004

Finally, the analysis taking into account clinico-pathological characteristics of patients as well as treatment-free survival and overall survival did not reveal any prognostic relevance (Supporting information; S1 Table).

Discussion

In the present study, we have analyzed the expression profiles of four genes of the RNP complex GAR1, NHP2, NOP10 and DKC1, as well as both telomerase units, hTR and hTERT, in CLL patients at diagnosis. It is important to pointing out that, to our knowledge, this is the first analysis of GAR1, NHP2, NOP10 and hTR in this pathology.

Our results showed significant decreased expression GAR1, NOP10, DKC1 and hTR, as well as increased mRNA levels of hTERT in CLL patients compared to controls, supporting the importance of changes in the expression of telomere-associated genes in this entity. There is scarce information about these genes in hematological malignancies. A report of our group [35] using microarray assay found increased expression of NOP10, GAR1 and NHP2 in myeloma cell lines and patients with plasma cell disorders, showing significant overexpression from monoclonal gammopathy of undetermined significance (MGUS) to plasma cell leukemia, indicating their association with tumor progression. On the contrary, there are some reports about DKC1 and hTERT expression in CLL. In concordance with our findings, two studies found decreased expression of DKC1 compared to controls [36, 37]. Sbarrato et al [37] suggested that the low DKC1 expression leads to an imbalance in ribosomal proteins that could influence the response of leukemic cells to the microenvironment, and proposed a general impact of DKC1 dysregulation on the translational machinery in CLL. Other authors suggest that DKC1 deficiency may contribute to tumorigenesis by altering the splicing of specific mRNAs, or by modulating the level of certain snoRNAs [13]. In vitro experiments showed that the loss of DKC1 function affects telomerase activity by reducing hTR levels, leading to premature telomere shortening that may result in chromosomal end-to-end fusions, breakage and rearrangements associated to tumor development [13].

In reference to hTERT expression, there are discordant results probably related with the number of patients and/or the clinical characteristics of the cohorts evaluated. Poncet et al [36] found hTERT down-regulation and Damle et al [38] did not observed significant differences in CLL patients compared to controls. On the contrary, Hoxha et al [3926] observed hTERT overexpression by microarray analysis, data that could not be validated by qRT-PCR. Interestingly, we found increased hTERT mRNA levels associated to UM-CLL and short TL. Conflicting results were reported about the association between hTERT mRNA expression and IGHV mutational status [21, 26, 38–40]. Our data, in concordance with those observed by Rampazzo et al [22], support the relationship between short TL, UM-IGHV status and increased hTERT expression in this pathology, as well as the importance of hTERT levels to maintain the replicative potential of tumor cells. In addition, a previous study of our group in patients with multiple myeloma (MM) also showed an inverse correlation between hTERT levels and TL, in which shortest telomeres had the highest hTERT expression [41]. Furthermore, it is important to pointing out that, besides maintaining the TL, hTERT is also involved in other cellular functions like cell survival and prevention of apoptosis [42, 43], which would be important to evaluate in CLL.

With regards to genomic aberrations, our data showed significantly increased mRNA levels of DKC1, NHP2 and NPO10 genes in CLL cases with two or more alterations. Using microarray assay, von Stedingk et al [44] found DKC1, NHP2 and GAR1 overexpression in high-stage neuroblastoma associated to poor prognosis and genomic complexity, as well as telomere dysfunction, proposing their relationship with tumor aggressiveness. More studies will be necessary to clarify the role of these genes in CLL progression. Simultaneously, we also observed increased DKC1 expression in patients with short TL, supporting the probable involvement of this gene in progressive telomere reduction, promoting genomic instability and immortalization of cancer cells. In concordance, a previous report [45] in lung cancer cells also found a strong association between Dyskerin and short telomeres and more recently our group showed DKC1 overexpression in MM patients with short TL [46]. In reference to GAR1, we did not find association with the different prognostic factors evaluated. A recent report on cell lines [47] found differential expression of this protein in response to different genotoxic agents and particularly to DNA damage response. The authors propose that this overexpression may be responsible, at least in part, to the survival and proliferations of tumor cells, suggesting that this protein may function independently of its role within the H/ACA RNP complex.

Finally, we were interested in correlating our data with clinical characteristics of patients as well as the evolution of the disease. This analysis did not reveal any prognostic relevance for the genes evaluated. The literature shows discordant results for hTERT expression, with authors observing strong impact of telomerase activity on overall survival [22, 25, 39] while others did not find any association [221, 276]. In addition, a recent report [37] did not find impact of DKC1 expression on overall survival or progression free survival, but they observed that cases with low DKC1 levels had a reduced survival following chemotherapy. There are no reports about this matter for the remaining genes.

The current study is subject to a number of limitations. Our analysis was performed in a retrospective cohort of CLL patients. Thus, as regards prognostic factors and outcome, the lack in statistical significance could be due to the limited number of patients in our series. Additionally, we have not investigated protein expression, results that would increase the significance of this study. Future analysis in larger cohorts will be useful to better define the prognostic value of these genes as well as functional modifications.

Concluding, we have observed modifications in the expression profiles of GAR1, NOP10, DKC1, hTERT and hTR genes in CLL patients. The correlation with prognostic factors of the disease showed associations among DKC1, NHP2 and NPO10 mRNA levels with the presence of two or more genetic alterations and of hTERT and DKC1 expression with short telomere length, suggesting a role for these telomere associated genes in genomic instability and telomere dysfunction in this pathology.

Supporting information

[Figure omitted. See PDF.]

S1 Fig. Correlation among gene expression profiles of NHP2, GAR1 and NOP10.

https://doi.org/10.1371/journal.pone.0179883.s001

(TIF)

S2 Fig. Heat map.

Telomere-associated gene expression profiles in CLL patients.

https://doi.org/10.1371/journal.pone.0179883.s002

(TIF)

S1 Table. mRNA expression levels of telomere-associated genes and clinical and laboratory features.

https://doi.org/10.1371/journal.pone.0179883.s003

(DOCX)

Author Contributions

1. Conceptualization: IS JP PCDS.

2. Data curation: PCDS VPN FS CS.

3. Formal analysis: PCDS CS JP.

4. Funding acquisition: IS.

5. Investigation: PCDS VPN FS.

6. Methodology: PCDS JP IS.

7. Project administration: IS.

8. Resources: RB IS.

9. Supervision: IS.

10. Visualization: PCDS JP IS.

11. Writing – original draft: PDS JP IS.

12. Writing – review & editing: PDS JP IS.

Citation: Dos Santos PC, Panero J, Stanganelli C, Palau Nagore V, Stella F, Bezares R, et al. (2017) Dysregulation of H/ACA ribonucleoprotein components in chronic lymphocytic leukemia. PLoS ONE 12(6): e0179883. https://doi.org/10.1371/journal.pone.0179883

References

1. de Lange T. Shelterin: the protein complex that shapes and safeguards human telomeres. Genes Dev 2005; 19: 2100–2110. pmid:16166375

2. Blasco MA. Telomeres and human disease: ageing, cancer and beyond, Nat. Rev.Genet. 2005; 6: 611–622. pmid:16136653

3. Ulaner GA. Telomere maintenance in clinical medicine. Am J Med. 2004; 117: 262–269. pmid:15308436

4. Collins K, Mitchell JR. Telomerase in the human organism. Oncogene 2002; 21:564–579. pmid:11850781

5. Forsyth NR, Wright WE, Shay JW. Telomerase and differentiation in multicellular organisms: turn it off, turn it on, and turn it off again. Differentiation 2002; 69:188–197. pmid:11841477

6. Meier UT. The many facets of H/ACA ribonucleoproteins. Chromosoma 2005; 114:1–14. pmid:15770508

7. Bachand F, Triki I, Autexier C. Human telomerase RNA-protein interactions. Nucleic Acids Res 2001; 29: 3385–3393. pmid:11504876

8. Meier UT. How a single protein complex accommodates many different H/ACA RNAs. Trends Biochem Sci 2006; 31:311–315. pmid:16647858

9. Egan ED, Collins K. An enhanced H/ACA RNP assembly mechanism for human telomerase RNA. Mol Cell Biol. 2012; 32: 2428–2439. pmid:22527283

10. Pogacić V, Dragon F, Filipowicz W. Human H/ACA small nucleolar RNPs and telomerase share evolutionarily conserved proteins NHP2 and NOP10. Mol Cell Biol. 2000; 20: 9028–9040. pmid:11074001

11. Collins K. The biogenesis and regulation of telomerase holoenzymes. Nat Rev Mol Cell Biol 2006; 7: 484–494. pmid:16829980

12. Darzacq X, Kittur N, Roy S, Shav-Tal Y, Singer RH, Meier UT. Stepwise RNP assembly at the site of H/ACA RNA transcription in human cells. J Cell Biol 2006; 173:207–218 pmid:16618814

13. Angrisani A, Vicidomini R, Turano M, Furia M. Human dyskerin: beyond telomeres. Biol Chem 2014; 395: 593–610. pmid:24468621

14. Mason PJ, Bessler M. The genetics of dyskeratosis congenita. Cancer Genet. 2011; 204: 635–645 pmid:22285015

15. Zenz T, Mertens D, Kuppers R, Dohner H, Stilgenbauer S. From pathogenesis to treatment of chronic lymphocytic leukaemia. Nature Rev Cancer 2010; 10: 37–50.

16. Damle RN, Wasil T, Fais F, Ghiotto F, Valetto A, Allen SL, et al. Ig V gene mutation status and CD38 expression as novel prognostic indicators in chronic lymphocytic leukemia. Blood 1999; 94: 1840–1847. pmid:10477712

17. Hamblin TJ, Davis Z, Gardiner A, Oscier DG, Stevenson FK. Unmutated Ig V(H) genes are associated with a more aggressive form of chronic lymphocytic leukemia. Blood 1999; 94:1848–1854. pmid:10477713

18. Döhner H, Stilgenbauer S, Benner A, Leupolt E, Kröber A, Bullinger L, et al. Genomic aberrations and survival in chronic lymphocytic leukemia. N Eng J Med 2000; 343:1910–1916.

19. Travella A, Ripollés L, Aventin A, Rodríguez A, Bezares RF, Caballín MR, et al. Structural alterations in chronic lymphocytic leukaemia. Cytogenetic and FISH analysis. Hematol Oncol 2013; 31: 339–347.

20. Roos G, Kröber A, Grabowski P, Kienle D, Bühler A, Döhner H, et al. Short telomeres are associated with genetic complexity, high-risk genomic aberrations, and short survival in chronic lymphocytic leukemia. Blood 2008; 111: 2246–2252. pmid:18045969

21. Sellmann L, de Beer D, Bartels M, Opalka B, Nückel H, Dührsen U, et al. Telomeres and prognosis in patients with chronic lymphocytic leukemia. Int J Hematol 2011; 93: 74–82. pmid:21203871

22. Rampazzo E, Bonaldi L, Trentin L, Visco C, Keppel S, Giunco S, et al. Telomere length and telomerase levels delineate subgroups of B-cell chronic lymphocytic leukemia with different biological characteristics and clinical outcomes. Haematologica 2012; 97: 56–63. pmid:21933855

23. Dos Santos P, Panero J, Palau Nagore V, Stanganelli C, Bezares RF, et al. Telomere shortening associated with increased genomic complexity in chronic lymphocytic leukemia. Tumour Biol 2015; 36: 8317–8324. pmid:26008147

24. Mansouri L,Grabowski P, Degerman S, Svenson U, Gunnarsson R, Cahill N, et al. Short telomere length is associated with NOTCH1/SF3B1/TP53 aberrations and poor outcome in newly diagnosed chronic lymphocytic leukemia patients. Am J Hematol 2013; 88: 647–651. pmid:23620080

25. 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: 4918–4922. pmid:9810000

26. Hoxha M, Fabris S, Agnelli L, Bollati V, Cutrona G, Matis S, et al. Relevance of telomere/telomerase system impairment in early stage chronic lymphocytic leukemia. Genes Chrom Cancer 2014; 53: 612–621. pmid:24706380

27. Hallek M, Cheson BD, Catovsky D, Caligaris-Cappio F, Dighiero G, Döhner H, et al. Guidelines for the diagnosis and treatment of chronic lymphocytic leukemia: a report from the International Workshop on Chronic Lymphocytic Leukemia updating the National Cancer Institute-Working Group 1996 guidelines. Blood 2008; 11: 446–456.

28. Rai KR, Sawitsky A, Cronkite EP, Chanana AD, Levy RN, Pasternack BS. Clinical staging of chronic lymphocytic leukemia. Blood 1975; 46: 219–234. pmid:1139039

29. Chiodi I, Belgiovine C, Zongaro S, Ricotti R, Horard B, Lossani A, et al. Super-telomeres in transformed human fibroblasts. Biochim Biophys Acta. 2013; 1833: 1885–1893. pmid:23570868

30. Bièche I, Noguès C, Paradis V, Olivi M, Bedossa P, Lidereau R, et al Quantitation of hTERT gene expression in sporadic breast tumors with a real-time reverse transcription-polymerase chain reaction assay. Clin Cancer Res 2000; 6: 452–459. pmid:10690523

31. Ferri CA, Bianchini M, Bengió RM, Moiraghi EB, Gonzalez MS, Noriega MF, et al. Clinical activity of ponatinib in one patient with chronic myeloid leukemia in chronic phase with e19a2 transcript and T315I mutation. Eur J Haematol 2015; 94: 270–272. pmid:24766374

32. Panero J, O'Callaghan NJ, Fenech M, Slavutsky I. Absolute qPCR for measuring telomere length in bone marrow samples of plasma cell disorders. Mol Biotechnol 2015; 57: 155–159. pmid:25311116

33. O’Callaghan N., Dhillon V., Thomas P., and Fenech M. A quantitative real-time PCR method for absolute telomere length. Biotechniques 2008; 44: 807–809. pmid:18476834

34. Stanganelli C, Travella A, Bezares R, Slavutsky I. Immunoglobulin gene rearrangements and mutational status in argentinian patients with chronic lymphocytic leukemia. Clin Lymphoma Myeloma Leuk 2013; 13: 447–457. pmid:23665144

35. Díaz de la Guardia R, Catalina P, Panero J, Elosua C, Pulgarin A, López MB, et al. Expression profile of telomere-associated genes in multiple myeloma. J Cell Mol Med 2012; 16: 3009–3021. pmid:22947336

36. Poncet D, Belleville A, t'kint de Roodenbeke C, Roborel de Climens A, Ben Simon E, Merle-Beral H, et al. Changes in the expression of telomere maintenance genes suggest global telomere dysfunction in B-chronic lymphocytic leukemia. Blood 2008; 111: 2388–2391. pmid:18077792

37. Sbarrato T, Horvilleur E, Pöyry T, Hill K, Chaplin LC, Spriggs RV, et al. A ribosome-related signature in peripheral blood CLL B cells is linked to reduced survival following treatment. Cell Death Dis 2016; 7: e2249 pmid:27253413

38. Damle RN, Batliwalla FM, Ghiotto F, Valetto A, Albesiano E, Sison C, et al. Telomere length and telomerase activity delineate distinctive replicative features of the B-CLL subgroups defined by immunoglobulin V gene mutations. Blood 2004; 103: 375–382. pmid:14504108

39. Terrin L, Trentin L, Degan M, Corradini I, Bertorelle R, Carli P, et al. Telomerase expression in B-cell chronic lymphocytic leukemia predicts survival and delineates subgroups of patients with the same IgVH mutation status and different outcome. Leukemia 2007; 21: 965–972. pmid:17344921

40. Tchirkov A, Chaleteix C, Magnac C, Vasconcelos Y, Davi F, Michel A et al. hTERT expression and prognosis in B-chronic lymphocytic leukemia. Ann Oncol 2004; 15: 1476–1480. pmid:15367406

41. Panero J, Arbelbide J, Fantl DB, Rivello HG, Kohan D, Slavutsky I. Altered mRNA expression of telomere-associated genes in monoclonal gammopathy of undetermined significance and multiple myeloma. Mol Med 2010; 16: 471–478. pmid:20644899

42. Cao Y, Li H, Deb S, Liu JP. TERT regulates cell survival independent of telomerase enzymatic activity. Oncogene 2002; 21: 3130–3138. pmid:12082628

43. Folini M, Brambilla C, Villa R, Gandellini P, Vignati S, Paduano F et al. Antisense oligonucleotide-mediated inhibition of hTERT, but not hTERC, induces rapid cell growth decline and apoptosis in the absence of telomere shortening in human prostate cancer cells. Eur J Cancer 2005; 41: 624–634. pmid:15737568

44. von Stedingk K, Koster J, Piqueras M, Noguera R, Navarro S, Pahlman S, et al. snoRNPs regulate telomerase activity in neuroblastoma and are associated with poor prognosis. Transl. Oncol 2013; 6: 447–457. pmid:23908688

45. Fernandez-Garcia I, Marcos T, Muñoz-Barrutia A, Serrano D, Pio R, Montuenga LM, et al. Multiscale in situ analysis of the role of dyskerin in lung cancer cells. Integr Biol (Camb) 2013; 5: 402–413.

46. Panero J, Stella F, Schutz N, Fantl DB, Slavutsky I. Differential expression of non-shelterin genes associated with high telomerase levels and telomere shortening in plasma cell disorders. PLoS One. 2015; 14; 10: e0137972. pmid:26366868

47. Lin P, Mobasher ME, Hakakian Y, Kakarla V, Naseem AF, Ziai H, et al. Differential requirements for H/ACA ribonucleoprotein components in cell proliferation and response to DNA damage. Histochem Cell Biol 2015; 144: 543–558. pmid:26265134