mir-16-5p as a Suitable Reference Gene for Normalization of Quantitative Real Time PCR in Acute Lymphoblastic Leukemia

Samiah Shahid1, Jawaria Shaheen1, Wajeehah Shahid2, M. Waheed Akhtar1,3 and Saima Sadaf1,*

1Institute of Biochemistry and Biotechnology, University of the Punjab, Lahore-54590

2Department of Physics, The University of Lahore, Lahore

3School of Biological Sciences, University of the Punjab, Lahore-54590

ABSTRACT

Quantitative real time polymerase chain reaction PCR (qPCR) may be utilized as a sensitive and reliable technique for the determination of circulating miRNA expression. Despite recent advancements, there is not a single consensus on the use of reference gene for quantification of circulating miRNAs through qPCR analyses in ALL. In the current study, we identified the reference gene that is the most suitable for qPCR normalization in patient and control plasma samples of ALL. Three highly reported reference genes namely RNU6-2, mir-16-5p and cel-mir-39-1 were selected as the candidate genes for normalization. Preliminary, quantification was performed through real time PCR in eight samples. The geNorm algorithm and comparative delta Cq method were used for selection of suitable reference gene out of the three candidate genes. The validation studies were, thereafter, performed on 112 samples including 87 patients and 25 normal healthy controls. The geNorm algorithm exhibited circulating mir-16-5p with the highest expression stability demonstrated by geNorm M value of 0.418. The comparative delta Cq method showed no significant difference in Cq values of ALL patients and healthy controls suggesting that mir-16-5p is the most stable normalizer for miRNA expression profiling assays in ALL. The results revealed that circulating mir-16-5p may be utilized as a suitable reference gene for normalization and improved quantification of miRNA.

Article Information

Received 26 August 2018

Revised 30 September 2018

Accepted 25 October 2018

Available online 01 March 2019

Authors’ Contribution

MWA and SS designed the experiment, supervised the research and reviewed the manuscript. SS and WS collected the clinical samples and statistically analyzed the study. SS and JS performed the experiments. SS prepared the manuscript.

Key words

mir-16, qPCR, Normalization, Acute lymphoblastic leukemia, Reference gene.

DOI: http://dx.doi.org/10.17582/journal.pjz/2019.51.2.747.754

* Corresponding author: saima.ibb@pu.edu.pk

0030-9923/2019/0002-0747 $ 9.00/0

Copyright 2019 Zoological Society of Pakistan

Introduction

MicroRNAs (miRNAs) are small non-coding RNAs, discovered first time in c. elegans in 1993 (Lee et al., 1993), play the major role in the regulation of gene expression either by suppressing their target messenger RNA (mRNA) and/or by up regulating the target genes in the 3′-untranslated region (3′-UTR) (Iorio and Croce, 2012; Reid et al., 2012). They are involved in different biological and molecular mechanisms like cell proliferation and stem cell formation (Hasani et al., 2014; Qin et al., 2018). Several miRNAs are engaged in the regulatory processes of hematopoiesis as they regulate differentiation of hematopoietic stem cells at different stages (Marcucci et al., 2011). A strong relationship between miRNA expression level and acute lymphoblastic leukemia (ALL) has been elucidated. ALL is the most common childhood malignancy worldwide including Pakistan with peak incidence at 6 years of age. The diagnosis of ALL involves bone marrow biopsy for which the children become reluctant due to its invasiveness. So the diagnosis of ALL demands the search for non-invasive techniques than can be utilized in developed and under-developed countries like Pakistan. Studies have demonstrated the potential of miRNAs in the diagnosis and prognosis of ALL. However, accurate miRNA quantification remained challenging in miRNA expression profiling studies. The miRNA profiling by quantitative real-time PCR (qPCR) is relatively more specific, reproducible, easily available and cost effective. It also requires very little amount of sample, even a single cell is sufficient for miRNA profiling. As miRNA expression levels are low in circulating fluid like plasma, so qPCR is well adapted for plasma miRNA profiling (El-Halawani et al., 2014). It has been considered a robust method for quantification of miRNA in ALL. However, quantification through qPCR demands normalization for accurate quantification of miRNA expression (Roberts et al., 2014; Schwarzenbach et al., 2015). Normalization controls for biological and technical variations in the qPCR assays. Several reports have been published on the selection of reference gene for different experiments; still there is no harmony on the use of reference gene for normalization of miRNA quantification using real-time PCR in ALL. The present research was established to determine suitable reference genes as a normalizer in qPCR assays for profiling circulating miRNA in ALL. We determined the stability of three reference genes including RNU6, mir-16-5p and c. eleg-mir-39 in plasma of ALL. We observed that circulating mir-16-5p exhibited highest stable expression in patients and control samples of ALL rendering it the reference gene most suitable for qPCR normalization in miRNA profiling.

 

