Skip to main content

Circulating microRNAs as minimal residual disease biomarkers in childhood acute lymphoblastic leukemia

Abstract

Background

Treatment stratification based on bone marrow minimal residual disease (MRD) at set time points has resulted in considerably improved survival in pediatric acute lymphoblastic leukemia (ALL). Treatment response is assessed using bone marrow samples. MicroRNAs (miRs) easily traffic among fluid spaces and are more stable than most other RNA classes. We examined the role of circulating miRs as putative less invasive MRD biomarkers.

Methods

In an exploratory experiment, expression of 46 preselected miRs was studied in platelet-free blood plasma samples of 15 de novo, 5 relapsed ALL patients and 10 controls by Custom TaqMan Array Advanced MicroRNA Card. Based on their high expression in ALL compared to controls, and on the reduction observed along the induction therapy, four miRs were selected for further analyses: miR-128-3p, -181a-5p, -181b-5p and 222-3p. Their expression was measured by qPCR at 4 time points in 27 de novo ALL patients treated in the ALL IC-BFM 2009 study.

Results

The expression of all 4 miRs significantly decreased over the first week of therapy (miR-128-3p: log2 fold change − 2.86; adjusted p 3.6 × 10−7; miR-181b-5p: log2 fold change − 1.75; adjusted p 1.48 × 10−2; miR-181a-5p: log2 fold change -1.33; adjusted p 3.12 × 10−2; miR-222-3p: log2 fold change − 1.25; adjusted p 1.66 × 10−2). However, no significant further reduction in miR expression was found after the 8th day of therapy. Measured drop in expression of 2 miRs at day 8 strongly correlated with day 15 bone marrow flow cytometry MRD results (miR-128-3p: Pearson’s r = 0.88, adjusted p = 2.71 × 10−4; miR-222-3p: r = 0.81, adjusted p = 2.99 × 10−3).

Conclusion

In conclusion, these circulating miRs might act as biomarkers of residual leukemia. MiR-128-3p and miR-222-3p in blood predict day 15 flow cytometry MRD results 7 days earlier. Although, their sensitivity falls behind that of bone marrow flow cytometry MRD at day 15.

Background

Treatment stratification based on early response to therapy greatly improved survival in pediatric acute lymphoblastic leukemia (ALL). Treatment response is assessed as quantified minimal residual disease (MRD) in bone marrow (BM) at set time points. For the last one or two decades, MRD assessment has been part of the routine clinical practice in the developed world. Flow cytometry and/or qPCR based methods are used in various countries, with sensitivities of 10−3–10−4 and 10−4–10−6, respectively [1]. Currently, application of high throughput sequencing methods are being studied to increase MRD sensitivity and to improve patient stratification [2]. Compared to individual IgH and TCR qPCR testing, next generation sequencing would also save working hours of laboratory staff. Whichever method is used, children are put under general anesthesia (GA) and undergo BM aspiration for MRD follow up 2–5 times following their diagnosis. This is because if little residual leukemia is left, it is most reliably and abundantly found in the BM [1]. GA and marrow aspiration, together with preparations and related issues (e.g. consenting, starving for anesthesia, checks, theatre lists, observations post procedure) put quite a burden on children, their families and on hospital staff.

MicroRNAs (miRs) are small non-coding RNAs, ranging from 19 to 25 nucleotides in size [3]. Their role is in posttranscriptional regulation of mRNA expression. MiRs regulate about half of the human genes [4]. Binding of miRs to their target sequence, with perfect complementarity in seed region and base paring in central region, results in mRNA cleavage and degradation. In case of less complementarity, translation repression is their typical function [5]. There are over 2000 miRs known in human [6]. The miR profile is distinct in every body fluid. In the blood, for example, platelets have a repertoire of approximately 750 miRs [7, 8], with the let-7 family contributing in almost 50% [8]. MiRs can be tissue specific regarding their origin, or there are examples when certain miRs are expressed at specific stages of ontogenesis or specific stages of cell differentiation [9, 10].

MiR profiles in normal tissues and in different diseases are intensively studied despite of numerous pitfalls [11]. MiRs are actively secreted into the extracellular space via extracellular vesicles [12]. The secreted miRs partly reflect the intracellular miR signature of the source cells, but certain miRs are withheld or selectively excreted in higher amounts [13]. Compared to mRNAs, miRs are much more stable. Some of them persist for 5 days [14, 15]. MiRs were found to be protected from RNases by encapsulation into extracellular vesicles, by complexes formed with Argonaute proteins or Nucleophosmin-1 and by their uptake into lipoproteins [16,17,18]. As a very universal phenomenon, malignant cells secrete higher amounts of miRs than healthy cells or tissues [19]. This, together with the selective over-secretion of some miRs, such as miR-128-3p [20,21,22], their easy distribution among fluid spaces and their stability makes miRs good candidates for cancer biomarkers. Indeed, very promising such results were published in solid tumors, e.g. in breast cancer and non-small cell lung cancer [23]. Some similar biomarker studies were published in leukemia, too. These focused on the miR signatures in BM, peripheral blood (PB) or plasma taken at the time of diagnosis as biomarkers of leukemia diagnosis, leukemia subtype, with disease characteristics or predictors of prognosis [24,25,26].

Our aim was to evaluate circulating miRs as biomarkers for MRD monitoring in pediatric ALL. Compared to the current practice of MRD monitoring, this approach would have the great advantage of less invasive sampling and no need for GA. Our strategy was to (1) get serial patient plasma samples along the induction cycle prospectively, (2) discard platelets very early as the largest source of miRs in the blood, (3) choose miRs for analyses based on the previously described miR profile of leukemic BM (4) for our first analysis on an exploratory cohort with array, and (5) proceed to validation of the selected best candidate biomarkers with qPCR on an extended patient cohort.

Methods

Patients and samples

