Multiple Transitions between Y Chromosome and Autosome in Tago's Brown Frog Species Complex.

Sex chromosome turnover is the transition between sex chromosomes and autosomes. Although many cases have been reported in poikilothermic vertebrates, their evolutionary causes and genetic mechanisms remain unclear. In this study, we report multiple transitions between the Y chromosome and autosome in the Japanese Tago's brown frog complex. Using chromosome banding and molecular analyses (sex-linked and autosomal single nucleotide polymorphisms, SNPs, from the nuclear genome), we investigated the frogs of geographic populations ranging from northern to southern Japan of two species, Rana tagoi and Rana sakuraii (2n = 26). Particularly, the Chiba populations of East Japan and Akita populations of North Japan in R. tagoi have been, for the first time, investigated here. As a result, we identified three different sex chromosomes, namely chromosomes 3, 7, and 13, in the populations of the two species. Furthermore, we found that the transition between the Y chromosome (chromosome 7) and autosome was repeated through hybridization between two or three different populations belonging to the two species, followed by restricted chromosome introgression. These dynamic sex chromosome turnovers represent the first such findings in vertebrates and imply that speciation associated with inter- or intraspecific hybridization plays an important role in sex chromosome turnover in frogs.

Parallel Evolution of Sex-Linked Genes across XX/XY and ZZ/ZW Sex Chromosome Systems in the Frog Glandirana rugosa.

Genetic sex-determination features male (XX/XY) or female heterogamety (ZZ/ZW). To identify similarities and differences in the molecular evolution of sex-linked genes between these systems, we directly compared the sex chromosome systems existing in the frog Glandirana rugosa. The heteromorphic X/Y and Z/W sex chromosomes were derived from chromosomes 7 (2n = 26). RNA-Seq, de novo assembly, and BLASTP analyses identified 766 sex-linked genes. These genes were classified into three different clusters (XW/YZ, XY/ZW, and XZ/YW) based on sequence identities between the chromosomes, probably reflecting each step of the sex chromosome evolutionary history. The nucleotide substitution per site was significantly higher in the Y- and Z-genes than in the X- and W- genes, indicating male-driven mutation. The ratio of nonsynonymous to synonymous nucleotide substitution rates was higher in the X- and W-genes than in the Y- and Z-genes, with a female bias. Allelic expression in gonad, brain, and muscle was significantly higher in the Y- and W-genes than in the X- and Z-genes, favoring heterogametic sex. The same set of sex-linked genes showed parallel evolution across the two distinct systems. In contrast, the unique genomic region of the sex chromosomes demonstrated a difference between the two systems, with even and extremely high expression ratios of W/Z and Y/X, respectively.

rBC2LCN-reactive SERPINA3 is a glycobiomarker candidate for pancreatic ductal adenocarcinoma.

Early detection is urgently needed to improve the patient's pancreatic ductal adenocarcinoma (PDAC) survival. Previously, we identified a novel tumor-associated glycan, H-type3, which is expressed on PDAC cells and is detected by rBC2LCN (recombinant N-terminal domain of BC2L-C identified from Burkholderia cenocepacia) lectin. Here, we identified that SERPINA3 is an rBC2LCN-reactive glycoprotein (BC2-S3) secreted from PDAC cells into the blood in patients with PDAC by liquid chromatography-tandem mass spectrometry analysis and lectin blotting. In immune staining, BC2-S3 was detected specifically in the tumor but not in normal tissues of PDAC. Lectin-ELISA was then developed to measure the serum level of BC2-S3 in healthy control (HC, n = 99) and patients with PDAC (n = 88). BC2-S3 exhibited higher in patients with PDAC than in those with HC. BC2-S3 showed similar diagnostic performance in all stages of PDAC (stages IA-IV, true positive rate = 76.1%, true negative rate = 81.8%) to CA19-9 (72.7%, 75.8%). Remarkably, BC2-S3 showed a significantly higher detection rate (89.7%) for early stage PDAC (IA-IIA) than CA19-9 (62.1%, P = 0.029). The combination of BC2-S3 and CA19-9 further improved the diagnostic ability for all stages of PDAC (81.8%, 87.9%). In conclusion, BC2-S3 is a glycobiomarker candidate for PDAC.

Evolutionary features of ligands and their receptors via protein-protein interactions and essentiality in primates.

