Thyroid hormones (THs) play crucial roles in various developmental and physiological processes in vertebrates, including squamate reptiles. The effect of THs on shedding frequency is interesting in Squamata, since the effects on lizards are quite the reverse of those in snakes: injection of thyroxine increases shedding frequency in lizards, but decreases it in snakes. However, the mechanism underlying this differential effect remains unclear. To facilitate the investigation of the molecular mechanism of the physiological functions of THs in Squamata, their two specific receptor (TRαand β) cDNAs, which are members of the nuclear hormone receptor superfamily, were cloned from a lizard, the leopard gecko, Eublepharis macularius. This is the first molecular cloning of thyroid hormone receptors (TRs) from reptiles. The deduced amino acid sequences showed high identity with those of other species, especially in the C and E/F domains, which are characteristic domains in nuclear hormone receptors. Expression analysis revealed that TRs were widely expressed in many tissues and organs, as in other animals. To analyze their role in the skin, temporal expression analysis was performed by RT-PCR, revealing that the two TRs had opposing expression patterns: TRαwas expressed more strongly after than before skin shedding, whereas TRβ was expressed more strongly before than after skin shedding. This provides good evidence that THs play important roles in the skin, and that the roles of their two receptor isoforms are distinct from each other.
The thyroid hormones (THs), thyroxine (T4) and thyronine (T3), are pleiotropic factors important for many developmental and physiological processes in vertebrates. There has been a lot of research into the physiological significance of THs in various vertebrates. For example, THs are known to be important for inner ear and retina development, liver metabolism in mice (Flamant and Samarut, 2003), metamorphosis in axolotl and Xenopus (Nakajima et al., 2005; Sachs et al., 2000; Safi et al., 2004), and embryogenesis and metamorphosis in many teleost fish (Power et al., 2001).
In reptiles, THs have been suggested to affect tail regeneration (Turner and Tipton, 1971), metabolic rate and metabolic enzyme activity (John-Alder, 1990; John-Alder and Joos, 1991), and shedding frequency (Chiu et al., 1967; Chiu and Lynn, 1970). Above all, their effect on shedding frequency is particularly interesting. In lizards, the injection of thyroxine increases shedding frequency, and thyroidectomy decreases it (Chiu et al., 1967). In contrast, in snakes, the injection of thyroxine decreases the shedding frequency, and thyroidectomy increases it (Chiu and Lynn, 1970). The mechanisms underlying these completely opposite phenomena have not been clarified, partly due to the lack of investigation of the molecular mechanism of THs in reptiles.
THs can regulate target genes by interacting with thyroid hormone receptors (TRs), which are members of the nuclear receptor superfamily. Two isoforms, TRα and TRβ, have been isolated from species of four classes of vertebrate, but not from reptiles (Forrest et al., 1990; Kawakami et al., 2003; Murray et al., 1988; Yaoita et al., 1990). These isoforms share high homology and have similar biochemical properties. However, they have distinct spatial and temporal expression profiles in overlapping patterns, suggesting that two genes mediate both individual and common biological functions.
Unlike other squamate animals, the leopard gecko, Eublepharis macularius, is easily maintained and bred in the laboratory. The leopard gecko is therefore expected to become an experimental model. Indeed, several molecular studies of the endocrine system have already been conducted on this species (Endo and Park, 2004; Endo and Park, 2005; Ikemoto and Park, 2003; Ikemoto et al., 2004; Kato et al., 2005; Valleley et al., 2001).
In this study, we cloned TRα and β from the leopard gecko to augment investigations on the molecular mechanisms of the physiological functions of THs in reptiles. In addition to identifying two isoforms of TR, we performed phylogenetic and expression analyses. We also demonstrated the differential expression of the TR isoforms, and herein discuss their possible roles in shedding.
MATERIALS AND METHODS
The leopard geckos (Eublepharis macularius) were treated according to the guidelines of the Bioscience Committee at the University of Tokyo. The animals were provided meal worms, crickets, water, and powdered calcium supplement ad libitum. Animals were anesthetized with sodium pentobarbital and killed by rapid decapitation, followed by complete bloodletting. Tissues and organs were immediately dissected, frozen in liquid nitrogen, and stored at −80°C until use.
RNA preparation and cDNA synthesis
Total RNA was extracted using ISOGEN (NIPPON GENE, Tokyo, Japan). The cDNAs used as templates for RT-PCR were synthesized from 3 μg of denatured total RNA using 5 μM oligo(dT) primer and 100 units of M-MLV Reverse Transcriptase (Promega, Madison, WI) in a 20 μl reaction volume with incubation at 42°C for 1.5 h. The cDNA used for rapid amplification of cDNA ends (RACE) was synthesized from 3 μg of total RNA using a SMART RACE cDNA Amplification Kit (BD Biosciences Clontech, Palo Alto, CA, USA) according to the manufacturer's instructions.
Molecular cloning of TR cDNAs by RT-PCR and RACE
RT-PCR was carried out to obtain partial TR cDNAs from skin cDNA using degenerate primers (Table 1). All of the following PCR amplifications were performed in 20 μl reaction volumes containing each primer at 1 μM, 0.25 unit of TaKaRa Ex Taq (TaKaRa, Shiga, Japan), each dNTP at 250 μM, and Ex Taq Buffer (TaKaRa). The PCR product was separated by electrophoresis, extracted using a QIAquick Gel Extraction Kit (QIAGEN K.K., Tokyo, Japan) and directly sequenced. After determination of the partial sequence, RACE was carried out to determine the complete sequence of TR cDNAs. PCR and nested PCR were performed with gene-specific primers (Table 1) in combination with Universal Primer A Mix (Clontech) or Nested Universal Primer A (Clontech). The amplified products were sequenced as described above. This procedure was repeated independently at least twice to avoid PCR amplification errors.
Oligonucleotide primers used for RACE, RT-PCR, and sequencing.
Comparison of the amino acid sequences of various TRs
CLUSTAL X software (version 1.81) (Thompson et al., 1997) was used with default settings to align the deduced amino acid sequences of the TRs of the leopard gecko and other species.
Amino acid identities were calculated for the C domain, E/F domain, and entire ORF.
Molecular phylogenetic analysis
The nucleotide sequences of the entire ORFs of the TRs from the leopard gecko and from several species representing all other vertebrate classes were aligned using CLUSTAL X with default settings. The alignment of the nucleotide sequences was used to generate a phylogenetic tree, using the neighbor-joining method (Saitou and Nei, 1987). Bootstrap values were calculated with 1000 replications to estimate the robustness of internal nodes. The Gen-Bank accession numbers of TRs used in the phylogenic analysis are as follows: Homo sapiens (human) TRα, M24748; Homo sapiens TRβ, X04707; Mus musculus (mouse) TRα, MMCERBA1; Mus musculus TRβ, S62756; Gallus gallus (chicken) TRα, Y00987; Gallus gallus TRβ, X17504; Xenopus laevis (African clawed frog) TRα,M35344; Xenopus laevis TRβ, M35361; Ambystoma mexicanum (axolotl) TRα, AY174871; Ambystoma mexicanum TRβ, AY174872; Salmo salar (salmon) TRα, AF146775; Salmo salar TRβ,AF302251; Conger myriaster (conger eel) TRαΑ;,AB183396; Conger myriaster TRβ1, AB183394.
Expression analysis of TRs
To identify the target organs of THs, the spatial expression pattern of the TRs was examined by RT-PCR. Twenty-five nanograms of cDNA from the whole brain, heart, liver, small intestine, large intestine, testis, ovary, and thymus, and 5 ng from the pituitary, were amplified using primers specific for TRs. The primer sets used were TRαSE03 and TRαAS02 for TRα, and TRβSE03 and TRβAS03 for TRβ (Table 1). The PCR products were visualized by electrophoresis on a 1.2% TAE agarose gel and stained with ethidium bromide. Each DNA fragment was extracted from the gel and directly sequenced to confirm its identity.
Temporal expression analysis of TRs in skin
Total RNA was extracted from the skin of three animals within 24 hours before or after shedding. cDNA from the skin (7.5 ng) was amplified using the specific primers described above. The PCR conditions were as follows: 94°C for 3 min; 30 cycles of 94°C for 40 s, 64°C for 25 s, and 72°C for 30 s; and 72°C for 2 min. The PCR products were analyzed by electrophoresis on a 1.2% TAE agarose gel.
Molecular cloning of TR cDNAs from the leopard gecko
Full-length TR cDNAs were isolated from the skin of the leopard gecko by RT-PCR and RACE. TRα cDNA comprised 1,744 bp, which included a 5′-UTR of 254 bp, an ORF of 1,227 bp encoding 408 amino acid residues, and a 3′-UTR of 263 bp. TRβcDNA comprised 1,481bp, including a 5′-UTR of 104 bp, an ORF of 1,110 bp encoding 369 amino acid residues, and a 3′-UTR of 267 bp. The domains are indicated in Fig. 1. Sequences of full-length cDNAs were deposited in GenBank (Accession Nos. AB204861 and AB204862).
Comparison of the amino acid sequences of various TRs
Alignments of the predicted amino acid sequence of leopard gecko TRs (lgTRs) with those of other species are shown in Fig. 2. Across the entire ORF, the lgTRs showed very high identity (87–96% for TRα and 77–98% for TRβ) with their homologs from other species. When the specific domains were considered, stronger conservation was observed in the C domain (91–97% between lgTRα and its vertebrate homologs; 95–99% for TRβ) and in the E/F domain (94–95% between lgTRα and its vertebrate homologs; 92–99% for TRβ).
Molecular phylogenetic analysis
A phylogenetic tree of the TRs was constructed using the entire ORF nucleotide sequences of selected species representing all classes of vertebrate (Fig. 3). The TRs formed two groups, TRα and TRβ, in accordance with these two isoforms being derived from distinct genes. As expected, both the leopard gecko TRs clustered with their chicken homologs.
As expected, a wide TR distribution was observed. RT-PCR products of the expected size were obtained from all tissues and organs examined (Fig. 4). The authenticity of the RT-PCR products was confirmed by direct sequencing. A control without RT was also used for 40 or 45 cycles of PCR, and no signal was detected (data not shown).
The THs are pleiotropic factors important for many functions in vertebrates. In reptiles, THs are suggested to affect tail regeneration, metabolic rate, metabolic enzyme activity, and shedding frequency. To augment the investigation of the molecular mechanism of THs in reptiles, we characterized their receptors at the molecular level.
In this study, we cloned two isoforms of TR from the leopard gecko, Eublepharis macularius. This is the first molecular identification of full-length TRs from reptiles. The deduced amino acid sequences of the cDNAs demonstrated the classic modular structure of members of the nuclear receptor superfamily, and exhibited high identity with their homologs in other species. The highest identities were shown with their chicken homologs (96 and 98%). There are eleven conserved Cys residues in the C domain of TRα. Although the seventh Cys in lgTRα, which does not contribute to the formation of a disulfide bond or zinc finger (Zhao et al., 1998), is substituted by Ser, other Cys residues are completely conserved in all TRαs. It is therefore conceivable that cloned lgTRαdoes not lose the ability to bind DNA by this substitution. In the future, ligand binding studies will be helpful. In the phylogenetic tree, TRs formed groups, TRαand TRβ. Both lgTRs clustered with their chicken counterparts, as expected.
mRNA expression of the lgTRs was detected in all the tissues and organs examined, indicating that thyroid hormones are pleiotropic factors important for many functions in the leopard gecko as well as other animals. In lizards, it has been reported that THs can regulate cardiac function (Venditti et al., 1996), enzyme activity in the liver (John-Alder, 1990), and testis activity (Cardone et al., 2000; Plowman and Lynn, 1973). It is therefore conceivable that TRs mediate such effects in these organs. In other species, the expression of TRs has been also demonstrated in various tissues, but significant isoform-specific functions are poorly understood. For instance, although the expression of TRs in the adult brain has been demonstrated in mammals, their specific roles have not yet been clarified (Schwartz et al., 1992). It is known that THs can regulate steroidogenesis; however, although the expression of TRs in the human ovary has been confirmed, their specific roles remain unknown (Zhang et al., 1997). There is less information for the testis, and the type of TR expressed there remains controversial (Maran, 2003).
THs are known to regulate the shedding frequency in Squamata. Intriguingly, the effect of thyroxine appears to be reversed between lizards and snakes. To obtain a better understanding of the potential role of THs in skin shedding, we analyzed the temporal expression of TRs in the skin. The expression of TRα was stronger after skin shedding than before, whereas the result was the opposite for TRβ. Although Chiu et al. 1967 have discussed the indirect effect of THs on shedding frequency, our results strongly suggest that THs can directly affect the skin, and that the two isoforms of TR play distinct roles in skin shedding.
The shedding cycle can be divided into two phases: resting and proliferation. As the skin is in the resting phase after shedding (Maderson and Licht, 1967), our result suggests that TRα plays a role in the resting phase, such as maintaining this condition so as not to enter the proliferation phase. Furthermore, we suspect that TRβ mediates the effect of THs in the proliferation phase, since TRβ was strongly expressed in the skin before shedding. It is conceivable that the condition of the skin taken before shedding in this experiment was at around the last stage of the proliferation phase (Maderson and Licht, 1967). This is supported by an in vitro study demonstrating that the epidermis can differentiate by itself but cannot shed (Flexman et al., 1968). This indicates that the capacity for the complex changing pattern of cell differentiation is intrinsic to the epidermis, but that shedding is not. Extrinsic factor(s) is/are necessary for shedding, and THs may be one such factor. It has also been reported that THs have no effect on the skin in the proliferation phase, either directly or indirectly, as even thyroidectomized animals shed (Chiu et al., 1967). This discrepancy may be because it is always difficult to completely remove the thyroid surgically (Chiu et al., 1967). The interpretation of the physiological significance of up- or down-regulation of mRNA expression of TRs needs further study, such as in situ hybridization to analyze where and when during the shedding cycle TRs are expressed. The differential expression of TR isoforms has also been reported for other species, such as frogs during development (Sachs et al., 2000). Although both isoforms appear to be involved in regulating metamorphosis, their functional differences are not yet clear.
In this report, we characterized leopard gecko TRs and demonstrated the possibility of their direct involvement in skin shedding. These results will facilitate investigation of the physiological significance of THs and of the molecular mechanisms by which they regulate shedding frequency in Squamata.
We are grateful to Prof. Y. Oka, Dr. Y. Akazome, Dr. H. Abe, Ms. M. Kyokuwa, Mr. M. Enomoto, Mr. T. Ikemoto, Mr. K. Kato, and Ms. M. Utsumi, Department of Biological Sciences, Graduate School of Science, The University of Tokyo, for valuable discussions throughout this study. This work was supported by Grantsin-Aid from the Japan Society for the Promotion of Science.
- A. Cardone, F. Angelini, T. Esposito, R. Comitato, and B. Varriale . 2000. The expression of androgen receptor messenger RNA is regulated by triiodothyronine in lizard testis. J Steroid Biochem Mol Biol 72:133–141. Google Scholar
- K. W. Chiu and W. G. Lynn . 1970. The role of thyroid in skin-shedding in the shovel-nosed snake, Chionactis occipitalis. Gen Comp Endocrinol 14:467–474. Google Scholar
- K. W. Chiu, J. G. Phillips, and P. F. A. Maderson . 1967. The role of the thyroid in the control of the sloughing cycle in the tokay (Gecko gecko, lacertilian). J Endocrinol 39:463–472. Google Scholar
- D. Endo and M. K. Park . 2004. Molecular characterization of the leopard gecko POMC gene and expressional change in the testis by acclimation to low temperature and with a short photoperiod. Gen Com Endocrinol 138:70–77. Google Scholar
- D. Endo and M. K. Park . 2005. Molecular cloning of P450 aromatase from the leopard gecko and its experssion in the ovary. J Steroid Biochem Mol Biol 96:131–140. Google Scholar
- F. Flamant and J. Samarut . 2003. Thyroid hormone receptors: lessons from knockout and knock-in mutant mice. Trends Endocrinol Metab 14:85–90. Google Scholar
- B. A. Flexman, P. F. A. Maderson, G. Szabo, and S. I. Roth . 1968. Control of cell differentiation in lizard epidermis in vitro. Develop Biol 18:354–374. Google Scholar
- D. Forrest, M. Sjöberg, and B. Vennström . 1990. Contrasting developmental and tissue-specific expression of alpha and beta thyroid hormone receptor genes. EMBO J 9:1519–1528. Google Scholar
- T. Ikemoto and M. K. Park . 2003. Identification and characterization of the reptilian GnRH-II gene in the leopard gecko, Eublepharis macularius, and its evolutionary consideration. Gene 31:157–165. Google Scholar
- T. Ikemoto, M. Enomoto, and M. K. Park . 2004. Identification and characterization of a reptilian GnRH receptor from the leopard gecko. Mol Cell Endocrinol 214:137–147. Google Scholar
- H. B. John-Alder 1990. Thyroid regulation of resting metabolic rate and intermediary metabolic enzymes in a lizard (Sceloporus occidentalis). Gen Comp Endocrinol 77:52–62. Google Scholar
- H. B. John-Alder and B. Joos . 1991. Interactive effects of thyroxine and experimental location on running endurance, tissue masses, and enzyme activities in captive versus field-active lizards (Sceloporus undulatus). Gen Comp Endocrinol 81:120–132. Google Scholar
- K. Kato, T. Ikemoto, and M. K. Park . 2005. Identification of the reptilian prolactin and its receptor cDNAs in the leopard gecko, Eublepharis macularius. Gene 346:267–276. Google Scholar
- Y. Kawakami, M. Tanda, S. Adachi, and K. Yamauchi . 2003. Characterization of thyroid hormone receptor alpha and beta in the metamorphosing Japanese conger eel, Conger myriaster. Gen Comp Endocrinol 132:321–332. Google Scholar
- P. F. A. Maderson and P. Licht . 1967. Epidermal morphology and sloughing frequency in normal and prolactin treated Anolis carolinensis (Iguanidae, Lacerilia). J Morph 123:157–172. Google Scholar
- R. R. M. Maran 2003. Thyroid hormones: their role in testicular steroidogenesis. Arch Androl 49:375–388. Google Scholar
- M. B. Murray, N. D. Zilz, N. L. McCreary, M. J. MacDonald, and H. C. Towle . 1988. Isolation and characterization of rat cDNA clones for two distinct thyroid hormone receptors. J Biol Chem 263:12770–12777. Google Scholar
- K. Nakajima, K. Fujimoto, and Y. Yaoita . 2005. Programmed cell death during amphibian metamorphosis. Semin Cell Dev Biol 16:271–280. Google Scholar
- M. M. Plowman and W. G. Lynn . 1973. The role of the thyriod in testicular function in the Gecko, Coleonys variegatus. Gen Comp Endocrinol 20:342–246. Google Scholar
- D. M. Power, L. Llewellyn, M. Faustino, M. A. Nowell, B. T. Bjornsson, I. E. Einarsdottir, A. V. Canario, and G. E. Sweeney . 2001. Thyroid hormones in growth and development of fish. Comp Biochem Physiol C Toxicol Pharmacol 130:447–459. Google Scholar
- L. M. Sachs, S. Damjanovski, P. L. Jones, Q. Li, T. Amano, S. Ueda, Y. B. Shi, and A. Ishizuya-Oka . 2000. Dual functions of thyroid hormone receptors during Xenopus development. Comp Biochem Physiol B Biochem Mol Biol 126:199–211. Google Scholar
- R. Safi, S. Bertrand, O. Marchand, M. Duffraisse, A. Luze, J. M. Vanacker, M. Maraninchi, A. Margotat, B. Demeneix, and V. Laudet . 2004. The axolotl (Ambystoma mexicanum), a neotenic amphibian, expresses functional thyroid hormone receptors. Endocrinology 145:760–772. Google Scholar
- N. Saitou and M. Nei . 1987. The neighbor-joining method: a new method for reconstructing phylogenic trees. Mol Biol Evol 4:406–425. Google Scholar
- H. L. Schwartz, K. A. Strait, N. C. Ling, and J. H. Oppenheimer . 1992. Quantitation of rat tissue thyroid hormone binding receptor isoforms by immunoprecipitation of nuclear triiodothyronine binding capacity. J Biol Chem 267:11794–11799. Google Scholar
- J. D. Thompson, T. J. Gibson, F. Plewniak, F. Jeanmouginm, and D. G. Higgins . 1997. The CLUSTAL X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25:4876–4882. Google Scholar
- J. E. Turner and S. R. Tipton . 1971. The role of the lizard thyroid gland in tail regeneration. J Exp Zool 178:63–86. Google Scholar
- E. M. Vallelay, E. J. Cartwright, N. J. Croft, A. F. Markham, and P. L. Coletta . 2001. Characterisation and expression of Sox9 in the Lepard gecko, Eublepharis macularius. J Exp Zool 291:85–91. Google Scholar
- P. Venditti, S. D. Meo, and P. M. Rosaroll . 1996. Effect of T3 administration on electrophysiological properties of lizard ventricular muscle fibers. J Comp Physiol B 165:552–557. Google Scholar
- Y. Yaoita, Y. B. Shi, and D. D. Brown . 1990. Xenopus laevis alpha and beta thyroid hormone receptors. Proc Natl Acad Sci USA 87:7090–7094. Google Scholar
- S. S. Zhang, A. J. Carrillo, and D. S. Darling . 1997. Expression of multiple thyroid hormone receptor mRNAs in human oocytes, cumulus cells, and granulosa cells. Mol Hum Reprod 3:555–562. Google Scholar
- Q. Zhao, S. Khorasanizadeh, Y. Miyoshi, M. A. Lazar, and F. Rastinejad . 1998. Structural elements of an orphan nuclear receptor-DNA complex. Mol Cell 1:849–861. Google Scholar