Transcriptional activation of candidate paused genes after partial hepatectomy.

Partial hepatectomy is a model of acute liver injury that is known to induce a strong reprogrammation of gene expression. Transcriptional induction of Immediate Early Genes is extremely fast and this would be due to the release of RNA Polymerase II poised for elongation at 'paused' genes. Using bioinformatic analysis, we identified 23 genes sharing features of paused genes before hepatectomy, and with predicted quick and strong expression induction after. This transcriptional dynamic, confirmed by RT-qPCR for Jun , Fos , Btg2, is very precocious. RNA Pol II CTD Ser2 hyperphosphorylation indicates a switch to productive elongation and release from transcriptional pause.

Xenopus tropicalis osteoblast-specific open chromatin regions reveal promoters and enhancers involved in human skeletal phenotypes and shed light on early vertebrate evolution.

While understanding the genetic underpinnings of osteogenesis has far-reaching implications for skeletal diseases and evolution, a comprehensive characterization of the osteoblastic regulatory landscape in non-mammalian vertebrates is still lacking. Here, we compared the ATAC-Seq profile of Xenopus tropicalis (Xt) osteoblasts to a variety of non mineralizing control tissues, and identified osteoblast-specific nucleosome free regions (NFRs) at 527 promoters and 6747 distal regions. Sequence analyses, Gene Ontology, RNA-Seq and ChIP-Seq against four key histone marks confirmed that the distal regions correspond to bona fide osteogenic transcriptional enhancers exhibiting a shared regulatory logic with mammals. We report 425 regulatory regions conserved with human and globally associated to skeletogenic genes. Of these, 35 regions have been shown to impact human skeletal phenotypes by GWAS, including one trps1 enhancer and the runx2 promoter, two genes which are respectively involved in trichorhinophalangeal syndrome type I and cleidocranial dysplasia. Intriguingly, 60 osteoblastic NFRs also align to the genome of the elephant shark, a species lacking osteoblasts and bone tissue. To tackle this paradox, we chose to focus on dlx5 because its conserved promoter, known to integrate regulatory inputs during mammalian osteogenesis, harbours an osteoblast-specific NFR in both frog and human. Hence, we show that dlx5 is expressed in Xt and elephant shark odontoblasts, supporting a common cellular and genetic origin of bone and dentine. Taken together, our work (i) unravels the Xt osteogenic regulatory landscape, (ii) illustrates how cross-species comparisons harvest data relevant to human biology and (iii) reveals that a set of genes including bnc2, dlx5, ebf3, mir199a, nfia, runx2 and zfhx4 drove the development of a primitive form of mineralized skeletal tissue deep in the vertebrate lineage.

Overlapping action of T3 and T4 during Xenopus laevis development.

Thyroid hormones are involved in many biological processes such as neurogenesis, metabolism, and development. However, compounds called endocrine disruptors can alter thyroid hormone signaling and induce unwanted effects on human and ecosystems health. Regulatory tests have been developed to detect these compounds but need to be significantly improved by proposing novel endpoints and key events. The Xenopus Eleutheroembryonic Thyroid Assay (XETA, OECD test guideline no. 248) is one such test. It is based on Xenopus laevis tadpoles, a particularly sensitive model system for studying the physiology and disruption of thyroid hormone signaling: amphibian metamorphosis is a spectacular (thus easy to monitor) life cycle transition governed by thyroid hormones. With a long-term objective of providing novel molecular markers under XETA settings, we propose first to describe the differential effects of thyroid hormones on gene expression, which, surprisingly, are not known. After thyroid hormones exposure (T3 or T4), whole tadpole RNAs were subjected to transcriptomic analysis. By using standard approaches coupled to system biology, we found similar effects of the two thyroid hormones. They impact the cell cycle and promote the expression of genes involves in cell proliferation. At the level of the whole tadpole, the immune system is also a prime target of thyroid hormone action.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs.

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Mitochondria Transfer from Mesenchymal Stem Cells Confers Chemoresistance to Glioblastoma Stem Cells through Metabolic Rewiring.

