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.

A new odorous frog species of Odorrana (Amphibia, Anura, Ranidae) from Guizhou Province, China

The frog genus Odorrana is distributed across east and southeastern Asia. Based on morphological differences and molecular phylogenetics, a new species of the genus occurring from Leigong Mountain in Guizhou Province, China is described. Phylogenetic analyses based on DNA sequences of the mitochondrial 12S rRNA, 16S rRNA, and ND2 genes supported the new species as an independent lineage. The uncorrected genetic distances between the 12S rRNA, 16S rRNA, and ND2 genes between the new species and its closest congener were 5.0%, 4.9%, and 16.3%, respectively. The new species is distinguished from its congeners by a combination of the following characters: body size moderate (SVL 39.1-49.4 mm in males, 49.7 mm in female); head width larger than head length; tympanum distinctly visible; small rounded granules scattered all over dorsal body and limbs; dorsolateral folds absent; heels overlapping when thighs are positioned at right angles to the body; tibiotarsal articulation reaching the level between eye to nostril when leg stretched forward; vocal sacs absent in male and nuptial pads present on the base of finger I.

Characterization and Comparison of the Two Mitochondrial Genomes in the Genus Rana

The mitochondrial genome (mitogenome) possesses several invaluable attributes, including limited recombination, maternal inheritance, a fast evolutionary rate, compact size, and relatively conserved gene arrangement, all of which make it particularly useful for applications in phylogenetic reconstruction, population genetics, and evolutionary research. In this study, we aimed to determine the complete mitogenomes of two morphologically similar Rana species (Rana hanluica and Rana longicrus) using next-generation sequencing. The entire circular mitogenome was successfully identified, with a length of 19,395 bp for R. hanluica and 17,833 bp for R. longicrus. The mitogenomes of both species contained 37 genes, including 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, and one control region; mitogenome size varied predominantly with the length of the control region. The two synonymous codon usages in 13 PCGs showed that T and A were used more frequently than G and C. The ratios of non-synonymous to synonymous substitutions of all 13 PCGs were <1 in the Rana species, indicating that the PCGs were under purifying selection. Finally, phylogenetic relationship analyses suggested that R. hanluica and R. longicrus were classified in the R. japonica group. Our study provides valuable reference material for the taxonomy of the genus Rana.

A synoptic review of the Amphibians of Iran: bibliography, taxonomy, synonymy, distribution, conservation status, and identification key to the eggs, larvae, and adults

This study provides an illustrated account, a comprehensive update of the systematics, and a bibliography of the 15 species of anurans in five families, eight genera; and of the six species of urodeles in two families, four genera in Iran. Bufonidae, with eight species, is the most diverse family; Salamandridae has five species and Ranidae has four species. This study also presents updated identification keys for the eggs, larvae, and metamorphosed amphibians of Iran. We designated specimen NMW 19855.1 as neotype of Pelophylax persicus (Schneider, 1799) comb. nov.. Along with distribution maps obtained from all the reliable localities and museum specimens known at this time, the modelled habitat of species, and for the first time, the National Red List of amphibians based on the IUCN red list categories and criteria. Based on our evaluation we propose to categorize Bufo eichwaldi, Paradactylodon persicus, Neurergus derjugini, and N. kaiseri as Vulnerable at National Red List, and to move Bufotes (Calliopersa) luristanicus, B. (C.) surdus, Firouzophrynus olivaceus, and Rana pseudodalmatina from the category of Least Concern (LC) to Near Threatened (NT). The National Red List of amphibians that we propose has significant implications for endangered species management and conservation. Forty-one percent of amphibian species in Iran are endemic to the country, and more than forty percent of the Iranian amphibians are at risk of extinction. Zagros Mountain forest and Hyrcaniain forests have more than 80% (i.e. 18 species) of the diversity of Iranian amphibians. A considerable amount of scientific literature published on Iranian amphibians in Persian language is not easily accessible to researchers outside Iran. This monograph attempts to remedy the situation and provides broader access to international herpetology. We recognize that taxonomy is always in a state of flux, and the names and synonymies used here reflect our current view.

A global survey of tapeworms (Cestoda: Proteocephalidae) of 'true' frogs (Amphibia: Ranidae), including a tabulated list of all proteocephalids parasitising amphibians

Proteocephalid tapeworms of frogs of the family Ranidae ('true' frogs) are reviewed with emphasis on their species diversity, host specificity and geographical distribution. New molecular data (nuclear lsrDNA and mitochondrial COI sequences) are presented for tapeworms of four species of ranid frogs in North America, including the poorly known Ophiotaenia saphena Osler, 1931 of Rana clamitans Latreille and R. catesbeiana (Shaw), which is redescribed using new material from Arkansas, USA. Tapeworms of R. sphenocephala (Cope) and R. pipiens Schreber, the latter previously identified as O. saphena, represent another, putative new species, but are not formally described due to insufficient available material. Proteocephalus papuensis Bursey, Goldberg et Kraus, 2008 from Sylvirana supragrisea (Menzies) is transferred to Ophiotaenia La Rue, 1911 as a new combination. After a critical review of the literature, only nine nominal species of Ophiotaenia are recognised as valid, which is in contrast to the large number of ranid frogs (> 440 spp.). The reasons for this striking disparity are briefly discussed, and a key based on morphology is presented for the identification of all species of Ophiotaenia from the Ranidae. Molecular data are available for only two taxa from North America that form a monophyletic group. The relationships among tapeworms of ranid frogs occurring in other zoogeographical regions are not yet known. The taxonomic status of Batrachotaenia Rudin, 1917, which was erected to accommodate proteocephalids from amphibians, is also discussed. To facilitate future studies, a tabulated summary of all 32 species of proteocephalids belonging to three genera reported from amphibians (frogs and salamanders) is presented, with information on their hosts, distribution, and taxonomically important characters, including key measurements.

2b or not 2b? 2bRAD is an effective alternative to ddRAD for phylogenomics

Restriction-site-associated DNA sequencing (RADseq) has become an accessible way to obtain genome-wide data in the form of single-nucleotide polymorphisms (SNPs) for phylogenetic inference. Nonetheless, how differences in RADseq methods influence phylogenetic estimation is poorly understood because most comparisons have largely relied on conceptual predictions rather than empirical tests. We examine how differences in ddRAD and 2bRAD data influence phylogenetic estimation in two non-model frog groups. We compare the impact of method choice on phylogenetic information, missing data, and allelic dropout, considering different sequencing depths. Given that researchers must balance input (funding, time) with output (amount and quality of data), we also provide comparisons of laboratory effort, computational time, monetary costs, and the repeatability of library preparation and sequencing. Both 2bRAD and ddRAD methods estimated well-supported trees, even at low sequencing depths, and had comparable amounts of missing data, patterns of allelic dropout, and phylogenetic signal. Compared to ddRAD, 2bRAD produced more repeatable datasets, had simpler laboratory protocols, and had an overall faster bioinformatics assembly. However, many fewer parsimony-informative sites per SNP were obtained from 2bRAD data when using native pipelines, highlighting a need for further investigation into the effects of each pipeline on resulting datasets. Our study underscores the importance of comparing RADseq methods, such as expected results and theoretical performance using empirical datasets, before undertaking costly experiments.

First evidence of ranavirus in native and invasive amphibians in Colombia

