Comparative transcriptomic analysis delineates adaptation strategies of Rana kukunoris toward cold stress on the Qinghai-Tibet Plateau.

BACKGROUND: Cold hardiness is fundamental for amphibians to survive during the extremely cold winter on the Qinghai-Tibet plateau. Exploring the gene regulation mechanism of freezing-tolerant Rana kukunoris could help us to understand how the frogs survive in winter. RESULTS: Transcriptome of liver and muscle of R. kukunoris collected in hibernation and spring were assisted by single molecule real-time (SMRT) sequencing technology. A total of 10,062 unigenes of R. kukunoris were obtained, and 9,924 coding sequences (CDS) were successfully annotated. Our examination of the mRNA response to whole body freezing and recover in the frogs revealed key genes concerning underlying antifreeze proteins and cryoprotectants (glucose and urea). Functional pathway analyses revealed differential regulated pathways of ribosome, energy supply, and protein metabolism which displayed a freeze-induced response and damage recover. Genes related to energy supply in the muscle of winter frogs were up-regulated compared with the muscle of spring frogs. The liver of hibernating frogs maintained modest levels of protein synthesis in the winter. In contrast, the liver underwent intensive high levels of protein synthesis and lipid catabolism to produce substantial quantity of fresh proteins and energy in spring. Differences between hibernation and spring were smaller than that between tissues, yet the physiological traits of hibernation were nevertheless passed down to active state in spring. CONCLUSIONS: Based on our comparative transcriptomic analyses, we revealed the likely adaptive mechanisms of R. kukunoris. Ultimately, our study expands genetic resources for the freezing-tolerant frogs.

Karyotypes of water frogs from the Pelophylax esculentus complex: results of cross-species chromosomal painting.

Amphibian species have the largest genome size enriched with repetitive sequences and relatively similar karyotypes. Moreover, many amphibian species frequently hybridize causing nuclear and mitochondrial genome introgressions. In addition, hybridization in some amphibian species may lead to clonality and polyploidization. All such events were found in water frogs from the genus Pelophylax. Among the species within the genus Pelophylax, P. esculentus complex is the most widely distributed and well-studied. This complex includes two parental species, P. ridibundus and P. lessonae, and their hybrids, P. esculentus, reproducing hemiclonally. Parental species and their hybrids have similar but slightly polymorphic karyotypes, so their precise identification is still required. Here, we have developed a complete set of 13 chromosome painting probes for two parental species allowing the precise identification of all chromosomes. Applying chromosomal painting, we identified homologous chromosomes in both parental species and orthologous chromosomes in their diploid hemiclonal hybrids. Comparative painting did not reveal interchromosomal exchanges between the studied water frog species and their hybrids. Using cross-specific chromosome painting, we detected unequal distribution of the signals along chromosomes suggesting the presence of species-specific tandem repeats. Application of chromosomal paints to the karyotypes of hybrids revealed differences in the intensity of staining for P. ridibundus and P. lessonae chromosomes. Thus, both parental genomes have a divergence in unique sequences. Obtained chromosome probes may serve as a powerful tool to unravel chromosomal evolution in phylogenetically related species, identify individual chromosomes in different cell types, and investigate the elimination of chromosomes in hybrid water frogs.

The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system.

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which species share a conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which mimicry is produced in poison frogs. We assembled a 6.02-Gbp genome with a contig N50 of 310 Kbp, a scaffold N50 of 390 Kbp and 85% of expected tetrapod genes. We leveraged this genome to conduct gene expression analyses throughout development of four colour morphs of Ranitomeya imitator and two colour morphs from both R. fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes differentially expressed throughout development, many of them related to melanophores/melanin, iridophore development and guanine synthesis. We also identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species, and identify several plausible candidates for colouration in vertebrates (e.g. cd36, ep-cadherin and perlwapin). Finally, we hypothesise that keratin genes (e.g. krt8) are important for producing different structural colours within these frogs.

Rapid molecular identification of Rana dybowskii by species-specific primers.

Oviductus Ranae is the dried oviduct from Rana dybowskii, a forest frog species with medicinal, tonic, and cosmetic properties. Due to the high price and resource shortage, counterfeit varieties of Oviductus Ranae often appear in the market. However, traditional identification methods cannot accurately differentiate between Oviductus Ranae and its adulterants. In this study, a rapid molecular identification method has been established. The method involves extracting genomic DNA in just 30 s using filter paper purification, species-specific rapid polymerase chain reaction (PCR) amplification, and finally, fluorescence detection of the products. It can accurately identify Oviductus Ranae and its three common adulterants in about 30 min, making the process simple, fast, and highly specific.