Materials and methods

Study patients, controls and plasma samples

In this study, blood samples were obtained from 112 study subjects including patients (n=87) and normal controls (n=25). Patients registered at INMOL (Institute of nuclear medicine and oncology), Lahore, Pakistan from January, 2013 to November, 2016 and newly diagnose with ALL were included in this study. Normal controls were healthy individuals who have not any malignancy or infectious disease like hepatitis, HIV and tuberculosis. ALL patients and controls those have gone through any medical or surgical treatment were excluded from the study. Plasma was isolated within two to three hours after sample collection and stored at -80°C in aliquots until further use. Informed written consent was taken from the enrolled subjects to draw blood and use for the current research. The research protocols were approved by the review board of the ethical committee of School of Biological Sciences, University of the Punjab, Lahore, Pakistan with reference number 284/15 according to the rules of the Helsinki declaration approving human related research.

miRNA extraction from plasma

The total RNA containing miRNA was extracted from plasma using miRNA easy kit for serum/plasma from Qiagen (Hilden, Germany) according to the instructions from the manufacturer. The total RNA containing miRNA was stored quickly at -20°C till further processing.

Quantification of miRNA

The quantity of eluted miRNA in the total RNA was evaluated by NanoDrop (ND-2000) from Thermo scientific (USA). The optical density (OD) was taken at 260 nm against RNase free water blank.

Selection of reference gene for normalization of qPCR experiments

Since normalization is the prerequisite for quantification of miRNAs in miRNA expression profiling, we designed our experiment to select for the suitable reference gene for miRNA profiling in ALL. For this purpose, preliminary experiment was performed on eight samples by taking three frequently reported reference genes including RNU6-2 (snRNA), mir-16-5p and cel-mir-39-1. The validation studies for the expression stability of the most stable gene were performed in 112 plasma samples of ALL. Stability of expression for all the candidate reference genes was measured by qPCR and the GeNorm M (average stability of expression) was calculated through qBase+ software (BioGazelle, Zwijnaarde, Belgium).

Primer assays of candidate genes

The sequences of primers for candidate genes were obtained from miRNA database (www.mirbase.org). Primer assays (10X) were obtained from Qiagen (Hilden, Germany). Details on primer sequences, type and their accession number are provided in Table I.

Reverse transcription

The cDNA of total RNA containing miRNA was synthesized by reverse transcription reactions utilizing miScript II RT Kit from Qiagen (Hilden, Germany) following the protocols established by the company. The cDNA synthesized was kept at -20ºC for further analysis.

miRNA expression analysis by relative qPCR

The quantification cycle (Cq) of three reference genes including RNU6, cel-mir-39-1 and mir-16-5p was determined through real time PCR using SYBR green dye in miScript PCR kit from Qiagen (Hilden, Germany). The template used for real time PCR was cDNA synthesized by reverse transcription. The reactions were performed using real time thermal cycler CFX96TM system from BioRad®. Exogenous control i.e. c. elegans mir-39 was spiked in during RNA extraction as a normalization control and NTC (no template control) was considered as the negative control. All reactions were performed in duplicate and to check inter-run variations, two independent samples were run in all qPCR experiments.

 

Table I.- Primers of candidate reference genes with their type, mature sequence and accession numbers.

Candidate gene	Type	Mature sequence	Accession number
miR-16-5p	Hsa-miRNA	5′-UAGCAGCACGUAAAUAUUGGCG-3′	MIMAT0000069
cel-miR-39-1	C. eleg-miRNA	5′-UCACCGGGUGUAAAUCAGCUUG-3′	MIMAT0000010
RNU6-2	SnRNA	5′-UAGCAG-CACGUAAAUAUUGGCG-3′	NR002752

 

