Phosphorylation states change Otx2 activity for cell proliferation and patterning in the Xenopus embryo.

The homeodomain transcription factor Otx2 has essential roles in head and eye formation via the negative and positive regulation of its target genes, but it remains elusive how this dual activity of Otx2 affects cellular functions. In the current study, we first demonstrated that both exogenous and endogenous Otx2 are phosphorylated at multiple sites. Using Xenopus embryos, we identified three possible cyclin-dependent kinase (Cdk) sites and one Akt site, and analyzed the biological activities of phosphomimetic (4E) and nonphosphorylatable (4A) mutants for those sites. In the neuroectoderm, the 4E but not the 4A mutant downregulated the Cdk inhibitor gene p27xic1 (cdknx) and posterior genes, and promoted cell proliferation, possibly forming a positive-feedback loop consisting of Cdk, Otx2 and p27xic1 for cell proliferation, together with anteriorization. Conversely, the 4A mutant functioned as an activator on its own and upregulated the expression of eye marker genes, resulting in enlarged eyes. Consistent with these results, the interaction of Otx2 with the corepressor Tle1 is suggested to be phosphorylation dependent. These data suggest that Otx2 orchestrates cell proliferation, anteroposterior patterning and eye formation via its phosphorylation state.

Sox13 is a novel early marker for hair follicle development.

Sox13, a group D member of the Sry-related high-mobility group box (Sox) transcription factor family, is expressed in various tissues including the hair follicle. However, its spatiotemporal expression patterns in the hair follicle and its role in hair development remain to be elucidated. To address these questions, we generated Sox13-LacZ-knock-in mice (Sox13LacZ/+), in which the LacZ reporter gene was inserted in-frame into exon 2, which contains the translation initiation codon. X-gal staining in Sox13LacZ/+ embryos revealed that Sox13 is initially expressed in the epithelial portion of the placode, and subsequently in the hair germ and the hair peg during early hair follicle development. In postnatal catagen and anagen, Sox13 was detected in the epithelial sheath, whereas in telogen, Sox13 was localized in the bulge region, where hair follicle stem cells reside. Immunohistochemistry with an anti-ß-galactosidase antibody and anti-hair keratin antibodies that specifically mark the different layers of the hair follicle revealed that Sox13 was predominantly expressed in the outer root sheath in anagen. However, the integumentary structures of Sox13LacZ/LacZ mice were grossly and histologically indistinguishable from those of wild type mice. These results suggest that although Sox13 is dispensable for epidermal and adnexal development, Sox13 is a useful marker for early hair follicle development.

Transcription factors Mix1 and VegT, relocalization of vegt mRNA, and conserved endoderm and dorsal specification in frogs.

Protein expression of the transcription factor genes mix1 and vegt characterized the presumptive endoderm in embryos of the frogs Engystomops randi, Epipedobates machalilla, Gastrotheca riobambae, and Eleutherodactylus coqui, as in Xenopus laevis embryos. Protein VegT was detected in the animal hemisphere of the early blastula in all frogs, and only the animal pole was VegT-negative. This finding stimulated a vegt mRNA analysis in X. laevis eggs and embryos. vegt mRNA was detected in the animal region of X. laevis eggs and early embryos, in agreement with the VegT localization observed in the analyzed frogs. Moreover, a dorso-animal relocalization of vegt mRNA occurred in the egg at fertilization. Thus, the comparative analysis indicated that vegt may participate in dorsal development besides its known roles in endoderm development, and germ-layer specification. Zygotic vegt (zvegt) mRNA was detected as a minor isoform besides the major maternal (mvegt) isoform of the X. laevis egg. In addition, α-amanitin-insensitive vegt transcripts were detected around vegetal nuclei of the blastula. Thus, accumulation of vegt mRNA around vegetal nuclei was caused by relocalization rather than new mRNA synthesis. The localization of vegt mRNA around vegetal nuclei may contribute to the identity of vegetal blastomeres. These and previously reportedly localization features of vegt mRNA and protein derive from the master role of vegt in the development of frogs. The comparative analysis indicated that the strategies for endoderm, and dorsal specification, involving vegt and mix1, have been evolutionary conserved in frogs.

