Validation of Reference Genes for Normalization Gene Expression in Reverse Transcription Quantitative PCR in Human Normal Thyroid and Goiter Tissue

Weber, Raquel; Bertoni, Ana Paula Santin; Bessestil, Laura Walter; Brasil, Beatriz Maria de Azevedo Assis; Brum, llma Simoni; Furlanetto, Tania Weber

doi:https://doi.org/10.1155/2014/198582

BioMed Research International

On this page

Abstract Introduction Methods Results Discussion Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2014 | Article ID 198582 | https://doi.org/10.1155/2014/198582

Validation of Reference Genes for Normalization Gene Expression in Reverse Transcription Quantitative PCR in Human Normal Thyroid and Goiter Tissue

Raquel Weber,¹Ana Paula Santin Bertoni,^1,2Laura Walter Bessestil,³Beatriz Maria de Azevedo Assis Brasil,⁴llma Simoni Brum,³and Tania Weber Furlanetto¹

Academic Editor: Muy-Teck Teh

Received28 Feb 2014

Accepted18 Apr 2014

Published11 May 2014

Abstract

Reverse transcription quantitative polymerase chain reaction (RT-qPCR) has been recognized as the most accurate method for quantifying mRNA transcripts, but normalization of samples is a prerequisite for correct data interpretation. So, this study aimed to evaluate the most stable reference gene for RT-qPCR in human normal thyroid and goiter tissues. Beta-actin (ACTB); glyceraldehyde-3-phosphate dehydrogenase (GAPDH); succinate dehydrogenase, subunit A, flavoprotein (Fp) (SDHA); hypoxanthine phosphoribosyltransferase I (HPRTI); tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide (YWHAZ); and beta-2-microglobulin (B2M) were evaluated in 14 thyroid tissue samples (7 normal and 7 goiter tissues) by RT-qPCR. The mean Cq and the maximum fold change (MFC) and NormFinder software were used to assess the stability of the genes. As a result, ACTB gene was more stable than GAPDH, SDHA, HPRTI, YWHAZ, and B2M. In conclusion, ACTB could be used to normalize RT-qPCR data in normal thyroid and goiter tissues.

1. Introduction

Reverse transcription quantitative polymerase chain reaction (RT-qPCR) has been recognized as the most accurate, sensitive, and easy method for quantifying mRNA transcripts in biological samples [1, 2]. Nevertheless, normalization of samples is a prerequisite for correct data interpretation, since its accuracy is significantly affected by sample quality, reagent, and operator technique [3].

Cell number, quantification of rRNA or total RNA, and reference genes are some standardization methods that have been proposed so far [2, 4, 5]. Normalization with another gene, which is assumed to have stable expression, is the method of choice because this reference gene is exposed to the same preparation steps as the gene of interest [6]. Therefore, a full procedure control is obtained during the whole assay.

A good reference gene needs to present low variation in its expression in the different samples which would be evaluated and demonstrated to be minimally regulated during the individual experiment or pathological condition. Therefore, there is a need to validate genes to allow the accuracy of RNA transcription analysis [3, 6]. Nowadays, there is emerging evidence that reference genes could vary their expression in different cellular development stages and in other conditions [5, 7, 8].

The thyroid normal gland is a fairly homogenous structure, but in goiter a genetic heterogeneity is described in follicular cells; some nodules have a polyclonal origin and others a monoclonal origin. These changes may be related to mutations in oncogenes which do not produce malignancy per se but would alter growth and function [14, 15]. These alterations could modify the expression of many thyroid genes, including the ones usually described as reference genes, so the evaluation of candidates to normalize RT-qPCR is necessary.

Little is known about the expression stability of reference genes in thyroid cells or tissue, as shown in Table 1.

Considering the species, there are many differences in reference genes stability evaluated in these studies. Although in human tissues the same software was used to analyze data, differences in RT-qPCR and cells status could have contributed to these findings. In addition, a stable expression of a reference gene in one tumor type does not predict a stable expression in another tumor type [7]. In this way, there is no universal reference gene being necessary to validate potential reference genes for any experimental condition. So, this study aimed to evaluate the most stable reference gene for RT-qPCR in human normal thyroid and goiter tissues.

2. Material and Methods

2.1. Tissue Acquisition