Table II.- Demographic and clinical characteristics of the enrolled study subjects (n=112).

Characteristics	ALL patients (n=87)	Normal controls (n=25)	P-value
Age (years)	18.18±10.75 (4.00-48.0)	22.96±14.01 (4.0-47.0)	0.1831
Height (cm)	149.4±27.59 (39.0-185)	148.8±23.45 (38.5-190)	>0.05
Weight (kg)	48.04±20.30 (13.0-100)	47.95±15.63 (13.5-98.2)	>0.05
Hemoglobin (g/dl)	7.45±2.09 (2.90-11.9)	13.5±0.79 (12.3-15.6)	<0.0001
RBC’s (106/µL)	2.85±0.87 (1.36-5.40)	4.36±0.51 (3.4-5.6)	<0.0001
HCT (%)	26.34±9.80 (10.9-86.8)	41.31±2.55 (37.6-45.8)	<0.0001
MCV (fl)	81.17±15.09 (21.1-110)	86.90±5.36 (75.3-98.1)	0.0282
MCH (pg)	23.62±5.10 (10.5-32.2)	30.3±2.67 (25.3-35.2)	<0.0001
MCHC (g/dL)	33.76±7.50 (17.0-90.0)	33.63±1.67 (30.9-36.2)	0.0489
WBC’s (109/L)	127.9±111.1 (6.0-481)	8.22±1.87 (5.3-11.1)	<0.0001
Lymphocytes (%)	77.12±9.08 (43.5-90.3)	31.07±8.70 (15.0-49.1)	<0.0001
Mixed cells (%)	4.82±2.69 (1.10-15.0)	18.62±10.42 (4.20-44.0)	<0.0001
Neutrophils (%)	18.06±8.56 (4.00-44.0)	50.30±10.05 (39.5-71.3)	<0.0001
Platelets (109/L)	47.23±43.89 (1.70-224)	299.7±81.45 (165-456)	<0.0001

The values are mean±SD (standard deviation). P-values are significant at <0.05.

 

Melt curve analysis

Melt curve analysis was performed to check the specificity of the amplification products. Melt curves were run from 65°C to 95°C using CFX96TM real time thermal cycler from BioRad®.

Statistical analyses

Raw data of fluorescent signals was analyzed by BioRad CFX96TM manager software (Version 3.1) with baseline subtracted curve fit and threshold set at 10.0 for Cq determination. The mean Cq values were calculated using duplicate Cq values. The data was analyzed using statistical tools including SPSS statistics 21 and prism software with GraphPad, version 7.0. Student’s t test was used for validation and Mann Whitney test was applied to compare two groups. All P values were two-sided, and the values less than 0.05 reflected statistically significant results.

 

Results

Study subjects characteristics

The study consisted of 112 subjects including 87 patients and 25 normal controls. There were 69 males and 18 females, and 16 males and 9 females in patients and control samples, respectively. The minimum age was 4 years in ALL patients and controls and the maximum age was 48 and 47 years in ALL patients and normal controls, respectively. There were 69% children and 31% adult cases in ALL patients. The basic demographic and clinical features of the enrolled subjects are presented in Table II. There was no significant difference in age, height and weight of patients and healthy controls (P>0.05). However, significant differences were observed in clinical features of the patients and controls (P<0.05).

Quantification of miRNA

The quantification of miRNA performed using NanoDrop showed concentration of miRNA at 260 nm and 280 nm. The optical density (OD) at A260/A280 gives an estimation of the RNA purity. The RNA with absorbance of A260/A280 ratio from 1.8 to 2.1 was used in the downstream processing.

 

Relative quantification revealed stable expression of circulating mir-16-5p

In our preliminary experiment, relative quantification was performed on 8 samples by taking three frequently reported reference genes including RNU6-2, and mir-16-5p and cel-mir-39-1. All samples were measured in the same run for the candidate reference targets to control for inter-run variation. Cq (quantification cycle) values of all three candidate genes were presented in Figure 1.

 

 

Table III.- Quantification cycle (Cq) values of candidate reference genes in ALL.

Candidate gene	Cq Min.	Cq Max.	Cq	Cq patient	Cq control	P-value
RNU6-2	22.78	23.76	23.16±0.29	23.14±0.18	23.18±0.41	0.8652
cel-miR-39-1	24.4	25.64	25.18±0.46	25.31±0.41	25.06±0.54	0.4857
miR-16-5p	25.01	25.4	25.13±0.11	25.17±0.16	25.10±0.06	0.8857