Ranaviruses can cause mass mortality events in amphibians, thereby becoming a threat to populations that are already facing dramatic declines. Ranaviruses affect all life stages and persist in multiple amphibian hosts. The detrimental effects of ranavirus infections to amphibian populations have already been observed in the UK and in North America. In Central and South America, the virus has been reported in several countries, but the presence of the genus Ranavirus (Rv) in Colombia is unknown. To help fill this knowledge gap, we surveyed for Rv in 60 species of frogs (including one invasive species) in Colombia. We also tested for co-infection with Batrachochytrium dendrobatidis (Bd) in a subset of individuals. For Rv, we sampled 274 vouchered liver tissue samples collected between 2014 and 2019 from 41 localities covering lowlands to mountaintop páramo habitat across the country. Using quantitative polymerase chain reaction (qPCR) and end-point PCR, we detected Rv in 14 individuals from 8 localities, representing 6 species, including 5 native frogs of the genera Osornophryne, Pristimantis and Leptodactylus, and the invasive American bullfrog Rana catesbeiana. Bd was detected in 7 of 140 individuals, with 1 co-infection of Rv and Bd in an R. catesbeiana specimen collected in 2018. This constitutes the first report of ranavirus in Colombia and should set off alarms about this new emerging threat to amphibian populations in the country. Our findings provide some preliminary clues about how and when Rv may have spread and contribute to understanding how the pathogen is distributed globally.

Population Subdivision in the Gopher Frog (Rana capito) across the Fragmented Longleaf Pine-Wiregrass Savanna of the Southeastern USA

Delineating genetically distinct population segments of threatened species and quantifying population connectivity are important steps in developing effective conservation and management strategies aimed at preventing extinction. The gopher frog (Rana capito) is a xeric-adapted, pond-breeding species endemic to the Gulf and Atlantic coastal plains of the southeastern United States. This species has experienced extensive habitat loss and fragmentation in the formerly widespread longleaf pine-wiregrass savanna where it lives, resulting in individual abundance declines and population extinctions throughout its range. We used individual-based clustering methods along with Bayesian inference of historical migration based on almost 1500 multilocus microsatellite genotypes to examine genetic structure in this taxon. Clustering analyses identified panhandle and peninsular populations in Florida as distinct genetic clusters separated by the Aucilla River, consistent with the division between the Coastal Plain and peninsular mitochondrial lineages, respectively. Analysis of historical migration indicated an east–west population divergence event followed by immigration to the east. Together, our results indicate that the genetically distinct Coastal Plain and peninsular Florida lineages should be considered separately for conservation and management purposes.

Comparative transcriptomes of nine tissues for the Heilongjiang brown frog (Rana amurensis)

The Heilongjiang brown frog (Rana amurensis) is widely used in traditional Chinese medicine. In particular, the oviduct and skin have been developed into various health products. However, limited numbers of complete genomes of amphibian species have been reported, excluding the Heilongjiang brown frog. Here, the transcriptomes of 45 samples from the liver, spleen, heart, ovaries, thigh muscles, skin, oviduct, stomach and intestine of five Heilongjiang brown frog were reassembled and analyzed. A total of 1,085,532 unigenes with an average length of 676.6 bp and N50 of 722 bp were obtained. Comparative transcriptomics of different tissues detected tissue-specific expression. There were 3248 differentially expressed genes (DEGs) in the ovary, and the number of unique DEGs between the ovary and spleen was the largest. The results of DEGs enrichment showed there were many pathways and items related to protein synthesis and metabolism in the oviduct. The DEGs of the skin were enriched with many bacterial defense items, indicating that there were a large number of antimicrobial peptides in the skin. Thus, these were suitable as biological sources for the development and extraction of antimicrobial peptides. Through the assembly of transcriptome sequencing data and functional annotation of the Heilongjiang brown frog genome, this study provides reference materials for further exploring and utilizing functional gene resources of frogs and lays a foundation for medical research and the development of new products.

Gene expression of secretory proteins in the nuptial pads of three Lithobates species (Anura: Ranidae)

Many amphibian species possess male-specific secretory breeding glands. In anurans, such specialized glands are often present as nuptial pads at the first digit of the hand, which are pressed against the female during amplexus. However, the role of nuptial pad secretions remains largely unknown. Here, we investigate the inner morphology as well as the RNA expression patterns of the nuptial pads of several Central American Lithobates species (Ranidae). As shown for the breeding glands of other amphibian species, the Lithobates nuptial pads are composed of large specialized mucus glands, excreting proteinaceous content to the surface. Whole-transcriptome sequencing revealed that for one of the species the most highly expressed transcripts encoding secretory proteins in the nuptial pads are sodefrin precursor-like factor (SPF) proteins, known to have a pheromone function in multiple amphibian species. The other two species, however, lack high expression of SPF transcripts but express other secretory proteins, whose roles in amphibian breeding glands remain unknown. Several of these proteins are highly expressed in all samples, suggesting a direct role during courtship or for the maintenance/development of the glands.

A regression-based QSAR-model to predict acute toxicity of aromatic chemicals in tadpoles of the Japanese brown frog (Rana japonica): Calibration, validation, and future developments to support risk assessment of chemicals in amphibians

Amphibian populations are undergoing a global decline worldwide. Such decline has been attributed to their unique physiology, ecology, and exposure to multiple stressors including chemicals, temperature, and biological hazards such as fungi of the Batrachochytrium genus, viruses such as Ranavirus, and habitat reduction. There are limited toxicity data for chemicals available for amphibians and few quantitative structure-activity relationship (QSAR) models have been developed and are publicly available. Such QSARs provide important tools to assess the toxicity of chemicals particularly in a data poor context. QSARs provide important tools to assess the toxicity of chemicals particularly when no toxicological data are available. This manuscript provides a description and validation of a regression-based QSAR model to predict, in a quantitative manner, acute lethal toxicity of aromatic chemicals in tadpoles of the Japanese brown frog (Rana japonica). QSAR models for acute median lethal molar concentrations (LC50-12 h) of waterborne chemicals using the Monte Carlo method were developed. The statistical characteristics of the QSARs were described as average values obtained from five random distributions into training and validation sets. Predictions from the model gave satisfactory results for the overall training set (R² = 0.72 and RMSE = 0.33) and were even more robust for the validation set (R² = 0.96 and RMSE = 0.11). Further development of QSAR models in amphibians, particularly for other life stages and species, are discussed.

Transcriptomic evidence for visual adaptation during the aquatic to terrestrial metamorphosis in leopard frogs

BACKGROUND

Differences in morphology, ecology, and behavior through ontogeny can result in opposing selective pressures at different life stages. Most animals, however, transition through two or more distinct phenotypic phases, which is hypothesized to allow each life stage to adapt more freely to its ecological niche. How this applies to sensory systems, and in particular how sensory systems adapt across life stages at the molecular level, is not well understood. Here, we used whole-eye transcriptomes to investigate differences in gene expression between tadpole and juvenile southern leopard frogs (Lithobates sphenocephalus), which rely on vision in aquatic and terrestrial light environments, respectively. Because visual physiology changes with light levels, we also tested the effect of light and dark exposure.

RESULTS

We found 42% of genes were differentially expressed in the eyes of tadpoles versus juveniles and 5% for light/dark exposure. Analyses targeting a curated subset of visual genes revealed significant differential expression of genes that control aspects of visual function and development, including spectral sensitivity and lens composition. Finally, microspectrophotometry of photoreceptors confirmed shifts in spectral sensitivity predicted by the expression results, consistent with adaptation to distinct light environments.

CONCLUSIONS

Overall, we identified extensive expression-level differences in the eyes of tadpoles and juveniles related to observed morphological and physiological changes through metamorphosis and corresponding adaptive shifts to improve vision in the distinct aquatic and terrestrial light environments these frogs inhabit during their life cycle. More broadly, these results suggest that decoupling of gene expression can mediate the opposing selection pressures experienced by organisms with complex life cycles that inhabit different environmental conditions throughout ontogeny.

Comparative Mitogenomics of True Frogs (Ranidae, Anura), and Its Implications for the Phylogeny and Evolutionary History of Rana

Simple Summary

