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.

Differential Mitochondrial Genome Expression of Four Hylid Frog Species under Low-Temperature Stress and Its Relationship with Amphibian Temperature Adaptation.

Extreme weather poses huge challenges for animals that must adapt to wide variations in environmental temperature and, in many cases, it can lead to the local extirpation of populations or even the extinction of an entire species. Previous studies have found that one element of amphibian adaptation to environmental stress involves changes in mitochondrial gene expression at low temperatures. However, to date, comparative studies of gene expression in organisms living at extreme temperatures have focused mainly on nuclear genes. This study sequenced the complete mitochondrial genomes of five Asian hylid frog species: Dryophytes japonicus, D. immaculata, Hyla annectans, H. chinensis and H. zhaopingensis. It compared the phylogenetic relationships within the Hylidae family and explored the association between mitochondrial gene expression and evolutionary adaptations to cold stress. The present results showed that in D. immaculata, transcript levels of 12 out of 13 mitochondria genes were significantly reduced under cold exposure (p < 0.05); hence, we put forward the conjecture that D. immaculata adapts by entering a hibernation state at low temperature. In H. annectans, the transcripts of 10 genes (ND1, ND2, ND3, ND4, ND4L, ND5, ND6, COX1, COX2 and ATP8) were significantly reduced in response to cold exposure, and five mitochondrial genes in H. chinensis (ND1, ND2, ND3, ND4L and ATP6) also showed significantly reduced expression and transcript levels under cold conditions. By contrast, transcript levels of ND2 and ATP6 in H. zhaopingensis were significantly increased at low temperatures, possibly related to the narrow distribution of this species primarily at low latitudes. Indeed, H. zhaopingensis has little ability to adapt to low temperature (4 °C), or maybe to enter into hibernation, and it shows metabolic disorder in the cold. The present study demonstrates that the regulatory trend of mitochondrial gene expression in amphibians is correlated with their ability to adapt to variable climates in extreme environments. These results can predict which species are more likely to undergo extirpation or extinction with climate change and, thereby, provide new ideas for the study of species extinction in highly variable winter climates.

Diversity and Molecular Evolution of Nonvisual Opsin Genes across Environmental, Developmental, and Morphological Adaptations in Frogs.

Nonvisual opsins are transmembrane proteins expressed in the eyes and other tissues of many animals. When paired with a light-sensitive chromophore, nonvisual opsins form photopigments involved in various nonvisual, light-detection functions including circadian rhythm regulation, light-seeking behaviors, and seasonal responses. Here, we investigate the molecular evolution of nonvisual opsin genes in anuran amphibians (frogs and toads). We test several evolutionary hypotheses including the predicted loss of nonvisual opsins due to nocturnal ancestry and potential functional differences in nonvisual opsins resulting from environmental light variation across diverse anuran ecologies. Using whole-eye transcriptomes of 81 species, combined with genomes, multitissue transcriptomes, and independently annotated genes from an additional 21 species, we identify which nonvisual opsins are present in anuran genomes and those that are also expressed in the eyes, compare selective constraint among genes, and test for potential adaptive evolution by comparing selection between discrete ecological classes. At the genomic level, we recovered all 18 ancestral vertebrate nonvisual opsins, indicating that anurans demonstrate the lowest documented amount of opsin gene loss among ancestrally nocturnal tetrapods. We consistently found expression of 14 nonvisual opsins in anuran eyes and detected positive selection in a subset of these genes. We also found shifts in selective constraint acting on nonvisual opsins in frogs with differing activity periods, habitats, distributions, life histories, and pupil shapes, which may reflect functional adaptation. Although many nonvisual opsins remain poorly understood, these findings provide insight into the diversity and evolution of these genes across anurans, filling an important gap in our understanding of vertebrate opsins and setting the stage for future research on their functional evolution across taxa.

Genome Assembly of the Dyeing Poison Frog Provides Insights into the Dynamics of Transposable Element and Genome-Size Evolution.

