Introduction
In order to be benchmarked, a novel biomarker should be: i) measurable in a specific clinical setting, ii) validated by multiple studies, iii) directly impacting on the medical decision (1). Although they do not yet satisfy all these criteria, microRNAs (miRNAs), and circulating miRNAs in particular, have gained a place among new-generation biomarkers. Circulating miRNAs can be measured in 12 types of body fluids (e.g., plasma and serum where they are particularly abundant, blood, urine, cerebrospinal fluid, saliva) and 61 different types of tissues (2, 3). Nonetheless, their validation and impact on medical decision criteria are still lacking. Validation is hindered by the lack of analytical (i.e., low inter-platform concordance) and postanalytical (i.e., normalization strategy) standardization (4, 5). Their impact on medical decision is weakened by uncertainty about the meaning of the changes in their expression: being fine-tuners of gene expression, miRNA expression is highly variable (6).
Despite these limitations, the abnormal expression of a single miRNA or miRNA families has been more or less specifically associated with disorders such as cancer, cardiovascular diseases, inflammatory and autoimmune diseases (7). Indeed, the diagnostic potential of the whole miRNAs expression profile (“miRNome”) could be even greater. Hence, given their key biological role, it is desirable to increase their implementation in routine diagnostics. To this end, standardization of preanalytical and postanalytical phases are needed (6, 8). Furthermore, current knowledge is still poor about the preanalytical warnings concerning samples intended for miRNAs analysis.
Circulating miRNAs can be passively released from cells because of physiological renewal or cell death, constituting free-circulating miRNAs associated with carrier proteins such as argonaute (Ago)-2 and nucleophosmin (NPM)-1, or they can be actively secreted encapsulated within microvesicles/exosomes targeting other cells/tissues (9). Their concentration in plasma, as estimated by quantitative real-time polymerase chain reaction (qPCR), is physiologically highly variable and ranges between 1 and 308 μg/L (10, 11). Therefore, any variable other than the intrinsic biological variables would impact on the outcome of miRNAs expression (12).
Among the still unsolved questions about the preanalytical management of samples for miRNA analysis, the choice of sample matrix poses multiple challenges. Although this has never been considered as a real issue in practical terms, it does have implications for the final analytical output (6). Previous studies have demonstrated the stability of miRNAs in serum and plasma, with a higher variability for the former than the latter, as also recently confirmed by Basso et al. (13-16). For plasma, the choice of the anticoagulant is pivotal. Because heparin and citrate interfere with the enzyme activities in PCR, they are not recommended for qPCR-based miRNAs quantitation (2). However, heparinase treatment improves the results in heparin-anticoagulated plasma, although Basso et al. reported a severe interference by heparin (15, 16). In addition, ethylendiaminetetraacetate (EDTA) potassium salts (K2 and K3), which can be removed from the PCR mastermix, are regarded as the best anticoagulants for assaying circulating miRNAs (17). Sodium fluoride/potassium oxalate (NaF/KOx), an alternative to EDTA, increases miRNA detection (15).
Against this background, we investigated the preanalytical variables, including sample matrix, time- and temperature-dependent stability, in the measurement of a wide panel of circulating miRNAs in EDTA-anticoagulated plasma. In addition, we checked for possible contamination from platelet-derived miRNAs. To this end we assayed: i) two different EDTA-anticoagulated sample matrices either with or without gel separator (plasma preparation tubes (PPT) and K2EDTA tubes, respectively); ii) two different centrifugation settings for plasma separation (single standard and single standard plus platelet-depleting centrifugation); iii) three different plasma storage conditions (immediate freezing at - 80°C, and storage for 24h at either room temperature or 4°C).
Materials and methods
Study design and subjects
This observational study was designed to determine the stability and quantification of circulating miRNAs in plasma subjected to different preanalytical variables, such as sample collection and storage conditions.
In this study, ten adult male volunteers (aged 25 – 40 years) were recruited among the personnel of the IRCCS Istituto Ortopedico Galeazzi, between November 2016 and January 2017. Eligibility criteria were: global health status with no acute or chronic disease and no chronic or current ongoing medical treatment, normal body-mass index (20 kg/m2 ≤ BMI < 25 kg/m2), and non-smoking. Subjects were asked to refrain from alcohol intake, physical activity, pro-inflammatory activities (e.g., shaving and waxing) in the 48h preceding blood drawing, and to keep a regular diet.
Venous blood was collected in K2EDTA and PPT (BD Vacutainer, Becton Dickinson, Milan, Italy). Samples were centrifuged at standard condition of 1300xg for 10 min at room temperature (RT) of 22°C. The plasma was collected and an aliquot from the K2EDTA tubes was further centrifuged at 2500xg for 15 min at RT to eliminate platelet residues and to obtain the platelet-poor plasma (PPP) (18). Plasma aliquots from K2EDTA, PPP, and PPT were frozen at - 80°C either immediately or after storage for 24 h at either RT or 4°C. The protocol was approved by the ethical committee (Comitato Etico Ospedale San Raffaele, Milano, on 4 February 2016) and registered with ClinicalTrials.gov (registration number NCT03386981, 4th January 2018). Written, informed consent to participate was obtained from all subjects after being informed about procedures, risks, and benefits associated with the research, as required by the Helsinki Declaration.
Methods
Total plasma RNA enriched in miRNAs was isolated with the miRCURY RNA Isolation Kit (Exiqon A/S, Vedbaek, Denmark) and reverse transcription was carried out using miRCURY LNA™ Universal reverse transcription miRNA PCR (Exiqon A/S, Vedbaek, Denmark). The spike-in controls UniSp2, UniSp4, UniSp5, and UniSp6 and cel-39 (Exiqon A/S, Vedbaek, Denmark), added to each sample at recommended concentration, were used to check the efficiency of RNA extraction and reverse transcription, respectively. Real-time qPCR (StepOne Plus, Applied Biosystems, Foster City, USA) was performed with ExiLENT SYBR Green Master Mix (Exiqon A/S, Vedbaek, Denmark) and a primer set of the 179 miRNAs most highly represented in the circulation (serum/plasma miRCURY LNA miRNA focus panel, Exiqon A/S, Vedbaek, Denmark). Three inter-plate calibrators (IPC), with specific primer pairs and DNA template, were provided with the panel kit. The PCR was conducted over 40 cycles (10 sec at 95°C, 60 sec at 60°C) preceded by 10-min activation of the polymerase at 95°C and followed the melting curves. Each sample was analysed in triplicate.
Reverse transcription-qPCR data analysis was conducted as previously described (8). Briefly, the threshold cycle (CT) values from qPCR were analysed with GenEx software v6 (Exiqon A/S, Vedbaek, Denmark). Threshold cycle values from the plate runs of different experiments were adjusted based on IPC CT values. Only those miRNAs with a CT < 37 were further considered in the analysis. The 2−ΔΔCT method was applied to calculate the relative expression of each miRNA and the global mean (i.e., the geometrical mean of the whole dataset CT values) was used to normalize the data, as this strategy is defined as being highly reliable in this specific experimental setting according to Faraldi et al. (8). The CT difference between hsa-miR-23a and hsa-miR-451 was used to check for micro-haemolysis, which can eventually affect the analysis (19). None of the samples exceeded the cut-off value of 7.
Comparisons of detected and undetected miRNAs in the different analysed conditions were performed to define miRNA detectability. Fold change (FC, the ratio between the normalized outputs from an “experimental” condition and “reference” condition) analyses were performed to determine the differences in circulating miRNA expression in all the different conditions. More specifically, we focused on the measurement of clinically relevant miRNAs: miRNAs associated with bone fracture risk; skeletal muscle modulated, both in tissue and circulation, during physical activity; and miRNAs associated with cancer diagnosis (with a reported area under curve (AUC) > 0.70) (20-23). Results are shown only for the miRNAs that were significantly ≥ 5-fold either up- or down-regulated, in at least one condition.
Statistical analysis
Detected and undetected miRNAs were compared within- and between-group considering the sample collection or the storage condition, using Pearson’s Chi-Square test with Yates’ continuity correction. Additionally, both Pearson’s Chi-Square test with Yates’ continuity correction and McNemar’s Chi-square test with continuity correction, for paired data, were performed comparing two conditions at a time. The non parametric Kruskal-Wallis test with Dunn’s multiple comparison test was used for within- and between-group comparison of the expression of each miRNA. Differences were considered statistically significant if P-values < 0.05. Statistical analysis was performed with Statistical 13.1 software (Statsoft, Tulsa, USA).
Results
Effect of sample matrix on samples immediately frozen at - 80°C
The number of undetected miRNAs in plasma collected in K2EDTA tubes was significantly higher compared to plasma in PPT tubes immediately frozen at - 80°C (as shown in Table 1).
Table 1
Because 15 miRNAs (hsa-miR-100-5p, hsa-miR-133a, hsa-miR136-3p, hsa-miR136-5p, hsa-miR-141-3p, hsa-miR-20b-5p, hsa-miR-200a-5p, hsa-miR-208-3p, hsa-miR22-5p, hsa-miR-33a-5p, hsa-miR-335-5p, hsa-miR-362, hsa-miR-365a, hsa-miR497-5p, and hsa-miR-7-1-3p) were undetectable in both conditions, all the following analyses were performed on 164 miRNAs. Two miRNAs (1% of the 164 miRNAs) were either up- or down-regulated in PPT compared to K2EDTA samples (Table 2).
Table 2
Effect of platelet depletion on samples immediately frozen at - 80°C
Platelet depletion was highlighted by the lower plasma expression of hsa-miR-126-3p, hsa-miR-223-3p, and hsa-miR-27b-3p, three of the most highly expressed miRNAs in platelets, in PPP samples compared to K2EDTA and PPT (Table 2, Supplementary Table 1) (24, 25).
Considerable differences were observed in the number of undetected miRNAs in K2EDTA, PPT, and PPP samples when immediately frozen at - 80°C (Table 3): specifically, the stepwise centrifugation significantly increased the frequency of undetected miRNAs compared to PPT samples (Supplementary Table 2 and 4). The 15 undetected miRNAs in the K2EDTA and PPT samples were undetected also in the PPP samples. The FC analyses showed that the percentage of up- and down-regulated miRNAs in K2EDTA as compared to PPP samples were 5% and 7%, respectively (Table 2), while in the PPT as compared to the PPP samples were 2% and 1%, respectively (Table 2).
Table 3
Storage temperature and delay: between-group analyses
The 15 miRNAs undetected in the samples immediately frozen at - 80°C were also undetected in all the samples stored for 24h at RT and 4°C.
Considerable differences in the number of undetected miRNAs were observed in the sample stored for 24h at either RT and 4°C considering the different collection conditions (Table 3, Supplementary Table 2 and 4), with the exception of samples stored for 24h at 4°C collected in K2EDTA and PPT tubes. Fold change analyses showed 2% up-regulated and 12% down-regulated miRNAs in PPP samples compared to the K2EDTA stored for 24h at RT and 1% up-regulated and 5% down-regulated miRNAs in PPP samples compared to the K2EDTA samples stored for 24h at 4°C (Table 2). No differences in miRNA expression were observed between PPT and K2EDTA samples stored for 24h at RT, whereas after storage for 24h at 4°C, 1% of miRNAs were upregulated in the PPT compared to the K2EDTA samples (Table 2).
Storage temperature and delay: within-group analyses
Analysis of miRNAs detectability revealed significant differences between samples immediately frozen at - 80°C, and stored for 24h at either RT or 4°C, collected in both K2EDTA and PPP (Table 4, Supplementary Tables 3 and 5). On the contrary, no differences in miRNA detectability emerged when comparing different storage conditions for samples collected in PPT (Table 4, Supplementary Tables 3 and 5).
Table 4
Analysis of FC revealed that the percentage of up- and down-regulated miRNAs after storage for 24h at RT, compared with samples immediately frozen at - 80°C, were, respectively, 5% and 4% for K2EDTA samples, 2% and 1% for PPP samples, 1% and 0% for PPT samples. Compared with samples immediately frozen at - 80°C, the results after plasma sample storage for 24h at 4°C showed 1% up-regulated and 1% down-regulated miRNAs in K2EDTA samples, 0% up-regulated and 1% down-regulated miRNA in PPP samples, and 0% up-regulated and 1% down-regulated miRNA in 1 PPT samples (Table 5).
Table 5
Circulating miRNA expression quantification
Among clinically relevant miRNAs associated with skeletal muscle, bone, and cancer from K2EDTA samples (hsa-miR-125b-5p, hsa-miR-29b-3p, and hsa-miR-376a-3p), PPP samples (hsa-miR-125a-5p), and PPT samples (hsa-miR-125b-5p), the circulating profile in the three storage conditions (- 80°C, 24h RT, and 24h 4°C) differed significantly (Table 5, Supplementary Table 6).
Some clinically relevant miRNAs were down-regulated in PPP samples compared to the other collection conditions. Specifically, hsa-miR-125b-5p, hsa-miR-152-3p, hsa-miR-375 were lower in PPP samples compared to K2EDTA, while hsa-miR-378a-3p compared to PPT samples, when immediately stored at - 80°C. When stored for 24h at RT, hsa-miR-1, hsa-miR-195a-5p, hsa-miR-375 were down-regulated in PPP samples compared to K2EDTA sample, and hsa-miR-378a-3p compared to PPT samples. Similarly, hsa-miR-1 was lower in PPP compared to K2EDTA samples, and hsa-miR-375 compared to PPT samples when stored for 24h at 4°C (Table 2, Supplementary Table 1). Despite platelet depletion, some miRNAs were actually up-regulated in the PPP as compared to the PPT (hsa-miR-152-3p and hsa-miR-378-3p) and to K2EDTA samples (hsa-miR-210-3p) (Table 2, Supplementary Table 1).
Noteworthy, the miRNA expression in the PPT samples were generally less affected by storage condition due to the lowest FC values for PPT samples (Table 5, Supplementary Table 6).
Discussion
Several preanalytical variables are thought to influence the quantitation of circulating miRNAs. Here, we provide an analysis of four preanalytical variables (collection matrix, platelet depletion, storage delay time and storage delay temperature) that could affect plasma circulating miRNA stability and, thus, miRNA detection and quantitation.
Regarding the sample matrix, EDTA has been demonstrated not to interfere with enzyme activities in PCR-based assays and is considered the best anticoagulant for circulating miRNA profiling (17). Hence, for our study we chose tubes spray-coated with EDTA either with or without a gel separator (PPT and K2EDTA, respectively).
Several studies have demonstrated that miRNA concentrations in EDTA-plasma samples are influenced by centrifugation time and speed, since platelets and microparticles represent a source of contamination for circulating miRNAs (24, 26). Recently, Basso et al. reported that EDTA-, citrate-, acid citrate dextrose-anticoagulated plasma and serum were all comparable and superior to heparin-plasma in terms of miRNA quality. In terms of miRNA yields, EDTA-plasma was affected by storage (over 9h at RT) and double centrifugation. Finally, EDTA-plasma miRNA expression profiles were strongly affected by the centrifugation protocol rather than temperature or time to separation (16). In another recent study, by Binderup, the impact of standard centrifugation, single prolonged centrifugation (3000xg for 30 min) and double centrifugation (2x3000xg for 15 min) on plasma miR-92, miR-126, and miR-16 expression was investigated. Differences in miR-92 and miR-126 expression, after single prolonged and double centrifugation, were evidenced when data were normalized on miR-16, miR-16-cel-39, but not when normalized on cel-39. Plasma expression of miR-92 and miR-126 fairly correlated when samples underwent to either standard (2000xg for 10 min) or double centrifugation and normalized on miR-16 (27). Besides the preanalytical-related information, this paper further highlights the issue of sample processing (5). McDonald and colleagues observed that differential processing alters miRNAs quantitation. An additional centrifugation step in plasma processing to reduce the effects of cellular contaminants significantly influenced miR-15 and miR-24 concentrations, while miR-26 concentration was only slightly reduced (28). Accordingly, for our K2EDTA samples, we applied stepwise centrifugation of plasma to assess the contribution of platelet-derived contaminants potentially affecting miRNA quantitation (PPP).
The effect of storage on plasma miRNAs stability remains debated (14, 29). Mitchell et al. demonstrated that miRNA expression (miR-15, miR-16, miR-16, and miR-24) were only minimally affected by incubation at RT for 24h, whereas Glinge demonstrated a reduction in plasma miR-1 and miR-21 expression after 24h at RT (14, 29).
To assess the effect of storage conditions on plasma circulating miRNA expression, for our study miRNA expression in plasma samples collected and immediately frozen at -80°C were compared to those collected and frozen with 24h delay stored at either RT or 4°C. This was done to simulate real-world situations in the clinical laboratory. By the end, 179 circulating miRNAs were assayed in nine different conditions: plasma from K2EDTA tubes or PPT tubes immediately frozen at - 80°C or stored for 24h at either RT or at 4°C and plasma from K2EDTA tubes depleted of platelets (PPP) stored as above.
Our results show that each condition was associated with a specific circulating miRNA profile, demonstrating that sample matrix, processing, and storage can all affect the analytical output in circulating miRNA measurement. Analyses of detectability showed that regardless of the storage condition, miRNA detectability in the PPP samples was notably affected (Table 3, Supplementary Tables 2 and 4). The frequency of undetected miRNAs in the K2EDTA compared to PPT samples was significantly reduced in each storage condition, except for 24h at 4°C in which miRNA detectability was comparable (Supplementary Tables 2 and 4). Regardless of the sample matrix, miRNA detectability was less influenced by storage conditions in the PPT samples: no differences of undetected miRNAs were observed (Table 4, Supplementary Table 3 and 5). Contrarily, K2EDTA samples displayed the greatest variability in miRNA detectability (Table 4, Supplementary Table 3 and 5): the percentage of detectable miRNAs increased during storage for 24h, particularly at 4°C. The amount of detectable miRNAs was higher in the PPT samples than in either the K2EDTA or the PPP samples, revealing that gel-mediated separation can actually increase the yield in miRNAs by contemporarily keeping cell contamination low.
According to previous studies on plasma processing and platelets contaminant removal (26-28), the stepwise centrifugation in K2EDTA samples led to reduced expression of platelet-derived miRNAs (24, 25), and also of other miRNAs of other cellular origin (16). This reduction was more evident after storage for 24h, highlighting a consistent instability of miRNAs in the PPP samples as compared to miRNA in the K2EDTA samples. The miRNA expression were higher in the PPT than the PPP tubes, while they were comparable in the PPT and the K2EDTA tubes when the plasma samples were immediately frozen at -80°C. Analyses of FC showed that miRNAs in the PPT samples were less influenced by temperature and length of storage: the PPT samples had the lowest number of modulated (up- and down-regulated) miRNAs. These results indicate that the PPT samples, possibly thanks to the low cellularity in the plasma due to gel-mediated separation, have greater stability of circulating miRNAs in different storage conditions, e.g., those found in routine laboratory practice. Hence, PPT plasma samples could be considered a suitable matrix for miRNA measurement.
Analyses of single miRNAs of clinical relevance (e.g., associated with skeletal muscle, bone, and cancer) shared the previous results about the greater stability of miRNA from samples collected in the PPT tubes, regardless of storage length and temperature. Notably, it emerged that each miRNA had a different circulating profile in each condition. This is of clinical relevance because a combination of various conditions could have a profound effect on analytical output and, hence, give rise to misinterpretation of results. Furthermore, our analysis showed that different miRNAs are differently affected by the same collection method or storage condition. This represents a great limitation for studies of pre-analytical variables, where few miRNAs are analysed, since a limited number of miRNAs cannot be representative of all circulating miRNAs. The analysis of a large amount of miRNAs, as in the present study, may offer a better overview and a more detailed description of the effects of different pre-analytical variables on the measurability of circulating miRNAs.
The main limitation of the present study is the limited number of variables that were considered. Indeed, given a potential list of preanalytical variables that can affect analytical output, here we just focused on only a few. Nonetheless, our study addresses part of most complex subject that needs further investigation. Furthermore, the rigorous experimental design provides reasonable certainty about the practical utility of the results reported here.
Circulating miRNAs are emerging as potential biomarkers for diagnostic, prognostic, and predictive purposes; meanwhile, there is a growing interest in the preanalytical variables that could affect their quantitation (2, 6). Standardization of these variables represents a critical step in the clinical implementation of miRNAs, since, as shown in this study, variables such as matrix, platelet depletion, and storage conditions can all alter miRNA measurement.