Genes coding for cyclin-dependent kinase inhibitors are fragile in Xenopus.

Cell proliferation is strictly regulated by the dosage balance among cell-cycle regulators such as CDK/cyclin complexes and CDK-Inhibitors. Even in the allotetraploid genome of Xenopus laevis, the dosage balance must be maintained for animals to stay alive, and the duplicated homeologous genes seem to have gradually changed, through evolution, resulting in the best genes for them to thrive. In the Xenopus laevis genome, while homeologous gene pairs of CDKs are fundamentally maintained and a few cyclin genes are amplified, homeologous gene pairs of the important CDK-Inhibitors, CDKn1c and CDKn2a, are deleted from chromosomes L and S. Although losses of CDKn1c and CDKn2a can lead to diseases in humans, their loss in X. laevis does not affect the animals' health. Also, another gene coding CDKn1b is lost besides CDKn1c and CDKn2a in the genome of Xenopus tropicalis. These findings suggest a high resistance of Xenopus to diseases. We also found that CDKn2c.S expression is higher than that of CDKn2c.L, and a conserved noncoding sequence (CNS) of CDKn2c genomic loci on X. laevis chromosome S and X. tropicalis has an enhancement activity in regulating the different expression. These findings together indicate a surprising fragility of CDK inhibitor gene loci in the Xenopus genome in spite of their importance, and may suggest that factors other than CDK-inhibitors decelerate cell-cycling in Xenopus.

Sex chromosome differentiation and the W- and Z-specific loci in Xenopus laevis.

Genetic sex-determining systems in vertebrates include two basic types of heterogamety; XX (female)/XY (male) and ZZ (male)/ZW (female) types. The African clawed frog Xenopus laevis has a ZZ/ZW-type sex-determining system. In this species, we previously identified a W-specific sex (female)-determining gene dmw, and specified W and Z chromosomes, which could be morphologically indistinguishable (homomorphic). In addition to dmw, we most recently discovered two genes, named scanw and ccdc69w, and one gene, named capn5z in the W- and Z-specific regions, respectively. In this study, we revealed the detail structures of the W/Z-specific loci and genes. Sequence analysis indicated that there is almost no sequence similarity between 278kb W-specific and 83kb Z-specific sequences on chromosome 2Lq32-33, where both the transposable elements are abundant. Synteny and phylogenic analyses indicated that all the W/Z-specific genes might have emerged independently. Expression analysis demonstrated that scanw and ccdc69w or capn5z are expressed in early differentiating ZW gonads or testes, thereby suggesting possible roles in female or male development, respectively. Importantly, the sex-determining gene (SDG) dmw might have been generated after allotetraploidization, thereby indicating the construction of the new sex-determining system by dmw after species hybridization. Furthermore, by direct genotyping, we confirmed that diploid WW embryos developed into normal female frogs, which indicate that the Z-specific region is not essential for female development. Overall, these findings indicate that sex chromosome differentiation has started, although no heteromorphic sex chromosomes are evident yet, in X. laevis. Homologous recombination suppression might have promoted the accumulation of mutations and transposable elements, and enlarged the W/Z-specific regions, thereby resulting in differentiation of the W/Z chromosomes.

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