The true frogs of the genus Rana are a complex and diverse group. Many new species have been discovered with the help of molecular markers and morphological traits. However, the evolutionary history in Rana were not well understood. In this study, we sequenced and annotated the complete mitochondrial genome of R. longicrus and R. zhenhaiensis. In 13 protein codon genes, the COI was the most conserved, and ATP8 had a fast rate of evolution. The Ka/Ks ratio analysis among Rana indicated the protein-coding genes were suffering purify selection. There were three kinds of gene arrangement patterns found. This study provides mitochondrial genetic information, improving our understanding of mitogenomic structure and evolution, and recognizes the phylogenetic relationship and taxonomy among Rana.

Abstract

The true frogs of the genus Rana are a complex and diverse group, containing approximately 60 species with wide distribution across Eurasia and the Americas. Recently, many new species have been discovered with the help of molecular markers and morphological traits. However, the evolutionary history in Rana was not well understood and might be limited by the absence of mitogenome information. In this study, we sequenced and annotated the complete mitochondrial genome of R. longicrus and R. zhenhaiensis, containing 22 tRNAs, 13 protein-coding genes, two ribosomal RNAs, and a non-coding region, with 17,502 bp and 18,006 bp in length, respectively. In 13 protein codon genes, the COI was the most conserved, and ATP8 had a fast rate of evolution. The Ka/Ks ratio analysis among Rana indicated the protein-coding genes were suffering purify selection. There were three kinds of gene arrangement patterns found. The mitochondrial gene arrangement was not related to species diversification, and several independent shifts happened in evolutionary history. Climate fluctuation and environmental change may have played an essential role in species diversification in Rana. This study provides mitochondrial genetic information, improving our understanding of mitogenomic structure and evolution, and recognizes the phylogenetic relationship and taxonomy among Rana.

The effects of intra- and interspecific competitions on personality and individual plasticity in two sympatric brown frogs

Abstract

We studied how individuals modify their behavior in response to inter- and intraspecific competitors and how these changes affected the pattern of variation between populations and species. As study models, we used tadpoles of two brown frogs, Rana latastei and R. dalmatina. Since R. latastei is always sympatric to R. dalmatina, whereas R. dalmatina is sympatric to R. latastei only in the periphery of its range, we predicted a stronger response to heterospecifics in R. latastei than in R. dalmatina and, within each species, in syntopic than in allotopic populations. To test these predictions, we raised tadpoles, from either syntopic or allotopic populations, in either syntopy or allotopy and repeatedly tested them in open field trials in the presence of a caged conspecific, a caged heterospecific, or an empty cage. As predicted, we found that, on average, R. latastei tadpoles modified their behavior across treatments more than R. dalmatina tadpoles and individuals from the syntopic population changed more than their conspecifics from the allotopic population. In both species, the pattern of variation at the individual level mirrored that at the population and species levels providing no evidence for an individual-by-environment interaction ($$I\times E$$ I × E ). Besides these differences, however, individuals of the two species also showed unpredicted and context-independent behavioral differences, suggesting that there might be more to interspecific behavioral variation than the effect of selection by heterospecific competitors.

Significance statement

Does the distribution range of a species influence the evolution of plastic behaviors to heterospecific competitors? And how do differences in plasticity affect animal personality? To answer these questions, we raised tadpoles of two brown frog species, Rana dalmatina and R. latastei, and studied how the amount and the type of their swimming varied with the presence of the other species. R. latastei, whose small distribution range fully overlaps with that of R. dalmatina, plastically responds to it, whereas R. dalmatina, which is sympatric to R. latastei only in the periphery of its broader range, does not. These interspecific differences mirrored those among individuals: tadpoles of both species show repeatable behaviors, but only those of R. latastei plastically changed their behavior with the presence of the other species; however, neither R. latastei nor R. dalmatina show among-individual variation in plasticity.

Host species is linked to pathogen genotype for the amphibian chytrid fungus (Batrachochytrium dendrobatidis)

Host-pathogen specificity can arise from certain selective environments mediated by both the host and pathogen. Therefore, understanding the degree to which host species identity is correlated with pathogen genotype can help reveal historical host-pathogen dynamics. One animal disease of particular concern is chytridiomycosis, typically caused by the global panzootic lineage of the amphibian chytrid fungus (Batrachochytrium dendrobatidis, Bd), termed Bd-GPL. This pathogen lineage has caused devastating declines in amphibian communities around the world. However, the site of origin for the common ancestor of modern Bd-GPL and the fine-scale transmission dynamics of this lineage have remained a mystery. This is especially the case in North America where Bd-GPL is widespread, but disease outbreaks occur sporadically. Herein, we use Bd genetic data collected throughout the United States from amphibian skin swabs and cultured isolate samples to investigate Bd genetic patterns. We highlight two case studies in Pennsylvania and Nevada where Bd-GPL genotypes are strongly correlated with host species identity. Specifically, in some localities bullfrogs (Rana catesbeiana) are infected with Bd-GPL lineages that are distinct from those infecting other sympatric amphibian species. Overall, we reveal a previously unknown association of Bd genotype with host species and identify the eastern United States as a Bd diversity hotspot and potential site of origin for Bd-GPL.

Competition between the tadpoles of Japanese toads versus frogs

Competition within and among species can play a key role in structuring the assemblages of anuran tadpoles. Previous studies have reported that tadpoles of the invasive cane toad (Rhinella marina) are more strongly disadvantaged by the presence of native frog tadpoles than by the same number of conspecific toad tadpoles. That effect might arise from a lack of coevolution of the invasive toad with its competitors; and/or from a generalized superiority of frog tadpoles over toad tadpoles. To clarify those possibilities, we conducted experimental trials using the larvae of a native rather than invasive toad (Bufo japonicus formosus in Japan) exposed to larvae of native anurans (the sympatric frogs Rana japonica and Rana ornativentris and the parapatric toad Bufo japonicus japonicus). In intraspecific competition trials, higher densities of B. j. formosus prolonged the larval period and reduced size at metamorphosis, but did not affect survival. In interspecific competition trials, the effects of the other anuran species on B. j. formosus were similar to the effects of the same number of conspecific larvae. This similarity in impact of interspecific versus intraspecific competition argues against any overall competitive superiority of frog larvae over toad larvae. Instead, the vulnerability of larval cane toads to frog tadpoles may result from a lack of coevolutionary history.

An integrative taxonomic revision of slug-eating snakes (Squamata: Pareidae: Pareineae) reveals unprecedented diversity in Indochina

Slug-eating snakes of the subfamily Pareinae are an insufficiently studied group of snakes specialized in feeding on terrestrial mollusks. Currently Pareinae encompass three genera with 34 species distributed across the Oriental biogeographic region. Despite the recent significant progress in understanding of Pareinae diversity, the subfamily remains taxonomically challenging. Here we present an updated phylogeny of the subfamily with a comprehensive taxon sampling including 30 currently recognized Pareinae species and several previously unknown candidate species and lineages. Phylogenetic analyses of mtDNA and nuDNA data supported the monophyly of the three genera Asthenodipsas, Aplopeltura, and Pareas. Within both Asthenodipsas and Pareas our analyses recovered deep differentiation with each genus being represented by two morphologically diagnosable clades, which we treat as subgenera. We further apply an integrative taxonomic approach, including analyses of molecular and morphological data, along with examination of available type materials, to address the longstanding taxonomic questions of the subgenus Pareas, and reveal the high level of hidden diversity of these snakes in Indochina. We restrict the distribution of P. carinatus to southern Southeast Asia, and recognize two subspecies within it, including one new subspecies proposed for the populations from Thailand and Myanmar. We further revalidate P. berdmorei, synonymize P. menglaensis with P. berdmorei, and recognize three subspecies within this taxon, including the new subspecies erected for the populations from Laos and Vietnam. Furthermore, we describe two new species of Pareas from Vietnam: one belonging to the P. carinatus group from southern Vietnam, and a new member of the P. nuchalis group from the central Vietnam. We provide new data on P. temporalis, and report on a significant range extension for P. nuchalis. Our phylogeny, along with molecular clock and ancestral area analyses, reveal a complex diversification pattern of Pareinae involving a high degree of sympatry of widespread and endemic species. Our analyses support the "upstream" colonization hypothesis and, thus, the Pareinae appears to have originated in Sundaland during the middle Eocene and then colonized mainland Asia in early Oligocene. Sundaland and Eastern Indochina appear to have played the key roles as the centers of Pareinae diversification. Our results reveal that both vicariance and dispersal are responsible for current distribution patterns of Pareinae, with tectonic movements, orogeny and paleoclimatic shifts being the probable drivers of diversification. Our study brings the total number of Pareidae species to 41 and further highlights the importance of comprehensive taxonomic revisions not only for the better understanding of biodiversity and its evolution, but also for the elaboration of adequate conservation actions.

