Transcriptomic profiling of Rana [Lithobates] catesbeiana back skin during natural and thyroid hormone-induced metamorphosis under different temperature regimes with particular emphasis on innate immune system components

As amphibians undergo thyroid hormone (TH)-dependent metamorphosis from an aquatic tadpole to the terrestrial frog, their innate immune system must adapt to the new environment. Skin is a primary line of defense, yet this organ undergoes extensive remodelling during metamorphosis and how it responds to TH is poorly understood. Temperature modulation, which regulates metamorphic timing, is a unique way to uncover early TH-induced transcriptomic events. Metamorphosis of premetamorphic tadpoles is induced by exogenous TH administration at 24 °C but is paused at 5 °C. However, at 5 °C a “molecular memory” of TH exposure is retained that results in an accelerated metamorphosis upon shifting to 24 °C. We used RNA-sequencing to identify changes in Rana (Lithobates) catesbeiana back skin gene expression during natural and TH-induced metamorphosis. During natural metamorphosis, significant differential expression (DE) was observed in >6500 transcripts including classic TH-responsive transcripts (thrb and thibz), heat shock proteins, and innate immune system components: keratins, mucins, and antimicrobial peptides (AMPs). Premetamorphic tadpoles maintained at 5 °C showed 83 DE transcripts within 48 h after TH administration, including thibz which has previously been identified as a molecular memory component in other tissues. Over 3600 DE transcripts were detected in TH-treated tadpoles at 24 °C or when tadpoles held at 5 °C were shifted to 24 °C. Gene ontology (GO) terms related to transcription, RNA metabolic processes, and translation were enriched in both datasets and immune related GO terms were observed in the temperature-modulated experiment. Our findings have implications on survival as climate change affects amphibia worldwide.