Nuclear gene introgressions in hybrid populations of water frog Pelophylax esculentus complex: geographical analysis of the phenomenon and its interpretation.

Reproduction of water frog hybrids Pelophylax esculentus (Pelophylax ridibundus x Pelophylax lessonae) is associated with hemiclonal reproduction and backcrossing. The hemiclonal mode of reproduction occurs within P. esculentus allodiploids. In this case, the unrecombined genome of one parental species is transmitted to the offspring after premeiotic elimination of the chromosome set of the second parental species. Usually, the chromosome set of P. lessonae is eliminated, and the altered genome of P. ridibundus is passed on to the progeny. The hemiclonal inheritance within diploid Pelophylax esculentus hybrids may be accompanied by certain aberrations of premeiotic elimination. As a result, the formation of P. ridibundus specimens with introgressions of the P. lessonae genetic material, or the formation of recombinant hybrids occurs, depending on which of the parental species backcrossing takes place. The aim of our study is to describe the aberration of premeiotic elimination within the water frog P. esculentus complex detected by the nuclear gene Ldh-B inheritance, with an attempt to find out the causes of this phenomenon. It has been established that aberrations of premeiotic elimination are widespread, but only within populations of water frog from the river system of Upper Dnieper within Ukraine. The highest level of introgression takes place in the water frog populations within Kiev metropolis under conditions of expressed anthropogenization, while the maximum frequency of recombinants was detected within populations from the basin of Desna River, that has preserved native ecosystems. It was demonstrated that the frequency of premeiotic aberrations does not correlate with the intensity of interspecific water frog hybridization. Populations with introgressions are more common than populations with recombinants, however, within the latter, the frequency of recombination events is higher. The primary factor of gametogenesis aberrations, most likely, is the genetic characteristics of the local populations of parental species, since unambiguous explanations of this phenomenon based on the action of environmental stress (pollution of water systems) are not obvious.

Microcystin-LR Induces Estrogenic Effects at Environmentally Relevant Concentration in Black-Spotted Pond Frogs (Pelophylax nigromaculatus): In Situ, In Vivo, In Vitro, and In Silico Investigations.

Harmful cyanobacterial blooms are frequent and intense worldwide, creating hazards for aquatic biodiversity. The potential estrogen-like effect of Microcystin-LR (MC-LR) is a growing concern. In this study, we assessed the estrogenic potency of MC-LR in black-spotted frogs through combined field and laboratory approaches. In 13 bloom areas of Zhejiang province, China, the MC-LR concentrations in water ranged from 0.87 to 8.77 µg/L and were correlated with sex hormone profiles in frogs, suggesting possible estrogenic activity of MC-LR. Tadpoles exposed to 1 µg/L, an environmentally relevant concentration, displayed a female-biased sex ratio relative to controls. Transcriptomic results revealed that MC-LR induces numerous and complex effects on gene expression across multiple endocrine axes. In addition, exposure of male adults significantly increased the estradiol (E2)/testosterone (T) ratio by 3.5-fold relative to controls. Downregulation of genes related to male reproductive endocrine function was also identified. We also showed how MC-LR enhances the expression of specific estrogen receptor (ER) proteins, which induce estrogenic effects by activating the ER pathway and hypothalamic-pituitary-gonadal (HPG) axis. In aggregate, our results reveal multiple lines of evidence demonstrating that, for amphibians, MC-LR is an estrogenic endocrine disruptor at environmentally relevant concentrations. The data presented here support the need for a shift in the MC-LR risk assessment. While hepatoxicity has historically been the focus of MC-LR risk assessments, our data clearly demonstrate that estrogenicity is a major mode of toxicity at environmental levels and that estrogenic effects should be considered for risk assessments on MC-LR going forward.

Effects of atrazine on movement, metabolism and gene expression in Pelophylax nigromaculatus larvae under global warming.

