Gonadal sex differentiation in frogs: how testes become shorter than ovaries

Gonadal morphogenesis in the South American toad Rhinella Arenarum (Anura, Bufonidae) unveils an extremely delayed rate of sex differentiation

Among anurans, Bufonids are recognized for their retarded sex differentiation. However, few studies have addressed gonadal morphogenesis in this family. Here, we analyzed the early gonadogenesis in laboratory-reared Rhinella arenarum. Few germ cells were identified in the genital ridge at Gosner stage 26. At metamorphosis, somatic cells and germ cells were observed in the outer region of the undifferentiated gonad, whereas the central region was occupied by stromal tissue. A cortico-medullary organization was first recognized on Day 7 postmetamorphosis. The cortex was composed of germ cells and encompassing epithelial cells, whereas the medulla contained cells presumptively derived from the coelomic epithelium. Medullary somatic cells formed metameric knots along the length of the undifferentiated gonad. Consequently, a series of 12-14 gonomeres became recognizable externally. The first sign of ovarian differentiation was observed on Day 15 postmetamorphosis, when a cavity was formed within each gonomere. In contrast, testes were recognized by a uniform distribution of germ cells and intermingled somatic cells, as the division into cortex and medulla was lost. By Day 50 postmetamorphosis, the gonadal metameric organization was still apparent both in the ovaries and testes. Follicles containing diplotene oocytes were observed within the ovary. In the testis, an incipient lobular architecture was recognized without initiation of meiosis within the seminiferous cords. These observations reveal an extremely delayed gonadal development in R. arenarum, not reported previously for other anuran species. In addition, the late differentiation of the gonads contrasted with the early appearance of follicles in the Bidder's organ. Lastly, we observed that delayed metamorphs exhibited an undifferentiated gonad, demonstrating that gonadogenesis in this species is more dependent on somatic development than on age.

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.

Bisphenol B disrupts testis differentiation partly via the estrogen receptor-mediated pathway and subsequently causes testicular dysgenesis in Xenopus laevis

There is growing concern about adverse effects of bisphenol A alternatives including bisphenol B (BPB) due to their estrogenic activity. However, limited data are available concerning the influences of BPB on male reproductive development in vertebrates, especially in amphibians, which are believed to be susceptible to estrogenic chemicals. The present study investigated the effects of 10, 100 and 1000 nM BPB (2.42, 24.2 and 242 μg/L) on testis development in Xenopus laevis, a model amphibian species for studying gonadal feminization. We found that exposure to BPB from stages 45/46 to 52 resulted in down-regulation of testis-biased gene expression and up-regulation of ovary-biased gene and vitellogenin (vtgb1) expression in gonad-mesonephros complexes (GMCs) of tadpoles at stage 52, coupled with suppressed cell proliferation in testes and reduced gonadal metameres, resembling the effects of 17ß-estradiol. Moreover, an estrogen receptor (ER) antagonist ICI 182780 antagonized BPB-caused up-regulation of ovary-biased gene and vtgb1 expression to some degree, indicating that the effects of BPB on X. laevis testis differentiation could be partly mediated by ER. All observations demonstrate that early exposure to BPB inhibited testis differentiation and exerted certain feminizing effects during gonadal differentiation. When exposure was extended to post-metamorphosis, testes exhibited histological and morphological abnormalities including segmented, discontinuous and fragmented shapes, besides altered sex-dimorphic gene expression. Notably, most of BPB-caused alterations were not concentration-dependent, but the lowest concentration indeed exerted significant effects. Overall, our study for the first time reveals that low concentrations of BPB can disrupt testis differentiation partly due to its estrogenic activity and subsequently cause testicular dysgenesis after metamorphosis, highlighting its reproductive risk to amphibians and other vertebrates including humans. Our finding also implies that estrogenic chemicals-caused testis differentiation inhibition at tadpole stages could predict later testicular dysgenesis after metamorphosis, meaning a possibility of early detection of abnormal testis development caused by estrogenic chemicals.

Analyzing the gonadal transcriptome of the frog Hoplobatrachus rugulosus to identify genes involved in sex development

BACKGROUND