Normal thyroid and goiter tissues were obtained from patients who underwent total thyroidectomy as part of treatment for differentiated thyroid cancer in the Hospital de Clínicas de Porto Alegre (HCPA). After the evaluation of macroscopic and frozen sections of surgical specimens by two pathologists, samples of normal thyroid or goiter tissue were frozen in liquid nitrogen and stored at −80°C until further processing. The pattern used to confirm the presence of goiter was a well-defined fibrous capsule with a mixture of macrofollicles and microfollicles and, in some cases, some degenerative changes such as fibrosis and hemorrhage [16]. This study was approved by the Ethics Committee of HCPA (GPPG: 12-0272).

2.2. Nucleic Acid Extraction and Reverse Transcription

Normal thyroid and goiter tissue were homogenized mechanically with Omnimix for 30 seconds in Trizol (Invitrogen, Life Technologies, Carlsbad, CA, USA), and total RNA was extracted with the same commercial kit, according to the manufacturer’s protocol. The RNA concentration and purity were assessed with the NanoDrop 1000 Spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA). Samples were then stored at −80°C. 1 μg total RNA was reverse-transcribed to produce cDNA using random and oligo-dT primers and Superscript II reverse transcriptase (Invitrogen Life Technologies) in total reaction volumes of 20 μL. All cDNA samples were diluted 10-fold with diethyl pyrocarbonate- (DEPC-) treated water and stored at −20°C.

2.3. Gene Selection and Quantitative PCR

Six genes commonly used as reference in RT-qPCR gene expression experiments were selected (Table 1): beta-actin (ACTB, related to cell structure); glyceraldehyde-3-phosphate dehydrogenase (GAPDH, related to carbohydrate metabolism); succinate dehydrogenase, subunit A, flavoprotein (Fp) (SDHA, related to energy metabolism); hypoxanthine phosphoribosyltransferase I (HPRTI, related to nucleotide metabolism); tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide (YWHAZ, related to cell growth and death); and beta-2-microglobulin (B2M, related to major histocompatibility complex). Primers sequence shown in Table 2 was described previously [17]. Quantification of amplified samples was performed based on a standard curve with five successive tenfold dilution points of a pool of cDNA samples.

RT-qPCR was performed using Applied Biosystems StepOne Real-Time PCR System using Kit Platinum SYBR Green qPCR SuperMix-UDG (Invitrogen Life Technologies, Carlsbad, CA, USA). The RT-qPCR was performed with an initial denaturation (10 min at 95°C), followed by 40 cycles of denaturation (15 s at 95°C), annealing, and extension (45 s at 60°C). Melting curve was acquired to ensure that a single product was amplified in each reaction.

2.4. Statistics

To evaluate the stability of the candidate reference genes in thyroid normal and goiter tissues, the NormFinder algorithm was used [18]. In addition, the analysis of raw quantification cycle (Cq) values of each gene was used to evaluate their stability. Mean Cq values, standard deviation (SD), coefficient of variation (CV), and maximum fold change (MFC, the ratio of the maximum and minimum values observed within the dataset) were calculated.

3. Results

Six reference genes were amplified in 14 thyroid samples: 7 normal and 7 goiter. All RT-qPCR assays produced a single peak in the melting curve. The mean Cq cycle values, SD, CV, and MFC obtained for each gene in goiter and normal thyroid are described in Table 3. Analyzing these data, none of the genes had a constant expression. Moreover, B2M showed higher variability in Cq values. Thereby, we proceeded the analysis of quantification data using the NormFinder algorithm. The values of stability of the candidate genes obtained from the NormFinder analysis are shown in Table 4. The most stable genes were ACTB for goiter and GAPDH for normal thyroid tissue. When all samples were analyzed together, the most stable genes were SDHA and ACTB. Stability values and intra/intergroup variations are shown, respectively, in Table 5 and Figure 1.

4. Discussion

In this study, we evaluated 6 reference genes in 14 thyroid specimens (7 normal thyroid and 7 goiter) using RT-qPCR. The genes with the lowest variations in expression for normal samples, goiter, and all samples were, respectively, GAPDH, ACTB, SDHA and ACTB, when evaluated by the NormFinder.

According to de Jonge et al. [19], a good candidate gene should present a small coefficient of variation and a MFC < 2. Moreover, a mean expression level lower than the maximum expression level subtracted with 2 SD is a prerequisite for a candidate reference gene. By these criteria, only B2M and GAPDH genes could be used as internal control for goiter, and none of these genes could be used for normal thyroid, although all had MFC < 2.