Cell signaling pathways, such as Wnt, Hedgehog (Hh), Notch, and Hippo, are essential for embryogenesis, organogenesis, and tissue homeostasis. In this study, we analyzed 415 genes involved in these pathways in the allotetraploid frog, Xenopus laevis. Most genes are retained in two subgenomes called L and S (193 homeologous gene pairs and 29 singletons). This conservation rate of homeologs is much higher than that of all genes in the X. laevis genome (86.9% vs 60.2%). Among singletons, 24 genes are retained in the L subgenome, a rate similar to the average for all genes (82.8% vs 74.6%). In addition, as general components of signal transduction, we also analyzed 32 heparan sulfate proteoglycan (HSPG)-related genes and eight TLE/Groucho transcriptional corepressors-related genes. In these gene sets, all homeologous pairs have been retained. Transcriptome analysis using RNA-seq data from developmental stages and adult tissues demonstrated that most homeologous pairs of signaling components have variable expression patterns, in contrast to the conservative expression profiles of homeologs for transcription factors. Our results indicate that homeologous gene pairs for cell signaling regulation have tended to become subfunctionalized after allotetraploidization. Diversification of signaling pathways by subfunctionalization of homeologs may enhance environmental adaptability. These results provide insights into the evolution of signaling pathways after polyploidization.

Genomic organization and modulation of gene expression of the TGF-ß and FGF pathways in the allotetraploid frog Xenopus laevis.

Inductive interactions mediated by the TGF-ß and FGF-MAPK pathways are essential for specification of the germ layers and embryonic body axes during early vertebrate embryogenesis. TGF-ß and FGF ligands signal through receptor Ser/Thr and Tyr kinases, respectively, and these signaling pathways cross-talk to regulate transcription and cell behavior. The allotetraploid Xenopus laevis and its ancestral diploid Xenopus tropicalis are versatile model organisms with which to study the inductive interactions and mechanisms of these signal transduction pathways. Here we have analyzed the draft genome of X. laevis with respect to the genomic organization and differential expression of genes in the TGF-ß and FGF pathways. Genomic structure and gene expression analyses of pathway components in X. laevis revealed that genetic modulations, including deletions resulting in singletons and differential expression of homeologs, have occurred frequently among extracellular regulatory factors of the TGF-ß pathway after allotetraploidization. Moreover, differential gene expression was found for factors regulating various cellular responses including co-receptors, decoy receptors, and intracellular negative regulators in both the TGF-ß and FGF-MAPK pathways. We summarize the patterns of genetic alterations in the allotetraploid frog X. laevis and discuss the importance of these changes with regard to developmental processes.

Identification and comparative analyses of Siamois cluster genes in Xenopus laevis and tropicalis.

Two siamois-related homeobox genes siamois (sia1) and twin (sia2), have been reported in Xenopus laevis. These genes are expressed in the blastula chordin- and noggin-expressing (BCNE) center and the Nieuwkoop center, and have complete secondary axis-inducing activity when over-expressed on the ventral side of the embryo. Using whole genome sequences of X. tropicalis and X. laevis, we identified two additional siamois-related genes, which are tandemly duplicated near sia1 and sia2 to form the siamois gene cluster. Four siamois genes in X. tropicalis are transcribed at blastula to gastrula stages. In X. laevis, the siamois gene cluster is present on both homeologous chromosomes, XLA3L and XLA3S. Transcripts from seven siamois genes (three on XLA3L and four on XLA3S) in X. laevis were detected at blastula to gastrula stages. A transcribed gene, sia1p. S, encodes an inactive protein without a homeodomain. When over-expressed ventrally, all siamois-related genes tested in this study except for sia1p. S induced a complete secondary axis, indicating that X. tropicalis and X. laevis have four and six active siamois-related genes, respectively. Of note, each gene required different amounts of mRNA for full activity. These results suggest the possibility that siamois cluster genes have functional redundancy to endow robustness and quickness to organizer formation in Xenopus species.

Conservatism and variability of gene expression profiles among homeologous transcription factors in Xenopus laevis.