Global warming and environmental pollutants both pose a threat to the behavior and physiology of animals, but research on the combined effects of the two is limited. Atrazine, a widely used herbicide, has toxic effects on organisms. In this study, the effects of environmental concentrations of atrazine exposure (100 µg/L) for seven days on the movement, metabolism and gene expression related to motility of Pelophylax nigromaculatus larvae (GS8) were investigated under global warming. The results showed that compared to the optimal growth temperature (18 °C), atrazine treatment under global warming (21 °C) significantly increased the average speed (about 11.2 times) and maximum acceleration (about 1.98 times) of P. nigromaculatus larvae, altered the relative abundance of 539 metabolites, including Formyl-5-hydroxykynurenamine, 2,4-Dihydroxybenzophenone, and FAPy-adenine, and changed the nucleotide metabolism, pyrimidine metabolism, glycerophospholipid metabolism, and purine metabolism, as well as increased the gene expression of SPLA2 (about 6.46 times) and CHK (about 3.25 times). In summary, atrazine treatment under global warming caused metabolic disorders in amphibian larvae and increased the expression of some movement-related genes in the brain, resulting in abnormally active.

Expression of cholesterol synthesis and steroidogenic markers in females of the Chinese brown frog (Rana dybowskii) during prespawning and prehibernation.

The oviduct of the Chinese brown frog (Rana dybowskii) expands in prehibernation rather than in prespawning, which is one of the physiological phenomena that occur in the preparation for hibernation. Steroid hormones are known to regulate oviductal development. Cholesterol synthesis and steroidogenesis may play an important role in the expansion of the oviduct before hibernation. In this study, we investigated the expression patterns of the markers that are involved in the de novo steroid synthesis pathway in the oviduct of R. dybowskii during prespawning and prehibernation. According to histological analysis, the oviduct of R. dybowskii contains epithelial cells, glandular cells, and tubule lumens. During prehibernation, oviductal pipe diameter and weight were significantly larger than during prespawning. 3-Hydroxy-3-methylglutaryl CoA reductase (HMGCR), low-density lipoprotein receptor (LDLR), steroidogenic acute regulatory protein (StAR), cytochrome P450 cholesterol side-chain cleavage enzyme (P450scc), and steroidogenic factor 1 (SF-1) were detected in epithelial cells in prehibernation and glandular cells during prespawning. HMGCR, LDLR, StAR, and P450scc protein expression levels were higher in prehibernation than during prespawning, but the SF-1 protein expression level did not significantly differ. HMGCR, LDLR, StAR, P450scc (CYP11A1), and SF-1 (NR5A1) mRNA expression levels were significantly higher in prehibernation compared with prespawning. The transcriptome results showed that the steroid synthesis pathway was highly expressed during prehibernation. Existing results indicate that the oviduct is able to synthesize steroid hormones using cholesterol, and that steroid hormones may affect the oviductal functions of R. dybowskii.

Large X-effects are absent in torrent frogs with nascent sex chromosomes.

Sex chromosomes are popularized as a special role in driving speciation. However, the empirical evidence from natural population processes has been limited to organisms with degenerated sex chromosomes, where hemizygosity is mainly considered to act as the driver of reproductive isolation. Here, we examined several hybrid zones of torrent frog Amolops mantzorum species complex, using an approach by mapping species-diagnostic loci onto the reference genome to compare sex-linked versus autosomal patterns of introgression. We find little support in sex-linked incompatibilities for large X-effects for these populations in hybrid zones with homomorphic sex chromosomes, due to the absence of the hemizygous effects. As expected, the large X-effects were not found in those with heteromorphic but newly evolved sex chromosomes, owing to the absence of strong genetic differences between X and Y chromosomes. The available data so far on amphibians suggest little role for sex-linked genes in speciation. The large X-effects in those with nascent sex chromosomes may not be as ubiquitous as presumed across the animal kingdom.

Intestinal response of Rana chensinensis larvae exposed to Cr and Pb, alone and in combination.

Although numerous investigations on the adverse impact of Cr and Pb have been performed, studies on intestinal homeostasis in amphibians are limited. Here, single and combined effects of Cr (104 µg/L) and Pb (50 µg/L) on morphological and histological features, bacterial community, digestive enzymes activities, as well as transcriptomic profile of intestines in Rana chensinensis tadpoles were assessed. Significant decrease in the relative intestine length (intestine length/snout-to-vent length, IL/SVL) was observed after exposure to Pb and Cr/Pb mixture. Intestinal histology and digestive enzymes activities were altered in metal treatment groups. In addition, treatment groups showed significantly increased bacterial richness and diversity. Tadpoles in treatment groups were observed to have differential gut bacterial composition from controls, especially for the abundance of phylum Proteobacteria, Firmicutes, Verrucomicrobia, Actinobacteria, and Fusobacteria as well as genus Citrobacter, Anaerotruncus, Akkermansia, and Alpinimonas. Moreover, transcriptomic analysis showed that the transcript expression profiles of GPx and SOD isoforms responded differently to Cr and/or Pb exposure. Besides, transcriptional activation of pro-apoptotic and glycolysis-related genes, such as Bax, Apaf 1, Caspase 3, PK, PGK, TPI, and GPI were detected in all treatment groups but downregulation of Bcl2 in Pb and Cr/Pb mixture groups. Collectively, these results suggested that Cr and Pb exposure at environmental relevant concentration, alone and in combination, could disrupt intestinal homeostasis of R. chensinensis tadpoles.

