High variability of expression profiles of homeologous genes for Wnt, Hh, Notch, and Hippo signaling pathways in Xenopus laevis

Temporal Transcriptomic Profiling of the Developing Xenopus laevis Eye

Retinal progenitor cells (RPCs) are a multipotent and highly proliferative population that give rise to all retinal cell types during organogenesis. Defining their molecular signature is a key step towards identifying suitable approaches to treat visual impairments. Here, we performed RNA-sequencing of whole eyes from Xenopus at three embryonic stages and used differential expression analysis to define the transcriptomic profiles of optic tissues containing proliferating and differentiating RPCs during retinogenesis. Gene Ontology and KEGG pathway analyses showed that genes associated with developmental pathways (including Wnt and Hedgehog signaling) were upregulated during the period of active RPC proliferation in early retinal development (Nieuwkoop Faber st. 24 and 27). Developing eyes had dynamic expression profiles and shifted to enrichment for metabolic processes and phototransduction during RPC progeny specification and differentiation (st. 35). Furthermore, conserved adult eye regeneration genes were also expressed during early retinal development including sox2, pax6, nrl, and Notch signaling components. The eye transcriptomic profiles presented here span RPC proliferation to retinogenesis and included regrowth-competent stages. Thus, our dataset provides a rich resource to uncover molecular regulators of RPC activity and will allow future studies to address regulators of RPC proliferation during eye repair and regrowth.

The heparan sulfate modification enzyme, Hs6st1, governs Xenopus neuroectodermal patterning by regulating distributions of Fgf and Noggin

The nervous system has various types of cells derived from three neuroectodermal regions: neural plate (NP), neural crest (NC), and preplacodal ectoderm (PPE). Differentiation of these regions is regulated by various morphogens. However, regulatory mechanisms of morphogen distribution in neural patterning are still debated. In general, an extracellular component, heparan sulfate (HS), is essential to regulate morphogen gradients by modulating morphogen binding. The present study focused on an HS modification enzyme, heparan sulfate 6-O-sulfotransferase 1 (Hs6st1), which is highly expressed during the neurula stage in Xenopus. Our present in situ hybridization analysis revealed that Hs6st1 is expressed in the lateral sensorial layer of neuroectoderm. Overexpression of Hs6st1 expands Sox3 (NP marker gene) expression, and slightly dampens FoxD3 (NC marker) expression. Hs6st1 knockout using the CRISPR/Cas9 system also expands the neural plate region, followed by retinal malformation. These results imply that 6-O sulfation, mediated by Hs6st1, selectively regulates morphogen distribution required for neuroectodermal patterning. Among morphogens required for patterning, Fgf8a accumulates on Hs6st1-expressing cells, whereas a secreted BMP antagonist, Noggin, diffuses away from those cells. Thus, cell-autonomous 6-O sulfation of HS at the sensorial layer of neuroectoderm also affects neuroectodermal patterning in neighboring regions, including neural plate and neural crest, not only through accumulation, but also through dispersal of specific morphogens.

Ndst1, a heparan sulfate modification enzyme, regulates neuroectodermal patterning by enhancing Wnt signaling in Xenopus

Neural tissue is derived from three precursor regions: neural plate, neural crest, and preplacodal ectoderm. These regions are determined by morphogen-mediated signaling. Morphogen distribution is generally regulated by binding to an extracellular matrix component, heparan sulfate (HS) proteoglycan. HS is modified by many enzymes, such as N-deacetyl sulfotransferase 1 (Ndst1), which is highly expressed in early development. However, functions of HS modifications in ectodermal patterning are largely unknown. In this study, we analyzed the role of Ndst1 using Xenopus embryos. We found that ndst1 was expressed in anterior neural plate and the trigeminal region at the neurula stage. ndst1 overexpression expanded the neural crest (NC) region, whereas translational inhibition reduced not only the trigeminal region, but also the adjacent NC region, especially the anterior part. At a later stage, ndst1 knocked-down embryos showed defects in cranial ganglion formation. We also found that Ndst1 activates Wnt signaling pathway at the neurula stage. Taken together, our results suggest that N-sulfonated HS accumulates Wnt ligand and activates Wnt signaling in ndst1-expressing cells, but that it inhibits signaling in non-ndst1-expressing cells, leading to proper neuroectodermal patterning.

Polyploidy as a Fundamental Phenomenon in Evolution, Development, Adaptation and Diseases