Xenopus laevis has an allotetraploid genome of 3.1Gb, in contrast to the diploid genome of a closely related species, Xenopus tropicalis. Here, we identified 412 genes (189 homeolog pairs, one homeologous gene cluster pair, and 28 singletons) encoding transcription factors (TFs) in the X. laevis genome by comparing them with their orthologs from X. tropicalis. Those genes include the homeobox gene family (Mix/Bix, Lhx, Nkx, Paired, POU, and Vent), Sox, Fox, Pax, Dmrt, Hes, GATA, T-box, and some clock genes. Most homeolog pairs for TFs are retained in two X. laevis subgenomes, named L and S, at higher than average rates (87.1% vs 60.2%). Among the 28 singletons, 82.1% were deleted from chromosomes of the S subgenome, a rate similar to the genome-wide average (82.1% vs 74.6%). Interestingly, nkx2-1, nkx2-8, and pax9, which reside consecutively in a postulated functional gene cluster, were deleted from the S chromosome, suggesting cluster-level gene regulation. Transcriptome correlation analysis demonstrated that TF homeolog pairs tend to have more conservative developmental expression profiles than most other types of genes. In some cases, however, either of the homeologs may show strongly different spatio-temporal expression patterns, suggesting neofunctionalization, subfunctionalization, or nonfunctionalization after allotetraploidization. Analyses of otx1 suggests that homeologs with much lower expression levels have undergone greater amino acid sequence diversification. Our comprehensive study implies that TF homeologs are highly conservative after allotetraploidization, possibly because the DNA sequences that they bind were also duplicated, but in some cases, they differed in expression levels or became singletons due to dosage-sensitive regulation of their target genes.

Genome organization of the vg1 and nodal3 gene clusters in the allotetraploid frog Xenopus laevis.

Extracellular factors belonging to the TGF-ß family play pivotal roles in the formation and patterning of germ layers during early Xenopus embryogenesis. Here, we show that the vg1 and nodal3 genes of Xenopus laevis are present in gene clusters on chromosomes XLA1L and XLA3L, respectively, and that both gene clusters have been completely lost from the syntenic S chromosome regions. The presence of gene clusters and chromosome-specific gene loss were confirmed by cDNA FISH analyses. Sequence and expression analyses revealed that paralogous genes in the vg1 and nodal3 clusters on the L chromosomes were also altered compared to their Xenopus tropicalis orthologs. X. laevis vg1 and nodal3 paralogs have potentially become pseudogenes or sub-functionalized genes and are expressed at different levels. As X. tropicalis has a single vg1 gene on chromosome XTR1, the ancestral vg1 gene in X. laevis appears to have been expanded on XLA1L. Of note, two reported vg1 genes, vg1(S20) and vg1(P20), reside in the cluster on XLA1L. The nodal3 gene cluster is also present on X. tropicalis chromosome XTR3, but phylogenetic analysis indicates that nodal3 genes in X. laevis and X. tropicalis were independently expanded and/or evolved in concert within each cluster by gene conversion. These findings provide insights into the function and molecular evolution of TGF-ß family genes in response to allotetraploidization.

Clustered Xenopus keratin genes: A genomic, transcriptomic, and proteomic analysis.

Keratin genes belong to the intermediate filament superfamily and their expression is altered following morphological and physiological changes in vertebrate epithelial cells. Keratin genes are divided into two groups, type I and II, and are clustered on vertebrate genomes, including those of Xenopus species. Various keratin genes have been identified and characterized by their unique expression patterns throughout ontogeny in Xenopus laevis; however, compilation of previously reported and newly identified keratin genes in two Xenopus species is required for our further understanding of keratin gene evolution, not only in amphibians but also in all terrestrial vertebrates. In this study, 120 putative type I and II keratin genes in total were identified based on the genome data from two Xenopus species. We revealed that most of these genes are highly clustered on two homeologous chromosomes, XLA9_10 and XLA2 in X. laevis, and XTR10 and XTR2 in X. tropicalis, which are orthologous to those of human, showing conserved synteny among tetrapods. RNA-Seq data from various embryonic stages and adult tissues highlighted the unique expression profiles of orthologous and homeologous keratin genes in developmental stage- and tissue-specific manners. Moreover, we identified dozens of epidermal keratin proteins from the whole embryo, larval skin, tail, and adult skin using shotgun proteomics. In light of our results, we discuss the radiation, diversification, and unique expression of the clustered keratin genes, which are closely related to epidermal development and terrestrial adaptation during amphibian evolution, including Xenopus speciation.