The NormFinder algorithm was used to rank the candidate reference genes based on their stability values. The gene expression is more stable when this value is closer to zero. The stability value less than 0.15 is the cutoff for an acceptable reference gene [4]. Considering these criteria, none of the selected genes could be used to normalize RT-qPCR data for goiter tissues, and GAPDH and SDHA could be used for normal thyroid. When all data were analyzed together, SDHA and ACTB genes had the lowest stability values, respectively, 0.076 and 0.082. We suggest using ACTB as reference gene, when comparing gene expression in normal thyroid and goiter, because this gene had the lowest intra- and intergroup variations. Our results differ somewhat from those of Chantawibul et al. [10] who found the highest stability value for HPRT1, while the GAPDH had the lowest stability value. These differences may partly be explained by the sample selection, which included goiters, adenomas, carcinomas, and lymphocytic thyroiditis.

ACTB mRNA is expressed at moderately abundant levels in most cell types and its protein plays a key role in the cytoskeleton maintenance. Recently, some studies have suggested that this gene does not satisfy the basic requirements for application as internal control due to change in expression under various biomedical stimuli [20, 21]. Nevertheless, in normal thyroid cells, ACTB expression appeared to be stable in response to different experimental treatments [9]. Other studies have demonstrated that ACTB can be used as reference gene [22–24] for different cell types.

One limitation of this study is the small number of samples of goiter; nevertheless, all samples were evaluated by two pathologists, who used the same criteria to define goiter.

In conclusion, the results of the present study suggest that ACTB gene is more stable than SDHA, GAPDH, HPRTI, YWHAZ, and B2M when evaluating human normal thyroid and goiter together.

Conflict of Interests

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

Acknowledgments

The authors gratefully acknowledge FIPE/HCPA, CAPES/PROF, and CNPq for financial support.

References