The tiger frog (Hoplobatrachus rugulosus) is listed as a national Class II protected species in China. In the context of global warming, the sex ratio of amphibians will be affected, and the development of the population will be limited. Therefore, considering the potential for a decrease in the number of amphibians, studying sex evolution and molecular regulation of gonadal development in H. rugulosus, phenomenon that are currently unclear, is of great significance.

RESULTS

Here, H. rugulosus was used to explore the mechanisms regulating gonadal development in amphibians. Illumina HiSeq 3000 was used to sequence the gonadal transcriptome of male and female H. rugulosus at two growth stages to identify genes related to gonadal development and analyze expression differences in the gonads. This analysis indicated that cyp17α, hsd3β, hsd11β1, cyp19α, and hsd17β12 perform vital functions in sex development in amphibians. Specifically, the expression of cyp3α, cyp17α, hsd3β, hsd11β1, sox2, sox9, sox30, soat, cyp19α, hsd17β12, and hspα1s was correlated with gonadal development and differentiation in H. rugulosus, as determined using the quantitative reverse transcriptase-polymerase chain reaction.

CONCLUSION

Significant differences were found in the gonadal gene expression levels in H. rugulosus of both sexes, and we identified a steroid hormone synthesis pathway in this species and analyzed related gene expression, but the changes during sex differentiation were still unclear. To our knowledge, this report presents the first analysis of the H. rugulosus gonadal transcriptome and lays the foundation for future research.

Testis Development and 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.

Development of Testis Cords and the Formation of Efferent Ducts in Xenopus laevis: Differences and Similarities with Other Vertebrates

The knowledge of testis development in amphibians relative to amniotes remains limited. Here, we used Xenopus laevis to investigate the process of testis cord development. Morphological observations revealed the presence of segmental gonomeres consisting of medullary knots in male gonads at stages 52-53, with no distinct gonomeres in female gonads. Further observations showed that cell proliferation occurs at specific sites along the anterior-posterior axis of the future testis at stage 50, which contributes to the formation of medullary knots. At stage 53, adjacent gonomeres become close to each other, resulting in fusion; then (pre-)Sertoli cells aggregate and form primitive testis cords, which ultimately become testis cords when germ cells are present inside. The process of testis cord formation in X. laevis appears to be more complex than in amniotes. Strikingly, steroidogenic cells appear earlier than (pre-)Sertoli cells in differentiating testes of X. laevis, which differs from earlier differentiation of (pre-)Sertoli cells in amniotes. Importantly, we found that the mesonephros is connected to the testis gonomere at a specific site at early larval stages and that these connections become efferent ducts after metamorphosis, which challenges the previous concept that the mesonephric side and the gonadal side initially develop in isolation and then connect to each other in amphibians and amniotes.

When Sex Chromosomes Recombine Only in the Heterogametic Sex: Heterochiasmy and Heterogamety in Hyla Tree Frogs

Sex chromosomes are classically predicted to stop recombining in the heterogametic sex, thereby enforcing linkage between sex-determining (SD) and sex-antagonistic (SA) genes. With the same rationale, a pre-existing sex asymmetry in recombination is expected to affect the evolution of heterogamety, for example, a low rate of male recombination might favor transitions to XY systems, by generating immediate linkage between SD and SA genes. Furthermore, the accumulation of deleterious mutations on nonrecombining Y chromosomes should favor XY-to-XY transitions (which discard the decayed Y), but disfavor XY-to-ZW transitions (which fix the decayed Y as an autosome). Like many anuran amphibians, Hyla tree frogs have been shown to display drastic heterochiasmy (males only recombine at chromosome tips) and are typically XY, which seems to fit the above expectations. Instead, here we demonstrate that two species, H. sarda and H. savignyi, share a common ZW system since at least 11 Ma. Surprisingly, the typical pattern of restricted male recombination has been maintained since then, despite female heterogamety. Hence, sex chromosomes recombine freely in ZW females, not in ZZ males. This suggests that heterochiasmy does not constrain heterogamety (and vice versa), and that the role of SA genes in the evolution of sex chromosomes might have been overemphasized.

Tissues injury and pathological changes in Hyla intermedia juveniles after chronic larval exposure to tebuconazole

