Identification of thyroid hormone response genes in the remodeling of dorsal muscle during Microhyla fissipes metamorphosis

Introduction Extensive morphological, biochemical, and cellular changes occur during anuran metamorphosis, which is triggered by a single hormone, thyroid hormone (TH). The function of TH is mainly mediated through thyroid receptor (TR) by binding to the specific thyroid response elements (TREs) of direct response genes, in turn regulating the downstream genes in the cascade. The remodeling of dorsal skeletal muscle during anuran metamorphosis provides the perfect model to identify the immediate early and direct response genes that are important during apoptosis, proliferation, and differentiation of the muscle. Methods In our current study, we performed Illumina sequencing combined with single-molecule real-time (SMRT) sequencing in the dorsal muscle of Microhyla fissipes after TH, cycloheximide (CHX), and TH_CHX treatment. Results and Discussion We first identified 1,245 differentially expressed transcripts (DETs) after TH exposure, many of which were involved in DNA replication, protein processing in the endoplasmic reticulum, cell cycle, apoptosis, p53 signaling pathway, and protein digestion and absorption. In the comparison of the TH group vs. control group and TH_CHX group vs. CHX group overlapping gene, 39 upregulated and 6 downregulated genes were identified as the TH directly induced genes. Further analysis indicated that AGGTCAnnTnAGGTCA is the optimal target sequence of target genes for TR/RXR heterodimers in M. fissipes . Future investigations on the function and regulation of these genes and pathways should help to reveal the mechanisms governing amphibian dorsal muscle remodeling. These full-length and high-quality transcriptomes in this study also provide an important foundation for future studies in M. fissipes metamorphosis.