Xenopus tropicalis Genome Re-Scaffolding and Re-Annotation Reach the Resolution Required for In Vivo ChIA-PET Analysis

Transcriptome and Methylome Analysis Reveal Complex Cross-Talks between Thyroid Hormone and Glucocorticoid Signaling at Xenopus Metamorphosis

Background: Most work in endocrinology focus on the action of a single hormone, and very little on the cross-talks between two hormones. Here we characterize the nature of interactions between thyroid hormone and glucocorticoid signaling during Xenopus tropicalis metamorphosis. Methods: We used functional genomics to derive genome wide profiles of methylated DNA and measured changes of gene expression after hormonal treatments of a highly responsive tissue, tailfin. Clustering classified the data into four types of biological responses, and biological networks were modeled by system biology. Results: We found that gene expression is mostly regulated by either T₃ or CORT, or their additive effect when they both regulate the same genes. A small but non-negligible fraction of genes (12%) displayed non-trivial regulations indicative of complex interactions between the signaling pathways. Strikingly, DNA methylation changes display the opposite and are dominated by cross-talks. Conclusion: Cross-talks between thyroid hormones and glucocorticoids are more complex than initially envisioned and are not limited to the simple addition of their individual effects, a statement that can be summarized with the pseudo-equation: TH ∙ GC > TH + GC. DNA methylation changes are highly dynamic and buffered from genome expression.

DNA methylation dynamics underlie metamorphic gene regulation programs in Xenopus tadpole brain

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.

A Pivotal Genetic Program Controlled by Thyroid Hormone during the Maturation of GABAergic Neurons

Mammalian brain development critically depends on proper thyroid hormone signaling, via the TRα1 nuclear receptor. The downstream mechanisms by which TRα1 impacts brain development are currently unknown. In order to investigate these mechanisms, we used mouse genetics to induce the expression of a dominant-negative mutation of TRα1 specifically in GABAergic neurons, the main inhibitory neurons in the brain. This triggered post-natal epileptic seizures and a profound impairment of GABAergic neuron maturation in several brain regions. Analysis of the transcriptome and TRα1 cistrome in the striatum allowed us to identify a small set of genes, the transcription of which is upregulated by TRα1 in GABAergic neurons and which probably plays an important role during post-natal maturation of the brain. Thus, our results point to GABAergic neurons as direct targets of thyroid hormone during brain development and suggest that many defects seen in hypothyroid brains may be secondary to GABAergic neuron malfunction.

Contaminant and Environmental Influences on Thyroid Hormone Action in Amphibian Metamorphosis

Aquatic and terrestrial environments are increasingly contaminated by anthropogenic sources that include pharmaceuticals, personal care products, and industrial and agricultural chemicals (i. e., pesticides). Many of these substances have the potential to disrupt endocrine function, yet their effect on thyroid hormone (TH) action has garnered relatively little attention. Anuran postembryonic metamorphosis is strictly dependent on TH and perturbation of this process can serve as a sensitive barometer for the detection and mechanistic elucidation of TH disrupting activities of chemical contaminants and their complex mixtures. The ecological threats posed by these contaminants are further exacerbated by changing environmental conditions such as temperature, photoperiod, pond drying, food restriction, and ultraviolet radiation. We review the current knowledge of several chemical and environmental factors that disrupt TH-dependent metamorphosis in amphibian tadpoles as assessed by morphological, thyroid histology, behavioral, and molecular endpoints. Although the molecular mechanisms for TH disruption have yet to be determined for many chemical and environmental factors, several affect TH synthesis, transport or metabolism with subsequent downstream effects. As molecular dysfunction typically precedes phenotypic or histological pathologies, sensitive assays that detect changes in transcript, protein, or metabolite abundance are indispensable for the timely detection of TH disruption. The emergence and application of 'omics techniques-genomics, transcriptomics, proteomics, metabolomics, and epigenomics-on metamorphosing tadpoles are powerful emerging assets for the rapid, proxy assessment of toxicant or environmental damage for all vertebrates including humans. Moreover, these highly informative 'omics techniques will complement morphological, behavioral, and histological assessments, thereby providing a comprehensive understanding of how TH-dependent signal disruption is propagated by environmental contaminants and factors.

Insufficiency of Thyroid Hormone in Frog Metamorphosis and the Role of Glucocorticoids

