Skip to main content
  • Research article
  • Open access
  • Published:

Hsa-miR-181a-5p, hsa-miR-182-5p, and hsa-miR-26a-5p as potential biomarkers for BCR-ABL1 among adult chronic myeloid leukemia treated with tyrosine kinase inhibitors at the molecular response

Abstract

Background

Tyrosine kinase inhibitors (TKIs) as first-line therapy for Chronic Myeloid Leukemia (CML) show a high success rate. However, a low number of patients with long-term treatment-free remission (TFR) were observed. Molecular relapse after imatinib discontinuation occurred at 50% at 24 months, with 80% occurrence within the first 6 months. One of the reasons for relapse is untimely TKIs discontinuation caused by large errors from estimates at very low-level or undetectable disease, thus warranting new biomarkers for CML.

Methods

Next Generation Sequencing (NGS) was used to identify microRNAs (miRNAs) at the molecular response in CML adult patients receiving TKIs treatment. A total of 86 samples were collected, 30 from CML patients responsive and 28 from non-responsive to imatinib therapy, and 28 from blood donors. NGS was conducted whereby 18 miRNAs were selected and validated by real-time RT-qPCR in triplicate.

Results

Hsa-miR-181a-5p was expressed significantly (p-value< 0.05) with 2.14 and 2.33-fold down-regulation in both patient groups, respectively meanwhile hsa-miR-182-5p and hsa-miR-26a-5p were significant only in the non-responsive group with 2.08 and 2.39 fold up-regulation. The down-regulation was consistent with decreased amounts of BCR-ABL1 in patients taking TKIs regardless of molecular responses. The up-regulation was consistent with the substantial presence of BCR-ABL1 in CML patients treated with TKIs at the molecular response.

Conclusions

Therefore, these miRNAs have potential as new therapeutic biomarkers for BCR-ABL1 status in adult CML patients treated with TKIs at molecular responses. These could improve current approaches and require further analysis to look for targets of these miRNAs in CML.

Peer Review reports

Background

CML occurs mainly in adults with Philadelphia chromosome (Ph) present in most patients [1] as a result of a reciprocal translocation between chromosome 9 and 22 that gives rise to BCR-ABL1 transcripts, t(9;22)(q34.1;q11.2) [2, 3]. There are 3 phases of CML: chronic, accelerated, and blast crisis; CML is treated using Tyrosine Kinase Inhibitors (TKIs). Imatinib (Glivec®, Novartis), a first-generation TKI, has been widely used since the Food and Drug Administration (FDA) approval in 2001, given orally (400 mg) once daily, and patients in the chronic phase usually respond well to treatment [1, 4, 5]. It is a targeted therapy whereby imatinib competes for the BCR-ABL1 tyrosine kinase site to prevent phosphorylation and inhibit proliferation [6]. Patients treated with imatinib are usually experiencing mild or moderate adverse effects [5].

Treatment progress is measured based on hematologic, cytogenetic, and molecular responses. Complete Hematologic Response (CHR) is when clinical presentations return to normal meanwhile, Complete Cytogenetic Response (CCyR) is when Ph is not detected. A molecular response is the major treatment endpoint with the optimal molecular response at BCR-ABL1 transcript level ≤ 10% by 3 months, ≤1% by 6 months is similar to CCyR and ≤ 0.1% by 12 months, known as Major Molecular Response (MMR). Beyond optimal response is known as warning or failure whereby each indicates treatment continuation, carefully considering or changing according to the patient’s conditions. Patients non-responsive or intolerant to imatinib will be administered second and third-generation (2G or 3G) TKIs as necessary, whereby nilotinib, dasatinib, and bosutinib, 2GTKIs are now choices of front-line treatment. These enable patients to not only achieve MMR faster but also much deeper molecular responses than imatinib. However, 2G and 3G TKIs showed harsher side-effects and clinical constraints without improved overall survival [4].

The BCR-ABL1 transcripts level quantification is by reverse transcription-quantitative polymerase chain reaction (RT-qPCR). It is a ratio of BCR-ABL1 transcripts to ABL1 transcripts or another internationally accepted control transcript regarding the International Randomized Study of Interferon and STI571 (IRIS) trial. It is referred to as the International Scale (IS). The IRIS standardized baseline is considered as 100% and below this or 1, 0.1, 0.01% (MR4), 0.0032% (MR4.5) and 0.001% (MR5) correspond to 2, 3, 4, 4.5, and 5-log reduction from the IRIS baseline respectively. MR4, MR4.5, and MR5 are known as deep molecular responses (DMR) [6]. A conversion factor converts laboratory ratio to an IS ratio and is determined through sample exchange with an established reference laboratory or calibrated to WHO International Genetic Reference Panel [7].

Patients who achieved DMR (MR4, MR4.5) have a quality of life near the general population; enable cessation studies in suitable patients. It is promising as losing MMR is uncommon after a year in treatment-free remission (TFR) otherwise indicates failure. Based on shared decision making, requirements for TKI discontinuation include first chronic phase, motivated patient, no prior treatment failure, typical e13a2 or e14a2 BCR-ABL1 transcripts, more than 5 years on TKI, and achieved DMR for more than 2 or 3 years for MR4.5 and MR4 respectively [4].

DMR of 4.5-log reduction or lower is preferable for discontinuation of imatinib by controlled studies due to fewer events observed and better clinical outcomes [8,9,10]. These low-level estimates are measured by using RT-qPCR which gives substantial errors, particularly from Poisson variation, with a limit of detection (LoD) of 10− 4, 10–4.5, and below [11,12,13]. These lead to disagreements in reporting as recommended by EUTOS, particularly in assigning targets to 3 for less than 3 copies [11, 12, 14]. Among suggestions for reporting of the result were to include confidence limit [12], average targets or as an individual result for undetectable replicate [11], and LoD was proposed [15]. In most centres, the routine LoD for RT-qPCR is MR4.5 [16]. All these have compromised the consistency of DMR which could hinder timely TKI discontinuation; therefore, a new therapeutic biomarker is necessary.