Glioblastomas (GBM) are heterogeneous tumors with high metabolic plasticity. Their poor prognosis is linked to the presence of glioblastoma stem cells (GSC), which support resistance to therapy, notably to temozolomide (TMZ). Mesenchymal stem cells (MSC) recruitment to GBM contributes to GSC chemoresistance, by mechanisms still poorly understood. Here, we provide evidence that MSCs transfer mitochondria to GSCs through tunneling nanotubes, which enhances GSCs resistance to TMZ. More precisely, our metabolomics analyses reveal that MSC mitochondria induce GSCs metabolic reprograming, with a nutrient shift from glucose to glutamine, a rewiring of the tricarboxylic acid cycle from glutaminolysis to reductive carboxylation and increase in orotate turnover as well as in pyrimidine and purine synthesis. Metabolomics analysis of GBM patient tissues at relapse after TMZ treatment documents increased concentrations of AMP, CMP, GMP, and UMP nucleotides and thus corroborate our in vitro analyses. Finally, we provide a mechanism whereby mitochondrial transfer from MSCs to GSCs contributes to GBM resistance to TMZ therapy, by demonstrating that inhibition of orotate production by Brequinar (BRQ) restores TMZ sensitivity in GSCs with acquired mitochondria. Altogether, these results identify a mechanism for GBM resistance to TMZ and reveal a metabolic dependency of chemoresistant GBM following the acquisition of exogenous mitochondria, which opens therapeutic perspectives based on synthetic lethality between TMZ and BRQ. Significance: Mitochondria acquired from MSCs enhance the chemoresistance of GBMs. The discovery that they also generate metabolic vulnerability in GSCs paves the way for novel therapeutic approaches.

Crosstalk between Thyroid Hormone and Corticosteroid Signaling Targets Cell Proliferation in Xenopus tropicalis Tadpole Liver.

Thyroid hormones (TH) and glucocorticoids (GC) are involved in numerous developmental and physiological processes. The effects of individual hormones are well documented, but little is known about the joint actions of the two hormones. To decipher the crosstalk between these two hormonal pathways, we conducted a transcriptional analysis of genes regulated by TH, GC, or both hormones together in liver of Xenopus tropicalis tadpoles using RNA-Seq. Among the differentially expressed genes (DE), 70.5% were regulated by TH only, 0.87% by GC only, and 15% by crosstalk between the two hormones. Gene ontology analysis of the crosstalk-regulated genes identified terms referring to DNA replication, DNA repair, and cell-cycle regulation. Biological network analysis identified groups of genes targeted by the hormonal crosstalk and corroborated the gene ontology analysis. Specifically, we found two groups of functionally linked genes (chains) mainly composed of crosstalk-regulated hubs (highly interactive genes), and a large subnetwork centred around the crosstalk-regulated genes psmb6 and cdc7. Most of the genes in the chains are involved in cell-cycle regulation, as are psmb6 and cdc7, which regulate the G2/M transition. Thus, the biological action of these two hormonal pathways acting together in the liver targets cell-cycle regulation.

Tetrabromobisphenol A effects on differentiating mouse embryonic stem cells reveals unexpected impact on immune system.

Tetrabromobisphenol A (TBBPA) is a potent flame retardant used in numerous appliances and a major pollutant in households and ecosystems. In vertebrates, it was shown to affect neurodevelopment, the hypothalamic-pituitary-gonadal axis and thyroid signaling, but its toxicity and modes of actions are still a matter of debate. The molecular phenotype resulting from exposure to TBBPA is only poorly described, especially at the level of transcriptome reprogramming, which further limits our understanding of its molecular toxicity. In this work, we combined functional genomics and system biology to provide a system-wide description of the transcriptomic alterations induced by TBBPA acting on differentiating mESCs, and provide potential new toxicity markers. We found that TBBPA-induced transcriptome reprogramming affect a large collection of genes loosely connected within the network of biological pathways, indicating widespread interferences on biological processes. We also found two hotspots of action: at the level of neuronal differentiation markers, and surprisingly, at the level of immune system functions, which has been largely overlooked until now. This effect is particularly strong, as terminal differentiation markers of both myeloid and lymphoid lineages are strongly reduced: the membrane T cell receptor (Cd79a, Cd79b), interleukin seven receptor (Il7r), macrophages cytokine receptor (Csf1r), monocyte chemokine receptor (Ccr2). Also, the high affinity IgE receptor (Fcer1g), a key mediator of allergic reactions, is strongly induced. Thus, the molecular imbalance induce by TBBPA may be stronger than initially realized.

Disruption by stealth - Interference of endocrine disrupting chemicals on hormonal crosstalk with thyroid axis function in humans and other animals.

Thyroid hormones (THs) are important regulators of growth, development, and homeostasis of all vertebrates. There are many environmental contaminants that are known to disrupt TH action, yet their mechanisms are only partially understood. While the effects of Endocrine Disrupting Chemicals (EDCs) are mostly studied as "hormone system silos", the present critical review highlights the complexity of EDCs interfering with TH function through their interactions with other hormonal axes involved in reproduction, stress, and energy metabolism. The impact of EDCs on components that are shared between hormone signaling pathways or intersect between pathways can thus extend beyond the molecular ramifications to cellular, physiological, behavioral, and whole-body consequences for exposed organisms. The comparatively more extensive studies conducted in mammalian models provides encouraging support for expanded investigation and highlight the paucity of data generated in other non-mammalian vertebrate classes. As greater genomics-based resources become available across vertebrate classes, better identification and delineation of EDC effects, modes of action, and identification of effective biomarkers suitable for HPT disruption is possible. EDC-derived effects are likely to cascade into a plurality of physiological effects far more complex than the few variables tested within any research studies. The field should move towards understanding a system of hormonal systems' interactions rather than maintaining hormone system silos.

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 T3 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.