Thyroid hormone (TH) is the most important hormone in frog metamorphosis, a developmental process which will not occur in the absence of TH but can be induced precociously by exogenous TH. However, such treatments including in-vitro TH treatments often do not replicate the events of natural metamorphosis in many organs, including lung, brain, blood, intestine, pancreas, tail, and skin. A potential explanation for the discrepancy between natural and TH-induced metamorphosis is the involvement of glucocorticoids (GCs). GCs are not able to advance development by themselves but can modulate the rate of developmental progress induced by TH via increased tissue sensitivity to TH. Global gene expression analyses and endocrine experiments suggest that GCs may also have direct actions required for completion of metamorphosis independent of their effects on TH signaling. Here, we provide a new review and analysis of the requirement and necessity of TH signaling in light of recent insights from gene knockout frogs. We also examine the independent and interactive roles GCs play in regulating morphological and molecular metamorphic events dependent upon TH.

Opposite T3 Response of ACTG1-FOS Subnetwork Differentiate Tailfin Fate in Xenopus Tadpole and Post-hatching Axolotl

Amphibian post-embryonic development and Thyroid Hormones (TH) signaling are deeply and intimately connected. In anuran amphibians, TH induce the spectacular and complex process known as metamorphosis. In paedomorphic salamanders, at similar development time, raising levels of TH fail to induce proper metamorphosis, as many "larval" tissues (e.g., gills, tailfin) are maintained. Why does the same evolutionary conserved signaling pathway leads to alternative phenotypes? We used a combination of developmental endocrinology, functional genomics and network biology to compare the transcriptional response of tailfin to TH, in the post-hatching paedormorphic Axolotl salamander and Xenopus tadpoles. We also provide a technological framework that efficiently reduces large lists of regulated genes down to a few genes of interest, which is well-suited to dissect endocrine regulations. We first show that Axolotl tailfin undergoes a strong and robust TH-dependent transcriptional response at post embryonic transition, despite the lack of visible anatomical changes. We next show that Fos and Actg1, which structure a single and dense subnetwork of cellular sensors and regulators, display opposite regulation between the two species. We finally show that TH treatments and natural variations of TH levels follow similar transcriptional dynamics. We suggest that, at the molecular level, tailfin fate correlates with the alternative transcriptional states of an fos-actg1 sub-network, which also includes transcription factors and regulators of cell fate. We propose that this subnetwork is one of the molecular switches governing the initiation of distinct TH responses, with transcriptional programs conducting alternative tailfin fate (maintenance vs. resorption) 2 weeks post-hatching.

Methods for Investigating the Larval Period and Metamorphosis in Xenopus

Anuran metamorphosis resembles postembryonic development in mammals, a period around birth when many organs/tissues mature into their adult form as circulating thyroid and stress hormone levels are high. Unlike uterus-enclosed mammalian embryos, tadpoles develop externally and undergo the dramatic changes of hormone-dependent development totally independent of maternal influence, making them a valuable model in which to study vertebrate postembryonic organ development and maturation. Various protocols have been developed and/or adapted for studying metamorphosis in Xenopus laevis and X. tropicalis, two highly related and well-studied frog species. Here, we introduce some of the methods for contemporary cell and molecular studies of gene function and regulation during metamorphosis.

Chromatin Immunoprecipitation for Chromatin Interaction Analysis Using Paired-End-Tag (ChIA-PET) Sequencing in Tadpole Tissues

Proper gene expression involves communication between the regulatory elements and promoters of genes. Because regulatory elements can be located over a large range of genomic distances (from as close as a few hundred bp to as much as several Mb away), contact and communication between regulators and the core transcriptional machinery at promoters are mediated through DNA looping. Today, chromosome conformation capture (3C)-based methods efficiently probe chromosome folding in the nucleus and thus provide a molecular description of physical proximity between enhancer(s) and their target promoter(s). One such method, chromatin interaction analysis using paired-end-tag (ChIA-PET) sequencing, is a leading high-throughput method for detection of genome wide chromatin interactions. Briefly, the method involves cross-linkage of chromatin (-DNA) fibers in cells in situ, fragmentation of the fixed chromatin-DNA complexes by sonication, followed by enrichment of the chromatin complexes with a dedicated antibody through the process of immunoprecipitation (IP). Next, application of the ChIA-PET protocol followed by deep sequencing and mapping of reads to the reference genome reveals both binding sites and remote chromatin interactions mediated by the protein factors of interest. The method detailed here focuses on ChIP sample preparation and can be completed in ∼5 d. The ChIA-PET method is detailed in an associated protocol. Because not all chromatin immunoprecipitation protocols are suitable for ChIA-PET, it is important to strictly follow this procedure before performing the ChIA-PET protocol.