DNA replication during cell proliferation is 'vertical' copying, which reproduces an initial amount of genetic information. Polyploidy, which results from whole-genome duplication, is a fundamental complement to vertical copying. Both organismal and cell polyploidy can emerge via premature cell cycle exit or via cell-cell fusion, the latter giving rise to polyploid hybrid organisms and epigenetic hybrids of somatic cells. Polyploidy-related increase in biological plasticity, adaptation, and stress resistance manifests in evolution, development, regeneration, aging, oncogenesis, and cardiovascular diseases. Despite the prevalence in nature and importance for medicine, agri- and aquaculture, biological processes and epigenetic mechanisms underlying these fundamental features largely remain unknown. The evolutionarily conserved features of polyploidy include activation of transcription, response to stress, DNA damage and hypoxia, and induction of programs of morphogenesis, unicellularity, and longevity, suggesting that these common features confer adaptive plasticity, viability, and stress resistance to polyploid cells and organisms. By increasing cell viability, polyploidization can provide survival under stressful conditions where diploid cells cannot survive. However, in somatic cells it occurs at the expense of specific function, thus promoting developmental programming of adult cardiovascular diseases and increasing the risk of cancer. Notably, genes arising via evolutionary polyploidization are heavily involved in cancer and other diseases. Ploidy-related changes of gene expression presumably originate from chromatin modifications and the derepression of bivalent genes. The provided evidence elucidates the role of polyploidy in evolution, development, aging, and carcinogenesis, and may contribute to the development of new strategies for promoting regeneration and preventing cardiovascular diseases and cancer.

Gene Structure Analysis of Chemokines and Their Receptors in Allotetraploid Frog, Xenopus laevis

Chemokines, relatively small secreted proteins, are involved in cell migration and function in various biological events, including immunity, morphogenesis, and disease. Due to their nature, chemokines tend to be a target of hijacking of immunity by virus and therefore show an exceptionally high mutation rate. Xenopus laevis is considered an excellent model to investigate the effect of whole-genome duplication for gene family evolution. Because its allotetraploidization occurred around 17–18 million years ago, ancestral subgenomes L and S were well conserved. Based on the gene model of human and diploid frog Xenopus tropicalis, we identified 52 chemokine genes and 26 chemokine receptors in X. laevis. The retention rate of the gene in the X. laevis L and S subgenomes was 96% (45/47) and 68% (32/47), respectively. We conducted molecular phylogenetic analysis and found clear orthologies in all receptor genes but not in the ligand genes, suggesting rapid divergences of the ligand. dN/dS calculation demonstrated that dN/dS ratio greater than one was observed in the four ligand genes, cxcl8b.1.S, cxcl18.S, ccl21.S, and xcl1.L, but nothing in receptor genes. These results revealed that the whole-genome duplication promotes diversification of chemokine ligands in X. laevis while conserving the genes necessary for homeostasis, suggesting that selective pressure also supports a rapid divergence of the chemokines in amphibians.

Enhanced Loss of Retinoic Acid Network Genes in Xenopus laevis Achieves a Tighter Signal Regulation

Retinoic acid (RA) is a major regulatory signal during embryogenesis produced from vitamin A (retinol) by an extensive, autoregulating metabolic and signaling network to prevent fluctuations that result in developmental malformations. Xenopus laevis is an allotetraploid hybrid frog species whose genome includes L (long) and S (short) chromosomes from the originating species. Evolutionarily, the X. laevis subgenomes have been losing either L or S homoeologs in about 43% of genes to generate singletons. In the RA network, out of the 47 genes, about 47% have lost one of the homoeologs, like the genome average. Interestingly, RA metabolism genes from storage (retinyl esters) to retinaldehyde production exhibit enhanced gene loss with 75% singletons out of 28 genes. The effect of this gene loss on RA signaling autoregulation was studied. Employing transient RA manipulations, homoeolog gene pairs were identified in which one homoeolog exhibits enhanced responses or looser regulation than the other, while in other pairs both homoeologs exhibit similar RA responses. CRISPR/Cas9 targeting of individual homoeologs to reduce their activity supports the hypothesis where the RA metabolic network gene loss results in tighter network regulation and more efficient RA robustness responses to overcome complex regulation conditions.

Activation of the Hippo Pathway in Rana sylvatica: Yapping Stops in Response to Anoxia

Wood frogs (Rana sylvatica) display well-developed anoxia tolerance as one component of their capacity to endure prolonged whole-body freezing during the winter months. Under anoxic conditions, multiple cellular responses are triggered to efficiently cope with stress by suppressing gene transcription and promoting activation of mechanisms that support cell survival. Activation of the Hippo signaling pathway initiates a cascade of protein kinase reactions that end with phosphorylation of YAP protein. Multiple pathway components of the Hippo pathway were analyzed via immunoblotting, qPCR or DNA-binding ELISAs to assess the effects of 24 h anoxia and 4 h aerobic recovery, compared with controls, on liver and heart metabolism of wood frogs. Immunoblot results showed significant increases in the relative levels of multiple proteins of the Hippo pathway representing an overall activation of the pathway in both organs under anoxia stress. Upregulation of transcript levels further confirmed this. A decrease in YAP and TEAD protein levels in the nuclear fraction also indicated reduced translocation of these proteins. Decreased DNA-binding activity of TEAD at the promoter region also suggested repression of gene transcription of its downstream targets such as SOX2 and OCT4. Furthermore, changes in the protein levels of two downstream targets of TEAD, OCT4 and SOX2, established regulated transcriptional activity and could possibly be associated with the activation of the Hippo pathway. Increased levels of TAZ in anoxic hearts also suggested its involvement in the repair mechanism for damage caused to cardiac muscles during anoxia. In summary, this study provides the first insights into the role of the Hippo pathway in maintaining cellular homeostasis in response to anoxia in amphibians.