Diversity, distribution and molecular species delimitation in frogs and toads from the Eastern Palaearctic

Biodiversity analyses can greatly benefit from coherent species delimitation schemes and up-to-date distribution data. In this article, we have made the daring attempt to delimit and map described and undescribed lineages of anuran amphibians in the Eastern Palaearctic (EP) region in its broad sense. Through a literature review, we have evaluated the species status considering reproductive isolation and genetic divergence, combined with an extensive occurrence dataset (nearly 85k localities). Altogether 274 native species from 46 genera and ten families were retrieved, plus eight additional species introduced from other realms. Independent hotspots of species richness were concentrated in southern Tibet (Medog County), the circum-Sichuan Basin region, Taiwan, the Korean Peninsula and the main Japanese islands. Phylogeographic breaks responsible for recent in situ speciation events were shared around the Sichuan Mountains, across Honshu and between the Ryukyu Island groups, but not across shallow water bodies like the Yellow Sea and the Taiwan Strait. Anuran compositions suggested to restrict the zoogeographical limits of the EP to East Asia. In a rapidly evolving field, our study provides a checkpoint to appreciate patterns of species diversity in the EP under a single, spatially explicit, species delimitation framework that integrates phylogeographic data in taxonomic research.

Molecular phylogeny and historical biogeography of the cave fish genus Sinocyclocheilus (Cypriniformes: Cyprinidae) in southwest China

Modern accumulations of genetic data offer unprecedented opportunities for understanding the systematic classification and origins of specific groups of organisms. The genus Sinocyclocheilus is among the most cave abundant genera in Cyprinidae, with 76 recognized species, belonging to four species groups. Recent phylogenetic studies have shown that the classification of species groups within the genus Sinocyclocheilus remains controversial. In this study, we constructed a sequence supermatrix of 26 species from four species groups of the genus Sinocyclocheilus using the mitochondrial genome to reveal phylogenetic relationships, historical biogeography and patterns of species diversification in the genus Sinocyclocheilus. Phylogenetic analysis strongly supports the monophyletic groups of the three species groups (S. jii, S. cyphotergous, and S. tingi groups) except the S. angularis group. Phylogenetic analysis showed that S. anshuiensis and S. microphthalmus, which were recognized as numbers of S. angularis group, formed a strongly supported independent clade. Therefore, we propose a new species group, the S. microphthalmus group, which contains S. anshuiensis and S. microphthalmus. Biogeographic reconstruction suggests that the living Sinocyclocheilus may have originated in north-central Guangxi at the late Eocene and dispersed outward after a vicariance at 32.31 Million years ago (Ma). Early diversification is focused on the late Oligocene (ca. 25 Ma), which is related to the second uplift of the Qinghai-Tibetan Plateau and the lateral extrusion of the Indochina at the Oligocene/Miocene boundary. Our results suggest that two uplifts of the Qinghai-Tibetan Plateau and climate change in the Miocene may have influenced the diversification of the Sinocyclocheilus lineage. This article is protected by copyright. All rights reserved.

An annotated checklist of the amphibians and reptiles of North Padre Island, Texas, USA, with comparisons to adjacent barrier island and mainland herpetofauna

Padre Island is the world's longest barrier island and includes the longest stretch of undeveloped barrier island in the world. Largely due to harsh environmental conditions and difficult access, only cursory and incomplete checklists and subjective estimates of abundance have been produced. The results of an inventory of amphibians and reptiles of North Padre Island conducted 2002-2020, including the results of extensive field surveys conducted 2002-2003, are reported herein. Natural history museum and iNaturalist records are summarized and compared among North and South Padre and Mustang islands and the mainland portion of the seven counties in which the islands occur. The conservation status of rare species and extirpation of others is noted. The morphology and taxonomic status of some unique occurrences are discussed. Eleven species of amphibians and 39 species of reptiles presently occur or have occurred naturally or as introduced or accidental species on North Padre Island. Twelve species of amphibians and 50 species of reptiles occur or have occurred on North Padre, South Padre, and Mustang islands. Thirty-one species of amphibians and 93 species of reptiles have been reported from the seven counties in which the islands occur.

An empirical test of the relationship between the bootstrap and likelihood ratio support in maximum likelihood phylogenetic analysis

In maximum likelihood (ML), the support for a clade can be calculated directly as the likelihood ratio (LR) or log-likelihood difference (S, LLD) of the best trees with and without the clade of interest. However, bootstrap (BS) clade frequencies are more pervasive in ML phylogenetics and are almost universally interpreted as measuring support. In addition to theoretical arguments against that interpretation, BS has several undesirable attributes for a support measure. For example, it does not vary in proportion to optimality or identify clades that are rejected by the evidence and can be overestimated due to missing data. Nevertheless, if BS is a reliable predictor of S, then it might be an efficient indirect method of measuring support-an attractive possibility, given the speed of many BS implementations. To assess the relationship between S and BS, we analyzed 106 empirical datasets retrieved from TreeBASE. Also, to evaluate the degree to which S and BS are affected by the number of replicates during suboptimal tree searches for S and pseudoreplicates during BS estimation, we randomly selected 5 of the 106 datasets and analyzed them using variable numbers of replicates and pseudoreplicates, respectively. The correlation between S and BS was extremely weak in the datasets we analyzed. Increasing the number of replicates during tree search decreased the estimated values of S for most clades, but the magnitude of change was small. In contrast, although increasing pseudoreplicates affected BS values for only approximately 40% of clades, values both increased and decreased, and they did so at much greater magnitudes. Increasing replicates/pseudoreplicates affected the rank order of clades in each tree for both S and BS. Our findings show decisively that BS is not an efficient indirect method of measuring support and suggest that even quite superficial searches to calculate S provide better estimates of support.

A new species of the genus Rana sensu lato Linnaeus, 1758 (Anura, Ranidae) from Wuyi Mountain, Fujian Province, China

A new species of the frog genus Rana sensu lato from Wuyi Mountain, Fujian Province, China is described. Molecular phylogenetic analyses clustered the new species into the R.johnsi group and indicated that it was genetically divergent from its closely related species. The new species could be distinguished from its congeners by a combination of the following characters: body size medium, SVL 41.4–45.6 mm (42.9 ± 1.9 mm, n = 4) in adult males and 47.6–50.3 mm (n = 2) in adult females; adult male with a pair of internal subgular vocal sacs; lateroventral grooves present on tip of toes; webbing on fourth toes reaching the tip of toe; transverse skin ridges distinctly present on the dorsal surface of thigh and tibia, the number large (mean 26.5 ± 2.7, range 22–29, n = 6); breeding males possess creamy white nuptial pad with tiny velvety spines on the dorsal surface of the first finger, divided into three parts.

GENETIC DIVERSITY OF CYSTODISCUS SPECIES IN AMPHIBIANS IN THE SOUTHERN UNITED STATES