The values are mean±SD (standard deviation). P-values are significant at <0.05.

 

The mean Cq values of circulating mir-16-5p showed insignificant difference in patients and control samples (P >0.05) (Fig. 2C). Similarly, mean Cq values for RNU6-2 and cel-mir-39-1 showed P >0.05, showing non-significant difference between two groups including patient and control (Fig. 2A, B).

Results of comparative delta Cq method

The mean Cq ± standard deviation of RNU6-2, mir-16-5p and cel-mir-39-1 were calculated using comparative delta Cq method that showed mean values of 23.16±0.2963, 25.13±0.1193 and 25.18±0.4669, respectively. The mir-16-5p represented lowest standard deviation, exhibiting highest suitability for normalization. There was no statistically significant difference in mean Cq values of patients and control samples (P>0.05) as described in Table III.

Specificity and efficiency of reference genes

Melt curve analysis showed a single peak for reference genes tested demonstrating specificity of the amplified products in the present assay with 100% efficiency (Fig. 3A). Relative fluorescence unit (RFU) f more than 103 was observed (Fig. 3B).

Genorm M values of candidate reference genes analyzed by qBase+ software

qBase+ software was used to analyze quantification results that uses geNorm algorithm to calculate expression stability (GeNorm M) and variation in using n+1 or n reference genes (GeNorm V) to evaluate the least number of reference genes required for appropriate normalization. The geNorm M of 1.5 was used as a threshold value and the lower M score identified the genes with the most stable expression. The mean expression levels of RNU6-2, mir-16-5p and cel-mir-39-1 were measured and average expression stability was calculated using M value that provides ranking of the candidate reference genes showing the lowest and highest stability in expression. Ranks of the three candidate genes according to expression stability were shown (Fig. 4). Mir-16-5p revealed most stable expression in plasma of ALL and healthy controls representing geNorm M score of 0.418 as compared to RNU6-2 and cel-mir-39-1 with M value of 0.437 and 0.483, respectively.

 

Circulating mir-16-5p as a suitable reference gene for normalization

The comparative delta cq method, geNorm algorithm and relative quantification demonstrated mirr-16-5p as the candidate reference gene showing potential to serve as the normalizer for qPCR assays. On the basis of these results, mir-16-5p was proposed as the suitable normalizer gene for miRNA profiling through qPCR analysis.

 

Discussion