Tebuconazole (TBZ), an azole pesticide, is one of the most frequently detected fungicides in surface water. Despite its harmful effects, mainly related to endocrine disturbance, the consequences of TBZ exposure in amphibians remain poorly understood. Here, we investigated the adverse and delayed effects of TBZ chronic exposure on a native anuran species, often inhabiting cultivated areas, the Italian tree frog (Hyla intermedia). To disclose the multiple mechanisms of action through which TBZ exerts its toxicity we exposed tadpoles over the whole larval period to two sublethal TBZ concentrations (5 and 50 μg/L), and we evaluated histological alterations in three target organs highly susceptible to xenobiotics: liver, kidney, and gonads. We also assessed morphometric and gravimetric parameters: snout-vent length (SVL), body mass (BM), liver somatic index (LSI), and gonad-mesonephros complex index (GMCI) and determined sex ratio, gonadal development, and differentiation. Our results show that TBZ induces irreversible effects on multiple target organs in H. intermedia, exerting its harmful effects through several pathological pathways, including a massive inflammatory response. Moreover, TBZ markedly affects sexual differentiation also by inducing the appearance of sexually undetermined individuals and a general delay of germ cell maturation. Given the paucity of data on the effects of TBZ in amphibians, our results will contribute to a better understanding of the environmental risk posed by this fungicide to the most endangered group of vertebrates.

Pattern of Gonadal Sex Differentiation in the Rice Field Frog Hoplobatrachus rugulosus (Anura: Dicroglossidae)

Sex differentiation during gonadal development is diversified among anuran amphibian species. In this study, the anuran experimental species Hoplobatrachus rugulosus was examined. The pattern of gonadal sex differentiation was observed by morphological and histological approaches. The gonad was observed morphologically at Gosner stage 33, while distinct testis and ovary were evident from 3-4 weeks after metamorphosis ended. Histological analysis showed that genital ridge formation began at stage 25 and ovarian differentiation began at stage 36. The developing ovary appeared with numerous primary oogonia, which developed into oocytes, while the medulla regressed to form an ovarian cavity. During metamorphosis, only an ovary was observed. Testicular differentiation seemed to begin later, during the first week after metamorphosis, and occurred via an intersex condition. The intersex gonads contained developing testicular tissue with both normal and atretic oocytes. The fully developed testis was first identified at 6 weeks after metamorphosis. Comparing the times of gonadal differentiation and somatic development revealed that the ovary exhibited a basic rate of differentiation while the testis exhibited a retarded one. These results establish that males of this species develop later than do females, and the testis develops through an intersex gonad, as is evident from its seminiferous cord formation, the presence of testis-ova, and atretic oocytes in the tissue. Thus, the pattern of gonadal sex differentiation in H. rugulosus is an undifferentiated type, in which only female gonads are observed during metamorphosis and intersex and male gonads are observed later. These results are crucial for further research on the sexual development of anurans.

Novel genetic sex markers reveal high frequency of sex reversal in wild populations of the agile frog (Rana dalmatina) associated with anthropogenic land use

Populations of ectothermic vertebrates are vulnerable to environmental pollution and climate change because certain chemicals and extreme temperatures can cause sex reversal during early ontogeny (i.e. genetically female individuals develop male phenotype or vice versa), which may distort population sex ratios. However, we have troublingly little information on sex reversals in natural populations, due to unavailability of genetic sex markers. Here, we developed a genetic sexing method based on sex-linked single nucleotide polymorphism loci to study the prevalence and fitness consequences of sex reversal in agile frogs (Rana dalmatina). Out of 125 juveniles raised in laboratory without exposure to sex-reversing stimuli, 6 showed male phenotype but female genotype according to our markers. These individuals exhibited several signs of poor physiological condition, suggesting stress-induced sex reversal and inferior fitness prospects. Among 162 adults from 11 wild populations in North-Central Hungary, 20% of phenotypic males had female genotype according to our markers. These individuals occurred more frequently in areas of anthropogenic land use; this association was attributable to agriculture and less strongly to urban land use. Female-to-male sex-reversed adults had similar body mass as normal males. We recorded no events of male-to-female sex reversal either in the laboratory or in the wild. These results support recent suspicions that sex reversal is widespread in nature, and suggest that human-induced environmental changes may contribute to its pervasiveness. Furthermore, our findings indicate that sex reversal is associated with stress and poor health in early life, but sex-reversed individuals surviving to adulthood may participate in breeding.

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.