Genome size varies greatly across the tree of life and transposable elements are an important contributor to this variation. Among vertebrates, amphibians display the greatest variation in genome size, making them ideal models to explore the causes and consequences of genome size variation. However, high-quality genome assemblies for amphibians have, until recently, been rare. Here, we generate a high-quality genome assembly for the dyeing poison frog, Dendrobates tinctorius. We compare this assembly to publicly available frog genomes and find evidence for both large-scale conserved synteny and widespread rearrangements between frog lineages. Comparing conserved orthologs annotated in these genomes revealed a strong correlation between genome size and gene size. To explore the cause of gene-size variation, we quantified the location of transposable elements relative to gene features and find that the accumulation of transposable elements in introns has played an important role in the evolution of gene size in D. tinctorius, while estimates of insertion times suggest that many insertion events are recent and species-specific. Finally, we carry out population-scale mobile-element sequencing and show that the diversity and abundance of transposable elements in poison frog genomes can complicate genotyping from repetitive element sequence anchors. Our results show that transposable elements have clearly played an important role in the evolution of large genome size in D. tinctorius. Future studies are needed to fully understand the dynamics of transposable element evolution and to optimize primer or bait design for cost-effective population-level genotyping in species with large, repetitive genomes.

Genetic Variants Underlying Plasticity in Natural Populations of Spadefoot Toads: Environmental Assessment versus Phenotypic Response.

Many organisms facultatively produce different phenotypes depending on their environment, yet relatively little is known about the genetic bases of such plasticity in natural populations. In this study, we describe the genetic variation underlying an extreme form of plasticity--resource polyphenism--in Mexican spadefoot toad tadpoles, Spea multiplicata. Depending on their environment, these tadpoles develop into one of two drastically different forms: a carnivore morph or an omnivore morph. We collected both morphs from two ponds that differed in which morph had an adaptive advantage and performed genome-wide association studies of phenotype (carnivore vs. omnivore) and adaptive plasticity (adaptive vs. maladaptive environmental assessment). We identified four quantitative trait loci associated with phenotype and nine with adaptive plasticity, two of which exhibited signatures of minor allele dominance and two of which (one phenotype locus and one adaptive plasticity locus) did not occur as minor allele homozygotes. Investigations into the genetics of plastic traits in natural populations promise to provide novel insights into how such complex, adaptive traits arise and evolve.

Understanding the genetic sex-determining mechanism in Hyla eximia treefrog inferred from H-Y antigen.

Genetic sex-determining mechanisms have been extensively elucidated in mammals; however, the sex chromosomes, sex-determining genes, and gene regulatory networks involved in sex differentiation remain poorly understood in amphibians. In this study, we investigated the sex-determining mechanism in the Hyla eximia treefrog based on karyotypic analysis and identification of H-Y antigen, a sex-linked peptide that is present in the gonads of the heterogametic sex (XY or ZW) in all vertebrates. Results show a diploid chromosome number 2n = 24 with homomorphic sex chromosomes. The heterogametic sex, ZW-female, were hypothesized based on H-Y antigen mRNA expression in female gonads (24,ZZ/24,ZW). The treefrog H-Y peptide exhibited a high percentage of identity with other vertebrate sequences uploaded to GenBank database. To obtain gene expression profiles, we also obtained the coding sequence of the housekeeping Actb gene. High H-Y antigen expression levels were further confirmed in ovaries using real-time polymerase chain reaction (RT-PCR) during non-breeding season, we noted a decrease in the expression of the H-Y antigen during breeding season. This study provides evidence that sex hormones might suppress H-Y antigen expression in the gonads of heterogametic females 24,ZW during the breeding season. These findings suggest that H-Y gene expression is a well-suited model for studying heterogametic sex by comparing the male and female gonads.

Effects of hibernation on two important contractile tissues in tibetan frogs, Nanorana parkeri: a perspective from transcriptomics and metabolomics approaches.

