Cytogenetics of the Hybridogenetic Frog Pelophylax grafi and Its Parental Species Pelophylax perezi.

Hybrid taxa from the genus Pelophylax can propagate themselves in a modified way of sexual reproduction called hybridogenesis ensuring the formation of clonal gametes containing the genome of only one parental (host) species. Pelophylax grafi from South-Western Europe is a hybrid composed of P. ridibundus and P. perezi genomes and it lives with a host species P. perezi (P-G system). Yet it is unknown, whether non-Mendelian inheritance is fully maintained in such populations. In this study, we characterize P. perezi and P. grafi somatic karyotypes by using comparative genomic hybridization, genomic in situ hybridization, fluorescent in situ hybridization, and actinomycin D-DAPI. Here, we show the homeology of P. perezi and P. grafi somatic karyotypes to other Pelophylax taxa with 2n = 26 and equal contribution of ridibundus and perezi chromosomes in P. grafi which supports F1 hybrid genome constitution as well as a hemiclonal genome inheritance. We show that ridibundus chromosomes have larger regions of interstitial (TTAGGG)n repeats flanking the nucleolus organizing region on chromosome no. 10 and a high quantity of AT pairs in the centromeric regions. In P. perezi, we found species-specific sequences in metaphase chromosomes and marker structures in lampbrush chromosomes. Pericentromeric RrS1 repeat sequence was present in perezi and ridibundus chromosomes, but the blocks were stronger in ridibundus. Various cytogenetic techniques applied to the P-G system provide genome discrimination between ridibundus and perezi chromosomal sets. They could be used in studies of germ-line cells to explain patterns of clonal gametogenesis in P. grafi and broaden the knowledge about reproductive strategies in hybrid animals.

Genome elimination from the germline cells in diploid and triploid male water frogs Pelophylax esculentus.

Hybridogenesis is a hemiclonal reproductive strategy in diploid and triploid hybrids. Our study model is a frog P. esculentus (diploid RL and triploids RLL and RRL), a natural hybrid between P. lessonae (LL) and P. ridibundus (RR). Hybridogenesis relies on elimination of one genome (L or R) from gonocytes (G) in tadpole gonads during prespermatogenesis, but not from spermatogonial stem cells (SSCs) in adults. Here we provide the first comprehensive study of testis morphology combined with chromosome composition in the full spectrum of spermatogenic cells. Using genomic in situ hybridization (GISH) and FISH we determined genomes in metaphase plates and interphase nuclei in Gs and SSCs. We traced genomic composition of SSCs, spermatocytes and spermatozoa in individual adult males that were crossed with females of the parental species and gave progeny. Degenerating gonocytes (24%-39%) and SSCs (18%-20%) led to partial sterility of juvenile and adult gonads. We conclude that elimination and endoreplication not properly completed during prespermatogenesis may be halted when gonocytes become dormant in juveniles. After resumption of mitotic divisions by SSCs in adults, these 20% of cells with successful genome elimination and endoreplication continue spermatogenesis, while in about 80% spermatogenesis is deficient. Majority of abnormal cells are eliminated by cell death, however some of them give rise to aneuploid spermatocytes and spermatozoa which shows that hybridogenesis is a wasteful process.

Maintenance of pure hybridogenetic water frog populations: Genotypic variability in progeny of diploid and triploid parents.