Evaluation of the effects of low concentrations of bisphenol AF on gonadal development using the Xenopus laevis model: A finding of testicular differentiation inhibition coupled with feminization

Developmental exposures to estrogenic chemicals possibly cause structural and functional abnormalities of reproductive organs in vertebrates. Bisphenol AF (BPAF), a bisphenol A (BPA) analogue, has been shown to have higher estrogenic activity than BPA, but little is known about the effects of BPAF on gonadal development, particularly gonadal differentiation. We aimed to determine whether low concentrations of BPAF could disrupt gonadal differentiation and subsequent development using Xenopus laevis, a model species for studying feminizing effects of estrogenic chemicals. X. laevis tadpoles were exposed to BPAF (1, 10, 100 nM) or 17β-estradiol (E2, positive control) from stages 45/46 to 53 and 66 in a semi-static exposure system, with a prolonged treatment with the highest concentration to the eighth week post-metamorphosis (WPM8). Gonadal morphology and histology as well as sexually dimorphic gene expression were examined to evaluate the effects of BPAF. All concentrations of BPAF caused changes in testicular morphology at different developmental stages compared with controls. Specifically, at stage 53, BPAF like E2 resulted in decreases in both the size and the number of gonadal metameres (gonomeres) in testes, looking like ovaries. Some of BPAF-treated testes remained segmented and even became discontinuous and fragmented at subsequent stages. Histological abnormalities were also observed in BPAF-treated testes, such as ovarian cavity at stages 53 and 66 and poorly developed seminiferous tubules on WPM8. At the molecular level, BPAF inhibited expression of male highly expressed genes in testes at stage 53. Correspondingly, BPAF, like E2, inhibited cell proliferation in testes at stage 50. All results show that low concentrations of BPAF inhibited testicular differentiation and subsequent development in X. laevis, along with feminizing effects to some degree. Our finding implies a risk of BPAF to the male reproductive system of vertebrates including humans.

Gonadal, body color, and genotoxic alterations in Lithobates catesbeianus tadpoles exposed to nonylphenol

Endocrine disrupting chemicals are one of the most important factors contributing to worldwide amphibian decline. The 4-nonylphenol (NP) is a degradation product of several compounds, such as detergents and pesticides, affecting the aquatic environment. Here, we test whether treatment with NP has an effect on developing ovarian tissue, nuclear abnormalities in erythrocytes, and body darkness in pre-metamorphic tadpoles of the bullfrog Lithobates catesbeianus. Tadpoles were exposed for 14 days to three different concentrations of NP (1, 10, and 100 μg/L) besides the control group, which was maintained only with water. After determining body coloration, animals were euthanized and gonads and blood were collected and processed for histology and genotoxic analysis. Even though most animals were females, intersex tadpoles were observed in control and treated groups and there were no males in any group. The highest concentration of NP showed an increase in atretic oocytes, but the area corresponding to somatic compartment and early and late germ cells were not affected. Furthermore, all treated groups presented higher amount of nuclear abnormalities in erythrocytes and body darkening when compared with the control group. These results suggest that NP causes genetic damage and morphological alterations in L. catesbeianus tadpoles by disrupting oogenesis, inducing genotoxicity and increasing body coloration. Its effects on gonadal development could cause future impairments in reproduction, while its deleterious effects on genotoxicity and body pigmentation could be used as a biomarker of effect to this compound.

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.

Evolutionary and developmental dynamics of sex-biased gene expression in common frogs with proto-Y chromosomes

BACKGROUND

The patterns of gene expression on highly differentiated sex chromosomes differ drastically from those on autosomes, due to sex-specific patterns of selection and inheritance. As a result, X chromosomes are often enriched in female-biased genes (feminization) and Z chromosomes in male-biased genes (masculinization). However, it is not known how quickly sexualization of gene expression and transcriptional degeneration evolve after sex-chromosome formation. Furthermore, little is known about how sex-biased gene expression varies throughout development.

