Intestinal remodeling during Xenopus metamorphosis as a model for studying thyroid hormone signaling and adult organogenesis

Intestinal development takes places in two phases, the initial formation of neonatal (mammals)/larval (anurans) intestine and its subsequent maturation into the adult form. This maturation occurs during postembryonic development when plasma thyroid hormone (T3) level peaks. In anurans such as the highly related Xenopus laevis and Xenopus tropicalis, the larval/tadpole intestine is drastically remodeled from a simple tubular structure to a complex, multi-folded adult organ during T3-dependent metamorphosis. This involved complete degeneration of larval epithelium via programmed cell death and de novo formation of adult epithelium, with concurrent maturation of the muscles and connective tissue. Here, we will summarize our current understanding of the underlying molecular mechanisms, with a focus on more recent genetic and genome-wide studies.

Revealing mitf functions and visualizing allografted tumor metastasis in colorless and immunodeficient Xenopus tropicalis

Transparent immunodeficient animal models not only enhance in vivo imaging investigations of visceral organ development but also facilitate in vivo tracking of transplanted tumor cells. However, at present, transparent and immunodeficient animal models are confined to zebrafish, presenting substantial challenges for real-time, in vivo imaging studies addressing specific biological inquiries. Here, we employed a mitf-/-/prkdc-/-/il2rg-/- triple-knockout strategy to establish a colorless and immunodeficient amphibian model of Xenopus tropicalis. By disrupting the mitf gene, we observed the loss of melanophores, xanthophores, and granular glands in Xenopus tropicalis. Through the endogenous mitf promoter to drive BRAFV600E expression, we confirmed mitf expression in melanophores, xanthophores and granular glands. Moreover, the reconstruction of the disrupted site effectively reinstated melanophores, xanthophores, and granular glands, further highlighting the crucial role of mitf as a regulator in their development. By crossing mitf-/- frogs with prkdc-/-/il2rg-/- frogs, we generated a mitf-/-/prkdc-/-/il2rg-/- Xenopus tropicalis line, providing a colorless and immunodeficient amphibian model. Utilizing this model, we successfully observed intravital metastases of allotransplanted xanthophoromas and migrations of allotransplanted melanomas. Overall, colorless and immunodeficient Xenopus tropicalis holds great promise as a valuable platform for tumorous and developmental biology research.

DNA methylation clocks for clawed frogs reveal evolutionary conservation of epigenetic aging

To address how conserved DNA methylation-based epigenetic aging is in diverse branches of the tree of life, we generated DNA methylation data from African clawed frogs (Xenopus laevis) and Western clawed frogs (Xenopus tropicalis) and built multiple epigenetic clocks. Dual species clocks were developed that apply to both humans and frogs (human-clawed frog clocks), supporting that epigenetic aging processes are evolutionary conserved outside mammals. Highly conserved positively age-related CpGs are located in neural-developmental genes such as uncx, tfap2d as well as nr4a2 implicated in age-associated disease. We conclude that signatures of epigenetic aging are evolutionary conserved between frogs and mammals and that the associated genes relate to neural processes, altogether opening opportunities to employ Xenopus as a model organism to study aging.

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.

Mutations in the histone methyltransferase Ezh2 drive context-dependent leukemia in Xenopus tropicalis

CRISPR-mediated simultaneous targeting of candidate tumor suppressor genes in Xenopus tropicalis allows fast functional assessment of co-driver genes for various solid tumors. Genotyping of tumors that emerge in the mosaic mutant animals rapidly exposes the gene mutations under positive selection for tumor establishment. However, applying this simple approach to the blood lineage has not been attempted. Multiple hematologic malignancies have mutations in EZH2, encoding the catalytic subunit of the Polycomb Repressive Complex 2. Interestingly, EZH2 can act as an oncogene or a tumor suppressor, depending on cellular context and disease stage. We show here that mosaic CRISPR/Cas9 mediated ezh2 disruption in the blood lineage resulted in early and penetrant acute myeloid leukemia (AML) induction. While animals were co-targeted with an sgRNA that induces notch1 gain-of-function mutations, sequencing of leukemias revealed positive selection towards biallelic ezh2 mutations regardless of notch1 mutational status. Co-targeting dnm2, recurrently mutated in T/ETP-ALL, induced a switch from myeloid towards acute T-cell leukemia. Both myeloid and T-cell leukemias engrafted in immunocompromised hosts. These data underline the potential of Xenopus tropicalis for modeling human leukemia, where mosaic gene disruption, combined with deep amplicon sequencing of the targeted genomic regions, can rapidly and efficiently expose co-operating driver gene mutations.

Age-associated DNA methylation changes in Xenopus frogs

Age-associated changes in DNA methylation have been characterized across various animals, but not yet in amphibians, which are of particular interest because they include widely studied model organisms. In this study, we present clear evidence that the aquatic vertebrate species Xenopus tropicalis displays patterns of age-associated changes in DNA methylation. We have generated whole-genome bisulfite sequencing (WGBS) profiles from skin samples of nine frogs representing young, mature, and old adults and characterized the gene- and chromosome-scale DNA methylation changes with age. Many of the methylation features and changes we observe are consistent with what is known in mammalian species, suggesting that the mechanism of age-related changes is conserved. Moreover, we selected a few thousand age-associated CpG sites to build an assay based on targeted DNA methylation analysis (TBSseq) to expand our findings in future studies involving larger cohorts of individuals. Preliminary results of a pilot TBSeq experiment recapitulate the findings obtained with WGBS setting the basis for the development of an epigenetic clock assay. The results of this study will allow us to leverage the unique resources available for Xenopus to study how DNA methylation relates to other hallmarks of ageing.

Phenotype-genotype relationships in Xenopus sox9 crispants provide insights into campomelic dysplasia and vertebrate jaw evolution

Since CRISPR-based genome editing technology works effectively in the diploid frog Xenopus tropicalis, a growing number of studies have successfully modeled human genetic diseases in this species. However, most of their targets were limited to non-syndromic diseases that exhibit abnormalities in a small fraction of tissues or organs in the body. This is likely because of the complexity of interpreting the phenotypic variations resulting from somatic mosaic mutations generated in the founder animals (crispants). In this study, we attempted to model the syndromic disease campomelic dysplasia (CD) by generating sox9 crispants in X. tropicalis. The resulting crispants failed to form neural crest cells at neurula stages and exhibited various combinations of jaw, gill, ear, heart, and gut defects at tadpole stages, recapitulating part of the syndromic phenotype of CD patients. Genotyping of the crispants with a variety of allelic series of mutations suggested that the heart and gut defects depend primarily on frame-shift mutations expected to be null, whereas the jaw, gill, and ear defects could be induced not only by such mutations but also by in-frame deletion mutations expected to delete part of the jawed vertebrate-specific domain from the encoded Sox9 protein. These results demonstrate that Xenopus crispants are useful for investigating the phenotype-genotype relationships behind syndromic diseases and examining the tissue-specific role of each functional domain within a single protein, providing novel insights into vertebrate jaw evolution.

Expanding EMC foldopathies: Topogenesis deficits alter the neural crest