MicroRNAs (miRNAs) are small non-protein-coding, endogenous RNAs around 22 nucleotides long and function as regulators in gene expression. The expression could be down-regulated or up-regulated. Mature miRNAs are incorporated in the RNA-induced silencing complex and annealed to mRNA. Regulation of gene expression is by mRNA cleavage or translational repression [17]. miRNAs have been studied extensively in cancer, and studies showed miRNAs have potential as a tumour suppressor or oncogene with potential use in the diagnosis, prognosis, and treatment of cancers [18,19,20,21,22,23,24,25].

NGS is suitable for detecting low-level variants as it allows deep sequencing with high sensitivity and high sequencing coverage. There are several platforms available with unique specifications [26]. Miseq system (illumina®), an easy-to-use NGS instrument that offers short sequencing run times, long read lengths, and high data quality, is suitable for targeted resequencing. The sequencing workflow involves library preparation, cluster generation, sequencing by synthesis, and data analysis. Except for library preparation, the rest are carried out by the instrument making it a reliable platform.

In this study, we investigate miRNAs of adult CML patients, responsive and non-responsive to imatinib therapy at the molecular response. Imatinib is the front-line treatment used to treat CML, and in the event of intolerance or failure, other approved TKIs are an option as appropriate. Current treatment success is measured by the minimal presence or undetectability of BCR-ABL1 transcripts by RT-qPCR. miRNAs from the whole blood of these patients were profiled using NGS and compared with normal control. miRNAs suitable for the study were selected and verified by real-time RT-qPCR, the known gold standard for quantifying gene expression. Therefore, this study aimed to identify miRNAs that have therapeutic potential concerning TKIs at the molecular response in CML. These miRNAs could lead to discovering conditions that favour TKI discontinuation accurately, thus making long-term TFR achievable.

Methods

Patients and samples

A cross-sectional study was conducted among adults aged 18 and above in Ampang Hospital and National Blood Centre (NBC) Malaysia between 2013 and 2014. Of these respondents, 30 were CML patients responsive to imatinib therapy, and 28 were CML patients non-responsive to imatinib therapy taken from Ampang Hospital, while 28 were blood donors mostly from NBC. Peripheral blood samples were collected from all the respondents. A standardized data collection form was used to collect information from the respondents. The collection form has three domains that were sociodemographic information, treatment history with imatinib, and treatment status responsive or non-responsive. Responsive is defined as a CML patient assessed and confirmed by the clinician through clinical diagnostics and evaluations, receiving imatinib not less than 18 months, the latest BCR-ABL1IS ratio was 0.1% or less and ABL1 control gene equalled 10,000 copies or more. Non-responsive, defined as any of the above criteria that are not met and treated with imatinib, other TKIs, or other drugs as necessary. Control, defined as normal people taken from healthy donors, mostly was from the NBC who came to donate blood. The lab used a conversion factor of 0.81 to convert the BCR-ABL1 ratio to the IS ratio. Blood samples and clinical data were collected from the clinic and sent to the research laboratory.

Preparation of samples

Peripheral blood (2.5 ml) was collected from each participant using BD® Vacutainer Safety-Lok™ Blood Collection Set and drawn into PAXgene® Blood RNA tube (PreAnalytiX). Purification of total RNA including miRNA was conducted using; PAXgene® Blood miRNA Kit (PreAnalytiX) with silica-membrane technology. The quantity of total RNA was measured by NanoDrop® ND-1000 spectrophotometer. Meanwhile, the quality of total RNA and small RNA were analysed by the 2100 Bioanalyzer (Agilent) using Eukaryote Total RNA Nano assay (Agilent) for RNA Integrity Number (RIN) and Small RNA assay (Agilent) for the amount (%) of miRNA.

Next generation sequencing