We studied 28 pediatric patients with de novo and 5 patients with 1st relapse of precursor-B cell ALL, age 1–18 years, diagnosed between 2016 and 2018 at the 2nd Department of Pediatrics, Semmelweis University or at Heim Pál Children’s Hospital, both in Budapest, Hungary. Informed consent was requested from legal guardians of the patients. The study was approved by the Ethics Committee of the Hungarian Medical Research Council (60106-1/2015/EKU and 6886/2019/EKU) and conducted according to the principles of the Declaration of Helsinki. De novo ALL patients from only three cytogenetic subgroups were included: those with hyperdiploidy (blast chromosome number ≥ 47), those with t(12;21) (p13;q22) alias ETV6/RUNX1 fusion and those with no known cytogenetic abnormality (normal findings on a set FISH panel targeting frequent ALL rearrangements, found to be diploid by flow cytometry DNA index, either normal karyotype or unsuccessful karyotyping). Karyotype of relapsed patients was not used as selection criteria or for grouping because of the limited number of available samples.

Peripheral blood (PB) (3.5 ml) of patients was collected at diagnosis and also at days 8, 15 and 33 of the induction cycle of the ALL IC-BFM 2009 trial protocol. Bone marrow (BM) aspirates (1 ml) were collected at diagnosis before any treatment and on days 15 and 33. Blood and bone marrow samples of relapsed patients were collected only at diagnosis. In addition, control PB samples were collected from non-leukemic subjects. Control patients had one of the following conditions: vitamin D deficiency, otitis media, impetigo, hypothyreosis, neurofibromatosis type 1, iron deficiency anemia, phimosis.

Preparation of platelet-free plasma (PFP) was carried out within 2 h of sampling and was based on the recommendation of the International Society of Thrombosis and Haemostasis [27]. Sodium citrate tubes were used for all sample types. Samples taken at diagnosis were collected before any steroid or cytostatic treatment would have been started for the patient. The use of PFP was chosen based on the fact that platelets have a large repertoire of miRs and are the major source of miRs in plasma [28]. Whole blood and bone marrow samples were centrifuged at 2500×g for 15 min at 16 °C. Supernatant was transferred into a new centrifuge tube and centrifuged again for 15 min at 16 °C 2500×g. Supernatant was aliquoted immediately and stored at − 80 °C until use.

Clinical data were collected during the first month of therapy: white blood count at diagnosis, absolute blast number at day 8 in peripheral blood, bone marrow flow cytometry (FC) MRD at day 15, absolute white blood cell count at day 33, absolute lymphocyte count at day 33 and FC MRD at day 33. FC MRD was measured in routine diagnostics as per the standards of ALL IC-BFM 2009 study.

Detailed description of patients and samples used is provided in Table 1, 2, and 3.

Table 1 Samples analyzed on TaqMan low density array and with advanced qPCR
Table 2 Characteristics of the population studied on custom TaqMan low density array
Table 3 Characteristics of the population studied by advanced qPCR method

RNA isolation, cDNA synthesis and miR profiling on array

With the aim to find MRD biomarkers, 46 miRs were selected for analysis. Only miRs found to be overexpressed in blasts in at least two different previously published studies were chosen (Table 4).

Table 4 The 46 microRNAs selected for custom TaqMan low density array

MiRNeasy Serum/Plasma Kit (Qiagen, Hilden, Germany) was used for RNA isolation from defrosted PFP samples. Synthesis of cDNA was carried out using TaqMan Advanced miRNA cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA) according to the protocol supplied by the manufacturer.

MiR profiling of 65 PFP samples from 20 ALL patients (15 de novo and 5 relapsed ALL) and 10 control samples was performed using Custom TaqMan Advanced Low Density miRNA Array (TLDA) Card (Thermo Fisher Scientific, Waltham, MA, USA). Measurements were carried out using QuantStudio 7 Flex Real-Time PCR System (Thermo Fisher Scientific, Waltham, MA, USA). GeNorm algorithm was used to calculate the gene expression stability measure for each potential reference miR, based on the average pairwise variation between all candidate reference miRs. As a validation, NormFinder algorithm was also used which takes intra- and inter-group variability into account. Based on these analyses, miR-484 was selected for normalizing, and the ΔCt value was calculated for each miR. Four samples measured in duplicates on TLDA Card and correlated very highly (Pearson’s r = 0.981 (CI 95% 0.974–0.986), p < 10−16).

Quantitative real-time PCR

Based on the results of TLDA cards, 4 candidate miRs were selected for further investigation on an extended population using Advanced qPCR methodology. Quantitative Real-time PCR TaqMan Advanced miRNA assays were used according to the manufacturer’s instruction (Thermo Fisher Scientific, Waltham, MA, USA). Thermo Fisher 7900HT Fast Real-Time PCR system was used for real-time quantitative PCR measurement (Thermo Fisher Scientific, Waltham, MA, USA). PFP samples originated from peripheral blood of 28 patients with ALL and 10 controls were measured, in total 110 samples (Table 3 for details). Detection of 5 miRs (miR-181a-5p, miR-81b-5p, miR-128-3p, miR-22-3p and control miR-484) were performed in duplicates (correlation between duplicates was calculated: Pearson’s r = 0.998, p < 2.2*10−16). Correlation between Advanced qPCR and TLDA measurements was: Pearson’s r = 0.873, p = 3.19*10−44.

Statistical analysis