Protein evolution rate is negatively correlated with several effectors, such as expression level, expression distribution, protein-protein interactions (PPIs), and essentiality for survival. These effectors can characterize the signaling pathways mediated by ligand-receptor binding. However, it is unclear whether these effectors are constraining factors on the pathway-specific evolution of ligands and receptors. To clarify the relation between the effectors and protein evolution (dN /dS ratio) in ligands and their receptors considering each signaling pathway, we investigated 377 proteins in 20 peptide/protein ligand groups and their receptor groups using 15 primate sequences. The dN /dS ratios between peptide/protein ligand groups and their receptor groups were positively correlated, suggesting the protein evolution under the influence of signaling pathway to which they belong. Comparing each signaling pathway, ligands and receptors mainly related to development and growth (FGF/Hedgehog/Notch/WNT groups) showed lower dN /dS ratios, higher PPI numbers, and higher essentiality, whereas those mainly related to immune process (CSF/IFN/IL/TNF groups) showed higher dN /dS ratios, lower PPI numbers, and lower essentiality. Most ligands and receptors were poorly expressed, and expression level was not a constraining factor on the protein evolution. These findings indicate that PPI and essentiality are constraining factors that characterize the pathway-specific evolution of ligands and receptors.

Highly rapid and diverse sex chromosome evolution in the Odorrana frog species complex.

Sex chromosomes in poikilothermal vertebrates are characterized by rapid and diverse evolution at the species or population level. Our previous study revealed that the Taiwanese frog Odorrana swinhoana (2n = 26) has a unique system of multiple sex chromosomes created by three sequential translocations among chromosomes 1, 3, and 7. To reveal the evolutionary history of sex chromosomes in the Odorrana species complex, we first identified the original, homomorphic sex chromosomes, prior to the occurrence of translocations, in the ancestral-type population of O. swinhoana. Then, we extended the investigation to a closely related Japanese species, Odorrana utsunomiyaorum, which is distributed on two small islands. We used a high-throughput nuclear genomic approach to analyze single-nucleotide polymorphisms and identify the sex-linked markers. Those isolated from the O. swinhoana ancestral-type population were found to be aligned to chromosome 1 and showed male heterogamety. In contrast, almost all the sex-linked markers isolated from O. utsunomiyaorum were heterozygous in females and homozygous in males and were aligned to chromosome 9. Morphologically, we confirmed chromosome 9 to be heteromorphic in females, showing a ZZ-ZW sex determination system, in which the W chromosomes were heterochromatinized in a stripe pattern along the chromosome axis. These results indicated that after divergence of the two species, the ancestral homomorphic sex chromosome 1 underwent highly rapid and diverse evolution, i.e., sequential translocations with two autosomes in O. swinhoana, and turnover to chromosome 9 in O. utsunomiyaorum, with a transition from XY to ZW heterogamety and change to heteromorphy.

Identification of ancestral sex chromosomes in the frog Glandirana rugosa bearing XX-XY and ZZ-ZW sex-determining systems.

Sex chromosomes constantly exist in a dynamic state of evolution: rapid turnover and change of heterogametic sex during homomorphic state, and often stepping out to a heteromorphic state followed by chromosomal decaying. However, the forces driving these different trajectories of sex chromosome evolution are still unclear. The Japanese frog Glandirana rugosa is one taxon well suited to the study on these driving forces. The species has two different heteromorphic sex chromosome systems, XX-XY and ZZ-ZW, which are separated in different geographic populations. Both XX-XY and ZZ-ZW sex chromosomes are represented by chromosome 7 (2n = 26). Phylogenetically, these two systems arose via hybridization between two ancestral lineages of West Japan and East Japan populations, of which sex chromosomes are homomorphic in both sexes and to date have not yet been identified. Identification of the sex chromosomes will give us important insight into the mechanisms of sex chromosome evolution in this species. Here, we used a high-throughput genomic approach to identify the homomorphic XX-XY sex chromosomes in both ancestral populations. Sex-linked DNA markers of West Japan were aligned to chromosome 1, whereas those of East Japan were aligned to chromosome 3. These results reveal that at least two turnovers across three different sex chromosomes 1, 3 and 7 occurred during evolution of this species. This finding raises the possibility that cohabitation of the two different sex chromosomes from ancestral lineages induced turnover to another new one in their hybrids, involving transition of heterogametic sex and evolution from homomorphy to heteromorphy.