Total RNA of a CML patient responsive to imatinib therapy with RIN 7.7, a CML patient non-responsive to imatinib therapy with RIN 7.4, and a normal control with RIN 8.1 were selected for NGS. TruSeq Small RNA kit (illumina®); was used to prepare libraries for subsequent cluster generation. A library was generated from each sample and labelled using a unique index sequence each. The kit contained adapters designed to directly and specifically ligate to a 5′-phosphate and a 3′-hydroxyl group of mature miRNAs. The ligated small RNA was reverse transcribed and amplified with a unique index per tube to create cDNA constructs according to the protocol (Catalogue # RS-930-1012, Part # 15004197 Rev. C, March 2011). The cycling conditions for PCR amplification were 30 s at 98 °C, 15 cycles of 10 s at 98 °C, 30 s at 60 °C, and 15 s at 72 °C, followed by 10 min at 72 °C and hold at 4 °C. Each sample was analysed using High Sensitivity DNA assay (Agilent) and obtained around 150 bp band and peak each at a 100x dilution.

The cDNA constructs (library) were purified using 6% Novex TBE PAGE Gel, 1.0 mm (Invitrogen). Libraries with unique indexes were pooled, added Novex DNA Loading Dye (Invitrogen), and loaded into 2 wells (25 μl each). The gel ran for 60 min at 145 V in the XCell SureLock™ Mini-Cell electrophoresis unit (Invitrogen™) containing 1X TBE buffer. Upon completion, the gel was removed and stained with 0.5 mg/ml ethidium bromide solution for 2 to 3 min. It was viewed under UV light, and bands between 160 and 145 bp were cut out using a razor blade. The excised gel fragment was passed through a 0.5 ml Gel Breaker tube by centrifugation, added 300 μl of Ultra-Pure water, passed through a 5 μm Filter tube by centrifugation, and concentrated by ethanol precipitation. The pellet was then resuspended in 10 μl 10 mM Tris-HCl, pH 8.5, and kept at − 150 to -25 °C. The purified sample library (different indexes) was validated using Agilent’s DNA 1000 assay (obtained a single sample peak around 150 bp) and High Sensitivity DNA assay (diluted 100x and obtained 152 bp peak, 1500 pmol/l) measured using the 2100 Bioanalyzer (Agilent).

The purified sample library was normalized to a 4 nM library using 10 mM Tris-Cl pH 8.5 with 0.1% tween 20, denatured to a 20 pM library in 1 mM NaOH using 0.2 N NaOH and pre-chilled HT1, and diluted to 13 pM using pre-chilled HT1. These were prepared according to Preparing Libraries for Sequencing on the Miseq® guide: Part # 15039740 Rev.D October 2013 (illumina®). The supplied 10 nM PhiX library was diluted to 4 nM using 10 mM Tris-Cl pH 8.5 with 0.1% tween 20, denatured to 20 pM, and further diluted to 12.5 pM. The denatured and diluted PhiX control was spiked in at least 5% to the purified sample library. The libraries were then dispensed into sample reservoir of reagent cartridge according to MiSeq® Reagent Kit v2 (50 cycles) Reagent Preparation Guide (illumina®) and ran on Miseq (illumina, USA) using Standard Flow Cell with 14 tiles imaged top and bottom for 15 M reads (approximately 5 M reads per sample). Clusters were generated and sequenced by synthesis with a single read, and then were imaged using LED and filter combinations meant for each fluorescently-labelled nucleotide. BaseSpace (illumina®) aligned, assembled, and analysed the reads. miRNAs of the responsive sample and miRNAs of the non-responsive sample were compared to miRNAs of the normal control. Eighteen miRNAs were selected mainly with regards to substantial expression level and clinical significance.

Real-time reverse transcriptase-quantitative PCR (RT-qPCR)

The selected miRNAs were ordered as ready-to-use 96-well plate Custom miScript™ miRNA PCR Array (Qiagen). The array plate consisted of the 18 miRNAs, normalisation controls (SNORD61, SNORD95, SNORD96A, RNU6–2), reverse transcription control (miRTC), and positive PCR control (PPC) in the 24 × 4 format. cDNA was prepared using miScript® II RT kit (Qiagen) for a 20 μl reaction. The reaction consisted of 4 μl 5x miScript HiSpec Buffer, 2 μl 10x miScript Nucleic Mix, 2 μl miScript Reverse Transcriptase Mix, RNase-free water, and template RNA. A starting input of 250 ng total RNA (containing miRNA) was used as recommended, which will result in 0.5 to 1 ng cDNA per array well. The reaction was incubated for 60 min at 37 °C, followed by 5 min at 95 °C, and placed on ice. The 5x miScript HiSpec buffer facilitated the conversion of mature miRNAs into cDNA while suppressing long RNAs.

miScript® SYBR® Green PCR kit (Qiagen); was used with our 96-well plate Custom miScript™ miRNA PCR Array. A reaction mixture consisted of 1375 μl 2x QuantiTect SYBR® Green PCR Master Mix, 275 μl 10x miScript Universal Primer, 1000 μl RNase-free water, and 100 μl template cDNA. Out of this, 25 μl per well was added to the verified array plate. Then the array plate was placed in Light Cycler 480 (Roche), a real-time cycler. The cycling conditions used were; initial activation for 15 min at 95 °C, 45 cycles of; denaturation for 15 s at 94 °C followed by annealing for 30 s at 55 °C, and extension for 30 s at 70 °C. PCR amplification was followed by dissociation curve analysis performed for 10 s at 95 °C, 1 min at 65 °C, and a continuous 95 °C. The miScript PCR system used; allowed sensitive and specific detection and quantification of the mature miRNAs in triplicate in all case and control samples.

A threshold value of 1.15 was used across all PCR runs and obtained threshold cycle (CT) of PPC (CTPPC) value equivalent to 19 + 2 and CTmiRTC- CTPPC value less than 7, thus allowing comparison of results. These indicated high-quality RNA samples, correct cycling program run, and the correct threshold was used, and no apparent inhibition of the reverse-transcription reaction occurred, respectively. In addition, a single melting peak was exhibited in each sample representing a specific PCR product. The Light Cycler uses a crossing point with the Fit Points method while Gene Globe Data Analysis Centre (Qiagen) uses CT, which is similar. Data was formatted in an Excel Spreadsheet using a provided template and uploaded into the Gene Globe Data Analysis Centre in a web browser (QIagen). Normalisation was conducted using the Global Threshold Cycle Mean of expressed miRNAs. Relative quantification was performed using ∆∆CT method; whereby ∆CT = CTmiRNA - CT normalisation factor, while ∆∆CT = ∆CT (experimental sample) – ∆CT (control sample). Fold change was calculated as 2-∆∆CT whereby greater than 1 was reported as fold up-regulation while less than 1 as fold down-regulation. P-values were calculated based on the Student’s t-test of the replicate 2(−∆CT) values for each miRNA. Any significant miRNA observed will be linked to TKIs received and the clinical outcomes of patients.

Results

The majority of the respondents were Malays (51%), followed by Chinese (34%) and Indians (15%). There were more males than females in this study. The hallmark of CML (BCR-ABL1 transcript) was present in all patients who were diagnosed with CML. In responsive patients, BCR-ABL1 ratios showed a declining pattern with imatinib therapy. Common chromosomal abnormalities; usually seen with CML were observed such as; trisomy 8, which was seen in responsive and non-responsive patients. Loss of chromosome Y, an age-related chromosomal abnormality, was seen in Ph-negative older male patients. This study showed that the number of responsive and non-responsive patients was similar in Malays. The responsiveness is markedly higher among the Chinese. However, the numbers of Indian respondents were too small for comparison (Table 1).

Table 1 Demographic characteristic of CML patients responsive and non-responsive to imatinib therapy and blood donors

In responsive patients, 53% achieved MMR of 7 Malays, 6 Chinese, and 3 Indians. The remaining 47% achieved DMR of 6 Malays, 6 Chinese, and 2 Indians, in which 40% were with 4-log reduction (MR4) and 7% with 4.5-log reduction (MR4.5) from the IRIS baseline. Longer duration on imatinib of around 6 years (mean 71.64 ± 34.50 months) was observed in DMR patients compared to about 5.5 years (mean 66 ± 30.43 months) in MMR patients. Overall mean duration on imatinib was 5.7 years with 53% being on imatinib for more than 5 years (Table 2). In non-responsive patients, imatinib was changed to other TKIs (82%) or other treatment options (11%) according to mutations in the kinase domain and patients’ conditions and tolerability. BCR-ABL1IS ratio was < 0.1% in 29% of patients, and the average of the latest available BCR-ABL1IS ratio was 13.7% (Table 3).

Table 2 Clinical findings of CML patients responsive to imatinib therapy and a potential biomarker at the molecular response
Table 3 Clinical findings of CML patients non-responsive to imatinib therapy and potential biomarkers at the molecular response

Profiles of miRNA expression by NGS showed a similar amount of mature miRNA was present in all samples, with 40.8% in normal control, 40.9% in responsive sample, and 40.3% in the non-responsive sample. A list of miRNAs with various expression levels was generated from each patient sample and compared to normal control. Among these, hsa-miR-181-5p expressed 0.45 and 0.51 times in responsive and non-responsive samples respectively.

These were consistent with real-time RT-qPCR data whereby significant down-regulations of hsa-miR-181a-5p were observed in both responsive (Group 1) and non-responsive (Group 2) groups when compared to the control group (Figs. 1 and 2 respectively). Fold change of hsa-miR-181a-5p was 0.47 (95% CI = 0.36, 0.58, p-value = 0.000006 that was < 0.05) in responsive group (Fig. 3) and 0.43 (95% CI = 0.31, 0.55, p-value = 0.000024 that was < 0.05) in non-responsive group (Fig. 4) which corresponded to 2.14 and 2.33 fold down-regulation respectively compared to control group.

Fig. 1
figure 1

Significant down-regulation of hsa-miR-181a-5p and RNU6-2 as shown by arrows when Group 1 (CML patients responsive to imatinib therapy); was compared with the Control Group (blood donors) using custom array plate and real-time RT-qPCR. The custom array plate consisted of 18 miRNAs selected from NGS. RNU6-2 is a reference gene

Fig. 2
figure 2

Significant down-regulation of hsa-miR-181a-5p as shown by an arrow when Group 2 (CML patients non-responsive to imatinib therapy); was compared with the Control Group (blood donors). In addition, the other two arrows showed significant up-regulation of hsa-miR-26a-5p and hsa-miR-182-5p

Fig. 3
figure 3

Lower expression of hsa-miR-181a-5p of 0.47 was observed in Group 1 (CML patients responsive to imatinib therapy) compared with the Control Group (blood donors) by real-time RT-qPCR (p-value < 0.05). RNU6–2 showed a similar expression

Fig. 4
figure 4

Lower expression of hsa-miR-181a-5p of 0.43 was also observed in Group 2 (CML patients non-responsive to imatinib therapy) compared with the Control Group (blood donors) by real-time RT-qPCR (p-value < 0.05). In addition, higher expressions were only observed from hsa-miR-182-5p and hsa-miR-26a-5p of 2.08 and 2.39, respectively, in this group

In addition, another 2 miRNAs were also significant but only in the non-responsive group. These were hsa-miR-182-5p with fold change of 2.08 (95% CI = 1.60, 2.57, p-value = 0.000012 that was < 0.05) and hsa-miR-26a-5p with fold change of 2.39 (95% CI = 1.82, 2.96, p-value = 0 that was < 0.05) as shown in Fig. 4 which corresponded to 2.08 and 2.39 fold up-regulation respectively when compared to the control group. The statistical analyses of real-time PCR data were as shown in Tables 4 and 5. Comparison of these biomarkers to BCR-ABL1 and IS levels were presented in Tables 2 and 3.

Table 4 Normalized expressions (2^(−Avg.(Delta(Ct)) of hsa-miR-181a-5p and RNU6–2 in Group 1 (CML patients responsive to imatinib therapy), and Control Group (blood donors) calculated by Gene Globe Data Analysis Centre (QIAGEN)
Table 5 Normalized miRNA expression (2^(−Avg.(Delta(Ct)) in Group 2 (CML patients non-responsive to imatinib therapy) and Control Group (blood donors) of significant miRNAs calculated by Gene Globe Data Analysis Centre (QIAGEN)

In the responsive group, a reference gene, RNU6–2, a small non-coding RNA, was also significant with fold change of 0.48 (95% CI = 0.37, 0.58, p-value = 0 that was < 0.05) corresponded to 2.10 fold down-regulation when compared to the control group (Fig. 3).

Discussion

We looked at our CML patients concerning imatinib therapy, a front-line treatment for CML in our government hospitals; because of treatment discontinuation for optimal CML responders as suggested by the European LeukemiaNet (ELN) 2013 and also ELN 2020. Data showed that half of our responsive patients received imatinib therapy for more than 5 years with a mean age of 48 years old, an average of 5.5 years for patients with MMR, and 6 years for DMR. These indicated on average, MMR achieved by 18 months extended survival for at least 4 years, and DMR is achievable with a longer time on therapy in our CML patients. DMR was observed in almost half of the responsive patients, and if consistent, these patients can be considered for imatinib discontinuation studies, as minimum requirements are at least 24 months of DMR [4], and previously 12 months of DMR [27]. DMR was reported suitable for imatinib discontinuation as it has survival benefits and lesser events [8]. However, determining the right time for imatinib discontinuation is critical. Long-term TFR is still low, with relapse after discontinuation being considerably high in remission patients. The weighted mean molecular relapse rate of CML was 42%, and mean molecular remission was more than 37 months before discontinuation of imatinib. A shorter duration of imatinib therapy and a shorter time with an undetectable level of BCR-ABL1 transcripts were reasons for relapse [28]. Therefore sufficient time on imatinib and consistency in DMR is crucial before imatinib discontinuation.

Management of CML has remarkably improved over the years with a high success rate for treatment using TKIs in chronic patients. This is assisted by new analysis techniques complementing existing ones for comprehensive detection and monitoring of CML. RT-qPCR, a sensitive technique has become a gold standard in the molecular response and is used to analyse treatment response, especially from TKIs as the first line. Since RT-qPCR has limitations at very low-levels thus new biomarkers are needed not only to assist analysis but also to understand better molecular and cellular interactions and associations to prolong TFR.

Numerous novel findings; were reported from CML studies using NGS, a new technology with relatively high sensitivity and high coverage. In CML patients treated with imatinib, NGS has enabled identification of novel BCR-ABL1 fusions gene comprising BCR intron 14 and ABL1 intron 2 breakpoints, giving partial deletion of SH3 domain [29]. Meanwhile, with whole-genome sequencing, NGS has enabled identification of e13a2-like BCR-ABL1 fusion comprising broken BCR exon 13 and 9 ABL1 intron 1 nucleotides forming a novel chimeric exon [30]. In our study, NGS has assisted in identifying miRNAs at the molecular response in local CML patients responsive and non-responsive to imatinib therapy. Validation by custom array real-time RT-qPCR has identified hsa-miR 181-5p, hsa-miR-182-5p, and hsa-miR-26a-5p as potential new therapeutic biomarkers in response to treatment using TKIs in CML. Hsa-miR 181-5p could be a tumour suppressor miRNA and might be involved in major or partial suppression of CML. Meanwhile, hsa-miR-182-5p and hsa-miR-26a-5p could be oncomiRs in which their expressions may indicate non-responsiveness to TKIs treatment or relapse.

Hsa-miR-181a-5p was significantly down-regulated in both the responsive and non-responsive to imatinib groups compared with the control group. The down-regulation indicated significant interactions between hsa-miR-181a-5p and TKI therapies. These showed that imatinib, nilotinib, dasatinib, and bosutinib were taken by CML patients, over time had caused inhibition to CML cell proliferation, led to apoptosis, and thus reduced the amount of BCR-ABL1 in CML patients. The down-regulation observed in the responsive to imatinib group represented 3 to 4.5-log reduction (MMR to MR4.5) of BCR-ABL1 levels from the IRIS baseline with an average of 5.7 years of intake. Similarly, down-regulation observed in the non-responsive to imatinib group represented less than 1 to 5-log reduction (approaching MMR to MR5) of BCR-ABL1 levels from the IRIS baseline. These down-regulations of hsa-miR-181a-5p indicated a positive correlation to log reduction or reduced levels of BCR-ABL1 in CML patients treated with TKIs. Therefore hsa-miR-181a-5p has the potential to be used as a therapeutic biomarker for the positive performance of TKIs shown by reduced levels of BCR-ABL1 indicating patients were responding to TKIs treatment in CML. Hsa-miR-181a-5p could be a tumour-suppressive miRNA that prevents mRNA from coding specific protein directly or associated with CML in patients treated with TKIs.

On the other hand, hsa-miR-182-5p and hsa-miR-26a-5p were significantly up-regulated in the non-responsive to imatinib group compared with the control group. The fold up-regulation was 2.08 and 2.39, respectively, and was observed only in this group. These were patients whereby most of their latest available BCR-ABL1IS ratios were more than 0.1% that was neither in MMR nor did DMR, with a substantial average, and only 7% [2] treated with imatinib. Fold regulation of these miRNAs in the imatinib responsive group, whereby BCR-ABL1 in these patients was minimal or absent (MMR or DMR), was not significant compared with the control group. It indicated imatinib had reduced the amount of BCR-ABL1 in these responsive patients to the extent not significantly different compared with healthy donors. Thus the significant up-regulation could indicate that hsa-miR-182-5p and hsa-miR-26a-5p have a positive correlation to the substantial presence of BCR-ABL1 in CML patients treated with TKIs. Therefore, these miRNAs could function as new therapeutic biomarkers for the substantial presence of BCR-ABL1 in CML patients treated with TKIs. Hsa-miR-182-5p and hsa-miR-26a-5p could be oncogenic miRNAs (oncomiRNAs) and might be associated with CML in patients treated with TKIs.

Studies have shown hsa-miR-181a-5p as a therapeutic biomarker, has potential use for clinical improvement in cancers. It is demonstrated in gastric cancer cells whereby hsa-miR-181a-5p has been extensively studied and reported as a potential regulator of MEG2, a tumour suppressor gene [31], regulated RASSF6, and in combination predicts poor prognosis in gastric cancer [32]. It also inhibited MTMR3 expression in AGS gastric cancer cells and was identified as novel autophagy [33]. In non-small-cell-lung cancer tissues and cell lines, hsa-miR-181a-5p was reported significantly reduced and has a potential role in tumour suppression by partially targeting Kras [34], a protein that regulates cell growth. In breast (BC) and colon cancers (CC), hsa-miR-181a-5p was reportedly down-regulated and inversely related to matrix metalloproteinase-14, which is elevated in tumours, in which, to prevent cancer metastasis in BC and CC, is by elevating hsa-miR-181a-5p [35]. In hepatocellular carcinoma: miR-181 was reported to significantly turn on the MAPK/JNK pathway, the regulator of cell proliferation, and by limiting it, would suppress the pathway [36]. Hsa-miR-181a-5p, a tumour suppressor is considered as a prognostic marker in Acute Myeloid Leukemia patients treated with intensive induction chemotherapy and autologous stem cell transplant [37].

Higher expression of hsa-miR-182-5p as observed in our study was also seen in other studies as in prostate cancer (PCa); lung squamous cell carcinoma tissues analysed from The Cancer Genome Atlas database, the Gene Expression Omnibus database, and real-time qPCR [38]; and also in colorectal cancer cells using tumorigenic variant cell line MICOL-14tum compared to MICOL-14h-tert cells [39]. In PCa, the biphasic role of hsa-miR-182-5p was observed, with a higher expression in localized PCa and contrarily lower expression in aggressive cancers [40]. Lower expression was also observed in renal cancer whereby hsa-miR-182-5p was down-regulated in tumour tissue compared with adjacent normal tissues, and overexpression decreased tumour growth in mice, demonstrating antitumor effect [41].

Hsa-miR-26a-5p expression in our study is consistent with a bigger-sample size study than ours on CML patients indicating the expression of miR-26a is influenced by the response to TKIs used [42]. The study showed miR-26a was differentially expressed; in HL-60.BCR-ABL cells when treated with tyrosine kinase inhibitors when compared to HL-60 cells. miR-26a expression levels increased gradually with both dasatinib and nilotinib treatments. The increased expression levels with dasatinib and nilotinib were consistent with our findings of hsa-miR-26a expression was up-regulated in CML patients molecularly non-responsive to imatinib compared with controls whereby most patients were treated with nilotinib or dasatinib. Thus, with dasatinib and nilotinib treatments, up-regulation of hsa-miR-26a was observed in CML patients’ non-responsiveness to imatinib. However, HL-60.BCR-ABL cells treated with imatinib mesylate showed increased miR-26a expression at 4-h but decreased expression at 8-h. It correlates with miR-26a expression was lower in chronic phase patients (62.5% Complete Cytogenetic Response (CCyR) and 37.5% imatinib-resistant) when compared with healthy individuals. miR-26a expression was also lower in imatinib-resistant patients when compared to patients who achieved CCyR (< 1% BCR-ABL1) [43]. Hsa-miR-26a-5p expression in our study was not significant in CML patients molecularly responsive to imatinib compared with controls as these were patients either in MMR or DMR (≤ 0.1% to undetectable BCR-ABL1). These showed that with imatinib, expressions of hsa-miR-26a were decreasing from responsiveness towards resistance. Thus these expressions indicated that imatinib, nilotinib, and dasatinib had altered BCR-ABL kinase activity differently.

The main limitation of this study was the low number of patient samples collected, which was due to only a small number of patients fulfilling all the required criteria within the 2-year duration. Hence, a longer time of sample collection is appropriate. In addition, our custom plate of 4 × 24 introduced variations between plates for replicates on different array plates. It made determining a suitable threshold challenging. Instead, using an array plate of 3 × 32 for triplicate is recommended.

These studies showed that hsa-miR-181a-5p, hsa-miR-182-5p, and hsa-miR-26a-5p are promising therapeutic biomarkers and can be manipulated to improve cancer treatment. In this study, only 18 miRNAs from NGS profiles were validated using real-time RT-qPCR; thus the significant miRNAs observed could be solely or partially involved in the decreasing and or increasing of BCR-ABL1 in CML patients. Further studies are needed to determine their targets and the effect of involvement in CML. Knowing their interactions and targets in CML will give their net effects that allow manipulation by inhibiting oncomiRs or stimulating tumour suppressor miRNAs. These will give new values for TKIs responses which could assist in determining the accurate time for TKIs discontinuation and better prediction of relapse in CML. Timely discontinuation of TKIs in CML remission patients is crucial for long-term TFR. Prolong TFR would subsequently allow CML patients to achieve a good quality of life as healthy people.

Conclusions

Hsa-miR-181a-5p has the potential as a therapeutic biomarker for BCR-ABL1 at below current LoD in CML adult patients treated with TKIs at the molecular response and could be a tumour-suppressive miRNA. Meanwhile, hsa-miR-182-5p and hsa-miR-26a-5p have potential as therapeutic biomarkers for the substantial presence of BCR-ABL1 in CML adults and could be oncomiRNAs.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

TKI:

Tyrosine Kinase Inhibitor

CML:

Chronic Myeloid Leukemia

TFR:

Treatment free remission

NGS:

Next Generation Sequencing

RT-qPCR:

Reverse transcriptase-quantitative PCR

miRNA:

MicroRNA

Ph:

Philadelphia chromosome

IRIS:

International Randomized Study of Interferon and STI571

IS:

International Scale

DMR:

Deep Molecular Response

CCyR:

Complete Cytogenetic Response

LoD:

Limit of Detection

RIN:

RNA Integrity Number

PPC:

Positive PCR Control

miRTC:

Reverse Transcription Control

CT :

Threshold Cycle

ELN:

European LeukemiaNet

BC:

Breast Cancer

CC:

Colon Cancer

PCa:

Prostate Cancer

References

  1. Alikian M, Gale RP, Apperley JF, Foroni L. Molecular techniques for the personalised management of patients with chronic myeloid leukaemia. Biomol Detect Quantif. 2017;11:4–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Heisterkamp N, Stam K, Groffen J, de Klein A, Grosveld G. Structural organization of the bcr gene and its role in the Ph’ translocation. Nature. 1985;315(6022):758–61.

    Article  CAS  PubMed  Google Scholar 

  3. Canaani E, Steiner-Saltz D, Aghai E, Gale RP, Berrebi A, Januszewicz E. Altered transcription of an oncogene in chronic myeloid Leukaemia. Lancet. 1984;323(8377):593–5.

    Article  Google Scholar 

  4. Hocchaus A. European LeukemiaNet 2020 recommendations for treating chronic myeloid leukemia. Leukemia. 2020;34:966–84.

    Article  Google Scholar 

  5. Pretel-Irazabal M, Tuneu-Valls A, Ormaechea-Pérez N. Adverse skin effects of Imatinib, a tyrosine kinase inhibitor. Actas Dermosifiliográficas. 2014;105(7):655–62.

    Article  CAS  Google Scholar 

  6. Druker BJ, Lydon NB. Lessons learned from the development of an abl tyrosine kinase inhibitor for chronic myelogenous leukemia. J Clin Invest. 2000;105(1):3–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Cross NCP, White HE, Colomer D, Ehrencrona H, Foroni L, Gottardi E, et al. Laboratory recommendations for scoring deep molecular responses following treatment for chronic myeloid leukemia. Leukemia. 2015;29(5):999–1003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Kaygusuz Atagunduz I, Toptas T, Deniz R, Kara O, Eser A, Sezgin A, et al. Effects of deeper molecular responses on outcomes in chronic myeloid leukemia patients in chronic phase treated with Imatinib Mesylate. Clin Lymphoma Myeloma Leuk. 2017;17(2):120–5.

    Article  PubMed  Google Scholar 

  9. Etienne G, Dulucq S, Nicolini F-E, Morisset S, Fort M-P, Schmitt A, et al. Achieving deeper molecular response is associated with a better clinical outcome in chronic myeloid leukemia patients on imatinib front-line therapy. Haematologica. 2014;99(3):458–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Rousselot P, Charbonnier A, Cony-Makhoul P, Agape P, Nicolini FE, Varet B, et al. Loss of major molecular response as a trigger for restarting tyrosine kinase inhibitor therapy in patients with chronic-phase chronic myelogenous leukemia who have stopped imatinib after durable undetectable disease. J Clin Oncol. 2014;32(5):424–30.

    Article  CAS  PubMed  Google Scholar 

  11. Young SS, Tucker T, Bosdet I, Karsan A. Estimating deep molecular responses in chronic myelogenous leukemia: a Bayesian approach. Leukemia. 2016;30(7):1628–9.

    Article  CAS  PubMed  Google Scholar 

  12. Morley AA. Reporting results for deep molecular responses in chronic myeloid leukemia. Leukemia. 2016;30(7):1630–1.

    Article  CAS  PubMed  Google Scholar 

  13. Latham S, Bartley PA, Budgen B, Ross DM, Hughes E, Branford S, et al. BCR-ABL1 expression, RT-qPCR and treatment decisions in chronic myeloid leukaemia. J Clin Pathol. 2016;69(9):817–21.

    Article  CAS  PubMed  Google Scholar 

  14. Ip H-W, So C-C. Deep molecular response in chronic myelogenous leukemia: ensuring accuracy and consistency. Leukemia. 2015;29(7):1620–1.

    Article  CAS  PubMed  Google Scholar 

  15. Cross NCP, Müller MC, Hochhaus A. Response to ho-wan Ip and chi-Chiu. Leukemia. 2015;29(7):1619.

    Article  CAS  PubMed  Google Scholar 

  16. Cross NC, Hochhaus A. Deep molecular response in chronic myeloid leukemia. Leukemia. 2016;30(7):1632.

    Article  CAS  PubMed  Google Scholar 

  17. Bartel DP. MicroRNAs. Cell. 2004;116(2):281–97.

    Article  CAS  PubMed  Google Scholar 

  18. Ma L, Reinhardt F, Pan E, Soutschek J, Bhat B, Marcusson EG, et al. Therapeutic silencing of miR-10b inhibits metastasis in a mouse mammary tumor model. Nat Biotechnol. 2010;28(4):341–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Jang E, Kim E, Son H-Y, Lim E-K, Lee H, Choi Y, et al. Nanovesicle-mediated systemic delivery of microRNA-34a for CD44 overexpressing gastric cancer stem cell therapy. Biomaterials. 2016;105:12–24.

    Article  CAS  PubMed  Google Scholar 

  20. Dong J, Liu Y, Liao W, Liu R, Shi P, Wang L. miRNA-223 is a potential diagnostic and prognostic marker for osteosarcoma. J Bone Oncol. 2016;5(2):74–9.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Chen S, Sun K-X, Liu B-L, Zong Z-H, Zhao Y. MicroRNA-505 functions as a tumor suppressor in endometrial cancer by targeting TGF-$α$. Mol Cancer. 2016;15(1):11.

  22. Wei Z, Liu Y, Wang Y, Zhang Y, Luo Q, Man X, et al. Downregulation of Foxo3 and TRIM31 by miR-551b in side population promotes cell proliferation, invasion, and drug resistance of ovarian cancer. Med Oncol. 2016;33(11):126.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Sato T, Shiba-Ishii A, Kim Y, Dai T, Husni RE, Hong J, et al. miR-3941: a novel microRNA that controls IGBP1 expression and is associated with malignant progression of lung adenocarcinoma. Cancer Sci. 2017;108(3):536–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. El Fatimy R, Subramanian S, Uhlmann EJ, Krichevsky AM. Genome editing reveals Glioblastoma addiction to MicroRNA-10b. Mol Ther. 2017;25(2):368–78.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Bandi V, Baluchamy S. miR-181a-2 downregulates the E3 ubiquitin ligase CUL4A transcript and promotes cell proliferation. Med Oncol. 2017;34(8):146.

    Article  PubMed  Google Scholar 

  26. Quail M, Smith ME, Coupland P, Otto TD, Harris SR, Connor TR, et al. A tale of three next generation sequencing platforms: comparison of ion torrent, pacific biosciences and illumina MiSeq sequencers. BMC Genomics. 2012;13(1):341.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Hughes TP, Ross DM. Moving treatment-free remission into mainstream clinical practice in CML. Blood. 2016;128(1):17–23.

    Article  CAS  PubMed  Google Scholar 

  28. Campiotti L, Basilio M, Guasti L, Piazza R, Gambacorti-passerini C, Maria A. Imatinib discontinuation in chronic myeloid leukaemia patients with undetectable BCR-ABL transcript level : a systematic review and a meta-analysis. Eur J Cancer. 2017;77:48–56.

    Article  CAS  PubMed  Google Scholar 

  29. Zhou F, Jin R, Hu Y, Mei H. A novel BCR-ABL1 fusion gene with genetic heterogeneity indicates a good prognosis in a chronic myeloid leukemia case. Mol Cytogenet. 2017;10(1):19.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Fu S, Hu Y, Fu Y, Chen F, Liu X, Zhang M, et al. Novel BCR-ABL1 fusion and leukemic mutations of SETBP1, PAX5, and TP53 detected by next generation sequencing in chronic myeloid leukemia. Cancer Biol Ther. 2016;17(10):1003–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Liu Z, Sun F, Hong Y, Liu Y, Fen M, Yin K, et al. MEG2 is regulated by miR-181a-5p and functions as a tumour suppressor gene to suppress the proliferation and migration of gastric cancer cells. Mol Cancer. 2017;16(1):133.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Mi Y, Zhang D, Jiang W, Weng J, Zhou C, Huang K, et al. miR-181a-5p promotes the progression of gastric cancer via RASSF6-mediated MAPK signalling activation. Cancer Lett. 2017;389:11–22.

    Article  CAS  PubMed  Google Scholar 

  33. Lin Y, Zhao J, Wang H, Cao J, Nie Y. miR-181a modulates proliferation, migration and autophagy in AGS gastric cancer cells and downregulates MTMR3. Mol Med Rep. 2017;15(5):2451–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Ma Z, Qiu X, Wang D, Li Y, Zhang B, Yuan T, et al. MiR-181a-5p inhibits cell proliferation and migration by targeting Kras in non-small cell lung cancer A549 cells. Acta Biochim Biophys Sin. 2015;47(8):630–8.

    Article  CAS  PubMed  Google Scholar 

  35. Li Y, Kuscu C, Banach A, Zhang Q, Pulkoski-Gross A, Kim D, et al. miR-181a-5p inhibits Cancer cell migration and angiogenesis via Downregulation of matrix Metalloproteinase-14. Cancer Res. 2015;75(13):2674–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Tan JYL, Habib NA, Chuah YW, Yau YH, Geifman-Shochat S, Chen WN. Identification of cellular targets of MicroRNA-181a in HepG2 cells: a new approach for functional analysis of MicroRNAs. PLoS One. 2015;10(4):1–17.

    Article  Google Scholar 

  37. Seipel K, Messerli C, Wiedemann G, Bacher U, Pabst T. MN1, FOXP1 and hsa-miR-181a-5p as prognostic markers in acute myeloid leukemia patients treated with intensive induction chemotherapy and autologous stem cell transplantation. Leuk Res. 2020;89:106296.

    Article  CAS  PubMed  Google Scholar 

  38. Luo J, Shi K, Yin S, Tang R, Chen W, Huang L, et al. Clinical value of miR-182-5p in lung squamous cell carcinoma: a study combining data from TCGA, GEO, and RT-qPCR validation. World J Surg Oncol. 2018;16(1):1–16.

    Article  Google Scholar 

  39. Perilli L, Tessarollo S, Albertoni L, Curtarello M, Pastò A, Brunetti E, et al. Silencing of miR-182 is associated with modulation of tumorigenesis through apoptosis induction in an experimental model of colorectal cancer. BMC Cancer. 2019;19(1):1–13.

    Article  CAS  Google Scholar 

  40. Baumann B, Acosta AM, Richards Z, Deaton R, Sapatynska A, Murphy A, et al. Association of High miR-182 levels with low-risk prostate Cancer. Am J Pathol. 2019;189(4):911–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Kulkarni P, Dasgupta P, Bhat NS, Shahryari V, Shiina M, Hashimoto Y, et al. Elevated miR-182-5p associates with renal Cancer cell mitotic arrest through diminished MALAT-1 expression. Mol Cancer Res. 2018;16(11):1750–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Ferreira AF, Moura LG, Tojal I, Ambrósio L, Pinto-Simões B, Hamerschlak N, et al. ApoptomiRs expression modulated by BCR–ABL is linked to CML progression and imatinib resistance. Blood Cells Mol Dis. 2014;53(1–2):47–55.

    Article  CAS  PubMed  Google Scholar 

  43. Ferreira AF, Moura LG, Tojal I, Ambrósio L, Pinto-Simões B, Hamerschlak N, et al. Corrigendum to “ApoptomiRs expression modulated by BCR-ABL is linked to CML progression and imatinib resistance” [Blood Cells, Mol Dis. 53, (1–2), (June–August 2014) 47–55]. Blood Cells Mol Dis. 2015;55(4):420.

    Article  Google Scholar 

Download references

Acknowledgments

We would like to thank the Director-General of Health Malaysia for his permission to publish this article, the Director of Institute for Medical Research for his technical support, respondents and assisting staff from Hospital Ampang and National Blood Centre, Malaysia; Asma Shahirah Azali and Muhammad Asyraf Noh for their comprehensive samples and clinical data collection.

Funding

This work was supported by a Research Grant from the Ministry of Health, Malaysia (NMRR 12–926-13248). The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

AMY was actively involved in all aspects of the study except for the clinical selection of patients, epidemiology, and biostatistics analysis, and was a major contributor in writing the manuscript, NAM involved in epidemiology and biostatistics research designed and critically reviewed the manuscript, CKM selected suitable CML patients for the study, HAH contributed in NGS preparation and reviewed the manuscript, YMY reviewed hematopathology analysis of the manuscript, and LI reviewed molecular analysis of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Aliza Mohd Yacob.

Ethics declarations

Ethics approval and consent to participate

All procedures performed involving human participants were approved by the Medical Research and Ethics Committee (Ministry of Health Malaysia, P12–837), in accordance with the ethical standards of the institutional and or national research committee, and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Written informed consent was obtained from all the respondents in the study before blood collection.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access 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/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mohd Yacob, A., Muhamad, N.A., Chang, K.M. et al. Hsa-miR-181a-5p, hsa-miR-182-5p, and hsa-miR-26a-5p as potential biomarkers for BCR-ABL1 among adult chronic myeloid leukemia treated with tyrosine kinase inhibitors at the molecular response. BMC Cancer 22, 332 (2022). https://doi.org/10.1186/s12885-022-09396-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-022-09396-5

Keywords