Myxosporean species in the genus Cystodiscus are parasites of amphibians and have been reported from several continents. Typically used for the identification of myxozoans, the spores produced by these species are similar to one another, possessing 2 polar capsules and being ovoid. The number of transverse depressions on the spore can be useful for delineating species, but these can sometimes be difficult to distinguish. In North America, Cystodiscus serotinus and Cystodiscus melleni have been described, and for C. serotinus in particular, numerous reports and a wide range of hosts have been associated with this species. Given the challenges of identifying some of these species, we questioned whether all encounters of Cystodiscus species can be attributed to these 2 described species, or if there may be additional undescribed species or cryptic species. Over 7 yr, 383 amphibians representing 13 species of toads, frogs, and salamanders were collected from sites in Oklahoma and Arkansas. Cystodiscus infections were found in 56 individuals (14.6%). Tissues from these infected individuals were preserved in alcohol for genetic analysis. The small subunit (SSU) and large subunit (LSU) ribosomal RNA genes were partially sequenced and analyzed phylogenetically. Nine distinct SSU sequence types and 7 distinct LSU sequence types were identified. Phylogenetically, sequence types were attributable to C. serotinus, C. melleni, Cystodiscus axonis, and an undescribed species. For the previously described species, there were multiple SSU sequence types: 4 for C. serotinus and 2 for both C. melleni and C. axonis. Phylogenetic patterns were similar for the LSU sequence analysis using a shorter sequence than the SSU, and we propose that the LSU is useful for initial barcoding of Cystodiscus species in any future surveys. In our qualitative assessment of sequence types compared to geography and host species, SSU types C1 and C2 (C. axonis) were only found in Union County, Arkansas, and McCurtain County, Oklahoma, respectively. Also, salamanders were only infected with SSU types B or D (C. melleni), and type B was only found in salamanders. Our finding of C. axonis in North America is notable because this species was described in Australia and is associated with host pathology. Our work reveals that there are cryptic species of Cystodiscus in the United States, one of which may be a pathogen, highlighting the importance of genetic analysis for future surveys of these species.

Clothianidin alters leukocyte profiles and elevates measures of oxidative stress in tadpoles of the amphibian, Rana pipiens

Neonicotinoid pesticide use is widespread and highly debated, as evidenced by recent attention received from the public, academics and pesticide regulatory agencies. However, relatively little is known about the physiological effects of neonicotinoid insecticides on aquatic vertebrates. Amphibians (larval stages in particular) are excellent vertebrate bioindicators in aquatic systems due to their risk of exposure and sensitivity to environmental stressors. Previous work with wood frog (Rana sylvatica) tadpoles exposed to formulated products containing thiamethoxam or clothianidin in outdoor mesocosms found significant shifts in leukocyte profiles, suggesting the tadpoles were physiologically stressed. The main objective of the present study was to characterize this stress response further using complementary measures of stress after exposure to clothianidin on northern leopard frogs (Rana pipiens) during their aquatic larval stages. Laboratory static-renewal exposures were conducted over eight weeks with the technical product clothianidin at 0, 0.23, 1, 10 and 100 μg/L, and diquat dibromide at 532 μg/L was used as a positive control. We assessed tadpole leukocyte profiles and measures of oxidative stress as these sub-lethal alterations could affect amphibian fitness. We found changes in several types of leukocytes at 1 and 10 μg/L, suggesting that these tadpoles exhibited signs of mild physiological stress. Clothianidin also induced an oxidative stress response at 0.23, 1 and 100 μg/L. However, we found no differences in survival, growth, development time or hepatosomatic index in frogs exposed to clothianidin. Our study indicates that tadpoles chronically exposed to clothianidin have increased stress responses, but in the absence of concentration-response relationships and effects on whole-organism endpoints, the implications on the overall health and fitness of these changes are unclear.

A new species of Gorgoderina (Digenea: Gorgoderidae) from Rana berlandieri in Los Tuxtlas tropical rainforest, Mexico, with the elucidation of its life cycle

Species in the genus Gorgoderina Looss, 1902 are parasites of the urinary bladder of amphibians and include around 50 species described globally. Molecular data on species of the genus are scarce, as is the information of their life-cycle patterns. During a survey on the genetic characterization of the frog trematodes in the tropical rain forest of Los Tuxtlas, in the Gulf of Mexico slope of Mexico, specimens of two morphotypes of Gorgoderina were sampled from the Rio Grande leopard frog, Rana berlandieri. One of them represented an undescribed species which is described herein as Gorgoderina rosamondae n. sp., whereas the other one was morphologically very similar to an apparently widely distributed North American species, G. attenuata, which has been previously reported in the same geographical area. Specimens of both morphotypes were sequenced for two nuclear and one mitochondrial genes. Phylogenetic trees corroborated the distinction of the new species, and data on the internal transcribed spacer 2 revealed genetic differences between G. attenuata sequenced from frogs in USA and specimens of Gorgoderina sp. from Los Tuxtlas, indicating the possibility that they also represent an undescribed species. COI sequences showed high genetic divergence values between the new species and Gorgoderina sp. from Los Tuxtlas (8.63-9.99%). Additionally, COI sequences of the larval forms (sporocyst, cercariae and metacercariae) sampled in the same locality from their first and second intermediate hosts (Pisidium sp. and Agriogomphus tumens, respectively) showed conspecificity, and the 3 host life-cycle of the new species was elucidated.

Spatiotemporal adaptive evolution of an MHC immune gene in a frog-fungus disease system

Genetic diversity of major histocompatibility complex (MHC) genes is linked to reduced pathogen susceptibility in amphibians, but few studies also examine broad spatial and temporal patterns of MHC and neutral genetic diversity. Here, we characterized range-wide MHC diversity in the Northern leopard frog, Rana pipiens, a species found throughout North America that is experiencing disease-related declines. We used previously sequenced neutral markers (mitochondrial DNA and microsatellites), sequenced an expressed MHC class IIß gene fragment, and measured infection prevalence and intensity of the global fungal pathogen Batrachochytrium dendrobatidis (Bd) across 14 populations. Four populations were sampled across two decades, enabling temporal comparisons of selection and demography. We recovered 37 unique MHC alleles, including 17 that were shared across populations. Phylogenetic and population genetic patterns between MHC and neutral markers were incongruent, and five MHC codon positions associated with peptide binding were under positive selection. MHC heterozygosity, but not neutral marker heterozygosity, was a significant factor explaining spatial patterns of Bd prevalence, whereas only environmental variables predicted Bd intensity. MHC allelic richness (AR) decreased significantly over time but microsatellite-based AR did not, highlighting a loss of functional immunogenetic diversity that may be associated with Bd selective pressures. MHC supertype 4 was significantly associated with an elevated risk of Bd infection, whereas one supertype 2 allele was associated with a nearly significant reduced risk of Bd. Taken together, these results provide evidence that positive selection contributes to MHC class IIß evolution in R. pipiens and suggest that functional MHC differences across populations may contribute to disease adaptation.

Stability in the Balkans: phylogeography of the endemic Greek stream frog, Rana graeca

We still have little knowledge concerning the phylogeography of amphibians and reptiles from the Balkan Peninsula compared with the other two Mediterranean peninsulas. This raises concerns for endemic taxa from these peninsulas, because it might interfere with further conservation efforts. Here we focus on the endemic Greek stream frog (Rana graeca) and reconstruct its biogeography and evolutionary history. Using four genetic markers (Cytb, 16S, COI and BDNF) in &gt; 350 sequences covering the whole distribution range, we conducted phylogenetic, demographic and ecological niche analyses, which revealed the phylogeography of this species. Surprisingly, this examination of R. graeca reveals a very shallow level of intraspecific genetic variability through the Balkans, with two main, statistically supported lineages having a partly sympatric distribution. The most variable marker was Cytb, which showed 19 haplotypes in 123 analysed sequences in the whole species distribution area. Here presented genetic data, together with the environmental niche projection and demographic analyses suggest that R. graeca was probably affected only marginally by climatic oscillations, with the Hellenides as the most suitable area for the occurrence of the species in different geological periods. This is consistent with the observed genetic diversity, which is mostly related to these mountains. Although the species shows a certain level of phenotypic variability and ecological preferences, this might be related to species plasticity affected by the micro-climatic conditions in small areas, which merits further research. Comparing phylogeography of other amphibian and reptile species in the Balkans, we showed that the observed pattern represents a new view on the phylogeography of the Balkan herpetofauna.