An intriguing outcome of hybridisation is the emergence of clonally and hemiclonally reproducing hybrids, that can sustain, reproduce, and lead to the emergence of polyploid forms. However, the maintenance of diploid and polyploid hybrid complexes in natural populations remains unresolved. We selected water frogs from the Pelophylax esculentus complex to study how diploid and triploid hybrids, which reproduce hemiclonally via hybridogenesis, are maintained in natural populations. During gametogenesis in diploid hybrids, one of the parental genomes is eliminated, and the remaining genome is endoreplicated. In triploid hybrids, the single-copy genome is typically eliminated, while genome endoreplication does not occur. To investigate how diploid and triploid hybrid frogs reproduce in populations without parental species, we crossed these hybrid animals from two separate pure hybrid populations located in Poland. Using cytogenetic analysis of tadpoles that emerged from the crosses, we established which gametes were produced by parental hybrids. The majority of hybrid females and hybrid males produced one type of gamete with the P. ridibundus genome. However, in both studied populations, approximately half of the diploid and triploid hybrids simultaneously produced gametes with different genome compositions and ploidy levels, specifically, the P. ridibundus and P. lessonae genomes, as well as diploid gametes with genomes of both parental species. Triploid hybrid males and females mostly produced haploid gametes with the P. lessonae genome; however, gametes with the P. ridibundus genome have also been observed. These results suggest that not all hybrids follow the classical hybridogenetic reproduction program and reveal a significant level of alterations in the gametogenesis pathways. In addition, we found a variable survival rate of particular progeny genotypes when we crossed hybrid females with different males suggesting the important role of postzygotic barriers on the maintenance of pure hybrid systems. We suggest that the observed variability in produced gametes and the different survival rate of the progeny with certain genotypes is crucial for the existence of pure hybrid systems.

Hybridization and introgression between toads with different sex chromosome systems.

The growing interest in the lability of sex determination in non-model vertebrates such as amphibians and fishes has revealed high rates of sex chromosome turnovers among closely related species of the same clade. Can such lineages hybridize and admix with different sex-determining systems, or could the changes have precipitated their speciation? We addressed these questions in incipient species of toads (Bufonidae), where we identified a heterogametic transition and characterized their hybrid zone with genome-wide markers (RADseq). Adult and sibship data confirmed that the common toad B. bufo is female heterogametic (ZW), while its sister species the spined toad B. spinosus is male heterogametic (XY). Analysis of a fine scale transect across their parapatric ranges in southeastern France unveiled a narrow tension zone (∼10 km), with asymmetric mitochondrial and nuclear admixture over hundreds of kilometers southward and northward, respectively. The geographic extent of introgression is consistent with an expansion of B. spinosus across B. bufo's former ranges in Mediterranean France, as also suggested by species distribution models. However, widespread cyto-nuclear discordances (B. spinosus backrosses carrying B. bufo mtDNA) run against predictions from the dominance effects of Haldane's rule, perhaps because Y and W heterogametologs are not degenerated. Common and spined toads can thus successfully cross-breed despite fundamental differences in their sex determination mechanisms, but remain partially separated by reproductive barriers. Whether and how the interactions of their XY and ZW genes contribute to these barriers shall provide novel insights on the debated role of labile sex chromosomes in speciation.

Impacts of the synthetic androgen Trenbolone on gonad differentiation and development - comparisons between three deeply diverged anuran families.

Using a recently developed approach for testing endocrine disruptive chemicals (EDCs) in amphibians, comprising synchronized tadpole exposure plus genetic and histological sexing of metamorphs in a flow-through-system, we tested the effects of 17ß-Trenbolone (Tb), a widely used growth promoter in cattle farming, in three deeply diverged anuran families: the amphibian model species Xenopus laevis (Pipidae) and the non-models Bufo(tes) viridis (Bufonidae) and Hyla arborea (Hylidae). Trenbolone was applied in three environmentally and/or physiologically relevant concentrations (0.027 µg/L (10-10 M), 0.27 µg/L (10-9 M), 2.7 µg/L (10-8 M)). In none of the species, Tb caused sex reversals or masculinization of gonads but had negative species-specific impacts on gonad morphology and differentiation after the completion of metamorphosis, independently of genetic sex. In H. arborea and B. viridis, mounting Tb-concentration correlated positively with anatomical abnormalities at 27 µg/L (10-9 M) and 2.7 µg/L (10-8 M), occurring in X. laevis only at the highest Tb concentration. Despite anatomical aberrations, histologically all gonadal tissues differentiated seemingly normally when examined at the histological level but at various rates. Tb-concentration caused various species-specific mortalities (low in Xenopus, uncertain in Bufo). Our data suggest that deep phylogenetic divergence modifies EDC-vulnerability, as previously demonstrated for Bisphenol A (BPA) and Ethinylestradiol (EE2).