BACKGROUND: In response to seasonal cold and food shortage, the Xizang plateau frogs, Nanorana parkeri (Anura: Dicroglossidae), enter a reversible hypometabolic state where heart rate and oxygen consumption in skeletal muscle are strongly suppressed. However, the effect of winter hibernation on gene expression and metabolic profiling in these two tissues remains unknown. In the present study, we conducted transcriptomic and metabolomic analyses of heart and skeletal muscle from summer- and winter-collected N. parkeri to explore mechanisms involved in seasonal hibernation. RESULTS: We identified 2407 differentially expressed genes (DEGs) in heart and 2938 DEGs in skeletal muscle. Enrichment analysis showed that shared DEGs in both tissues were enriched mainly in translation and metabolic processes. Of these, the expression of genes functionally categorized as "response to stress", "defense mechanisms", or "muscle contraction" were particularly associated with hibernation. Metabolomic analysis identified 24 and 22 differentially expressed metabolites (DEMs) in myocardium and skeletal muscle, respectively. In particular, pathway analysis showed that DEMs in myocardium were involved in the pentose phosphate pathway, glycerolipid metabolism, pyruvate metabolism, citrate cycle (TCA cycle), and glycolysis/gluconeogenesis. By contrast, DEMs in skeletal muscle were mainly involved in amino acid metabolism. CONCLUSIONS: In summary, natural adaptations of myocardium and skeletal muscle in hibernating N. parkeri involved transcriptional alterations in translation, stress response, protective mechanisms, and muscle contraction processes as well as metabolic remodeling. This study provides new insights into the transcriptional and metabolic adjustments that aid winter survival of high-altitude frogs N. parkeri.

Contrasting nidification behaviors facilitate diversification and colonization of the Music frogs under a changing paleoclimate.

In order to cope with the complexity and variability of the terrestrial environment, amphibians have developed a wide range of reproductive and parental behaviors. Nest building occurs in some anuran species as parental care. Species of the Music frog genus Nidirana are known for their unique courtship behavior and mud nesting in several congeners. However, the evolution of these frogs and their nidification behavior has yet to be studied. With phylogenomic and phylogeographic analyses based on a wide sampling of the genus, we find that Nidirana originated from central-southwestern China and the nidification behavior initially evolved at ca 19.3 Ma but subsequently lost in several descendants. Further population genomic analyses suggest that the nidification species have an older diversification and colonization history, while N. adenopleura complex congeners that do not exhibit nidification behavior have experienced a recent rapid radiation. The presence and loss of the nidification behavior in the Music frogs may be associated with paleoclimatic factors such as temperature and precipitation. This study highlights the nidification behavior as a key evolutionary innovation that has contributed to the diversification of an amphibian group under past climate changes.

Diversity and Evolution of Frog Visual Opsins: Spectral Tuning and Adaptation to Distinct Light Environments.

Visual systems adapt to different light environments through several avenues including optical changes to the eye and neurological changes in how light signals are processed and interpreted. Spectral sensitivity can evolve via changes to visual pigments housed in the retinal photoreceptors through gene duplication and loss, differential and coexpression, and sequence evolution. Frogs provide an excellent, yet understudied, system for visual evolution research due to their diversity of ecologies (including biphasic aquatic-terrestrial life cycles) that we hypothesize imposed different selective pressures leading to adaptive evolution of the visual system, notably the opsins that encode the protein component of the visual pigments responsible for the first step in visual perception. Here, we analyze the diversity and evolution of visual opsin genes from 93 new eye transcriptomes plus published data for a combined dataset spanning 122 frog species and 34 families. We find that most species express the four visual opsins previously identified in frogs but show evidence for gene loss in two lineages. Further, we present evidence of positive selection in three opsins and shifts in selective pressures associated with differences in habitat and life history, but not activity pattern. We identify substantial novel variation in the visual opsins and, using microspectrophotometry, find highly variable spectral sensitivities, expanding known ranges for all frog visual pigments. Mutations at spectral-tuning sites only partially account for this variation, suggesting that frogs have used tuning pathways that are unique among vertebrates. These results support the hypothesis of adaptive evolution in photoreceptor physiology across the frog tree of life in response to varying environmental and ecological factors and further our growing understanding of vertebrate visual evolution.

Marine introgressions and Andean uplift have driven diversification in neotropical Monkey tree frogs (Anura, Phyllomedusinae).