An alternative D. melanogaster 7SK snRNP.

BACKGROUND: The 7SK small nuclear RNA (snRNA) found in most metazoans is a key regulator of P-TEFb which in turn regulates RNA polymerase II elongation. Although its primary sequence varies in protostomes, its secondary structure and function are conserved across evolutionary distant taxa. RESULTS: Here, we describe a novel ncRNA sharing many features characteristic of 7SK RNAs, in D. melanogaster. We examined the structure of the corresponding gene and determined the expression profiles of the encoded RNA, called snRNA:7SK:94F, during development. It is probably produced from the transcription of a lncRNA which is processed into a mature snRNA. We also addressed its biological function and we show that, like dm7SK, this alternative 7SK interacts in vivo with the different partners of the P-TEFb complex, i.e. HEXIM, LARP7 and Cyclin T. This novel RNA is widely expressed across tissues. CONCLUSION: We propose that two distinct 7SK genes might contribute to the formation of the 7SK snRNP complex in D. melanogaster.

Two repeated motifs enriched within some enhancers and origins of replication are bound by SETMAR isoforms in human colon cells.

Setmar is a gene specific to simian genomes. The function(s) of its isoforms are poorly understood and their existence in healthy tissues remains to be validated. Here we profiled SETMAR expression and its genome-wide binding landscape in colon tissue. We found isoforms V3 and V6 in healthy and tumour colon tissues as well as incell lines. In two colorectal cell lines SETMAR binds to several thousand Hsmar1 and MADE1 terminal ends, transposons mostly located in non-genic regions of active chromatin including in enhancers. It also binds to a 12-bp motifs similar to an inner motif in Hsmar1 and MADE1 terminal ends. This motif is interspersed throughout the genome and is enriched in GC-rich regions as well as in CpG islands that contain constitutive replication origins. It is also found in enhancers other than those associated with Hsmar1 and MADE1. The role of SETMAR in the expression of genes, DNA replication and in DNA repair are discussed.

Are paedomorphs actual larvae?

Amphibians display very diverse life cycles and development can be direct, where it occurs in ovo and a juvenile hatches directly, or biphasic, where an aquatic larva hatches and later undergoes metamorphosis followed by sexual maturation. In both cases, metamorphosis, corresponds to the post embryonic transition (PETr). A third strategy, only found in Urodeles, is more complex as larvae reach sexual maturity before metamorphosis, which can become accessory. The resulting paedomorphs retain their larval characters and keep their aquatic habitat. Does it mean that paedomorphs do not undergo PETr? Recent work using high throughput technologies coupled to system biology and developmental endocrinology revisited this question and provided novel datasets indicating that a paedomorph's "larval" tissue undergoes a proper developmental transition. Together with historical data, we propose that this transition is a marker of the PETr, which would be distinct from metamorphosis. This implies that (a) complex life cycles would result from the uncoupling of PETr and metamorphosis, and (b) biphasic life cycles would be a special cases where they occur simultaneously.

The C-terminal Domain of piggyBac Transposase Is Not Required for DNA Transposition.

PiggyBac(PB)-like elements (pble) are members of a eukaryotic DNA transposon family. This family is of interest to evolutionary genomics because pble transposases have been domesticated at least 9 times in vertebrates. The amino acid sequence of pble transposases can be split into three regions: an acidic N-terminal domain (~100 aa), a central domain (~400 aa) containing a DD[D/E] catalytic triad, and a cysteine-rich domain (CRD; ~90 aa). Two recent reports suggested that a functional CRD is required for pble transposase activity. Here we found that two CRD-deficient pble transposases, a PB variant and an isoform encoded by the domesticated PB-derived vertebrate transposase gene 5 (pgbd5) trigger transposition of the Ifp2 pble. When overexpressed in HeLa cells, these CRD-deficient transposases can insert Ifp2 elements with proper and improper transposon ends, associated with deleterious effects on cells. Finally, we found that mouse CRD-deficient transposase Pgbd5, as well as PB, do not insert pbles at random into chromosomes. Transposition events occurred more often in genic regions, in the neighbourhood of the transcription start sites and were often found in genes predominantly expressed in the human central nervous system.

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.

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.

De Novo Transcriptomic Approach to Study Thyroid Hormone Receptor Action in Non-mammalian Models.

Thyroid hormones are pleiotropic hormones involved in chordates physiology. Understanding their functions and mechanisms is also instrumental to diagnose dys-regulations and get a predictive power that can applied to medicine, ecology, etc. Today, high-throughput sequencing technologies offer the opportunity to address this issue not only in model organisms but also in non-model organisms. Here, we describe a method that makes use of RNA-seq to address differential expression analysis in non-model organism.

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.