Chromatin Interaction Analysis Using Paired-End-Tag (ChIA-PET) Sequencing in Tadpole Tissues

Proper gene expression involves communication between the regulatory elements and promoters of genes. Today, chromosome conformation capture (3C)-based methods efficiently probe chromosome folding in the nucleus and thus provide a molecular description of physical proximity through DNA looping between enhancer(s) and their target promoter(s). One such method, chromatin interaction analysis using paired-end-tag (ChIA-PET) sequencing is a powerful high-throughput method for detection of genome-wide chromatin interactions. Following enrichment of the chromatin complexes with a dedicated antibody, through a process of immunoprecipitation (IP), DNA fragments are end-joined with specifically designed DNA-linkers through proximity ligation. The DNA-linkers contain the binding site for the type II restriction enzyme MmeI, which cleaves 20 bp from each end of the ligated fragments, thus releasing a "paired end tag" (PET): [20 bp tag]-[linker]-[20 bp tag]. The PETs are then deep-sequenced and reads are mapped to the reference genome, revealing both binding sites, as well as remote chromatin interactions mediated by the protein factors of interest. The method detailed here focuses on ChIA-PET library construction and can be completed in 2 wk.

The North American bullfrog draft genome provides insight into hormonal regulation of long noncoding RNA

Frogs play important ecological roles, and several species are important model organisms for scientific research. The globally distributed Ranidae (true frogs) are the largest frog family, and have substantial evolutionary distance from the model laboratory Xenopus frog species. Unfortunately, there are currently no genomic resources for the former, important group of amphibians. More widely applicable amphibian genomic data is urgently needed as more than two-thirds of known species are currently threatened or are undergoing population declines. We report a 5.8 Gbp (NG50 = 69 kbp) genome assembly of a representative North American bullfrog (Rana [Lithobates] catesbeiana). The genome contains over 22,000 predicted protein-coding genes and 6,223 candidate long noncoding RNAs (lncRNAs). RNA-Seq experiments show thyroid hormone causes widespread transcriptional change among protein-coding and putative lncRNA genes. This initial bullfrog draft genome will serve as a key resource with broad utility including amphibian research, developmental biology, and environmental research.

Genome-wide identification of thyroid hormone receptor targets in the remodeling intestine during Xenopus tropicalis metamorphosis

Thyroid hormone (T3) affects development and metabolism in vertebrates. We have been studying intestinal remodeling during T3-dependent Xenopus metamorphosis as a model for organ maturation and formation of adult organ-specific stem cells during vertebrate postembryonic development, a period characterized by high levels of plasma T3. T3 is believed to affect development by regulating target gene transcription through T3 receptors (TRs). While many T3 response genes have been identified in different animal species, few have been shown to be direct target genes in vivo, especially during development. Here we generated a set of genomic microarray chips covering about 8000 bp flanking the predicted transcription start sites in Xenopus tropicalis for genome wide identification of TR binding sites. By using the intestine of premetamorphic tadpoles treated with or without T3 and for chromatin immunoprecipitation assays with these chips, we determined the genome-wide binding of TR in the control and T3-treated tadpole intestine. We further validated TR binding in vivo and analyzed the regulation of selected genes. We thus identified 278 candidate direct TR target genes. We further provided evidence that these genes are regulated by T3 and likely involved in the T3-induced formation of adult intestinal stem cells during metamorphosis.

Using Xenopus to understand human disease and developmental disorders

Model animals are crucial to biomedical research. Among the commonly used model animals, the amphibian, Xenopus, has had tremendous impact because of its unique experimental advantages, cost effectiveness, and close evolutionary relationship with mammals as a tetrapod. Over the past 50 years, the use of Xenopus has made possible many fundamental contributions to biomedicine, and it is a cornerstone of research in cell biology, developmental biology, evolutionary biology, immunology, molecular biology, neurobiology, and physiology. The prospects for Xenopus as an experimental system are excellent: Xenopus is uniquely well-suited for many contemporary approaches used to study fundamental biological and disease mechanisms. Moreover, recent advances in high throughput DNA sequencing, genome editing, proteomics, and pharmacological screening are easily applicable in Xenopus, enabling rapid functional genomics and human disease modeling at a systems level.

Frogs model man: In vivo thyroid hormone signaling during development

