Homoeologous chromosomes of Xenopus laevis are highly conserved after whole-genome duplication

Using Xenopus to discover new candidate genes involved in BOR and other congenital hearing loss syndromes

Hearing in infants is essential for brain development, acquisition of verbal language skills, and development of social interactions. Therefore, it is important to diagnose hearing loss soon after birth so that interventions can be provided as early as possible. Most newborns in the United States are screened for hearing deficits and commercially available next-generation sequencing hearing loss panels often can identify the causative gene, which may also identify congenital defects in other organs. One of the most prevalent autosomal dominant congenital hearing loss syndromes is branchio-oto-renal syndrome (BOR), which also presents with defects in craniofacial structures and the kidney. Currently, mutations in three genes, SIX1, SIX5, and EYA1, are known to be causative in about half of the BOR patients that have been tested. To uncover new candidate genes that could be added to congenital hearing loss genetic screens, we have combined the power of Drosophila mutants and protein biochemical assays with the embryological advantages of Xenopus, a key aquatic animal model with a high level of genomic similarity to human, to identify potential Six1 transcriptional targets and interacting proteins that play a role during otic development. We review our transcriptomic, yeast 2-hybrid, and proteomic approaches that have revealed a large number of new candidates. We also discuss how we have begun to identify how Six1 and co-factors interact to direct developmental events necessary for normal otic development.

Unleashing diversity through flexibility: The evolutionary journey of sex chromosomes in amphibians and reptiles

Sex determination systems have greatly diversified between amphibians and reptiles, with such as the different sex chromosome compositions within a single species and transition between temperature-dependent sex determination (TSD) and genetic sex determination (GSD). In most sex chromosome studies on amphibians and reptiles, the whole-genome sequence of Xenopous tropicalis and chicken have been used as references to compare the chromosome homology of sex chromosomes among each of these taxonomic groups, respectively. In the present study, we reviewed existing reports on sex chromosomes, including karyotypes, in amphibians and reptiles. Furthermore, we compared the identified genetic linkages of sex chromosomes in amphibians and reptiles with the chicken genome as a reference, which is believed to resemble the ancestral tetrapod karyotype. Our findings revealed that sex chromosomes in amphibians are derived from genetic linkages homologous to various chicken chromosomes, even among several frogs within single families, such as Ranidae and Pipidae. In contrast, sex chromosomes in reptiles exhibit conserved genetic linkages with chicken chromosomes, not only across most species within a single family, but also within closely related families. The diversity of sex chromosomes in amphibians and reptiles may be attributed to the flexibility of their sex determination systems, including the ease of sex reversal in these animals.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Evolution of casein kinase 1 and functional analysis of new doubletime mutants in Drosophila

Circadian clocks are timing devices that rhythmically adjust organism's behavior, physiology, and metabolism to the 24-h day-night cycle. Eukaryotic circadian clocks rely on several interlocked transcription-translation feedback loops, where protein stability is the key part of the delay between transcription and the appearance of the mature proteins within the feedback loops. In bilaterian animals, including mammals and insects, the circadian clock depends on a homologous set of proteins. Despite mostly conserved clock components among the fruit fly Drosophila and mammals, several lineage-specific differences exist. Here we have systematically explored the evolution and sequence variability of insect DBT proteins and their vertebrate homologs casein kinase 1 delta (CKIδ) and epsilon (CKIε), dated the origin and separation of CKIδ from CKIε, and identified at least three additional independent duplications of the CKIδ/ε gene in Petromyzon, Danio, and Xenopus. We determined conserved regions in DBT specific to Diptera, and functionally tested a subset of those in D. melanogaster. Replacement of Lysine K224 with acidic residues strongly impacts the free-running period even in heterozygous flies, whereas homozygous mutants are not viable. K224D mutants have a temperature compensation defect with longer free-running periods at higher temperatures, which is exactly the opposite trend of what was reported for corresponding mammalian mutants. All DBTs of dipteran insects contain the NKRQK motif at positions 220-224. The occurrence of this motif perfectly correlates with the presence of BRIDE OF DOUBLETIME, BDBT, in Diptera. BDBT is a non-canonical FK506-binding protein that physically interacts with Drosophila DBT. The phylogeny of FK506-binding proteins suggests that BDBT is either absent or highly modified in non-dipteran insects. In addition to in silico analysis of DBT/CKIδ/ε evolution and diversity, we have identified four novel casein kinase 1 genes specific to the Drosophila genus.

Evolutionary diversification of the canonical Wnt signaling effector TCF/LEF in chordates

Wnt signaling is essential during animal development and regeneration, but also plays an important role in diseases such as cancer and diabetes. The canonical Wnt signaling pathway is one of the most conserved signaling cascades in the animal kingdom, with the T‐cell factor/lymphoid enhancer factor (TCF/LEF) proteins being the major mediators of Wnt/β‐catenin‐regulated gene expression. In comparison with invertebrates, vertebrates possess a high diversity of TCF/LEF family genes, implicating this as a possible key change to Wnt signaling at the evolutionary origin of vertebrates. However, the precise nature of this diversification is only poorly understood. The aim of this study is to clarify orthology, paralogy, and isoform relationships within the TCF/LEF gene family within chordates via in silico comparative study of TCF/LEF gene structure, molecular phylogeny, and gene synteny. Our results support the notion that the four TCF/LEF paralog subfamilies in jawed vertebrates (gnathostomes) evolved via the two rounds of whole‐genome duplications that occurred during early vertebrate evolution. Importantly, gene structure comparisons and synteny analysis of jawless vertebrate (cyclostome) TCFs suggest that a TCF7L2‐like form of gene structure is a close proxy for the ancestral vertebrate structure. In conclusion, we propose a detailed evolutionary path based on a new pre‐whole‐genome duplication vertebrate TCF gene model. This ancestor gene model highlights the chordate and vertebrate innovations of TCF/LEF gene structure, providing the foundation for understanding the role of Wnt/β‐catenin signaling in vertebrate evolution.

Pharmacological modulation of the cAMP signaling of two isoforms of melanocortin-3 receptor by melanocortin receptor accessory proteins in the tetrapod Xenopus laevis