Genomics-informed delineation of conservation units in a desert amphibian.

Delineating conservation units (CUs, e.g., evolutionarily significant units, ESUs, and management units, MUs) is critical to the recovery of declining species because CUs inform both listing status and management actions. Genomic data have strengths and limitations in informing CU delineation and related management questions in natural systems. We illustrate the value of using genomic data in combination with landscape, dispersal and occupancy data to inform CU delineation in Nevada populations of the Great Basin Distinct Population Segment of the Columbia spotted frog (Rana luteiventris). R. luteiventris occupies naturally fragmented aquatic habitats in this xeric region, but beaver removal, climate change and other factors have put many of these populations at high risk of extirpation without management intervention. We addressed three objectives: (i) assessing support for ESUs within Nevada; (ii) evaluating and revising, if warranted, the current delineation of MUs; and (iii) evaluating genetic diversity, effective population size, adaptive differentiation and functional connectivity to inform ongoing management actions. We found little support for ESUs within Nevada but did identify potential revisions to MUs based on unique landscape drivers of connectivity that distinguish these desert populations from those in the northern portion of the species range. Effective sizes were uniformly small, with low genetic diversity and weak signatures of adaptive differentiation. Our findings suggest that management actions, including translocations and genetic rescue, might be warranted. Our study illustrates how a carefully planned genetic study, designed to address priority management goals that include CU delineation, can provide multiple insights to inform conservation action.

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

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

Giant Tree Frog diversification in West and Central Africa: Isolation by physical barriers, climate, and reproductive traits.

Secondary sympatry amongst sister lineages is strongly associated with genetic and ecological divergence. This pattern suggests that for closely related species to coexist in secondary sympatry, they must accumulate differences in traits that mediate ecological and/or reproductive isolation. Here, we characterized inter- and intraspecific divergence in three giant tree frog species whose distributions stretch across West and Central Africa. Using genome-wide single-nucleotide polymorphism data, we demonstrated that species-level divergence coincides temporally and geographically with a period of large-scale forest fragmentation during the late Pliocene. Our environmental niche models further supported a dynamic history of climatic suitability and stability, and indicated that all three species occupy distinct environmental niches. We found modest morphological differentiation amongst the species with significant divergence in tympanum diameter and male advertisement call. In addition, we confirmed that two species occur in secondary sympatry in Central Africa but found no evidence of hybridization. These patterns support the hypothesis that cycles of genetic exchange and isolation across West and Central Africa have contributed to globally significant biodiversity. Furthermore, divergence in both ecology and reproductive traits appear to have played important roles in maintaining distinct lineages. At the intraspecific level, we found that climatic refugia, precipitation gradients, marine incursions, and potentially riverine barriers generated phylogeographic structure throughout the Pleistocene and into the Holocene. Further studies examining phenotypic divergence and secondary contact amongst these geographically structured populations may demonstrate how smaller scale and more recent biogeographic barriers contribute to regional diversification.