Shedding Light on a Secretive Tertiary urodelean Relict: Hynobiid salamanders (Paradactylodon persicus s.l.) from Iran, Illuminated by Phylogeographic, Developmental and Transcriptomic Data.

The Hyrcanian Forests present a unique Tertiary relict ecosystem, covering the northern Elburz and Talysh Ranges (Iran, Azerbaijan), a poorly investigated, unique biodiversity hotspot with many cryptic species. Since the 1970s, two nominal species of Urodela, Hynobiidae, Batrachuperus (later: Paradactylodon) have been described: Paradactylodon persicus from northwestern and P. gorganensis from northeastern Iran. Although P. gorganensis has been involved in studies on phylogeny and development, there is little data on the phylogeography, systematics, and development of the genus throughout the Hyrcanian Forests; genome-wide resources have been entirely missing. Given the huge genome size of hynobiids, making whole genome sequencing hardly affordable, we aimed to publish the first transcriptomic resources for Paradactylodon from an embryo and a larva (9.17 Gb RNA sequences; assembled to 78,918 unigenes). We also listed 32 genes involved in vertebrate sexual development and sex determination. Photographic documentation of the development from egg sacs across several embryonal and larval stages until metamorphosis enabled, for the first time, comparison of the ontogeny with that of other hynobiids and new histological and transcriptomic insights into early gonads and timing of their differentiation. Transcriptomes from central Elburz, next-generation sequencing (NGS) libraries of archival DNA of topotypic P. persicus, and GenBank-sequences of eastern P. gorganensis allowed phylogenetic analysis with three mitochondrial genomes, supplemented by PCR-amplified mtDNA-fragments from 17 museum specimens, documenting <2% uncorrected intraspecific genetic distance. Our data suggest that these rare salamanders belong to a single species P. persicus s.l. Humankind has a great responsibility to protect this species and the unique biodiversity of the Hyrcanian Forest ecosystems.

Common occurrence of Sharpey's fibres in amphibian phalanges.

Sharpey's fibres are known mainly as providing anchorage between tooth and the periodontal ligament but they occur also in other types of bones. In the postcranial skeleton these fibres are usually present at the muscle or tendon attachment sites. They were reported in all major groups of extant vertebrates, as well as in putative lissamphibian ancestors-temnospondyls and lepospondyls. However, it was recently stated that their presence was very rarely described in extant amphibians. In limbs, they were reported predominantly from proximal bones. They have not yet been reported from phalanges, which are the most commonly sectioned amphibian bones. Here, we describe phalangeal histology of nine species representing most major clades of lissamphibians. These results show that Sharpey's fibres occur commonly in lissamphibian phalanges. In shaft, they are radially oriented and occur in the periosteal bone, at sites of tendon attachment. They can also occur in the metaphysis and contact the cartilage. This may provide a basis for foot muscle reconstructions in fossil amphibians.

The programmed DNA elimination and formation of micronuclei in germ line cells of the natural hybridogenetic water frog Pelophylax esculentus.