Independent pseudogenizations and losses of sox15 during amniote diversification following asymmetric ohnolog evolution.

BACKGROUND: Four ohnologous genes (sox1, sox2, sox3, and sox15) were generated by two rounds of whole-genome duplication in a vertebrate ancestor. In eutherian mammals, Sox1, Sox2, and Sox3 participate in central nervous system (CNS) development. Sox15 has a function in skeletal muscle regeneration and has little functional overlap with the other three ohnologs. In contrast, the frog Xenopus laevis and zebrafish orthologs of sox15 as well as sox1-3 function in CNS development. We previously reported that Sox15 is involved in mouse placental development as neofunctionalization, but is pseudogenized in the marsupial opossum. These findings suggest that sox15 might have evolved with divergent gene fates during vertebrate evolution. However, knowledge concerning sox15 in other vertebrate lineages than therian mammals, anuran amphibians, and teleost fish is scarce. Our purpose in this study was to clarify the fate and molecular evolution of sox15 during vertebrate evolution. RESULTS: We searched for sox15 orthologs in all vertebrate classes from agnathans to mammals by significant sequence similarity and synteny analyses using vertebrate genome databases. Interestingly, sox15 was independently pseudogenized at least twice during diversification of the marsupial mammals. Moreover, we observed independent gene loss of sox15 at least twice during reptile evolution in squamates and crocodile-bird diversification. Codon-based phylogenetic tree and selective analyses revealed an increased dN/dS ratio for sox15 compared to the other three ohnologs during jawed vertebrate evolution. CONCLUSIONS: The findings revealed an asymmetric evolution of sox15 among the four ohnologs during vertebrate evolution, which was supported by the increased dN/dS values in cartilaginous fishes, anuran amphibians, and amniotes. The increased dN/dS value of sox15 may have been caused mainly by relaxed selection. Notably, independent pseudogenizations and losses of sox15 were observed during marsupial and reptile evolution, respectively. Both might have been caused by strong relaxed selection. The drastic gene fates of sox15, including neofunctionalization and pseudogenizations/losses during amniote diversification, might be caused by a release from evolutionary constraints.

SSEA-1-positive fibronectin is secreted by cells deviated from the undifferentiated state of human induced pluripotent stem cells.

Quality control for human induced pluripotent stem cells (hiPSCs) is important for efficient and stable production of hiPSC-derived cell therapy products to be used for transplantation. During cell culture, hiPSCs spontaneously undergo morphological changes and lose pluripotent properties. Such cells are termed deviated cells, which are altered from the undifferentiated state of hiPSCs, and express the early differentiation marker stage-specific embryonic antigen 1 (SSEA-1). In this study, we searched for soluble SSEA-1+ glycoproteins secreted from deviated cells generated by culturing hiPSCs in cell culture medium containing heat-inactivated supplements. Glycoproteins obtained from cell culture supernatants of SSEA-1+ deviated cells were enriched by an O-glycan binding lectin and blotted with anti-SSEA-1 antibody. A single protein band at >250 kDa specifically detected by anti-SSEA-1 antibody was identified as fibronectin (FN) by LC-MS/MS analysis and immunoprecipitation combined with western blotting, indicating that FN is a carrier protein of SSEA-1. We then constructed a sandwich enzyme-linked immunosorbent assay to detect SSEA-1+ FN secreted from deviated cells. This FN-SSEA-1 test proved to be both sensitive and specific, allowing for non-destructive detection of SSEA-1+ deviated cells within mixed cell population, with a lower limit of detection of 100 cells/mL. The developed assay may provide a standard technology for quality control of hiPSCs used for regenerative medicine.

A technique for removing tumourigenic pluripotent stem cells using rBC2LCN lectin.