In this study, three candidate genes, RNU6-2, mir-16-5p and cel-mir-39-1 were analyzed for normalization of miRNA quantification in plasma samples of ALL and healthy controls as normal subjects. It was revealed that mir-16-5p showed non–significant difference in its expression in ALL patients and healthy controls (P>0.05) that demonstrated its stable expression in all samples (Table III). Further, geNorm analysis showed mir-16-5p lowest M score of 0.418 with highest stability. It has been reported that a reference gene must show stable expression in all the samples analyzed with no significant difference in expression (Das et al., 2016). We identified mir-16-5p as a suitable reference gene out of three candidate genes for normalization of qPCR analysis. Different strategies have been established to determine expression of circulating miRNAs including microarray, quantitative real-time PCR (q-PCR) (Yang et al., 2009), northern blotting (Lu et al., 2005), flow cytometry (Fichtlscherer et al., 2010), and deep sequencing (Wu et al., 2011). qPCR is well adapted for miRNA expression levels due to its high specificity, sensitivity and reproducibility (Mestdagh et al., 2008). As miRNA is present in a very little amount in circulation, so qPCR is best utilized for analyzing miRNA profiling in plasma. It requires a small amount of RNA sample as little as a single cell, while greater amount of sample (>1 μg) is required for profiling with other methods including microarray (El-Halawani et al., 2014). The qPCR has been reported as the gold standard method for quantitative expression profiling of miRNAs in serum/plasma and tissues as well. Normalization of qPCR analysis is mandatory to acquire accurate results in miRNA expression profiling. Although, exogenous miRNA spiked-in during extraction procedure endows check for technical variations, however, it is unable to normalize miRNA concentration in plasma due to sample to sample variation in patients and healthy controls. Therefore, selection of a suitable reference gene is mandatory to accurately interpret miRNA expression and improve its reliability and reproducibility. There is extensive literature on selection of reference genes for normalization of qPCR experiments in cells and tissues (Han et al., 2014; Rihani et al., 2013; Wotschofsky et al., 2011). Previously, RNU6-2 and RNU43 have been reported as normalizer in miRNA expression studies (Zhu et al., 2012). The candidate genes selected as reference for normalization as highly similar mature paralogous sequences do not exist for these particular miRNAs. A study on circulating miRNA profiling in ALL reported U6 as a reference gene (Swellam and El-Khazragy, 2016). However, there are certain limitations with respect to stability and abundance, in the use of snRNA and rRNA for expression profiling of circulating miRNAs. Recent studies reported that RNU6 is not a suitable reference gene for normalization of miRNA expression and the use of stably expressed miRNA has been suggested (Mase et al., 2017; Xiang et al., 2014). Our results also showed that stability of RNU6-2 was lower than mir-16-5p. A number of miRNAs have been reported as reference genes for normalization of qPCR analysis of circulating miRNA including miR-16, mir-93 in gastric carcinoma (Song et al., 2012), mir-191 and mir-148 in breast cancer (Hu et al., 2012). A study on plasma miRNA in B-ALL reported the use of multiple reference genes including mir-106a, mir-26b, mir-30a and mir-30d (Luna-Aguirre et al., 2015). Despite recent improvements in miRNA quantification strategies, there is no study reporting the use of reference gene for circulating miRNA profiling in ALL We designed our experiment to choose appropriate reference gene in qPCR normalization for miRNA expression studies in ALL. We identified mir-16-5p as the most suitable reference gene out of three candidate reference genes including RNU6-2, mir-16-5p and cel-mir-39-1 for qPCR normalization. Comparative delta Cq method identified mir-16-5p as the most stable normalizer gene for miRNA profiling in ALL. Furthermore, qBase+ was used to analyze quantification results (Hellemans et al., 2007) that uses geNorm algorithm to calculate expression stability (GeNorm M) (St-Pierre et al., 2017). The mir-16-5p showed lowest M value representing highest stability. The use of mir-16 has previously been described as a suitable reference genes for miRNA expression profiling in prostate cancer (Zhao et al., 2018) and cancer cell lines (Rihani et al., 2013). The mir-16 has also been reported as a reference gene to calculate normalized expression of plasma samples in chronic kidney disease (Lange et al., 2017). The present results proved the potential of mir-16-5p as a suitable reference gene to normalize qPCR analyses for determination of miRNA expression levels in ALL.

 

Conclusion

In our study, mir-16-5p showed highest expression stability in plasma of ALL and normal controls out of three previously reported reference genes, RNU6-2, cel-mir-39-1 and mir-16-5p. This is the first study that reported mir-16-5p with highest suitability to function as reference gene for qPCR normalization experiments in miRNA expression profiling of ALL.

 

Acknowledgements

The authors acknowledge the higher education commission (HEC), Pakistan for providing partial funding support to accomplish the present research.

 

Statement of conflict of interest

The authors declare that there is no conflict of interests regarding the publication of this article.

 

References

Das, M.K., Andreassen, R., Haugen, T.B. and Furu, K., 2016. Identification of endogenous controls for use in miRNA quantification in human cancer cell lines. Cancer Genom. Proteom., 13: 63-68.

El-Halawani, N., Hamed, N., Eldafrawi, M., Mourad, Z., Sorour, A. and Ghallab, O., 2014. Diagnostic and prognostic value of plasma level of microRNA-92a in acute myeloid leukemia. Am. J. mol. Biol., 4: 1-10. https://doi.org/10.4236/ajmb.2014.41001

Fichtlscherer, S., de Rosa, S., Fox, H., Schwietz, T., Fischer, A., Liebetrau, C., Weber, M., Hamm, C.W., Röxe, T. and Müller-Ardogan, M., 2010. Circulating microRNAs in patients with coronary artery diseasenovelty and significance. Cir. Res., 107: 677-684. https://doi.org/10.1161/CIRCRESAHA.109.215566