RESULTS

We sample a population of common frogs (Rana temporaria) with limited sex-chromosome differentiation (proto-sex chromosome), leaky genetic sex determination evidenced by the occurrence of XX males, and delayed gonadal development, meaning that XY individuals may first develop ovaries before switching to testes. Using high-throughput RNA sequencing, we investigate the dynamics of gene expression throughout development, spanning from early embryo to froglet stages. Our results show that sex-biased expression affects different genes at different developmental stages and increases during development, reaching highest levels in XX female froglets. Additionally, sex-biased gene expression depends on phenotypic, rather than genotypic sex, with similar expression in XX and XY males; correlates with gene evolutionary rates; and is not localized to the proto-sex chromosome nor near the candidate sex-determining gene Dmrt1.

CONCLUSIONS

The proto-sex chromosome of common frogs does not show evidence of sexualization of gene expression, nor evidence for a faster rate of evolution. This challenges the notion that sexually antagonistic genes play a central role in the initial stages of sex-chromosome evolution.

A rapid rate of sex-chromosome turnover and non-random transitions in true frogs

The canonical model of sex-chromosome evolution predicts that, as recombination is suppressed along sex chromosomes, gametologs will progressively differentiate, eventually becoming heteromorphic. However, there are numerous examples of homomorphic sex chromosomes across the tree of life. This homomorphy has been suggested to result from frequent sex-chromosome turnovers, yet we know little about which forces drive them. Here, we describe an extremely fast rate of turnover among 28 species of Ranidae. Transitions are not random, but converge on several chromosomes, potentially due to genes they harbour. Transitions also preserve the ancestral pattern of male heterogamety, in line with the ‘hot-potato’ model of sex-chromosome transitions, suggesting a key role for mutation-load accumulation in non-recombining genomic regions. The importance of mutation-load selection in frogs might result from the extreme heterochiasmy they exhibit, making frog sex chromosomes differentiate immediately from emergence and across their entire length.

The evolutionary forces that favour transitions in sex chromosomes are not well understood. Here, Jeffries and colleagues show a very high rate of sex chromosome turnover in true frogs, which may be driven by rapid mutation-load accumulation due to the low recombination rate in males.

Tissue Specificity and Dynamics of Sex-Biased Gene Expression in a Common Frog Population with Differentiated, Yet Homomorphic, Sex Chromosomes

Sex-biased genes are central to the study of sexual selection, sexual antagonism, and sex chromosome evolution. We describe a comprehensive de novo assembled transcriptome in the common frog Rana temporaria based on five developmental stages and three adult tissues from both sexes, obtained from a population with karyotypically homomorphic but genetically differentiated sex chromosomes. This allows the study of sex-biased gene expression throughout development, and its effect on the rate of gene evolution while accounting for pleiotropic expression, which is known to negatively correlate with the evolutionary rate. Overall, sex-biased genes had little overlap among developmental stages and adult tissues. Late developmental stages and gonad tissues had the highest numbers of stage- or tissue-specific genes. We find that pleiotropic gene expression is a better predictor than sex bias for the evolutionary rate of genes, though it often interacts with sex bias. Although genetically differentiated, the sex chromosomes were not enriched in sex-biased genes, possibly due to a very recent arrest of XY recombination. These results extend our understanding of the developmental dynamics, tissue specificity, and genomic localization of sex-biased genes.

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.

Mutual maintenance of di- and triploid Pelophylax esculentus hybrids in R-E systems: results from artificial crossings experiments

Background

Interspecies animal hybrids can employ clonal or hemiclonal reproduction modes where one or all parental genomes are transmitted to the progeny without recombination. Nevertheless, some interspecies hybrids retain strong connection with the parental species needed for successful reproduction. Appearance of polyploid hybrid animals may play an important role in the substitution of parental species and in the speciation process.

Results