INTRODUCTION: Tumourigenesis attributed to residual undifferentiated cells in a graft is considered to be a significant issue in cell therapy using human pluripotent stem cells. To ensure the safety of regenerative medicine derived from pluripotent stem cells, residual undifferentiated cells must be eliminated in the manufacturing process. We previously described the lectin probe rBC2LCN, which binds harmlessly and specifically to the cell surface of human pluripotent stem cells. We report here a technique using rBC2LCN to remove pluripotent cells from a heterogenous population to reduce the chance of teratoma formation. METHODS: We demonstrate a method for separating residual tumourigenic cells using rBC2LCN-bound magnetic beads. This technology is a novel use of their previous discovery that rBC2LCN is a lectin that selectively binds to pluripotent cells. We optimize and validate a method to remove hPSCs from a mixture with human fibroblasts using rBC2LCN-conjugated magnetic beads. RESULTS: Cells with the potential to form teratoma could be effectively eliminated from a heterogeneous cell population with biotin-labelled rBC2LCN and streptavidin-bound magnetic beads. The efficiency was measured by FACS, ddPCR, and animal transplantation, suggesting that magnetic cell separation using rBC2LCN is quite efficient for eliminating hPSCs from mixed cell populations. CONCLUSIONS: The removal of residual tumourigenic cells based on rBC2LCN could be a practical option for laboratory use and industrialisation of regenerative medicine using human pluripotent stem cells.

rBC2LCN lectin as a potential probe of early-stage HER2-positive breast carcinoma.

The recombinant N-terminal domain of BC2L-C lectin (rBC2LCN) is useful for detecting not only human pluripotent stem cells but also some cancers. However, the cancer types and stages that can be detected by rBC2LCN remain unclear. In this study, we identified the human breast carcinoma subtypes and stages that can be detected by rBC2LCN. Compared with rBC2LCN-negative breast carcinoma cell lines, the rBC2LCN-positive cells expressed higher levels of human epidermal growth factor receptor 2 (HER2) and epithelial marker genes. Importantly, rBC2LCN histochemical staining of human breast carcinoma tissues demonstrated the utility of rBC2LCN in detecting breast carcinoma types that express HER2 and have not spread much in the early phase of growth. We conclude that rBC2LCN may have potential as a detection probe and a drug delivery vehicle to identify and treat early-stage HER2-positive breast carcinoma.

Parallel Evolution of Two dmrt1-Derived Genes, dmy and dm-W, for Vertebrate Sex Determination.

Animal sex-determining genes, which bifurcate for female and male development, are diversified even among closely related species. Most of these genes emerged independently from various sex-related genes during species diversity as neofunctionalization-type genes. However, the common mechanisms of this divergent evolution remain poorly understood. Here, we compared the molecular evolution of two sex-determining genes, the medaka dmy and the clawed frog dm-W, which independently evolved from the duplication of the transcription factor-encoding masculinization gene dmrt1. Interestingly, we detected parallel amino acid substitutions, from serine (S) to threonine (T), on the DNA-binding domains of both ancestral DMY and DM-W, resulting from positive selection. Two types of DNA-protein binding experiments and a luciferase reporter assay demonstrated that these S-T substitutions could strengthen the DNA-binding abilities and enhance the transcriptional regulation function. These findings suggest that the parallel S-T substitutions may have contributed to the establishment of dmy and dm-W as sex-determining genes.

mRNA and miRNA expression profiles in an ectoderm-biased substate of human pluripotent stem cells.

The potential applications of human pluripotent stem cells, embryonic stem (ES) cells, and induced pluripotent stem (iPS) cells in cell therapy and regenerative medicine have been widely studied. The precise definition of pluripotent stem cell status during culture using biomarkers is essential for basic research and regenerative medicine. Culture conditions, including extracellular matrices, influence the balance between self-renewal and differentiation. Accordingly, to explore biomarkers for defining and monitoring the pluripotent substates during culture, we established different substates in H9 human ES cells by changing the extracellular matrix from vitronectin to Matrigel. The substate was characterised by low and high expression of the pluripotency marker R-10G epitope and the mesenchymal marker vimentin, respectively. Immunohistochemistry, induction of the three germ layers, and exhaustive expression analysis showed that the substate was ectoderm-biased, tended to differentiate into nerves, but retained the potential to differentiate into the three germ layers. Further integrated analyses of mRNA and miRNA microarrays and qPCR analysis showed that nine genes (COL9A2, DGKI, GBX2, KIF26B, MARCH1, PLXNA4, SLC24A4, TLR4, and ZHX3) were upregulated in the ectoderm-biased cells as ectoderm-biased biomarker candidates in pluripotent stem cells. Our findings provide important insights into ectoderm-biased substates of human pluripotent stem cells in the fields of basic research and regenerative medicine.

Independent evolution for sex determination and differentiation in the DMRT family in animals.

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.

The rBC2LCN-positive subpopulation of PC-3 cells exhibits cancer stem-like properties.