The species richness in the Neotropics has been linked to environmental heterogeneity and a complex geological history. We evaluated which biogeographic processes were associated with the diversification of Monkey tree frogs, an endemic clade from the Neotropics. We tested two competing hypotheses: the diversification of Phyllomedusinae occurred either in a "south-north" or a "north-south" direction in the Neotropics. We also hypothesized that marine introgressions and Andean uplift had a crucial role in promoting their diversification. We used 13 molecular markers in a Bayesian analysis to infer phylogenetic relationships among 57 species of Phyllomedusinae and to estimate their divergence times. We estimated ancestral ranges based on 12 biogeographic units considering the landscape modifications of the Neotropical region. We found that the Phyllomedusinae hypothetical ancestor range was probably widespread throughout South America, from Western Amazon to Southern Atlantic Forest, at 29.5 Mya. The Phyllomedusines' ancestor must have initially diverged through vicariance, generally followed by jump-dispersals and sympatric speciation. Dispersal among areas occurred mostly from Western Amazonia towards Northern Andes and the South American diagonal of dry landscapes, a divergent pattern from both "south-north" and "north-south" diversification hypotheses. Our results revealed a complex diversification process of Monkey tree frogs, occurring simultaneously with the orogeny of Northern Andes and the South American marine introgressions in the last 30 million years.

Distribution and diversification of Adelphobates, emblematic poison frogs from Brazilian Amazonia.

Adelphobates contains three species, and the inaccurate identification of A. quinquevittatus and the scarcity of records of A. castaneoticus complicate inference of their distributions; the latter species occurs in sympatry with A. galactonotus. Our objective was to revise the distributions of Adelphobates by compiling data and modeling habitat suitability, as range limits may be shaped by landscape features and biotic interactions. We initially analyzed the existence of operational taxonomic units within the nominal species and subsequently inferred the observed and potential distributions, taking into account the possible independent lineages for the three species, and we also generated a molecular timetree to understand the chronology of interspecific diversification events. Adelphobates quinquevittatus was found to have a more easterly distribution than previously described, and specimens with phenotypic variation were found to occur in areas inconsistent with the modeling, and A. castaneoticus was concentrated in the Tapajós-Xingu interfluve, surrounded by A. galactonotus. Models indicated that the right bank of the Xingu River is suitable for both species, indeed, both were found there. Despite Adelphobates species having their distributions delimited by major Amazonian rivers, estimated divergence times predate the formation of the modern river network, suggesting that other mechanisms were involved in their diversification.

Toxic to the touch: The makings of lethal mantles in pitohui birds and poison dart frogs.

How do chemically defended animals resist their own toxins? This intriguing question on the concept of autotoxicity is at the heart of how species interactions evolve. In this issue of Molecular Ecology (Molecular Ecology, 2024, 33), Bodawatta and colleagues report on how Papua New Guinean birds coopted deadly neurotoxins to create lethal mantles that protect against predators and parasites. Combining chemical screening of the plumage of a diverse collection of passerine birds with genome sequencing, the researchers unlocked a deeper understanding of how some birds sequester deadly batrachotoxin (BTX) from their food without poisoning themselves. They identified that birds impervious to BTX bear amino acid substitutions in the toxin-binding site of the voltage-gated sodium channel Nav1.4, whose function is essential for proper contraction and relaxation of vertebrate muscles. Comparative genetic and molecular docking analyses show that several of the substitutions associated with insensitivity to BTX may have become prevalent among toxic birds through positive selection. Intriguingly, poison dart frogs that also co-opted BTX in their lethal mantles were found to harbour similar toxin insensitivity substitutions in their Nav1.4 channels. Taken together, this sets up a powerful model system for studying the mechanisms behind convergent molecular evolution and how it may drive biological diversity.

Multiple Transitions between Y Chromosome and Autosome in Tago's Brown Frog Species Complex.

Sex chromosome turnover is the transition between sex chromosomes and autosomes. Although many cases have been reported in poikilothermic vertebrates, their evolutionary causes and genetic mechanisms remain unclear. In this study, we report multiple transitions between the Y chromosome and autosome in the Japanese Tago's brown frog complex. Using chromosome banding and molecular analyses (sex-linked and autosomal single nucleotide polymorphisms, SNPs, from the nuclear genome), we investigated the frogs of geographic populations ranging from northern to southern Japan of two species, Rana tagoi and Rana sakuraii (2n = 26). Particularly, the Chiba populations of East Japan and Akita populations of North Japan in R. tagoi have been, for the first time, investigated here. As a result, we identified three different sex chromosomes, namely chromosomes 3, 7, and 13, in the populations of the two species. Furthermore, we found that the transition between the Y chromosome (chromosome 7) and autosome was repeated through hybridization between two or three different populations belonging to the two species, followed by restricted chromosome introgression. These dynamic sex chromosome turnovers represent the first such findings in vertebrates and imply that speciation associated with inter- or intraspecific hybridization plays an important role in sex chromosome turnover in frogs.