To establish the mechanisms that enable parental species, diploid and polyploid hybrids coexist we have performed artificial crossing experiments of water frogs of Pelophylax esculentus complex. We identified tadpole karyotypes and oocyte genome composition in all females involved in the crossings. The majority of diploid and triploid hybrid frogs produced oocytes with 13 bivalents leading to haploid gametes with the same genome as parental species hybrids usually coexist with. After fertilization of such gametes only diploid animals appeared. Oocytes with 26 bivalents produced by some diploid hybrid frogs lead to diploid gametes, which give rise to triploid hybrids after fertilization. In gonads of all diploid and triploid hybrid tadpoles we found DAPI-positive micronuclei (nucleus-like bodies) involved in selective genome elimination. Hybrid male and female individuals produced tadpoles with variable karyotype and ploidy even in one crossing owing to gametes with various genome composition.

Conclusions

We propose a model of diploid and triploid hybrid frog reproduction in R-E population systems. Triploid Pelophylax esculentus hybrids can transmit genome of parental species they coexist with by producing haploid gametes with the same genome composition. Triploid hybrids cannot produce triploid individuals after crossings with each other and depend on diploid hybrid females producing diploid eggs. In contrast to other population systems, the majority of diploid and triploid hybrid females unexpectedly produced gametes with the same genome as parental species hybrids coexist with.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1063-3) contains supplementary material, which is available to authorized users.

Chronic exposures to fungicide pyrimethanil: multi-organ effects on Italian tree frog (Hyla intermedia)

Amphibian habitats are easily contaminated by several pollutants, and in agricultural landscapes the likely exposure scenario is represented by pesticides. Many of these substances are known or suspected to act as endocrine disrupting chemicals (EDCs). The goal of the present study was to assess the effects of pyrimethanil, a common-used but also overlooked fungicide, on liver, kidney and gonadal differentiation of Hyla intermedia. Through a multi-organ evaluation, we demonstrated that a long term exposure to two environmentally relevant concentrations of pyrimethanil (5 and 50 µg/L) elicits a range of toxic responses. First we showed that pyrimethanil induces underdevelopment of ovaries and interferes with normal sexual differentiation, thus revealing the endocrine disruption potential of this fungicide. Moreover we revealed that all considered organs are seriously affected by this fungicide and both necrosis and apoptosis contribute to the histological response. This is the first report on the effects of pyrimethanil on gonads, liver and kidney histology of a non-model species and it demonstrates that the hazardous properties of this fungicide can result from several pathological processes affecting different key compartments of amphibian.

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 (A_s). 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 A_s 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.

Lethal and sublethal effects of a methomyl-based insecticide in Hoplobatrachus rugulosus

The aim of this study was to investigate the lethal and sublethal toxicity of a methomyl-based insecticide in Hoplobatrachus rugulosus, as methomyl-based insecticides are applied in massive amounts and agrochemicals have effects on the decline in amphibian populations. To evaluate the toxic effects of methomyl from agricultural application, a methomyl-based insecticide containing 40% methomyl was selected. The median lethal concentration of 96 hours of methomyl exposure was 8.69 ppm for H. rugulosus tadpoles. The lethal concentration also produced severe histological damage in the livers and kidneys of the exposed tadpoles. The sublethal concentration used for methomyl was 144 ppb during the metamorphosis period. It was found that the sublethal concentration of the methomyl compound could decrease growth, metamorphosis time, and size, disturb biochemical parameters, and produce histological damage. In livers, methomyl effects increased oxidative stress and dramatically decreased the glycogen level of the treated froglets. Mononuclear infiltration, blood congestion, amorphous substances, and hepatocytes vacuolization were observed throughout liver tissue. The methomyl-based insecticide also increased oxidative stress and decreased nitric oxide levels in the kidneys of the exposed froglets. Renal tissue damage including blood congestion, amorphous substances, and Bowman’s capsule spaces reduction were found in the methomyl exposure group. The methomyl compound also produced vacuoles in various stages of oocytes, but no histological damage was found in testicular tissue. Our results indicated strong toxic effects of the methomyl-based insecticide on H. rugulosus, a broadly tolerant anuran.

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.

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.

Dmrt1 polymorphism covaries with sex-determination patterns in Rana temporaria