DNA elimination is a radical form of gene silencing and occurs both in somatic and germ cells. The programmed DNA elimination occurs during gametogenesis in interspecies hybrids that reproduce by hybridogenesis (stick insects, fishes, and amphibians) and concerns removal of whole genomes of one of the parental species and production of clonal gametes propagating the genome of the other species. The cellular mechanisms differ considerably in hybridogenetic insects and fishes but remains unknown in edible frogs Pelophylax esculentus, natural hybrids between Pelophylax lessonae and Pelophylax ridibundus. Here we report DNA elimination mechanism in early developing gonads of diploid and triploid hybrid frogs, studied by TEM, immunofluorescence, and cytochemistry. In gonocytes of both sexes (primary oogonia and prespermatogonia), micronuclei emerge as detached nuclear buds formed during interphase. We found depletion of nuclear pore complexes in micronuclear membrane and chromatin inactivation via heterochromatinization followed by degradation of micronuclei by autophagy. Micronuclei formation does not lead to apoptotic cell death showing that genome elimination is a physiological process. Chromatin elimination via micronuclei in P. esculentus is unique among hybridogenetic animals and contributes to broadening the knowledge about reproductive modes in animals.

Prespermatogenesis and early spermatogenesis in frogs.

Spermatogenesis in frogs was for the first time divided into two phases: prespermatogenesis, when gonocytes proliferate in developing tadpole testes, and active spermatogenesis when spermatogonial stem cells (i.e. descendants of gonocytes), either self-renew or enter into meiotic cycles within cysts formed by Sertoli cells. We argue that amphibian larval gonocytes are homologues to mammalian gonocytes, whereas spermatogonial stem cells (SSCs) in adult frogs are homologous to mammalian single spermatogonia (As). Gonocytes constitute sex cords, i.e. the precursors of seminiferous tubules; they are bigger than SSCs and differ in morphology and ultrastructure. The nuclear envelope in gonocytes formed deep finger-like invaginations absent in SSCs. All stages of male germ cells contained lipid droplets, which were surrounded by glycogen in SSCs, but not in gonocytes. Mitochondria in gonocytes had enlarged edges of cristae, and in SSCs also lamellar mitochondria appeared. Minimal duration of prespermatogenesis was 46days after gonadal sex differentiation, but usually it lasted longer. SSCs give rise to secondary spermatogonia (equal to mammalian A, In, and B). Their lowest number inside a cyst was eight and this indicated the minimal number of cell cycles (three) of secondary spermatogonia necessary to enter meiosis. We sorted them according to the number of cell cycles (from 8 to 256 cells). This number is similar to that recorded for mammals as the result of a single As proliferation. The number of secondary spermatogonia correlates with the volume of a cyst. The general conclusion is that spermatogenesis in amphibians and mammals follows basically the same scheme.

Phalangeal bone anomalies in the European common toad Bufo bufo from polluted environments.

Every spring, many of amphibians are killed by motor vehicles on roads. These road-killed animals can be used as valuable material for non-invasive studies showing the effect of environmental pollution on amphibian populations. The aims of our research were to check whether the phalanges of road-killed toads may be useful as material for histological analysis, and whether various degrees of human impact influence the level in bone abnormalities in the common toad. We also examined whether the sex and age structure of toads can differ significantly depending in the different sites. We chose three toad breeding sites where road-killed individuals had been observed: near the centre of a city, the outskirts of a city, and a rural site. We collected dead individuals during spring migration in 2013. The sex of each individual was determined and the toes were used to determine age using the skeletochronology method. While performing age estimates, we looked for abnormalities in relation to normal bone tissue structure. In urban site, females dominate males (sex ratio 2.6:1), but in populations from rural and semi-urban sites, sex ratio was reverse (1:2.2 and 1:1.4, respectively). However, we did not find any significant differences between age structure of all populations (average age of each population: approximately 4 years). We observed abnormalities in more than 80 % of all toads from the city, compared to approximately 20 % from the rural and semi-urban sites. In particular, we found hypertrophic bone cells, misaligned intercellular substance, and irregular outer edges of bones. We suggest that these malformations are caused by different pollution, e.g. with heavy metals.

Impaired gonadal and somatic development corroborate vulnerability differences to the synthetic estrogen ethinylestradiol among deeply diverged anuran lineages.