As a member of the seven-transmembrane rhodopsin-like G protein-coupled receptor superfamily, the melanocortin-3 receptor (MC3R) is vital for the regulation of energy homeostasis and rhythms synchronizing in mammals, and its pharmacological effect could be directly influenced by the presence of melanocortin receptor accessory proteins (MRAPs), MRAP1 and MRAP2. The tetrapod amphibian Xenopus laevis (xl) retains higher duplicated genome than extant teleosts and serves as an ideal model system for embryonic development and physiological studies. However, the melanocortin system of the Xenopus laevis has not yet been thoroughly evaluated. In this work, we performed sequence alignment, phylogenetic tree, and synteny analysis of two xlMC3Rs. Co-immunoprecipitation and immunofluorescence assay further confirmed the co-localization and in vitro interaction of xlMC3Rs with xlMRAPs on the plasma membrane. Our results demonstrated that xlMRAP2.L/S could improve α-MSH-stimulated xlMC3Rs signaling and suppress their surface expression. Moreover, xlMC3R.L showed a similar profile on the ligands and surface expression in the presence of xlMRAP1.L. Overall, the distinct pharmacological modulation of xlMC3R.L and xlMC3R.S by dual MRAP2 proteins elucidated the functional consistency of melanocortin system during genomic duplication of tetrapod vertebrates.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

Comparative Distribution of Repetitive Sequences in the Karyotypes of Xenopus tropicalis and Xenopus laevis (Anura, Pipidae)

Xenopus laevis and its diploid relative, Xenopus tropicalis, are the most used amphibian models. Their genomes have been sequenced, and they are emerging as model organisms for research into disease mechanisms. Despite the growing knowledge on their genomes based on data obtained from massive genome sequencing, basic research on repetitive sequences in these species is lacking. This study conducted a comparative analysis of repetitive sequences in X. laevis and X. tropicalis. Genomic in situ hybridization (GISH) and fluorescence in situ hybridization (FISH) with Cot DNA of both species revealed a conserved enrichment of repetitive sequences at the ends of the chromosomes in these Xenopus species. The repeated sequences located on the short arm of chromosome 3 from X. tropicalis were not related to the sequences on the short arm of chromosomes 3L and 3S from X. laevis, although these chromosomes were homoeologous, indicating that these regions evolved independently in these species. Furthermore, all the other repetitive sequences in X. tropicalis and X. laevis may be species-specific, as they were not revealed in cross-species hybridizations. Painting experiments in X. laevis with chromosome 7 from X. tropicalis revealed shared sequences with the short arm of chromosome 3L. These regions could be related by the presence of the nucleolus organizer region (NOR) in both chromosomes, although the region revealed by chromosome painting in the short arm of chromosome 3L in X. laevis did not correspond to 18S + 28S rDNA sequences, as they did not colocalize. The identification of these repeated sequences is of interest as they provide an explanation to some problems already described in the genome assemblies of these species. Furthermore, the distribution of repetitive DNA in the genomes of X. laevis and X. tropicalis might be a valuable marker to assist us in understanding the genome evolution in a group characterized by numerous polyploidization events coupled with hybridizations.

The repertoire of vertebrate STAT transcription factors: Origin and variations in fish

The stat gene family diversified during early vertebrate evolution thanks to two rounds of whole genome duplication (WGD) to produce a typical repertoire composed of 6 STAT factors (named 1-6). In contrast, only one or two stat genes have been reported in C. elegans and in D. melanogaster. The main types of STAT found from bony fish to mammals are present in Agnathan genomes, but a typical STAT1-6 repertoire is only observed in jawed vertebrates. Comparative syntenies showed that STAT6 was the closest to the ancestor of the family. An extensive survey of stat genes across fish including polyploid species showed that whole genome duplications did not lead to a uniform expansion of stat genes. While 2 to 5 stat1 are present in salmonids, whose genome duplicated about 35My ago, only one copy of stat2 and stat6 is retained. In contrast, common carp, with a recent whole genome duplication (5-10My), possesses a doubled stat repertoire indicating that the elimination of stat2 and stat6 additional copies is not immediate. Altogether our data shed light on the multiplicity of evolutionary pathways followed by key components of the canonical cytokine receptor signalling pathway, and point to differential selective constraints exerted on these factors.

Heterogeneity of synonymous substitution rates in the Xenopus frog genome

With the increasing availability of high quality genomic data, there is opportunity to deeply explore the genealogical relationships of different gene loci between closely related species. In this study, we utilized genomes of Xenopus laevis (XLA, a tetraploid species with (L) and (S) sub-genomes) and X. tropicalis (XTR, a diploid species) to investigate whether synonymous substitution rates among orthologous or homoeologous genes displayed any heterogeneity. From over 1500 orthologous/homoeologous genes collected, we calculated proportion of synonymous substitutions between genomes/sub-genomes (k) and found variation within and between chromosomes. Within most chromosomes, we identified higher k with distance from the centromere, likely attributed to higher substitution rates and recombination in these regions. Using maximum likelihood methods, we identified further evidence supporting rate heterogeneity, and estimated species divergence times and ancestral population sizes. Estimated species divergence times (XLA.L-XLA.S: ~25.5 mya; XLA-XTR: ~33.0 mya) were slightly younger compared to a past study, attributed to consideration of population size in our study. Meanwhile, we found very large estimated population size in the ancestral populations of the two species (N_A = 2.55 x 10⁶). Local hybridization and population structure, which have not yet been well elucidated in frogs, may be a contributing factor to these possible large population sizes.

Activation of mouse NBCe1-B by Xenopus laevis and mouse IRBITs: Role of the variable Nt appendage of IRBITs

The IP3 receptor binding protein released with inositol 1,4,5-trisphosphate (IRBIT) plays important roles in the regulation of intracellular Ca²⁺ signaling and intracellular pH. The mammals express two IRBIT paralogs, i.e., IRBIT1 (encoded by AHCYL1) and IRBIT2 (encoded by AHCYL2). The clawed frog Xenopus laevis oocyte is widely used for biophysical studies on ion channels and transporters. It remains unknown whether endogenous IRBIT is expressed in Xenopus oocytes. Here, we cloned from frog oocyte irbit2.L and irbit2.S, orthologs of mammalian IRBIT2. When over-expressed, the frog IRBITs powerfully stimulate the electrogenic Na⁺/HCO₃^- cotransporter NBCe1-B as mouse IRBIT2-V2 does. Expression of an isolated Nt fragment of NBCe1-B containing the IRBIT-binding domain greatly decreases NBCe1-B activity in oocytes, suggesting that the basal activity of NBCe1-B contains a large component derived from the stimulation by endogenous frog IRBIT. The frog IRBITs are highly homologous to the mammalian ones in the carboxyl-terminal region, but varies greatly in the amino-terminal (Nt) appendage. Interestingly, truncation study showed that the Nt appendage of IRBIT1 and the long Nt appendage of IRBIT2-V2 modestly enhance, whereas the short Nt appendage of IRBIT2-V4 greatly inhibits the functional interaction between IRBIT and NBCe1-B. Finally, Ala-substitution of Ser68, a key phosphorylation site in the PEST domain of IRBIT, causes distinct functional consequences depending on the structural context of the Nt appendage in different IRBIT isoforms. We conclude that the Nt appendage of IRBITs is not necessary for, but plays an important regulatory role in the functional interaction between IRBIT and NBCe1-B.