All statistical analyses were performed using R statistical software (R Foundation for Statistical Computing, Vienna, Austria; version 3.5.1). To select suitable reference miRs for normalization, the geNorm and NormFinder algorithms of the NormqPCR [41] package was used. Statistical differential expression of miRs was determined by the Limma [42] package. For that, a linear model was fitted for each gene based on the subgroup of the patients or the different sampling times. The gender of the patients was used as a coefficient in both cases. Then, moderated t-statistics and log-odds of differential expression were calculated by empirical Bayes moderation of the standard errors towards a common value. The resulting p-values obtained were corrected for multiple testing using the Benjamini–Hochberg procedure. Genes were considered to be differentially expressed when the adjusted p-value false discovery rate was below 0.05. Pearson correlation coefficients were calculated to evaluate the correlation between the flow cytometry MRD results and the log2 transformed expression level or fold change of miRs, respectively. In order to deal with multiple comparisons, the p-values were adjusted by the Benjamini–Hochberg false discovery rate method with type I error rate of 5%.

Results

MiR expression pattern of the exploratory cohort on array

Expressions of 46 selected miRs were measured in PFP samples of de novo and relapsed pediatric precursor-B acute lymphoblastic leukemia patients and controls on custom TaqMan low density array (TLDA). See the results illustrated in Figs. 1, 2, and 3. There was no significant correlation between miR expressions of the bone marrow and the peripheral blood (adjusted p > 0.05, data not shown).

Fig. 1
figure 1

TaqMan low density array miR expression, comparison of ALL subgroups and controls. The heatmap displays the expression pattern of microRNAs measured in platelet-free plasma derived from bone marrow and peripheral blood samples of pediatric ALL patients. The horizontal axis shows the origin of samples in different subgroups (normal karyotype: no alteration on FISH, DNA index diploid and cytogenetics normal or unsuccessful; HD: hyperdiploid subgroup), the day of sampling in relation to start of treatment. The vertical axis shows the 44 previously selected microRNA repertoire. miR-484 was used as normalizing miR. Expression of further 2 miRs was not detected in any samples. The scale shows the normalized log2 transformed expression levels. Yellow indicates microRNA down-expression, red indicates overexpression compared to the normalizing miR

Fig. 2
figure 2

MiR expression pattern in ALL and control PB PFP samples measured by TLDA. Case indicates de novo and relapsed ALL day 0 PB PFP samples, control indicates non-leukemic PB PFP samples. Data obtained on TLDA measurements. Box: the 2nd and 3rd quartiles; thick line in the box: median; lower whiskers: minimal value if there are no low outrange values, or Q1 − 1.5*IQR; upper whiskers: maximum value if there is no upper outlier, or Q3 + 1.5*IQR; dots: outliners, lines above the boxes: significant correlation (p < 0.05); ***p value: 0–0.001; **p value: 0.001–0.01; *p value: 0.01–0.05

Fig. 3
figure 3

MiR expression in de novo and relapsed ALL PB PFP at the day of diagnosis (by TLDA). Box: the 2nd and 3rd quartiles; thick line in the box: median; lower whiskers: minimal value if there are no low outrange values, or Q1 − 1.5*IQR; upper whiskers: maximum value if there is no upper outlier, or Q3 + 1.5. No significant differences were found

The expressions of 19 miRs in ALL PB PFP at diagnosis were significantly different compared to the control samples (p < 0.05). In line with our selection criteria, 18 miRs showed higher expression, while one miR showed downregulation in leukemic patients (Fig. 2). Different karyotype (hyperdiploid, t(12, 21), normal) subgroups were separately compared to the control PFP samples. Five miRs were significantly upregulated in every subgroup: miR-128-3p, miR-146a-5p, miR-181b-5p, miR-222-3p and miR-532-5p (p < 0.05). There was no unique miR expression pattern in subgroups of patients based on our TLDA measurements. Platelet originated miRs (let-7a-5p and let-7f-5p) were not significantly upregulated or downregulated compared to the control samples.

When miR expression differences were analysed between day 0 and day 33 of therapy in PB PFP, 4 miRs showed significantly lower expression at day 33: miR-128-3p, miR-181a-5p, miR-181b-5p, miR-222-3p. These 4 miRs were further investigated by Advanced qPCR on the extended cohort.

MiR expression pattern of the extended cohort by advanced qPCR

Relative expression levels of the 4 selected miRs were first compared between ALL day 0 PB PFP samples and PB PFP of non-leukemic controls. ALL subgroups of different karyotypes were examined separately as well as the whole ALL group (Table 5).

Table 5 miR expressions in diagnostic ALL samples compared to controls (PB PFP samples, advanced qPCR results)

Changes of PB PFP miR expression along the induction cycle were analyzed in paired tests. Higher expression of all 4 miRs identified on TLDA (miR-128-3p, miR-181a-5p, miR-181b-5p and miR-222-3p), was validated in the extended study population.

Measured by qPCR in this larger cohort, the expression of all four miRs decreased significantly after the start of treatment. However, no further significant drop could be demonstrated after day 8 with regards to any of the miRs (Fig. 4; Table 6). This question was also explored in the three ALL cytogenetic subgroups separately, and similar trends were found (Fig. 5). The best correlations and tendencies were seen with miR-181b-5p, miR-128-3p expression within the normal karyotype group, but no significant change was found after day 8 in any karyotype subgroup or any miR.

Fig. 4
figure 4

MiR expressions in PB PFP during the first month of therapy (by TaqMan Advanced qPCR). Box: the 2nd and 3rd quartiles; thick line in the box: median; whiskers: minimum and maximum values if there are no outrange values, or Q1 − 1.5*IQR; dots: outliers, lines above the boxes: significant correlation (p < 0.05); ***p value: 0 ≤ p < 0.001, **p value: 0.001 ≤ p < 0.01, *p value: 0.01 ≤ p < 0.05

Table 6 miR expression changes in peripheral blood PFP during the 1st month of chemotherapy
Fig. 5
figure 5

