- Research article
- Open Access
Identification of a set of endogenous reference genes for miRNA expression studies in Parkinson’s disease blood samples
BMC Research Notes volume 7, Article number: 715 (2014)
Research on microRNAs (miRNAs) is becoming an increasingly attractive field, as these small RNA molecules are involved in several physiological functions and diseases. To date, only few studies have assessed the expression of blood miRNAs related to Parkinson’s disease (PD) using microarray and quantitative real-time PCR (qRT-PCR). Measuring miRNA expression involves normalization of qRT-PCR data using endogenous reference genes for calibration, but their choice remains a delicate problem with serious impact on the resulting expression levels. The aim of the present study was to evaluate the suitability of a set of commonly used small RNAs as normalizers and to identify which of these miRNAs might be considered reliable reference genes in qRT-PCR expression analyses on PD blood samples.
Commonly used reference genes snoRNA RNU24, snRNA RNU6B, snoRNA Z30 and miR-103a-3p were selected from the literature. We then analyzed the effect of using these genes as reference, alone or in any possible combination, on the measured expression levels of the target genes miR-30b-5p and miR-29a-3p, which have been previously reported to be deregulated in PD blood samples.
We identified RNU24 and Z30 as a reliable and stable pair of reference genes in PD blood samples.
MicroRNAs (miRNAs) are small non-coding RNAs of 20–22 nucleotides involved in transcriptional and post-transcriptional regulation of gene expression. MiRNAs function via base-pairing with complementary sequences within target mRNA molecules, usually resulting in gene silencing via translational repression or target degradation . MiRNAs are involved in several physiological functions, such as cell cycle, apoptosis, proliferation, differentiation and development . In recent years the research on miRNAs has increased in intensity because of their involvement not only in physiological processes but also in different diseases. Several studies have shown abnormal expression levels of miRNAs in different pathologies, including Parkinson’s disease . The main clinical hallmarks of this neurodegenerative disorder are resting tremor, muscular rigidity, bradykinesia, impaired balance and coordination. Pathologically PD is characterized by the loss of the dopaminergic neurons in the motor control region substantia nigra and an accumulation of protein-filled structures called Lewy bodies . Lewy bodies and other brain specific diagnostic signs cannot be observed until after death of the PD affected individuals. Biomarkers from easy accessible resources, such as peripheral blood, plasma, serum, urine, saliva, and eventually cerebrospinal fluid could be used to detect and monitor the disease much earlier, even before symptoms appear. Unfortunately to date no proven biomarkers are available for the diagnosis of PD . MiRNAs are strong and specific gene regulators and therefore promising candidates to be diagnostic markers and therapeutic targets. Peripheral blood can be considered as a potential diagnostic tool because it is readily obtainable and reflects dynamically a system-wide biology .
To date, only few studies have assessed the expression of blood miRNAs related to Parkinson’s disease using microarray and quantitative real-time PCR [7, 8]. Several methods have been developed to measure miRNA expression . Of these methods qRT-PCR is superior due to its sensitivity, specificity and linear dynamic range of quantification [10, 11] and is therefore generally accepted as the gold standard for accurate gene quantification and small RNA profiling. To distinguish true biological changes from technical variations, several variables must be considered and controlled for gene-expression analysis, such as the amount of the starting material and differences between tissues in overall transcriptional activity. The use of multiple stable reference genes is generally necessary as the method of choice for qRT-PCR data normalization. The normalization procedure is adapted to compare the expression levels of target genes with the expression of a reference gene to obtain a relative quantification of the investigated target gene . The expression of an ideal reference gene should not vary in the tissues or cells under investigation and in response to treatment or presence of a pathology . Accurate expression profiling is critically dependent on technical variations due to the applied normalization strategies. A misleading estimation of the expression analysis data can arise if only a single or improper reference genes are used . These stable reference genes can be identified from a set of candidate reference genes in a pilot experiment on a selection of samples that are representative for the experimental conditions under investigation. Different algorithms, such as geNorm and Normfinder allow the ranking of candidate reference genes according to their stability and indicate the optimal number and combination of reference genes required for accurate normalization of gene expression . In the case of miRNA expression profiling, only few candidate reference miRNAs, including miRNA 103, have been identified [15, 16]. Typically, other small endogenous noncoding RNAs such as small nuclear (snRNAs) and small nucleolar RNAs (snoRNAs) have been used  as they share similar properties, such as RNA stability and size and are expressed abundantly. Unfortunately their expression is not always stable . As no universal reference gene has so far been identified, there is no substitute for empirical validation of normalization that is appropriate to the particular experimental design and goals . Recently expression studies of circulating miRNAs in serum or plasma as potential biomarkers for PD are emerging [18, 19]. In these studies synthetic, non-human spike-in miRNAs, used frequently to monitor RNA purification and retrotranscription efficiencies, are almost exclusively used as well for normalization in expression analyses, because of the lack of established reference genes . Moreover, new emerging technologies including Droplet Digital PCR technology, which provides an absolute quantification of nucleic acids, are trying to overcome the problem of the normalization, but for now they are available only in few laboratories.
The aim of the present study was to evaluate the suitability of a set of commonly used small RNAs as normalizers and to identify which of these miRNAs might be considered reliable reference genes in qRT-PCR expression analyses on PD blood samples. We compared the potential reference genes, snoRNA RNU24 (RNU24), snRNA RNU6B (RNU6B), snoRNA Z30 (Z30) and miR-103a-3p to find the best set of reference genes to use for normalization using (i) comparative delta-Ct , (ii) the NormFinder algorithm  and (iii) qbasePLUS, an improved version of the geNorm algorithm . The effect of different reference genes used on the relative quantification was shown using the target genes miR-30b-5p and miR-29a-3p, previously shown to be differentially expressed in PD blood samples [7, 8].
Patient and control samples collection
38 patients affected by idiopathic Parkinson’s disease (PD) attending the Movement Disorders outpatient clinic of the Bolzano Hospital (Italy) were enrolled. Diagnosis of PD was made according to Gelb PD criteria . All the patients were treated with Levodopa. The subjects with cognitive impairment or unable to sign informed consent or affected by atypical Parkinsonism were excluded. 38 disease-free controls matched for sex and age (range ± 3 years), collected among spouses or unrelated companions of the patients, were also enrolled. The study was approved by the ethics committee of the Bolzano Hospital and all of the enrolled subjects provided written informed consent to participate.
RNA isolation and quality control
Whole peripheral blood from patients and controls was collected per participant in a Na3-Citrate buffered Venosafe® Plastic Tube (VF-054SBCS07, Terumo) and stored at room temperature. Within 6 hours of blood collection the peripheral blood mononuclear cells (PBMCs) were isolated after lysis of red blood cells by Red Blood Cell (RBC) Lysing Buffer. Total RNA, including small RNAs, was extracted from white blood cells using TRIzol® reagent (Invitrogen) according to the manufacturer’s instructions and stored at -80°C.
The quality and quantity of the extracted RNA was assessed by agarose gel electrophoresis and by Experion™ Automated Electrophoresis System (Bio-Rad Laboratories s.r.l., Milano, Italy), using Experion™ StdSens Analysis Kit (Bio-Rad Laboratories s.r.l., Milano, Italy). The RNA Quality Indicator (RQI) feature of Experion software allows estimating the level of the integrity of eukaryotic total RNA samples. All 76 collected samples showed an a RQI ≥ 8.0 and were included in the expression analyses.
Reverse Transcription and gene specific quantitative real time PCR
Expression analyses were performed by single gene specific stem loop RT-PCR. Reverse transcription (RT) reactions were performed on 7 ng of total RNA using the TaqMan miRNA Reverse Transcription Kit (Applied Biosystem; Part.No 4366597) and miRNA-specific stem-loop primers (Table 1) in a scaled down volume of 10 μl RT reaction, according to the manufacturer’s instructions.
The thermal cycling parameters of reverse transcription were 30 min at 16°C, 30 min at 42°C, and 5 min at 85°C. The cDNA samples were diluted in nuclease-free water and stored at -20°C.
All quantitative real-time PCR reactions were performed in triplicate on Bio-Rad 96CFX instrument (Bio-Rad Laboratories s.r.l., Milano, Italy), in scaled down 10 μl reaction volumes using 0.7 μl of RT product per reaction using Universal MasterMixII, no UNG and gene specific TaqMan® small RNA Assays (Applied Biosystem®; Table 1) according to the manufacturer’s instructions. All reactions were performed on hard shell PCR plates (Bio-Rad Laboratories Inc. Part.No HSP9645), sealed using adhesive Microseal ‘B’ Films (Bio-Rad Laboratories Inc. Part.No RSN 102595). One reference sample and one NTC (blank) were included in triplicate on each plate. The thermal cycling parameters were 10 min at 95°C followed by 40 cycles of 15 sec at 95°C and 1 min at 60°C.
PCR amplification efficiencies for all target and reference genes were determined from the slope of the log-linear portion of calibration curves, generated from a two-fold dilution series of one reference sample at five dilution points for three technical replicates for each gene assay (Table 1). PCR reactions were performed as described previously for the quantitative real-time PCR reactions. The relative expression of miR-30b-5p and miR-29a-3p were assessed using the Bio-Rad CFX Manager v1.6 (Bio-Rad Laboratories s.r.l., Milano, Italy) software. For normalization, different combinations of reference genes among RNU24, RNU6B, Z30 and miR-103a-3p were used.
In order to determine the expression stability of the four candidate reference genes, statistical analyses of their expression across all samples were performed using three different algorithms: NormFinder , GeNorm , and comparative delta-Ct method . NormFinder computes a stability measure and the lowest value indicates the most stable gene expression. Samples are grouped to allow direct estimation of expression variation, ranking genes according to the similarity of their expression profiles by using a model-based approach. NormFinder takes into account inter- and intragroup variation for normalization factor calculations avoiding misinterpretation caused by artificial selection of co-regulated genes. GeNorm computes an M value describing the variation of a gene compared to all other candidate genes. Similarly to the NormFinder stability value, lower M values indicate stable gene expression. Specifically, stable reference genes are supposed to have an M value smaller than 1.5.
Comparative delta-Ct compares relative expression of pairs of genes within each sample set and ranks stability of reference genes according to the repeatability of the gene expression difference.
Scatter plots and box plots were used to visualize the presence of outliers in the expression Ct values; outlier values, as evidenced by the box plot, were labelled for possible subsequent removal. Medians and interquartile ranges were produced to summarize data. The medians were calculated on the differences of the expression between each case and its matched control.
For miR-29a-3p and miR-30b-5p, the distribution of the difference of the expression levels between cases and controls, within each matched pair and normalized with different combinations of reference genes, was checked using a Skewness-Kurtosis test for normality. A subset of the ladder of powers for variable transformation in case of departure from normality was evaluated .
Relative expression of miR-29a-3p and miR-30b-5p was analysed using a Wilcoxon matched-pairs signed-ranks test. Although 15 tests for each microRNA were performed, no adjustment for multiple comparisons was done. This is because in practice the only test that would have been performed is the one based on data a priori normalized for the best reference gene combination. Spearman’s rank correlation coefficients were computed on Ct reference gene data and on normalized miR-29a-3p and miR-30b-5p data. GeNorm analysis was performed using Biogazelle’s qbasePLUS software (Bio-Rad Laboratories Inc.) and NormFinder analysis was performed using “NormFinder.xla”, a Microsoft Excel-based Visual Basic application. All other analyses were performed using Stata 12 .
The selection of the best set of endogenous reference genes for gene expression studies in Parkinson’s disease blood samples was based on the efficiency of the TaqMan® MicroRNA Assays, the quality of the related expression data, and on the expression stability analysis.
Evaluation of the amplification efficiencies of the TaqMan® MicroRNA Assays
The hallmarks for an accurate and optimized qRT-PCR assay are a linear correlation coefficient (r2) equal or greater than 0.98 and a PCR amplification efficiency from 90% to 110% . The r2 value indicates the quality of the fit of the standard curve to the plotted data points. Primer efficiency indicates the amplicon doubling rate of a specific primer pair during a PCR. An efficiency of 100% indicates that the cDNA target is duplicated at every PCR cycle during the exponential phase. The efficiencies of the TaqMan® MicroRNA Assays of the four candidate reference genes (RNU24, RNU6B, Z30 and miR-103a-3p), and the two target genes (miR-30b-5p and miR-29a-3p) were calculated from the slope of the log-linear portion of calibration curves as described (Materials and Methods; Reverse Transcription and quantitative real time PCR) and are reported in Table 1. RNU24, Z30, miR-103a-3p, miR-30b-5p, and miR-29a-3p showed amplification efficiencies and r2 values ranging respectively between 90.4% and 97.5% and between 0.92 and 0.99 demonstrating a high performance. In contrast, RNU6B showed a low amplification efficiency of 79.9% with a correlation coefficient of 0.967.
Data quality control
The analysis of the generated Ct data revealed the presence of a few outliers, as evidenced in Figure 1. Because no evident technical reason could be found to exclude them, given the high quality of the extracted RNA as evaluated by Experion RNA chip electrophoresis, and the overall high performance of the qRT-PCR reactions, they were not removed from the analyses, but they were rather considered an expression of normal biological variability. The differences of the relative gene expression in cases and in controls, within each matched pair, did not always follow a normal distribution. Since no unique transformation was able to restore deviations from normality, a non-parametric Wilcoxon matched-pairs signed-ranks test was subsequently used to analyse data.
Expression stability analysis
The results of the expression stability assessment of Z30, RNU24, RNU6B and miR-103a-3p genes among all analysed samples are reported in Tables 2 and 3, showing the results of the comparative delta-Ct method, and the comparison with the NormFinder and geNorm results. Independently from the method used, the gene ranking remained unchanged, indicating Z30 and RNU24 as the best reference genes. In particular, the combination of Z30 and RNU24 was selected as the best among all the others, as underlined by the stability value of 0.011 calculated by NormFinder, the M value of 0.489 calculated by geNorm (see Additional file 1: Figure S1), and the delta-Ct values (Table 2). This combination showed the smallest and therefore best stability values among all the other combinations of reference genes (see Additional file 1: Figure S1) whereas the worst combination was attributed to RNU6B and miR-103a-3p.
Unexpectedly, miR-103a-3p, previously indicated as a suitable reference gene for quantification of microRNAs in peripheral blood of PD patients through qRT-PCR , was consistently found to be the worst reference gene in the analysed sample (Tables 2 and 3).
Furthermore, since the expression of reference genes is supposed to remain constant among different samples, we would expect that the relationship between the expression of two reference genes is linear.
To estimate the inter-gene relations of the four putative reference genes, correlation analysis was performed. Results indicate that the strength of the linear relationship is good for Z30 and RNU24 expression (r = 0.8004, p < 0.0001), which are those candidates ranked as the best reference genes. The correlation of the expression of RNU6B and miR-103a-3p with Z30 and RNU24 is smaller, with the lowest correlations observed between the expression levels of miR-103a-3p and the remaining reference genes (see Additional file 1: Table S1).
The same trend is also indicated by the scatter plot matrix (see Additional file 1: Figure S2), where the most linear relationship could be observed between Z30 and RNU24.
All the analyses were also repeated excluding outlier samples (sensitivity analyses), but no difference in the rank of the reference genes was observed (data not shown).
Consequence of normalizing expression of miR-29a-3p and miR-30b-5p in PD cases and matched controls with different combinations of reference genes
We tested the effect of using different combinations of reference genes on the relative expression values using miR-29a-3p and miR-30b-5p as targets, which were found to be deregulated in PD in previous studies [7, 8].
The results of the Wilcoxon matched-pairs signed-ranks test (Table 4) show that there is an increased relative expression of miR-29a-3p and miR-30b-5p in patients with PD in comparison to controls.
Normalizing the expression values of miR-29a-3p and miR-30b-5p for the best reference gene combination of Z30 and RNU24 results in a median difference of the relative expression of 0.68 (IQR = -0.17; 1.74) for miR-29a-3p (p value = 0.0032), and equal to 1.74 (IQR = 0.05; 4.26) for miR-30b-5p (p value =0.0009). This trend was consistently observed using Z30, RNU24 and RNU6B as reference genes, either alone or in different combinations (Table 4).
However, the normalization of the microRNA expression using miR-103a-3p as reference gene results in an opposite effect, with a lower expression of miR-29a-3p and miR-30b-5p in cases, in comparison with controls. The usage of miR-103a-3p alone as a reference gene was sufficient to completely reverse the direction of the expression difference between cases and controls of the investigated miRNAs (miR-29a-3p median difference = -0.58, IQR = -1.33; -0.05, p = 0.0003 and miR-30b-5p median difference = -1.16, IQR = -2.53; 0.23, p = 0.0005).
The inclusion of miR-103a-3p in any other combination of reference genes systematically worsened the observed p values. This effect was clearly visible with the addition of miR-103a-3p to the combination of Z30 and RNU24, which had a strong impact on the results observed for both miR-29a-3p and miR-30b-5p, with the p value of miR-29a-3p changing from p = 0.0032 to p = 0.5473 and the p value of miR-30b-5p changing from p = 0.0009 to p = 0.3422. Indeed, the relative expression levels of miR-29a-3p/miR-30b-5p were not always well correlated with each other. But some differences could be observed, depending on the reference gene set used for normalization, as shown in Additional file 1: Tables S2 and S3. Smaller correlation coefficients could be observed when miR-103a-3p alone was used as a reference gene, while correlations were improved when it was present in other combinations.
Sensitivity analyses were finally repeated excluding outlier samples, but no difference in the results was observed (data not shown).
Quantitative real-time PCR is the most commonly available and reliable method used in expression profile studies aimed at detecting changes in gene expression, because alternative approaches do not yet match the sensitivity and specificity of PCR-based approach. The choice of housekeeping genes in a qRT-PCR study is a critical step before starting any experiment, in order to avoid inaccurate interpretation of the data that may lead to biased results. Indeed, it has been shown in several publications that traditional reference genes used in qRT-PCR studies do not always show a stable expression pattern. The same gene revealed as almost invariant for certain tissues or cell types, could present highly variable expression levels in other tissues or experimental conditions [20, 25, 26]. Thus, suitable control genes are extremely specific for particular sample sets and experimental models, being a crucial component in assessing gene expression patterns with confidence. Striking examples are provided by the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and beta-actin genes, which have been extensively used for normalization of gene expression data in a broad range of different tissues and pathologies, including the nervous system [25–28]. It has been shown however, that GAPDH and beta-actin genes are indeed direct targets of miR-644a in prostate cancer cell lines, demonstrating the unsuitability of GAPDH and beta-actin as internal controls in miR-644a functional studies and emphasizing the need to carefully consider the choice of a reference gene according to the specific study design .
Small nuclear RNAs as well as miRNAs and snoRNAs are widely used as reference genes for miRNA studies in tissue samples and blood. To date, in the field of the Parkinson's disease there is no work reporting the validation of a set of reliable reference genes for a miRNA profiling study in blood PD samples.
For that reason, the present study aims not only to identify and to validate specific endogenous reference genes as reported in the literature, but also to evaluate the impact of the use of different combinations of these reference genes on the interpretation of expression profile differences between cases and controls of certain target miRNAs.
We selected from the literature the most commonly used reference genes in qRT-PCR blood studies or in neurological diseases studies. We then assessed their suitability as normalizers, and their effect, alone or in combination, on the target genes miR-30b-5p and miR-29a-3p in PD blood samples, previously shown to have altered expression profiles in PD [7, 8]. The selection of the best set of endogenous reference genes for gene expression studies in Parkinson’s disease blood samples from those tested was based on the efficiency of the TaqMan® MicroRNA Assays, the quality of the related expression data, on the expression stability analysis, and biological data available in the literature.
Using these metrics of efficiency, quality and stability, RNU24 ranks top of the list, followed by Z30 as the best reference genes in our study. Moreover, neither RNU24 nor Z30 have been related so far to neurodegeneration and/or Parkinson’s disease. On the contrary, miR-103a-3p is ranked as the worst reference gene, so much so that in combination with other reference genes it is able to bias results.
It is important to underline that miR-103a-3p alone or in combination with the other reference genes reverses the direction of the expression of miR-29a-3p and miR-30b-5p. Our results are consistent with the available literature evidence, which support the unsuitability of miR-103a-3p as a reference gene in PD studies for a number of reasons. MiR-103a-3p and miR-107 have been recently shown to target human CDK5R1, coding for the regulatory subunit 1 of cyclin-dependent kinase 5, a protein found to be hyperactivated in different neurodegenerative diseases, including Parkinson’s and Alzheimer's diseases [31, 32]. The same two miRNAs have also been shown to regulate insulin sensitivity by targeting caveolin-1 . The connection between miR-103a-3p and insulin is relevant, since insulin is known to regulate dopamine release. Moreover, in a recent study a correlation between PD and a hyperglycaemic status was shown in rats . Last, but not of minor importance when studying neurodegenerative disease, is the evidence that the expression of miR-103a-3p varies depending on the age of the patients . In our study, the paired case–control experimental design allows correction for the confounding effect of age, but adjustment for age would be required for unmatched designs. This may be unfeasible for small qRT-PCR studies, therefore the confounding effect would not be controlled and the results would be biased.
Contrary to miR-103a-3p, the insertion of RNU6B in combinations with other reference genes has a lower impact on the stability values, but we hypothesize that this effect can be attributed to the high stability of RNU24 and Z30. RNU6B was not considered to be a reliable reference gene for PD blood samples in our study, because the efficiency, the r2 and the stability values were too low. Furthermore, in support of this conclusion, RNU6B has been shown to be stable under some but not all conditions of neuronal differentiation .
In conclusion, the results of this study indicate that relative expression differences of the selected microRNAs between Parkinson’s cases and controls were profoundly affected by the choice of the reference gene set, thus demonstrating the complexity of the choice of reference gene set. A real endogenous reference gene should be stably expressed, regardless of physiological conditions, including age and disease. In our hands, using miR-103a-3p as a reference gene, alone or in combination with other suitable genes, produced contrasting results within our dataset. Our data suggest that RNU24 and Z30 constitute a reference gene set that will lead to reliable data for normalization of miRNA expression profiles in blood of PD patients versus controls, at least until new technologies for direct PCR quantification such as droplet-digital PCR become more widespread and affordable. The results confirm the importance of a careful choice of the reference genes, which should be evaluated depending on the specific on-going study, not relying solely on evidence already published.
Reverse transcription and quantitative real time polymerase chain reaction
Peripheral blood mononucleated cells
RNA quality indicator
Saini HK, Griffiths-Jones S, Enright AJ: Genomic analysis of human microRNA transcripts. Proc Natl Acad Sci U S A. 2007, 104 (45): 17719-17724. 10.1073/pnas.0703890104.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.
Heyer MP, Pani AK, Smeyne RJ, Kenny PJ, Feng G: Normal midbrain dopaminergic neuron development and function in miR-133b mutant mice. J Neurosci. 2012, 32 (32): 10887-10894. 10.1523/JNEUROSCI.1732-12.2012.
Bras J, Singleton A, Cookson MR, Hardy J: Emerging pathways in genetic Parkinson’s disease: Potential role of ceramide metabolism in Lewy body disease. FEBS J. 2008, 275 (23): 5767-5773. 10.1111/j.1742-4658.2008.06709.x.
Shtilbans A, Henchcliffe C: Biomarkers in Parkinson’s disease: an update. Curr Opin Neurol. 2012, 25 (4): 460-465. 10.1097/WCO.0b013e3283550c0d.
Liew CC, Ma J, Tang HC, Zheng R, Dempsey AA: The peripheral blood transcriptome dynamically reflects system wide biology: a potential diagnostic tool. J Lab Clin Med. 2006, 147 (3): 126-132. 10.1016/j.lab.2005.10.005.
Margis R, Margis R, Rieder CR: Identification of blood microRNAs associated to Parkinsonis disease. J Biotechnol. 2011, 152 (3): 96-101. 10.1016/j.jbiotec.2011.01.023.
Martins M, Rosa A, Guedes LC, Fonseca BV, Gotovac K, Violante S, Mestre T, Coelho M, Rosa MM, Martin ER, Vance JM, Outeiro TF, Wang L, Borovecki F, Ferreira JJ, Oliveira SA: Convergence of miRNA expression profiling, alpha-synuclein interacton and GWAS in Parkinson’s disease. PLoS One. 2011, 6 (10): e25443-10.1371/journal.pone.0025443.
Kang K, Peng X, Luo J, Gou D: Identification of circulating miRNA biomarkers based on global quantitative real-time PCR profiling. J Anim Sci Biotechnol. 2012, 3 (1): 1-9. 10.1186/2049-1891-3-1.
D’haene B, Mestdagh P, Hellemans J, Vandesompele J: miRNA expression profiling: from reference genes to global mean normalization. Methods Mol Biol. 2012, 822: 261-272. 10.1007/978-1-61779-427-8_18.
Hellemans J, Vandesompele J: qPCR data analysis - unlocking the secret to successful results. PCR Troubleshooting and Optimization: The Essential Guide. Edited by: Suzanne K, Nick O. 2011, Poole (UK): MO BIO Laboratories. Caister Academic Press, 1-13.
Wotschofsky Z, Meyer HA, Jung M, Fendler A, Wagner I, Stephan C, Busch J, Erbersdobler A, Disch AC, Mollenkopf HJ, Jung K: Reference genes for the relative quantification of microRNAs in renal cell carcinomas and their metastases. Anal Biochem. 2011, 417 (2): 233-241. 10.1016/j.ab.2011.06.009.
Sanders I, Holdenrieder S, Walgenbach-Brunagel G, von Ruecker A, Kristiansen G, Muller SC, Ellinger J: Evaluation of reference genes for the analysis of serum miRNA in patients with prostate cancer, bladder cancer and renal cell carcinoma. Int J Urol. 2012, 19 (11): 1017-1025. 10.1111/j.1442-2042.2012.03082.x.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3 (7): 1-12.
Peltier HJ, Latham GJ: Normalization of microRNA expression levels in quantitative RT-PCR assays: identification of suitable reference RNA targets in normal and cancerous human solid tissues. RNA. 2008, 14 (5): 844-852. 10.1261/rna.939908.
Latham G: Normalization of MicroRNA Quantitative RT-PCR Data in Reduced Scale Experimental Designs. MicroRNAs and the Immune System: Methods and Protocols, Methods in Molecular Biology. Volume 667. Edited by: Springer Science+Business Media, Monticelli S. 2010, New York: Humana Press, 19-31.
Gee HE, Buffa FM, Camps C, Ramachandran A, Leek R, Taylor M, Patil M, Sheldon H, Betts G, Homer J, West C, Ragoussis J, Harris AL: The small-nucleolar RNAs commonly used for microRNA normalisation correlate with tumour pathology and prognosis. Br J Cancer. 2011, 104 (7): 1168-1177. 10.1038/sj.bjc.6606076.
Khoo SK, Petillo D, Kang UJ, Resau JH, Berryhill B, Linder J, Forsgren L, Neuman LA, Tan AC: Plasma-Based Circulating MicroRNA Biomarkers for Parkinson’s Disease. J Parkinson’s Dis. 2012, 2: 321-331.
Cardo LF, Coto E, de Mena L, Ribacoba R, Moris G, Menendez M, Alvarez V: Profile of microRNAs in the plasma of Parkinson’s disease patients and healthy controls. J Neurol. 2013, 260 (5): 1420-1422. 10.1007/s00415-013-6900-8.
Silver N, Best S, Jiang J, Thein SL: Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol Biol. 2006, 7: 33-10.1186/1471-2199-7-33.
Andersen CL, Jensen JL, Orntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64 (15): 5245-5250. 10.1158/0008-5472.CAN-04-0496.
Gelb DJ, Oliver E, Gilman S: Diagnostic criteria for Parkinson disease. Arch Neurol. 1999, 56 (1): 33-39. 10.1001/archneur.56.1.33.
Tukey JW: Exploratory Data Analysis. 1977, Boston: Anonymous Addison-Wesley Publishing Company, Reading, MA
StataCorp: Stata Statistical Software: Release 12. 2011, College Station, TX: StataCorp LP
Dheda K, Huggett JF, Bustin SA, Johnson MA, Rook G, Zumla A: Validation of housekeeping genes for normalizing RNA expression in real-time PCR. Biotechniques. 2004, 37 (1): 112-–4, 116, 118–9
Valente V, Teixeira SA, Neder L, Okamoto OK, Oba-Shinjo SM, Marie SK, Scrideli CA, Paco-Larson ML, Carlotti CG: Selection of suitable housekeeping genes for expression analysis in glioblastoma using quantitative RT-PCR. BMC Mol Biol. 2009, 10: 17-2199. 10.1186/1471-2199-10-17. 10-17
Benn CL, Fox H, Bates GP: Optimisation of region-specific reference gene selection and relative gene expression analysis methods for pre-clinical trials of Huntington’s disease. Mol Neurodegener. 2008, 3: 1-17. 10.1186/1750-1326-3-1.
Glare EM, Divjak M, Bailey MJ, Walters EH: beta-Actin and GAPDH housekeeping gene expression in asthmatic airways is variable and not suitable for normalising mRNA levels. Thorax. 2002, 57 (9): 765-770. 10.1136/thorax.57.9.765.
Sikand K, Singh J, Ebron JS, Shukla GC: Housekeeping gene selection advisory: glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and beta-actin are targets of miR-644a. PLoS One. 2012, 7 (10): e47510-10.1371/journal.pone.0047510.
Moncini S, Salvi A, Zuccotti P, Viero G, Quattrone A, Barlati S, De Petro G, Venturin M, Riva P: The role of miR-103 and miR-107 in regulation of CDK5R1 expression and in cellular migration. PLoS One. 2011, 6 (5): e20038-10.1371/journal.pone.0020038.
Smith PD, Crocker SJ, Jackson-Lewis V, Jordan-Sciutto KL, Hayley S, Mount MP, O’Hare MJ, Callaghan S, Slack RS, Przedborski S, Anisman H, Park DS: Cyclin-dependent kinase 5 is a mediator of dopaminergic neuron loss in a mouse model of Parkinson’s disease. Proc Natl Acad Sci U S A. 2003, 100 (23): 13650-13655. 10.1073/pnas.2232515100.
Wang WX, Rajeev BW, Stromberg AJ, Ren N, Tang G, Huang Q, Rigoutsos I, Nelson PT: The expression of microRNA miR-107 decreases early in Alzheimer’s disease and may accelerate disease progression through regulation of beta-site amyloid precursor protein-cleaving enzyme 1. J Neurosci. 2008, 28 (5): 1213-1223. 10.1523/JNEUROSCI.5065-07.2008.
Trajkovski M, Hausser J, Soutschek J, Bhat B, Akin A, Zavolan M, Heim MH, Stoffel M: MicroRNAs 103 and 107 regulate insulin sensitivity. Nature. 2011, 474 (7353): 649-653. 10.1038/nature10112.
Herrera BM, Lockstone HE, Taylor JM, Ria M, Barrett A, Collins S, Kaisaki P, Argoud K, Fernandez C, Travers ME, Grew JP, Randall JC, Gloyn AL, Gauguier D, McCarthy MI, Lindgren CM: Global microRNA expression profiles in insulin target tissues in a spontaneous rat model of type 2 diabetes. Diabetologia. 2010, 53 (6): 1099-1109. 10.1007/s00125-010-1667-2.
Noren Hooten N, Abdelmohsen K, Gorospe M, Ejiogu N, Zonderman AB, Evans MK: microRNA expression patterns reveal differential expression of target genes with age. PLoS One. 2010, 5 (5): e10724-10.1371/journal.pone.0010724.
Lim QE, Zhou L, Ho YK, Wan G, Too HP: snoU6 and 5S RNAs are not reliable miRNA reference genes in neuronal differentiation. Neuroscience. 2011, 199: 32-43.
The authors are grateful to the study participants for their participation and collaboration in this research project. We thank Diana Riekschnitz and Clemens Egger for technical support. The authors acknowledge the help of Giulia Giannini, Agatha Eisendle, Edith Kompatscher, and Monika Mair in recruiting the study participants. The authors are grateful to Deborah Mascalzoni, Fabiola del Greco, Pietro Cortelli and Marika Falla for their support in ethical, statistical and medical issues. This work was supported by the Department for Promotion of Educational Policies, Universities and Research of the Autonomous Province of Bolzano, South Tyrol and the South Tyrolean Sparkasse Foundation.
The authors declare that they have no competing interests.
AS designed and performed the experiments and wrote the manuscript. LF performed the statistical analyses and wrote the manuscript. HB contributed to the expression analysis and revised the manuscript. SZ designed and contributed to the sample collection and clinical characterization of the patients. AZ contributed to the sample collection and experimental work. AP contributed to the expression stability analysis and revised the manuscript. PPP designed and contributed to the sample collection and clinical characterization of the patients. AAH critically revised the manuscript and contributed to the design of the sample collection. CS designed and supervised the study, contributed to the design of the sample collection, and wrote the manuscript. All authors read and approved the final manuscript.
Alice Serafin, Luisa Foco contributed equally to this work.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Serafin, A., Foco, L., Blankenburg, H. et al. Identification of a set of endogenous reference genes for miRNA expression studies in Parkinson’s disease blood samples. BMC Res Notes 7, 715 (2014). https://doi.org/10.1186/1756-0500-7-715
- geNorm algorithm
- Normfinder algorithm
- Comparative delta-Ct