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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Evolutionary and functional conservation of IRF7 in the Tibetan frog Nanorana parkeri.

BACKGROUND: The production of interferons (IFNs) is essential for the control of viral infections, and interferon regulatory factor 7 (IRF7) is considered as a vital regulator for the transcription of type I IFNs. Amphibians appear to possess a highly expanded type I IFN repertoire, consisting of intron-containing genes as observed in fish, and intronless genes as in other higher vertebrates. However, the knowledge on transcriptional regulatory mechanism of these two types of type I IFN genes is rather scarce in amphibians. METHODS AND RESULTS: A IRF7 gene named as Np-IRF7 was identified in Tibetan frog (Nanorana parkeri), and bioinformatic analysis revealed that the predicted protein of Np-IRF7 contains several important structural features known in IRF7. Expression analysis showed that Np-IRF7 gene was widely expressed and rapidly induced by poly(I:C) in different organs/tissues. Interestingly, luciferase reporter assay revealed that intronless IFN promoters were more effectively activated than intron-containing IFN promoter in Np-IRF7-transfected cells. Moreover, the overexpression of Np-IRF7 could induce the expression of ISGs and suppress the replication of FV3 in A6 cells. CONCLUSION: Np-IRF7 is indeed the ortholog of known IRF7, and IRF7 is structurally conserved in different lineages of vertebrates. Np-IRF7 played distinct roles in the activation of intron-containing and intronless type I IFN promoters, thus inducing the expression of interferon-stimulated antiviral effectors and providing a protection against ranavirus infection. The present research thus contributes to a better understanding of regulatory function of IRF7 in the IFN-mediated antiviral response of anuran amphibians.

Multiple contact zones and karyotypic evolution in a neotropical frog species complex.

Previous studies of DNA sequence and karyotypic data have revealed high genetic diversity in the Physalaemus cuvieri - Physalaemus ephippifer species complex-a group of small leptodactylid frogs in South America. To date, seven major genetic lineages have been recognized in this group, with species delimitation tests supporting four to seven of them as valid species. Among these, only P. ephippifer shows heteromorphic sex chromosomes, but the implications of cytogenetic divergence for the evolution of this group are unknown. We analyzed karyotypic, mitochondrial DNA, and 3RAD genomic data to characterize a putative contact zone between P. ephippifer and P. cuvieri Lineage 1, finding evidence for admixture and karyotypic evolution. We also describe preliminary evidence for admixture between two other members of this species complex-Lineage 1 and Lineage 3 of P. cuvieri. Our study sheds new light on evolutionary relationships in the P. cuvieri - P. ephippifer species complex, suggesting an important role of karyotypic divergence in its evolutionary history and underscoring the importance of hybridization as a mechanism of sex chromosome evolution in amphibians.

Unexpected diversity and co-occurrence of phytotelmic frogs (Guibemantis) around Andasibe, one of the most intensively surveyed amphibian hotspots of Madagascar, and descriptions of three new species.

The area around the Malagasy village of Andasibe, which includes Analamazaotra-Mantadia National Park as well as other protected areas, is characterized by very species-rich and well-studied communities of animals and plants, but new species are still regularly discovered. Three species of phytotelmic frogs of the subgenus Pandanusicola in the genus Guibemantis are known from this area, G. flavobrunneus, G. liber, and G. pulcher. Further Pandanusicola frogs from this area have been provisionally assigned to G. bicalcaratus or G. albolineatus, pending detailed taxonomic review. During preliminary exploration of the ecology of these specialized frogs that live and reproduce in the leaf axils of Pandanus screw pines, we noticed the syntopic presence of two differently colored and differently sized Pandanusicola in Andasibe that could not be unambiguously assigned to any known species. A genetic screening revealed that these correspond to yet two further species in the area. Based on our data, seven species of Pandanusicola occur in Andasibe and nearby forests: (1) G. liber, the only non-phytotelmic species of the subgenus in the region; (2) G. flavobrunneus which is the largest species and characterized by a diagnostic yellowish brown dorsal pattern; (3) G. pulcher, characterized by translucent-green color with purplish brown spotting not observed in any other species in the area; (4) G. methueni, a brownish species usually lacking contrasted dorsolateral bands that differs from the other species in the area by emitting a characteristic trill-like advertisement call series (rather than clicks or chirps) and according to our data is widespread along Madagascars east coast; as well as three new species: (5) G. ambakoana sp. nov., a brownish species, typically with contrasted incomplete light dorsolateral bands and with single click-like advertisement calls; (6) G. vakoa sp. nov., a species that is equally brownish but lacks contrasted light dorsolateral bands and that has single click-like advertisement calls of very short duration; and (7) G. rianasoa sp. nov., a species that is smaller sized and has less distinct femoral glands than all the others, and emits a short series of soft chirp-like advertisement calls. All these species are genetically highly distinct, with >5% uncorrected pairwise distances in the mitochondrial 16S rRNA gene, and lack of haplotype sharing in two nuclear-encoded genes. The co-occurrence of seven Pandanusicola frogs in a relatively small geographic area is unprecedented in Madagascar and calls for in-depth studies of a possible differentiation in habitat use and life history.