Evaluation of RNA purification methods by using different blood stabilization tubes: identification of key features for epidemiological studies
BMC Research Notes volume 13, Article number: 77 (2020)
Peripheral blood is the most promising source of RNA biomarkers for diagnostic and epidemiological studies, because the presence of disease and prognostic information is reflected in the gene expression pattern. Quality RNA is used by a number of different downstream applications, so the selection of the most appropriate RNA stabilization and purification method is important. We have analyzed the RNA purified from 300 blood samples from 25 donors processed by two technicians using three methodologies with Tempus and PaxGene tubes.
The best quality sample results were obtained with the Tempus Spin RNA Isolation Kit and the PaxGene Blood miRNA Kit, although larger amounts of RNA were obtained with the Tempus Spin RNA Isolation Kit. Lower Cq values were observed for RNA and miRNA genes in samples that were tested with PaxGene Blood miRNA Kit and Tempus Spin RNA Isolation Kit respectively. We identify the Tempus Spin RNA Isolation Kit as the most robust methodology, whilst the MagMax for Stabilized Blood Tubes RNA Isolation Kit showed the most instability. For biobanks, which process a large cohort and conduct epidemiological studies, the Tempus Spin RNA Isolation Kit is the most appropriate methodology. The study demonstrates the robustness of real-life procedures.
The major limitation associated with RNA is the phenomenon of RNA degradation during blood collection and storage . Both the sample collection and the RNA purification procedures influence the data outcome . To minimize the degradation, different RNA stabilization technologies have been developed [6, 7]. Tempus and PaxGene tubes are commonly-used for blood collection. The use of a stabilizer is mandatory if one wishes to maintain the gene expression profile. Different suppliers have released kits based that can be coupled to Tempus and/or PaxGene tubes . A number of previous studies have evaluated different RNA extraction methodologies [9,10,11], but the results are difficult to interpret because replicas of the same samples aren’t used during the comparison of methodologies. Different tubes demonstrate different levels of stabilization efficiency depending on the specific target, so the selection of the best methodology is fundamental to any reliable further analysis .
In this paper, we present an unpublished comparison of three blood RNA extraction methodologies. Although there exist research studies that have analyzed the main characteristics of these methodologies [10, 11], there are currently no studies that compare these methodologies together, and analyzed the methods reproducibility and robustness. The results must be taken into account for large scale studies and biobanking purposes.
The study design
Blood samples were collected from 25 healthy donors. The donors were recruited by the Andalusian Public Health System Biobank (Granada, Spain). Twelve samples from each donor were collected: four blood samples that were collected in PaxGene tubes (2.5 ml; Preanalytix, Switzerland), and eight samples in Tempus tubes (3 ml; by Life Technologies, Germany).
The potential influence that individual technicians may have over the outcome of the sample analysis was evaluated. Consequently, two technicians processed 150 blood tubes, each using the three extraction methodologies mentioned above. Duplicate samples were processed for the 25 donors (n = 50), to check the reproducibility of each technicians’ with regards to the three tested extraction methodologies (n = 150) (Fig. 1).
Three extraction methodologies were evaluated. For each donation, the RNA was extracted from four blood tubes by each methodology, so that duplicate samples could be processed for the 25 donors (n = 50) by each technician.
Methodology A: PaxGene® tubes and the PaxGene® Blood miRNA Kit (Preanalytix, Switzerland). The PaxGene tubes were semi-automated processed using the Qiacube robot (Qiagen) .
Methodology B: Tempus® tubes and the MagMax for Stabilized Blood Tubes RNA Isolation Kit (ThermoFisher, USA)  and a DynaMag-2 stand (Life Technologies, USA).
Methodology C: Tempus® tubes and the Tempus™ Spin RNA Isolation Kit (Applied Biosystems, USA) . No optional DNase treatment was performed.
RNA quantification and quality analysis
The RNA yield and its purity were measured using an Infinite F200 spectrophotometer (Tecan, Switzerland). The concentration was measured against the A260 nm value, while A260/A280 nm and A260/A230 ratios were used to calculating the purity.
Forty-two samples from 7 donors were measured by fluorimetry technique using the Quant-iT™ RiboGreen™ RNA Assay Kit (ThermoFisher Scientific, EEUU) (Fig. 1).
RNA integrity analysis
Samples RNA integrity was measured using the Bioanalyzer 2100 (Agilent Technologies Ind., USA) and the Agilent RNA 6000 Nano Kit (Agilent Technologies Inc.). The RNA integrity (RIN value) was calculated.
Samples from 7 donors (n = 42) were randomly selected for cDNA synthesis and expression analysis. Samples were extracted by the two technicians using the three methodologies. RNA was reverse-transcribed using miScript II RT Kit (Qiagen) (Additional file 1).
Expression assays were generated in triplicate for three miRNAs (miR16, miR26, and miR30), with the miScript SYBR Green PCR Kit (Qiagen) (Additional file 1).
Expression assays were generated in triplicate for three mRNAs (18S rRNA , ACTB , and GAPDH) ). The FastStart Essential DNA Green Master Kit (Roche LifeScience) was used for amplification (Additional file 1).
An extraordinary DNase treatment was carried out on 20 µl of selected RNAs using the RNase-Free DNase Set (Qiagen).
A statistical analysis was performed using the IBM SPSS Statistics program 17.0.3 (SPSS, Inc., Chicago, IL). A descriptive statistic analysis was also performed with the XLSTAT v2017.6 program. The ANOVA test (p < 0.05) and the ‘post hoc test after ANOVA’—Scheffé’s method—(p < 0.05) was applied to compare the mean differences between methodologies. Grubbs test was applied for outliers identification. Box plots were constructed with the XLSTAT v2017.6 program for the visualisation of replicas variability.
Samples were extracted using three methodologies. The main characteristics are summarized in Additional file 2: Table A.
A comparison of the RNA yield between methodologies and technicians’ were performed. The average yields, are indicated in Additional file 2: Table B and represented in Fig. 2a. The best results were obtained with Methodology C. Significant differences were observed between methodologies and between the two technicians with Methodologies B and C, but not with Methodology A.
The A260/A280 and A260/A230 ratios were analyzed for the three methodologies (Additional file 2: Table B). Ratios were obtained for all of the samples, except for 5 samples extracted with Methodology B, due to the low RNA concentration of the samples (2.56–16.08 ng/µl). The A260/A280 and A260/A230 values obtained are represented in Fig. 2b, c. Methodologies A and C showed the best A260/A280 ratio. No significant differences were observed between methodologies and technicians.
Regarding the A260/A230 ratio, large differences were observed between methodologies. The best values were obtained with Methodology C and the worst values were obtained with Methodology B. No significant differences were observed between the technicians for Methodologies A and C, but significant differences were observed for Methodology B (Additional file 2: Table B, Fig. 2b, c).
RIN values were calculated for all of the samples except for 13 samples (6 samples from Methodology A and 7 samples from Methodology B). The 6 samples from Methodology A but for which no RIN values was calculated correspond to samples with high RNA concentrations (118.16–172.16 ng/µl), for which the Bioanalyzer showed an anomalous profile. Those 6 RNA samples were subjected to an extraordinary DNase treatment. The RIN was analyzed again obtaining a good result (Additional file 3). In addition, the 6 original and the 6 aliquots treated with DNase were checked by PCR for the GAPDH gene . PCR amplification was observed in the 6 original samples, whilst the 6 purified samples showed no amplification (data not shown).
The average RIN value and standard deviation are indicated in Additional file 2: Table B and Fig. 2d. The best results were obtained with Methodologies A and C, and no significant differences were observed between methodologies and technicians.
Analysis of the reproducibility and robustness of the three methodologies
Variability between replicas from the same donor with respect to the three methodologies and both technicians was analysed and the percentage of variation between each pair of samples from each donor that was processed by each technician with each methodology was calculated. The best reproducibility was observed with Methodology C (Additional file 2: Table C and Fig. 3). Only significant differences were observed between technicians for variable “yield” with Methodology C. Technician A obtained better reproducibility results than Technician B with Methodology A. The lowest reproducibility was obtained for both technicians with Methodology B.
RNA expression analysis
Forty-two selected samples were used to check the utility of the extracted RNA, considering Ribogreen quantification values. The average Cq and standard deviation values were calculated (Additional file 4). Gene expression was detected in every case. The worst results for mRNA were obtained with methodology A, with significant differences for the 18S rRNA gene and ACTB gene fragments. However, no significant differences were observed for Methodology B and Methodology C.
The best results for miRNA were obtained with Methodology A (especially for mir-30 miRNA), while higher Cq values were observed for Methodology C.
Previous comparative studies have been performed on different methodologies . However, a simultaneous comparison of the proposed methodologies has not been published before. In our study, Methodology C was the most effective with respect to yield, so it provided a yield between 38% and a 56% greater in comparison to Methodologies B and A, respectively. Other studies have shown a higher RNA yield with Tempus tubes [10,11,12, 19], associated with a smaller volume of blood in the case of PaxGene tubes . However we note that Methodology B also uses Tempus tube and the obtained RNA yield was lower than for Methodology C, so it’s not only important the type of tube, but also the RNA extraction kit. The high degree of heterogeneity that was observed between the donors, with values ranging between 2.88–22.39, 0.43–19.70, and 0.76–24.98 µg for methodologies A, B, and C respectively, is mainly due to the intrinsic characteristics of the samples.
With respect to the purity of the RNA, the best results were obtained with Methodology C. Overall, ‘very good’ and comparable purity values were observed for Methodology A and Methodology C, in accordance with previous research [9, 12, 20]. However, we observed more unstable A260/A230 ratios with Methodology A and especially for Methodology B than these previous studies. The presence of impurities with Methodology B has been described previously . The differences observed in the A260/A230 ratio between the two technicians were very marked for Methodology B, possibly due to small differences in the protocol implementation. It has been reported that the observed differences in A260/A230 ratio with Methodology A, are due to the high salt content of PaxGene kit elution buffer .
Methodologies A and C showed extraordinary RNA integrity, not observed in other studies [10,11,12]. It was not possible to obtain the RIN value for seven samples extracted by using Methodology B due to the low concentration of the samples. The RIN value was also not available for 6 samples that were extracted with Methodology A, because of the anomalous profile. The fact that this occurs with some frequency with paired samples from the same donor suggests some peculiarity of the sample, and Methodology A wasn’t able to eliminate the contaminant. The bioanalyzer fluorogram only shows nucleic acids those are present in the sample . Using an extraordinary DNase treatment and PCR we demonstrated that DNA traces are present, and the kit hasn’t been able to eliminate in during RNA purification procedure, probably because of the high RNA concentration, due to a high cellularity count and a consequent large amount of DNA.
Additionally, Methodology C show greatest robustness between replicas and technicians. The differences observed between technicians, may be attributed to a better understanding and/or execution of a particular protocol. Clearly, there was a huge variability in the purity and integrity measures with Methodology B, possibly due to a non-efficient contaminants removal.
Slightly worse expression results were obtained for mRNAs with Methodology A, although this methodology showed a higher detection capacity for miRNAs. The purification capacity of miRNA for Methodology C should be noted, although the manufacturer doesn’t guarantee the purification of miRNA. Skogholt et al. suggest that differences in expression between methodologies are due to sampling systems, influenced by physical properties of the transcripts . Although 6 transcripts were studied in the present study, we observe correlation between methodologies that use PaxGene and the Tempus tubes with respect to the expression levels of mRNA. However, there is an effect of the extraction protocol on the expression of miRNA, so not only the type of tube is important [11, 19], and the selected methodology influences the mRNA and miRNA expression profiles. It is essential to establish the standardization of sample processing for large cohort studies, such data from mixed methodologies should never be analyzed together [11, 19, 23].
To summarize, superior results were obtained with Methodologies A and C. For Biobank large-cohort and epidemiological studies, Methodology C is the most appropriate methodology. However, the reduced hands-on time that is needed for Methodology A must also be considered.
DNase treatment was included in Methodologies A and B; however the additional treatment with DNase was not performed with Methodology C.
Availability of data and materials
The datasets used and/or analyzed during this study are available from the corresponding author on request.
Human glyceraldehyde-3-phosphate dehydrogenase
- 18S rRNA:
Ribosomic RNA 18S
Hebels DG, Georgiadis P, Keun HC, Athersuch TJ, Vineis P, Vermeulen R, Portengen L, Bergdahl IA, Hallmans G, Palli D, et al. Performance in omics analyses of blood samples in long-term storage: opportunities for the exploitation of existing biobanks in environmental health research. Environ Health Perspect. 2013;121:480–7.
Skogholt AH, Ryeng E, Erlandsen SE, Skorpen F, Schonberg SA, Saetrom P. Gene expression differences between PAXgene and Tempus blood RNA tubes are highly reproducible between independent samples and biobanks. BMC Res Notes. 2017;10:136.
Pantel K, Alix-Panabieres C. Circulating tumour cells in cancer patients: challenges and perspectives. Trends Mol Med. 2010;16:398–406.
Kim JH, Jin HO, Park JA, Chang YH, Hong YJ, Lee JK. Comparison of three different kits for extraction of high-quality RNA from frozen blood. Springerplus. 2014;3:76.
Meyer A, Paroni F, Gunther K, Dharmadhikari G, Ahrens W, Kelm S, Maedler K. Evaluation of existing methods for human blood mRNA isolation and analysis for large studies. PLoS ONE. 2016;11:e0161778.
Chai V, Vassilakos A, Lee Y, Wright JA, Young AH. Optimization of the PAXgene blood RNA extraction system for gene expression analysis of clinical samples. J Clin Lab Anal. 2005;19:182–8.
Rainen L, Oelmueller U, Jurgensen S, Wyrich R, Ballas C, Schram J, Herdman C, Bankaitis-Davis D, Nicholls N, Trollinger D, Tryon V. Stabilization of mRNA expression in whole blood samples. Clin Chem. 2002;48:1883–90.
Tan SC, Yiap BC. DNA, RNA, and protein extraction: the past and the present. J Biomed Biotechnol. 2009;2009:574398.
Asare AL, Kolchinsky SA, Gao Z, Wang R, Raddassi K, Bourcier K, Seyfert-Margolis V. Differential gene expression profiles are dependent upon method of peripheral blood collection and RNA isolation. BMC Genomics. 2008;9:474.
Duale N, Brunborg G, Ronningen KS, Briese T, Aarem J, Aas KK, Magnus P, Stoltenberg C, Susser E, Lipkin WI. Human blood RNA stabilization in samples collected and transported for a large biobank. BMC Res Notes. 2012;5:510.
Hantzsch M, Tolios A, Beutner F, Nagel D, Thiery J, Teupser D, Holdt LM. Comparison of whole blood RNA preservation tubes and novel generation RNA extraction kits for analysis of mRNA and MiRNA profiles. PLoS ONE. 2014;9:e113298.
Menke A, Rex-Haffner M, Klengel T, Binder EB, Mehta D. Peripheral blood gene expression: it all boils down to the RNA collection tubes. BMC Res Notes. 2012;5:1.
PAXgene Blood miRNA Kit Handbook. Qiagen-Preanalytix ed.; 2015.
MagMAX™ for Stabilized Blood Tubes RNA Isolation Kit Protocol. In: Part Number 4452007 Rev A (Technologies A-L ed., vol. Part Number 4452007 Rev. A; 2010.
Tempus™ Blood RNA Tube and Tempus™ Spin RNA Isolation Kit Protocol. D PNR ed. Applied Biosystems; 2008.
Lopez P, Wagner KD, Hofman P, Van Obberghen E. RNA activation of the vascular endothelial growth factor gene (VEGF) promoter by double-stranded rna and hypoxia: role of noncoding VEGF promoter transcripts. Mol Cell Biol. 2016;36:1480–93.
Boutz PL, Chawla G, Stoilov P, Black DL. MicroRNAs regulate the expression of the alternative splicing factor nPTB during muscle development. Genes Dev. 2007;21:71–84.
Kang H, Rho JG, Kim C, Tak H, Lee H, Ji E, Ahn S, Shin AR, Cho HI, Huh YH, et al. The miR-24-3p/p130Cas: a novel axis regulating the migration and invasion of cancer cells. Sci Rep. 2017;7:44847.
Nikula T, Mykkanen J, Simell O, Lahesmaa R. Genome-wide comparison of two RNA-stabilizing reagents for transcriptional profiling of peripheral blood. Transl Res. 2013;161:181–8.
Feddersen S, Bastholt L, Pedersen SM. Stabilization of circulating thyroglobulin mRNA transcripts in patients treated for differentiated thyroid carcinoma. Ann Clin Biochem. 2017;54:558–66.
Richards J, Unger ER, Rajeevan MS. Simultaneous extraction of mRNA and microRNA from whole blood stabilized in tempus tubes. BMC Res Notes. 2019;12:39.
Grissom SF, Lobenhofer EK, Tucker CJ. A qualitative assessment of direct-labeled cDNA products prior to microarray analysis. BMC Genomics. 2005;6:36.
Yip L, Fuhlbrigge R, Atkinson MA, Fathman CG. Impact of blood collection and processing on peripheral blood gene expression profiling in type 1 diabetes. BMC Genomics. 2017;18:636.
We would like to thank Susana Alonso Fernández for her economic management work and for her help in buying the necessary materials and in publishing the laboratory results. We would also like to thank Dr. Robert Ryan for the thorough review of the language in the article.
This work was supported by funds from Consejería de Salud (Junta de Andalucía) (Grant Number: 0001) which supported the hiring of the personnel in charge of the study, and from the Instituto de Salud, Carlos III (PT13/0001/0040–PT13/0010/0050, both co-financed by FEDER) that supported the acquisition of the necessary consumable materials. The funding bodies had no role in the design of the study, data collection, analysis, interpretation of the data, or writing of the manuscript.
Ethics approval and consent to participate
The study protocol was carried out in accordance with the Helsinki Declaration. The study was explained to the donors, who participated voluntarily. The handling of human biological samples was carried out according to the national legal framework (Spanish Law on Biomedical Research (July 2007)) following informed consent of the donors and following Biobank procedures. The Scientific Committee of the Andalusian Public Health System Biobank gave written consent for the transfer of the samples that were used in the project (Date 27/11/2015). The local ethics committee “Comité de ética para la investigación biomedica”, granted written consent and approved the procedures performed in this work (Project code: 37150070. Date 28/09/2015).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
(i) Reversotranscription, (ii) mRNA and miRNA real time PCR conditions, (iii) and oligo sequences details.
Table A. Summary of the technical characteristics, technical requirements and cost analysis for the three methodologies tested. Table B. Average and standard deviation obtained for RNA yield, A260/A280, A260/A230 and RIN obtained by both technicians with the three methodologies. The average values encompassing the values obtained by both technicians are also indicated. Table C. Robustness of each methodology analysis. Average percentage and standard deviation for variation between methodologies and technicians for variables yield, A260/A280, A260/A230 and RIN.
Diagram of an RNA electropherogram for six samples purified with Methodology A. Images a–f correspond to the RNA purified by Methodology A, following the manufacturer’s instructions. The images a’–f’ correspond to the same RNA samples after additional treatment with DNase. a and a’: Donor 16 Sample 1, b and b’: Donor 16 Sample 2, c and c’: Donor 17 sample 1, d and d’: Donor 23 Sample 1, e and e’: Donor 23 Sample 2, f and f’: Donor 23 Sample 3. The anomalous peaks of the electropherograms are shaded in the picture.
Representation of the average Cq value for the RNA obtained by the three methodologies. mRNA genes (18S rRNA, ACTB, and GAPDH) and miRNA (mir-26, mir-30, and mir-16) were tested. Error bars are indicated.
About this article
Cite this article
Carrillo-Ávila, J.A., de la Puente, R., Catalina, P. et al. Evaluation of RNA purification methods by using different blood stabilization tubes: identification of key features for epidemiological studies. BMC Res Notes 13, 77 (2020). https://doi.org/10.1186/s13104-020-04943-4
- RNA quality
- Liquid biopsy