Impact of the Miocene orogenesis on Kaloula spp. radiation and implication of local refugia on genetic diversification

The phylogeography of the Kaloula genus in East Asia is still poorly understood. One of the difficulties is the absence of fossils to corroborate molecular dating estimates. Here, we examined the mitochondrial structure of Kaloula spp. in East Asia and focused on the impact of glaciations on the northernmost species: Kaloula borealis. We determined the phylogenetic relationships, molecular dating, and genetic connectivity assessments within the genus from 1211 bp of concatenated mitochondrial 12S and 16S. The relaxed clock analyses reveal the emergence of Kaloula spp. common ancestor in East and Southeast Asia between the Eocene and Oligocene, c. 38.47 Ma (24.69-53.65). The genetic diversification of lineages then increased on the East Asian Mainland during the Lower Miocene, c. 20.10 (8.73-30.65), most likely originating from the vicariance and radiation triggered by the orogeny of the Qinghai-Tibetan Plateau. Later, the dispersal towards the North East Asian Mainland during the Upper Miocene drove the population diversification of K. borealis c. 9.01 Ma (3.66-15.29). Finally, the central mainland population became isolated following orogenesis events and diverged into K. rugifera during the Pliocene, c. 3.06 Ma (0.02-10.90). The combination of population genetic and barrier analyses revealed a significant genetic isolation between populations of Kaloula spp. matching with the massive Qinling-Daba Mountain chain located in south-central China. Finally, we highlight a young divergence within the Eastern Mainland population of K. borealis, possibly attributed to refugia in south eastern China from which populations later expanded.

Chytridiomycosis-induced mortality in a threatened anuran

Effectively planning conservation introductions involves assessing the suitability of both donor and recipient populations, including the landscape of disease risk. Chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), has caused extensive amphibian declines globally and may hamper reintroduction attempts. To determine Bd dynamics in potential source populations for conservation translocations of the threatened California red-legged frog (Rana draytonii) to Yosemite National Park, we conducted Bd sampling in two populations in the foothills of the Sierra Nevada Mountains, California, U.S.A. At one of two sites, we observed lethally high Bd loads in early post-metamorphic life stages and confirmed one chytridiomycosis-induced mortality, the first such report for this species. These results informed source population site selection for subsequent R. draytonii conservation translocations. Conservation efforts aimed at establishing new populations of R. draytonii in a landscape where Bd is ubiquitous can benefit from an improved understanding of risk through disease monitoring and ex situ infection studies.

Unequal sisters - Past and potential future range development of Anatolian and Hyrcanian brown frogs

Phylogeography can reconstruct historical evolutionary processes by comparing historical patterns of gene flow, divergence among species and by using species distribution models (SDM) upon geographic distribution. We investigate the phylogeographic patterns of Anatolian brown frogs including R. macrocnemis and R. tavasensis as well as the Hyrcanian brown frog, R. pseudodalmatina, using a fragment of the mitochondrial 16S rRNA gene for 145 specimens across the entire range of these frogs. We calculate parameters of molecular diversity, such as the number of variable sites (S), the number of haplotypes (h), haplotype diversity (Hd) and nucleotide diversity (π). We generated a haplotype network and used three methods (Neutrality tests, mismatch distributions and Bayesian skyline plots) to reconstruct the demographic histories of R. macrocnemis and R. pseudodalmatina. Finally, we used SDMs to predict the habitat suitability for three periods: The Present Day, the Last Glacial Maximum (LGM) and the future until 2070 for R. macrocnemis and R. pseudodalmatina. Our phylogenetic analyses support a late Miocene origin of Anatolian and Hyrcanian lineages. Hyrcanian brown frogs were enclosed in lowlands of the southern coast of the Caspian Sea after the uplift of the Elburz range and the Armenian plateau. The formation of a salinity belt from the north Aegean corridor (the south western Turkey) to northward during the Late Tortonian led to the subdivision of ancestor of the Anatolian lineage into today isolated western and eastern populations. The salinity belt had a considerable impact on the divergence of R. tavasensis from R. macrocnemis. Combined historical demographic analyses and SDMs revealed a rapid expansion occurring during the Pleistocene in R. macrocnemis and R. pseudodalmatina. Currently, suitable habitat for R. macrocnemis has declined compared to the LGM, and the species is predicted to do even worse under future climatic conditions. In contrast, R. pseudodalmatina found suitable habitat from the LGM to present within its restricted distribution area; it is predicted to do fine even under future climate.

Overwintering and cold tolerance in the Moor Frog (Rana arvalis) across its range

Only two species of boreal Holarctic frogs (genus Rana Linnaeus, 1758) can survive freezing and overwinter on land; they are found in the subarctic and cold regions of North America (Wood Frog, Rana sylvatica LeConte, 1825) and Eurasia (Moor Frog, Rana arvalis Nilsson, 1842) and are an example of an unusual adaptive strategy of overwintering. Freeze tolerance (down to –16 °C) of R. sylvatica has been thoroughly studied; however, little is known about cold resistance of R. arvalis in cold regions. We found that R. arvalis from European Russia and from West Siberia tolerate freezing down to –12 or –16 °C, whereas frogs from the Danish population survived freezing only to –4 °C (Y. Voituron et al. 2009b; J. Comp. Physiol. B, 179: 223–230). All of these populations, according to mitochondrial DNA markers, are closely related. We suggest that the observed differences in cold tolerance (–4 °C vs. –12 or –16 °C) could be caused either by adaptations to climatic factors or by differences in experimental protocols. The northeastern boundary of the geographic range of R. arvalis in Yakutia coincides with the transitional area between discontinuous and continuous permafrost; beyond this area, winter soil temperature sharply declines. The lower lethal temperature and overwintering ecology of R. arvalis in Siberia are similar to those of the North American R. sylvatica.

A new species of the Rana japonica group (Anura, Ranidae, Rana) from China, with a taxonomic proposal for the R. johnsi group

Ranajiulingensissp. nov., a new species from Hunan and Jiangxi, southeastern China, is described. The new species is assigned to the R.japonica group. The clade comprising R.jiulingensissp. nov. and R.dabieshanensis from Anhui is the sister taxon of R.omeimontis from Sichuan. Ranajiulingensissp. nov. can be distinguished by the significant divergences in the 16S and COI genes, and the combination of following morphological characters: body size medium, SVL 48.3–57.8 mm in adult males and 48.2–57.5 mm in adult females; dorsolateral fold straight; digits without circummarginal grooves; dorsal skin smooth; tibio-tarsal articulation reaching forward beyond the tip of snout; heels overlapping; webbing formula of toes: I 1⅓ – 2 II 1⅓ – 2⅓ III 1½ – 2⅔ IV 3 – 1⅓ V; absence of vocal sacs in males; and presence of creamy white nuptial pad with tiny hoar spines on the finger I and reddish tubercles on loreal and temporal regions in breeding males. Furthermore, based on our results and the previous literature, R.zhengi is synonymized with R.sangzhiensis, and a new species group, the Ranajohnsi group, is proposed for the R.johnsi and R.sangzhiensis. Currently, the Rana contains 41 recognized species, and the phylogenetic placements of several species remain unresolved.