Inferring the "Primordial Immune Complex": Origins of MHC Class I and Antigen Receptors Revealed by Comparative Genomics

Comparative analyses suggest that the MHC was derived from a prevertebrate "primordial immune complex" (PIC). PIC duplicated twice in the well-studied two rounds of genome-wide duplications (2R) early in vertebrate evolution, generating four MHC paralogous regions (predominantly on human chromosomes [chr] 1, 6, 9, 19). Examining chiefly the amphibian Xenopus laevis, but also other vertebrates, we identified their MHC paralogues and mapped MHC class I, AgR, and "framework" genes. Most class I genes mapped to MHC paralogues, but a cluster of Xenopus MHC class Ib genes (xnc), which previously was mapped outside of the MHC paralogues, was surrounded by genes syntenic to mammalian CD1 genes, a region previously proposed as an MHC paralogue on human chr 1. Thus, this gene block is instead the result of a translocation that we call the translocated part of the MHC paralogous region (MHCtrans) Analyses of Xenopus class I genes, as well as MHCtrans, suggest that class I arose at 1R on the chr 6/19 ancestor. Of great interest are nonrearranging AgR-like genes mapping to three MHC paralogues; thus, PIC clearly contained several AgR precursor loci, predating MHC class I/II. However, all rearranging AgR genes were found on paralogues derived from the chr 19 precursor, suggesting that invasion of a variable (V) exon by the RAG transposon occurred after 2R. We propose models for the evolutionary history of MHC/TCR/Ig and speculate on the dichotomy between the jawless (lamprey and hagfish) and jawed vertebrate adaptive immune systems, as we found genes related to variable lymphocyte receptors also map to MHC paralogues.

A chromosome-scale genome assembly and dense genetic map for Xenopus tropicalis

The Western clawed frog Xenopus tropicalis is a diploid model system for both frog genetics and developmental biology, complementary to the paleotetraploid X. laevis. Here we report a chromosome-scale assembly of the X. tropicalis genome, improving the previously published draft genome assembly through the use of new assembly algorithms, additional sequence data, and the addition of a dense genetic map. The improved genome enables the mapping of specific traits (e.g., the sex locus or Mendelian mutants) and the characterization of chromosome-scale synteny with other tetrapods. We also report an improved annotation of the genome that integrates deep transcriptome sequence from diverse tissues and stages. The exon-intron structures of these genes are highly conserved relative to both X. laevis and human, as are chromosomal linkages ("synteny") and local gene order. A network analysis of developmental gene expression will aid future studies.

Silurana Chromosomal Evolution: A New Piece to the Puzzle

The African clawed frogs of the subgenus Silurana comprise both diploid and tetraploid species. The root of the polyploidization event leading to the extant Xenopus calcaratus, X. mellotropicalis, and X. epitropicalis is not fully understood so far. In X. mellotropicalis, we previously proposed 2 evolutionary scenarios encompassing complete (scenario A) or incomplete (scenario B) translocation of a heterochromatic block from chromosome 9 to 2 in a diploid ancestor. To resolve this puzzle, we performed FISH coupled with tyramide signal amplification (FISH-TSA) using 5 X. tropicalis and X. mellotropicalis single copy gene probes (gyg2, cept1, fn1, ndufs1, and sf3b1) reflecting borders of the heterochromatic blocks in X. tropicalis chromosome 9 (XTR 9) and X. mellotropicalis chromosome 9b (XME 9b) and XME 2a. cDNA sequencing recognized both homoeologous genes in X. mellotropicalis. Comparison of gene physical mapping between X. tropicalis and X. mellotropicalis clearly confirmed complete rather than incomplete translocation t(9;2) of the heterochromatic block in the diploid predecessor and thus favored scenario A regarding the formation of an ancestral allotetraploid karyotype.

Expanding the genetic toolkit in Xenopus: Approaches and opportunities for human disease modeling

The amphibian model Xenopus, has been used extensively over the past century to study multiple aspects of cell and developmental biology. Xenopus offers advantages of a non-mammalian system, including high fecundity, external development, and simple housing requirements, with additional advantages of large embryos, highly conserved developmental processes, and close evolutionary relationship to higher vertebrates. There are two main species of Xenopus used in biomedical research, Xenopus laevis and Xenopus tropicalis; the common perception is that both species are excellent models for embryological and cell biological studies, but only Xenopus tropicalis is useful as a genetic model. The recent completion of the Xenopus laevis genome sequence combined with implementation of genome editing tools, such as TALENs (transcription activator-like effector nucleases) and CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR associated nucleases), greatly facilitates the use of both Xenopus laevis and Xenopus tropicalis for understanding gene function in development and disease. In this paper, we review recent advances made in Xenopus laevis and Xenopus tropicalis with TALENs and CRISPR-Cas and discuss the various approaches that have been used to generate knockout and knock-in animals in both species. These advances show that both Xenopus species are useful for genetic approaches and in particular counters the notion that Xenopus laevis is not amenable to genetic manipulations.

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.

Congruent population structure across paralogous and nonparalogous loci in Salish Sea chum salmon (Oncorhynchus keta)

Whole-genome duplications are major evolutionary events with a lasting impact on genome structure. Duplication events complicate genetic analyses as paralogous sequences are difficult to distinguish; consequently, paralogs are often excluded from studies. The effects of an ancient whole-genome duplication (approximately 88 MYA) are still evident in salmonids through the persistence of numerous paralogous gene sequences and partial tetrasomic inheritance. We use restriction site-associated DNA sequencing on 10 collections of chum salmon from the Salish Sea in the USA and Canada to investigate genetic diversity and population structure in both tetrasomic and rediploidized regions of the genome. We use a pedigree and high-density linkage map to identify paralogous loci and to investigate genetic variation across the genome. By applying multivariate statistical methods, we show that it is possible to characterize paralogous loci and that they display similar patterns of population structure as the diploidized portion of the genome. We find genetic associations with the adaptively important trait of run-timing in both sets of loci. By including paralogous loci in genome scans, we can observe evolutionary signals in genomic regions that have routinely been excluded from population genetic studies in other polyploid-derived species.

Chromosome divergence during evolution of the tetraploid clawed frogs, Xenopus mellotropicalis and Xenopus epitropicalis as revealed by Zoo-FISH

