Expressed Alu repeats as a novel, reliable tool for normalization of real-time quantitative RT-PCR data
- Manuela Marullo†1,
- Chiara Zuccato†1,
- Caterina Mariotti2,
- Nayana Lahiri3,
- Sarah J Tabrizi3,
- Stefano Di Donato2 and
- Elena Cattaneo1Email author
© Marullo et al.; licensee BioMed Central Ltd. 2010
Received: 28 July 2009
Accepted: 28 January 2010
Published: 28 January 2010
We describe a novel strategy for mRNA normalization in quantitative real-time PCR that is based on expressed Alu repeat amplification as a measure for the mRNA fraction. We show that expressed Alu repeat amplification is a fast, accurate normalization tool that can be successfully used for quantification of selected mRNA in the human transcriptome. This result is particularly important for clinical diagnosis and biomarker validation studies based on mRNA detection in human blood.
Several factors hinder the reliable quantification of mRNA in gene expression studies performed in human samples, including sample availability, handling, and storage and the strategy used to normalize selected mRNA over the total mRNA fraction. With real-time quantitative RT-PCR (RT-qPCR), quantification errors depend on variations in the amount of total RNA or cDNA between samples, which remain difficult to determine. Various strategies have been applied in efforts to normalize these variations. Under controlled conditions of reproducible extraction of good-quality RNA, one possible strategy relies on total RNA quantification (RNA mass quantity) and/or the use of synthetic internal standards. To date, the amplification of a reference gene as the internal standard is the most frequently used approach to normalizing the mRNA fraction. This strategy is best applied when a large number of target genes are to be screened and makes inter-laboratory standardization easier than using total RNA content, which might depend on the quantification method. However, the reliable quantification of selected targets requires identification of a proper internal control gene for normalization, one that exhibits stable expression in the given cell or tissue under investigation.
A first accurate strategy for normalization was based on the GeNorm algorithm, which allows identification of the most stably expressed control genes in a human tissue of interest. The expression of ten common reference genes is analyzed in the tissue investigated by RT-qPCR, and the GeNorm algorithm is applied to calculate the gene expression stability (M) of the different control genes. The level of the target gene is normalized over the normalization factor (NF) obtained from the geometric mean of expression value of three or more stable reference genes .
Although the GeNorm algorithm is a robust normalization method, the reference gene validation requires extensive experimental work and a high quantity of sample material. These requirements may represent a problem in the context of gene expression studies involving the human transcriptome for which a limited amount of human tissue samples are available. Moreover, the choice of stable control genes becomes particularly difficult and is an expensive and lengthy process in pathological conditions characterized by transcriptional dysregulation. Altered activity of transcription factors and DNA target sequences may also affect the stability of mRNA expression levels of the selected reference gene, thus producing an erroneous quantification of the selected mRNA. This possibility leads to the requirement of selecting many genes and extensive experimental work before choosing a panel of stable reference genes.
To avoid these problems, we propose a new strategy for mRNA normalization in RT-qPCR that is based on expressed Alu repeat (EAR) amplification as a measure for the total mRNA fraction. More than one million copies of the approximately 300-bp Alu element are interspersed throughout the human genome, with up to 75% of all known genes containing Aluinsertions within their introns and/or untranslated regions (UTRs) . Therefore, the differential expression of a number of genes in the tissues or cells under investigation will not influence EAR abundance in the transcriptome.
In this study, we set up a standardized RT-qPCR protocol for EAR amplification in the human transcriptome. We provide evidence that normalization based on EAR amplification is a suitable, fast, and precise tool for quantification of selected mRNA by RT-qPCR in human biological samples.
There is a long-standing interest in brain-derived neurotrophic factor (BDNF) because of its implications in neurodegenerative diseases . BDNF is a neurotrophin that is crucial for the survival, maintenance, and differentiation of subpopulations of neurons in the central and peripheral nervous systems. Alterations in BDNF brain levels and activity have been described in various neurodegenerative disorders, most notably Huntington's disease (HD) . Although BDNF is highly concentrated in the nervous system, it can also be detected in human blood and the blood of other mammals [4, 5]. Because of that, several efforts have targeted testing BDNF levels in blood with the aim of assessing its potential value as a disease biomarker. For HD and other diseases, such as depression, schizophrenia, and Alzheimer's disease, its reliable detection in peripheral tissues may be crucial for understanding biological and pathogenic processes in humans and for assessing the potential activity of candidate therapeutic agents. Unfortunately, in human blood, this measurement has always proved problematic, especially for the BDNF protein. In fact, many factors influence this measurement, including the way the blood is collected, its storage, and the methodology employed to prepare the samples .
Set up of a new normalization strategy: EAR amplification
Comparison of reference gene stability in human blood versus EARs
Primers sequences for GeNorm
NCBI accession number
Hypoxanthine purine synthesis phosphoribosyl-transferase 1
Tyrosine 3-monooxygenase/S tryptophan 5-monooxygenase activation protein, zeta polypeptide
Guanine nucleotide binding protein, beta polypeptide 2-like 1
Glyceraldehyde-3-oxidoreductase phosphate dehydrogenase
To assess the validity of EARs as a normalization tool versus the GeNorm strategy, we evaluated BDNF mRNA levels normalized to EARs and to NF in blood from ten healthy volunteers. Figure 3b shows that the quantification of BDNF mRNA in human blood is comparable after normalization with the two different strategies. The median value shown in Figure 3c confirms this result: the median values of blood BDNF mRNA levels normalized versus EARs and versus NF are similar to each other (0.95 versus 1.1), while the median value of BDNF mRNA normalized versus GAPDH mRNA is much lower (0.65). These data indicate that EARs represent a good tool for normalization and that normalization over a single reference gene provides incorrect results.
Subsequently, we applied the EAR normalization tool for the quantification of BDNF mRNA in the eight control participants analyzed in Figure 1, in which the messenger level of the neurotrophin is represented relative to GAPDH or β2M mRNA content. Figure 3d shows the median value of BDNF mRNA normalized over GAPDH (1), β2M (1.6), or EARs (0.55). As shown, the distribution of the BDNF mRNA values normalized to EARs is tighter with respect to the distribution of BDNF mRNA normalized to GAPDH or β2M mRNA, indicating that the EAR normalization strategy removes most of the non-specific variations arising from the use of an inappropriate internal reference gene.
The experiments described above prove the reliability of the EAR normalization strategy with respect to GeNorm as well as the efficiency of this new method with respect to normalization to a single housekeeping gene to measure BDNF mRNA levels in human blood.
EAR normalization applies to the quantification of mRNAs specifically expressed in blood
To evaluate the reliability of the EAR normalization strategy with other, more blood-specific genes, we analyzed IL8 (interleukin-8) mRNA levels in blood from control participants. Validation of EAR normalization for cytokine quantification in blood might be particularly useful for measuring the balance between pro- and anti-inflammatory cytokines in chronic inflammatory disorders or autoimmune disease. Cytokine mRNA quantification is therefore of clinical relevance to immuno-monitoring and for the evaluation of the disease status in patients.
We concluded that EARs are a suitable normalization tool for gene expression study of cytokines in human blood.
EAR normalization for human tissues other than blood
Fast, sensitive, and accurate quantification in RT-qPCR analysis is a crucial tool not only for assessing mRNA in human blood but also for evaluating the levels of given molecules in tissues such as post-mortem human brain. In a previous study, we analyzed BDNF mRNA levels in 38 post-mortem human cortices from control participants and HD patients using GAPDH as the reference gene . The GAPDH gene is, in fact, one of the more commonly used reference genes in such analyses .
The data described above show that, in addition to peripheral blood, the EAR normalization strategy can be applied to samples from human brain and skeletal muscle, implying that this method could be used for assessing different human samples.
Accurate normalization of mRNA levels is an absolute prerequisite for obtaining reliable results in gene expression studies. Recently, considerable attention has focused on the systematic study of normalization procedures. Up to the year 2000, a number of articles discussed the problem of normalizing mRNA values over a reference gene, indicating that the most commonly used reference genes, such as GAPDH and βACT, exhibit variation under different experimental conditions [9, 11, 12]. In addition, finding suitable reference genes is not straightforward, and different internal control genes should be tested in the context of a specific experimental design. Different groups have faced this problem [13, 14], and interesting strategies have been proposed. In particular, Vandesompele et al.  elucidated the errors related to the common practice of single control normalization. The results showed that a conventional normalization strategy based on a single reference gene leads to erroneous normalization and also demonstrated that ideal and universal control genes do not exist. A proposed alternative strategy is based on a gene stability measure (M) as the average pair-wise variation between a particular gene and all other control genes. This approach has led to the GeNorm algorithm, which evaluates the M of common reference genes and allows determination of the most stably expressed genes in the tissue analyzed during stepwise exclusion of the worst scoring gene points . The rationale for this analysis is that a normalization factor based on proper internal control genes should remove all non-specific variation. This robust procedure has led to more accurate gene expression studies; however, the method must be set up for different tissues analyzed and requires a high number of replicates, involving extensive experimental work and twice the amount of sample material. Despite being an important development, the GeNorm strategy does not completely satisfy the need for a universal method for normalization in RT-qPCR.
We therefore focused on the EARs. The location of EARs in UTRs of 1,500 genes makes the development of this novel strategy an important tool for normalization. In fact, even the transcriptional dysregulation that characterizes several pathological conditions does not significantly alter the total amount of EARs in the whole human transcriptome. This scenario applies in HD, in which the activity of several transcription factors and related DNA target sequences is impaired , thus specifically complicating the choice of a suitably stable reference gene both in the control and HD population.
One of the selected mRNA targets in this study, BDNF, is reported to be involved in the pathogenesis of HD [3, 16]. Previous studies have shown that BDNF mRNA is consistently reduced in cells and in vivo brain tissues from HD transgenic mice and patients . In addition, blood BDNF levels correlate with brain BDNF levels and disease progression in HD mice . The protocol described here for the rapid and accurate evaluation of BDNF mRNA in human blood samples represents a step forward in efforts to develop a reliable measurement of key mRNA in blood. Moreover, we clearly demonstrated that mRNA levels normalized over EARs are comparable to those obtained after normalization using the most robust normalization strategies available in both analyzed tissues.
More important, given the widespread distribution of EARs in the human genome and the absence of any tissue specificity in their expression, the EAR normalization tool described here can be universally applied to human biological samples (blood, skeletal muscle and post-mortem brain tissues), avoiding the time-consuming and expensive experimental work needed for selection and analysis of the most stable reference genes in the chosen tissue. It should be noted that the EAR normalization strategy minimizes the amount of RNA used for the analysis, an important advantage for expression studies on samples characterized by a limited amount of RNA where it would otherwise be very difficult to use alternative methods such as GeNorm.
A reliable and flexible normalization strategy has been an urgent need in the context of gene expression studies in human samples. We have identified a suitable novel normalization strategy that can be applied in different types of human biological material. EAR normalization enables the validation of an appropriate control gene for any specific experimental design to be avoided. More important, the strategy of reducing non-specific variation in samples by normalizing to EAR enables a straightforward assessment of relevant biological variation from mRNA levels. This possibility is particularly important in the context of clinical studies for evaluating the alteration of mRNA levels as reliable markers of disease progression.
Materials and methods
Recruitment of healthy volunteers
Characteristics of the study
University College London (UK)
Institute-IRCCS "Carlo Besta" (Italy)
Institute-IRCCS "Carlo Besta" (Italy)
Human tissue collection
A total of 2.5 ml whole venous blood was collected in an evacuated blood collection tube (PAXgene™ Blood RNA Tube, BD Vacutainer™, Plymouth, UK) containing a stabilizing additive. After sample collection, the PAX tubes were incubated and stabilized at room temperature for 2 hours, followed by immediate storage at -20°C until RNA extraction.
Skeletal muscle collection
Biopsies of vastus lateralis muscle were obtained after informed consent from National Neurological Institute-IRCCS "Carlo Besta" (Italy). Biopsies were obtained under local anesthesia and immediately frozen in dry ice or liquid nitrogen.
RNA was extracted from whole blood using the PAXgene™ Blood RNA System Kit (PreAnalytiX, QIAGEN) following the manufacturer's guidelines. Total RNA concentration and purity were measured using a NanoDrop 1000 Spectrophotometer (Thermo Scientific, Waltham, MA, USA), and quality was verified by means of agarose gel electrophoresis of 1 μg of each sample (Additional data file 2). The total RNA was stored in aliquots at -80°C.
The muscle samples were homogenized in TRIZOL (Invitrogen, Carlsbad, CA, USA) using a Ultraturex homogeniser. Total RNA was isolated according to the manufacturer's protocol and stored in aliquots at -80°C.
Using Superscript III RNaseH reverse transcriptase (Invitrogen) and random primers in a volume of 20 μl, 250 ng total RNA was reverse-transcribed to single-stranded cDNA, according to the manufacturer's instructions.
Three independent PCR analyses were performed of each of two independent RT reactions, for a total of six independent measurements for each of the analyzed mRNA samples.
BDNF, IL8, TNNC1 and housekeeping genes
Using an iCycler Thermal Cycler with a Multicolor Real-time PCR Detection System (Bio-Rad, Hercules, CA, USA), all of the reactions were performed in a total volume of 25 μl containing 5 μl of cDNA from a 1:10 RT dilution, 50 mM KCl, 20 mM Tris-HCI, pH 8.4, 0.2 mM dNTPs, iTaq DNA polymerase, 25 units/ml, 3 mM MgCl2, SYBR Green I, 10 nM fluorescein, stabilizers (iQTM SYBR Green Supermix, Bio-Rad), and 0.3 μM of forward and reverse primers. The amplification cycles consisted of an initial denaturing cycle at 95°C for 3 minutes, followed by 45 cycles of 30 s at 95°C, 30 s at 60°C, and 30 s at 72°C. Fluorescence was quantified during the 60°C annealing step, and product formation was confirmed by means of a melting curve analysis (55 to 94°C). The amounts of BDNF and IL8 mRNA were normalized to the geometric mean of the three most stable housekeeping genes (HPRT1, YWHAZ, and GNB2L; Table 1 for primer sequences).
The BDNF-specific primers set up on the BDNF coding sequence (GenBank accession number AF411339) were BDNF 5'-TAACGGCGGCAGACAAAAAGA-3' and BDNF 5'-GAAGTATTGCTTCAGTTGGCCT-3'. The obtained amplification product was 101 bp long. The IL8-specific primers were set up referring to GenBank accession number NM_000584.2; the primer sequences were IL8 5'-CCATCTCACTGTGTGTAAACATGAC-3' and IL8 5'-TCCACTCTCAATCACTCTCAGTTCT-3'. The obtained amplification product was 194 bp long. The TNNC1-specific primers were set up referring to GenBank accession number NM_003280.2; the primer sequences were TNNC1 5'-TGCAGGAGATGATCGATGAGGTG-3' and TNNC1 5'-TGCGGAAGAGGTCAGACAGCTC-3'. The obtained amplification product was 138 bp long.
Using an iCycler Thermal Cycler with a Multicolor Real-time PCR Detection System (Bio-Rad), all of the reactions were performed in a total volume of 25 μl containing 5 μl of cDNA from a 1:200 RT dilution, 50 mM KCl, 20 mM Tris-HCI, pH 8.4, 0.2 mM dNTPs, iTaq DNA polymerase, 25 units/ml, 3 mM MgCl2, SYBR Green I, 10 nM fluorescein, stabilizers (iQTM SYBR Green Supermix, Bio-Rad), and 0.3 μM of forward and reverse primers. The amplification cycles consisted of an initial denaturing cycle at 95°C for 3 minutes, followed by 45 cycles of 30 s at 95°C, 30 s at 60°C, and 30 s at 72°C. Fluorescence was quantified during the 60°C annealing step, and product formation was confirmed by means of a melting curve analysis (55 to 94°C).
The EAR-specific primers were set up on an Alu human interspersed repetitive sequence (Figure 2b). The primer sequences were EAR 5'-GAGGCTGAGGCAGGAGAATCG-3' and EAR 5'-GTCGCCCAGGCTGGAGTG-3'. The obtained amplification product was 87 bp long.
Statistical analyses and data representation
Because BDNF values did not have a normal distribution, the non-parametric Mann-Whitney two-tailed U-test was used for the statistical analyses, with significance set at P < 0.05. Median values were chosen because of the skewed distribution of the data. In each figure, the boundary of the box closest to zero indicates the 25th percentile, the line within the box marks the median or 50th percentile, and the boundary of the box farthest from zero indicates the 75th percentile. When ten or more samples were analyzed, whiskers above and below the box indicate the 90th and 10th percentiles.
The validation of the EAR normalization strategy with respect to the GeNorm algorithm was measured using Pearson correlations (r).
The materials and methods section and nomenclature follows the MIQE reporting guidelines .
brain-derived neurotrophic factor
expressed Alu repeat
guanine nucleotide binding protein (G protein), beta polypeptide 2-like 1
hypoxanthine phosphoribosyltransferase 1
gene expression stability measure
real-time quantitative RT-PCR
troponin C type 1
tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide.
We gratefully acknowledge Marco Bianchi from Bio-Rad Laboratories (Segrate, MI, Italy) for technical advice and for helpful comments on the manuscript. This work was supported by NeuroNe (sixth Framework Programme, European Union, LsHM-CT-2004-512039) and Stem-HD (Framework Programme, Union, LsHB-CT-2006-037349) (E Cattaneo). This paper and protocol was developed based on an original idea and proof of concept of the method presented by J Vandesompele at the 2nd International qPCR Symposium, Freising-Weihenstephan, Germany, September 6, 2005.
- 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: RESEARCH0034.1-RESEARCH0034.11. 10.1186/gb-2002-3-7-research0034.View ArticleGoogle Scholar
- Yulug IG, Yulug A, Fisher EM: The frequency and position of Alu repeats in cDNAs, as determined by database searching. Genomics. 1995, 27: 544-548. 10.1006/geno.1995.1090.PubMedView ArticleGoogle Scholar
- Zuccato C, Cattaneo E: Brain-derived neurotrophic factor in neurodegenerative diseases. Nat Rev Neurol. 2009, 5: 311-322. 10.1038/nrneurol.2009.54.PubMedView ArticleGoogle Scholar
- Fujimura H, Altar CA, Chen R, Nakamura T, Nakahashi T, Kambayashi J, Sun B, Tandon NN: Brain-derived neurotrophic factor is stored in human platelets and released by agonist stimulation. Thromb Haemost. 2002, 87: 728-734.PubMedGoogle Scholar
- Radka SF, Holst PA, Fritsche M, Altar CA: Presence of brain-derived neurotrophic factor in brain and human and rat but not mouse serum detected by a sensitive and specific immunoassay. Brain Res. 1996, 709: 122-130. 10.1016/0006-8993(95)01321-0.PubMedView ArticleGoogle Scholar
- Trajkovska V, Marcussen AB, Vinberg M, Hartvig P, Aznar S, Knudsen GM: Measurements of brain-derived neurotrophic factor: methodological aspects and demographical data. Brain Res Bull. 2007, 73: 143-149. 10.1016/j.brainresbull.2007.03.009.PubMedView ArticleGoogle Scholar
- Batzer MA, Deininger PL: Alu repeats and human genomic diversity. Nat Rev Genet. 2002, 3: 370-379. 10.1038/nrg798.PubMedView ArticleGoogle Scholar
- Zuccato C, Marullo M, Conforti P, MacDonald ME, Tartari M, Cattaneo E: Systematic assessment of BDNF and its receptor levels in human cortices affected by Huntington's disease. Brain Pathol. 2008, 18: 225-238. 10.1111/j.1750-3639.2007.00111.x.PubMedView ArticleGoogle Scholar
- Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation. Biotechniques. 2000, 29: 332-337.PubMedGoogle Scholar
- Strand AD, Aragaki AK, Shaw D, Bird T, Holton J, Turner C, Tapscott SJ, Tabrizi SJ, Schapira AH, Kooperberg C, Olson JM: Gene expression in Huntington's disease skeletal muscle: a potential biomarker. Hum Mol Genet. 2005, 14: 1863-1876. 10.1093/hmg/ddi192.PubMedView ArticleGoogle Scholar
- Bustin SA: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrinol. 2000, 25: 169-193. 10.1677/jme.0.0250169.PubMedView ArticleGoogle Scholar
- Bustin SA: Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J Mol Endocrinol. 2002, 29: 23-39. 10.1677/jme.0.0290023.PubMedView ArticleGoogle Scholar
- Overbergh L, Giulietti A, Valckx D, Decallonne R, Bouillon R, Mathieu C: The use of real-time reverse transcriptase PCR for the quantification of cytokine gene expression. J Biomol Tech. 2003, 14: 33-43.PubMedPubMed CentralGoogle Scholar
- Pachot A, Blond JL, Mougin B, Miossec P: Peptidylpropyl isomerase B (PPIB): a suitable reference gene for mRNA quantification in peripheral whole blood. J Biotechnol. 2004, 114: 121-124. 10.1016/j.jbiotec.2004.07.001.PubMedView ArticleGoogle Scholar
- Cha JH: Transcriptional signatures in Huntington's disease. Prog Neurobiol. 2007, 83: 228-248. 10.1016/j.pneurobio.2007.03.004.PubMedPubMed CentralView ArticleGoogle Scholar
- Zuccato C, Cattaneo E: Role of brain-derived neurotrophic factor in Huntington's disease. Prog Neurobiol. 2007, 81: 294-330. 10.1016/j.pneurobio.2007.01.003.PubMedView ArticleGoogle Scholar
- Conforti P, Ramos C, Apostol BL, Simmons DA, Nguyen HP, Riess O, Thompson LM, Zuccato C, Cattaneo E: Blood level of brain-derived neurotrophic factor mRNA is progressively reduced in rodent models of Huntington's disease: restoration by the neuroprotective compound CEP-1347. Mol Cell Neurosci. 2008, 39: 1-7. 10.1016/j.mcn.2008.04.012.PubMedView ArticleGoogle Scholar
- Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009, 55: 611-622. 10.1373/clinchem.2008.112797.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.