MiR expression during the first month of the therapy in cytogenetic subgroups of ALL (by qPCR). a Normal karyotype (no alteration on FISH, DNA index diploid and cytogenetics normal or unsuccessful); b hyperdiploid subgroup; c t(12;21). Box: the 2nd and 3rd quartiles; thick line: median; lower whiskers: minimal value if there are no low outlier values, or Q1 − 1.5*IQR; upper whiskers: maximum value if there is no upper outlier, or Q3 + 1.5*IQR; dots: outliners. Lines above the boxes: significant correlation (p < 0.05); ***p value: 0–0.001, **p value: 0.001–0.01, *p value: 0.01–0.05

Association of miR qPCR results with other MRD methods and prognostic factors

In order to further examine the value of the 4 selected miRs as MRD biomarkers, we tested their associations with results of other MRD measurement methods and with various risk factors: induction day 15 bone marrow flow cytometry MRD (the gold standard method to measure early treatment response in ALL-IC studies), day 8 absolute blast count in peripheral blood film by microscopic cytomorphology (also recorded in the ALL-IC and BFM studies), risk group stratification (based on age, initial white blood cell count (WBC), cytogenetics and early treatment response as per ALL-IC 2009) as well as end of induction absolute lymphocyte count (recently suggested as a prognostic marker) [43, 44]. All these questions were analyzed in the whole ALL group and in cytogenetic subgroups. See the significant associations highlighted in Table 7 and further results and details in Additional file 1.

Table 7 Associations between miR expression changes and MRD parameters or risk factors

All patients were MRD negative on day 33 bone marrow flow cytometry, hence this association could not be analyzed. To test if the blood WBC influences PB PFP miR expression causing bias, we checked this association on day 0 and day 33 and found no correlation.

The predictive accuracy of microRNAs to differentiate patients with high MRD (defined as induction day 15 bone marrow flow cytometry MRD greater than 1%) versus low MRD (< 1%) was also assessed by ROC analysis and by calculating the area under the receiver operating characteristic curve (AUC). MiR-128-3p had the highest predictive accuracy (AUC = 0.91) and miR-222-3p showed also high predictive accuracy (AUC = 0.79) (Fig. 6). See the full result set of ROC analyses in Additional file 2.

Fig. 6
figure 6

ROC curve analyses of miR expression changes between day 0 and 8. ROC curves showing the predictive accuracy of miR expression changes between day 0 and 8 to differentiate patients with high and low MRD. Area under the ROC curve (AUC) and its confidence interval are indicated. Patients with high MRD are defined as induction day 15 bone marrow flow cytometry MRD ≥ 1%. a miR-181b-5p; b miR-181a-5p; c miR-128-3p; d miR-222-3p

Discussion

We hypothesized that circulating miRs may be used as biomarkers of residual leukemia with the advantage of less invasive sampling for follow-up, compared to bone marrow aspirate based tests in present practice. We focused on pediatric precursor-B cell ALL patients.

Platelet-free plasma (PFP) derived from peripheral blood (PB) was chosen for analyses. In the first step, 46 candidate miRs were examined on a discovery population of de novo plus relapsed ALL patients and controls. Four miRs were selected for further analyses based on higher expression compared to controls and based on their decreasing expression by the end of induction: miR-128-3p, miR-181a-5p, miR-181b-5p and miR-222-3p. These four miRs were then tested by advanced qPCR at 4 time points along the induction cycle of an extended cohort of de novo ALL patients. Over the first week, the expression of all four miRs dropped significantly but there was no relevant further decrease at later time points. Correlation of miR expression with standard MRD results and risk factors could be demonstrated. On ROC analyses, the drop of normalized miR expression (day 8 or day 15 compared to day 0) performed better than the normalized expressions themselves at any given time points. Day 8 miR expressions performed somewhat better than day 15 miR expressions. After all, the day 0 to day 8 change in normalized expression of miR-128-3p and miR-181b-5p came out as best MRD-biomarkers.

To summarize our results, the selected circulating miRs can indeed be used as biomarkers of residual leukemia burden but are less sensitive than the day 15 bone marrow flow cytometry MRD measurement, the standard response measurement in ALL-IC. We know it from flow cytometry (as well as IgH and TCR based qPCR) MRD measurements that in reality, a monotonous multiple-log reduction occurs in the bone marrow blast percentage from day 0, through day 8 and day 15 to day 33 along the induction cycle [45,46,47]. By blood PFP miR measurement in our study, only a day 0 to day 8 reduction could be demonstrated. The most plausible explanation is that the overexpression of miRs originated from residual leukemic cells is measurable in the blood plasma but the noise (baseline miR expressions and their variability) is too high to allow for a sensitive leukemia follow-up. Indeed, the expression of the same miRs could be detected in the serum of healthy adults [48]. It is difficult to tell what this background noise is related to. The general understanding about miRs is that these participate in many regulation pathways in normal physiology. In microRNA Target Prediction Database (miRDB.org), an online database miR-128-3p has 1254 predicted targets while miR-222-3p has 619, miR-181a-5p and miR-181b-5p have 1408 predicted targets [49].

Our results are in parallel with findings of multiple studies which examined miR expression in leukemic cells [3, 34, 35, 50,51,52,53,54,55]. In these studies, miRs were extracted from whole bone marrow or from whole blood of pediatric- or partly pediatric-ALL patients, from mononuclear cells separated of the above samples, or from ALL cell lines. Very different sets of miRs were investigated, all partially overlapping with our selection. The 4 miRs chosen in our study were unanimously described as being overexpressed in ALL blasts. The reported rate of overexpression compared to healthy controls or to ALL patients in remission spread between 1.3 and 500 times, with very diverse normalizing RNAs used. The findings supported by most evidence are the overexpression of miR-128 and miR-181-family members. Many other overexpressed miRs were reported in ALL, e.g. let-7b, miR-222, miR-181b [31, 34].