The endoplasmic reticulum (ER) membrane protein complex (EMC) is essential for the insertion of a wide variety of transmembrane proteins into the plasma membrane across cell types. Each EMC is composed of Emc1-7, Emc10, and either Emc8 or Emc9. Recent human genetics studies have implicated variants in EMC genes as the basis for a group of human congenital diseases. The patient phenotypes are varied but appear to affect a subset of tissues more prominently than others. Namely, craniofacial development seems to be commonly affected. We previously developed an array of assays in Xenopus tropicalis to assess the effects of emc1 depletion on the neural crest, craniofacial cartilage, and neuromuscular function. We sought to extend this approach to additional EMC components identified in patients with congenital malformations. Through this approach, we determine that EMC9 and EMC10 are important for neural crest development and the development of craniofacial structures. The phenotypes observed in patients and our Xenopus model phenotypes similar to EMC1 loss of function likely due to a similar mechanism of dysfunction in transmembrane protein topogenesis.

Identification of tumor-related genes via RNA sequencing of tumor tissues in Xenopus tropicalis

Cancer treatment is still challenging because the disease is often caused by multiple mutations. Although genomic studies have identified many oncogenes and tumor suppressor genes, gene sets involved in tumorigenesis remain poorly understood. Xenopus, a genus of aquatic frogs, is a useful model to identify gene sets because it can be genetically and experimentally analyzed. Here, we analyzed gene expression in tumor tissues of three individuals in Xenopus tropicalis and identified 55 differentially expressed genes (DEGs). Gene ontology (GO) analysis showed that the upregulated genes in the tumor tissues were enriched in GO terms related to the extracellular matrix and collagen fibril organization. Hierarchical clustering showed that the gene expression patterns of tumor tissues in X. tropicalis were comparable to those of human connective, soft, and subcutaneous tissue-derived cancers. Additionally, pathway analysis revealed that these DEGs were associated with multiple pathways, including the extracellular matrix, collagen fibril organization, MET signaling, and keratan sulfate. We also found that the expression tendency of some DEGs that have not been well analyzed in the cancer field clearly determines the prognosis of human cancer patients. This study provides a remarkable reference for future experimental work on X. tropicalis to identify gene sets involved in human cancer.

Evolutionary and functional conservation of myeloid differentiation factor 88 (MyD88) in amphibian Xenopus tropicalis

As a universal adaptor used by most TLR members, the myeloid differentiation factor 88 (MyD88) plays essential roles in TLR-mediated inflammatory response of invertebrate and vertebrate animals, and functional features of MyD88 remain largely unknown in amphibians. In this study, a MyD88 gene named Xt-MyD88 was characterized in the Western clawed frog (Xenopus tropicalis). Xt-MyD88 and MyD88 in other species of vertebrates share similar structural characteristics, genomic structures, and flanking genes, suggesting that MyD88 is structurally conserved in different phyla of vertebrates ranging from fish to mammals. Moreover, Xt-MyD88 was widely expressed in different organs/tissues, and was induced by poly(I:C) in spleen, kidney, and liver. Importantly, overexpression of Xt-MyD88 triggered a marked activation of both NF-κB promoter and interferon-stimulated response elements (ISREs), implying that it may be play important roles in inflammatory responses of amphibians. The research represents the first characterization on the immune functions of amphibian MyD88, and reveals considerable functional conservation of MyD88 in early tetrapods.

Modification of maternally defined H3K4me3 regulates the inviability of interspecific Xenopus hybrids

Increasing evidence suggests that interspecific hybridization is crucial to speciation. However, chromatin incompatibility during interspecific hybridization often renders this process. Genomic imbalances such as chromosomal DNA loss and rearrangements leading to infertility have been commonly noted in hybrids. The mechanism underlying reproductive isolation of interspecific hybridization remains elusive. Here, we identified that modification of maternally defined H3K4me3 in Xenopus laevis and Xenopus tropicalis hybrids determines the different fates of the two types of hybrids as te×ls with developmental arrest and viable le×ts. Transcriptomics highlighted that the P53 pathway was overactivated, and the Wnt signaling pathway was suppressed in te×ls hybrids. Moreover, the lack of maternal H3K4me3 in te×ls disturbed the balance of gene expression between the L and S subgenomes in this hybrid. Attenuation of p53 can postpone the arrested development of te×ls. Our study suggests an additional model of reproductive isolation based on modifications of maternally defined H3K4me3.

Characterization of a novel corticosterone response gene in Xenopus tropicalis tadpole tails

Corticosteroids are critical for development and for mediating stress responses across diverse vertebrate taxa. Study of frog metamorphosis has made significant breakthroughs in our understanding of corticosteroid signaling during development in non-mammalian vertebrate species. However, lack of adequate corticosterone (CORT) response genes in tadpoles make identification and quantification of CORT responses challenging. Here, we characterized a CORT-response gene frzb (frizzled related protein) previously identified in Xenopus tropicalis tadpole tail skin by an RNA-seq study. We validated the RNA-seq results that CORT and not thyroid hormone induces frzb in the tails using quantitative PCR. Further, maximum frzb expression was achieved by 100-250 nM CORT within 12-24 hours. frzb is not significantly induced in the liver and brain in response to 100 nM CORT. We also found no change in frzb expression across natural metamorphosis when endogenous CORT levels peak. Surprisingly, frzb is only induced by CORT in X. tropicalis tails and not in Xenopus laevis tails. The exact downstream function of increased frzb expression in tails in response to CORT is not known, but the specificity of hormone response and its high mRNA expression levels in the tail render frzb a useful marker of exogenous CORT-response independent of thyroid hormone for exogenous hormone treatments and in-vivo endocrine disruption studies.

An Efficient CRISPR-Cas9 Method to Knock Out MiRNA Expression in Xenopus Tropicalis

In recent years CRISPR-Cas9 knockouts (KO) have become increasingly utilized to study gene function. MicroRNAs (miRNAs) are short noncoding RNAs, 20-25 nucleotides long, which affect gene expression through posttranscriptional repression. As miRNAs are so small and due to the limitations of known PAM sequences, it is difficult to design CRISPR sgRNAs that reproducibly lead to a KO. We have therefore developed a novel approach using two guide RNAs to effectively "drop out" a miRNA. Validation of efficient CRISPR miRNA KO and phenotype analysis included use of q-RT-PCR and Sanger sequencing. To show specificity of the phenotype, we provide a protocol to use miRNA mimics to rescue the KO phenotype.

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.

Upregulation of proto-oncogene ski by thyroid hormone in the intestine and tail during Xenopus metamorphosis

Thyroid hormone (T3) is important for adult organ function and vertebrate development, particularly during the postembryonic period when many organs develop/mature into their adult forms. Amphibian metamorphosis is totally dependent on T3 and can be easily manipulated, thus offering a unique opportunity for studying how T3 controls postembryonic development in vertebrates. Numerous early studies have demonstrated that T3 affects frog metamorphosis through T3 receptor (TR)-mediated regulation of T3 response genes, where TR forms a heterodimer with RXR (9-cis retinoic acid receptor) and binds to T3 response elements (TREs) in T3 response genes to regulate their expression. We have previously identified many candidate direct T3 response genes in Xenopus tropicalis tadpole intestine. Among them is the proto-oncogene Ski, which encodes a nuclear protein with complex function in regulating cell fate. We show here that Ski is upregulated in the intestine and tail of premetamorphic tadpoles upon T3 treatment and its expression peaks at stage 62, the climax of metamorphosis. We have further discovered a putative TRE in the first exon that can bind to TR/RXR in vitro and mediate T3 regulation of the promoter in vivo. These data demonstrate that Ski is activated by T3 through TR binding to a TRE in the first exon during Xenopus tropicalis metamorphosis, implicating a role of Ski in regulating cell fate during metamorphosis.

