Homomorphic sex chromosomes and rapid turnover of sex-determining genes can complicate establishing the sex chromosome system operating in a given species. This difficulty exists in Xenopus tropicalis, an anuran quickly becoming a relevant model for genetic, genomic, biochemical, and ecotoxicological research. Despite the recent interest attracted by this species, little is known about its sex chromosome system. Direct evidence that females are the heterogametic sex, as in the related species Xenopus laevis, has yet to be presented. Furthermore, X. laevis' sex-determining gene, DM-W, does not exist in X. tropicalis, and the sex chromosomes in the two species are not homologous. Here we identify X. tropicalis' sex chromosome system by integrating data from (i) breeding sex-reversed individuals, (ii) gynogenesis, (iii) triploids, and (iv) crosses among several strains. Our results indicate that at least three different types of sex chromosomes exist: Y, W, and Z, observed in YZ, YW, and ZZ males and in ZW and WW females. Because some combinations of parental sex chromosomes produce unisex offspring and other distorted sex ratios, understanding the sex-determination systems in X. tropicalis is critical for developing this flexible animal model for genetics and ecotoxicology.

The study of sex determination and sex chromosome organisation in non-model species has long been technically challenging, but new sequencing methodologies are now enabling precise and high-throughput identification of sex-specific genomic sequences. In particular, Restriction Site-Associated DNA Sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software designed to specifically search for sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analyzed a published dataset of Japanese medaka, Oryzias latipes , where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyze new RAD-Seq datasets from 15 fish species spanning multiple systematic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in non-model species and outperforms the commonly used RAD-Seq analysis software STACKS.RADSex in speed, resource usage, ease of application, and visualization options. Furthermore, our analysis of new datasets from 15 species provides new insights on sex determination in fish.

The study of sex determination and sex chromosome organisation in non-model species has long been technically challenging, but new sequencing methodologies are now enabling precise and high-throughput identification of sex-specific genomic sequences. In particular, Restriction Site-Associated DNA Sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software designed to specifically search for sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analyzed a published dataset of Japanese medaka, Oryzias latipes , where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyze new RAD-Seq datasets from 15 fish species spanning multiple systematic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in non-model species and outperforms the commonly used RAD-Seq analysis software STACKS.RADSex in speed, resource usage, ease of application, and visualization options. Furthermore, our analysis of new datasets from 15 species provides new insights on sex determination in fish.

The role of environmental factors in sexual differentiation in amphibians is not new. The effect of hormones or hormone-like compounds is widely demonstrated. However, the effect of temperature has traditionally been regarded as something anecdotal that occurs in extreme situations and not as a factor to be considered. The data currently available reveal a different situation. Sexual differentiation in some amphibian species can be altered even by small changes in temperature. On the other hand, although not proven, it is possible that temperature is related to the appearance of sex-reversed individuals in natural populations under conditions unrelated to environmental contaminants. According to this, temperature, through sex reversal (phenotypic sex opposed to genetic sex), could play an important role in the turnover of sex-determining genes and in the maintenance of homomorphic sex chromosomes in this group. Accordingly, and given the expected increase in global temperatures, growth and sexual differentiation in amphibians could easily be affected, altering the sex ratio in natural populations and posing major conservation challenges for a group in worldwide decline. It is therefore particularly urgent to understand the mechanism by which temperature affects sexual differentiation in amphibians.

Sex is determined genetically in amphibians; however, little is known about the sex chromosomes, testis-determining genes, and the genes involved in testis differentiation in this class. Certain inherent characteristics of the species of this group, like the homomorphic sex chromosomes, the high diversity of the sex-determining mechanisms, or the existence of polyploids, may hinder the design of experiments when studying how the gonads can differentiate. Even so, other features, like their external development or the possibility of inducing sex reversal by external treatments, can be helpful. This review summarizes the current knowledge on amphibian sex determination, gonadal development, and testis differentiation. The analysis of this information, compared with the information available for other vertebrate groups, allows us to identify the evolutionarily conserved and divergent pathways involved in testis differentiation. Overall, the data confirm the previous observations in other vertebrates—the morphology of the adult testis is similar across different groups; however, the male-determining signal and the genetic networks involved in testis differentiation are not evolutionarily conserved.

Amphibians have some of the most variable genome sizes among vertebrates. Genome size variation has been attributed to repetitive and noncoding DNA, including satellite repeats, transposable elements, introns, and nuclear insertions of viral and organelle DNA. In vertebrates, satellite DNAs have been widely described in mammals, but few molecular studies have been carried out in amphibians. Here, we provide a detailed characterization of a new family of satellite DNA, present in all 15 examined species of the family Bufonidae. Southern-blot analysis and PCR reveal that this satellite is formed by monomers of 807 bp, is organized in tandem arrays, and has an AT-content of 57.4%. Phylogenetic analyses show that most clades exhibit species-specific variances, indicating that this satellite DNA has evolved by concerted evolution. The homogenization/fixation process is heterogeneous in Bufonidae, where the genera Bufo and Bufotes do not show species-specific differences, while populations from Rhinella marina exhibit population-specific changes. Additionally, variants of this satellite DNA have been identified in Duttaphrynus melanostictus and R. marina, supporting the ‘library hypothesis’ (a set, ‘library’, of satellite DNAs is shared by a species group). Physical mapping in Bufo bufo, Bufo spinosus, Epidalea calamita and Bufotes viridis provides evidence that this repetitive DNA is not dispersed in the karyotype, but accumulated in pericentromeric regions of some chromosomal pairs. This location, together with its presence in the transcriptomes of bufonids, could indicate a role in centromere function or heterochromatin formation and maintenance.

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.

Taxonomy in Bufonidae witnessed notable transformations. Bufotes viridis and Epidalea calamita, previously included in genus Bufo, were relocated in other genera, while the genus Bufo was restricted to members of the earlier Bufo bufo group. On the other hand, Bufo bufo sensu lato now includes four species: Bufo bufo, Bufo spinosus, Bufo verricossissimus and Bufo eichwaldi. In this study, we examined three species of three Bufonidae genera (B. spinosus, B. viridis and E. calamita) by conventional (C-banding and Ag-NOR staining) and molecular (in situ hybridization with probes for telomeric repeats and rDNA loci, and genomic in situ hybridization (GISH)) cytogenetic methods. C-banding patterns are reported for the first time for B. spinosus and E. calamita populations from Iberian Peninsula and for B. viridis from Greece, and reveal several differences with the reported C-banded karyotypes described for other European populations of these species. Silver staining shows size heteromorphisms of the signals at the Nucleolar Organizing Region (NOR). By contrast, FISH with ribosomal probes only reveal size heteromorphism of rDNA sequences in E. calamita, suggesting that the differences observed after silver staining in B. spinosus and B. viridis should be attributed to differences in chromosomal condensation and/or gene activity rather than to differences in the copy number for ribosomal genes. Regarding telomeric repeats, E. calamita is the only species with interstitial telomeric sequences (ITS) located on centromeric regions, probably originated by accumulation of telomeric sequences in the centromeric heterochromatin. Finally, we analyzed the composition and distribution of repetitive sequences by genome in situ hybridization. These experiments reveal the accumulation of repetitive sequences in centromeric regions of the three species, although these sequences are not conserved when species from different genera are compared.