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.