Larger, unfiltered datasets are more effective at resolving phylogenetic conflict: Introns, exons, and UCEs resolve ambiguities in Golden-backed frogs (Anura: Ranidae; genus Hylarana)

Using FrogCap, a recently-developed sequence-capture protocol, we obtained >12,000 highly informative exons, introns, and ultraconserved elements (UCEs), which we used to illustrate variation in evolutionary histories of these classes of markers, and to resolve long-standing systematic problems in Southeast Asian Golden-backed frogs of the genus-complex Hylarana. We also performed a comprehensive suite of analyses to assess the relative performance of different genetic markers, data filtering strategies, tree inference methods, and different measures of branch support. To reduce gene tree estimation error, we filtered the data using different thresholds of taxon completeness (missing data) and parsimony informative sites (PIS). We then estimated species trees using concatenated datasets and Maximum Likelihood (IQ-TREE) in addition to summary (ASTRAL-III), distance-based (ASTRID), and site-based (SVDQuartets) multispecies coalescent methods. Topological congruence and branch support were examined using traditional bootstrap, local posterior probabilities, gene concordance factors, quartet frequencies, and quartet scores. Our results did not yield a single concordant topology. Instead, introns, exons, and UCEs clearly possessed different phylogenetic signals, resulting in conflicting, yet strongly-supported phylogenetic estimates. However, a combined analysis comprising the most informative introns, exons, and UCEs converged on a similar topology across all analyses, with the exception of SVDQuartets. Bootstrap values were consistently high despite high levels of incongruence and high proportions of gene trees supporting conflicting topologies. Although low bootstrap values did indicate low heuristic support, high bootstrap support did not necessarily reflect congruence or support for the correct topology. This study reiterates findings of some previous studies, which demonstrated that traditional bootstrap values can produce positively misleading measures of support in large phylogenomic datasets. We also showed a remarkably strong positive relationship between branch length and topological congruence across all datasets, implying that very short internodes remain a challenge to resolve, even with orders of magnitude more data than ever before. Overall, our results demonstrate that more data from unfiltered or combined datasets produced superior results. Although data filtering reduced gene tree incongruence, decreased amounts of data also biased phylogenetic estimation. A point of diminishing returns was evident, at which higher congruence (from more stringent filtering) at the expense of amount of data led to topological error as assessed by comparison to more complete datasets across different genomic markers. Additionally, we showed that applying a parameter-rich model to a partitioned analysis of concatenated data produces better results compared to unpartitioned, or even partitioned analysis using model selection. Despite some lingering uncertainties, a combined analysis of our genomic data and sequences supplemented from GenBank (on the basis of a few gene regions) revealed highly supported novel systematic arrangements. Based on these new findings, we transfer Amnirana nicobariensis into the genus Indosylvirana; and I. milleti and Hylarana celebensis to the genus Papurana. We also provisionally place H. attigua in the genus Papurana pending verification from positively identified (voucher substantiated) samples.

Vernal pool amphibian breeding ecology monitoring from 1931 to present: A harmonised historical and ongoing observational ecology dataset

Background

For 88 years (1931-present), the Mohonk Preserve's Daniel Smiley Research Center has been collecting data on occupancy and reproductive success of amphibian species, as well as associated water quality of 11 vernal pools each spring (February to May). Though sampling effort has varied over the dataset range, the size of the dataset is unprecedented within the field of amphibian ecology. With more than 2,480 individual species sampling dates and more than 151,701 recorded individual occurrences of the nine amphibian species, the described dataset represents the longest and largest time-series of herpetological sampling with paired water quality data.

New information

We describe the novel publication of a paired dataset of amphibian occurrence with environmental indicators spanning nearly 90 years of data collection. As of February 2020, the dataset includes 2,480 sampling dates across eleven vernal pools and 151,701 unique occurrences of egg masses or individuals recorded across nine species of amphibian. The dataset also includes environmental conditions associated with the species occurrences with complete coverage for air temperature and precipitation records and partial coverage for a variety of other weather and water quality measures. Data collection has included species, egg mass and tadpole counts; weather conditions including precipitation, sky and wind codes; water quality measurements including water temperature and pH; and vernal pool assessment including depth and surface vegetation coverage. Collection of data was sporadic from 1931–1991, but data have been collected consistently from 1991 to present. We also began monitoring dissolved oxygen, nitrate concentrations and conductivity of the vernal pools using a YSI Sonde Professional Plus Instrument and turbidity using a turbidity tube in February 2018. The dataset (and periodic updates), as well as metadata in the EML format, are available in the Environmental Data Initiative Repository under package edi.398.

The amphibians and reptiles of Colima, Mexico, with a summary of their conservation status

Colima is the fourth smallest Mexican state, covering only 0.3% of the surface area of Mexico, but due to the remarkable diversity of physiographic and environmental conditions present in Colima it contains a high biological diversity. We generated an up-to-date herpetofaunal checklist for Colima, with a summary of the conservation status of Colima's amphibians and reptiles. Our checklist contains a total of 153 species of amphibians and reptiles (three introduced). Thirty-nine are amphibians and 114 are reptiles. More than half of Colima's herpetofauna are Mexican endemics (66.7% of amphibians, 67.5% of reptiles). Less than 25% of the amphibian and reptile species in Colima are in protected categories according to the IUCN Red List and SEMARNAT. The reptiles in the Marine and Revillagigedo Archipelago regions are the most threatened taxa of the Colima herpetofauna. Colima shares > 80% of its herpetofauna with its neighboring states, Jalisco and Michoacán.

Antimicrobial peptide and sequence variation along a latitudinal gradient in two anurans

BACKGROUND

While there is evidence of both purifying and balancing selection in immune defense genes, large-scale genetic diversity in antimicrobial peptides (AMPs), an important part of the innate immune system released from dermal glands in the skin, has remained uninvestigated. Here we describe genetic diversity at three AMP loci (Temporin, Brevinin and Palustrin) in two ranid frogs (Rana arvalis and R. temporaria) along a 2000 km latitudinal gradient. We amplified and sequenced part of the Acidic Propiece domain and the hypervariable Mature Peptide domain (~ 150-200 bp) in the three genes using Illumina Miseq and expected to find decreased AMP genetic variation towards the northern distribution limit of the species similarly to studies on MHC genetic patterns.

RESULTS

We found multiple loci for each AMP and relatively high gene diversity, but no clear pattern of geographic genetic structure along the latitudinal gradient. We found evidence of trans-specific polymorphism in the two species, indicating a common evolutionary origin of the alleles. Temporin and Brevinin did not form monophyletic clades suggesting that they belong to the same gene family. By implementing codon evolution models we found evidence of strong positive selection acting on the Mature Peptide. We also found evidence of diversifying selection as indicated by divergent allele frequencies among populations and high Theta k values.

CONCLUSION

Our results suggest that AMPs are an important source of adaptive diversity, minimizing the chance of microorganisms developing resistance to individual peptides.

Circumventing surface tension: tadpoles suck bubbles to breathe air

The surface tension of water provides a thin, elastic membrane upon which many tiny animals are adapted to live and move. We show that it may be equally important to the minute animals living beneath it by examining air-breathing mechanics in five species (three families) of anuran (frog) tadpoles. Air-breathing is essential for survival and development in most tadpoles, yet we found that all tadpoles at small body sizes were unable to break through the water's surface to access air. Nevertheless, by 3 days post-hatch and only 3 mm body length, all began to breathe air and fill the lungs. High-speed macrovideography revealed that surface tension was circumvented by a novel behaviour we call 'bubble-sucking': mouth attachment to the water's undersurface, the surface drawn into the mouth by suction, a bubble 'pinched off' within the mouth, then compressed and forced into the lungs. Growing tadpoles transitioned to air-breathing via typical surface breaching. Salamander larvae and pulmonate snails were also discovered to 'bubble-suck', and two insects used other means of circumvention, suggesting that surface tension may have a broader impact on animal phenotypes than hitherto appreciated.

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as "refugia within refugia", provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and "ghost" mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD-sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply-diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation-by-distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.