Comprehensive analyses of hox gene expression in Xenopus laevis embryos and adult tissues.

From whole genome sequencing of an allotetraploid frog, Xenopus laevis, two homeologous sets (L and S) of four Hox clusters A through D (HoxA.L/S, HoxB.L/S, HoxC.L/S, and HoxD.L/S) and 13 paralogous groups (PGs) with 76 genes were identified, allowing us to carry out the first comprehensive analyses of hox gene expression in vertebrates. Expression of all hox genes during development and in adult tissues was analyzed by RNA-sequencing. The expression levels of most hox genes were similar between homeologs, but in some pairs, large differences were observed and several of these were confirmed by RT-PCR and whole mount in situ hybridization experiments. These results indicate that subfunctionalization of hox genes may have occurred since allotetraploidization. Furthermore, comprehensive analysis of hox gene expression during early development did not agree with the hypothesis of temporal collinearity especially in genes belonging to PG2 to PG10.

The medaka dhc2 mutant reveals conserved and distinct mechanisms of Hedgehog signaling in teleosts.

BACKGROUND: Primary cilia are essential for Hedgehog (Hh) signal transduction in vertebrates. Although the core components of the Hh pathway are highly conserved, the dependency on cilia in Hh signaling is considered to be lower in fish than in mice, suggesting the presence of species-specific mechanisms for Hh signal transduction. RESULTS: To precisely understand the role of cilia in Hh signaling in fish and explore the evolution of Hh signaling, we have generated a maternal-zygotic medaka (Oryzias latipes) mutant that lacks cytoplasmic dynein heavy chain 2 (dhc2; MZdhc2), a component required for retrograde intraflagellar transport. We found that MZdhc2 exhibited the shortened cilia and partial defects in Hh signaling, although the Hh defects were milder than zebrafish mutants which completely lack cilia. This result suggests that Hh activity in fish depends on the length of cilium. However, the activity of Hh signaling in MZdhc2 appeared to be higher than that in mouse Dnchc2 mutants, suggesting a lower requirement for cilia in Hh signaling in fish. We confirmed that Ptch1 receptor is exclusively localized on the cilium in fish as in mammals. Subsequent analyses revealed that Fused, an essential mediator for Hh signaling in Drosophila and fish but not in mammals, augments the activity of Hh signaling in fish as a transcriptional target of Hh signaling. CONCLUSIONS: Ciliary requirement for Hh signaling in fish is lower than that in mammals, possibly due to fused-mediated positive feedback in Hh signaling. The finding of this fish-specific augmentation provides a novel insight into the evolution of Hh signaling.

Dynamic in vivo binding of transcription factors to cis-regulatory modules of cer and gsc in the stepwise formation of the Spemann-Mangold organizer.

How multiple developmental cues are integrated on cis-regulatory modules (CRMs) for cell fate decisions remains uncertain. The Spemann-Mangold organizer in Xenopus embryos expresses the transcription factors Lim1/Lhx1, Otx2, Mix1, Siamois (Sia) and VegT. Reporter analyses using sperm nuclear transplantation and DNA injection showed that cerberus (cer) and goosecoid (gsc) are activated by the aforementioned transcription factors through CRMs conserved between X. laevis and X. tropicalis. ChIP-qPCR analysis for the five transcription factors revealed that cer and gsc CRMs are initially bound by both Sia and VegT at the late blastula stage, and subsequently bound by all five factors at the gastrula stage. At the neurula stage, only binding of Lim1 and Otx2 to the gsc CRM, among others, persists, which corresponds to their co-expression in the prechordal plate. Based on these data, together with detailed expression pattern analysis, we propose a new model of stepwise formation of the organizer, in which (1) maternal VegT and Wnt-induced Sia first bind to CRMs at the blastula stage; then (2) Nodal-inducible Lim1, Otx2, Mix1 and zygotic VegT are bound to CRMs in the dorsal endodermal and mesodermal regions where all these genes are co-expressed; and (3) these two regions are combined at the gastrula stage to form the organizer. Thus, the in vivo dynamics of multiple transcription factors highlight their roles in the initiation and maintenance of gene expression, and also reveal the stepwise integration of maternal, Nodal and Wnt signaling on CRMs of organizer genes to generate the organizer.