Thyroid hormone (TH) signaling comprises TH transport across cell membranes, metabolism by deiodinases, and molecular mechanisms of gene regulation. Proper TH signaling is essential for normal perinatal development, most notably for neurogenesis and fetal growth. Knowledge of perinatal TH endocrinology needs improvement to provide better treatments for premature infants and endocrine diseases during gestation and to counteract effects of endocrine disrupting chemicals. Studies in amphibians have provided major insights to understand in vivo mechanisms of TH signaling. The frog model boasts dramatic TH-dependent changes directly observable in free-living tadpoles with precise and easy experimental control of the TH response at developmental stages comparable to fetal stages in mammals. The hormones, their receptors, molecular mechanisms, and developmental roles of TH signaling are conserved to a high degree in humans and amphibians, such that with respect to developmental TH signaling "frogs are just little people that hop." The frog model is exceptionally illustrative of fundamental molecular mechanisms of in vivo TH action involving TH receptors, transcriptional cofactors, and chromatin remodeling. This review highlights the current need, recent successes, and future prospects using amphibians as a model to elucidate molecular mechanisms and functional roles of TH signaling during post-embryonic development.

Genetic Investigation of Thyroid Hormone Receptor Function in the Developing and Adult Brain

Thyroid hormones exert a broad influence on brain development and function, which has been extensively studied over the years. Mouse genetics has brought an important contribution, allowing precise analysis of the interplay between TRα1 and TRβ1 nuclear receptors in neural cells. However, the exact contribution of each receptor, the possible intervention of nongenomic signaling, and the nature of the genetic program that is controlled by the receptors remain poorly understood.

Systems Biology Analyses in Chicken: Workflow for Transcriptome and ChIP-Seq Analyses Using the Chicken Skin Paradigm

With advances in molecular biology, various biological phenomena can now be explored at higher resolution using mRNA sequencing (RNA-Seq) and chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-Seq), two powerful high-throughput next-generation sequencing (NGS) technologies. While methods are used widely in mouse, human, etc., less information is available in other animals, such as the chicken. Here we assemble a workflow of the RNA-Seq and ChIP-Seq analyses for the chicken studies using chicken skin appendage tissue as an example. We present guidelines for RNA-Seq quality control, alignment, quantification, normalization, and differentially expressed gene analysis. In the meantime, we outline a bioinformatics pipeline for ChIP-Seq quality control, alignment, peak calling, super-enhancer identification, and differential enrichment analysis.

Developmental and Thyroid Hormone Regulation of the DNA Methyltransferase 3a Gene in Xenopus Tadpoles

Thyroid hormone is essential for normal development in vertebrates. In amphibians, T3 controls metamorphosis by inducing tissue-specific gene regulation programs. A hallmark of T3 action is the modification of chromatin structure, which underlies changes in gene transcription. We found that mRNA for the de novo DNA methyltransferase (DNMT) dnmt3a, but not dnmt1, increased in the brain of Xenopus tadpoles during metamorphosis in parallel with plasma [T3]. Addition of T3 to the rearing water caused a time-dependent increase in dnmt3a mRNA in tadpole brain, tail, and hind limb. By analyzing data from a genome-wide analysis of T3 receptor (TR) binding in tadpole tail, we identified several putative T3 response elements (TREs) within the dnmt3a locus. Using in vitro DNA binding, transient transfection-reporter, and chromatin immunoprecipitation assays for TRs, we identified two functional TREs at -7.1 kb and +5.1 kb relative to the dnmt3a transcription start site. Sequence alignment showed that these TREs are conserved between two related frog species, X. laevis and X. tropicalis, but not with amniotes. Our previous findings showed that this gene is directly regulated by liganded TRs in mouse brain, and whereas the two mouse TREs are conserved among Eutherian mammals, they are not conserved in Xenopus species. Thus, although T3 regulation of dnmt3a may be an ancient pathway in vertebrates, the genomic sites responsible for hormone regulation may have diverged or arisen by convergent evolution. We hypothesize that direct T3 regulation of dnmt3a may be an important mechanism for modulating global changes in DNA methylation.

Futures Challenges in Thyroid Hormone Signaling Research

The canonical pathway of thyroid hormone signaling involves its binding to nuclear receptors (TRs) acting directly on the transcription of a number of genes. Recent genome-wide studies revealed that chromatin occupancy by TR is not sufficient for transactivation of gene expression. Reciprocally, in some cases, DNA binding by TR may not be required for cellular response. This leaves many new questions to be addressed in future research.