Among the cellular miR studies, some reported different miR expression patterns in ALL cytogenetic subgroups [31, 35, 36, 50], which could not be demonstrated in our present study in PFP. Interestingly, authors disagree if the miR expression pattern of bone marrow and peripheral blood are similar [50, 51] or different [52] in intrapersonal comparisons. We did not find strong correlation between marrow PFP and blood PFP miR expressions in our TLDA results (adjusted p > 0.05, data not shown).

Three publications are available that evaluated circulating miRs in pediatric ALL. These studied plasma or serum unlike our present study (PFP samples). The identified miRs are different from ours and from each other, however, all three papers conclude that circulating miRs offer a very sensitive and specific diagnostic tool to differentiate pediatric precursor-B cell ALL cases from healthy controls [31, 56, 57]. We are not aware of any previous study which would have evaluated PFP miR expression in pediatric ALL, or would have evaluated circulating miRs as ALL MRD biomarkers. There are encouraging such publications on patient cohorts of other malignancies, though [58, 59].

Our study has special strengths when compared to other papers. A possible bias by platelets—as massive sources of circulating miRs—could be ruled out because no association of leukemia MRD and platelet originated miRs was found despite the wide range of platelet count of our patients. As no correlation could be detected between PFP miR expression levels and WBC either on day 0 (all patients with full blown ALL, before treatment) or day 33 (all the patients were MRD negative with bone marrow flow cytometry), we are convinced that the measured miR expressions were not primarily influenced by healthy leukocytes or blasts in the peripheral blood.

Mir-181 might contribute to leukemogenesis and relapse [35], it is the inducer of Wnt/β-catenin signaling pathway. The uncontrolled pathway causes leukemia by stimulating cell proliferation. It has been demonstrated using ALL-derived cell lines and bone marrow samples of patients with ALL that miR-181a-5p targets and inactivates Wnt antagonist WIF1 (Wnt inhibitory factor-1) [60].

Regarding the biology in the background, miR-128 was found to inhibit or promote tumor proliferation in several tumors and its proved target genes alter cancer-related biological processes [61]. Its expression was higher in ALL cell lines [51], in patients with poor-prognosis ALL [51], in bone marrow samples of children with ALL at diagnosis, at relapse, but not in remission [35]. Expression of miR-128 decreased during therapy [51]. One of the predicted putative targets of miR-128b is BMI1, a transcriptional factor which is important in hematopoietic stem cells and leukemia stem-cell self-renewal [35, 61]. In silico predictions suggest that miR-128 and miR-181 are among the major regulators of lymphoid differentiation. These miRs may uphold early progenitor cells at early stem-progenitor stage, prevent their further differentiation towards more mature cells by regulating molecules having crucial role in early steps of hematopoiesis. MiR-222, most likely, blocks hematopoietic differentiation at terminal stages [62]. Chromosomal location of miR-222 is close to Xp11.3 and its target is C-KIT’s 3′UTR. C-kit is a tyrosine kinase receptor with a role in cell differentiation and growth. It is expressed in hematopoietic stem cells and progenitor cells [63].

Lacking standards regarding normalizing miRs and only partial correlation of various PCR based methods which are widely known difficulties of RNA and miR lab methodologies contribute to the limitations of this paper. Still, we were able to prove our points via both methods applied. Our sample collection started in 2016, so it was not possible to correlate miR expressions to patient survival yet, and the cohort would also be small for that analysis. Given the candidate miR approach used, the existence of a better miR biomarker of leukemia MRD cannot be ruled out.

Conclusion

In summary, circulating miR-128-3p and miR-181b-5p expression behave as biomarkers of residual leukemia in pediatric ALL. However, these are less sensitive than current standards of patient follow-up.

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

ABC:

absolute blast count

ALC:

absolute lymphocyte count

ALL:

acute lymphoblastic leukemia

BM:

bone marrow

FC:

flow cytometry

FISH:

fluorescent in situ hybridisation

GA:

general anaesthesia

HD:

hyperdiploid

IgH:

immunoglobulin heavy chain

miR:

microRNA

MRD:

minimal residual disease

mRNA:

messenger RNA

PB:

peripheral blood

PCR:

polymerase chain reaction

PFP:

platelet-free plasma

RNA:

ribonucleic acid

RNase:

ribonuclease

TCR:

T cell receptor

TLDA:

TaqMan low density array

WBC:

white blood cell

WIF1:

Wnt inhibitor factor-1

