Cooperative, synergistic gene regulation by nuclear hormone receptors can increase sensitivity and amplify cellular responses to hormones. We investigated thyroid hormone (TH) and glucocorticoid (GC) synergy on the Krüppel-like factor 9 (Klf9) gene, which codes for a zinc finger transcription factor involved in development and homeostasis of diverse tissues. We identified regions of the Xenopus and mouse Klf9 genes 5-6 kb upstream of the transcription start sites that supported synergistic transactivation by TH plus GC. Within these regions, we found an orthologous sequence of approximately 180 bp that is highly conserved among tetrapods, but absent in other chordates, and possesses chromatin marks characteristic of an enhancer element. The Xenopus and mouse approximately 180-bp DNA element conferred synergistic transactivation by hormones in transient transfection assays, so we designate this the Klf9 synergy module (KSM). We identified binding sites within the mouse KSM for TH receptor, GC receptor, and nuclear factor κB. TH strongly increased recruitment of liganded GC receptor and serine 5 phosphorylated (initiating) RNA polymerase II to chromatin at the KSM, suggesting a mechanism for transcriptional synergy. The KSM is transcribed to generate long noncoding RNAs, which are also synergistically induced by combined hormone treatment, and the KSM interacts with the Klf9 promoter and a far upstream region through chromosomal looping. Our findings support that the KSM plays a central role in hormone regulation of vertebrate Klf9 genes, it evolved in the tetrapod lineage, and has been maintained by strong stabilizing selection.

Methylation of cytosine residues in DNA influences chromatin structure and gene transcription, and its regulation is crucial for brain development. There is mounting evidence that DNA methylation can be modulated by hormone signaling. We analyzed genome-wide changes in DNA methylation and their relationship to gene regulation in the brain of Xenopus tadpoles during metamorphosis, a thyroid hormone-dependent developmental process. We studied the region of the tadpole brain containing neurosecretory neurons that control pituitary hormone secretion, a region that is highly responsive to thyroid hormone action. Using Methylated DNA Capture sequencing (MethylCap-seq) we discovered a diverse landscape of DNA methylation across the tadpole neural cell genome, and pairwise stage comparisons identified several thousand differentially methylated regions (DMRs). During the pre-to pro-metamorphic period, the number of DMRs was lowest (1,163), with demethylation predominating. From pre-metamorphosis to metamorphic climax DMRs nearly doubled (2,204), with methylation predominating. The largest changes in DNA methylation were seen from metamorphic climax to the completion of metamorphosis (2960 DMRs), with 80% of the DMRs representing demethylation. Using RNA sequencing, we found negative correlations between differentially expressed genes and DMRs localized to gene bodies and regions upstream of transcription start sites. DNA demethylation at metamorphosis revealed by MethylCap-seq was corroborated by increased immunoreactivity for the DNA demethylation intermediates 5-hydroxymethylcytosine and 5-carboxymethylcytosine, and the methylcytosine dioxygenase ten eleven translocation 3 that catalyzes DNA demethylation. Our findings show that the genome of tadpole neural cells undergoes significant changes in DNA methylation during metamorphosis, and these changes likely influence chromatin architecture, and gene regulation programs occurring during this developmental period.

Thyroid hormone (T3) plays pivotal roles in vertebrate development, acting via nuclear receptors (TRs) that regulate gene transcription by promoting posttranslational modifications to histones. Methylation of cytosine residues in DNA also modulates gene transcription, and our recent finding of predominant DNA demethylation in the brain of Xenopus tadpoles at metamorphosis, a T3-dependent developmental process, caused us to hypothesize that T3 induces these changes in vivo. Treatment of pre-metamorphic tadpoles with T3 for 24 or 48 hr increased immunoreactivity in several brain regions for the DNA demethylation intermediates 5-hydroxymethylcytosine (5-hmC) and 5-carboxylcytosine, and the methylcytosine dioxygenase ten-eleven translocation 3 (TET3). Thyroid hormone treatment induced locus-specific DNA demethylation in proximity to known T3 response elements within the DNA methyltransferase 3a and Krüppel-like factor 9 genes, analyzed by 5-hmC immunoprecipitation and methylation sensitive restriction enzyme digest. Chromatin-immunoprecipitation (ChIP) assay showed that T3 induced TET3 recruitment to these loci. Furthermore, the mRNAs for several genes encoding DNA demethylation enzymes were induced by T3 in a time-dependent manner in tadpole brain. A TR ChIP-sequencing experiment identified putative TR binding sites at several of these genes, and we provide multiple lines of evidence to support that tet2 contains a bona fide T3 response element. Our findings show that T3 can promote DNA demethylation in developing tadpole brain, in part by promoting TET3 recruitment to discrete genomic regions, and by inducing genes that encode DNA demethylation enzymes.

Amphibian metamorphosis is controlled by thyroid hormone (TH), which binds TH receptors (TRs) to regulate gene expression programs that underlie morphogenesis. Gene expression screens using tissues from premetamorphic tadpoles treated with TH identified some TH target genes, but few studies have analyzed genome-wide changes in gene regulation during spontaneous metamorphosis. We analyzed RNA sequencing data at four developmental stages from the beginning to the end of spontaneous metamorphosis, conducted on the neuroendocrine centers of Xenopus tropicalis tadpole brain. We also conducted chromatin immunoprecipitation sequencing (ChIP-seq) for TRs, and we compared gene expression changes during metamorphosis with those induced by exogenous TH. The mRNA levels of 26% of protein coding genes changed during metamorphosis; about half were upregulated and half downregulated. Twenty four percent of genes whose mRNA levels changed during metamorphosis had TR ChIP-seq peaks. Genes involved with neural cell differentiation, cell physiology, synaptogenesis and cell-cell signaling were upregulated, while genes involved with cell cycle, protein synthesis, and neural stem/progenitor cell homeostasis were downregulated. There is a shift from building neural structures early in the metamorphic process, to the differentiation and maturation of neural cells and neural signaling pathways characteristic of the adult frog brain. Only half of the genes modulated by treatment of premetamorphic tadpoles with TH for 16 h changed expression during metamorphosis; these represented 33% of the genes whose mRNA levels changed during metamorphosis. Taken together, our results provide a foundation for understanding the molecular basis for metamorphosis of tadpole brain, and they highlight potential caveats for interpreting gene regulation changes in premetamorphic tadpoles induced by exogenous TH.