A New Nomenclature of Xenopus laevis Chromosomes Based on the Phylogenetic Relationship to Silurana/Xenopus tropicalis.

Xenopus laevis (XLA) is an allotetraploid species which appears to have undergone whole-genome duplication after the interspecific hybridization of 2 diploid species closely related to Silurana/Xenopus tropicalis (XTR). Previous cDNA fluorescence in situ hybridization (FISH) experiments have identified 9 sets of homoeologous chromosomes in X. laevis, in which 8 sets correspond to chromosomes 1-8 of X. tropicalis (XTR1-XTR8), and the last set corresponds to a fusion of XTR9 and XTR10. In addition, recent X. laevis genome sequencing and BAC-FISH experiments support this physiological relationship and show no gross chromosome translocation in the X. laevis karyotype. Therefore, for the benefit of both comparative cytogenetics and genome research, we here propose a new chromosome nomenclature for X. laevis based on the phylogenetic relationship and chromosome length, i.e. XLA1L, XLA1S, XLA2L, XLA2S, and so on, in which the numbering of XLA chromosomes corresponds to that in X. tropicalis and the postfixes 'L' and 'S' stand for 'long' and 'short' chromosomes in the homoeologous pairs, which can be distinguished cytologically by their relative size. The last chromosome set is named XLA9L and XLA9S, in which XLA9 corresponds to both XTR9 and XTR10, and hence, to emphasize the phylogenetic relationship to X. tropicalis, XLA9_10L and XLA9_10S are also used as synonyms.

The Helicobacter pylori CagA oncoprotein disrupts Wnt/PCP signaling and promotes hyperproliferation of pyloric gland base cells.

Helicobacter pylori strains that deliver the oncoprotein CagA into gastric epithelial cells are the major etiologic agents of upper gastric diseases including gastric cancer. CagA promotes gastric carcinogenesis through interactions with multiple host proteins. Here, we show that CagA also disrupts Wnt-dependent planar cell polarity (Wnt/PCP), which orients cells within the plane of an epithelium and coordinates collective cell behaviors such as convergent extension to enable epithelial elongation during development. Ectopic expression of CagA in Xenopus laevis embryos impaired gastrulation, neural tube formation, and axis elongation, processes driven by convergent extension movements that depend on the Wnt/PCP pathway. Mice specifically expressing CagA in the stomach epithelium had longer pyloric glands and mislocalization of the tetraspanin proteins VANGL1 and VANGL2 (VANGL1/2), which are critical components of Wnt/PCP signaling. The increased pyloric gland length was due to hyperproliferation of cells at the gland base, where Lgr5+ stem and progenitor cells reside, and was associated with fewer differentiated enteroendocrine cells. In cultured human gastric epithelial cells, the N terminus of CagA interacted with the C-terminal cytoplasmic tails of VANGL1/2, which impaired Wnt/PCP signaling by inducing the mislocalization of VANGL1/2 from the plasma membrane to the cytoplasm. Thus, CagA may contribute to the development of gastric cancer by subverting a Wnt/PCP-dependent mechanism that restrains pyloric gland stem cell proliferation and promotes enteroendocrine differentiation.