Amphibians are undergoing a global decline. One poorly investigated reason could be the pollution of aquatic habitats by endocrine disrupting compounds (EDCs). We tested the susceptibility to the synthetically stabilized estrogen 17α-ethinylestradiol (EE2) in three deeply diverged anuran species, differing in sex determination systems, types of gonadogenesis and larval ecologies. To understand whether data from the amphibian model Xenopus laevis (Pipidae) are analogous and applicable to only distantly related non-model amphibians, tadpoles of X. laevis, Hyla arborea (Hylidae) and Bufo viridis (Bufonidae) were simultaneously exposed to 50, 500 and 5000ng/L EE2 from hatching until completion of metamorphosis, using a flow-through-system under identical experimental conditions. Comparing molecularly established genetic with histologically assessed phenotypic sex in all species, we have recently shown that EE2 provoked numerous genetic-male-to-phenotypic-female sex reversals and mixed sex individuals, confirming overall its expected feminizing effect. In the present study, we focus on the influence of EE2 on gonadal and somatic development. Anatomy and histology revealed several species-specific effects. In both non-model species, H. arborea and B. viridis, high numbers of anatomically impaired gonads were observed. In H. arborea, exposed to 5000ng/L EE2, numerous underdeveloped gonads were detected. Whereas EE2 did not alter snout-to-vent length and body weight of X. laevis metamorphs, H. arborea showed a treatment-dependent decrease, while B. viridis exhibited an increase in body weight and snout-to-vent length. Apart from a concentration-dependent occurrence of yellowish skin color in several H. arborea, no organ-specific effects were detected. Since EE2 ubiquitously occurs in many aquatic ecosystems and affects sexual and somatic development, among EDCs, it may indeed contribute to amphibian decline. The inter-species variation in developmental EE2-effects corroborates species-specific vulnerability differences towards EDCs between deeply diverged amphibian groups.

The plasticizer bisphenol A affects somatic and sexual development, but differently in pipid, hylid and bufonid anurans.

Due to their terrestrial habitats and aquatic reproduction, many amphibians are both very vulnerable and highly suitable bioindicators. The plasticizer bisphenol A (BPA) is one of the most produced chemical substances worldwide, and knowledge on its impacts on humans and animals is mounting. BPA is used for the industrial production of polycarbonate plastics and epoxy resins and found in a multitude of consumer products. Studies on BPA have involved mammals, fish and the fully aquatic anuran model Xenopus laevis. However, our knowledge about the sexual development of non-model, often semi-terrestrial anuran amphibians remains poor. Using a recently developed experimental design, we simultaneously applied BPA to two non-model species (Hyla arborea, Hylidae; Bufo viridis, Bufonidae) and the model X. laevis (Pipidae), compared their genetic and phenotypic sex for detection of sex reversals, and studied sexual development, focusing on anatomical and histological features of gonads. We compared three concentrations of BPA (0.023, 2.28 and 228 µg/L) to control groups in a high-standard flow-through-system, and tested whether conclusions, drawn from the model species, can be extrapolated to non-model anurans. In contrast to previous studies on fish and Xenopus, often involving dosages much higher than most environmental pollution data, we show that BPA causes neither the development of mixed sex nor of sex-reversed individuals (few, seemingly BPA-independent sex reversals) in all focal species. However, environmentally relevant concentrations, as low as 0.023 µg/L, were sufficient to provoke species-specific anatomically and histologically detectable impairments of gonads, and affected morphological traits of metamorphs. As the intensity of these effects differed between the three species, our data imply that BPA diversely affects amphibians with different evolutionary history, sex determination systems and larval ecologies. These results highlight the role of amphibians as a sensitive group that is responsive to environmental pollution.

Sex reversal assessments reveal different vulnerability to endocrine disruption between deeply diverged anuran lineages.