Han, H.S., Jo, Y.N., Lee, J.Y., Choi, S.Y., Jeong, Y., Yun, J. and Lee, O.J., 2014. Identification of suitable reference genes for the relative quantification of microRNAs in pleural effusion. Oncol. Lett., 8: 1889-1895. https://doi.org/10.3892/ol.2014.2404

Hasani, S.S., Hashemi, M., Eskandari-Nasab, E., Naderi, M., Omrani, M. and Sheybani-Nasab, M., 2014. A functional polymorphism in the miR-146a gene is associated with the risk of childhood acute lymphoblastic leukemia: A preliminary report. Tumor Biol., 35: 219-225. https://doi.org/10.1007/s13277-013-1027-1

Hellemans, J., Mortier, G., De Paepe, A., Speleman, F. and Vandesompele, J., 2007. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol., 8: R19. https://doi.org/10.1186/gb-2007-8-2-r19

Hu, Z., Dong, J., Wang, L.E., Ma, H., Liu, J., Zhao, Y., Tang, J., Chen, X., Dai, J. and Wei, Q., 2012. Serum microRNA profiling and breast cancer risk: The use of miR-484/191 as endogenous controls. Carcinogenesis, 33: 828-834. https://doi.org/10.1093/carcin/bgs030

Iorio, M.V. and Croce, C.M., 2012. microRNA involvement in human cancer. Carcinogenesis, 33: 1126-1133. https://doi.org/10.1093/carcin/bgs140

Lange, T., Stracke, S., Rettig, R., Lendeckel, U., Kuhn, J., Schlüter, R., Rippe, V., Endlich, K. and Endlich, N., 2017. Identification of miR-16 as an endogenous reference gene for the normalization of urinary exosomal miRNA expression data from CKD patients. PLoS One, 12: e0183435. https://doi.org/10.1371/journal.pone.0183435

Lee, R.C., Feinbaum, R.L. and Ambros, V., 1993. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell, 5: 843-854. https://doi.org/10.1016/0092-8674(93)90529-Y

Lu, J., Getz, G., Miska, E.A., Alvarez-Saavedra, E., Lamb, J., Peck, D., Sweet-Cordero, A., Ebert, B.L., Mak, R.H. and Ferrando, A.A., 2005. MicroRNA expression profiles classify human cancers. Nature, 435: 834. https://doi.org/10.1038/nature03702

Luna-Aguirre, C.M., de la Luz Martinez-Fierro, M., Mar-Aguilar, F., Garza-Veloz, I., Treviño-Alvarado, V., Rojas-Martinez, A., Jaime-Perez, J.C., Malagon-Santiago, G.I., Gutierrez-Aguirre, C.H. and Gonzalez-Llano, O., 2015. Circulating microRNA expression profile in B-cell acute lymphoblastic leukemia. Cancer Biomark, 15: 299-310. https://doi.org/10.3233/CBM-150465

Marcucci, G., Mrozek, K., Radmacher, M.D., Garzon, R. and Bloomfield, C.D., 2011. The prognostic and functional role of microRNAs in acute myeloid. Blood, 117: 1121-1129. https://doi.org/10.1182/blood-2010-09-191312

Mase, M., Grasso, M., Avogaro, L., D’Amato, E., Tessarolo, F., Graffigna, A., Denti, M.A. and Ravelli, F., 2017. Selection of reference genes is critical for miRNA expression analysis in human cardiac tissue. A focus on atrial fibrillation. Scient. Rep., 7: 41127. https://doi.org/10.1038/srep41127

Mestdagh, P., Feys, T., Bernard, N., Guenther, S., Chen, C., Speleman, F. and Vandesompele, J., 2008. High-throughput stem-loop RT-qPCR miRNA expression profiling using minute amounts of input RNA. Nucl. Acids Res., 36: e143. https://doi.org/10.1093/nar/gkn725

Qin, L., Zhang, G., Zhu, C., Wu, J., Zhao, Z., and Zhao, Y., 2018. The predicted target genes of miR-122/449a validation and effects on sertoli cells proliferation. Pakistan J. Zool., 50: 273-281. http://dx.doi.org/10.17582/journal.pjz/2018.50.1.273.281