References

  1. van Dongen JJM, van der Velden VHJ, Brüggemann M, Orfao A. Minimal residual disease diagnostics in acute lymphoblastic leukemia: need for sensitive, fast, and standardized technologies. Blood. 2015;125(26):3996–4009.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Wood B, Wu D, Crossley B, Dai Y, Williamson D, Gawad C, et al. Measurable residual disease detection by high-throughput sequencing improves risk stratification for pediatric B-ALL. Blood. 2018;131(12):1350–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Zhang H, Luo X-Q, Zhang P, Huang L-B, Zheng Y-S, Wu J, et al. MicroRNA patterns associated with clinical prognostic parameters and CNS relapse prediction in pediatric acute leukemia. PLoS ONE. 2009;4(11):e7826.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Umerez M, Garcia-Obregon S, Martin-Guerrero I, Astigarraga I, Gutierrez-Camino A, Garcia-Orad A. Role of miRNAs in treatment response and toxicity of childhood acute lymphoblastic leukemia. Pharmacogenomics. 2018;19:361–73.

    Article  CAS  PubMed  Google Scholar 

  5. Mohr AM, Mott JL. Overview of microRNA biology. Semin Liver Dis. 2015;35(1):3–11.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Hammond SM. An overview of microRNAs. Adv Drug Deliv Rev. 2015;87:3–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Bray PF, McKenzie SE, Edelstein LC, Nagalla S, Delgrosso K, Ertel A, et al. The complex transcriptional landscape of the anucleate human platelet. BMC Genomics. 2013;14:1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Ple H, Landry P, Benham A, Coarfa C, Gunaratne PH, Provost P. The repertoire and features of human platelet microRNAs. PLoS ONE. 2012;7(12):e50746.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, et al. A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007;129(7):1401–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, et al. The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature. 2000;403(6772):901–6.

    Article  CAS  PubMed  Google Scholar 

  11. Witwer KW. Circulating MicroRNA biomarker studies: pitfalls and potential solutions. Clin Chem. 2015;61(1):56–63.

    Article  CAS  PubMed  Google Scholar 

  12. Thery C. Exosomes: secreted vesicles and intercellular communications. F1000 Biol Rep. 2011;3:15.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Yanez-Mo M, Siljander PR, Andreu Z, Zavec AB, Borras FE, Buzas EI, et al. Biological properties of extracellular vesicles and their physiological functions. J Extracell Vesicles. 2015;4:27066.

    Article  PubMed  Google Scholar 

  14. Gantier MP, McCoy CE, Rusinova I, Saulep D, Wang D, Xu D, et al. Analysis of microRNA turnover in mammalian cells following Dicer1 ablation. Nucleic Acids Res. 2011;39(13):5692–703.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Leidinger P, Backes C, Meder B, Meese E, Keller A. The human miRNA repertoire of different blood compounds. BMC Genomics. 2014;15:474.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Wang K, Zhang S, Weber J, Baxter D, Galas DJ. Export of microRNAs and microRNA-protective protein by mammalian cells. Nucleic Acids Res. 2010;38(20):7248–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT. MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. 2011;13(4):423–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Mittelbrunn M, Gutierrez-Vazquez C, Villarroya-Beltri C, Gonzalez S, Sanchez-Cabo F, Gonzalez MA, et al. Unidirectional transfer of microRNA-loaded exosomes from T cells to antigen-presenting cells. Nat Commun. 2011;2:282.

    Article  PubMed  CAS  Google Scholar 

  19. Schwarzenbach H, Nishida N, Calin GA, Pantel K. Clinical relevance of circulating cell-free microRNAs in cancer. Nat Rev Clin Oncol. 2014;11:145.

    Article  CAS  PubMed  Google Scholar 

  20. Mets E, Van Peer G, Van der Meulen J, Boice M, Taghon T, Goossens S, et al. MicroRNA-128-3p is a novel oncomiR targeting PHF6 in T-cell acute lymphoblastic leukemia. Haematologica. 2014;99(8):1326–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Zhu YD, Wang L, Sun C, Fan L, Zhu DX, Fang C, et al. Distinctive microRNA signature is associated with the diagnosis and prognosis of acute leukemia. Med Oncol (Northwood, London, England). 2012;29(4):2323–31.

    Article  CAS  Google Scholar 

  22. Fulci V, Colombo T, Chiaretti S, Messina M, Citarella F, Tavolaro S, et al. Characterization of B- and T-lineage acute lymphoblastic leukemia by integrated analysis of MicroRNA and mRNA expression profiles. Genes Chromosom Cancer. 2009;48(12):1069–82.

    Article  CAS  PubMed  Google Scholar 

  23. Wang H, Peng R, Wang J, Qin Z, Xue L. Circulating microRNAs as potential cancer biomarkers: the advantage and disadvantage. Clin Epigenet. 2018;10(1):59.

    Article  CAS  Google Scholar 

  24. Yingchun L, Rong Z, Kun Y, Ying Y, Zhuogang L. Bone marrow MicroRNA-335 level predicts the chemotherapy response and prognosis of adult acute myeloid leukemia. Medicine. 2015;94(33):e0986.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Fathullahzadeh S, Mirzaei H, Honardoost MA, Sahebkar A, Salehi M. Circulating microRNA-192 as a diagnostic biomarker in human chronic lymphocytic leukemia. Cancer Gene Ther. 2016;23(10):327–32.

    Article  CAS  PubMed  Google Scholar 

  26. Papageorgiou SG, Kontos CK, Diamantopoulos MA, Bouchla A, Glezou E, Bazani E, et al. MicroRNA-155-5p overexpression in peripheral blood mononuclear cells of chronic lymphocytic leukemia patients is a novel, independent molecular biomarker of poor prognosis. Dis Markers. 2017;2017:2046545.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Lacroix R, Judicone C, Mooberry M, Boucekine M, Key NS, Dignat-George F, et al. Standardization of pre-analytical variables in plasma microparticle determination: results of the International Society on Thrombosis and Haemostasis SSC Collaborative workshop. J Thromb Haemost. 2013;11(6):1190–3.

    Article  Google Scholar 

  28. Mitchell AJ, Gray WD, Hayek SS, Ko Y-A, Thomas S, Rooney K, et al. Platelets confound the measurement of extracellular miRNA in archived plasma. Sci Rep. 2016;6:32651.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Marabita F, de Candia P, Torri A, Tegnér J, Abrignani S, Rossi RL. Normalization of circulating microRNA expression data obtained by quantitative real-time RT-PCR. Brief Bioinform. 2016;17(2):204–12.

    Article  PubMed  Google Scholar 

  30. Schwarzenbach H, da Silva AM, Calin G, Pantel K. Data normalization strategies for MicroRNA quantification. Clin Chem. 2015;61(11):1333–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Luna-Aguirre CM, de la Luz Martinez-Fierro M, Mar-Aguilar F, Garza-Veloz I, Trevino-Alvarado V, Rojas-Martinez A, et al. Circulating microRNA expression profile in B-cell acute lymphoblastic leukemia. Cancer Biomarkers. 2015;15(3):299–310.

    Article  CAS  PubMed  Google Scholar 

  32. Vilella F, Moreno-Moya JM, Balaguer N, Grasso A, Herrero M, Martinez S, et al. Hsa-miR-30d, secreted by the human endometrium, is taken up by the pre-implantation embryo and might modify its transcriptome. Development (Cambridge, England). 2015;142(18):3210–21.

    Article  CAS  Google Scholar 

  33. Rinnerthaler G, Hackl H, Gampenrieder SP, Hamacher F, Hufnagl C, Hauser-Kronberger C, et al. miR-16-5p is a stably-expressed housekeeping MicroRNA in breast cancer tissues from primary tumors and from metastatic sites. Int J Mol Sci. 2016;17(2):e156.

    Article  PubMed  CAS  Google Scholar 

  34. de Oliveira JC, Scrideli CA, Brassesco MS, Morales AG, Pezuk JA, Queiroz Rde P, et al. Differential miRNA expression in childhood acute lymphoblastic leukemia and association with clinical and biological features. Leuk Res. 2012;36(3):293–8.

    Article  PubMed  CAS  Google Scholar 

  35. Han B-W, Feng D-D, Li Z-G, Luo X-Q, Zhang H, Li X-J, et al. A set of miRNAs that involve in the pathways of drug resistance and leukemic stem-cell differentiation is associated with the risk of relapse and glucocorticoid response in childhood ALL. Hum Mol Genet. 2011;20(24):4903–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Schotte D, De Menezes RX, Akbari Moqadam F, Khankahdani LM, Lange-Turenhout E, Chen C, et al. MicroRNA characterize genetic diversity and drug resistance in pediatric acute lymphoblastic leukemia. Haematologica. 2011;96(5):703–11.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Harada M, Pokrovskaja-Tamm K, Soderhall S, Heyman M, Grander D, Corcoran M. Involvement of miR17 pathway in glucocorticoid-induced cell death in pediatric acute lymphoblastic leukemia. Leukemia Lymphoma. 2012;53(10):2041–50.

    Article  CAS  PubMed  Google Scholar 

  38. Diakos C, Zhong S, Xiao Y, Zhou M, Vasconcelos GM, Krapf G, et al. TEL-AML1 regulation of survivin and apoptosis via miRNA-494 and miRNA-320a. Blood. 2010;116(23):4885–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Neilson JR, Zheng GX, Burge CB, Sharp PA. Dynamic regulation of miRNA expression in ordered stages of cellular development. Genes Dev. 2007;21(5):578–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Avigad S, Verly IR, Lebel A, Kordi O, Shichrur K, Ohali A, et al. miR expression profiling at diagnosis predicts relapse in pediatric precursor B-cell acute lymphoblastic leukemia. Genes Chromosom Cancer. 2016;55(4):328–39.

    Article  CAS  PubMed  Google Scholar 

  41. Perkins JR, Dawes JM, McMahon SB, Bennett DL, Orengo C, Kohl M. ReadqPCR and NormqPCR: r packages for the reading, quality checking and normalisation of RT-qPCR quantification cycle (Cq) data. BMC Genomics. 2012;13:296.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Sun D, Elson P, Liedtke M, Medeiros BC, Earl M, Alizadeh A, et al. Absolute lymphocyte count at day 28 independently predicts event-free and overall survival in adults with newly diagnosed acute lymphoblastic leukemia. Am J Hematol. 2012;87(10):957–60.

    Article  PubMed  Google Scholar 

  44. Farkas T, Muller J, Erdelyi DJ, Csoka M, Kovacs GT. Absolute Lymphocyte Count (ALC) after induction treatment predicts survival of pediatric patients with acute lymphoblastic leukemia. Pathol Oncol Res. 2017;23(4):889–97.

    Article  CAS  PubMed  Google Scholar 

  45. Schumich A, Maurer-Granofszky M, Attarbaschi A, Potschger U, Buldini B, Gaipa G, et al. Flow-cytometric minimal residual disease monitoring in blood predicts relapse risk in pediatric B-cell precursor acute lymphoblastic leukemia in trial AIEOP-BFM-ALL 2000. Pediatr Blood Cancer. 2019;66(5):e27590.

    Article  CAS  PubMed  Google Scholar 

  46. Conter V, Bartram CR, Valsecchi MG, Schrauder A, Panzer-Grumayer R, Moricke A, et al. Molecular response to treatment redefines all prognostic factors in children and adolescents with B-cell precursor acute lymphoblastic leukemia: results in 3184 patients of the AIEOP-BFM ALL 2000 study. Blood. 2010;115(16):3206–14.

    Article  CAS  PubMed  Google Scholar 

  47. Basso G, Veltroni M, Valsecchi MG, Dworzak MN, Ratei R, Silvestri D, et al. Risk of relapse of childhood acute lymphoblastic leukemia is predicted by flow cytometric measurement of residual disease on day 15 bone marrow. J Clin Oncol. 2009;27(31):5168–74.

    Article  PubMed  Google Scholar 

  48. Juzenas S, Venkatesh G, Hubenthal M, Hoeppner MP, Du ZG, Paulsen M, et al. A comprehensive, cell specific microRNA catalogue of human peripheral blood. Nucleic Acids Res. 2017;45(16):9290–301.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Wong N, Wang X. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 2015;43(Database issue):D146–52.

    Article  CAS  PubMed  Google Scholar 

  50. Mi S, Lu J, Sun M, Li Z, Zhang H, Neilly MB, et al. MicroRNA expression signatures accurately discriminate acute lymphoblastic leukemia from acute myeloid leukemia. Proc Natl Acad Sci USA. 2007;104(50):19971–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Nemes K, Csóka M, Nagy N, Márk Á, Váradi Z, Dankó T, et al. Expression of certain leukemia/lymphoma related microRNAs and its correlation with prognosis in childhood acute lymphoblastic leukemia. Pathol Oncol Res. 2015;21(3):597–604.

    Article  CAS  PubMed  Google Scholar 

  52. Duyu M, Durmaz B, Gunduz C, Vergin C, Yilmaz Karapinar D, Aksoylar S, et al. Prospective evaluation of whole genome microRNA expression profiling in childhood acute lymphoblastic leukemia. BioMed Res Int. 2014;2014:967585.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Schotte D, Chau JC, Sylvester G, Liu G, Chen C, van der Velden VH, et al. Identification of new microRNA genes and aberrant microRNA profiles in childhood acute lymphoblastic leukemia. Leukemia. 2009;23(2):313–22.

    Article  CAS  PubMed  Google Scholar 

  54. Shafik RE, Abd El Wahab N, Senoun SA, Ebeid E, El Taweel MA. Expression of micro-RNA 128 and let-7b in Pediatric acute lymphoblastic leukemia cases. Asian Pac J Cancer Prev. 2018;19(8):2263–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Zhou J, Wang S, Sun K, Chng WJ. The emerging roles of exosomes in leukemogeneis. Oncotarget. 2016;7(31):50698–707.

    PubMed  PubMed Central  Google Scholar 

  56. Swellam M, Hashim M, Mahmoud MS, Ramadan A, Hassan NM. Aberrant expression of some circulating miRNAs in childhood acute lymphoblastic leukemia. Biochem Genet. 2018;56(4):283–94.

    Article  CAS  PubMed  Google Scholar 

  57. Jiang Q, Lu X, Huang P, Gao C, Zhao X, Xing T, et al. Expression of miR-652-3p and effect on apoptosis and drug sensitivity in pediatric acute lymphoblastic leukemia. Biomed Res Int. 2018;2018:5724686.

    PubMed  PubMed Central  Google Scholar 

  58. McDonald AC, Vira M, Shen J, Sanda M, Raman JD, Liao J, et al. Circulating microRNAs in plasma as potential biomarkers for the early detection of prostate cancer. Prostate. 2018;78(6):411–8.

    Article  CAS  PubMed  Google Scholar 

  59. Ostenfeld MS, Jensen SG, Jeppesen DK, Christensen LL, Thorsen SB, Stenvang J, et al. miRNA profiling of circulating EpCAM(+) extracellular vesicles: promising biomarkers of colorectal cancer. J Extracell Vesicles. 2016;5:31488.

    Article  PubMed  CAS  Google Scholar 

  60. Lyu X, Li J, Yun X, Huang R, Deng X, Wang Y, et al. miR-181a-5p, an inducer of Wnt-signaling, facilitates cell proliferation in acute lymphoblastic leukemia. Oncol Rep. 2017;37(3):1469–76.

    Article  CAS  PubMed  Google Scholar 

  61. Li M, Fu W, Wo L, Shu X, Liu F, Li C. miR-128 and its target genes in tumorigenesis and metastasis. Exp Cell Res. 2013;319(20):3059–64.

    Article  CAS  PubMed  Google Scholar 

  62. Georgantas RW 3rd, Hildreth R, Morisot S, Alder J, Liu CG, Heimfeld S, et al. CD34+ hematopoietic stem-progenitor cell microRNA expression and function: a circuit diagram of differentiation control. Proc Natl Acad Sci USA. 2007;104(8):2750–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Ju X, Li D, Shi Q, Hou H, Sun N, Shen B. Differential microRNA expression in childhood B-cell precursor acute lymphoblastic leukemia. Pediatr Hematol Oncol. 2009;26(1):1–10.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