Deep divergences among inconspicuously colored clades of Epipedobates poison frogs.

Poison frogs (Dendrobatidae) are famous for their aposematic species, having a combination of diverse color patterns and defensive skin toxins, yet most species in this family are inconspicuously colored and considered non-aposematic. Epipedobates is among the youngest genus-level clades of Dendrobatidae that includes both aposematic and inconspicuous species. Using Sanger-sequenced mitochondrial and nuclear markers, we demonstrate deep genetic divergences among inconspicuous species of Epipedobates but relatively shallow genetic divergences among conspicuous species. Our phylogenetic analysis includes broad geographic sampling of the inconspicuous lineages typically identified as E. boulengeri and E. espinosai, which reveals two putative new species, one in west-central Colombia (E. sp. 1) and the other in north-central Ecuador (E. aff. espinosai). We conclude that E. darwinwallacei is a junior subjective synonym of E. espinosai. We also clarify the geographic distributions of inconspicuous Epipedobates species including the widespread E. boulengeri. We provide a qualitative assessment of the phenotypic diversity in each nominal species, with a focus on the color and pattern of inconspicuous species. We conclude that Epipedobates contains eight known valid species, six of which are inconspicuous. A relaxed molecular clock analysis suggests that the most recent common ancestor of Epipedobates is âˆ¼11.1 million years old, which nearly doubles previous estimates. Last, genetic information points to a center of species diversity in the Chocó at the southwestern border of Colombia with Ecuador. A Spanish translation of this text is available in the supplementary materials.

Piecing the barcoding puzzle of Palearctic water frogs (Pelophylax) sheds light on amphibian biogeography and global invasions.

Palearctic water frogs (genus Pelophylax) are an outstanding model in ecology and evolution, being widespread, speciose, either threatened or threatening to other species through biological invasions, and capable of siring hybrid offspring that escape the rules of sexual reproduction. Despite half a century of genetic research and hundreds of publications, the diversity, systematics and biogeography of Pelophylax still remain highly confusing, in no small part due to a lack of correspondence between studies. To provide a comprehensive overview, we gathered >13,000 sequences of barcoding genes from >1700 native and introduced localities and built multigene mitochondrial (~17 kb) and nuclear (~10 kb) phylogenies. We mapped all currently recognized taxa and their phylogeographic lineages (>40) to get a grasp on taxonomic issues, cyto-nuclear discordances, the genetic makeup of hybridogenetic hybrids, and the origins of introduced populations. Competing hypotheses for the molecular calibration were evaluated through plausibility tests, implementing a new approach relying on predictions from the anuran speciation continuum. Based on our timetree, we propose a new biogeographic paradigm for the Palearctic since the Paleogene, notably by attributing a prominent role to the dynamics of the Paratethys, a vast paleo-sea that extended over most of Europe. Furthermore, our results show that distinct marsh frog lineages from Eastern Europe, the Balkans, the Near East, and Central Asia (P. ridibundus ssp.) are naturally capable of inducing hybridogenesis with pool frogs (P. lessonae). We identified 14 alien lineages (mostly of P. ridibundus) over ~20 areas of invasions, especially in Western Europe, with genetic signatures disproportionally pointing to the Balkans and Anatolia as the regions of origins, in line with exporting records of the frog leg industry and the stocks of pet sellers. Pelophylax thus emerges as one of the most invasive amphibians worldwide, and deserves much higher conservation concern than currently given by the authorities fighting biological invasions.

Lost in time: Unraveling the identity of Vitreorana parvula (Boulenger, 1895) (Anura: Centrolenidae).

Vitreorana parvula was the first glassfrog described for the Atlantic Forest. The species, however, has become a taxonomic puzzle as the only known individual is the lectotype from the 19th century, which is not particularly well-preserved or accompanied by a detailed original description. To solve this problem, we collected topotypic specimens, as well as advertisement calls, tissue samples, and natural history data, and compared them to other Vitreorana species. Our results show clear morphological, acoustic, and genetic differences between V. parvula and other species of Vitreorana, except for V. uranoscopa. Following our results, we consider V. uranoscopa as a junior synonym of V. parvula and redescribe the species based on topotypic material, while summarizing relevant variation from across its distribution.