Whole genome duplication (WGD) generates new species and genomic redundancy. In African clawed frogs of the genus Xenopus, this phenomenon has been especially important in that (i) all but one extant species are polyploid and (ii) whole genome sequences of some species provide an evidence for genomic rearrangements prior to or after WGD. Within Xenopus in the subgenus Silurana, at least one allotetraploidization event gave rise to three extant tetraploid (2n = 4x = 40) species–Xenopus mellotropicalis, X. epitropicalis, and X. calcaratus–but it is not yet clear the degree to which these tetraploid genomes experienced rearrangements prior to or after allotetraploidization. To explore genome evolution during diversification of these species, we performed cytogenetic analyses of X. mellotropicalis, including assessment of the localization of nucleolar organizer region, chromosome banding, and determination of the p/q arm ratios for each chromosome pair. We compared these data to a previously characterized karyotype of X. epitropicalis. Morphometric, C-banding and Zoo-FISH data support a previously hypothesized common allotetraploid predecessor of these species. Zoo-FISH with whole chromosome painting (WCP) probes derived from the closely related diploid species X. tropicalis confirmed the existence of ten chromosomal quartets in X. mellotropicalis somatic cells, as expected by its ploidy level and tetraploid ancestry. The p/q arm ratio of chromosome 2a was found to be substantially different between X. mellotropicalis (0.81) and X. epitropicalis (0.67), but no substantial difference between these two species was detected in this ratio for the homoeologous chromosome pair 2b, or for other chromosome pairs. Additionally, we identified variation between these two species in the locations of a heterochromatic block on chromosome pair 2a. These results are consistent with a dynamic history of genomic rearrangements before and/or after genome duplication, a surprising finding given the otherwise relatively conserved genomic structure of most frogs.

Co-accumulation of cis-regulatory and coding mutations during the pseudogenization of the Xenopus laevis homoeologs six6.L and six6.S

Common models for the evolution of duplicated genes after genome duplication are subfunctionalization, neofunctionalization, and pseudogenization. Although the crucial roles of cis-regulatory mutations in subfunctionalization are well-documented, their involvement in pseudogenization and/or neofunctionalization remains unclear. We addressed this issue by investigating the evolution of duplicated homeobox genes, six6.L and six6.S, in the allotetraploid frog Xenopus laevis. Based on a comparative expression analysis, we observed similar eye-specific expression patterns for the two loci and their single ortholog in the ancestral-type diploid species Xenopus tropicalis. However, we detected lower levels of six6.S expression than six6.L expression. The six6.S enhancer sequence was more highly diverged from the orthologous enhancer of X. tropicalis than the six6.L enhancer, and showed weaker activity in a transgenic reporter assay. Based on a phylogenetic analysis of the protein sequences, we observed greater divergence between X. tropicalis Six6 and Six6.S than between X. tropicalis Six6 and Six6.L, and the observed mutations were reminiscent of a microphthalmia mutation in human SIX6. Misexpression experiments showed that six6.S has weaker eye-enlarging activity than six6.L, and targeted disruption of six6.L reduced the eye size more significantly than that of six6.S. These results suggest that enhancer attenuation stimulates the accumulation of hypomorphic coding mutations, or vice versa, in one duplicated gene copy and facilitates pseudogenization. We also underscore the value of the allotetraploid genome of X. laevis as a resource for studying latent pathogenic mutations.

Xenopus Piwi proteins interact with a broad proportion of the oocyte transcriptome

Piwi proteins utilize small RNAs (piRNAs) to recognize target transcripts such as transposable elements (TE). However, extensive piRNA sequence diversity also suggests that Piwi/piRNA complexes interact with many transcripts beyond TEs. To determine Piwi target RNAs, we used ribonucleoprotein-immunoprecipitation (RIP) and cross-linking and immunoprecipitation (CLIP) to identify thousands of transcripts associated with the Piwi proteins XIWI and XILI (Piwi-protein-associated transcripts, PATs) from early stage oocytes of X. laevis and X. tropicalis Most PATs associate with both XIWI and XILI and include transcripts of developmentally important proteins in oogenesis and embryogenesis. Only a minor fraction of PATs in both frog species displayed near perfect matches to piRNAs. Since predicting imperfect pairing between all piRNAs and target RNAs remains intractable, we instead determined that PAT read counts correlate well with the lengths and expression levels of transcripts, features that have also been observed for oocyte mRNAs associated with Drosophila Piwi proteins. We used an in vitro assay with exogenous RNA to confirm that XIWI associates with RNAs in a length- and concentration-dependent manner. In this assay, noncoding transcripts with many perfectly matched antisense piRNAs were unstable, whereas coding transcripts with matching piRNAs were stable, consistent with emerging evidence that Piwi proteins both promote the turnover of TEs and other RNAs, and may also regulate mRNA localization and translation. Our study suggests that Piwi proteins play multiple roles in germ cells and establishes a tractable vertebrate system to study the role of Piwi proteins in transcript regulation.

Sequential Turnovers of Sex Chromosomes in African Clawed Frogs (Xenopus) Suggest Some Genomic Regions Are Good at Sex Determination

Sexual differentiation is fundamentally important for reproduction, yet the genetic triggers of this developmental process can vary, even between closely related species. Recent studies have uncovered, for example, variation in the genetic triggers for sexual differentiation within and between species of African clawed frogs (genus Xenopus). Here, we extend these discoveries by demonstrating that yet another sex determination system exists in Xenopus, specifically in the species Xenopus borealis. This system evolved recently in an ancestor of X. borealis that had the same sex determination system as X. laevis, a system which itself is newly evolved. Strikingly, the genomic region carrying the sex determination factor in X. borealis is homologous to that of therian mammals, including humans. Our results offer insights into how the genetic underpinnings of conserved phenotypes evolve, and suggest an important role for cooption of genetic building blocks with conserved developmental roles.

Genome evolution in the allotetraploid frog Xenopus laevis

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.

Mechanisms of transcription factor evolution in Metazoa

Transcriptions factors (TFs) are pivotal for the regulation of virtually all cellular processes, including growth and development. Expansions of TF families are causally linked to increases in organismal complexity. Here we study the evolutionary dynamics, genetic causes and functional implications of the five largest metazoan TF families. We find that family expansions dominate across the whole metazoan tree; however, some branches experience exceptional family-specific accelerated expansions. Additionally, we find that such expansions are often predated by modular domain rearrangements, which spur the expansion of a new sub-family by separating it from the rest of the TF family in terms of protein-protein interactions. This separation allows for radical shifts in the functional spectrum of a duplicated TF. We also find functional differentiation inside TF sub-families as changes in expression specificity. Furthermore, accelerated family expansions are facilitated by repeats of sequence motifs such as C2H2 zinc fingers. We quantify whole genome duplications and single gene duplications as sources of TF family expansions, implying that some, but not all, TF duplicates are preferentially retained. We conclude that trans-regulatory changes (domain rearrangements) are instrumental for fundamental functional innovations, that cis-regulatory changes (affecting expression) accomplish wide-spread fine tuning and both jointly contribute to the functional diversification of TFs.