The authors are thankful to all the patients and control subjects, nurses, physicians who took part in this study. We also thank Mónika Sándorné Vángor for the technical assistance.

Funding

This research was supported by the National Research, Development and Innovation Office (NKFIH, Grant No. K115861, D.J. Erdélyi, and Grant No. K112872, C. Szalai), the Hungarian Pediatric Oncology Network (Grants 09/MGYH-MGYGYT/2018 and 01/MGYH-MGYGYT/2018), the BME-Biotechnology FIKP Grant of EMMI (BME FIKP-BIO) and the Hungarian Society of Haematology and Transfusion (Grant to D.J.E. in 2017). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

AR performed the laboratory work, collected clinical data and wrote the manuscript. NK performed the laboratory work. AG performed statistical analysis. JCS performed the laboratory work. BE performed bioinformatics. GP provided patient samples and clinical data. HB performed bioinformatics. GN performed laboratory work. JM provided patient samples and clinical data. GTK supervised the research and provided patient samples. CS supervised the research and contributed to experimental discussion. AFS supervised the research, contributed to experimental discussion and reviewed the manuscript. DJE designed the study, analyzed the data and wrote the paper.

Corresponding author

Correspondence to Dániel J. Erdélyi.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee of the Hungarian Medical Research Council (60106-1/2015/EKU and 6886/2019/EKU) and conducted according to the principles of the Declaration of Helsinki.

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.

Supplementary information

Additional file 1.

Associations between miR expression changes and MRD parameters or risk factors. In these analyses, the dependent variable is the normalized miR expression change between the time points compared, all measured in PB PFP. Fold change indicates change in the normalized miR expression between the two time points, the median of these values are provided. Pearson’s r indicates the correlation between dependent and independent variables. Adjusted p refers to the statistical significance of the Pearson’s correlation. *Normal karyotype: no alteration on FISH, DNA index diploid and cytogenetics normal or unsuccessful.

Additional file 2.

ROC curve analysis. The predictive accuracy of miRs to differentiate patients with high MRD (defined as induction day 15 bone marrow flow cytometry MRD greater than 1%) versus low MRD (< 1%) was assessed by ROC analysis and by calculating the area under the receiver operating characteristic curve (AUC).

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rzepiel, A., Kutszegi, N., Gézsi, A. et al. Circulating microRNAs as minimal residual disease biomarkers in childhood acute lymphoblastic leukemia. J Transl Med 17, 372 (2019). https://doi.org/10.1186/s12967-019-2114-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-019-2114-x

Keywords