A temporally resolved transcriptome for developing "Keller" explants of the Xenopus laevis dorsal marginal zone

BACKGROUND

Explanted tissues from vertebrate embryos reliably develop in culture and have provided essential paradigms for understanding embryogenesis, from early embryological investigations of induction, to the extensive study of Xenopus animal caps, to the current studies of mammalian gastruloids. Cultured explants of the Xenopus dorsal marginal zone ("Keller" explants) serve as a central paradigm for studies of convergent extension cell movements, yet we know little about the global patterns of gene expression in these explants.

RESULTS

In an effort to more thoroughly develop this important model system, we provide here a time-resolved bulk transcriptome for developing Keller explants.

CONCLUSIONS

The dataset reported here provides a useful resource for those using Keller explants for studies of morphogenesis and provide genome-scale insights into the temporal patterns of gene expression in an important tissue when explanted and grown in culture.

Molecular markers for corneal epithelial cells in larval vs. adult Xenopus frogs

Corneal Epithelial Stem Cells (CESCs) and their proliferative progeny, the Transit Amplifying Cells (TACs), are responsible for maintaining the integrity and transparency of the cornea. These stem cells (SCs) are widely used in corneal transplants and ocular surface reconstruction. Molecular markers are essential to identify, isolate and enrich for these cells, yet no definitive CESC marker has been established. An extensive literature survey shows variability in the expression of putative CESC markers among vertebrates; being attributed to species-specific variations, or other differences in developmental stages of these animals, approaches used in these studies and marker specificity. Here, we expanded the search for CESC markers using the amphibian model Xenopus laevis. In previous studies we found that long-term label retaining cells (suggestive of CESCs and TACs) are present throughout the larval basal corneal epithelium. In adult frogs, these cells become concentrated in the peripheral cornea (limbal region). Here, we used immunofluorescence to characterize the expression of nine proteins in the corneas of both Xenopus larvae and adults (post-metamorphic). We found that localization of some markers change between larval and adult stages. Markers such as p63, Keratin 19, and β1-integrin are restricted to basal corneal epithelial cells of the larvae. After metamorphosis their expression is found in basal and intermediate layer cells of the adult frog corneal epithelium. Another protein, Pax6 was expressed in the larval corneas, but surprisingly it was not detected in the adult corneal epithelium. For the first time we report that Tcf7l2 can be used as a marker to differentiate cornea vs. skin in frogs. Tcf7l2 is present only in the frog skin, which differs from reports indicating that the protein is expressed in the human cornea. Furthermore, we identified the transition between the inner, and the outer surface of the adult frog eyelid as a key boundary in terms of marker expression. Although these markers are useful to identify different regions and cellular layers of the frog corneal epithelium, none is unique to CESCs or TACs. Our results confirm that there is no single conserved CESC marker in vertebrates. This molecular characterization of the Xenopus cornea facilitates its use as a vertebrate model to understand the functions of key proteins in corneal homeostasis and wound repair.

Developmental toxicity and apoptosis in zebrafish embryos induced by low-dose γ-ray irradiation

In this paper, the developmental toxicity and apoptosis in zebrafish (Danio rerio) embryos induced by 0.01, 0.05, and 0.10-Gy γ-ray irradiation were investigated and verified by single cell gel electrophoresis, acridine orange staining, flow cytometry, transmission electron microscopy, digital gene expression sequencing, and Western blot analysis. DNA damage, deformity rates, and apoptosis of zebrafish embryos were found to increase significantly with the increase of irradiation dose, and survival and hatching rates significantly decreased when the irradiation dose exceeds 0.10 and 0.05 Gy, respectively. Exposure to 0.10-Gy γ-ray irradiation resulted in the swelling of cell mitochondria of zebrafish embryos and changes in their intracellular vacuoles. mRNA and protein expression levels of Shh (sonic hedgehog 19 KDa) and Smo (smoothened 86 KDa) of Hh signaling pathway associated with the development of early embryos significantly increased with the increase of irradiation dose. Expression of the AKT (56 KDa) and PiK3r3 (55 KDa) genes, which are anti-apoptotic and involved with the PI3K/Akt signaling pathway, significantly decreased, while expression of the bada gene, which is pro-apoptotic, significantly increased. The results show that γ-ray irradiations of 0.01 and 0.05 Gy can induce developmental toxicity and apoptosis in zebrafish embryos via Hh and PI3K/Akt signaling pathways, respectively.