La sympatrie secondaire parmi les espèces sœurs est fortement associée à la divergence génétique et écologique. Ce modèle suggère que pour que des espèces étroitement liées coexistent en sympatrie secondaire, elles doivent accumuler des différences dans les traits qui contribuent à l'isolement écologique ou reproductif. Ici, nous avons caractérisé la divergence inter- et intra-spécifique chez trois espèces de grenouilles arboricoles géantes dont les distributions s'étendent à travers l'Afrique de l'Ouest et Centrale. Avec des données génétiques, nous avons démontré que la divergence au niveau des espèces coïncide temporellement et géographiquement avec une période de fragmentation forestière à la fin du Pliocène. Nos modèles de niches environnementales ont soutenu une histoire dynamique de stabilité climatique, et ont indiqué que les trois espèces occupent des niches environnementales distinctes. Nous avons trouvé une différenciation morphologique modeste parmi les trois espèces mais une divergence significative dans le diamètre du tympan et les cris des mâles. De plus, nous avons confirmé que deux espèces sont présentes en sympatrie secondaire en Afrique Centrale mais n'avons trouvé aucune preuve d'hybridation. Ces résultats soutiennent l'hypothèse que les cycles d'échange génétique et d'isolement à travers l'Afrique de l'Ouest et Centrale ont contribué à une profonde concentration de biodiversité dans la région. De plus, la divergence des traits écologiques et reproducteurs semble avoir joué un rôle important dans le maintien de lignées distinctes. Au niveau intra-spécifique, nous avons constaté que les refuges climatiques, les gradients de précipitation, les incursions marines et potentiellement les barrières fluviales ont généré une structure phylogéographique pendant le Pléistocène et jusqu'à l'Holocène. Des études examinant la divergence phénotypique et le contact secondaire entre ces populations géographiquement structurées pourraient démontrer comment des barrières biogéographiques à échelle plus petite et plus récentes contribuent à la diversification régionale.

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

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

Sequencing of laser captured Z and W chromosomes of the tocantins paradoxical frog (Pseudis tocantins) provides insights on repeatome and chromosomal homology.

Pseudis tocantins is the only frog species of the hylid genus Pseudis that possesses highly heteromorphic sex chromosomes. Z and W chromosomes of Ps. tocantins differ in size, morphology, position of the nucleolar organizer region (NOR) and the amount and distribution of heterochromatin. A chromosomal inversion and heterochromatin amplification on the W chromosome were previously inferred to be involved in the evolution of this sex chromosome pair. Despite these findings, knowledge related to the molecular composition of the large heterochromatic band of this W chromosome is restricted to the PcP190 satellite DNA, and no data are available regarding the gene content of either the W or the Z chromosome of Ps. tocantins. Here, we sequenced microdissected Z and W chromosomes of this species to further resolve their molecular composition. Comparative genomic analysis suggests that Ps. tocantins sex chromosomes are likely homologous to chromosomes 4 and 10 of Xenopus tropicalis. Analyses of the repetitive DNA landscape in the Z and W assemblies allowed for the identification of several transposable elements and putative satellite DNA sequences. Finally, some transposable elements from the W assembly were found to be highly diverse and divergent from elements found elsewhere in the genome, suggesting a rapid amplification of these elements on the W chromosome.

Heterogeneous Evolution of Sex Chromosomes in the Torrent Frog Genus Amolops.

In sharp contrast to birds and mammals, in numerous cold-blooded vertebrates, sex chromosomes have been described as homomorphic. This sex chromosome homomorphy has been suggested to result from the high turnovers often observed across deeply diverged clades. However, little is known about the tempo and mode of sex chromosome evolution among the most closely related species. Here, we examined the evolution of sex chromosome among nine species of the torrent frog genus Amolops. We analyzed male and female GBS and RAD-seq from 182 individuals and performed PCR verification for 176 individuals. We identified signatures of sex chromosomes involving two pairs of chromosomes. We found that sex-chromosome homomorphy results from both turnover and X-Y recombination in the Amolops species, which simultaneously exhibits heterogeneous evolution on homologous and non-homologous sex chromosomes. A low turnover rate of non-homologous sex chromosomes exists in these torrent frogs. The ongoing X-Y recombination in homologous sex chromosomes will act as an indispensable force in preventing sex chromosomes from differentiating.

First natural breeding of the endangered dusky gopher frog (Lithobates sevosus) in captivity.

To date, all captive breeding of the dusky gopher frog, Lithobates sevosus, a federally listed endangered species, has been accomplished using in vitro fertilization (IVF). Here, we describe multievent and highly fecund captive reproduction of dusky gopher frogs driven solely by natural environmental factors. Six pairs of L. sevosus were kept in a 3.7 × 4.4 m2 outdoor enclosure designed to resemble their natural breeding habitat, which included a pool and three artificial burrows. Modifications to the enclosure that simulated temperatures and conditions within their natural range during winter were added in October and removed in late February. Following a warm, rainy period, five egg masses were laid between March 5 and 11, 2020. The number of oocytes per egg mass was 2300 ± 409 (range = 1341-3565), with the total across all five egg masses being 11,501. Of these oocytes, the hatching rate was 68.58 ± 10.05% (range = 37.53%-95.59%), with a total of 7887 successful hatchlings overall. Clutch sizes were similar to those in the wild and greater than those typically produced using IVF; thus, natural breeding can substantially increase the number of frogs available for reintroduction programs. Although assisted reproductive technologies such as IVF will continue to be useful for ensuring the success of strategic genetic pairings of captive L. sevosus, the new tool of nonassisted reproduction in specifically designed outdoor enclosures is an important advancement for the conservation and recovery of this endangered species.