Whole-genome sequencing identifies I-SceI-mediated transgene integration sites in Xenopus tropicalis snai2:eGFP line

Transgenesis with the meganuclease I-SceI is a safe and efficient method, but the underlying mechanisms remain unclear due to the lack of information on transgene localization. Using I-SceI, we previously developed a transgenic Xenopus tropicalis line expressing enhanced green fluorescent protein driven by the neural crest-specific snai2 promoter/enhancer, which is a powerful tool for studying neural crest development and craniofacial morphogenesis. Here, we carried out whole-genome shotgun sequencing for the snai2:eGFP embryos to identify the transgene integration sites. With a 19x sequencing coverage, we estimated that 6 copies of the transgene were inserted into the Xenopus tropicalis genome in the hemizygous transgenic embryos. Two transgene integration loci adjacent to each other were identified in a noncoding region on chromosome 1, possibly as a result of duplication after a single transgene insertion. Interestingly, genomic DNA at the boundaries of the transgene integration loci contains short sequences homologous to the I-SceI recognition site, suggesting that the integration was not random but probably mediated by sequence homology. To our knowledge, our work represents the first genome-wide sequencing study on a transgenic organism generated with I-SceI, which is useful for evaluating the potential genetic effects of I-SceI-mediated transgenesis and further understanding the mechanisms underlying this transgenic method.

Differential maturation and chaperone dependence of the paralogous protein kinases DYRK1A and DYRK1B

The HSP90/CDC37 chaperone system not only assists the maturation of many protein kinases but also maintains their structural integrity after folding. The interaction of mature kinases with the HSP90/CDC37 complex is governed by the conformational stability of the catalytic domain, while the initial folding of the protein kinase domain is mechanistically less well characterized. DYRK1A (Dual-specificity tyrosine (Y)-phosphorylation Regulated protein Kinase 1A) and DYRK1B are closely related protein kinases with discordant HSP90 client status. DYRK kinases stoichiometrically autophosphorylate on a tyrosine residue immediately after folding, which served us as a traceable marker of successful maturation. In the present study, we used bacterial expression systems to compare the capacity of autonomous maturation of DYRK1A and DYRK1B in the absence of eukaryotic cofactors or chaperones. Under these conditions, autophosphorylation of human DYRK1B was severely compromised when compared with DYRK1A or DYRK1B orthologs from zebrafish and Xenopus. Maturation of human DYRK1B could be restored by bacterial expression at lower temperatures, suggesting that folding was not absolutely dependent on eukaryotic chaperones. The differential folding properties of DYRK1A and DYRK1B were largely due to divergent sequences of the C-terminal lobes of the catalytic domain. Furthermore, the mature kinase domain of DYRK1B featured lower thermal stability than that of DYRK1A when exposed to heat challenge in vitro or in living cells. In summary, our study enhances the mechanistic understanding of the differential thermodynamic properties of two closely related protein kinases during initial folding and as mature kinases.

Echinobase: leveraging an extant model organism database to build a knowledgebase supporting research on the genomics and biology of echinoderms

Echinobase (www.echinobase.org) is a third generation web resource supporting genomic research on echinoderms. The new version was built by cloning the mature Xenopus model organism knowledgebase, Xenbase, refactoring data ingestion pipelines and modifying the user interface to adapt to multispecies echinoderm content. This approach leveraged over 15 years of previous database and web application development to generate a new fully featured informatics resource in a single year. In addition to the software stack, Echinobase uses the private cloud and physical hosts that support Xenbase. Echinobase currently supports six echinoderm species, focused on those used for genomics, developmental biology and gene regulatory network analyses. Over 38 000 gene pages, 18 000 publications, new improved genome assemblies, JBrowse genome browser and BLAST + services are available and supported by the development of a new echinoderm anatomical ontology, uniformly applied formal gene nomenclature, and consistent orthology predictions. A novel feature of Echinobase is integrating support for multiple, disparate species. New genomes from the diverse echinoderm phylum will be added and supported as data becomes available. The common code development design of the integrated knowledgebases ensures parallel improvements as each resource evolves. This approach is widely applicable for developing new model organism informatics resources.

Sperm associated antigen 7 is activated by T3 during Xenopus tropicalis metamorphosis via a thyroid hormone response element within the first intron

Thyroid hormone (T3) affects many diverse physiological processes such as metabolism, organogenesis, and growth. The two highly related frog species, diploid Xenopus tropicalis and pseudo tetraploid Xenopus laevis, have been used as models for analyzing the effects of T3 during vertebrate development. T3 regulates T3-inducible gene transcription through T3 receptor (TR)-binding to T3-response elements (TREs). We have previously identified sperm associated antigen 7 (spag7) as a candidate T3 target gene that is potentially involved in adult stem cell development and/or proliferation during intestinal metamorphosis. To investigate whether T3 regulates spag7 directly at the transcriptional level via TR, we first conducted qRT-PCR to analyze its expression during natural and T3-induced metamorphosis and found that spag7 was up-regulated during natural metamorphosis in the intestine, tail, brain and hindlimb, peaking at the climax of metamorphosis in all those organs, and upon T3 treatment of premetamorphic tadpoles. Next, we demonstrated that an intronic TRE in spag7, first identified through bioinformatic analysis, could bind to TR in vitro and in vivo during metamorphosis. A dual luciferase assay utilizing a reconstituted frog oocyte transcription system showed that the TRE could mediate promoter activation by liganded TR. These results indicate that spag7 expression is directly regulated by T3 through the TRE in the first intron during metamorphosis, implicating a role for spag7 early during T3-regulated tissue remodeling and resorption.

Interaction between sex-determining genes from two species: clues from Xenopus hybrids

Hybrids provide an interesting model to study the evolution of sex-determining genes and sex chromosome systems as they offer the opportunity to see how independently evolving sex-determining pathways interact in vivo. In this context, the genus Xenopus represents a stimulating model, since species with non-homologous sex chromosomes and different sex-determining genes have been identified. In addition, the possibility of interspecies breeding is favoured in this group, which arose by alloploidization events, with species ploidy ranging from 2n = 2x = 20 in X. tropicalis (the only diploid representative of the genus) to 2n = 12x = 108 in X. ruwenzoriensis. To study how two sex-determining genes interact in vivo, X. laevis × X. tropicalis hybrids were produced. Gonadal differentiation in these hybrids revealed that the dm-w gene is dominant over X. tropicalis male-determining sex chromosomes (Y or Z), even though the Y chromosome is dominant in X. tropicalis (Y > W>Z). In the absence of the dm-w gene (the Z chromosome from X. laevis is present), the W chromosome from X. tropicalis is able to trigger ovarian development. Testicular differentiation will take place in the absence of W chromosomes from any of the parental species. The dominance/recessivity relationships between these sex-determining loci in the context of either parental genome remains unknown. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Spatio-temporal mRNA tracking in the early zebrafish embryo

Early stages of embryogenesis depend on subcellular localization and transport of maternal mRNA. However, systematic analysis of these processes is hindered by a lack of spatio-temporal information in single-cell RNA sequencing. Here, we combine spatially-resolved transcriptomics and single-cell RNA labeling to perform a spatio-temporal analysis of the transcriptome during early zebrafish development. We measure spatial localization of mRNA molecules within the one-cell stage embryo, which allows us to identify a class of mRNAs that are specifically localized at an extraembryonic position, the vegetal pole. Furthermore, we establish a method for high-throughput single-cell RNA labeling in early zebrafish embryos, which enables us to follow the fate of individual maternal transcripts until gastrulation. This approach reveals that many localized transcripts are specifically transported to the primordial germ cells. Finally, we acquire spatial transcriptomes of two xenopus species and compare evolutionary conservation of localized genes as well as enriched sequence motifs.