Reid, J.F., Sokolova, V., Zoni, E., Lampis, A., Pizzamiglio, S., Bertan, C., Zanutto, S., Perrone, F., Camerini, T. and Gallino, G., 2012. MicroRNA profiling in colorectal cancer highlights miR-1 involvement in MET-dependent proliferation. Mol. Cancer Res., 10: 504-515. https://doi.org/10.1158/1541-7786.MCR-11-0342

Rihani, A., Van Maerken, T., Pattyn, F., Van Peer, G., Beckers, A., De Brouwer, S., Kumps, C., Mets, E., Van der Meulen, J. and Rondou, P., 2013. Effective Alu repeat based RT-Qpcr normalization in cancer cell perturbation experiments. PLoS One, 8: e71776. https://doi.org/10.1371/journal.pone.0071776

Roberts, T.C., Coenen-Stass, A.M. and Wood, M.J., 2014. Assessment of RT-qPCR normalization strategies for accurate quantification of extracellular microRNAs in murine serum. PLoS One, 9: e89237. https://doi.org/10.1371/journal.pone.0089237

Schwarzenbach, H., da Silva, A.M., Calin, G. and Pantel, K., 2015. Data normalization strategies for microRNA quantification. Clin. Chem., 61: 1333-1342. https://doi.org/10.1373/clinchem.2015.239459

Song, J., Bai, Z., Han, W., Zhang, J., Meng, H., Bi, J., Ma, X., Han, S. and Zhang, Z., 2012. Identification of suitable reference genes for qPCR analysis of serum microRNA in gastric cancer patients. Dig. Dis. Sci., 57: 897-904. https://doi.org/10.1007/s10620-011-1981-7

St-Pierre, J., Grégoire, J.-C. and Vaillancourt, C., 2017. A simple method to assess group difference in RT-qPCR reference gene selection using GeNorm: The case of the placental sex. Scient. Rep., 7: 16923. https://doi.org/10.1038/s41598-017-16916-y

Swellam, M.and El-Khazragy, N., 2016. Clinical impact of circulating microRNAs as blood-based marker in childhood acute lymphoblastic leukemia. Tumor Biol., 37: 10571-10576. https://doi.org/10.1007/s13277-016-4948-7

Wotschofsky, Z., Meyer, H.-A., Jung, M., Fendler, A., Wagner, I., Stephan, C., Busch, J., Erbersdobler, A., Disch, A.C. and Mollenkopf, H.J., 2011. Reference genes for the relative quantification of microRNAs in renal cell carcinomas and their metastases. Anal. Biochem., 417: 233-241. https://doi.org/10.1016/j.ab.2011.06.009

Wu, Q., Lu, Z., Li, H., Lu, J., Guo, L. and Ge, Q., 2011. Next-generation sequencing of microRNAs for breast cancer detection. J. BioMed. Biotechnol., 2011: 597145. https://doi.org/10.1155/2011/597145

Xiang, M., Zeng, Y., Yang, R., Xu, H., Chen, Z., Zhong, J., Xie, H., Xu, Y. and Zeng, X., 2014. U6 is not a suitable endogenous control for the quantification of circulating microRNAs. Biochem. biophys. Res. Commun., 454: 210-214. https://doi.org/10.1016/j.bbrc.2014.10.064

Yang, H., Schmuke, J.J., Flagg, L.M., Roberts, J.K., Allen, E.M., Ivashuta, S., Gilbertson, L.A., Armstrong, T.A. and Christian, A.T., 2009. A novel real-time polymerase chain reaction method for high throughput quantification of small regulatory RNAs. Pl. Biotechnol. J., 7: 621-630. https://doi.org/10.1111/j.1467-7652.2009.00429.x

Zhao, H., Ma, T.F., Lin, J., Liu, L.L., Sun, W.J., Guo, L.X., Wang, S.Q., Otecko, N.O. and Zhang, Y.P., 2018. Identification of valid reference genes for mRNA and microRNA normalisation in prostate cancer cell lines. Scient. Rep., 8: 1949. https://doi.org/10.1038/s41598-018-19458-z

Zhu, Y.D., Wang, L., Sun, C., Fan, L., Zhu, D.X., Fang, C., Wang, Y.H., Zou, Z.J., Zhang, S.J. and Li, J.Y., 2012. Distinctive microRNA signature is associated with the diagnosis and prognosis of acute leukemia. Med. Oncol., 29: 2323-2331. https://doi.org/10.1007/s12032-011-0140-5