Evolution of Sex Chromosome Heteromorphy in Geographic Populations of the Japanese Tago's Brown Frog Complex.

The sex chromosomes of most anuran amphibians are characterized by homomorphy in both sexes, and evolution to heteromorphy rarely occurs at the species or geographic population level. Here, we report sex chromosome heteromorphy in geographic populations of the Japanese Tago's brown frog complex (2n = 26), comprising Rana sakuraii and R. tagoi. The sex chromosomes of R. sakuraii from the populations in western Japan were homomorphic in both sexes, whereas chromosome 7 from the populations in eastern Japan were heteromorphic in males. Chromosome 7 of R. tagoi, which is distributed close to R. sakuraii in eastern Japan, was highly similar in morphology to the Y chromosome of R. sakuraii. Based on this and on mitochondrial gene sequence analysis, we hypothesize that in the R. sakuraii populations from eastern Japan the XY heteromorphic sex chromosome system was established by the introduction of chromosome 7 from R. tagoi via interspecies hybridization. In contrast, chromosome 13 of R. tagoi from the 2 large islands in western Japan, Shikoku and Kyushu, showed a heteromorphic pattern of constitutive heterochromatin distribution in males, while this pattern was homomorphic in females. Our study reveals that sex chromosome heteromorphy evolved independently at the geographic lineage level in this species complex.

Ecological adaptation drives wood frog population divergence in life history traits.

Phenotypic variation among populations is thought to be generated from spatial heterogeneity in environments that exert selection pressures that overcome the effects of gene flow and genetic drift. Here, we tested for evidence of isolation by distance or by ecology (i.e., ecological adaptation) to generate variation in early life history traits and phenotypic plasticity among 13 wood frog populations spanning 1200 km and 7° latitude. We conducted a common garden experiment and related trait variation to an ecological gradient derived from an ecological niche model (ENM) validated to account for population density variation. Shorter larval periods, smaller body weight, and relative leg lengths were exhibited by populations with colder mean annual temperatures, greater precipitation, and less seasonality in precipitation and higher population density (high-suitability ENM values). After accounting for neutral genetic variation, the QST-FST analysis supported ecological selection as the key process generating population divergence. Further, the relationship between ecology and traits was dependent upon larval density. Specifically, high-suitability/high-density populations in the northern part of the range were better at coping with greater conspecific competition, evidenced by greater postmetamorphic survival and no difference in body weight when reared under stressful conditions of high larval density. Our results support that both climate and competition selection pressures drive clinal variation in larval and metamorphic traits in this species. Range-wide studies like this one are essential for accurate predictions of population's responses to ongoing ecological change.

Small-scale population divergence is driven by local larval environment in a temperate amphibian.

Genomic variation within and among populations is shaped by the interplay between natural selection and the effects of genetic drift and gene flow. Adaptive divergence can be found in small-scale natural systems even when population sizes are small, and the potential for gene flow is high, suggesting that local environments exert selection pressures strong enough to counteract the opposing effects of drift and gene flow. Here, we investigated genomic differentiation in nine moor frog (Rana arvalis) populations in a small-scale network of local wetlands using 16,707 ddRAD-seq SNPs, relating levels of differentiation with local environments, as well as with properties of the surrounding landscape. We characterized population structure and differentiation, and partitioned the effects of geographic distance, local larval environment, and landscape features on total genomic variation. We also conducted gene-environment association studies using univariate and multivariate approaches. We found small-scale population structure corresponding to 6-8 clusters. Local larval environment was the most influential component explaining 2.3% of the total genetic variation followed by landscape features (1.8%) and geographic distance (0.8%), indicative of isolation-by-environment, -by-landscape, and -by-distance, respectively. We identified 1000 potential candidate SNPs putatively under divergent selection mediated by the local larval environment. The candidate SNPs were involved in, among other biological functions, immune system function and development. Our results suggest that small-scale environmental differences can exert selection pressures strong enough to counteract homogenizing effects of gene flow and drift in this small-scale system, leading to observable population differentiation.