Mapping single-cell atlases throughout Metazoa unravels cell type evolution

Comparing single-cell transcriptomic atlases from diverse organisms can elucidate the origins of cellular diversity and assist the annotation of new cell atlases. Yet, comparison between distant relatives is hindered by complex gene histories and diversifications in expression programs. Previously, we introduced the self-assembling manifold (SAM) algorithm to robustly reconstruct manifolds from single-cell data (Tarashansky et al., 2019). Here, we build on SAM to map cell atlas manifolds across species. This new method, SAMap, identifies homologous cell types with shared expression programs across distant species within phyla, even in complex examples where homologous tissues emerge from distinct germ layers. SAMap also finds many genes with more similar expression to their paralogs than their orthologs, suggesting paralog substitution may be more common in evolution than previously appreciated. Lastly, comparing species across animal phyla, spanning sponge to mouse, reveals ancient contractile and stem cell families, which may have arisen early in animal evolution.

Comparative Distribution of Repetitive Sequences in the Karyotypes of Xenopus tropicalis and Xenopus laevis (Anura, Pipidae)

Xenopus laevis and its diploid relative, Xenopus tropicalis, are the most used amphibian models. Their genomes have been sequenced, and they are emerging as model organisms for research into disease mechanisms. Despite the growing knowledge on their genomes based on data obtained from massive genome sequencing, basic research on repetitive sequences in these species is lacking. This study conducted a comparative analysis of repetitive sequences in X. laevis and X. tropicalis. Genomic in situ hybridization (GISH) and fluorescence in situ hybridization (FISH) with Cot DNA of both species revealed a conserved enrichment of repetitive sequences at the ends of the chromosomes in these Xenopus species. The repeated sequences located on the short arm of chromosome 3 from X. tropicalis were not related to the sequences on the short arm of chromosomes 3L and 3S from X. laevis, although these chromosomes were homoeologous, indicating that these regions evolved independently in these species. Furthermore, all the other repetitive sequences in X. tropicalis and X. laevis may be species-specific, as they were not revealed in cross-species hybridizations. Painting experiments in X. laevis with chromosome 7 from X. tropicalis revealed shared sequences with the short arm of chromosome 3L. These regions could be related by the presence of the nucleolus organizer region (NOR) in both chromosomes, although the region revealed by chromosome painting in the short arm of chromosome 3L in X. laevis did not correspond to 18S + 28S rDNA sequences, as they did not colocalize. The identification of these repeated sequences is of interest as they provide an explanation to some problems already described in the genome assemblies of these species. Furthermore, the distribution of repetitive DNA in the genomes of X. laevis and X. tropicalis might be a valuable marker to assist us in understanding the genome evolution in a group characterized by numerous polyploidization events coupled with hybridizations.

A Role of Endogenous Histone Acetyltransferase Steroid Hormone Receptor Coactivator 3 in Thyroid Hormone Signaling During Xenopus Intestinal Metamorphosis

Background: Thyroid hormone (triiodothyronine [T3]) plays an important role in regulating vertebrate developmental, cellular, and metabolic processes via T3 receptor (TR). Liganded TR recruit coactivator complexes that include steroid receptor coactivators (SRC1, SRC2 or SRC3), which are histone acetyltransferases, to T3-responsive promoters. The functions of endogenous coactivators during T3-dependent mammalian adult organ development remain largely unclear, in part, due to the difficulty to access and manipulate late-stage embryos and neonates. We use Xenopus metamorphosis as a model for postembryonic development in vertebrates. This process is controlled by T3, involves drastic changes in every organ/tissue, and can be easily manipulated. We have previously found that SRC3 was upregulated in the intestine during amphibian metamorphosis. Methods: To determine the function of endogenous SRC3 during intestinal remodeling, we have generated Xenopus tropicalis animals lacking a functional SRC3 gene and analyzed the resulting phenotype. Results: Although removing SRC3 had no apparent effect on external development and animal gross morphology, the SRC3 (-/-) tadpoles displayed a reduction in the acetylation of histone H4 in the intestine compared with that in wild-type animals. Further, the expression of TR target genes was also reduced in SRC3 (-/-) tadpoles during intestinal remodeling. Importantly, SRC3 (-/-) tadpoles had inhibited/delayed intestinal remodeling during natural and T3-induced metamorphosis, including reduced adult intestinal stem cell proliferation and apoptosis of larval epithelial cells. Conclusion: Our results, thus, demonstrate that SRC3 is a critical component of the TR-signaling pathway in vivo during intestinal remodeling.

Thyroid Hormone Receptor Is Essential for Larval Epithelial Apoptosis and Adult Epithelial Stem Cell Development but Not Adult Intestinal Morphogenesis during Xenopus tropicalis Metamorphosis

Vertebrate postembryonic development is regulated by thyroid hormone (T3). Of particular interest is anuran metamorphosis, which offers several unique advantages for studying the role of T3 and its two nuclear receptor genes, TRα and TRβ, during postembryonic development. We have recently generated TR double knockout (TRDKO) Xenopus tropicalis animals and reported that TR is essential for the completion of metamorphosis. Furthermore, TRDKO tadpoles are stalled at the climax of metamorphosis before eventual death. Here we show that TRDKO intestine lacked larval epithelial cell death and adult stem cell formation/proliferation during natural metamorphosis. Interestingly, TRDKO tadpole intestine had premature formation of adult-like epithelial folds and muscle development. In addition, T3 treatment of premetamorphic TRDKO tadpoles failed to induce any metamorphic changes in the intestine. Furthermore, RNA-seq analysis revealed that TRDKO altered the expression of many genes in biological pathways such as Wnt signaling and the cell cycle that likely underlay the inhibition of larval epithelial cell death and adult stem cell development caused by removing both TR genes. Our data suggest that liganded TR is required for larval epithelial cell degeneration and adult stem cell formation, whereas unliganded TR prevents precocious adult tissue morphogenesis such as smooth-muscle development and epithelial folding.

ATP and large signaling metabolites flux through caspase-activated Pannexin 1 channels

Pannexin 1 (Panx1) is a membrane channel implicated in numerous physiological and pathophysiological processes via its ability to support release of ATP and other cellular metabolites for local intercellular signaling. However, to date, there has been no direct demonstration of large molecule permeation via the Panx1 channel itself, and thus the permselectivity of Panx1 for different molecules remains unknown. To address this, we expressed, purified, and reconstituted Panx1 into proteoliposomes and demonstrated that channel activation by caspase cleavage yields a dye-permeable pore that favors flux of anionic, large-molecule permeants (up to ~1 kDa). Large cationic molecules can also permeate the channel, albeit at a much lower rate. We further show that Panx1 channels provide a molecular pathway for flux of ATP and other anionic (glutamate) and cationic signaling metabolites (spermidine). These results verify large molecule permeation directly through caspase-activated Panx1 channels that can support their many physiological roles.

Crosstalk between codon optimality and cis-regulatory elements dictates mRNA stability

BACKGROUND