Patterns of sex-chromosome differentiation and gonadal development have been shown to vary among populations of Rana temporaria along a latitudinal transect in Sweden. Frogs from the northern-boreal population of Ammarnäs displayed well-differentiated X and Y haplotypes, early gonadal differentiation, and a perfect match between phenotypic and genotypic sex. In contrast, no differentiated Y haplotypes could be detected in the southern population of Tvedöra, where juveniles furthermore showed delayed gonadal differentiation. Here, we show that Dmrt1, a gene that plays a key role in sex determination and sexual development across all metazoans, displays significant sex differentiation in Tvedöra, with a Y-specific haplotype distinct from Ammarnäs. The differential segment is not only much shorter in Tvedöra than in Ammarnäs, it is also less differentiated and associates with both delayed gonadal differentiation and imperfect match between phenotypic and genotypic sex. Whereas Tvedöra juveniles with a local Y haplotype tend to ultimately develop as males, those without it may nevertheless become functional XX males, but with strongly female-biased progeny. Our findings suggest that the variance in patterns of sex determination documented in common frogs might result from a genetic polymorphism within a small genomic region that contains Dmrt1. They also substantiate the view that recurrent convergences of sex determination toward a limited set of chromosome pairs may result from the co-option of small genomic regions that harbor key genes from the sex-determination pathway.

The genetic contribution to sex determination and number of sex chromosomes vary among populations of common frogs (Rana temporaria)

The patterns of sex determination and sex differentiation have been shown to differ among geographic populations of common frogs. Notably, the association between phenotypic sex and linkage group 2 (LG2) has been found to be perfect in a northern Swedish population, but weak and variable among families in a southern one. By analyzing these populations with markers from other linkage groups, we bring two new insights: (1) the variance in phenotypic sex not accounted for by LG2 in the southern population could not be assigned to genetic factors on other linkage groups, suggesting an epigenetic component to sex determination; (2) a second linkage group (LG7) was found to co-segregate with sex and LG2 in the northern population. Given the very short timeframe since post-glacial colonization (in the order of 1000 generations) and its seemingly localized distribution, this neo-sex chromosome system might be the youngest one described so far. It does not result from a fusion, but more likely from a reciprocal translocation between the original Y chromosome (LG2) and an autosome (LG7), causing their co-segregation during male meiosis. By generating a strict linkage between several important genes from the sex-determination cascade (Dmrt1, Amh and Amhr2), this neo-sex chromosome possibly contributes to the 'differentiated sex race' syndrome (strictly genetic sex determination and early gonadal development) that characterizes this northern population.

Sex-chromosome differentiation and 'sex races' in the common frog (Rana temporaria)

Sex-chromosome differentiation was recently shown to vary among common frog populations in Fennoscandia, suggesting a trend of increased differentiation with latitude. By rearing families from two contrasted populations (respectively, from northern and southern Sweden), we show this disparity to stem from differences in sex-determination mechanisms rather than in XY-recombination patterns. Offspring from the northern population display equal sex ratios at metamorphosis, with phenotypic sexes that correlate strongly with paternal LG2 haplotypes (the sex chromosome); accordingly, Y haplotypes are markedly differentiated, with male-specific alleles and depressed diversity testifying to their smaller effective population size. In the southern population, by contrast, a majority of juveniles present ovaries at metamorphosis; only later in development do sex ratios return to equilibrium. Even at these later stages, phenotypic sexes correlate only mildly with paternal LG2 haplotypes; accordingly, there are no recognizable Y haplotypes. These distinct patterns of gonadal development fit the concept of 'sex races' proposed in the 1930s, with our two populations assigned to the 'differentiated' and 'semi-differentiated' races, respectively. Our results support the suggestion that 'sex races' differ in the genetic versus epigenetic components of sex determination. Analysing populations from the 'undifferentiated race' with high-density genetic maps should help to further test this hypothesis.

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.

Pattern of gonadal differentiation and development up to sexual maturity in the frogs, Microhyla ornata and Hylarana malabarica: A comparative study