The recombinant lectin rBC2LCN is a useful marker for discriminating the undifferentiated status of human induced or embryonic stem cells. Recently, rBC2LCN has also been used for detecting some cancers and niche cells. However, the generality of which types of cells are detected by rBC2LCN is unclear. In this study, we demonstrated the potential of rBC2LCN as a probe for detecting and isolating cancer stem-like cells. Interestingly, flow cytometric analysis of various human cell lines indicated that the human prostate cancer cell line PC-3 consisted of rBC2LCN-positive and -negative subpopulations. Compared with the rBC2LCN-negative subpopulation, the rBC2LCN-positive subpopulation possessed representative features of cancer stem cells and malignancy, such as slow proliferation, increased cell motility, anchorage-independent growth, and drug resistance. The comprehensive expression profiles revealed that the rBC2LCN-positive subpopulation expressed higher levels of cancer stem cell markers. These findings indicate that rBC2LCN is useful for detecting not only pluripotent stem cells but also the cancer stem-like subpopulation of PC-3 cells. Pluripotent and cancer cells with rBC2LCN positivity would be important for future stem cell research.

Masculinization-Related Genes and Cell-Mass Structures During Early Gonadal Differentiation in the African Clawed Frog Xenopus laevis.

The African clawed frog Xenopus laevis has a female heterogametic ZZ/ZW-type sex-determining system. We previously discovered a W-linked female sex-determining gene dm-W that is involved in ovary formation, probably through the up-regulation of the estrogen synthesis genes cyp19a1 and foxl2. We also reported that a unique "mass-in-line structure", which disappears from ZZ gonads during early testicular development, might serve as the basis for ovary differentiation in ZW gonads. However, the molecular mechanisms underlying early masculinization are poorly understood. To elucidate the development of bipotential gonads into testes after sex determination in this species, we focused on the orthologs of five mammalian sex-related genes: three nuclear factor genes, dax1, sf1 (also known as ad4bp), and sox9, and two genes encoding members of the tumor growth factor-ß (TGF-ß) family, anti-Müllerian hormone (amh) and inhibin ßb (inhbb). Quantitative RT-PCR analysis revealed that the expression of dax1, sox9, amh, and inhbb or sf1 was greatly or slightly higher in ZZ than in ZW gonads during early sex development. In situ hybridization analysis revealed that amh and inhbb mRNAs were expressed in somatic cells on the inner and outer sides of cell masses in the mass-in-line structure, respectively, in the developing ZZ gonads. Interestingly, estrogen exposure prevented the disappearance of the mass-in-line structure in early developing ZZ tadpoles. These findings suggest that TGF-ß signaling is involved in the destruction of the mass-in-line structure, which may be maintained by estrogen.

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.

Molecular Evolution of Two Distinct dmrt1 Promoters for Germ and Somatic Cells in Vertebrate Gonads.

The transcription factor DMRT1 has important functions in two distinct processes, somatic-cell masculinization and germ-cell development in mammals. However, it is unknown whether the functions are conserved during evolution, and what mechanism underlies its expression in the two cell lineages. Our analysis of the Xenopus laevis and Silurana tropicalis dmrt1 genes indicated the presence of two distinct promoters: one upstream of the noncoding first exon (ncEx1), and one within the first intron. In contrast, only the ncEx1-upstream promoter was detected in the dmrt1 gene of the agnathan sand lamprey, which expressed dmrt1 exclusively in the germ cells. In X. laevis, the ncEx1- and exon 2-upstream promoters were predominantly used for germ-cell and somatic-cell transcription, respectively. Importantly, knockdown of the ncEx1-containing transcript led to reduced germ-cell numbers in X. laevis gonads. Intriguingly, two genetically female individuals carrying the knockdown construct developed testicles. Analysis of the reptilian leopard gecko dmrt1 revealed the absence of ncEx1. We propose that dmrt1 regulated germ-cell development in the vertebrate ancestor, then acquired another promoter in its first intron to regulate somatic-cell masculinization during gnathostome evolution. In the common ancestor of reptiles and mammals, only one promoter got function for both the two cell lineages, accompanied with the loss of ncEx1. In addition, we found a conserved noncoding sequence (CNS) in the dmrt1 5'-flanking regions only among amniote species, and two CNSs in the introns among most vertebrates except for agnathans. Finally, we discuss relationships between these CNSs and the promoters of dmrt1 during vertebrate evolution.

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.