Comparison of transcriptome responses of the liver, tail fin, and olfactory epithelium of Rana [Lithobates] catesbeiana tadpoles disrupted by thyroid hormones and estrogen

Abstract

Thyroid hormones (THs) are important developmental regulators in vertebrates, including during the metamorphosis of a tadpole into a frog. Metamorphosis is a post-embryonic developmental period initiated by TH production in the tadpole thyroid gland. The two main bioactive forms of TH are L-thyroxine (T4) and 3,5,3’-triiodothyronine (T3); these hormones have overlapping but distinct mechanisms of action. Premetamorphic tadpoles are highly responsive to TH and can be induced to metamorphose through exogenous TH exposure, making them an important model for both the study of vertebrate TH signaling and endocrine disrupting chemicals (EDCs). It is important to differentiate TH-mediated responses from estrogenic responses in premetamorphic tadpoles when assessing dysregulation by EDCs as crosstalk between the two endocrine systems is well-documented. Herein, we compare the RNA-sequencing-derived transcriptomic profiles of three TH-responsive tissues (liver, olfactory epithelium, and tail fin) in premetamorphic bullfrog (Rana [Lithobates] catesbeiana) tadpoles exposed to T3, T4, and estradiol (E2). These profiles were generated using the latest available genome assembly for the species. The data indicate that there is a clear distinction, and little overlap, between the transcriptomic responses elicited by E2 and the THs. In contrast, within the THs, the T3- and T4-induced transcriptomic profiles generally show considerable overlap; however, the degree of overlap is highly tissue-dependent, illustrating the importance of distinguishing the two THs and the affected signaling pathways within the target tissue type when evaluating hormone active agents. The data herein also show that E2 and TH treatment can uniquely induce significant changes in expression of their respective “classic” bioindicator transcripts vtg (E2) and thra, thrb, and thibz (THs). However, care must be taken in the interpretation of increased vep or esr1 transcripts as a change in transcript levels can be induced by THs rather than solely E2.