The regulation of messenger RNA (mRNA) stability has a profound impact on gene expression dynamics during embryogenesis. For example, in animals, maternally deposited mRNAs are degraded after fertilization to enable new developmental trajectories. Regulatory sequences in 3' untranslated regions (3'UTRs) have long been considered the central determinants of mRNA stability. However, recent work indicates that the coding sequence also possesses regulatory information. Specifically, translation in cis impacts mRNA stability in a codon-dependent manner. However, the strength of this mechanism during embryogenesis, as well as its relationship with other known regulatory elements, such as microRNA, remains unclear.

RESULTS

Here, we show that codon composition is a major predictor of mRNA stability in the early embryo. We show that this mechanism works in combination with other cis-regulatory elements to dictate mRNA stability in zebrafish and Xenopus embryos as well as in mouse and human cells. Furthermore, we show that microRNA targeting efficacy can be affected by substantial enrichment of optimal (stabilizing) or non-optimal (destabilizing) codons. Lastly, we find that one microRNA, miR-430, antagonizes the stabilizing effect of optimal codons during early embryogenesis in zebrafish.

CONCLUSIONS

By integrating the contributions of different regulatory mechanisms, our work provides a framework for understanding how combinatorial control of mRNA stability shapes the gene expression landscape.

Analysis of Thyroid Hormone Receptor α-Knockout Tadpoles Reveals That the Activation of Cell Cycle Program Is Involved in Thyroid Hormone-Induced Larval Epithelial Cell Death and Adult Intestinal Stem Cell Development During Xenopus tropicalis Metamorphosis

Background: There are two highly conserved thyroid hormone (triiodothyronine [T3]) receptor (TR) genes, TRα and TRβ, in all vertebrates, and the expression of TRα but not TRβ is activated earlier than T3 synthesis during development. In human, high levels of T3 are present during the several months around birth, and T3 deficiency during this period causes severe developmental abnormalities including skeletal and intestinal defects. It is, however, difficult to study this period in mammals as the embryos and neonates depend on maternal supply of nutrients for survival. However, Xenopus tropicalis undergoes a T3-dependent metamorphosis, which drastically changes essentially every organ in a tadpole. Of interest is intestinal remodeling, which involves near complete degeneration of the larval epithelium through apoptosis. Concurrently, adult intestinal stem cells are formed de novo and subsequently give rise to the self-renewing adult epithelial system, resembling intestinal maturation around birth in mammals. We have previously demonstrated that T3 signaling is essential for the formation of adult intestinal stem cells during metamorphosis. Methods: We studied the function of endogenous TRα in the tadpole intestine by using knockout animals and RNA-seq analysis. Results: We observed that removing endogenous TRα caused defects in intestinal remodeling, including drastically reduced larval epithelial cell death and adult intestinal stem cell proliferation. Using RNA-seq on intestinal RNA from premetamorphic wild-type and TRα-knockout tadpoles treated with or without T3 for one day, before any detectable T3-induced cell death and stem cell formation in the tadpole intestine, we identified more than 1500 genes, which were regulated by T3 treatment of the wild-type but not TRα-knockout tadpoles. Gene Ontology and biological pathway analyses revealed that surprisingly, these TRα-regulated genes were highly enriched with cell cycle-related genes, in addition to genes related to stem cells and apoptosis. Conclusions: Our findings suggest that TRα-mediated T3 activation of the cell cycle program is involved in larval epithelial cell death and adult epithelial stem cell development during intestinal remodeling.

Molecular and functional characterization of interferon regulatory factor 1 (IRF1) in amphibian Xenopus tropicalis

Interferon regulatory factor 1 (IRF1) is an important regulator in controlling the transcription of type I interferon genes, and its functions have been well-characterized in mammals, birds and fish. However, little information is available regarding the function of amphibian IRF1. In this study, an IRF1 gene homolog named as Xt-IRF1 was identified in the Western clawed frog (Xenopus tropicalis), an amphibian model specie widely used for comparative immunology research. Xt-IRF1 and IRF1 in other vertebrates possess similar genomic structure and flanking genes, and were grouped together to form a separate clade in phylogenetic tree. In addition, Xt-IRF1 gene was constitutively expressed in all tissues examined, with the highest expression level observed in spleen, and was inducible after poly(I:C) stimulation. Importantly, the expression of Xt-IRF1 was markedly induced by recombinant type I interferon, and Xt-IRF1 induced a strong activation of both IFNβ and ISRE promoters. The present study opens the door to investigate the roles of IRF1 in amphibians, and thus contributes to a better understanding of the functional evolution of IRFs in lower tetrapods.

Axial Skeletal Malformations in Genetically Modified Xenopus laevis and Xenopus tropicalis

Skeletal malformations in captive-bred, adult Xenopus spp., have not previously been reported. Here we describe 10 sexually mature, genetically modified laboratory frogs (6 Xenopus laevis and 4 Xenopus tropicalis) with axial skeletal abnormalities. The young adult frogs were described by veterinary staff as presenting with "hunchbacks," but were otherwise considered to be in good health. All affected frogs were genetically engineered using various techniques: transcription activator-like effector nucleases (TALEN) editing using thyroid hormone receptor α TALEN mRNA, restriction enzyme-mediated integration methods involving insertion of the inducible transgene pCAR/TRDN, or via I-SceI meganuclease transgenesis using either pDRTREdpTR-HS4 or pDPCrtTA-TREG-HS4 plasmid sequences. Radiographic findings (6 frogs) and gross necropsy (10 frogs) revealed vertebral column malformations and sacroiliac deformities that resulted in moderate to severe kyphosis and kyphoscoliosis. These findings were confirmed and additional skeletal abnormalities were identified using computed tomography to create a 3D reconstruction of 4 frogs. Additional findings visible on the 3D reconstructions included incomplete vertebral segmentation, malformed transverse processes, and a short and/or curved urostyle. Histopathologic findings included misshapen intervertebral joints with nonconforming articular surfaces, narrowed joint cavities, flattened or irregularly-formed articular cartilage, irregular maturation lines and nonpolarized chondrocytes, excess fibrocartilage, and evidence of irregular bone resorption and growth. While the specific etiology of the vertebral skeletal abnormalities remains unclear, possibilities include: 1) egg/oocyte physical manipulation (dejellying, microinjection, fertilization, etc.), 2) induction and expression of the transgenes, 3) inactivation (knockout) of existing genes by insertional mutagenesis, or 4) a combination of the above. Furthermore, the possibility of undetected changes in the macro or microenvironment, or a feature of the genetic background of the affected frogs cannot be ruled out.

Thyroid Hormone Induces DNA Demethylation in Xenopus Tadpole Brain