Gonadogenesis was studied in Microhyla ornata (Family: Microhylidae) and Hylarana malabarica (Family: Ranidae) up to sexual maturity. Indifferent gonads of M. ornata directly differentiated into either testes or ovaries while those of H. malabarica differentiated into ovaries in all the individuals followed by testicular differentiation in males through an ovarian phase. In some tadpoles of M. ornata, formation of a central cavity at Gosner stage 27 marked ovary differentiation while meiosis was initiated at stage 29. Folliculogenesis was evident at stage 39. Vitellogenesis was initiated in females 9 months post-metamorphosis that attained maturity around 11 months after the completion of metamorphosis. Gonads of males with uniformly distributed germ and somatic cells remained undifferentiated until stage 41. Germ and somatic cells reorganized into seminiferous cords at stage 42. One month after completing metamorphosis, testes contained seminiferous tubules while those of 3 months old males exhibited all spermatogenic stages. In H. malabarica, germ cells entry into meiosis marked ovary differentiation at stage 29 while, ovarian cavity was discernable around stage 35. Post-metamorphosis, ovaries of 1-6 month old females contained pre-diplotene oocytes. Females were immature even 1 year after the completion of metamorphosis. In all the tadpoles, ovaries with central cavity and meiocytes were present up to the completion of metamorphosis. Gonads of prospective males displayed an obliterating ovarian cavity along with degenerating oocytes at the end of metamorphosis. Germ and somatic cells reorganized into seminiferous cords in males 3 months after the completion of metamorphosis. Testes of 4 months old males exhibited distinct seminiferous tubules while those of 6 months old exhibited meiosis. All spermatogenic stages were observed in testes of 9 months old males indicating maturity.

Sex-Chromosome Homomorphy in Palearctic Tree Frogs Results from Both Turnovers and X-Y Recombination

Contrasting with birds and mammals, poikilothermic vertebrates often have homomorphic sex chromosomes, possibly resulting from high rates of sex-chromosome turnovers and/or occasional X-Y recombination. Strong support for the latter mechanism was provided by four species of European tree frogs, which inherited from a common ancestor (∼ 5 Ma) the same pair of homomorphic sex chromosomes (linkage group 1, LG1), harboring the candidate sex-determining gene Dmrt1. Here, we test sex linkage of LG1 across six additional species of the Eurasian Hyla radiation with divergence times ranging from 6 to 40 Ma. LG1 turns out to be sex linked in six of nine resolved cases. Mapping the patterns of sex linkage to the Hyla phylogeny reveals several transitions in sex-determination systems within the last 10 My, including one switch in heterogamety. Phylogenetic trees of DNA sequences along LG1 are consistent with occasional X-Y recombination in all species where LG1 is sex linked. These patterns argue against one of the main potential causes for turnovers, namely the accumulation of deleterious mutations on nonrecombining chromosomes. Sibship analyses show that LG1 recombination is strongly reduced in males from most species investigated, including some in which it is autosomal. Intrinsically low male recombination might facilitate the evolution of male heterogamety, and the presence of important genes from the sex-determination cascade might predispose LG1 to become a sex chromosome.

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).

Maintenance of ancestral sex chromosomes in Palearctic tree frogs: direct evidence from Hyla orientalis

Contrasting with the situation found in birds and mammals, sex chromosomes are generally homomorphic in poikilothermic vertebrates. This homomorphy was recently shown to result from occasional X-Y recombinations (not from turnovers) in several European species of tree frogs (Hyla arborea, H. intermedia and H. molleri). Because of recombination, however, alleles at sex-linked loci were rarely diagnostic at the population level; support for sex linkage had to rely on multilocus associations, combined with occasional sex differences in allelic frequencies. Here, we use direct evidence, obtained from anatomical and histological analyses of offspring with known pedigrees, to show that the Eastern tree frog (H. orientalis) shares the same pair of sex chromosomes, with identical patterns of male heterogamety and complete absence of X-Y recombination in males. Conservation of an ancestral pair of sex chromosomes, regularly rejuvenated via occasional X-Y recombination, seems thus a widespread pattern among Hyla species. Sibship analyses also identified discrepancies between genotypic and phenotypic sex among offspring, associated with abnormal gonadal development, suggesting a role for sexually antagonistic genes on the sex chromosomes.