Biophysical studies on the antimicrobial activity of linearized esculentin 2EM

Linearized esculentin 2 EM (E2EM-lin) from the frog, Glandirana emeljanovi was highly active against Gram-positive bacteria (minimum lethal concentration ≤ 5.0 μM) and strongly α-helical in the presence of lipid mimics of their membranes (>55.0%). The N-terminal α-helical structure adopted by E2EM-lin showed the potential to form a membrane interactive, tilted peptide with an hydrophobicity gradient over residues 9 to 23. E2EM-lin inserted strongly into lipid mimics of membranes from Gram-positive bacteria (maximal surface pressure changes ≥5.5 mN m^-1), inducing increased rigidity (C_s^-1 ↑), thermodynamic instability (ΔG_mix < 0 → ΔG_mix > 0) and high levels of lysis (>50.0%). These effects appeared to be driven by the high anionic lipid content of membranes from Gram-positive bacteria; namely phosphatidylglycerol (PG) and cardiolipin (CL) species. The high levels of α-helicity (60.0%), interaction (maximal surface pressure change = 6.7 mN m^-1) and lysis (66.0%) shown by E2EM-lin with PG species was a major driver in the ability of the peptide to lyse and kill Gram-positive bacteria. E2EM-lin also showed high levels of α-helicity (62.0%) with CL species but only low levels of interaction (maximal surface pressure change = 2.9 mN m^-1) and lysis (21.0%) with the lipid. These combined data suggest that E2EM-lin has a specificity for killing Gram-positive bacteria that involves the formation of tilted structure and appears to be primarily driven by PG-mediated membranolysis. These structure/function relationships are used to help explain the pore forming process proposed to describe the membranolytic, antibacterial action of E2EM-lin.

Spatiotemporal heterogeneity decouples infection parameters of amphibian chytridiomycosis

Emerging infectious diseases are responsible for declines in wildlife populations around the globe. Mass mortality events associated with emerging infectious diseases are often associated with high number of infected individuals (prevalence) and high pathogen loads within individuals (intensity). At the landscape scale, spatial and temporal variation in environmental conditions can alter the relationship between these infection parameters and blur the overall picture of disease dynamics. Quantitative estimates of how infection parameters covary with environmental heterogeneity at the landscape scale are scarce. If we are to identify wild populations at risk of disease epidemics, we must elucidate the factors that shape, and potentially decouple, the link between pathogen prevalence and intensity of infection over complex ecological landscapes. Using a network of 41 populations of the amphibian host Rana pipiens in Ontario, Canada, we present the spatial and temporal heterogeneity in pathogen prevalence and intensity of infection of the chytrid fungus Batrachochytrium dendrobatidis (Bd), across a 3-year period. We then quantify how covariation between both infection parameters measured during late summer is modified by previously experienced spatiotemporal environmental heterogeneity across 14 repeat sampled populations. Late summer Bd infection parameters are governed, at least in part, by different environmental factors operating during separate host life-history events. Our results provide evidence for a relationship between Bd prevalence and thermal regimes prior to host breeding at the site level, and a relationship between intensity of infection and aquatic conditions (precipitation, hydroshed size and river density) throughout host breeding period at the site level. This demonstrates that microclimatic variation within temporal windows can drive divergent patterns of pathogen dynamics within and across years, by effecting changes in host behaviour which interfere with the pathogen's ability to infect and re-infect hosts. A clearer understanding of the role that spatiotemporal heterogeneity has upon infection parameters will provide valuable insights into host-pathogen epidemiology, as well as more fundamental aspects of the ecology and evolution of interspecific interactions.

Parachute geckos free fall into synonymy: Gekko phylogeny, and a new subgeneric classification, inferred from thousands of ultraconserved elements

Recent phylogenetic studies of gekkonid lizards have revealed unexpected, widespread paraphyly and polyphyly among genera, unclear generic boundaries, and a tendency towards the nesting of taxa exhibiting specialized, apomorphic morphologies within geographically widespread "generalist" clades. This is especially true in Australasia, where monophyly of Gekko proper has been questioned with respect to phenotypically ornate flap-legged geckos of the genus Luperosaurus, the Philippine false geckos of the genus Pseudogekko, and even the elaborately "derived" parachute geckos of the genus Ptychozoon. Here we employ sequence capture targeting 5060 ultraconserved elements (UCEs) to infer phylogenomic relationships among 42 representative ingroup gekkonine lizard taxa. We analyze multiple datasets of varying degrees of completeness (10, 50, 75, 95, and 100 percent complete with 4715, 4051, 3376, 2366, and 772 UCEs, respectively) using concatenated maximum likelihood and multispecies coalescent methods. Our sampling scheme addresses four persistent systematic questions in this group: (1) Are Luperosaurus and Ptychozoon monophyletic, and are any of these named species truly nested within Gekko? (2) Are prior phylogenetic estimates of Sulawesi's L. iskandari as the sister taxon to Melanesian G. vittatus supported by our genome-scale dataset? (3) Is the high-elevation L. gulat of Palawan Island correctly placed within Gekko? (4) And, finally, where do the enigmatic taxa P. rhacophorus and L. browni fall in a higher-level gekkonid phylogeny? We resolve these issues; confirm with strong support some previously inferred findings (placement of Ptychozoon taxa within Gekko; the sister taxon relationship between L. iskandari and G. vittatus); resolve the systematic position of unplaced taxa (L. gulat, and L. browni); and transfer L. iskandari, L. gulat, L. browni, and all members of the genus Ptychozoon to the genus Gekko. Our unexpected and novel systematic inference of the placement of Ptychozoon rhacophorus suggests that this species is not grouped with Ptychozoon or even Luperosaurus (as previously expected) but may, in fact, be most closely related to several Indochinese species of Gekko. With our resolved and strongly supported phylogeny, we present a new classification emphasizing the most inclusive, original generic name (Gekko) for these ~60 taxa, arranged into seven subgenera.

An integrative phylogenomic approach illuminates the evolutionary history of Old World tree frogs (Anura: Rhacophoridae)

Rhacophoridae are one of the most speciose and ecologically diverse families of amphibians. Resolution of their evolutionary relationships is key to understanding the accumulation of biodiversity, yet previous hypotheses based on Sanger sequencing exhibit much discordance amongst generic relationships. This conflict precludes the making of sound macroevolutionary conclusions. Herein, we conduct the first phylogenomic study using broad-scale sampling and sequences of 352 nuclear DNA loci obtained using anchored hybrid enrichment targeted sequencing. The robust time-calibrated phylogenetic hypothesis clarifies several long-disputed relationships and facilitates the testing of evolutionary hypotheses on spatiotemporal diversification and reproductive modes. The major extant lineages of Rhacophoridae appear to have radiated in mainland Asia, and the spatiotemporal process corresponds with several common accumulations of biodiversity in Asia. Analyses do not detect any case of "Out of Himalaya" in Rhacophoridae. All transitions of reproductive modes appear to have evolved in an ordered, gradual sequence associated with gaining independence of standing water for larval development. The different reproductive modes are phylogenetically conserved and the completion of their transitions appear to have occurred over a period of ~30 Ma, which does not fit a pattern of a rapid burst of diversification. Innovations in reproductive modes associate statistically with the uneven distribution of species-richness between clades, where higher diversification is linked to increased terrestrial modes of reproduction. These results strengthen the hypothesis that breeding innovations drive diversification by providing new opportunities for ecological release and dispersion.