Abstract
Background
Hypertrophy of the nucleolus is a distinctive cytological feature of malignant cells and corresponds to aggressive behaviour. This study aimed to identify the key gene associated with nucleolar prominence (NP) in breast cancer (BC) and determine its prognostic significance.
Methods
From The Cancer Genome Atlas (TCGA) cohort, digital whole slide images identified cancers having NP served as label and an information theory algorithm was applied to find which mRNA gene best explained NP. Dyskerin Pseudouridine Synthase 1 (DKC1) was identified. DKC1 expression was assessed using mRNA data of Molecular Taxonomy of Breast Cancer International Consortium (METABRIC, n = 1980) and TCGA (n = 855). DKC1 protein expression was assessed using immunohistochemistry in Nottingham BC cohort (n = 943).
Results
Nuclear and nucleolar expressions of DKC1 protein were significantly associated with higher tumour grade (p < 0.0001), high nucleolar score (p < 0.001) and poor Nottingham Prognostic Index (p < 0.0001). High DKC1 expression was associated with shorter BC-specific survival (BCSS). In multivariate analysis, DKC1 mRNA and protein expressions were independent risk factors for BCSS (p < 0.01).
Conclusion
DKC1 expression is strongly correlated with NP and its overexpression in BC is associated with unfavourable clinicopathological characteristics and poor outcome. This has been a detailed example in the correlation of phenotype with genotype.
Similar content being viewed by others
Background
Breast cancer (BC) is the most common cancer diagnosed in women worldwide, accounting for ~1 in 10 new cancer diagnoses each year and is the second most common cause of death1,2 due to cancer. BC is a heterogeneous disease with variable morphologies and response to therapy. Some morphological features, especially histological grade, have been well validated to have a strong prognostic value and their assessment helps in prognostic stratification of BC patients for treatment decisions.3
In the Nottingham cohort, nucleolar prominence (NP) has recently been shown to be a significant predictor for patient outcome as well as of being highly correlated with tumour grade. Since the NP is a distinctive morphological attribute, it is hypothesised to possibly serve as a substitute for the highly subjective pleomorphism component score of the Nottingham BC grading.4 Consequently, it is deemed imperative to explore the correlations between the nucleolar phenotype and genotype.
The major function of the nucleolus is synthesis and assembly of ribosomes5 where both are associated with malignant transformation and cancer progression.6,7 Indeed, ribosome biogenesis depends on the cancer growth rate, which is directly related to nucleolar size of malignant cells. Nucleolar size and cell kinetics’ parameters are interrelated because of the increasing rate of ribosome biogenesis in proliferating cells.8 In some solid cancer and haematological malignancies, the ribosome biogenesis rate increases as a consequence of overexpression of the oncogene c-Myc, which controls all the steps of ribosome biogenesis.9 Despite the biological and clinical significance of NP in BC, the key gene associated with NP and its prognostic significance remains to be defined.
Dyskerin Pseudouridine Synthase 1 (DKC1) is a predominantly nucleolar protein encoded by DKC1 gene and mapped at Xq28.10 DKC1 is a crucial component of the telomerase complex and is required for normal telomere maintenance and post-transcriptional processing of precursor rRNA. Therefore, DKC1 is necessary for tumour cell progression through mechanisms related to its function in the processing of rRNA precursor.11 Usually, clinically indolent and slow-growing tumours express lower levels of DKC1 and its inhibition slows or hinders the proliferation in most cell types.11,12 Through various deprivation of function approaches, emerging evidence suggests that DCK1 may regulate other cellular processes, including IRES-mediated translation, telomere maintenance independent of telomere length regulation, mitosis, transcription and possibly microRNA processing.13,14 Upregulation of DKC1 expression has been reported in several human cancers including hepatocellular carcinoma,15 neuroblastoma,16 lymphoma,17 melanoma,18 prostate cancer,19 colorectal cancer20 and ovarian carcinoma.21
Methods
Principle of DKC1 selection
We have applied an information theory (IT) approach to The Cancer Genome Atlas (TCGA) breast cancer dataset. The IT approach was used for feature selection to identify the key gene associated with NP, which was assessed morphologically in full face invasive BC sections stained with haematoxylin and eosin (H&E) using digital whole slide images (WSI) as explained in our previous study.4 The TCGA BC cohort was employed since it contains satisfactory whole slide images and mRNA-seq2 data present in 743 cancers. In the IT approach utilised, the nucleolar score served as a label and 20,339 mRNA transcripts served as predictor variables. The IT algorithm is a ‘greedy’ algorithm and reduces the number of features selected. “Greedy” is the term used in the machine learning community to describe an algorithm, which selects the optimal feature at each step and does not alter any choices already made based on findings from future choices.22 The analysis showed that the attribute, which exhibited the highest mutual information (information gain) with NP, was DKC1. Moreover, the detection of DKC1 was supported by LASSO regression feature selection. The required functions for LASSO were obtained from R library Glmnet.23 LASSO regression is also capable of reducing the number of predictors and thereby allowing for a focused study of a few attributes.24 Therefore, by applying the IT approach and LASSO regression, the selection was limited to a single gene (DKC1). This was followed by evaluating DKC1 mRNA and DKC1 protein expression in large clinically annotated cohorts of BC to evaluate its clinicopathological and prognostic value in invasive BC as described below.
Study cohorts for transcriptomic analysis
The discovery of DKC1 was by study of TCGA cohort. The TCGA was also used to assess the possible correlation between DKC1 mRNA expression and the variables recorded in this cohort.25 The Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) cohort (n = 1980) was used to evaluate DKC1 gene copy number (CN) aberrations and gene expression.26 Genomic and transcriptomic data for the METABRIC cohort had been obtained using the Affymetrix SNP 6.0 and Illumina HT-12v3 platforms, respectively.26 The association between DKC1 mRNA expression, copy number aberrations and clinicopathological parameters, molecular subtypes and patient outcome was investigated. Breast Cancer Gene Expression Miner online dataset v4.3 (http://bcgenex.centregauducheau.fr/BC-GEM/GEM-requete.php) was also used as external validation of DKC1 mRNA expression.
Study cohort for protein expression
The Nottingham BC patient cohort was used to evaluate the immunohistochemical (IHC) expression of DKC1. This cohort is a well-characterised large series (n = 943) of invasive BC patients aged ≤70 years and presented at Nottingham City Hospital between 1999 and 2006. The cohort has long-term clinical follow-up and clinicopathological data included patient’s age at diagnosis, histological tumour type, tumour grade, tumour size, lymph node status, Nottingham Prognostic Index (NPI) and lymphovascular invasion (LVI). Patient outcome data were obtained including BC-specific survival (BCSS), defined as the time (in months) from the date of primary surgical treatment to the time of death from BC and distant metastasis free survival (DMFS) defined as time (in months) from primary surgical treatment until the first event of distant metastasis. Patients were treated based on tumour features, NPI and hormone receptor status. Endocrine therapy was given to patients who had oestrogen receptor positivity (ER+) tumours with high NPI scores (>3.4), whereas no adjuvant therapy was given to patients with ‘good’ NPI scores (≤3.4). Premenopausal patients with moderate and poor NPI scores were candidates for chemotherapy, while postmenopausal patients with ‘moderate’ or ‘poor’ NPI scores were given hormonal therapy only. Classical treatment of cyclophosphamide, methotrexate and fluorouracil (CMF) was used as a therapy for patients presented with absence of ER expression and clinically fit to receive chemotherapy. None of the patients in the current study cohort received neoadjuvant therapy. Data related to the expression of ER, progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER2) as well as Ki67 were available.27,28,29,30 Molecular subtypes were based on tumour immunohistochemical (IHC) profile and the Elston–Ellis31 mitotic score as ER+/HER2−; low proliferation (mitotic score 1), ER+/HER2− high proliferation (mitotic scores 2 and 3), HER2-positive class: HER2+ regardless of ER status, TN: ER−, PR− and HER2−.32 There was no significant differences in the distribution of the clinicopathological parameters between the Nottingham and the METABRIC cohorts (all correlation coefficients ≥0.948, all p < 0.0001)33 (Supplementary Table 1).
DKC1 validation by western blotting
The antibody specificity of anti-DKC1 antibody (EPR10399, Abcam, UK) was validated using western blotting (WB) performed on cell lysates of a wild and transfected MDA-MB-231 human breast cancer cell line (American Type Culture Collection; Rockville, MD, USA). The forward transfection of siRNA procedure was followed according to DKC1 siRNAs manufacturer’s instructions. In brief, the cells were seeded in 6-well plate at a cell density of 3 × 105 cells per well and incubated overnight in 37 °C 5% CO2 incubator. The following day, the cells reached about 40% confluence and were transfected with 10 nM of three different IDs of DKC1 siRNA (Cat#:4392420, ThermoFisher Scientific, UK). A transfection with 10 nM scrambled siRNA sequence (Cat#:4390843, ThermoFisher Scientific, UK) was carried out in the experiment and considered as a negative control. DKC1 protein expression of untransfected & transfected cells was then determined by the Western blotting. Briefly, after collecting the cell lysates, a dilution of 1:1000 of the primary antibody and 1:15000 IRDye 800CW Donkey anti-rabbit secondary antibody (LI-COR Biosciences) were applied, and 5% milk /PBS-Tween (0.1%) (Marvel Original Dried Skimmed Milk, Premier Food Groups Ltd., UK) was used for blocking and antibodies incubation. Mouse monoclonal anti-β-actin primary antibody (1:5000) (Sigma–Aldrich, UK) with IRDye 800CW Donkey anti-mouse fluorescent secondary antibody (LI-COR Biosciences) were used to visualise a marker of endogenous control. Visualisation of DKC1 band was done by using the Odyssey Fc with Image Studio 4.0 (LI-COR Biosciences).
Tissue microarrays and immunohistochemical analysis
Invasive BC tissues were previously arrayed as tissue microarrays (TMA) using the Grand Master® (3D HISTECH®, Budapest, Hungary).34 IHC staining was performed on 4 μm TMA thick sections using the Novocastra Novolink™ Polymer Detection Systems kit (Code: RE7280-K, Leica, Biosystems, Newcastle, UK). Antigen retrieval was performed in citrate buffer pH 6.0 using a microwave (Whirlpool JT359 Jet Chef 1000 W) for 20 min. Rabbit monoclonal DKC1 was diluted at 1:50 in Leica antibody diluent (RE AR9352, Leica, Biosystems, UK) and incubated with the sections for 60 min at room temperature. A negative control was obtained by omitting the incubation with primary antibody while formalin fixed placenta tissue was used as a positive control according to manufacturer’s datasheet.
Assessment of DKC1 protein expression
Scanning of TMA stained sections into high-resolution digital images was performed by using a NanoZoomer scanner (NanoZoomer; Hamamatsu Photonics, Welwyn Garden City, UK) at ×20 magnification. Scoring of DKC1 nuclear and nucleolar4 expression was evaluated based on a semi-quantitative scoring using modified histochemical score (H-score), where the intensity of staining was multiplied by the percentage of positive cells in the tissue for each intensity, producing a score ranging between 0 and 300.35 A score index of 0, 1, 2 and 3 corresponding to negative, weak, moderate and strong respectively were used for intensity. The percentage (0–100) of positive cells for each intensity was evaluated subjectively. All non-representative cores including folded tissue during processing and staining, cores with only normal breast tissue and cores with invasive tumour <15% of core surface area were excluded from scoring. All the cores were scored by a trained observer (K. Elsharawy) blinded of histopathological and patient outcome data. Further, to test the interobserver’s reproducibility of the scoring, a subset of TMA cores (10%) was randomly selected and double scored by a second trained observer (M. Aleskandarany). Moreover, for further evaluation of scoring reproducibility, 20% of the cases were double scored by the main observer (K. Elsharawy) after 5 months washout period blind from the first scores.
Statistical analysis
IBM-SPSS statistical software 24.0 (SPSS, Chicago, IL, USA) was used in statistical analysis. Interobserver agreement in DKC1 IHC scoring was assessed using intraclass correlation coefficient (ICC). Dichotomisation of DKC1 proteomic and transcriptomic levels expression was determined based on the prediction of BCSS using X-tile bioinformatics software version 3.6.1 (School of Medicine, Yale University, New Haven, USA).36 The H-scores of 110 and 10 were the optimal cut-off values of DKC1 nuclear and nucleolar protein expression. Continuous data of DKC1 mRNA and DKC1 protein expression were used to assess the association with clinicopathological parameters. Differences between three or more groups were investigated using one-way analysis of variance (ANOVA) with the post-hoc Tukey multiple comparison test (for parametric data) or Kruskal–Wallis test (for non-parametric distribution). Student t-test (parametric data) or Mann–Whitney test (non-parametric distribution) were used to evaluate the differences between two groups. Spearman’s correlation coefficient was calculated to examine the association between continuous variables. Univariate analysis was visualised using Kaplan–Meier curves and significance was assessed by log-rank test. Cox’s proportional hazard regression models were built for the multivariate survival analysis to adjust for confounding factors. P values were adjusted by using Bonferroni correction for multiple testing. For all tests, p value < 0.05 was considered as statistically significant. This study followed the reporting recommendations for tumour markers prognostic studies (REMARK) criteria.37
Results
In this study, we have applied the bespoke bioinformatics tools to identify the key genes associated with NP and this identified DKC1 as the target gene. Then, DKC1 was investigated at the transcriptomic, genetic and protein levels.
DKC1 mRNA expression and CN aberrations
High DKC1 mRNA expression (log2 intensity >9.4) was observed in 709/1970 (36%) of the METABRIC cases. In all, 77/1980 (4%) of cases showed DKC1 CN gain, whereas 115/1980 (6%) showed a CN loss. A significant association was observed between DKC1 CN variation and DKC1 mRNA expression (p < 0.0001) (Fig. 1a).
DKC1 protein expression in breast cancer
Prior to IHC staining, the specificity of the antibody used was validated using WB performed on a DKC1 siRNA transfected BC cell line. A specific band at the predicted DKC1 molecular weight (58 kDa) was detected for proteins extracted from untransfected cells and those transfected with scrambled siRNA sequences. In addition, the DKC1 band intensity was significantly reduced with proteins extracted from DKC1 siRNAs transfected cells, confirming the specificity of the antibody utilised. A single band was observed at β-actin molecular weight (42 kDa) demonstrating the uniformity of loaded protein quantities (Fig. 2a).
DKC1 protein expression was observed in the nucleus and nucleoli of invasive BC cells, with expression levels varying from absent to strong (Fig. 2b–d). Strong concordance was observed between the two observers in DKC1 immuno-scoring in 10% of the cases (ICC = 0.864, p < 0.0001 for nuclear expression and ICC = 0.781, p < 0.001 for nucleolar expression). Moreover, second scoring of 20% of cases after 5 months washout period confirmed concordance ((ICC = 0.822, p < 0.0001 for nuclear expression and ICC = 0.804, p < 0.0001 for nucleolar expression). At the optimal DKC1 cut-off values (H-score 110 and 10, respectively), High DKC1 nuclear and nucleolar expression were observed in 574/942 (61%) and 153/942 (16%) of the informative tumours, respectively. There was a significant positive correlation between DKC1 nuclear and nucleolar expression (n = 429) (correlation coefficient = 0.143, p < 0.0001).
Correlation of DKC1 mRNA and protein expression with clinicopathological parameters
High DKC1 mRNA expression was significantly associated with younger patient age, larger tumour size, higher tumour grade and poorer NPI (p < 0.001, p = 0.024, p < 0.0001 and p < 0.0001) as shown in Fig. 1b-e, respectively. These associations were confirmed using the Breast Cancer Gene-Expression Miner v4.3 (Supplementary Fig. 1A-C)
In the TCGA BC dataset, similar associations, as described above, were observed with clinicopathological parameters. In particular, high DKC1 mRNA expression was significantly associated with high nucleolar score 34 (p < 0.0001, Supplementary Table 2).
High expression of DKC1 protein whether in the nucleus and/or nucleoli was associated with aggressive features of BC including higher tumour grade (p < 0.0001), larger tumour size (p = 0.04 only with nucleolar expression), higher mitotic scores (p < 0.0001), increased nuclear pleomorphism (p < 0.0001), higher scores of nucleolar prominence (p < 0.001), poor NPI (p < 0.0001) and the invasive ductal no special histological type (NST) (p < 0.0001), Table 1.
DKC1 expression and other markers
The correlation of DKC1 mRNA with other relevant genes was investigated using the METABRIC and TCGA datasets. The genes were chosen based on published information, being either regulatory genes or those that share or support DKC1 biological function especially those primarily involved in the ribosomal biogenesis. DKC1 was positively associated with GAR1 (p < 0.0001), NOP10 (p < 0.001) and NHP2 (p < 0.0001). Moreover, there was a significant association between DKC1 MKI67 and the regulatory genes c-Myc (all p < 0.001) (Supplementary Table 3). High DKC1 mRNA expression was associated with those tumours, which showed TP53 mutations (p < 0.0001, Supplementary Table 4). Also, the statistical analysis showed a significant positive association of high DKC1 protein expression with high Ki67 (χ2 = 8.815, p = 0.003).
DKC1 mRNA and protein expression in BC molecular subtypes
At the transcriptomic level in METABRIC cohort, high DKC1 expression was significantly associated with hormone receptor negative (ER− and PR−), HER2+ tumours and TNBC (all p < 0.0001) as shown in Supplementary Table 4. Similar results were observed upon analysing the publicly available gene-expression data available on the Breast Cancer Gene-Expression Miner v4.3 online platform (Supplementary Fig. 1D–G) and TCGA datasets (Supplementary Table 4).
Regarding the association with the intrinsic PAM50 subtypes,38 high expression of DKC1 mRNA was observed in basal-like, HER2+ and Luminal B tumours (Fig. 1f, p < 0.0001). These findings were confirmed using the Breast Cancer Gene-Expression Miner v4.3 (Supplementary Fig. 1H). In the SCMGENE subtypes, high expression of DKC1 mRNA was observed in the ER−/HER2− cases followed by ER+/HER2− high proliferation class (p < 0.0001, Fig. 1g).
DKC1 nuclear and nucleolar protein expression was associated with negative ER status (p = 0.04 and p < 0.0001 respectively). Moreover, DKC1 nucleolar protein showed a significant correlation within HER2+ and triple negative (TN) tumours (both p < 0.0001), Table 2.
There was a higher protein expression of DKC1 (nuclear & nucleolar) in the ER + high proliferative tumours than in the other molecular subtypes (p < 0.0001) as shown in Table 1.
Correlation of DKC1 mRNA and protein expression with patient outcome
In METABRIC cohort, high DKC1 mRNA expression was associated with poor BCSS in all cases (HR = 1.5, 95% CI = 1.3–1.8; p < 0.0001). Moreover, DKC1 mRNA expression was predictive of BCSS only in luminal B cases (HR = 1.5, 95%CI = 1.1–2.1; p = 0.015) as shown in Supplementary Fig. 2A–E. The relationship between high DKC1 mRNA expression and poor patient outcome in ER+ disease, but not ER− disease, was shown using the TCGA cohort (Supplementary Fig. 3).
Both high DKC1 nuclear and nucleolar protein expressions, when assessed individually, were associated with poor outcome (HR = 2.5, 95%CI = 1.7–3.7; p < 0.0001 and HR = 1.5, 95%CI = 1.1–2.2; p = 0.038, respectively) Fig. 3a, b. When the analysis was limited to molecular subtypes, high expression of DKC1 nuclear protein was significantly associated with poor outcome in ER+ high proliferation tumours (HR = 4.4, 95% CI = 1.6–12.3; p = 0.002), HER2+ tumours (HR = 2.6, 95% CI = 1.1–6.7; p = 0.039) and TNBC (HR = 1.5, 95% CI = 1.1–6.2; p = 0.035) Fig. 3d–g. However, no significant association of DKC1 nucleolar protein expression was identified with outcome in BC subtypes (p > 0.05).
For further analysis, combinatorial DKC1 protein expression groups were created [i.e. low nuclear/low nucleolar, high nuclear/low nucleolar expression, low nuclear/high nucleolar and high nuclear/high nucleolar]. A significant difference in patient survival was observed between these four groups, where tumours with low nuclear and low nucleolar DKC1 expression showed the best outcome, whereas the tumours with high nuclear and high nucleolar expression showed the worst outcome (p < 0.0001) Fig. 3c.
The multivariate Cox-proportional models, including other prognostic covariates such as tumour size, grade and nodal stage, showed that DKC1 nuclear, combinatorial protein expression and DKC1 mRNA in the METABRIC dataset were independent predictors for poor prognosis in whole cases (p = 0.001, HR 2.037, 95% CI = 1.373–3.023, p = 0.003, HR 2.746, 95% CI = 1.484–5.083 and p = 0.006, HR 1.316, 95% CI = 1.097–1.579) as shown in Table 3 and Supplementary Table 5.
In addition, high DKC1 nuclear protein expression was significantly associated with shorter distant metastases-free survival (DMFS) (HR = 2.1, 95% CI = 1.5–2.9; p < 0.0001). Likewise, combinatorial protein expression was associated with shorter DMFS (HR = 1.3, 95% CI = 1.1–1.4; p = 0.001) (Supplementary Fig. 4).
Discussion
In malignant tumours, the number and size of nucleoli are usually an indication of the rate of ribosome production, which is regarded as a major metabolic requisite for cell growth and proliferation.39,40 Nucleolar function and size are directly related to cell doubling time in cancer cells and quantitative morphometric evaluation of nucleolar size was considered as a cytological parameter of the tumour cells proliferation rate.4,41 In this study, we assessed NP in the TCGA breast cancer dataset as previously described,4 and used two greedy algorithms information theory and validated it using LASSO regression test22,23,24 to identify genes driving NP. This demonstrated that out of the 20,339 genes investigated, DKC1 was the top differentially expressed gene. Then DKC1 expression was evaluated at the proteomic, transcriptomic and genomic levels in large cohorts of invasive BC.
There were significant associations between high DKC1 mRNA in TCGA breast cancer dataset and DKC1 protein expression in Nottingham cohort with high nucleolar scoring.4 These findings supported our hypothesis that DKC1 plays a role in the nucleoli appearance and size likely through its mechanism in ribosomal biogenesis.
Our results also showed positive correlations between nuclear and nucleolar DKC1 protein expression in the breast tumour cells. It was reported that newly synthesised DKC1 initially localises to the nucleoplasm, followed by consecutive translocation to the nucleoli and the nuclear Cajal bodies. Usually, colocalisation of DKC1 on the nuclear Cajal bodies occurred only when it had already accumulated in the nucleoli.42 DKC1 is involved in the pseudouridination and processing of small spliceosomal RNAs through its binding to H/ACA small Cajal body RNAs.43
The current study confirms the significant association between the high expression of DKC1, at both protein and mRNA levels, and clinicopathological parameters characteristics of poor prognosis and with shorter survival; findings which are in-line with other studies.15,44,45 Some studies have also confirmed that DKC1 overexpression is involved in tumorigenic processes and has prognostic value in numerous types of cancer.19,21,46 The association between DKC1 mRNA and shorter survival was identified in both METABRIC and TCGA cohorts. Moreover, our analysis revealed that the prognostic significance of DKC1 protein and mRNA in BC was independent of other variables, demonstrating its potential clinical relevance in improving survival rate prediction.
When BC molecular subtypes were considered, the significant association between DKC1 protein and poor patient outcome was observed in the ER+ high proliferation (i.e. luminal B), HER2+ and TNBC classes whereas the high mRNA expression was only limited to the luminal B subtype. The most common type of BC constituting nearly 55–70% is the ER+/luminal tumour, and those tumours are variable in terms of recurrence, mortality rates and disease prognosis.47,48 These observations further endorse DKC1 functions in playing crucial roles in tumour growth and progression
DKC1 performs two fundamental functions for cell proliferation. First, DKC1 is a component of the H/ACA small nucleolar ribonucleoprotein particles (snoRNPs) involved the pseudouridylation of ribosomal RNA (rRNA) molecules and necessary for their processing. Second, it is required for telomerase activity by stabilising the telomerase RNA component.11 The faster the rate of cell proliferation, the higher the demand for protein production, which is compatible with increased rRNA synthesis.49,50 It has been reported that any dysregulation of DKC1 levels results in defects of ribosome biogenesis and a reduction of rRNA pseudouridylation, which in turn hinders the normal ribosome rRNA processing rate.51 For instance, Montanaro et al. have demonstrated that reduced DKC1 gene expression by specific RNA interference in BC cell lines resulted in a reduction of rRNA pseudouridylation, which subsequently effected the survival of proliferating cells.52 The role of DKC1 in mitosis was also confirmed,11 where dyskerin was identified as one of seventy genes which correlated with the development of aneuploidy. Alawi et al. have demonstrated that dyskerin expression peaks during G2/M and loss of dyskerin function has a widely disruptive effect on mitosis and triggers the spindle-assembly checkpoint.13 Our findings showed that high DKC1 expression was significantly associated with proliferation as assessed by Ki67 labelling index, which was also observed in other studies in BC52 hepatocellular carcinoma15 and prostate cancer,19 confirming that DKC1 is critical for mitotic progression and proliferation in these cancers.
DKC1 is the direct and conserved transcriptional target of c-Myc,53 which explains the strong correlation between its upregulation and active cell proliferation.54 In our study, we observed a significant positive association between DKC1 and c-Myc in mRNA expression. Previous studies have demonstrated that tumour oncogene c-Myc controls the transcription of DKC1 gene in addition to other proteins, which are required for rRNA processing.9,55 TP53 mutations were also highly prevalent in breast tumours with high DKC1 mRNA expression in METABRIC. There is mounting evidence that the usual increase of ribosome biogenesis (one of the main functions of DKC1) in cancer cells is the consequence of frequent alterations of two major tumour suppressors, TP53 and retinoblastoma (RB) genes.55 In addition, tumours with altered p53 and/or retinoblastoma protein pRb functions are characterised by significantly larger/more conspicuous nucleoli than tumours with normally functioning p53 and pRb.56
We further investigated the association of DKC1 expression with other H/ACA ribonucleoproteins, NHP2, NOP10 and GAR1, which play important roles in disease progression. DKC1, NHP2 and NOP10 form a core trimer that directly binds to H/ACA RNAs. The three proteins are interdependent with each other for stability and also regulate constancy of the bound RNAs.57 GAR1 binds only to DKC1 and is needed for a proper functioning of the H/ACA RNPs, but its absence does not reduce the stability of the rRNA.58 These findings confirmed the significant positive correlation between H/ACA ribonucleoproteins and DKC1 in BC.43 Alterations in DKC1 expression will potentially disrupt the biogenesis of H/ACA pathway and consequently affect ribosome synthesis and impair cell proliferation.
In the last decade, there have been a few attempts to construct DKC1 inhibitors. However, one in silico study successfully determined a small molecule inhibitor (Pyrazofurin) that exerted an ability to weaken the pharmacological and physiological activities of DKC1 through inhibiting its function in pseudouridylation of rRNA. Although Pyrazofurin failed to progress Phase 2 clinical trials; however, its chemical structure should continue to be exploited as a pharmacokinetic model to develop a potent, effective and safe DKC1 inhibitor that may eventually be used for BC highly expressing DKC1.59
A few limitations of this study findings are worth mentioning. On one hand, the semi-quantitative H-score method used to evaluate the immunohistochemical protein expression in the Nottingham cohort, might be regarded as having substantial subjectivity. This was addressed by double scoring a subset of the cancers to ensure the reproducibility and liability of the procedure. On the other hand, only one representative TMA core from each tumour tissue was arrayed and scored instead of considering replicates to express the tumour heterogeneity. This was due to the limited tissue resources in our biobank. However, to overcome the issue, 20 full face sections of randomly selected breast cancer cases were stained, prior to TMA application, with DKC1 antibody to assess the staining homogeneity and to evaluate the pertinence of using tissue microarrays (TMAs). These showed homogeneously distributed DKC1 expression, deeming the use of TMA to assess DKC1 an appropriate tissue platform, to mitigate the limited resources as well as testing the hypothesis a large BC cohort. Finally, the ‘weak’ but significant correlation between nuclear and nucleolar expression of DKC1 might also be regarded as a weakness point in the study. This might be due to the subjectivity of the method, which has been used in determining NP in our previous study.4
Conclusion
This study reveals a significant correlation between the morphological features of NP and an underlying molecular and protein description (DKC1). The importance of DKC1 was demonstrated in three independent datasets where each dataset contributed to the description of DKC1 from different perspectives. DKC1 is significantly associated with high nucleolar score and with poor prognostic characteristics and poor patients’ outcome. Overexpression of DKC1 appears to play a role in the proliferation and progression of the aggressive BC subtypes including the luminal B, TNBC and HER2 molecular subtypes. Findings here encourage further investigation of DKC1 as it might relate to guiding targeted therapies and to evaluate its role in response to chemotherapy.
References
Parada, H. Jr, Sun, X., Tse, C. K., Olshan, A. F. & Troester, M. A. Lifestyle patterns and survival following breast cancer in the Carolina Breast Cancer Study. Epidemiol. (Camb., Mass) 30, 83–92 (2019).
Andisha, N. M., McMillan, D. C., Gujam, F. J. A., Roseweir, A. & Edwards, J. The relationship between phosphorylation status of focal adhesion kinases, molecular subtypes, tumour microenvironment and survival in patients with primary operable ductal breast cancer. Cell. Signal. https://doi.org/10.1016/j.cellsig.2019.04.006 (2019).
Rakha, E. A., El-Sayed, M. E., Lee, A. H., Elston, C. W., Grainge, M. J., Hodi, Z. et al. Prognostic significance of Nottingham histologic grade in invasive breast carcinoma. J. Clin. Oncol. 26, 3153–3158 (2008).
Elsharawy, K. A., Toss, M. S., Abuelmaaty, S. R., Ball, G., Green, A. R., Aleskandarany, M. A. et al. Prognostic significance of nucleolar assessment in invasive breast cancer. Histopathology https://doi.org/10.1111/his.14036 (2019).
Lam, Y. W., Trinkle-Mulcahy, L. & Lamond, A. I. The nucleolus. J. Cell Sci. 118, 1335 (2005).
Montanaro, L., Trere, D. & Derenzini, M. Nucleolus, ribosomes, and cancer. Am. J. Pathol. 173, 301–310 (2008). e-pub ahead of print 2008/06/28.
Bywater, M. J., Poortinga, G., Sanij, E., Hein, N., Peck, A., Cullinane, C. et al. Inhibition of RNA polymerase I as a therapeutic strategy to promote cancer-specific activation of p53. Cancer cell 22, 51–65 (2012).
Derenzini, M., Trerè, D., Pession, A., Govoni, M., Sirri, V. & Chieco, P. Nucleolar size indicates the rapidity of cell proliferation in cancer tissues. J. Pathol. 191, 181–186 (2000).
van Riggelen, J., Yetil, A. & Felsher, D. W. MYC as a regulator of ribosome biogenesis and protein synthesis. Nat. Rev. Cancer 10, 301–309 (2010).
Heiss, N. S., Knight, S. W., Vulliamy, T. J., Klauck, S. M., Wiemann, S., Mason, P. J. et al. X-linked dyskeratosis congenita is caused by mutations in a highly conserved gene with putative nucleolar functions. Nat. Genet. 19, 32–38 (1998).
Alawi, F. & Lin, P. Dyskerin is required for tumor cell growth through mechanisms that are independent of its role in telomerase and only partially related to its function in precursor rRNA processing. Mol. Carcinog. 50, 334–345 (2011).
Lin, P., Mobasher, M. E. & Alawi, F. Acute dyskerin depletion triggers cellular senescence and renders osteosarcoma cells resistant to genotoxic stress-induced apoptosis. Biochem. Biophys. Res. Commun. 446, 1268–1275 (2014).
Alawi, F. & Lin, P. Dyskerin localizes to the mitotic apparatus and is required for orderly mitosis in human cells. PLoS ONE 8, e80805 (2013).
Fong, Y. W., Ho, J. J., Inouye, C. & Tjian, R. The dyskerin ribonucleoprotein complex as an OCT4/SOX2 coactivator in embryonic stem cells. eLife 2014 https://doi.org/10.7554/eLife.03573 (2014).
Liu, B., Zhang, J., Huang, C. & Liu, H. Dyskerin overexpression in human hepatocellular carcinoma is associated with advanced clinical stage and poor patient prognosis. PLoS ONE 7, e43147 (2012).
Westermann, F., Henrich, K.-O., Wei, J. S., Lutz, W., Fischer, M., König, R. et al. High Skp2 expression characterizes high-risk neuroblastomas independent of MYCN status. Clin. Cancer Res. 13, 4695–4703 (2007).
Piva, R., Pellegrino, E., Mattioli, M., Agnelli, L., Lombardi, L., Boccalatte, F. et al. Functional validation of the anaplastic lymphoma kinase signature identifies CEBPB and Bcl2A1 as critical target genes. J. Clin. Invest. 116, 3171–3182 (2006).
McDonald, S. L., Edington, H. D., Kirkwood, J. M. & Becker, D. Expression analysis of genes identified by molecular profiling of VGP melanomas and MGP melanoma-positive lymph nodes. Cancer Biol. Ther. 3, 110–120 (2004).
Sieron, P., Hader, C., Hatina, J., Engers, R., Wlazlinski, A., Müller, M. et al. DKC1 overexpression associated with prostate cancer progression. Br. J. Cancer (Mol. Diagnostics) 101, 1410 (2009).
Witkowska, A., Gumprecht, J., Glogowska-Ligus, J., Wystrychowski, G., Owczarek, A., Stachowicz, M. et al. Expression profile of significant immortalization genes in colon cancer. Int. J. Mol. Med. 25, 321–329 (2010).
Schaner, M. E., Ross, D. T., Ciaravino, G., Sorlie, T., Troyanskaya, O., Diehn, M. et al. Gene expression patterns in ovarian carcinomas. Mol. Biol. Cell 14, 4376–4386 (2003).
Brown, G., Pocock, A., Zhao, M.-J. & Luján, M. Conditional likelihood maximisation: a unifying framework for information theoretic feature selection. J. Mach. Learn Res. 13, 27–66 (2012).
Friedman, J., Hastie, T. & Tibshirani, R. Regularization paths for generalized linear models via coordinate descent. J. Stat. Softw. 33, 1–22 (2010).
Tibshirani, R., Bien, J., Friedman, J., Hastie, T., Simon, N., Taylor, J. et al. Strong rules for discarding predictors in lasso-type problems. J. R. Stat. Soc. 74, 245–266 (2012).
Ciriello, G., Gatza, M. L., Beck, A. H., Wilkerson, M. D., Rhie, S. K., Pastore, A. et al. Comprehensive molecular portraits of invasive lobular breast cancer. Cell 163, 506–519 (2015).
Curtis, C., Shah, S. P., Chin, S. F., Turashvili, G., Rueda, O. M., Dunning, M. J. et al. The genomic and transcriptomic architecture of 2000 breast tumours reveals novel subgroups. Nature 486, 346–352 (2012).
Aleskandarany, M. A., Abduljabbar, R., Ashankyty, I., Elmouna, A., Jerjees, D., Ali, S. et al. Prognostic significance of androgen receptor expression in invasive breast cancer: transcriptomic and protein expression analysis. Breast Cancer Res. Treat. 159, 215–227 (2016).
Rakha, E. A., Agarwal, D., Green, A. R., Ashankyty, I., Ellis, I. O., Ball, G. et al. Prognostic stratification of oestrogen receptor-positive HER2-negative lymph node-negative class of breast cancer. Histopathology 70, 622–631 (2017).
Rakha, E. A., Elsheikh, S. E., Aleskandarany, M. A., Habashi, H. O., Green, A. R., Powe, D. G. et al. Triple-negative breast cancer: distinguishing between basal and nonbasal subtypes. Clin. Cancer Res. 15, 2302–2310 (2009).
Muftah, A. A., Aleskandarany, M. A., Al-Kaabi, M. M., Sonbul, S. N., Diez-Rodriguez, M., Nolan, C. C. et al. Ki67 expression in invasive breast cancer: the use of tissue microarrays compared with whole tissue sections. Breast Cancer Res. Treat. 164, 341–348 (2017).
Elston, C. W. & Ellis, I. O. Pathological prognostic factors in breast cancer. I. The value of histological grade in breast cancer: experience from a large study with long-term follow-up. Histopathology 19, 403–410 (1991).
Senkus, E., Kyriakides, S., Ohno, S., Penault-Llorca, F., Poortmans, P., Rutgers, E. et al. Primary breast cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann. Oncol. 26, v8–v30 (2015).
El Ansari, R., Craze, M. L., Miligy, I., Diez-Rodriguez, M., Nolan, C. C., Ellis, I. O. et al. The amino acid transporter SLC7A5 confers a poor prognosis in the highly proliferative breast cancer subtypes and is a key therapeutic target in luminal B tumours. Breast Cancer Res. 20, 21 (2018).
Abd El-Rehim, D. M., Ball, G., Pinder, S. E., Rakha, E., Paish, C., Robertson, J. F. et al. High-throughput protein expression analysis using tissue microarray technology of a large well-characterised series identifies biologically distinct classes of breast cancer confirming recent cDNA expression analyses. Int. J. Cancer 116, 340–350 (2005).
McCarty, K. S. Jr & KS, Mc. Carty Sr Histochemical approaches to steroid receptor analyses. Semin. Diagn. Pathol. 1, 297–308 (1984).
Camp, R. L., Dolled-Filhart, M. & Rimm, D. L. X-tile. A new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin. Cancer Res. 10, 7252–7259 (2004).
McShane, L. M., Altman, D. G., Sauerbrei, W., Taube, S. E., Gion, M., Clark, G. M. et al. REporting recommendations for tumour MARKer prognostic studies (REMARK). Br. J. Cancer 93, 387–391 (2005).
Parker, J. S., Mullins, M., Cheang, M. C., Leung, S., Voduc, D., Vickery, T. et al. Supervised risk predictor of breast cancer based on intrinsic subtypes. J. Clin. Oncol. 27, 1160–1167 (2009).
Fischer, A. H., Bardarov, S. Jr & Jiang, Z. Molecular aspects of diagnostic nucleolar and nuclear envelope changes in prostate cancer. J. Cell. Biochem. 91, 170–184 (2004).
Hernandez-Verdun, D. The nucleolus: a model for the organization of nuclear functions. Histochem. Cell Biol. 126, 135–148 (2006).
Derenzini, M., Trere, D., Pession, A., Montanaro, L., Sirri, V. & Ochs, R. L. Nucleolar function and size in cancer cells. Am. J. Pathol. 152, 1291–1297 (1998).
Heiss, N. S., Girod, A., Salowsky, R., Wiemann, S., Pepperkok, R. & Poustka, A. Dyskerin localizes to the nucleolus and its mislocalization is unlikely to play a role in the pathogenesis of Dyskeratosis congenita. Hum. Mol. Genet. 8, 2515–2524 (1999).
Ye, K. H/ACA guide RNAs, proteins and complexes. Curr. Opin. Struct. Biol. 17, 287–292 (2007).
Zhang, M., Pan, Y., Jiang, R., Hou, P., Shan, H., Chen, F. et al. DKC1 serves as a potential prognostic biomarker for human clear cell renal cell carcinoma and promotes its proliferation, migration and invasion via the NFkappaB pathway. Oncol. Rep. 40, 968–978 (2018).
von Stedingk, K., Koster, J., Piqueras, M., Noguera, R., Navarro, S., Påhlman, S. et al. snoRNPs regulate telomerase activity in neuroblastoma and are associated with poor prognosis. Transl. Oncol. 6, 447–457 (2013).
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 111, 2388–2391 (2008).
Dawson, S. J., Rueda, O. M., Aparicio, S. & Caldas, C. A new genome-driven integrated classification of breast cancer and its implications. EMBO J. 32, 617–628 (2013).
Rakha, E. A., El-Sayed, M. E., Green, A. R., Paish, E. C., Powe, D. G., Gee, J. et al. Biologic and clinical characteristics of breast cancer with single hormone receptor positive phenotype. J. Clin. Oncol. 25, 4772–4778 (2007).
Mochizuki, Y., He, J., Kulkarni, S., Bessler, M. & Mason, P. J. Mouse dyskerin mutations affect accumulation of telomerase RNA and small nucleolar RNA, telomerase activity, and ribosomal RNA processing. Proc. Natl Acad. Sci. USA 101, 10756–10761 (2004).
Filipowicz, W. & Pogacic, V. Biogenesis of small nucleolar ribonucleoproteins. Curr. Opin. Cell Biol. 14, 319–327 (2002).
Ruggero, D., Grisendi, S., Piazza, F., Rego, E., Mari, F., Rao, P. H. et al. Dyskeratosis congenita and cancer in mice deficient in ribosomal RNA modification. Science 299, 259–262 (2003).
Montanaro, L., Brigotti, M., Clohessy, J., Barbieri, S., Ceccarelli, C., Santini, D. et al. Dyskerin expression influences the level of ribosomal RNA pseudo-uridylation and telomerase RNA component in human breast cancer. J. Pathol. 210, 10–18 (2006).
Alawi, F. & Lee, M. N. DKC1 is a direct and conserved transcriptional target of c-MYC. Biochem. Biophys. Res. Commun. 362, 893–898 (2007).
Alawi, F., Lin, P., Ziober, B. & Patel, R. Correlation of dyskerin expression with active proliferation independent of telomerase. Head Neck 33, 1041–1051 (2011).
Penzo, M., Montanaro, L., Treré, D. & Derenzini, M. The ribosome biogenesis-cancer connection. Cells 8, 55 (2019).
Trere, D., Ceccarelli, C., Montanaro, L., Tosti, E. & Derenzini, M. Nucleolar size and activity are related to pRb and p53 status in human breast cancer. J. Histochem. Cytochem. 52, 1601–1607 (2004).
Grozdanov, P. N., Roy, S., Kittur, N. & Meier, U. T. SHQ1 is required prior to NAF1 for assembly of H/ACA small nucleolar and telomerase RNPs. RNA 15, 1188–1197 (2009).
Darzacq, X., Kittur, N., Roy, S., Shav-Tal, Y., Singer, R. H. & Meier, U. T. Stepwise RNP assembly at the site of H/ACA RNA transcription in human cells. J. Cell Biol. 173, 207–218 (2006).
Rocchi, L., Barbosa, A. J., Onofrillo, C., Del Rio, A. & Montanaro, L. Inhibition of human dyskerin as a new approach to target ribosome biogenesis. PLoS ONE 9, e101971 (2014).
Acknowledgements
We thank Breast Cancer Now and PathLAKE project for supporting this study. These new Centres are supported by the Data to Early Diagnosis and Precision Medicine strand of the government’s Industrial Strategy Challenge, managed and delivered by UK Research and Innovation (UKRI).
Author information
Authors and Affiliations
Contributions
K.E.: scored all the cases, took the lead in writing the manuscript, data analysis and interpretation. O.M. and M.A.: helped in double scoring, data interpretation and reviewing the article. A.H., H.G. and M.I.A.: contributed in data analysis, interpretation, writing and reviewing the article. A.G. and L.D.: contributed in data analysis, study design and reviewing the article. E.R.: conceived and planned the presented idea, data interpretation and reviewing the article. All authors reviewed and approved the final version of the manuscript. All authors declared their contribution(s) to the study
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This work obtained ethics approval by the North West—Greater Manchester Central Research Ethics Committee under the title; Nottingham Health Science Biobank (NHSB), reference number 15/NW/0685. All patients included were consented to participate in the study and to use their materials in research. All samples from Nottingham used in this study were pseudo-anonymised and stored in compliance with the UK Human Tissue Act. The study was performed in accordance with the Declaration of Helsinki.
Data availability
The authors confirm the data that have been used in this work are available on reasonable request.
Competing interests
The authors declare no competing interests.
Funding information
This research was supported and funded by Damietta University and Egyptian Ministry of Higher education and Scientific Research.
Additional information
Note This work is published under the standard license to publish agreement. After 12 months the work will become freely available and the license terms will switch to a Creative Commons Attribution 4.0 International (CC BY 4.0).
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Elsharawy, K.A., Mohammed, O.J., Aleskandarany, M.A. et al. The nucleolar-related protein Dyskerin pseudouridine synthase 1 (DKC1) predicts poor prognosis in breast cancer. Br J Cancer 123, 1543–1552 (2020). https://doi.org/10.1038/s41416-020-01045-7
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41416-020-01045-7
This article is cited by
-
The role of RNA-modifying proteins in renal cell carcinoma
Cell Death & Disease (2024)
-
A pan-cancer analysis of Dyskeratosis congenita 1 (DKC1) as a prognostic biomarker
Hereditas (2023)
-
Dyskerin and telomerase RNA component are sex-differentially associated with outcomes and Sunitinib response in patients with clear cell renal cell carcinoma
Biology of Sex Differences (2023)
-
SUMO specific peptidase 3 halts pancreatic ductal adenocarcinoma metastasis via deSUMOylating DKC1
Cell Death & Differentiation (2023)
-
The importance of pseudouridylation: human disorders related to the fifth nucleoside
Biologia Futura (2023)