Thyroid hormone (T3) plays pivotal roles in vertebrate development, acting via nuclear T3 receptors (TRs) that regulate gene transcription by promoting post-translational modifications to histones. Methylation of cytosine residues in deoxyribonucleic acid (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 premetamorphic tadpoles with T3 for 24 or 48 hours 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 messenger ribonucleic acid 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.

A frog with three sex chromosomes that co-mingle together in nature: Xenopus tropicalis has a degenerate W and a Y that evolved from a Z chromosome

In many species, sexual differentiation is a vital prelude to reproduction, and disruption of this process can have severe fitness effects, including sterility. It is thus interesting that genetic systems governing sexual differentiation vary among-and even within-species. To understand these systems more, we investigated a rare example of a frog with three sex chromosomes: the Western clawed frog, Xenopus tropicalis. We demonstrate that natural populations from the western and eastern edges of Ghana have a young Y chromosome, and that a male-determining factor on this Y chromosome is in a very similar genomic location as a previously known female-determining factor on the W chromosome. Nucleotide polymorphism of expressed transcripts suggests genetic degeneration on the W chromosome, emergence of a new Y chromosome from an ancestral Z chromosome, and natural co-mingling of the W, Z, and Y chromosomes in the same population. Compared to the rest of the genome, a small sex-associated portion of the sex chromosomes has a 50-fold enrichment of transcripts with male-biased expression during early gonadal differentiation. Additionally, X. tropicalis has sex-differences in the rates and genomic locations of recombination events during gametogenesis that are similar to at least two other Xenopus species, which suggests that sex differences in recombination are genus-wide. These findings are consistent with theoretical expectations associated with recombination suppression on sex chromosomes, demonstrate that several characteristics of old and established sex chromosomes (e.g., nucleotide divergence, sex biased expression) can arise well before sex chromosomes become cytogenetically distinguished, and show how these characteristics can have lingering consequences that are carried forward through sex chromosome turnovers.

Structure of nucleosome-bound DNA methyltransferases DNMT3A and DNMT3B

CpG methylation by de novo DNA methyltransferases (DNMTs) 3A and 3B is essential for mammalian development and differentiation and is frequently dysregulated in cancer1. These two DNMTs preferentially bind to nucleosomes, yet cannot methylate the DNA wrapped around the nucleosome core2, and they favour the methylation of linker DNA at positioned nucleosomes3,4. Here we present the cryo-electron microscopy structure of a ternary complex of catalytically competent DNMT3A2, the catalytically inactive accessory subunit DNMT3B3 and a nucleosome core particle flanked by linker DNA. The catalytic-like domain of the accessory DNMT3B3 binds to the acidic patch of the nucleosome core, which orients the binding of DNMT3A2 to the linker DNA. The steric constraints of this arrangement suggest that nucleosomal DNA must be moved relative to the nucleosome core for de novo methylation to occur.

Paired Box 9 (PAX9), the RNA polymerase II transcription factor, regulates human ribosome biogenesis and craniofacial development

Dysregulation of ribosome production can lead to a number of developmental disorders called ribosomopathies. Despite the ubiquitous requirement for these cellular machines used in protein synthesis, ribosomopathies manifest in a tissue-specific manner, with many affecting the development of the face. Here we reveal yet another connection between craniofacial development and making ribosomes through the protein Paired Box 9 (PAX9). PAX9 functions as an RNA Polymerase II transcription factor to regulate the expression of proteins required for craniofacial and tooth development in humans. We now expand this function of PAX9 by demonstrating that PAX9 acts outside of the cell nucleolus to regulate the levels of proteins critical for building the small subunit of the ribosome. This function of PAX9 is conserved to the organism Xenopus tropicalis, an established model for human ribosomopathies. Depletion of pax9 leads to craniofacial defects due to abnormalities in neural crest development, a result consistent with that found for depletion of other ribosome biogenesis factors. This work highlights an unexpected layer of how the making of ribosomes is regulated in human cells and during embryonic development.

CRISPR/Cas9 mediated mutation of the mtnr1a melatonin receptor gene causes rod photoreceptor degeneration in developing Xenopus tropicalis

Nighttime surges in melatonin levels activate melatonin receptors, which synchronize cellular activities with the natural light/dark cycle. Melatonin receptors are expressed in several cell types in the retina, including the photon-sensitive rods and cones. Previous studies suggest that long-term photoreceptor survival and retinal health is in part reliant on melatonin orchestration of circadian homeostatic activities. This scenario would accordingly envisage that disruption of melatonin receptor signaling is detrimental to photoreceptor health. Using in vivo CRISPR/Cas9 genomic editing, we discovered that a small deletion mutation of the Mel1a melatonin receptor (mtnr1a) gene causes a loss of rod photoreceptors in retinas of developing Xenopus tropicalis heterozygous, but not homozygous mutant tadpoles. Cones were relatively spared from degeneration, and the rod loss phenotype was not obvious after metamorphosis. Localization of Mel1a receptor protein appeared to be about the same in wild type and mutant retinas, suggesting that the mutant protein is expressed at some level in mutant retinal cells. The severe impact on early rod photoreceptor viability may signify a previously underestimated critical role in circadian influences on long-term retinal health and preservation of sight. These data offer evidence that disturbance of homeostatic, circadian signaling, conveyed through a mutated melatonin receptor, is incompatible with rod photoreceptor survival.

The neurodevelopmental disorder risk gene DYRK1A is required for ciliogenesis and control of brain size in Xenopus embryos

DYRK1A [dual specificity tyrosine-(Y)-phosphorylation-regulated kinase 1 A] is a high-confidence autism risk gene that encodes a conserved kinase. In addition to autism, individuals with putative loss-of-function variants in DYRK1A exhibit microcephaly, intellectual disability, developmental delay and/or congenital anomalies of the kidney and urinary tract. DYRK1A is also located within the critical region for Down syndrome; therefore, understanding the role of DYRK1A in brain development is crucial for understanding the pathobiology of multiple developmental disorders. To characterize the function of this gene, we used the diploid frog Xenopus tropicalis We discover that Dyrk1a is expressed in ciliated tissues, localizes to ciliary axonemes and basal bodies, and is required for ciliogenesis. We also demonstrate that Dyrk1a localizes to mitotic spindles and that its inhibition leads to decreased forebrain size, abnormal cell cycle progression and cell death during brain development. These findings provide hypotheses about potential mechanisms of pathobiology and underscore the utility of X. tropicalis as a model system for understanding neurodevelopmental disorders.

Chromatin accessibility dynamics and single cell RNA-Seq reveal new regulators of regeneration in neural progenitors

Vertebrate appendage regeneration requires precisely coordinated remodeling of the transcriptional landscape to enable the growth and differentiation of new tissue, a process executed over multiple days and across dozens of cell types. The heterogeneity of tissues and temporally-sensitive fate decisions involved has made it difficult to articulate the gene regulatory programs enabling regeneration of individual cell types. To better understand how a regenerative program is fulfilled by neural progenitor cells (NPCs) of the spinal cord, we analyzed pax6-expressing NPCs isolated from regenerating Xenopus tropicalis tails. By intersecting chromatin accessibility data with single-cell transcriptomics, we find that NPCs place an early priority on neuronal differentiation. Late in regeneration, the priority returns to proliferation. Our analyses identify Pbx3 and Meis1 as critical regulators of tail regeneration and axon organization. Overall, we use transcriptional regulatory dynamics to present a new model for cell fate decisions and their regulators in NPCs during regeneration.

Cooperation of the NEIL3 and Fanconi anemia/BRCA pathways in interstrand crosslink repair

The NEIL3 DNA glycosylase is a base excision repair enzyme that excises bulky base lesions from DNA. Although NEIL3 has been shown to unhook interstrand crosslinks (ICL) in Xenopus extracts, how NEIL3 participants in ICL repair in human cells and its corporation with the canonical Fanconi anemia (FA)/BRCA pathway remain unclear. Here we show that the NEIL3 and the FA/BRCA pathways are non-epistatic in psoralen-ICL repair. The NEIL3 pathway is the major pathway for repairing psoralen-ICL, and the FA/BRCA pathway is only activated when NEIL3 is not present. Mechanistically, NEIL3 is recruited to psoralen-ICL in a rapid, PARP-dependent manner. Importantly, the NEIL3 pathway repairs psoralen-ICLs without generating double-strand breaks (DSBs), unlike the FA/BRCA pathway. In addition, we found that the RUVBL1/2 complex physically interact with NEIL3 and function within the NEIL3 pathway in psoralen-ICL repair. Moreover, TRAIP is important for the recruitment of NEIL3 but not FANCD2, and knockdown of TRAIP promotes FA/BRCA pathway activation. Interestingly, TRAIP is non-epistatic with both NEIL3 and FA pathways in psoralen-ICL repair, suggesting that TRAIP may function upstream of the two pathways. Taken together, the NEIL3 pathway is the major pathway to repair psoralen-ICL through a unique DSB-free mechanism in human cells.

Early Xenopus gene regulatory programs, chromatin states, and the role of maternal transcription factors

For decades, the early development of the Xenopus embryo has been an essential model system to study the gene regulatory mechanisms that govern cellular specification. At the top of the hierarchy of gene regulatory networks, maternally deposited transcription factors initiate this process and regulate the expression of zygotic genes that give rise to three distinctive germ layer cell types (ectoderm, mesoderm, and endoderm), and subsequent generation of organ precursors. The onset of germ layer specification is also closely coupled with changes associated with chromatin modifications. This review will examine the timing of maternal transcription factors initiating the zygotic genome activation, the epigenetic landscape of embryonic chromatin, and the network structure that governs the process.

Comprehensive RNA-Seq analysis of notochord-enriched genes induced during Xenopus tropicalis tail resorption

Anuran metamorphosis is perhaps the most dramatic developmental process regulated by thyroid hormone (TH). One of the unique processes that occur during metamorphosis is the complete resorption of the tail, including the notochord. Interestingly, recent gene knockout studies have shown that of the two known vertebrate TH receptors, TRα and TRβ, TRβ appears to be critical for notochord regression during tail resorption in Xenopus tropicalis. To determine the mechanisms underlying notochord regression, we carried out a comprehensive gene expression analysis in the notochord during metamorphosis by using RNA-Seq analyses of whole tail at stage 60 before any noticeable tail length reduction, whole tail at stage 63 when the tail length is reduced by about one half, and the rest of the tail at stage 63 after removing the notochord. This allowed us to identify many notochord-enriched, metamorphosis-induced genes at stage 63. Future studies on these genes should help to determine if they are regulated by TRβ and play any roles in notochord regression.

The histone methyltransferase KMT2D, mutated in Kabuki syndrome patients, is required for neural crest cell formation and migration

Kabuki syndrome is an autosomal dominant developmental disorder with high similarities to CHARGE syndrome. It is characterized by a typical facial gestalt in combination with short stature, intellectual disability, skeletal findings and additional features like cardiac and urogenital malformations, cleft palate, hearing loss and ophthalmological anomalies. The major cause of Kabuki syndrome are mutations in KMT2D, a gene encoding a histone H3 lysine 4 (H3K4) methyltransferase belonging to the group of chromatin modifiers. Here we provide evidence that Kabuki syndrome is a neurocrestopathy, by showing that Kmt2d loss-of-function inhibits specific steps of neural crest (NC) development. Using the Xenopus model system, we find that Kmt2d loss-of-function recapitulates major features of Kabuki syndrome including severe craniofacial malformations. A detailed marker analysis revealed defects in NC formation as well as migration. Transplantation experiments confirm that Kmt2d function is required in NC cells. Furthermore, analyzing in vivo and in vitro NC migration behavior demonstrates that Kmt2d is necessary for cell dispersion but not protrusion formation of migrating NC cells. Importantly, Kmt2d knockdown correlates with a decrease in H3K4 monomethylation and H3K27 acetylation supporting a role of Kmt2d in the transcriptional activation of target genes. Consistently, using a candidate approach, we find that Kmt2d loss-of-function inhibits Xenopus Sema3F expression, and overexpression of Sema3F can partially rescue Kmt2d loss-of-function defects. Taken together, our data reveal novel functions of Kmt2d in multiple steps of NC development and support the hypothesis that major features of Kabuki syndrome are caused by defects in NC development.

Morpholino Studies in Xenopus Brain Development

Antisense morpholino oligonucleotides (MOs) have become a valuable method to knockdown protein levels, to block with mRNA splicing and to interfere with miRNA function. MOs are widely used to alter gene expression in development of Xenopus and Zebrafish, where they are typically injected into the fertilized egg or blastomeres. Here we present methods to use electroporation to target delivery of MOs to the central nervous system of Xenopus laevis or Xenopus tropicalis tadpoles. Briefly, MO electroporation is accomplished by injecting MO solution into the brain ventricle and driving the MOs into cells of the brain with current passing between 2 platinum plate electrodes, positioned on either side of the target brain area. The method is relatively straightforward and uses standard equipment found in many neuroscience labs. A major advantage of electroporation is that it allows spatial and temporal control of MO delivery and therefore knockdown. Co-electroporation of MOs with cell type-specific fluorescent protein expression plasmids allows morphological analysis of cellular phenotypes. Furthermore, co-electroporation of MOs with rescuing plasmids allows assessment of specificity of the knockdown and phenotypic outcome. By combining MO-mediated manipulations with sophisticated assays of neuronal function, such as electrophysiological recording, behavioral assays, or in vivo time-lapse imaging of neuronal development, the functions of specific proteins and miRNAs within the developing nervous system can be elucidated. These methods can be adapted to apply antisense morpholinos to study protein and RNA function in a variety of complex tissues.

Maternal pluripotency factors initiate extensive chromatin remodelling to predefine first response to inductive signals

Embryonic development yields many different cell types in response to just a few families of inductive signals. The property of signal-receiving cells that determines how they respond to inductive signals is known as competence, and it differs in different cell types. Here, we explore the ways in which maternal factors modify chromatin to specify initial competence in the frog Xenopus tropicalis. We identify early-engaged regulatory DNA sequences, and infer from them critical activators of the zygotic genome. Of these, we show that the pioneering activity of the maternal pluripotency factors Pou5f3 and Sox3 determines competence for germ layer formation by extensively remodelling compacted chromatin before the onset of inductive signalling. This remodelling includes the opening and marking of thousands of regulatory elements, extensive chromatin looping, and the co-recruitment of signal-mediating transcription factors. Our work identifies significant developmental principles that inform our understanding of how pluripotent stem cells interpret inductive signals.

Identifying Common Features in the Activation of Melanocortin-2 Receptors: Studies on the Xenopus tropicalis Melanocortin-2 Receptor

The interaction between the pituitary hormone, adrenocorticotropin (ACTH), and melanocortin-2 receptor (MC2R) orthologs involves the H6 F7 R8 W9 and R/K15 K16 R17 R18 motifs in ACTH making contact with corresponding contact sites on MC2R. Earlier studies have localized the common HFRW binding site of all melanocortin receptors to residues in TM2, TM3, and TM6 that are located close to the extracellular space. The current study has identified residues in Xenopus tropicalis (xt) MC2R in TM4 (I158, F161), in EC2 (M166), and in TM5 (V172) that also are involved in activation of xtMC2R, and may be in the R/KKRR contact site of xtMC2R. These results are compared to earlier studies on the corresponding domains of human MC2R and rainbow trout MC2R in an effort to identify common features in the activation of teleost and tetrapod MC2R orthologs following stimulation with ACTH.

ML277 specifically enhances the fully activated open state of KCNQ1 by modulating VSD-pore coupling

Upon membrane depolarization, the KCNQ1 potassium channel opens at the intermediate (IO) and activated (AO) states of the stepwise voltage-sensing domain (VSD) activation. In the heart, KCNQ1 associates with KCNE1 subunits to form IKs channels that regulate heart rhythm. KCNE1 suppresses the IO state so that the IKs channel opens only to the AO state. Here, we tested modulations of human KCNQ1 channels by an activator ML277 in Xenopus oocytes. It exclusively changes the pore opening properties of the AO state without altering the IO state, but does not affect VSD activation. These observations support a distinctive mechanism responsible for the VSD-pore coupling at the AO state that is sensitive to ML277 modulation. ML277 provides insights and a tool to investigate the gating mechanism of KCNQ1 channels, and our study reveals a new strategy for treating long QT syndrome by specifically enhancing the AO state of native IKs currents.

Generation of Xenopus Haploid, Triploid, and Hybrid Embryos

Frog species of the genus Xenopus are widely used for studies of cell and developmental biology, and recent genome sequencing has revealed interesting phylogenetic relationships. Here we describe methods to generate haploid, triploid, and hybrid species starting from eggs and sperm of Xenopus laevis and Xenopus tropicalis that enable investigation of how genome size and content affect physiology at the organismal, cellular, and subcellular levels.

Functional integrity of honeybee (Apis mellifera L.) resistant to dieldrin γ-aminobutyric acid receptor channels conjugated with three fluorescent proteins

To generate an efficient tool used in Xenopus oocyte expression for in situ investigation of channel receptor expression, distribution and function, the C-terminus of the honeybee (Apis mellifera L.) resistant to dieldrin (RDL) subunit was fused with *FP, including monomeric red, enhanced yellow or enhanced green fluorescent protein (referred to as mRFP, EYFP and EGFP, respectively). In the present study, all fused *FP-AmRDLs could be visualized using fluorescence and laser confocal microscopy in cRNA-injected oocytes. Fluorescence was distributed isotropically in the cellular membrane. The potencies of the agonist γ-aminobutyric acid (GABA), but not β-alanine, and the test antagonists (fipronil, flufiprole, dieldrin, α-endosulfan, bifenazate and avermectin B1a) in the *FP-AmRDL receptor did not significantly differ from that of the untagged receptor with two-electrode voltage clamp detection. The half maximal effective concentrations (EC50 s) of GABA in AmRDL, EGFP-AmRDL, EYFP-AmRDL and mRFP-AmRDL receptors were 11.98, 12.61, 18.92 and 22.11 μM, respectively, and those of β-alanine were 651.6, 629.6, 1643.0 and 2146.0 μM, respectively. Inhibition percentages of test antagonists against *FP-AmRDL and AmRDL were not significantly different from each other. Overall, the consistency in functional properties between *FP-AmRDL and AmRDL receptors makes pGH19-*FP a promising tool for further in situ investigation of GABA receptors.

Regulation of a plant aquaporin by a Casparian strip membrane domain protein-like

The absorption of soil water by roots allows plants to maintain their water status. At the endodermis, water transport can be affected by initial formation of a Casparian strip and further deposition of suberin lamellas and regulated by the function of aquaporins. Four Casparian strip membrane domain protein-like (CASPL; CASPL1B1, CASPL1B2, CASPL1D1, and CASPL1D2) were previously shown to interact with PIP2;1. The present work shows that CASPL1B1, CASPL1B2, and CASPL1D2 are exclusively expressed in suberized endodermal cells, suggesting a cell-specific role in suberization and/or water transport regulation. When compared with wild-type plants, and by contrast to caspl1b1*caspl1b2 double loss of function, caspl1d1*caspl1d2 double mutants showed, in some control or NaCl stress experiments and not upon abscisic acid (ABA) treatment, a weak enlargement of the continuous suberization zone. None of the mutants showed root hydraulic conductivity (Lp_r ) phenotype, whether in control, NaCl, or ABA treatment conditions. The data suggest a slight negative role for CASPL1D1 and CASPL1D2 in suberization under control or salt stress conditions, with no major impact on whole root transport functions. At the molecular level, CASPL1B1 was able to physically interact with PIP2;1 and potentially could influence the regulation of aquaporins by acting on their phosphorylated form.

Cytochrome P450 1A transcript is a suitable biomarker of both exposure and response to diluted bitumen in developing frog embryos

In order for Alberta's thick bitumen to be transported through pipelines, condensates are added creating a diluted bitumen (dilbit) mixture. Recent pipeline expansion projects have generated concern about potential dilbit spills on aquatic wildlife health. Studies have suggested that polycyclic aromatic compounds (PACs) are toxic to aquatic vertebrates and could potentially also interfere with their endocrine system. The research objectives of this study were to investigate the toxicity of dilbit to developing frog embryos and to identify the molecular mechanisms of action involved. Fertilized embryos of Western clawed frog (Silurana tropicalis) were exposed for 72 h to water accommodated fractions (WAF; 0.7-8.9 μg/L TPACs) and chemically-enhanced WAFs (CEWAF; 0.09-56.7 μg/L TPACs) of Access Western Blend (AWB) and Cold Lake Blend (CLB) dilbits. Both dilbit's CEWAFs significantly increased embryonic mortality and malformation incidence in the highest treatments tested, while WAF treatments led to no visible toxic effects. Increases of the cytochrome P450 1A (cyp1a) mRNA levels were observed for all WAF and CEWAF dilbit treatments suggesting that phase I detoxification is activated in the dilbit-exposed larvae. When exposed to PAC concentrations ranging from 0.09 to 8.9 μg/L, the frogs displayed no observable malformations, but expressed significant increases of cyp1a mRNA levels (2- to 25-fold; indicating a suitable biomarker of exposure); however, when concentrations were of 46.6 μg/L or higher, both malformed frog phenotype and induction of cyp1a mRNA level (>250-fold) were measured (indicating a suitable biomarker of response). The expression of several genes related to cellular detoxification and endocrine disruption were also measured, but were not significantly altered by the treatments. In sum, cyp1a mRNA level is a highly sensitive endpoint to measure subtle molecular changes induced by PAC exposure in the frog embryos and larvae, and data suggest that PAC concentration higher than 46 μg/L would be toxic to the developing S. tropicalis.

CRISPR/Cas9-Mediated Knockout of Rb1 in Xenopus tropicalis

At this time, no molecular targeted therapies exist for treatment of retinoblastoma. This can be, in part, attributed to the lack of animal models that allow for both rapid identification of novel therapeutic targets and hypothesis driven drug testing. Within this scope, we have recently reported the first genuine genetic nonmammalian retinoblastoma cancer model within the aquatic model organism Xenopus tropicalis (Naert et al., Sci Rep 6: 35263, 2016). Here we describe the methods to generate rb1 mosaic mutant Xenopus tropicalis by employing the CRISPR/Cas9 technology. In depth, we discuss short guide RNA (sgRNA) design parameters, generation, quality control, quantification, and delivery followed by several methods for assessing genome editing efficiencies. As such the reader should be capable, by minor changes to the methods described here, to (co-) target rb1 or any one or multiple gene(s) within the Xenopus tropicalis genome by multiplex CRISPR/Cas9 methodology.