Expansion of amphibian intronless interferons revises the paradigm for interferon evolution and functional diversity

Interferons (IFNs) are key cytokines identified in vertebrates and evolutionary dominance of intronless IFN genes in amniotes is a signature event in IFN evolution. For the first time, we show that the emergence and expansion of intronless IFN genes is evident in amphibians, shown by 24–37 intronless IFN genes in each frog species. Amphibian IFNs represent a molecular complex more complicated than those in other vertebrate species, which revises the established model of IFN evolution to facilitate re-inspection of IFN molecular and functional diversity. We identified these intronless amphibian IFNs and their intron-containing progenitors, and functionally characterized constitutive and inductive expression and antimicrobial roles in infections caused by zoonotic pathogens, such as influenza viruses and Listeria monocytogenes. Amphibians, therefore, may serve as overlooked vectors/hosts for zoonotic pathogens, and the amphibian IFN system provides a model to study IFN evolution in molecular and functional diversity in coping with dramatic environmental changes during terrestrial adaption.

A linkage map for the Newt Notophthalmus viridescens: Insights in vertebrate genome and chromosome evolution

Genetic linkage maps are fundamental resources that enable diverse genetic and genomic approaches, including quantitative trait locus (QTL) analyses and comparative studies of genome evolution. It is straightforward to build linkage maps for species that are amenable to laboratory culture and genetic crossing designs, and that have relatively small genomes and few chromosomes. It is more difficult to generate linkage maps for species that do not meet these criteria. Here, we introduce a method to rapidly build linkage maps for salamanders, which are known for their enormous genome sizes. As proof of principle, we developed a linkage map with thousands of molecular markers (N=2349) for the Eastern newt (Notophthalmus viridescens). The map contains 12 linkage groups (152.3-934.7cM), only one more than the number of chromosome pairs. Importantly, this map was generated using RNA isolated from a single wild caught female and her 28 offspring. We used the map to reveal chromosome-scale conservation of synteny among N. viridescens, A. mexicanum (Urodela), and chicken (Amniota), and to identify large conserved segments between N. viridescens and Xenopus tropicalis (Anura). We also show that met1, a major effect QTL that regulates the expression of alternate metamorphic and paedomorphic modes of development in Ambystoma, associates with a chromosomal fusion that is not found in the N. viridescens map. Our results shed new light on the ancestral amphibian karyotype and reveal specific fusion and translocation events that shaped the genomes of three amphibian model taxa. The ability to rapidly build linkage maps for large salamander genomes will enable genetic and genomic analyses within this important vertebrate group, and more generally, empower comparative studies of vertebrate biology and evolution.

Open and closed evolutionary paths for drastic morphological changes, involving serial gene duplication, sub-functionalization, and selection

Twin-tail goldfish strains are examples of drastic morphological alterations that emerged through domestication. Although this mutation is known to be caused by deficiency of one of two duplicated chordin genes, it is unknown why equivalent mutations have not been observed in other domesticated fish species. Here, we compared the chordin gene morphant phenotypes of single-tail goldfish and common carp (close relatives, both of which underwent chordin gene duplication and domestication). Morpholino-induced knockdown depleted chordin gene expression in both species; however, while knockdown reproduced twin-tail morphology in single-tail goldfish, it had no effect on common carp morphology. This difference can be explained by the observation that expression patterns of the duplicated chordin genes overlap completely in common carp, but are sub-functionalized in goldfish. Our finding implies that goldfish drastic morphological changes might be enhanced by the subsequent occurrence of three different types of evolutionary event (duplication, sub-functionalization, and selection) in a certain order.

Goldfish morphology as a model for evolutionary developmental biology

Morphological variation of the goldfish is known to have been established by artificial selection for ornamental purposes during the domestication process. Chinese texts that date to the Song dynasty contain descriptions of goldfish breeding for ornamental purposes, indicating that the practice originated over one thousand years ago. Such a well-documented goldfish breeding process, combined with the phylogenetic and embryological proximities of this species with zebrafish, would appear to make the morphologically diverse goldfish strains suitable models for evolutionary developmental (evodevo) studies. However, few modern evodevo studies of goldfish have been conducted. In this review, we provide an overview of the historical background of goldfish breeding, and the differences between this teleost and zebrafish from an evolutionary perspective. We also summarize recent progress in the field of molecular developmental genetics, with a particular focus on the twin-tail goldfish morphology. Furthermore, we discuss unanswered questions relating to the evolution of the genome, developmental robustness, and morphologies in the goldfish lineage, with the goal of blazing a path toward an evodevo study paradigm using this teleost species as a new model species. For further resources related to this article, please visit the WIREs website.

Towards the bridging of molecular genetics data across Xenopus species

BACKGROUND

The clawed African frog Xenopus laevis has been one of the main vertebrate models for studies in developmental biology. However, for genetic studies, Xenopus tropicalis has been the experimental model of choice because it shorter life cycle and due to a more tractable genome that does not result from genome duplication as in the case of X. laevis. Today, although still organized in a large number of scaffolds, nearly 85% of X. tropicalis and 89% of X. laevis genomes have been sequenced. There is expectation for a comparative physical map that can be used as a Rosetta Stone between X. laevis genetic studies and X. tropicalis genomic research.

RESULTS

In this work, we have mapped using coarse-grained alignment the 18 chromosomes of X. laevis, release 9.1, on the 10 reference scaffolds representing the haploid genome of X. tropicalis, release 9.0. After validating the mapping with theoretical data, and estimating reference averages of genome sequence identity, 37 to 44% between the two species, we have carried out a synteny analysis for 2,112 orthologous genes. We found that 99.6% of genes are in the same organization.

CONCLUSIONS

Taken together, our results make possible to establish the correspondence between 62 and 65.5% of both genomes, percentage of identity, synteny and automatic annotation of transcripts of both species, providing a new and more comprehensive tool for comparative analysis of these two species, by allowing to bridge molecular genetics data among them.

Nucleolar dominance and maternal control of 45S rDNA expression

Using a system of interspecies hybrids, trihybrids, and recombinants with varying proportions of genomes from three distinct Xenopus species, we provide evidence for de novo epigenetic silencing of paternal 45 S ribosomal ribonucleic acid (rRNA) genes and their species-dependent expression dominance that escapes transcriptional inactivation after homologous recombination. The same pattern of imprinting is maintained in the offspring from mothers being genetic males (ZZ) sex-reversed to females, indicating that maternal control of ribosomal deoxyribonucleic acid (rDNA) expression is not sex-chromosome linked. Nucleolar dominance (nucleolus underdevelopment) in Xenopus hybrids appears to be associated with a major non-Mendelian reduction in the number of 45 S rDNA gene copies rather than a specific pattern of their expression. The loss of rRNA gene copies in F1 hybrids was non-random with respect to the parental species, with the transcriptionally dominant variant preferentially removed from hybrid zygotes. This dramatic disruption in the structure and function of 45 S rDNA impacts transcriptome patterns of small nucleolar RNAs and messenger RNAs, with genes from the ribosome and oxidative stress pathways being among the most affected. Unorthodoxies of rDNA inheritance and expression may be interpreted as hallmarks of genetic conflicts between parental genomes, as well as defensive epigenetic mechanisms employed to restore genome integrity.