Discovery and description of a new species of burrowing frog Sphaerotheca Gnther, 1859 (Anura: Dicroglossidae) from the suburban landscapes of Bengaluru, India.

Recent resolution of prevailing taxonomic ambiguities in the genus Sphaerotheca and new species discoveries from urban/suburban landscapes highlight the need for attention to non-forested habitats for amphibian conservation. In this paper, we review the status of the members of the genus Sphaerotheca and justify the synonymy of Sphaerotheca magadha as a junior synonym of Sphaerotheca swani. The prospects of resurrection of Sphaerotheca swani (herein preliminarily referred to as Sphaerotheca cf. breviceps [swani]) are discussed. In addition, we describe a new species Sphaerotheca varshaabhu sp. nov. from the suburban region of Bengaluru, India. We employ an integrative taxonomic approach to characterize the new species using molecular phylogeny, genetic distance, morphological characters, and geographical isolation as lines of evidence. We also provide a description of vocal repertoire of Sphaerotheca varshaabhu sp. nov. and provide comparative bioacoustics data for four species. This previously undescribed species from the suburban areas of Bengaluru described herein as Sphaerotheca varshaabhu sp. nov. forms a genetically divergent lineage and its genetic distance varied from 3.6% to 12.2% for 16S rRNA with respect to other species of Sphaerotheca. Our phylogenetic analysis for the genus including the new species confirms the synonymy of one recently described species, resulting in 10 valid species in the genus Sphaerotheca. These results emphasize the need for utilizing an integrative taxonomic approach for uncovering hidden diversity of suburban areas. Given these recent discoveries, we advocate for more robust surveys in human dominated areas, so that these amphibians may receive more attention.

The phylogeography of Middle Eastern tree frogs in Israel.

Western Palearctic treefrogs of the genus Hyla provide an example of a morphologically and ecologically cryptic group. Up to three distinct Hyla species have been proposed as resident in Israel and this number has consistently been subject to taxonomical debates. Here, we analyzed 16S rRNA and COI gene fragments of 658 individuals sampled at 47 pools in nine regions across Israel and the West Bank in order to resolve the taxonomic status of Hyla frogs. We generated both Bayesian and Maximum Likelihood phylogenies, and constructed time-calibrated trees to provide an evolutionary and historical context of sequence variations. We further applied SAMOVA as well as Monmonier's maximum-difference algorithm to study the genetic structure among populations and to identify potential zones acting as barriers to gene flow across locations. Our results revealed two distinct haplogroups for each gene fragment, with 95% CI divergence times dated from 8.9-17.1 Mya (16S) and 7.1-23.6 Mya (COI), respectively. SAMOVA and barrier analyses partitioned the populations into three groups. Our results highlight that, while there are probably only two Hyla species in Israel, one population of one of the species might qualify as a separate evolutionarily significant unit. Our findings elucidate the taxonomic status of Hyla frogs in Israel and provide the basis for determining appropriate management and conservation priorities.

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.

Multi-omics resources for the Australian southern stuttering frog (Mixophyes australis) reveal assorted antimicrobial peptides.

The number of genome-level resources for non-model species continues to rapidly expand. However, frog species remain underrepresented, with up to 90% of frog genera having no genomic or transcriptomic data. Here, we assemble the first genomic and transcriptomic resources for the recently described southern stuttering frog (Mixophyes australis). The southern stuttering frog is ground-dwelling, inhabiting naturally vegetated riverbanks in south-eastern Australia. Using PacBio HiFi long-read sequencing and Hi-C scaffolding, we generated a high-quality genome assembly, with a scaffold N50 of 369.3 Mb and 95.1% of the genome contained in twelve scaffolds. Using this assembly, we identified the mitochondrial genome, and assembled six tissue-specific transcriptomes. We also bioinformatically characterised novel sequences of two families of antimicrobial peptides (AMPs) in the southern stuttering frog, the cathelicidins and ß-defensins. While traditional peptidomic approaches to peptide discovery have typically identified one or two AMPs in a frog species from skin secretions, our bioinformatic approach discovered 12 cathelicidins and two ß-defensins that were expressed in a range of tissues. We investigated the novelty of the peptides and found diverse predicted activities. Our bioinformatic approach highlights the benefits of multi-omics resources in peptide discovery and contributes valuable genomic resources in an under-represented taxon.