Gli2 is required for the induction and migration of Xenopus laevis neural crest

The neural crest (NC) is a multipotent migratory embryonic population that is formed during late gastrulation and gives rise to a wide array of derivatives, including cells from the peripheral nervous system (PNS), the craniofacial bones and cartilages, peripheral glial cells, and melanocyte cells, among others. In this work we analyzed the role of the Hedgehog signaling pathway effector gli2 in Xenopus NC. We provide evidence that the gli2 gene is expressed in the prospective, premigratory and migratory NC. The use of a specific morpholino against gli2 and the pharmacological specific inhibitor GANT61 in different experimental approaches allowed us to determine that gli2 is required for the induction and specification of NC cells as a transcriptional activator. Moreover, gli2 also acts by reducing apoptosis in the NC without affecting its cell proliferation status. We also demonstrated that gli2 is required cell-autonomously for NC migration, and for the formation of NC derivatives such as the craniofacial cartilages, melanocytes and the cranial ganglia. Altogether, our results showed that gli2 is a key transcriptional activator to accomplish the proper specification and development of Xenopus NC cells.

Roles of two types of heparan sulfate clusters in Wnt distribution and signaling in Xenopus

Wnt proteins direct embryonic patterning, but the regulatory basis of their distribution and signal reception remain unclear. Here, we show that endogenous Wnt8 protein is distributed in a graded manner in Xenopus embryo and accumulated on the cell surface in a punctate manner in association with “N-sulfo-rich heparan sulfate (HS),” not with “N-acetyl-rich HS”. These two types of HS are differentially clustered by attaching to different glypicans as core proteins. N-sulfo-rich HS is frequently internalized and associated with the signaling vesicle, known as the Frizzled/Wnt/LRP6 signalosome, in the presence of Wnt8. Conversely, N-acetyl-rich HS is rarely internalized and accumulates Frzb, a secreted Wnt antagonist. Upon interaction with Frzb, Wnt8 associates with N-acetyl-rich HS, suggesting that N-acetyl-rich HS supports Frzb-mediated antagonism by sequestering Wnt8 from N-sulfo-rich HS. Thus, these two types of HS clusters may constitute a cellular platform for the distribution and signaling of Wnt8.

Wnt proteins mediate embryonic development but how protein localization and patterning is regulated is unclear. Here, the authors show that distinct structures with different heparan sulfate modifications (‘N-sulfo-rich’ and ‘N-acetyl-rich’) regulate cellular localization and signal transduction of Wnt8 in Xenopus.

A molecular atlas of the developing ectoderm defines neural, neural crest, placode, and nonneural progenitor identity in vertebrates

During vertebrate neurulation, the embryonic ectoderm is patterned into lineage progenitors for neural plate, neural crest, placodes and epidermis. Here, we use Xenopus laevis embryos to analyze the spatial and temporal transcriptome of distinct ectodermal domains in the course of neurulation, during the establishment of cell lineages. In order to define the transcriptome of small groups of cells from a single germ layer and to retain spatial information, dorsal and ventral ectoderm was subdivided along the anterior-posterior and medial-lateral axes by microdissections. Principal component analysis on the transcriptomes of these ectoderm fragments primarily identifies embryonic axes and temporal dynamics. This provides a genetic code to define positional information of any ectoderm sample along the anterior-posterior and dorsal-ventral axes directly from its transcriptome. In parallel, we use nonnegative matrix factorization to predict enhanced gene expression maps onto early and mid-neurula embryos, and specific signatures for each ectoderm area. The clustering of spatial and temporal datasets allowed detection of multiple biologically relevant groups (e.g., Wnt signaling, neural crest development, sensory placode specification, ciliogenesis, germ layer specification). We provide an interactive network interface, EctoMap, for exploring synexpression relationships among genes expressed in the neurula, and suggest several strategies to use this comprehensive dataset to address questions in developmental biology as well as stem cell or cancer research.

Vertebrate embryo germ layers become progressively regionalized by evolutionarily conserved molecular processes. Catching the early steps of this dynamic spatial cell diversification at the scale of the transcriptome was challenging, even with the advent of efficient RNA sequencing. We have microdissected complementary and defined areas of a single germ layer, the developing ectoderm, and explored how the transcriptome changes over time and space in the ectoderm during the differentiation of frog epidermis, neural plate, and neural crest. We have created EctoMap, a searchable interface using these regional transcriptomes, to predict the expression of the 31 thousand genes expressed in neurulae and their networks of co-expression, predictive of functional relationships. Through several examples, we illustrate how these data provide insights in development, cancer, evolution and stem cell biology.