The chromosomal constitution of fish hybrid lineage revealed by 5S rDNA FISH

Background

The establishment of the bisexual fertile fish hybrid lineage including the allodiploid and allotetraploid hybrids, from interspecific hybridization of red crucian carp (Carassius auratus red var. 2n = 100, 2n = AA) (♀) × common carp (Cyprinus carpio L. 2n = 100, 2n = BB) (♂), provided a good platform to investigate genetic relationship between the parents and their hybrid progenies.

Results

The chromosomal inheritance of diploid and allotetraploid hybrid progenies in successive generations, was studied by applying 5S rDNA fluorescence in situ hybridization. Signals of 5S rDNA distinguished the chromosomal constitution of common carp (B-genome) from red crucian carp (A-genome), in which two strong signals were observed on the first submetacentric chromosome, while no major signal was found in common carp. After fish hybridization, one strong signal of 5S rDNA was detected in the same locus on the chromosome of diploid hybrids. As expected, two strong signals were observed in 4nF₃ tetraploid hybrids offspring and it is worth mentioning that two strong signals were detected in a separating bivalent of a primary spermatocyte in 4nF₃. Furthermore, the mitosis of heterozygous chromosomes was shown normal and stable with blastular tissue histological studies.

Conclusions

We revealed that 5S rDNA signal can be applied to discern A-genome from B-genome, and that 5S rDNA bearing chromosomes can be stably passed down in successive generations. Our work provided a significant method in fish breeding and this is important for studies in fish evolutionary biology.

The Xenopus ORFeome: A resource that enables functional genomics

Functional characterisation of proteins and large-scale, systems-level studies are enabled by extensive sets of cloned open reading frames (ORFs) in an easily-accessible format that enables many different applications. Here we report the release of the first stage of the Xenopus ORFeome, which contains 8673 ORFs from the Xenopus Gene Collection (XGC) for Xenopus laevis, cloned into a Gateway® donor vector enabling rapid in-frame transfer of the ORFs to expression vectors. This resource represents an estimated 7871 unique genes, approximately 40% of the non-redundant X. laevis gene complement, and includes 2724 genes where the human ortholog has an association with disease. Transfer into the Gateway system was validated by 5' and 3' end sequencing of the entire collection and protein expression of a set of test clones. In a parallel process, the underlying ORF predictions from the original XGC collection were re-analysed to verify quality and full-length status, identifying those proteins likely to exhibit truncations when translated. These data are integrated into Xenbase, the Xenopus community database, which associates genomic, expression, function and human disease model metadata to each ORF, enabling end-users to search for ORFeome clones with links to commercial distributors of the collection. When coupled with the experimental advantages of Xenopus eggs and embryos, the ORFeome collection represents a valuable resource for functional genomics and disease modelling.

Genome-Wide Identification and Expression of Xenopus F-Box Family of Proteins

Protein degradation via the multistep ubiquitin/26S proteasome pathway is a rapid way to alter the protein profile and drive cell processes and developmental changes. Many key regulators of embryonic development are targeted for degradation by E3 ubiquitin ligases. The most studied family of E3 ubiquitin ligases is the SCF ubiquitin ligases, which use F-box adaptor proteins to recognize and recruit target proteins. Here, we used a bioinformatics screen and phylogenetic analysis to identify and annotate the family of F-box proteins in the Xenopus tropicalis genome. To shed light on the function of the F-box proteins, we analyzed expression of F-box genes during early stages of Xenopus development. Many F-box genes are broadly expressed with expression domains localized to diverse tissues including brain, spinal cord, eye, neural crest derivatives, somites, kidneys, and heart. All together, our genome-wide identification and expression profiling of the Xenopus F-box family of proteins provide a foundation for future research aimed to identify the precise role of F-box dependent E3 ubiquitin ligases and their targets in the regulatory circuits of development.

The sea lamprey meiotic map improves resolution of ancient vertebrate genome duplications

It is generally accepted that many genes present in vertebrate genomes owe their origin to two whole-genome duplications that occurred deep in the ancestry of the vertebrate lineage. However, details regarding the timing and outcome of these duplications are not well resolved. We present high-density meiotic and comparative genomic maps for the sea lamprey (Petromyzon marinus), a representative of an ancient lineage that diverged from all other vertebrates ∼550 million years ago. Linkage analyses yielded a total of 95 linkage groups, similar to the estimated number of germline chromosomes (1n ∼ 99), spanning a total of 5570.25 cM. Comparative mapping data yield strong support for the hypothesis that a single whole-genome duplication occurred in the basal vertebrate lineage, but do not strongly support a hypothetical second event. Rather, these comparative maps reveal several evolutionarily independent segmental duplications occurring over the last 600+ million years of chordate evolution. This refined history of vertebrate genome duplication should permit more precise investigations of vertebrate evolution.

Inbreeding Ratio and Genetic Relationships among Strains of the Western Clawed Frog, Xenopus tropicalis

The Western clawed frog, Xenopus tropicalis, is a highly promising model amphibian, especially in developmental and physiological research, and as a tool for understanding disease. It was originally found in the West African rainforest belt, and was introduced to the research community in the 1990s. The major strains thus far known include the Nigerian and Ivory Coast strains. However, due to its short history as an experimental animal, the genetic relationship among the various strains has not yet been clarified, and establishment of inbred strains has not yet been achieved. Since 2003 the Institute for Amphibian Biology (IAB), Hiroshima University has maintained stocks of multiple X. tropicalis strains and conducted consecutive breeding as part of the National BioResource Project. In the present study we investigated the inbreeding ratio and genetic relationship of four inbred strains at IAB, as well as stocks from other institutions, using highly polymorphic microsatellite markers and mitochondrial haplotypes. Our results show successive reduction of heterozygosity in the genome of the IAB inbred strains. The Ivory Coast strains clearly differed from the Nigerian strains genetically, and three subgroups were identified within both the Nigerian and Ivory Coast strains. It is noteworthy that the Ivory Coast strains have an evolutionary divergent genetic background. Our results serve as a guide for the most effective use of X. tropicalis strains, and the long-term maintenance of multiple strains will contribute to further research efforts.

