Figures
Abstract
In order to identify the optimal internal control for relative real-time PCR when studying target gene expression in the red alga Porphyra yezoensis, we quantified the expression of seven housekeeping genes (18S ribosomal RNA, 30S ribosomal protein S8, Polyubiquitin-2, Glyceraldehyde-3-phosphate dehydrogenase, Elongation factor 1-alpha, Beta-tubulin and Actin 3) at different life-history stages. Absolute quantification was done by normalization to total RNA quantity and by normalization to genomic DNA quantity. We used these two normalization approaches, comparing the differences of expression levels of all candidate housekeeping genes between any two generations and across three life-history stages (filamentous sporophytes, leafy gametophytes and conchospores). We found GAPDH had the best stability in all cases and we recommend that GAPDH be considered as a potential internal control for gene expression studies at different life-history stages in P. yezoensis.
Citation: Wu X, Huang A, Xu M, Wang C, Jia Z, Wang G, et al. (2013) Variation of Expression Levels of Seven Housekeeping Genes at Different Life-History Stages in Porphyra yezoensis. PLoS ONE 8(4): e60740. https://doi.org/10.1371/journal.pone.0060740
Editor: Debashish Bhattacharya, Rutgers University, United States of America
Received: June 21, 2011; Accepted: March 4, 2013; Published: April 18, 2013
Copyright: © 2013 Wu et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by National Natural Science Foundation of China (B49082401; 30830015)(http://www.nsfc.gov.cn/Portal0/default124.htm). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Porphyra (Bangiales, Rhodophyta) includes more than 100 species distributed worldwide [1]. Owing to its economic importance and health benefits [2]–[4], Porphyra is treated as essential material in practical application research [5], [6]. Porphyra grows on intertidal marine rocks, surviving under extreme conditions such as high temperature, strong light and desiccation for many hours each day [7] and is, therefore, suitable experimental material for understanding how intertidal organisms respond to various forms of stress [8], [9]. Moreover, the multiphasic life history of Porphyra can be observed in the laboratory within a few months [10], accelerating research on growth, reproduction and photosynthetic mechanisms [11]–[13]. These characteristics have prompted increasing attention to Porphyra, which has been proposed as a model system for physiological and genetic investigation [14], [15].
The red alga Porphyra yezoensis is one of the most extensively cultivated species in Japan, China and South Korea. It has a complex heteromorphic life cycle with an alternation between macroscopic haploid leafy gametophytes and microscopic filamentous sporophytes. Shell-boring filamentous sporophytes release conchospores that can germinate and grow into leafy gametophytes [14], [16], [17]. It is apparent that different life-history stages of P. yezoensis exhibit different morphological traits and occur under distinct growth conditions; the gametophytes live freely in seawater and the sporophytes penetrate the shells. One question that needs to be answered is whether expression levels of housekeeping genes are constant in the various life-history stages?
The analysis of target gene expression at different life-history stages is an increasingly significant field of life science research [11], [18], [19]. In order to conduct a gene expression assay it is crucial to choose a suitable internal control gene that is expressed stably and, thus, the selection of this internal control gene is an absolute prerequisite for accurate quantification of target gene expression [20]–[22]. In fact, it has become clear that there is no gene that can be used as internal control for all species [23]–[26]. However, no ideal internal control for gene expression analysis in P. yezoensis has been identified. In addition, it is worth noting that results obtained by normalization to a single internal control are generally less rigorous or even inaccurate, and that results obtained by normalization to two or more different internal controls are contradictory in some cases. Owing to the complex life-history of P. yezoensis, this problem is particularly pertinent to the investigation of gene expression. The requirement for suitable internal control genes for normalization is increasingly stringent for in-depth research into P. yezoensis, for which the mechanism underlying gene regulation and expression across the entire life-history is extremely important.
Our objective was to identify internal control genes suitable for the expression profile of a target at different life-history stages. We chose seven housekeeping genes commonly used as likely candidates and the transcript numbers of these genes were determined with absolute quantification. The results given here serve as a reliable reference for selection of internal controls in future experimental designs.
Materials and Methods
P. yezoensis culture
The filamentous sporophytes were cultured in sterile seawater with added nutrients with constant aeration at 20°C with 50°μmol of photons m−2 s−1 (12 h light/12 h dark). The nutrient-enriched seawater was changed every week. Basically, the shell-boring sporophytes were incubated in fresh filtered seawater at 24°C with PES medium under 20 μmol of photons m−2 s−1 with a 12 h light/12 h dark cycle. After a few days, the temperature was reduced gradually to 20°C and we started to collect conchospores from the shells. Before collection, the shells were brushed with filtered seawater to remove any contamination. The seawater containing conchospores was poured slowly into a glass beaker, then transferred to a 1.5 ml microcentrifuge tube and centrifuged (Eppendorf Centrifuge 5804R; Eppendorf AG, Hamburg, Germany) at 4000 g for 5 min at room temperature [27]. The conchospores that had already been isolated were immediately frozen in liquid nitrogen and stored at –80°C. The leafy gametophytes were cultivated at 15°C with 50 μmol of photons m−2 s−1 in 12 h light/12 h dark cycle.
RNA extraction and cDNA synthesis
Total RNA was extracted from samples at different life-history stages using a Tiangen RNAprep pure plant kit (Tiangen, Beijing) according to the manufacturer's instructions. It should be noted that RNA extractions from sporophytes and gametophytes were replicated three times. However, although collected every day, the amount of conchospores was insufficient to complete replicate experiments. For this reason, RNA extraction from conchospores was conducted only one time.
Before the extraction of RNA, each sample was dried with sterile filter paper to avoid degradation, ground to a powder in liquid nitrogen and then weighed. Purified RNA was eluted with RNase-free water and immediately stored at –80°C. RNA integrity was checked by electrophoresis in a 2% (w/v) agarose gel, RNA purity was estimated by measuring OD 260/280 absorbance ratio, and the concentration of RNA was determined with a Qubit® Fluorometer and the corresponding RNA Assay Kit (Invitrogen, USA). RNA yield was calculated as RNA quantity/sample weight. RNase-free DNase I (Promega, USA) was used to eliminate genomic DNA contamination before RT-PCR. First-strand cDNA was synthesized using M-MLV reverse transcriptase (Promega, USA) with random primers on a PCR apparatus (Eppendorf, Genman) following the protocols suggested by the manufacturer. The cDNA product was stored at –20°C.
Genomic DNA extraction
Genomic DNA was extracted from samples at different life-history stages using the Universal Genomic DNA Extraction Kit (TaKaRa Biotechnology, Dalian, China) according to the manufacturer's instructions. Similarly, DNA extractions from sporophytes and gametophytes were replicated three times, and DNA extraction from conchospores was conducted only one time. In order to avoid error caused by grinding, samples of the same powder were used for extraction of RNA and DNA. In addition, the amount of materials used to extract RNA and DNA also was the same in terms of weight. DNA quality was evaluated by measuring OD 260/280 absorbance ratio. The concentration of genomic DNA was determined with a Qubit® Fluorometer and the corresponding DNA Assay Kit (Invitrogen, USA). DNA yield was calculated as DNA quantity/sample weight.
Primer design
The sequences of commonly used housekeeping genes 18S, RPS8, PUB-2, GAPDH, EF1alpha, TubB and Act3 were downloaded from the Genbank database (http://www.ncbi.nlm.nih.gov). Primer design and evaluation were completed with Primer Premier 5.0 (Premier Biosoft International, Palo Alto, CA) and Oligo Primer analysis software (Wojciech and Piotr Rychlik Copyright, version 6.31). The criteria used for primer design included a primer length of 17–21 bp, a PCR product size of 104–187 bp and other general rules (Table 1).
Construction of DNA standards
The fragment of each housekeeping gene was amplified by PCR, using the specific primers (Table 2). PCR amplicon standard was purified using the E.Z.N.A TM Gel Extraction Kit (OMEGA Bio-Tek, Doraville, GA, USA) and then quantified with a Qubit® fluorometer and the corresponding dsDNA HS Assay Kit (Invitrogen Crop., Carlsbad, CA, USA). Transcript number was calculated as:(1)
Quantification of the tested genes in various samples
Real-time PCR was done with a BIO-RAD IQ5 real-time PCR detection system (Bio-Rad, USA) containing 9 μl of 2× SYBR® Green Master Mix (Tiangen, Beijing), 2 μl of template (standard diluted in series or cDNA diluted 10-fold), 5 μl of sense and antisense primer (2 μM) and made to a final volume of 20 μl with RNase-free water. Conditions were: 3 min at 94°C, followed by 40 cycles of 15 s at 94°C, 40 s at 59°C and 30 s at 72°C. The melting curve was generated by heating for 30 s from 55°C to 95°C, with a ramp speed of 0.5°C/cycle. A 10-fold series dilution of the standard was used to construct the standard curves [28]–[30] by plotting the logarithm of the threshold cycles against the logarithm of the transcript number of the template. The levels of all unknown samples should be distributed within the range of the standard curve. Each sample was amplified in triplicate. Amplification efficiency (E) and the linear correlation coefficient (R2) were used to evaluate primer efficiency:(2)
On the basis of the Ct values obtained, the transcript numbers of all tested genes in different samples were calculated and normalized to total RNA quantity and to genomic DNA quantity. The expression stability of all candidate housekeeping genes was evaluated by comparing the differences of observed transcript numbers between any two generations and across three life-history stages [31]. The entire protocol, from nucleic acid extraction to Real-time PCR, was replicated.
Results
Quality and yields of total RNA and genomic DNA
An OD 260/280 ratio of 1.8–2.0 is usually considered as an acceptable indicator of good nucleic acid quality. In this study, OD 260/280 ratios of total RNA and genomic DNA extracted from different samples varied between 1.8 and 2.0 (Table 3).
The yields of total RNA and genomic DNA were calculated by dividing quantity by the corresponding sample weight. The mean RNA yield was highest in the sporophytes (0.158±0.059 µg/mg) and lowest in the conchospores (0.069 µg/mg), whereas the mean DNA yield was highest in the conchospores (0.088 µg/mg) and lowest in the gametophytes (0.028±0.002 µg/mg) (Figure 1).
Primers specificity and efficiency
In this study, the amplification efficiency of all PCRs was between 90% and 105%, and R2 was >0.980 (Figure 2). In the melting curve analysis, the amplicons of all seven genes revealed a single product with a melting temperature in accord with the expected value (Figure 3).
The results showed that amplification efficiency was between 96% and 103%, and linear correlation coefficient was >0.99.
Melting peaks were examined with standard samples and unkown samples (sporophytes, gametophytes and conchospores). The melting curve for each gene had only one peak.
Expression profiling and stability of candidate housekeeping genes: normalization to total RNA
On the basis of standard curves, the transcript number of the housekeeping gene in different samples was obtained by normalizing to total RNA. The expression levels of the genes were 18S > GAPDH > RPS8 > > EF1alpha > PUB-2> TubB > Act3, range 2.25×103–3.47×107 (based on the expression level in the sporophytes), and the expression of 18S was the highest at all life-history stages. The trend in expression levels was the same for all candidate housekeeping genes at different life-history stages: conchospores > sporophytes > gametophytes (Table 4).
We ranked the tested genes by the difference between the smallest and the largest transcript number of each gene across all samples (Table 5). GAPDH had the smallest difference (2.2-fold), and there was little variation across all samples. The difference for 18 S varied greatly, ranging from 1.11×107 to7.77×109 copies. Act3, EF1alpha and RPS8 showed huge variance with >100-fold differences. The difference was 21.5-fold for TubB and 91.1-fold for PUB-2.
We calculated the differences of expression values between any two generations. The differences between the expression levels of all tested genes in the sporophytes and in the gametophytes were relatively small, ranging from 1.1-fold for GAPDH to 8.1-fold for PUB-2 (Table 6). According to the difference between the sporophytes and the conchospores, we found that GAPDH was still the most satisfactory of the tested genes; values for the other tested genes ranged from 9.9-fold to 224.0-fold (Table 7). Ranking by the differences between the gametophytes and the conchospores was identical with that based on differences across all samples (Table 8).
Expression profiling and stability of candidate housekeeping genes: normalization to genomic DNA
The expression levels of the tested genes were 18S > GAPDH > RPS8 > > EF1alpha > PUB-2> TubB > Act3, ranging from 5.86×103 to 8.06×107 (based on the expression level in the sporophytes). The transcript abundance of all the tested genes was greatest for the conchospores, except GAPDH, which was expressed highest in the gametophytes (Table 9).
We compared the difference of expression levels of the tested genes between any two generations and across three life-history stages. The expression levels of all the tested genes were similar between the sporophytes and the gametophytes, ranging from 1.4-fold to 5.4-fold (Table 6). By comparing the transcript profile of sporophytes or gametophytes against conconspores and analyzing the difference of transcript numbers across three life-history stages, we found that 18S, Act3, EF1alpha, PUB-2, RPS8 and TubB were not good choices as internal controls in these cases. However, GAPDH displayed the greatest level of stability (Table 5, Table 7, Table 8).
Discussion
It is essential to identify an appropriate internal control for accurate analysis of gene expression at different life-history stages and in various tissues [32], [33]. At present, 18S ribosomal RNA, glyceraldehyde-3-phosphate dehydrogenase and β-actin are most commonly selected as internal controls for higher animals and plants [34]–[38]. Other housekeeping genes, such as Elongation factor and 30S Ribosomal Protein Subunit, have been proposed to be suitable for normalization in real-time PCR experiments under certain conditions [23], [38]. Generally, in such studies, differences in gene expression have been compared in different organs and tissues or under various experimental conditions but at the life-history stage had the same ploidy. Although expression profiling can be done by normalizing against the above-mentioned internal control genes, there are still researchers advocating the use of two or more internal control genes, in order to ensure accuracy [21], [39], [40]. Thus, for P. yezoensis with a complex life history involving changes in ploidy, it is necessary to choose optimal internal controls for gene expression studies with caution.
Absolute quantification requires standardization, and several approaches have been suggested [30], [41]. Ideally, the gene transcript number is standardized to the number of cells [21], [28]; however, P. yezoensis is a multicellular organism and accurate enumeration of cells is impossible. A strategy based on total RNA mass quantity is widely used [42]–[44] and we suggest normalization against the genomic DNA quantity of tested samples could be used. In order to identify a suitable internal control for P. yezoensis, we have quantified the expression of seven housekeeping genes at different life-history stages with absolute quantification by normalization to total RNA quantity and to genomic DNA quantity.
Besides, absolute quantification requires the construction of a standard curve for each target gene, using DNA standard. Although it had been demonstrated that the amplification efficiency was unchanged by using circular or linear DNA as the standard [29], the types of DNA standard (circular plasmid, linearized plasmid and linear PCR amplicon) affected quantification accuracy. A recent study showed that a circular plasmid was unsuitable as a standard, by which gene transcript number was overestimated seriously, whereas the linear standards gave highly accurate estimates [45]. In this study, a linear PCR amplicon was used as the standard for constructing standard curve due to its high flexibility and sensitivity. An ideally suited internal control gene is one that shows no, or only a limited, variation of expression across the sample set [46]. Consequently, we chose the relatively stably expressed gene GAPDH as a suitable internal control for relative quantification under certain experimental designs. No matter which normalization method was applied, when differences of gene expression levels between any two generations or across three generations were compared, the ranking of housekeeping genes was consistent, identifying GAPDH as the most stable of the tested genes with differences ranging from 1.4- fold to 2.0-fold. However, the expression levels of the other six tested genes varied greatly and to different degrees. Based on the above results, we believed it was appropriate to use GAPDH as the internal control in this study.
It should be pointed out, however, that the results showed enormous expression differences exist in housekeeping genes between vegetative cells and conchospores. Released conchospores are looking for matrix to attach rapidly by amoeboid movement. In amoeboid movement, rich expression of actin is significant for cytoskeletal polymerization and depolymerization. Simultaneously, ribosome and some related functional proteins are synthesized abundantly to meet the need of the movement. Therefore, it is understandable that there is large expression differences between vegetative and conchospores.
In conclusion, we recommend that GAPDH be used as the internal control for gene expression studies at different life-history stages of P. yezoensis, allowing accurate quantification of target gene expression by real-time quantitative PCR experimental designs.
Author Contributions
Conceived and designed the experiments: GCW JFN. Performed the experiments: XJW JFN MLX CW ZJJ. Analyzed the data: XJW AYH JFN. Contributed reagents/materials/analysis tools: JFN. Wrote the paper: XJW JFN GCW.
References
- 1. Yoshida T, Notoya M, Kikuchi N, Miyata M (1997) Catalogue of species of Porphyra in the world, with special reference to the type locality and bibliography. Nat Hist Res, Special Issue 3: 5–18.
- 2. Fleurence J (1999) Seaweed proteins: biochemical, nutritional aspects and potential uses. Trends Food Sci Tech 10: 25–28.
- 3.
FAO (2006) State of World Aquaculture 2006. Rome, Italy: Food & Agriculture Org.
- 4. MacArtain P, Gill CIR, Brooks M, Campbell R, Rowland IR (2007) Nutritional value of edible seaweeds. Nutr Rev 65: 535–543.
- 5. Daggett TL, Pearce CM, Tingley M, Robinson SMC, Chopin T (2005) Effect of prepared and macroalgal diets and seed stock source on somatic growth of juvenile green sea urchins (Strongylocentrotus droebachiensis). Aquaculture 244: 263–281.
- 6. De la Coba F, Aguilera J, De Gálvez MV, Álvarez M, Gallego E, et al. (2009) Prevention of the ultraviolet effects on clinical and histopathological changes, as well as the heat shock protein-70 expression in mouse skin by topical application of algal UV-absorbing compounds. J Dermatol Sci 55: 161–169.
- 7. Blouin NA, Brodie JA, Grossman AC, Xu P, Brawley SH (2010) Porphyra: a marine crop shaped by stress. Trends Plant Sci 16: 29–37.
- 8. Reed RH, Collins JC, Russell G (1980) The Effects of Salinity upon Galactosyl-Glycerol Content and Concentration of the Marine Red Alga Porphyra purpurea (Roth) C. Ag. J Exp Bot 31: 1539–1554.
- 9. Wiencke C, Läuchli A (1980) Growth, cell volume, and fine structure of Porphyra umbilicalis in relation to osmotic tolerance. Planta 150: 303–311.
- 10. SAHOO D, TANG X, YARISH C (2002) Porphyra: the economic seaweed as a new experimental system. Curr Sci 83: 1313–1316.
- 11. Asamizu E, Nakajima M, Kitade Y, Saga N, Nakamura Y, et al. (2003) Comparison of RNA expression profiles between the two generations of Porphyra yezoensis (Rhodophyta), based on expressed sequence tag frequency analysis. J Phycol 39: 923–930.
- 12. Sidirelli-Wolff M (1992) The influence of temperature, irradiance and photoperiod on the reproductive life history of Porphyra leucosticta (Bangiales, Rhodophyta) in laboratory culture. Bot Mar 35: 251–257.
- 13. Sheath RG, Hellebust JA, Sawa T (1977) Changes in plastid structure, pigmentation and photosynthesis of the conchocelis stage of Porphyra leucosticta (Rhodophyta, Bangiophyceae) in response to low light and darkness. Phycologia 16: 265–275.
- 14. Kitade Y, Taguchi G, Shin JA, Saga N (1998) Porphyra monospore system (Bangiales, Rhodophyta): A model for the developmental biology of marine plants. Phycol Res 46: 17–20.
- 15. Fukuda S, Mikami K, Uji T, Park EJ, Ohba T, et al. (2008) Factors influencing efficiency of transient gene expression in the red macrophyte Porphyra yezoensis. Plant Sci 174: 329–339.
- 16. Yamazaki A, Nakanishi K, Saga N (1998) Axenic tissue culture and morphogenesis in Porphyra yezoensis (Bangiales, Rhodophyta). J Phycol 34: 1082–1087.
- 17. Kitade Y, Fukuda S, Saga N (1999) Preliminary study on early development of monospore of Porphyra yezoensis (Bangiales, Rhodophyta). Fish Genet Breed Sci 28: 27–34.
- 18. Fukuda S, Kitade Y, Miyamoto H, Nagashima S, Takahashi S, et al. (2003) Isolation and characterization of an elongation factor-1α gene in Porphyra yezoensis (Rhodophyta). J Appl Phycol 15: 81–86.
- 19. Kitade Y, Fukuda S, Nakajima M, Watanabe T, Saga N (2002) Isolation of a cDNA encoding a homologue of actin from Porphyra yezoensis (Rhodophyta). J Appl Phycol 14: 135–141.
- 20.
Pfaffl MW (2004) Quantification strategies in real-time PCR. In: Bustin SA, editor. A-Z of quantitative PCR. La Jolla, California, USA: International University Line. 89–113.
- 21. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3: research0034–research0034.11.
- 22. Bustin SA (2002) Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J Mol Endocrinol 29: 23–39.
- 23. Ransbotyn V, Reusch TBH (2006) Housekeeping gene selection for quantitative real-time PCR assays in the seagrass Zostera marina subjected to heat stress. Limnol Oceanogr Meth 4: 367–373.
- 24. Jorgensen SM, Kleveland EJ, Grimholt U, Gjoen T (2006) Validation of reference genes for real-time polymerase chain reaction studies in Atlantic salmon. Mar Biotechnol 8: 398–408.
- 25. Nicot N, Hausman JF, Hoffmann L, Evers D (2005) Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress. J Exp Bot 56: 2907–2914.
- 26. Savli H, Karadenizli A, Kolayli F, Gundes S, Ozbek U, et al. (2003) Expression stability of six housekeeping genes: a proposal for resistance gene quantification studies of Pseudomonas aeruginosa by real-time quantitative RT-PCR. J Med Microbiol 52: 403–408.
- 27. Zhang BY, Yang F, Wang GC, Peng G (2010) CLONING AND QUANTITATIVE ANALYSIS OF THE CARBONIC ANHYDRASE GENE FROM PORPHYRA YEZOENSIS. J Phycol 46: 290–296.
- 28. Wong ML, Medrano JF (2005) Real-time PCR for mRNA quantitation. Biotechniques 39: 75–85.
- 29. Whelan JA, Russell NB, Whelan MA (2003) A method for the absolute quantification of cDNA using real-time PCR. J Immunol Meth 278: 261–269.
- 30. ustin SA (2000) Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrinol 25: 169–193.
- 31. Giricz O, Lauer-Fields JL, Fields GB (2008) The normalization of gene expression data in melanoma: investigating the use of glyceraldehyde 3-phosphate dehydrogenase and 18S ribosomal RNA as internal reference genes for quantitative real-time PCR. Anal Biochem 380: 137–139.
- 32. Jain M, Nijhawan A, Tyagi AK, Khurana JP (2006) Validation of housekeeping genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochem Biophys Res Comm 345: 646–651.
- 33. Brunner AM, Yakovlev IA, Strauss SH (2004) Validating internal controls for quantitative plant gene expression studies. BMC Plant Biol 4: 14.
- 34. Stordeur P, Poulin LF, Craciun L, Zhou L, Schandené L, et al. (2002) Cytokine mRNA quantification by real-time PCR. J Immunol Meth 259: 55–64.
- 35. Huang T, Böhlenius H, Eriksson S, Parcy F, Nilsson O (2005) The mRNA of the Arabidopsis gene FT moves from leaf to shoot apex and induces flowering. Science 309: 1694–1696.
- 36. Iskandar HM, Simpson RS, Casu RE, Bonnett GD, Maclean DJ, et al. (2004) Comparison of reference genes for quantitative real-time polymerase chain reaction analysis of gene expression in sugarcane. Plant Mol Biol Rep 22: 325–337.
- 37. Jung M, Ramankulov A, Roigas J, Johannsen M, Ringsdorf M, et al. (2007) In search of suitable reference genes for gene expression studies of human renal cell carcinoma by real-time PCR. BMC Mol Biol 8: 47.
- 38. Siaut M, Heijde M, Mangogna M, Montsant A, Coesel S, et al. (2007) Molecular toolbox for studying diatom biology in Phaeodactylum tricornutum. Gene 406: 23–35.
- 39. Remans T, Smeets K, Opdenakker K, Mathijsen D, Vangronsveld J, et al. (2008) Normalisation of real-time RT-PCR gene expression measurements in Arabidopsis thaliana exposed to increased metal concentrations. Planta 227: 1343–1349.
- 40. Tricarico C, Pinzani P, Bianchi S, Paglierani M, Distante V, et al. (2002) Quantitative real-time reverse transcription polymerase chain reaction: normalization to rRNA or single housekeeping genes is inappropriate for human tissue biopsies. Anal Biochem 309: 293–300.
- 41. Huggett J, Dheda K, Bustin S, Zumla A (2005) Real-time RT-PCR normalisation; strategies and considerations. Gene Immun 6: 279–284.
- 42. Tichopad A, Pfaffl MW, Didier A (2003) Tissue-specific expression pattern of bovine prion gene: quantification using real-time RT-PCR. Mol Cell Probes 17: 5–10.
- 43. Marelli SP, Terova G, Cozzi MC, Lasagna E, Sarti FM, et al. (2010) GENE EXPRESSION OF HEPATIC GLUCOCORTICOID RECEPTOR NR3C1 AND CORRELATION WITH PLASMATIC CORTICOSTERONE IN ITALIAN CHICKENS. Anim Biotechnol 21: 140–148.
- 44. Lončar J, Popović M, Zaja R, Smital T (2010) Gene expression analysis of the ABC efflux transporters in rainbow trout (Oncorhynchus mykiss). Comp Biochem Physiol C Toxicol Pharmacol 151: 209–215.
- 45. Hou Y, Zhang H, Miranda L, Lin S (2010) Serious overestimation in quantitative PCR by circular (supercoiled) plasmid standard: microalgal pcna as the model gene. Plos One 5: e9545.
- 46. Andersen CL, Jensen JL, Ørntoft TF (2004) 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. Canc Res 64: 5245–5250.