To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We demonstrate the allotetraploid origin of X. laevis by partitioning its genome into two homeologous subgenomes, marked by distinct families of “fossil” transposable elements. Based on the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged ~34 million years ago (Mya) and combined to form an allotetraploid ~17–18 Mya. 56% of all genes are retained in two homeologous copies. Protein function, gene expression, and the amount of flanking conserved sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

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.

Y-linked Dmy (also called dmrt1bY) in the teleost fish medaka, W-linked Dm-W in the African clawed frog (Xenopus laevis), and Z-linked Dmrt1 in the chicken are all sex chromosome-linked Dmrt1 homologues required for sex determination. Dmy and Dm-W both are Dmrt1 palalogues evolved through Dmrt1 duplication, while chicken Dmrt1 is a Z-linked orthologue. The eutherian sex-determining gene, Sry, evolved from an allelic gene, Sox3. Here we analyzed the exon–intron structures of the Dmrt1 homologues of several vertebrate species through information from databases and by determining the transcription initiation sites in medaka, chicken, Xenopus, and mouse. Interestingly, medaka Dmrt1 and Dmy and Xenopus Dm-W and Dmrt1 have a noncoding-type first exon, while mouse and chicken Dmrt1 do not. We next compared the 5′-flanking sequences of the Dmrt1 noncoding and coding exons 1 of several vertebrate species and found conservation of the presumptive binding sites for some transcription factors. Importantly, based on the phylogenetic trees for Dmrt1 and Sox3 homologues, it was implied that the sex-determining gene Dmy, Dm-W, and Sry have a higher substitution rate than thier prototype genes. Finally, we discuss the evolutionary relationships between vertebrate sex chromosomes and the sex-determining genes Dmy/Dm-W and Sry, which evolved by neofunctionalization of Dmrt1 and Sox3, respectively, for sex determining function. We propose a coevolution model of sex determining gene and sex chromosome, in which undifferentiated sex chromosomes easily allow replacement of a sex-determining gene with another new one, while specialized sex chromosomes are restricted a particular sex-determining gene.Electronic supplementary materialThe online version of this article (doi:10.1007/s10577-011-9265-9) contains supplementary material, which is available to authorized users.

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.

Some DMRT family genes including arthropod dsx, nematode mab-3, and vertebrate dmrt1 are involved in sex determination and/or differentiation in bilaterian animals. Although there have been some reports about evolutionary analyses of the family by using its phylogenetic trees, it is still undecided as to whether these three sex determination-related genes share orthologous relationships or not. To clarify this question, we analyzed evolutional relationships among the family members in various bilaterians by using not only phylogenetic tree analysis, but also synteny analysis. We found that only four genes, dmrt2a/2b, dmrt3, dmrt4/5 and dmrt93B were commonly present in invertebrate bilateria. The syntenies of dmrt2a/2b-dmrt3 and dmrt4/5-dmrt93B are conserved before and after two rounds of whole genome duplication in the ancestral vertebrate. Importantly, this indicates that dmrt1 must have appeared in the common vertebrate ancestor. In addition, dmrt1, dsx, or mab-3 formed each different cluster at a distance in our phylogenetic tree. From these findings, we concluded that the three sex determination-related genes, dmrt1, dsx, and mab-3 have no orthologous relationships, and suggested independent evolution for sex determination and differentiation in the DMRT gene family. Our results may supply clues about why sex-determining systems have diverged during animal evolution.

Amphibian metamorphosis induced by T(3) involves programmed cell death and the differentiation of various types of cells in degenerated and reconstructed tissues. However, the signaling pathway that directs the T(3)-dependent cell-fate determinations remains unclear. TNF-alpha is a pleiotropic cytokine that affects diverse cellular responses. Engagement of TNF-alpha with its receptor (TNFR1) causes intracellular apoptotic and/or survival signaling. To investigate TNF signaling functions during anuran metamorphosis, we first identified Xenopus laevis orthologs of TNF (xTNF)-alpha and its receptor. We found that xTNF-alpha activated nuclear factor-kappaB in X. laevis A6 cells through the Fas-associated death domain and receptor-interacting protein 1. Interestingly, xTNF-alpha mRNA in blood cells showed prominent expression at prometamorphosis during metamorphosis. Next, to elucidate the apoptotic and/or survival signaling induced by xTNF-alpha in an in vitro model of metamorphosis, we established a vascular endothelial cell line, XLgoo, from X. laevis tadpole tail. XLgoo cells formed actin stress fibers and elongated in response to xTNF-alpha. T(3) induced apoptosis in these cells, but the addition of xTNF-alpha blocked the T(3)-induced apoptosis. In addition, treatment of the cells with T(3) for 2 d induced the expression of thyroid hormone receptor-beta and caspase-3, and this thyroid hormone receptor-beta induction was drastically repressed by xTNF-alpha. Furthermore, in organ culture of the tail, xTNF-alpha significantly attenuated the tail degeneration induced by T(3). These findings suggested that xTNF-alpha could protect vascular endothelial cells from apoptotic cell death induced by T(3) during metamorphosis and thereby participate in the regulation of cell fate.

The African clawed frog, Xenopus laevis, has a ZZ/ZW-type sex-determination system. We previously reported that a W-linked gene, Dm-W, can determine development as a female. However, the mechanisms of early sex differentiation remain unclear. We used microarrays to screen for genes with sexually dimorphic expression in ZZ and ZW gonads during early sex differentiation in X laevis and found several steroidogenic genes. Importantly, the steroid 17α-hydroxylase gene Cyp17a1 and the aromatase gene Cyp19a1 were highly expressed in ZZ and ZW gonads, respectively, just after sex determination. At this stage, we found that Cyp17a1, Cyp19a1, or both were expressed in the ZZ and ZW gonads in a unique mass-in-line structure, in which several masses of cells, each surrounded by a basement membrane, were aligned along the anteroposterior axis. In fact, during sex differentiation, ovarian cavities formed inside each mass of Cyp17a1- and Cyp19a1-positive cells in the ZW gonads. However, the mass-in-line structure disappeared during testicular development in the ZZ testes. These results suggested that the mass-in-line structure found in both ZZ and ZW gonads just after sex determination might be formed in advance to produce ovarian cavities and then oocytes. Consequently, we propose a view that the default sex may be female in the morphological aspect of gonads in X laevis.

No abstract