Chromosome Banding in Amphibia. XXXII. The Genus Xenopus (Anura, Pipidae)

Mitotic chromosomes of 16 species of the frog genus Xenopus were prepared from kidney and lung cell cultures. In the chromosomes of 7 species, high-resolution replication banding patterns could be induced by treating the cultures with 5-bromodeoxyuridine (BrdU) and deoxythymidine (dT) in succession, and in 6 of these species the BrdU/dT-banded chromosomes could be arranged into karyotypes. In the 3 species of the clade with 2n = 20 and 4n = 40 chromosomes (X. tropicalis, X. epitropicalis, X. new tetraploid 1), as well as in the 3 species with 4n = 36 chromosomes (X. laevis, X. borealis, X. muelleri), the BrdU/dT-banded karyotypes show a high degree of homoeology, though differences were detected between these groups. Translocations, inversions, insertions or sex-specific replication bands were not observed. Minor replication asynchronies found between chromosomes probably involve heterochromatic regions. BrdU/dT replication banding of Xenopus chromosomes provides the landmarks necessary for the exact physical mapping of genes and repetitive sequences. FISH with an X. laevis 5S rDNA probe detected multiple hybridization sites at or near the long-arm telomeric regions in most chromosomes of X. laevis and X. borealis, whereas in X. muelleri, the 5S rDNA sequences are located exclusively at the long-arm telomeres of a single chromosome pair. Staining with the AT base pair-specific fluorochrome quinacrine mustard revealed brightly fluorescing heterochromatic regions in the majority of X. borealis chromosomes which are absent in other Xenopus species.

Polyploidy in Amphibia

This review summarizes the current status of the known extant genuine polyploid anuran and urodelan species, as well as spontaneously originated and/or experimentally produced amphibian polyploids. The mechanisms by which polyploids can originate, the meiotic pairing configurations, the diploidization processes operating in polyploid genomes, the phenomenon of hybridogenesis, and the relationship between polyploidization and sex chromosome evolution are discussed. The polyploid systems in some important amphibian taxa are described in more detail.

Benchmarking Transcriptome Quantification Methods for Duplicated Genes in Xenopus laevis

Xenopus is an important model organism for the study of genome duplication in vertebrates. With the full genome sequence of diploid Xenopus tropicalis available, and that of allotetraploid X. laevis close to being finished, we will be able to expand our understanding of how duplicated genes have evolved. One of the key features in the study of the functional consequence of gene duplication is how their expression patterns vary across different conditions, and RNA-seq seems to have enough resolution to discriminate the expression of highly similar duplicated genes. However, most of the current RNA-seq analysis methods were not designed to study samples with duplicate genes such as in X. laevis. Here, various computational methods to quantify gene expression in RNA-seq data were evaluated, using 2 independent X. laevis egg RNA-seq datasets and 2 reference databases for duplicated genes. The fact that RNA-seq can measure expression levels of similar duplicated genes was confirmed, but long paired-end reads are more informative than short single-end reads to discriminate duplicated genes. Also, it was found that bowtie, one of the most popular mappers in RNA-seq analysis, reports significantly smaller numbers of unique hits according to a mapping quality score compared to other mappers tested (BWA, GSNAP, STAR). Calculated from unique hits based on a mapping quality score, both expression levels and the expression ratio of duplicated genes can be estimated consistently among biological replicates, demonstrating that this method can successfully discriminate the expression of each copy of a duplicated gene pair. This comprehensive evaluation will be a useful guideline for studying gene expression of organisms with genome duplication using RNA-seq in the future.

Xenopus Cytogenetics and Chromosomal Evolution

The genus Xenopus represents important model organisms in the field of developmental biology and chromosomal evolution. Developmental processes are tightly coupled with the analysis of gene function via genetic linkage and mapping. Cytogenetic techniques such as chromosome banding or FISH are essential tools for the determination of gene position and subsequently for the construction of linkage and physical maps. Here, we present a summary of key achievements in X. tropicalis and X. laevis cytogenetics with emphasis on the gene localization to chromosomes. The second part of this review is focused on the chromosomal evolution regarding both above-mentioned species. With respect to methodology, hybridization techniques such as FISH and chromosome-specific painting FISH are highlighted.

Homeolog-specific targeted mutagenesis in Xenopus laevis using TALENs

Transcription activator-like effector nucleases (TALENs) have previously been used for targeted genome editing in various organisms including Xenopus laevis. However, because of genomic polyploidization, X. laevis usually possess homeologous genes (homeologs) with quite similar sequences that make the analysis of gene function difficult. In the present study, we show methodological examples of targeted gene modification of X. laevis homeologs. The X. laevis cytoglobin gene (cygb) consists of two homeologs (xlcygba and xlcygbb), and molecular phylogenetic analysis suggested that they have potentially different functions. Thus, there is a need to establish a method of homeolog-specific gene disruption to clarify gene functions in detail. Here, we show successful examples of homeolog-specific and simultaneous gene disruption for xlcygba and xlcygbb. We found that selective digestion can be performed with at least three mismatches in TALEN target sites in both homeologs. This report paves the way for the functional analyses of X. laevis homeologs, even those containing nearly identical sequences.

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.

Xenopus as a Model for GI/Pancreas Disease

Diseases affecting endodermal organs like the pancreas, lung and gastrointestinal (GI) tract have a substantial impact on human welfare. Since many of these are congenital defects that arise as a result of defects during development broad efforts are focused on understanding the development of these organs so as to better identify risk factors, disease mechanisms and therapeutic targets. Studies implementing model systems, like the amphibian Xenopus, have contributed immensely to our understanding of signaling (e.g. Wnt, FGF, BMP, RA) pathways and gene regulation (e.g. hhex, ptf1a, ngn3) that underlie normal development as well as disease progression. Recent advances in genome engineering further enhance the capabilities of the Xenopus model system for pursuing biomedical research, and will undoubtedly result in a boom of new information underlying disease mechanisms ultimately leading to advancements in diagnosis and therapy.

Biased gene conversion and GC-content evolution in the coding sequences of reptiles and vertebrates

Mammalian and avian genomes are characterized by a substantial spatial heterogeneity of GC-content, which is often interpreted as reflecting the effect of local GC-biased gene conversion (gBGC), a meiotic repair bias that favors G and C over A and T alleles in high-recombining genomic regions. Surprisingly, the first fully sequenced nonavian sauropsid (i.e., reptile), the green anole Anolis carolinensis, revealed a highly homogeneous genomic GC-content landscape, suggesting the possibility that gBGC might not be at work in this lineage. Here, we analyze GC-content evolution at third-codon positions (GC3) in 44 vertebrates species, including eight newly sequenced transcriptomes, with a specific focus on nonavian sauropsids. We report that reptiles, including the green anole, have a genome-wide distribution of GC3 similar to that of mammals and birds, and we infer a strong GC3-heterogeneity to be already present in the tetrapod ancestor. We further show that the dynamic of coding sequence GC-content is largely governed by karyotypic features in vertebrates, notably in the green anole, in agreement with the gBGC hypothesis. The discrepancy between third-codon positions and noncoding DNA regarding GC-content dynamics in the green anole could not be explained by the activity of transposable elements or selection on codon usage. This analysis highlights the unique value of third-codon positions as an insertion/deletion-free marker of nucleotide substitution biases that ultimately affect the evolution of proteins.