Multiple anthropogenic stressors cause worldwide amphibian declines. Among several poorly investigated causes is global pollution of aquatic ecosystems with endocrine disrupting compounds (EDCs). These substances interfere with the endocrine system and can affect the sexual development of vertebrates including amphibians. We test the susceptibility to an environmentally relevant contraceptive, the artificial estrogen 17α-ethinylestradiol (EE2), simultaneously in three deeply divergent systematic anuran families, a model-species, Xenopus laevis (Pipidae), and two non-models, Hyla arborea (Hylidae) and Bufo viridis (Bufonidae). Our new approach combines synchronized tadpole exposure to three EE2-concentrations (50, 500, 5,000 ng/L) in a flow-through-system and pioneers genetic and histological sexing of metamorphs in non-model anurans for EDC-studies. This novel methodology reveals striking quantitative differences in genetic-male-to-phenotypic-female sex reversal in non-model vs. model species. Our findings qualify molecular sexing in EDC-analyses as requirement to identify sex reversals and state-of-the-art approaches as mandatory to detect species-specific vulnerabilities to EDCs in amphibians.

Sex-chromosome differentiation parallels postglacial range expansion in European tree frogs (Hyla arborea).

Occasional XY recombination is a proposed explanation for the sex-chromosome homomorphy in European tree frogs. Numerous laboratory crosses, however, failed to detect any event of male recombination, and a detailed survey of NW-European Hyla arborea populations identified male-specific alleles at sex-linked loci, pointing to the absence of XY recombination in their recent history. Here, we address this paradox in a phylogeographic framework by genotyping sex-linked microsatellite markers in populations and sibships from the entire species range. Contrasting with postglacial populations of NW Europe, which display complete absence of XY recombination and strong sex-chromosome differentiation, refugial populations of the southern Balkans and Adriatic coast show limited XY recombination and large overlaps in allele frequencies. Geographically and historically intermediate populations of the Pannonian Basin show intermediate patterns of XY differentiation. Even in populations where X and Y occasionally recombine, the genetic diversity of Y haplotypes is reduced below the levels expected from the fourfold drop in copy numbers. This study is the first in which X and Y haplotypes could be phased over the distribution range in a species with homomorphic sex chromosomes; it shows that XY-recombination patterns may differ strikingly between conspecific populations, and that recombination arrest may evolve rapidly (<5000 generations).

Sex chromosome conservation, DMRT1 phylogeny and gonad morphology in diploid Palearctic green toads (Bufo viridis subgroup).

Due to the prevailing sex chromosome homomorphy and large genome size, the knowledge on sex determination systems, sex chromosomes and sex-determining genes in amphibians remains scarce. Using 3 cross-amplifying sex-linked microsatellite markers, we uncover sex determination systems and sex chromosomes in purebred, diploid Palearctic green toads (Bufo viridis subgroup), which had so far only been characterized in laboratory-bred hybrids. Our data support an XY system in B. balearicus, B. viridis and B. variabilis. While females show recombination, it is strongly suppressed (or not detectable) in males. Markers corroborate the largest chromosome pair 1 (homologous to linkage group 1 of Xenopus tropicalis) to represent the sex chromosomes in diploid species of the B. viridis subgroup (B. siculus, B. shaartusiensis, B. balearicus, B. turanensis, B. variabilis, B. viridis, and probably B. boulengeri). This chromosome harbors DMRT1, a key gene of the sexual pathway in deeply divergent animal taxa. However, our phylogenetic analysis of a 600-bp fragment of that gene in diploid green toad taxa reveals that X and Y alleles cluster by species and not by gametolog. This suggests that XY-sequence similarity stems from occasional XY recombination within DMRT1, and we preliminarily reject its role as the master sex determination gene, pending future extension of this evidence to the entire DMRT1 gene. We further create a chain of evidence, which supports the hypothesis that linkage group 1 of X. tropicalis appears to be maintained as the largest chromosome (1), and thus is homologous in anuran karyotype evolution from pipid to hylid, bufonid and ranid anurans.