[en] Reference genes are essential for gene expression analysis when using real-time quantitative PCR (RT-qPCR). Xenopus laevis is a popular amphibian model for studying vertebrate embryogenesis and development. Further, X. laevis is ideal for studying thyroid signaling due to its thyroid dependent metamorphosis, a stage comparable to birth in humans. When using PCR based studies, a primary concern is the choice of reference genes. Commonly used references are eef1a1, odc1, rpl8, and actnB, although there is a lack of ad hoc reference genes for X. laevis. Here, we used previously published RNA-seq data on different X. laevis stages and identified the top 14 candidate genes with respect to their expression levels as a function of developmental stage and degree of variation. We further evaluated the stability of these and other candidate genes using RT-qPCR on various stages including the unfertilised eggs, whole embryos during early development and brains during late development. We used four different statistical software packages: deltaCT, geNorm, NormFinder and BestKeeper. We report optimized reference gene pair combinations for studying development (early whole embryos), brains at later stages (metamorphosis and adult), and thyroid signalling. These reference gene pairs are suitable for studying different aspects of X. laevis development and organogenesis.
Gurdon, J. B. & Hopwood, N. The introduction of Xenopus laevis into developmental biology: Of empire, pregnancy testing and ribosomal genes. Int. J. Dev. Biol. 44, 43-50 (2000).
Hellsten, U. et al. The Genome of the Western Clawed Frog Xenopus tropicalis. Science (80-). 328, 633-636 (2010).
Session, A. M. et al. Genome evolution in the allotetraploid frog Xenopus laevis. Nature 538, 336-343 (2016).
Chervoneva, I. et al. Selection of optimal reference genes for normalization in quantitative RT-PCR. BMC Bioinformatics 11, 253 (2010).
Coulson, D. T. et al. Identification of valid reference genes for the normalization of RT qPCR gene expression data in human brain tissue. BMC Mol. Biol. 9, 46 (2008).
Casadei, R. et al. Identification of housekeeping genes suitable for gene expression analysis in the zebrafish. Gene Expr. Patterns 11, 271-276 (2011).
Šindelka, R., Ferjentsik, Z. & Jonák, J. Developmental expression profiles of Xenopus laevis reference genes. Dev. Dyn. 235, 754-758 (2006).
Dhorne-Pollet, S., Thélie, A. & Pollet, N. Validation of novel reference genes for RT-qPCR studies of gene expression in Xenopus tropicalis during embryonic and post-embryonic development. Dev. Dyn. 242, 709-717 (2013).
Moe, T. K., Ziliang, J., Barathi, A. & Beuerman, R. W. Differential expression of glyceraldehyde-3-phosphate dehydrogenase (GAPDH), beta actin and hypoxanthine phosphoribosyltransferase (HPRT) in postnatal rabbit sclera. Curr. Eye Res. 23, 44-50 (2001).
Sarropoulou, E., Nousdili, D., Kotoulas, G. & Magoulas, A. Functional Divergences of GAPDH Isoforms During Early Development in Two Perciform Fish Species. Mar. Biotechnol. 13, 1115-1124 (2011).
Chambers, S. A. M. & Shaw, B. R. Levels of histone H4 diacetylation decrease dramatically during sea urchin embryonic development and correlate with cell doubling rate. J. Biol. Chem. 259, 13458-13463 (1984).
Hu, Y., Xie, S. & Yao, J. Identification of novel reference genes suitable for qRT-PCR normalization with respect to the zebrafish developmental stage. PLoS One 11, 1-13 (2016).
Kudo, T. et al. Identification of reference genes for quantitative expression analysis using large-scale RNA-seq data of <i>Arabidopsis thaliana</i> and model crop plants. Genes Genet. Syst. 91, 111-125 (2016).
Yim, A. K.-Y. et al. Using RNA-Seq Data to Evaluate Reference Genes Suitable for Gene Expression Studies in Soybean. PLoS One 10, e0136343 (2015).
Cheung, T. T., Weston, M. K. & Wilson, M. J. Selection and evaluation of reference genes for analysis of mouse (Mus musculus) sex-dimorphic brain development. PeerJ 5, e2909 (2017).
Silver, N., Best, S., Jiang, J. & Thein, S. L. Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol. Biol. 7, 33 (2006).
Vandesompele, J. et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 3, RESEARCH0034 (2002).
Pfaffl, M. W., Tichopad, A., Prgomet, C. & Neuvians, T. P. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper-Excel-based tool using pair-wise correlations. Biotechnol. Lett. 26, 509-515 (2004).
Andersen, C. L., Jensen, J. L. & Ørntoft, T. F. 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. 64, 5245-5250 (2004).
Nieuwkoop, P. D. & Faber, J. Normal Table of Xenopus Laevis (Daudin): A Systematical & Chronological Survey of the Development from the Fertilized Egg till the End of Metamorphosis. Http://Www.Xenbase.Org/Anatomy/Alldev.Do 33, (Amsterdam: North-Holland Publishing Company. Guilders, 1994).
Lim, W., Nguyen, N.-H., Yang, H. Y., Scanlan, T. S. & Furlow, J. D. A Thyroid Hormone Antagonist That Inhibits Thyroid Hormone Action in Vivo. J. Biol. Chem. 277, 35664-35670 (2002).
Fini, J.-B. et al. Human amniotic fluid contaminants alter thyroid hormone signalling and early brain development in Xenopus embryos. Sci. Rep. 7, 43786 (2017).
Untergasser, A. et al. Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 35, W71-4 (2007).
Dubois, G. M. et al. Deiodinase activity is present in Xenopus laevis during early embryogenesis. Endocrinology 147, 4941-4949 (2006).
Karpinka, J. B. et al. Xenbase, the Xenopusmodel organism database; new virtualized system, data types and genomes. Nucleic Acids Res. 43, D756-D763 (2015).
Veldhoen, N. et al. Corrigendum to 'The bactericidal agent triclosan modulates thyroid hormone-associated gene expression and disrupts postembryonic anuran development' [Aquat. Toxicol. 80 (2006) 217-227] (https://doi.org/10.1016/j.aquatox.2006.08.010). Aquatic Toxicology 83, 84 (2007).
Crofton, K. M., Paul, K. B., DeVito, M. J. & Hedge, J. M. Short-term in vivo exposure to the water contaminant triclosan: Evidence for disruption of thyroxine. Environ. Toxicol. Pharmacol. 24, 194-197 (2007).
Koeppe, E. S., Ferguson, K. K., Colacino, J. A. & Meeker, J. D. Relationship between urinary triclosan and paraben concentrations and serum thyroid measures in NHANES 2007-2008. Sci. Total Environ. 445-446, 299-305 (2013).
Fini, J. B. et al. Thyroid hormone signaling in the Xenopus laevis embryo is functional and susceptible to endocrine disruption. Endocrinology 153, 5068-81 (2012).
Burmeister, L. A., Pachucki, J. & St. Germain, D. L. Thyroid hormones inhibit type 2 iodothyronine deiodinase in the rat cerebral cortex by both pre-and posttranslational mechanisms. Endocrinology 138, 5231-5237 (1997).
Gereben, B. The mRNA Structure Has Potent Regulatory Effects on Type 2 Iodothyronine Deiodinase Expression. Mol. Endocrinol. 16, 1667-1679 (2002).
Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2-ΔΔCT Method.