Gtpbp2 is required for BMP signaling and mesoderm patterning in Xenopus embryos

Smad proteins convey canonical intracellular signals for activated receptors in the TGFβ superfamily, but the activity of Smads and their impact on target genes are further regulated by a wide variety of cofactors and partner proteins. We have identified a new Smad1 partner, a GTPase named Gtpbp2 that is a distant relative of the translation factor eEf1a. Gtpbp2 affects canonical signaling in the BMP branch of the TGFβ superfamily, as morpholino knockdown of Gtpbp2 decreases, and overexpression of Gtpbp2 enhances, animal cap responses to BMP4. During Xenopus development, gtpbp2 transcripts are maternally expressed and localized to the egg animal pole, and partitioned into the nascent ectodermal and mesodermal cells during cleavage and early gastrulation stages. Subsequently, gtpbp2 is expressed in the neural folds, and in early tadpoles undergoing organogenesis gtpbp2 is expressed prominently in the brain, eyes, somites, ventral blood island and branchial arches. Consistent with its expression, morpholino knockdown of Gtpbp2 causes defects in ventral-posterior germ layer patterning, gastrulation and tadpole morphology. Overexpressed Gtpbp2 can induce ventral-posterior marker genes and localize to cell nuclei in Xenopus animal caps, highlighting its role in regulating BMP signaling in the early embryo. Here, we introduce this large GTPase as a novel factor in BMP signaling and ventral-posterior patterning.

Engineering Xenopus embryos for phenotypic drug discovery screening

Many rare human inherited diseases remain untreatable despite the fact that the disease causing genes are known and adequate mouse disease models have been developed. In vivo phenotypic drug screening relies on isolating drug candidates by their ability to produce a desired therapeutic phenotype in whole organisms. Embryos of zebrafish and Xenopus frogs are abundant, small and free-living. They can be easily arrayed in multi-well dishes and treated with small organic molecules. With the development of novel genome modification tools, such a zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and CRISPR/Cas, it is now possible to efficiently engineer non-mammalian models of inherited human diseases. Here, we will review the rapid progress made in adapting these novel genome editing tools to Xenopus. The advantages of Xenopus embryos as in vivo models to study human inherited diseases will be presented and their utility for drug discovery screening will be discussed. Being a tetrapod, Xenopus complements zebrafish as an indispensable non-mammalian animal model for the study of human disease pathologies and the discovery of novel therapeutics for inherited diseases.

Targeted mutagenesis of multiple and paralogous genes in Xenopus laevis using two pairs of transcription activator-like effector nucleases

Transcription activator-like effector nucleases (TALENs) have been extensively used in genome editing in various organisms. In some cases, however, it is difficult to efficiently disrupt both paralogous genes using a single pair of TALENs in Xenopus laevis because of its polyploidy. Here, we report targeted mutagenesis of multiple and paralogous genes using two pairs of TALENs in X. laevis. First, we show simultaneous targeted mutagenesis of three genes, tyrosinase paralogues (tyra and tyrb) and enhanced green fluorescent protein (egfp) by injection of two TALENs pairs in transgenic embryos carrying egfp. Consistent with the high frequency of both severe phenotypic traits, albinism and loss of GFP fluorescence, frameshift mutation rates of tyr paralogues and egfp reached 40-80%. Next, we show early introduction of TALEN-mediated mutagenesis of these target loci during embryogenesis. Finally, we also demonstrate that two different pairs of TALENs can simultaneously introduce mutations to both paralogues encoding histone chaperone with high efficiency. Our results suggest that targeted mutagenesis of multiple genes using TALENs can be applied to analyze the functions of paralogous genes with redundancy in X. laevis.

Cytogenetic analysis of Phyllomedusa distincta Lutz, 1950 (2n = 2x = 26), P. tetraploidea Pombal and Haddad, 1992 (2n = 4x = 52), and their natural triploid hybrids (2n = 3x = 39) (Anura, Hylidae, Phyllomedusinae)

Background

Natural polyploidy has played an important role during the speciation and evolution of vertebrates, including anurans, with more than 55 described cases. The species of the Phyllomedusa burmeisteri group are mostly characterized by having 26 chromosomes, but a karyotype with 52 chromosomes was described in P. tetraploidea. This species was found in sintopy with P. distincta in two localities of São Paulo State (Brazil), where triploid animals also occur, as consequence of natural hybridisation. We analyse the chromosomes of P. distincta, P. tetraploidea, and their triploid hybrids, to enlighten the origin of polyploidy and to obtain some evidence on diploidisation of tetraploid karyotype.

Results

Phyllomedusa distincta was 2n = 2x = 26, whereas P. tetraploidea was 2n = 4x = 52, and the hybrid individuals was 2n = 3x = 39. In meiotic phases, bivalents were observed in the diploid males, whereas both bivalents and tetravalents were observed in the tetraploid males. Univalents, bivalents or trivalents; metaphase II cells carrying variable number of chromosomes; and spermatids were detected in the testis preparations of the triploid males, indicating that the triploids were not completely sterile. In natural and experimental conditions, the triploids cross with the parental species, producing abnormal egg clutches and tadpoles with malformations. The embryos and tadpoles exhibited intraindividual karyotype variability and all of the metaphases contained abnormal constitutions. Multiple NORs, detected by Ag-impregnation and FISH with an rDNA probe, were observed on chromosome 1 in the three karyotypic forms; and, additionally, on chromosome 9 in the diploids, mostly on chromosome 8 in the tetraploids, and on both chromosome 8 and 9 in the triploids. Nevertheless, NOR-bearing chromosome 9 was detected in the tetraploids, and chromosome 9 carried active or inactive NORs in the triploids. C-banding, base-specific fluorochrome stainings with CMA₃ and DAPI, FISH with a telomeric probe, and BrdU incorporation in DNA showed nearly equivalent patterns in the karyotypes of P. distincta, P. tetraploidea, and the triploid hybrids.

Conclusions

All the used cytogenetic techniques have provided strong evidence that the process of diploidisation, an essential step for stabilising the selective advantages produced by polyploidisation, is under way in distinct quartets of the tetraploid karyotype.