M. W. Pfaffl, A. Tichopad, C. Prgomet, and T. P. Neuvians, “Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper—Excel-based tool using pair-wise correlations,” Biotechnology Letters, vol. 26, no. 6, pp. 509–515, 2004.
View at: Publisher Site | Google Scholar
S. A. Bustin, “Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems,” Journal of Molecular Endocrinology, vol. 29, no. 1, pp. 23–39, 2002.
View at: Publisher Site | Google Scholar
S. A. Bustin, V. Benes, J. A. Garson et al., “The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments,” Clinical Chemistry, vol. 55, no. 4, pp. 611–622, 2009.
View at: Publisher Site | Google Scholar
J. Vandesompele, K. De Preter, F. Pattyn et al., “Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes,” Genome Biology, vol. 3, no. 7, Article ID RESEARCH0034, 2002.
View at: Publisher Site | Google Scholar
O. Thellin, W. Zorzi, B. Lakaye et al., “Housekeeping genes as internal standards: use and limits,” Journal of Biotechnology, vol. 75, no. 2-3, pp. 291–295, 1999.
View at: Publisher Site | Google Scholar
A. Radonić, S. Thulke, I. M. Mackay, O. Landt, W. Siegert, and A. Nitsche, “Guideline to reference gene selection for quantitative real-time PCR,” Biochemical and Biophysical Research Communications, vol. 313, no. 4, pp. 856–862, 2004.
View at: Publisher Site | Google Scholar
F. Haller, B. Kulle, S. Schwager et al., “Equivalence test in quantitative reverse transcription polymerase chain reaction: confirmation of reference genes suitable for normalization,” Analytical Biochemistry, vol. 335, no. 1, pp. 1–9, 2004.
View at: Publisher Site | Google Scholar
J. Dundas and M. Ling, “Reference genes for measuring mRNA expression,” Theory in Biosciences, vol. 131, no. 4, pp. 215–223, 2012.
View at: Google Scholar
A. P. Santin, A. F. Souza, L. S. Brum, and T. W. Furlanetto, “Validation of reference genes for normalizing gene expression in real-time quantitative reverse transcription PCR in human thyroid cells in primary culture treated with progesterone and estradiol,” Molecular Biotechnology, vol. 54, no. 2, pp. 278–282, 2013.
View at: Publisher Site | Google Scholar
S. Chantawibul, A. Anuwong, and K. Leelawat, “Validation of appropriate reference genes for gene expression studies in human thyroid gland using real-time RT-PCR,” Journal of the Medical Association of Thailand, vol. 95, supplement 3, pp. S36–S40, 2012.
View at: Google Scholar
C. Lecchi, F. Dilda, P. Sartorelli, and F. Ceciliani, “Widespread expression of SAA and Hp RNA in bovine tissues after evaluation of suitable reference genes,” Veterinary Immunology and Immunopathology, vol. 145, no. 1-2, pp. 556–562, 2012.
View at: Publisher Site | Google Scholar
Q. Li, K. J. Domig, T. Ettle, W. Windisch, C. Mair, and K. Schedle, “Evaluation of potential reference genes for relative quantification by RT-qPCR in different porcine tissues derived from feeding studies,” International Journal of Molecular Sciences, vol. 12, no. 3, pp. 1727–1734, 2011.
View at: Publisher Site | Google Scholar
P. Lisowski, M. Pierzchała, J. Gościk, C. S. Pareek, and L. Zwierzchowski, “Evaluation of reference genes for studies of gene expression in the bovine liver, kidney, pituitary, and thyroid,” Journal of Applied Genetics, vol. 49, no. 4, pp. 367–372, 2008.
View at: Google Scholar
H. Studer, H. J. Peter, and H. Gerber, “Natural heterogeneity of thyroid cells: the basis for understanding thyroid function and nodular goiter growth,” Endocrine Reviews, vol. 10, no. 2, pp. 125–135, 1989.
View at: Google Scholar
K. Krohn, D. Führer, Y. Bayer et al., “Molecular pathogenesis of euthyroid and toxic multinodular goiter,” Endocrine Reviews, vol. 26, no. 4, pp. 504–524, 2005.
View at: Publisher Site | Google Scholar
Z. W. Baloch and V. A. LiVolsi, “Follicular-patterned lesions of the thyroid: the bane of the pathologist,” American Journal of Clinical Pathology, vol. 117, no. 1, pp. 143–150, 2002.
View at: Publisher Site | Google Scholar
A. F. Souza, I. S. Brum, B. S. Neto, M. Berger, and G. Branchini, “Reference gene for primary culture of prostate cancer cells,” Molecular Biology Reports, vol. 40, no. 4, pp. 2955–2962, 2012.
View at: Publisher Site | Google Scholar
C. L. Andersen, J. L. Jensen, and T. F. Orntoft, “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 Research, vol. 64, no. 15, pp. 5245–5250, 2004.
View at: Publisher Site | Google Scholar
H. J. de Jonge, R. S. N. Fehrmann, E. S. J. M. de Bont et al., “Evidence based selection of housekeeping genes,” PLoS ONE, vol. 2, no. 9, article e898, 2007.
View at: Publisher Site | Google Scholar
S. Selvey, E. W. Thompson, K. Matthaei, R. A. Lea, M. G. Irving, and L. R. Griffiths, “β-actin—an unsuitable internal control for RT-PCR,” Molecular and Cellular Probes, vol. 15, no. 5, pp. 307–311, 2001.
View at: Publisher Site | Google Scholar
W. Ruan and M. Lai, “Actin, a reliable marker of internal control?” Clinica Chimica Acta, vol. 385, no. 1-2, pp. 1–5, 2007.
View at: Publisher Site | Google Scholar
R. Mori, Q. Wang, K. D. Danenberg, J. K. Pinski, and P. V. Danenberg, “Both β-actin and GAPDH are useful reference genes for normalization of quantitative RT-PCR in human FFPE tissue samples of prostate cancer,” Prostate, vol. 68, no. 14, pp. 1555–1560, 2008.
View at: Publisher Site | Google Scholar
R. Mehta, A. Birerdinc, N. Hossain et al., “Validation of endogenous reference genes for qRT-PCR analysis of human visceral adipose samples,” BMC Molecular Biology, vol. 11, article 39, 2010.
View at: Publisher Site | Google Scholar
R. Hayashi, S. Hayashi, K. Arai, M. Chikuda, and Y. Obara, “Effects of antioxidant supplementation on mRNA expression of glucose-6-phosphate dehydrogenase, β-actin and 18S rRNA in the anterior capsule of the lens in cataract patients,” Experimental Eye Research, vol. 96, no. 1, pp. 48–54, 2012.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2014 Raquel Weber et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

3993

Downloads

1212

Citations