Evolutionary and biogeographic origins of high tropical diversity in old world frogs (Ranidae)

Complete Mitogenome of Oreolalax Omeimontis Reveals Phylogenetic Status and Novel Gene Arrangement of Archaeobatrachia

Species of the genus Oreolalax displayed crucial morphological characteristics of vertebrates transitioning from aquatic to terrestrial habitats; thus, they can be regarded as a representative vertebrate genus for this landing phenomenon. But the present phylogenetic status of Oreolalax omeimontis has been controversial with morphological and molecular approaches, and specific gene rearrangements were discovered in all six published Oreolalax mitogenomes, which are rarely observed in Archaeobatrachia. Therefore, this study determined the complete mitogenome of O. omeimontis with the aim of identifying its precise phylogenetic position and novel gene arrangement in Archaeobatrachia. Phylogenetic analysis with Bayesian inference and maximum likelihood indicates O. omeimontis is a sister group to O. lichuanensis, which is consistent with previous phylogenetic analysis based on morphological characteristics, but contrasts with other studies using multiple gene fragments. Moreover, although the duplication of trnM occurred in all seven Oreolalax species, the translocation of trnQ and trnM occurred differently in O. omeimontis to the other six, and this unique rearrangement would happen after the speciation of O. omeimontis. In general, this study sheds new light on the phylogenetic relationships and gene rearrangements of Archaeobatrachia.

Hidden in plain sight: novel molecular data reveal unexpected genetic diversity among paramphistome parasites (Digenea: Paramphistomoidea) of European water frogs

Parasites might represent a helpful tool in understanding the historical dispersion and phylogeography of their hosts. In order to reveal whether the migration routes and diversification of hosts can be traceable in the genetic structures of their parasites, we investigated the diversity of paramphistomoid trematodes of Pelophylax frogs in 2 geographically distant European regions. Water frogs belonging to the genus Pelophylax represent a striking example of a species with a high variety of ecological adaptations and a rich evolutionary history. The parasites were collected from 2 Balkan endemic species, P. epeiroticus and P. kurtmuelleri, and 2 species in Slovakia, P. esculentus and P. ridibundus. While in Slovakia, Pelophylax frogs harboured 2 species, the diplodiscid Diplodiscus subclavatus and the cladorchiid Opisthodiscus diplodiscoides, only the former was recorded in the south-western Balkans. Remarkably high genetic diversity (16 unique mitochondrial cox1 haplotypes, recognized among 60 novel sequences) was observed in D. subclavatus, and subsequent phylogenetic analyses revealed a strong population-genetic structure associated with geographical distribution. We also evidenced the existence of 2 divergent D. subclavatus cox1 haplogroups in the south-western Balkans, which might be associated with the historical diversification of endemic water frogs in the regional glacial microrefugia.

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.

Conserved ancestral tropical niche but different continental histories explain the latitudinal diversity gradient in brush-footed butterflies

The global increase in species richness toward the tropics across continents and taxonomic groups, referred to as the latitudinal diversity gradient, stimulated the formulation of many hypotheses to explain the underlying mechanisms of this pattern. We evaluate several of these hypotheses to explain spatial diversity patterns in a butterfly family, the Nymphalidae, by assessing the contributions of speciation, extinction, and dispersal, and also the extent to which these processes differ among regions at the same latitude. We generate a time-calibrated phylogeny containing 2,866 nymphalid species (~45% of extant diversity). Neither speciation nor extinction rate variations consistently explain the latitudinal diversity gradient among regions because temporal diversification dynamics differ greatly across longitude. The Neotropical diversity results from low extinction rates, not high speciation rates, and biotic interchanges with other regions are rare. Southeast Asia is also characterized by a low speciation rate but, unlike the Neotropics, is the main source of dispersal events through time. Our results suggest that global climate change throughout the Cenozoic, combined with tropical niche conservatism, played a major role in generating the modern latitudinal diversity gradient of nymphalid butterflies.

Relict groups of spiny frogs indicate Late Paleogene-Early Neogene trans-Tibet dispersal of thermophile faunal elements

BACKGROUND

The Himalaya-Tibet orogen (HTO) presents an outstanding geologically active formation that contributed to, and fostered, modern Asian biodiversity. However, our concepts of the historical biogeography of its biota are far from conclusive, as are uplift scenarios for the different parts of the HTO. Here, we revisited our previously published data set of the tribe Paini extending it with sequence data from the most western Himalayan spiny frogs Allopaa and Chrysopaa and using them as an indirect indicator for the potential paleoecological development of Tibet.

METHODS

We obtained sequence data of two mitochondrial loci (16S rRNA, COI) and one nuclear marker (Rag1) from Allopaa samples from Kashmir Himalaya as well as Chrysopaa sequence data from the Hindu Kush available from GenBank to complement our previous data set. A Maximum likelihood and dated Bayesian gene tree were generated based on the concatenated data set. To resolve the inconsistent placement of Allopaa, we performed different topology tests.

RESULTS

Consistent with previous results, the Southeast Asian genus Quasipaa is sister to all other spiny frogs. The results further reveal a basal placement of Chrysopaa relative to Allopaa and Nanorana with an estimated age of ca. 26 Mya. Based on the topology tests, the phylogenetic position of Allopaa as a sister clade to Chaparana seems to be most likely, resulting in a paraphyletic genus Nanorana and a separation from the latter clade around 20 Mya, although a basal position of Allopaa to the genus Nanorana cannot be entirely excluded. Both, the placements of Chrysopaa and Allopaa support the presence of basal Paini lineages in the far northwestern part of the HTO, which is diametrically opposite end of the HTO with respect to the ancestral area of spiny frogs in Southeast Asia. These striking distributional patterns can be most parsimoniously explained by trans-Tibet dispersal during the late Oligocene (subtropical Chrysopaa) respectively early Miocene (warm temperate Allopaa). Within spiny frogs, only members of the monophyletic Nanorana+Paa clade are adapted to the colder temperate climates, indicating that high-altitude environments did not dominate in the HTO before ca. 15 Mya. Our results are consistent with fossil records suggesting that large parts of Tibet were characterized by subtropical to warm temperate climates at least until the early Miocene. They contradict prevalent geological models of a highly uplifted late Paleogene proto-Plateau.

Unexpectedly high levels of lineage diversity in Sundaland puddle frogs (Dicroglossidae: Occidozyga Kuhl and van Hasselt, 1822)

One of the most urgent contemporary tasks for taxonomists and evolutionary biologists is to estimate the number of species on earth. Recording alpha diversity is crucial for protecting biodiversity, especially in areas of elevated species richness, which coincide geographically with increased anthropogenic environmental pressures - the world's so-called biodiversity hotspots. Although the distribution of Puddle frogs of the genus Occidozyga in South and Southeast Asia includes five biodiversity hotspots, the available data on phylogeny, species diversity, and biogeography are surprisingly patchy. Samples analyzed in this study were collected throughout Southeast Asia, with a primary focus on Sundaland and the Philippines. A mitochondrial gene region comprising ~ 2000 bp of 12S and 16S rRNA with intervening tRNA Valine and three nuclear loci (BDNF, NTF3, POMC) were analyzed to obtain a robust, time-calibrated phylogenetic hypothesis. We found a surprisingly high level of genetic diversity within Occidozyga, based on uncorrected p-distance values corroborated by species delimitation analyses. This extensive genetic diversity revealed 29 evolutionary lineages, defined by the > 5% uncorrected p-distance criterion for the 16S rRNA gene, suggesting that species diversity in this clade of phenotypically homogeneous forms probably has been underestimated. The comparison with results of other anuran groups leads to the assumption that anuran species diversity could still be substantially underestimated in Southeast Asia in general. Many genetically divergent lineages of frogs are phenotypically similar, indicating a tendency towards extensive morphological conservatism. We present a biogeographic reconstruction of the colonization of Sundaland and nearby islands which, together with our temporal framework, suggests that lineage diversification centered on the landmasses of the northern Sunda Shelf. This remarkably genetically structured group of amphibians could represent an exceptional case for future studies of geographical structure and diversification in a widespread anuran clade spanning some of the most pronounced geographical barriers on the planet (e.g., Wallace's Line). Studies considering gene flow, morphology, ecological and bioacoustic data are needed to answer these questions and to test whether observed diversity of Puddle frog lineages warrants taxonomic recognition.

Incorporating the speciation process into species delimitation

The "multispecies" coalescent (MSC) model that underlies many genomic species-delimitation approaches is problematic because it does not distinguish between genetic structure associated with species versus that of populations within species. Consequently, as both the genomic and spatial resolution of data increases, a proliferation of artifactual species results as within-species population lineages, detected due to restrictions in gene flow, are identified as distinct species. The toll of this extends beyond systematic studies, getting magnified across the many disciplines that rely upon an accurate framework of identified species. Here we present the first of a new class of approaches that addresses this issue by incorporating an extended speciation process for species delimitation. We model the formation of population lineages and their subsequent development into independent species as separate processes and provide for a way to incorporate current understanding of the species boundaries in the system through specification of species identities of a subset of population lineages. As a result, species boundaries and within-species lineages boundaries can be discriminated across the entire system, and species identities can be assigned to the remaining lineages of unknown affinities with quantified probabilities. In addition to the identification of species units in nature, the primary goal of species delimitation, the incorporation of a speciation model also allows us insights into the links between population and species-level processes. By explicitly accounting for restrictions in gene flow not only between, but also within, species, we also address the limits of genetic data for delimiting species. Specifically, while genetic data alone is not sufficient for accurate delimitation, when considered in conjunction with other information we are able to not only learn about species boundaries, but also about the tempo of the speciation process itself.

Revisiting the phylogenetic predicament of the genus Huia (Amphibia: Ranidae) using molecular data and tadpole morphology

Despite a considerable amount of research, the systematics of the ranid genus Huia have remained unresolved, mostly owing to insufficient sampling and morphological similarities. As currently circumscribed, Huia consists of five species, but multiple studies have consistently demonstrated that it is not a monophyletic genus. However, no study has approached the problem with adequate data and provided a systematically sound solution, leaving the genus to languish in a classification that is phylogenetically incoherent. We generated the most comprehensive sampling of Huia to date, based on extensive fieldwork in Java and Sumatra. Using an integrative taxonomy framework, we analysed four mitochondrial and two nuclear markers and, in conjunction with tadpole morphology, investigated the phylogenetics of Huia and its congeners. Corroborating previous studies, Huia is recovered as a paraphyletic group. Huia cavitympanum emerges as the sister taxon to Meristogenys. The remaining members of Huia form a monophyletic group, sister to the H. cavitympanum + Meristogenys clade. Our extensive geographical sampling in Sumatra and Java reveals multiple highly divergent lineages that potentially represent undescribed species. Using our expanded molecular and morphological dataset, we resolve the paraphyly of Huia by restricting the genus to its type species and propose a new genus to accommodate the strongly supported clade of Sumatran and Javan populations previously belonging to Huia.

Characteristics of the mitochondrial genome of Rana omeimontis and related species in Ranidae: Gene rearrangements and phylogenetic relationships

The Omei wood frog (Rana omeimontis), endemic to central China, belongs to the family Ranidae. In this study, we achieved detail knowledge about the mitogenome of the species. The length of the genome is 20,120 bp, including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a noncoding control region. Similar to other amphibians, we found that only nine genes (ND6 and eight tRNA genes) are encoded on the light strand (L) and other genes on the heavy strand (H). Totally, The base composition of the mitochondrial genome included 27.29% A, 28.85% T, 28.87% C, and 15.00% G, respectively. The control regions among the Rana species were found to exhibit rich genetic variability and A + T content. R. omeimontis was clustered together with R. chaochiaoensis in phylogenetic tree. Compared to R. amurensis and R. kunyuensi, it was more closely related to R. chaochiaoensis, and a new way of gene rearrangement (ND6-trnE-Cytb-D-loop-trnL2 (CUN)-ND5-D-loop) was also found in the mitogenome of R. amurensis and R. kunyuensi. Our results about the mitochondrial genome of R. omeimontis will contribute to the future studies on phylogenetic relationship and the taxonomic status of Rana and related Ranidae species.

Phylogeographic analysis and species distribution modelling of the wood frog Batrachyla leptopus (Batrachylidae) reveal interglacial diversification in south western Patagonia

Background

The evolutionary history of southern South American organisms has been strongly influenced by Pleistocene climate oscillations. Amphibians are good models to evaluate hypotheses about the influence of these climate cycles on population structure and diversification of the biota, because they are sensitive to environmental changes and have restricted dispersal capabilities. We test hypotheses regarding putative forest refugia and expansion events associated with past climatic changes in the wood frog Batrachyla leptopus distributed along ∼1,000 km of length including glaciated and non-glaciated areas in southwestern Patagonia.

Methods

Using three mitochondrial regions (D-loop, cyt b, and coI) and two nuclear loci (pomc and crybA1), we conducted multilocus phylogeographic analyses and species distribution modelling to gain insights of the evolutionary history of this species. Intraspecific genealogy was explored with maximum likelihood, Bayesian, and phylogenetic network approaches. Diversification time was assessed using molecular clock models in a Bayesian framework, and demographic scenarios were evaluated using approximate Bayesian computation (ABC) and extended Bayesian skyline plot (EBSP). Species distribution models (SDM) were reconstructed using climatic and geographic data.

Results

Population structure and genealogical analyses support the existence of four lineages distributed north to south, with moderate to high phylogenetic support (Bootstrap > 70%; BPP > 0.92). The diversification time of B. leptopus’ populations began at ∼0.107 mya. The divergence between A and B lineages would have occurred by the late Pleistocene, approximately 0.068 mya, and divergence between C and D lineages was approximately 0.065 mya. The ABC simulations indicate that lineages coalesced at two different time periods, suggesting the presence of at least two glacial refugia and a postglacial colonization route that may have generated two southern lineages (p = 0.93, type I error: <0.094, type II error: 0.134). EBSP, mismatch distribution and neutrality indexes suggest sudden population expansion at ∼0.02 mya for all lineages. SDM infers fragmented distributions of B. leptopus associated with Pleistocene glaciations. Although the present populations of B. leptopus are found in zones affected by the last glacial maximum (∼0.023 mya), our analyses recover an older history of interglacial diversification (0.107–0.019 mya). In addition, we hypothesize two glacial refugia and three interglacial colonization routes, one of which gave rise to two expanding lineages in the south.

Diversity of Antimicrobial Peptides in Three Partially Sympatric Frog Species in Northeast Asia and Implications for Evolution

Antimicrobial peptides (AMPs) are evolutionarily ancient molecules that play an essential role in innate immunity across taxa from invertebrates to vertebrates. The evolution system of AMP system has not been well explained in the literature. In this study, we cloned and sequenced AMP transcriptomes of three frog species, namely Rana dybowskii, Rana amurensis, and Pelophylax nigromaculatus, which are partially sympatric in northeast Asia, but show different habitat preferences. We found that each species contained 7 to 14 families of AMPs and the diversity was higher in species with a large geographic range and greater habitat variation. All AMPs are phylogenetically related but not associated with the speciation process. Most AMP genes were under negative selection. We propose that the diversification and addition of novel functions and improvement of antimicrobial efficiency are facilitated by the expansion of family members and numbers. We also documented significant negative correlation of net charges and numbers of amino acid residues between the propiece and mature peptide segments. This supports the Net Charge Balance Hypothesis. We propose the Cut Point Sliding Hypothesis as a novel diversification mechanism to explain the correlation in lengths of the two segments.

Phylogeny of spiny frogs Nanorana (Anura: Dicroglossidae) supports a Tibetan origin of a Himalayan species group

Recent advances in the understanding of the evolution of the Asian continent challenge the long-held belief of a faunal immigration into the Himalaya. Spiny frogs of the genus Nanorana are a characteristic faunal group of the Himalaya-Tibet orogen (HTO). We examine the phylogeny of these frogs to explore alternative biogeographic scenarios for their origin in the Greater Himalaya, namely, immigration, South Tibetan origin, strict vicariance. We sequenced 150 Nanorana samples from 62 localities for three mitochondrial (1,524 bp) and three nuclear markers (2,043 bp) and complemented the data with sequence data available from GenBank. We reconstructed a gene tree, phylogenetic networks, and ancestral areas. Based on the nuDNA, we also generated a time-calibrated species tree. The results revealed two major clades (Nanorana and Quasipaa), which originated in the Lower Miocene from eastern China and subsequently spread into the HTO (Nanorana). Five well-supported subclades are found within Nanorana: from the East, Central, and Northwest Himalaya, the Tibetan Plateau, and the southeastern Plateau margin. The latter subclade represents the most basal group (subgenus Chaparana), the Plateau group (Nanorana) represents the sister clade to all species of the Greater Himalaya (Paa). We found no evidence for an east-west range expansion of Paa along the Himalaya, nor clear support for a strict vicariance model. Diversification in each of the three Himalayan subclades has probably occurred in distinct areas. Specimens from the NW Himalaya are placed basally relative to the highly diverse Central Himalayan group, while the lineage from the Tibetan Plateau is placed within a more terminal clade. Our data indicate a Tibetan origin of Himalayan Nanorana and support a previous hypothesis, which implies that a significant part of the Himalayan biodiversity results from primary diversification of the species groups in South Tibet before this part of the HTO was uplifted to its recent heights.

The complete mitochondrial genome of Pyxicephalus adspersus: high gene rearrangement and phylogenetics of one of the world's largest frogs

The family Pyxicephalidae including two subfamilies (Cacosterninae and Pyxicephalinae) is an ecologically important group of frogs distributed in sub-Saharan Africa. However, its phylogenetic position among the Anura has remained uncertain. The present study determined the complete mitochondrial genome sequence of Pyxicephalus adspersus, the first representative mitochondrial genome from the Pyxicephalinae, and reconstructed the phylogenetic relationships within Ranoidae using 10 mitochondrial protein-coding genes of 59 frog species. The P. adspersus mitochondrial genome showed major gene rearrangement and an exceptionally long length that is not shared with other Ranoidae species. The genome is 24,317 bp in length, and contains 15 protein-coding genes (including extra COX3 and Cyt b genes), four rRNA genes (including extra 12S rRNA and 16S rRNA genes), 29 tRNA genes (including extra tRNALeu (UAG), tRNALeu (UUR), tRNAThr , tRNAPro , tRNAPhe , tRNAVal , tRNAGln genes) and two control regions (CRs). The Dimer-Mitogenome and Tandem duplication and random loss models were used to explain these gene arrangements. Finally, both Bayesian inference and maximum likelihood analyses supported the conclusion that Pyxicephalidae was monophyletic and that Pyxicephalidae was the sister clade of (Petropedetidae + Ptychadenidae).

Identification of leopard frogs (Anura: Ranidae: Lithobates) distributed in some localities of the Southern Mexican Plateau using mitochondrial DNA sequences

Species of the genus Lithobates are difficult to identify, especially on the 'Rana pipiens' complex, because the morphological differences among the species are often subtle. In this study, we used 12S ribosomal RNA gene partial sequences to identify specimens of leopard frogs from some new localities (not sampled on previous studies) of the Southern Mexican Plateau, to carry out a phylogenetic analysis, and also a morphometric analysis of some morphological features were analyzed to evaluate their morphological variation. A phylogenetic analysis using partial sequences of 12S rDNA mitochondrial gene for 31 samples was carried out using Bayesian Inference, Maximum Parsimony, and Maximum Likelihood. In addition, 20 morphological lineal measurements from 97 specimens were analyzed by principal component analysis (PCA) to evaluate if the frogs studied are conspecific. Partial sequences of the 12S rDNA obtained from frogs distributed in the Southern Mexican Plateau show two haplotypes with a divergence of 0.4%. Phylogenetic hypothesis shows an exclusive group with a previous sequence of Lithobates montezumae. The PCA indicates that variables are not linearly correlated and specimens belong to a single group. Evidence found, let us consider that specimens studied belong to one species of the Lithobates montezumae subgroup, distributed principally in the Southern Mexican Plateau.

Timing and ecological priority shaped the diversification of sedges in the Himalayas

BACKGROUND

Diversification patterns in the Himalayas have been important to our understanding of global biodiversity. Despite recent broad-scale studies, the most diverse angiosperm genus of the temperate zone-Carex L. (Cyperaceae), with ca. 2100 species worldwide-has not yet been studied in the Himalayas, which contains 189 Carex species. Here the timing and phylogenetic pattern of lineage and ecological diversification were inferred in this ecologically significant genus. We particularly investigated whether priority, adaptation to ecological conditions, or both explain the highly successful radiation of the Kobresia clade (ca. 60 species, of which around 40 are present in the Himalayas) of Himalayan Carex.

METHODS

Phylogenetic relationships were inferred using maximum likelihood analysis of two nuclear ribosomal DNA (nrDNA) regions (ITS and ETS) and one plastid gene (matK); the resulting tree was time-calibrated using penalized likelihood and a fossil calibration at the root of the tree. Biogeographical reconstruction for estimation of historical events and ancestral ranges was performed using the dispersal-extinction-cladogenesis (DEC) model, and reciprocal effects between biogeography and diversification were inferred using the geographic state speciation and extinction (GeoSSE) model. Climatic envelopes for all species for which mapped specimen data available were estimated using climatic data from WORLDCLIM, and climatic niche evolution was inferred using a combination of Ornstein-Uhlenbeck models of shifting adaptive optima and maximum likelihood inference of ancestral character states under a Brownian motion model.

RESULTS

The Himalayan Carex flora represents three of the five major Carex clades, each represented by multiple origins within the Himalayas. The oldest Carex radiation in the region, dating to ca. 20 Ma, near the time of Himalayan orogeny, gave rise to the now abundant Kobresia clade via long-distance dispersal from the Nearctic. The Himalayan Carex flora comprises a heterogeneous sample of diversifications drawn from throughout the cosmopolitan, but mostly temperate, Carex radiation. Most radiations are relatively recent, but the widespread and diverse Himalayan Kobresia radiation arose at the early Miocene. The timing and predominance of Kobresia in high-elevation Himalayan meadows suggests that Kobresia may have excluded other Carex lineages: the success of Kobresia in the Himalayas, in other words, appears to be a consequence largely of priority, competitive exclusion and historical contingency.

Sequence capture using AFLP-generated baits: A cost-effective method for high-throughput phylogenetic and phylogeographic analysis

Target sequence capture is an efficient technique to enrich specific genomic regions for high-throughput sequencing in ecological and evolutionary studies. In recent years, many sequence capture approaches have been proposed, but most of them rely on commercial synthetic baits which make the experiment expensive. Here, we present a novel sequence capture approach called AFLP-based genome sequence capture (AFLP Capture). This method uses the AFLP (amplified fragment length polymorphism) technique to generate homemade capture baits without the need for prior genome information, thus is applicable to any organisms. In this approach, biotinylated AFLP fragments representing a random fraction of the genome are used as baits to capture the homologous fragments from genomic shotgun sequencing libraries. In a trial study, by using AFLP Capture, we successfully obtained 511 orthologous loci (>700,000 bp in total length) from 11 Odorrana species and more than 100,000 single nucleotide polymorphisms (SNPs) in four analyzed individuals of an Odorrana species. This result shows that our method can be used to address questions of various evolutionary depths (from interspecies level to intraspecies level). We also discuss the flexibility in bait preparation and how the sequencing data are analyzed. In summary, AFLP Capture is a rapid and flexible tool and can significantly reduce the experimental cost for phylogenetic studies that require analyzing genome-scale data (hundreds or thousands of loci).

A new ancient lineage of frog (Anura: Nyctibatrachidae: Astrobatrachinae subfam. nov.) endemic to the Western Ghats of Peninsular India

The Western Ghats (WG) is an escarpment on the west coast of Peninsular India, housing one of the richest assemblages of frogs in the world, with three endemic families. Here, we report the discovery of a new ancient lineage from a high-elevation massif in the Wayanad Plateau of the southern WG. Phylogenetic analysis reveals that the lineage belongs to Natatanura and clusters with Nyctibatrachidae, a family endemic to the WG/Sri Lanka biodiversity hotspot. Based on geographic distribution, unique morphological traits, deep genetic divergence, and phylogenetic position that distinguishes the lineage from the two nyctibatrachid subfamilies Nyctibatrachinae Blommers-Schlösser, 1993 and Lankanectinae Dubois & Ohler, 2001, we erect a new subfamily Astrobatrachinae subfam. nov. (endemic to the WG, Peninsular India), and describe a new genus Astrobatrachus gen. nov. and species, Astrobatrachus kurichiyana sp. nov. The discovery of this species adds to the list of deeply divergent and monotypic or depauperate lineages with narrow geographic ranges in the southern massifs of the WG. The southern regions of the WG have long been considered geographic and climatic refugia, and this new relict lineage underscores their evolutionary significance. The small range of this species exclusively outside protected areas highlights the significance of reserve forest tracts in the WG in housing evolutionary novelty. This reinforces the need for intensive sampling to uncover new lineages and advance our understanding of the historical biogeography of this ancient landmass.

Genetic variation and cryptic lineage diversity of the Nigerian red-headed rock agama Agama agama associate with eco-geographic zones

Nigeria is an Afrotropical region with considerable ecological heterogeneity and levels of biotic endemism. Among its vertebrate fauna, reptiles have broad distributions, thus, they constitute a compelling system for assessing the impact of ecological variation and geographic isolation on species diversification. The red-headed rock agama, Agama agama, lives in a wide range of habitats and, thus, it may show genetic structuring and diversification. Herein, we tested the hypothesis that ecology affects its genetic structure and population divergence. Bayesian inference phylogenetic analysis of a mitochondrial DNA (mtDNA) gene recovered four well-supported matrilines with strong evidence of genetic structuring consistent with eco-geographic regions. Genetic differences among populations based on the mtDNA also correlated with geographic distance. The ecological niche model for the matrilines had a good fit and robust performance. Population divergence along the environmental axes was associated with climatic conditions, and temperature ranked highest among all environmental variables for forest specialists, while precipitation ranked highest for the forest/derived savanna, and savanna specialists. Our results cannot reject the hypothesis that niche conservatism promotes geographic isolation of the western populations of Nigerian A. agama. Thus, ecological gradients and geographic isolation impact the genetic structure and population divergence of the lizards. This species might be facing threats due to recent habitat fragmentation, especially in western Nigeria. Conservation actions appear necessary.

Evolutionary biogeography of the centipede genus Ethmostigmus from Peninsular India: testing an ancient vicariance hypothesis for Old World tropical diversity

BACKGROUND

Understanding the relative influence of vicariance and dispersal in shaping Old World tropical biodiversity remains a challenge. We aimed to infer the roles of these alternative biogeographic processes using a species time-tree for the centipede genus Ethmostigmus from the Old World tropics. Additionally, we explored fine-scale biogeographic patterns for an endemic radiation of Ethmostigmus from the peninsular Indian Plate (PIP), an area with complex geological and climatic history.

RESULTS

Divergence time estimates suggest that Ethmostigmus began diversifying in the Late Cretaceous, 99 (± 25) million years ago (Ma), its early biogeographic history shaped by vicariance. Members of Ethmostigmus in PIP form a monophyletic group that underwent endemic radiation in the Late Cretaceous, 72 (± 25) Ma. In contrast, a new species of Ethmostigmus from north-east India formed a clade with African/Australian species. Fine-scale biogeographic analyses in PIP predict that Indian Ethmostigmus had an ancestor in southern-central parts of the Western Ghats. This was followed by four independent dispersal events from the southern-central Western Ghats to the Eastern Ghats, and between different parts of the Western Ghats in the Cenozoic.

CONCLUSIONS

Our results are consistent with Gondwanan break-up driving the early evolutionary history of the genus Ethmostigmus. Multiple dispersal events coinciding with geo-climatic events throughout the Cenozoic shaped diversification in PIP. Ethmostigmus species in PIP are restricted to wet forests and have retained that niche throughout their diversification.

Spatial Phylogenetics of Florida Vascular Plants: The Effects of Calibration and Uncertainty on Diversity Estimates

Recent availability of biodiversity data resources has enabled an unprecedented ability to estimate phylogenetically based biodiversity metrics over broad scales. Such approaches elucidate ecological and evolutionary processes yielding a biota and help guide conservation efforts. However, the choice of appropriate phylogenetic resources and underlying input data uncertainties may affect interpretation. Here, we address how differences among phylogenetic source trees and levels of phylogenetic uncertainty affect these metrics and test existing hypotheses regarding geographic biodiversity patterns across the diverse vascular plant flora of Florida, US. Ecological niche models for 1,490 Florida species were combined with a "purpose-built" phylogenetic tree (phylogram and chronogram), as well as with trees derived from community resources (Phylomatic and Open Tree of Life). There were only modest differences in phylodiversity metrics given the phylogenetic source tree and taking into account the level of phylogenetic uncertainty; we identify similar areas of conservation interest across Florida regardless of the method used.

A new species of the odorous frog genus Odorrana (Amphibia, Anura, Ranidae) from southwestern China

The genus Odorrana is widely distributed in the mountains of East and Southeastern Asia. An increasing number of new species in the genus have been recognized especially in the last decade. Phylogenetic studies of the O. schmackeri species complex with wide distributional range also revealed several cryptic species. Here, we describe a new species in the species complex from Guizhou Province of China. Phylogenetic analyses based on mitochondrial DNA indicated the new species as a monophyly clustered into the Odorrana clade and sister to O. schmackeri, and nuclear DNA also indicated it as an independent lineage separated from its related species. Morphologically, the new species can be distinguished from its congeners based on a combination of the following characters: (1) having smaller body size in males (snout-vent length (SVL) <43.3 mm); (2) head longer than wide; (3) dorsolateral folds absent; (4) tympanum of males large and distinct, tympanum diameter twice as long as width of distal phalanx of finger III; (5) two metacarpal tubercles; (6) relative finger lengths: II < I < IV < III; (7) tibiotarsal articulation reaching to the level between eye to nostril when leg stretched forward; (8) disks on digits with circum-marginal grooves; (9) toes fully webbed to disks; (10) the first subarticular tubercle on fingers weak; (11) having white pectoral spinules, paired subgular vocal sacs located at corners of throat, light yellow nuptial pad on the first finger in males.

Evolutionary origin of the latitudinal diversity gradient in liverworts

A latitudinal diversity gradient towards the tropics appears as one most recurrent patterns in ecology, but the mechanisms underlying this pattern remain an area of controversy. In angiosperms, the tropical conservatism hypothesis proposes that most groups originated in the tropics and are adapted to a tropical climatic regime, and that relatively few species have evolved physiological adaptations to cold, dry or unpredictable climates. This mechanism is, however, unlikely to apply across land plants, and in particular, to liverworts, a group of about 7500 species, whose ability to withstand cold much better than their tracheophyte counterparts is at odds with the tropical conservatism hypothesis. Molecular dating, diversification rate analyses and ancestral area reconstructions were employed to explore the evolutionary mechanisms that account for the latitudinal diversity gradient in liverworts. As opposed to angiosperms, tropical liverwort genera are not older than their extra-tropical counterparts (median stem age of tropical and extra-tropical liverwort genera of 24.35 ± 39.65 Ma and 39.57 ± 49.07 Ma, respectively), weakening the 'time for speciation hypothesis'. Models of ancestral area reconstructions with equal migration rates between tropical and extra-tropical regions outperformed models with asymmetrical migration rates in either direction. The symmetry and intensity of migrations between tropical and extra-tropical regions suggested by the lack of resolution in ancestral area reconstructions towards the deepest nodes are at odds with the tropical niche conservatism hypothesis. In turn, tropical genera exhibited significantly higher net diversification rates than extra-tropical ones, suggesting that the observed latitudinal diversity gradient results from either higher extinction rates in extra-tropical lineages or higher speciation rates in the tropics. We discuss a series of experiments to help deciphering the underlying evolutionary mechanisms.

Phylogenetic patterns of trait and trait plasticity evolution: Insights from amphibian embryos

Environmental variation favors the evolution of phenotypic plasticity. For many species, we understand the costs and benefits of different phenotypes, but we lack a broad understanding of how plastic traits evolve across large clades. Using identical experiments conducted across North America, we examined prey responses to predator cues. We quantified five life-history traits and the magnitude of their plasticity for 23 amphibian species/populations (spanning three families and five genera) when exposed to no cues, crushed-egg cues, and predatory crayfish cues. Embryonic responses varied considerably among species and phylogenetic signal was common among the traits, whereas phylogenetic signal was rare for trait plasticities. Among trait-evolution models, the Ornstein-Uhlenbeck (OU) model provided the best fit or was essentially tied with Brownian motion. Using the best fitting model, evolutionary rates for plasticities were higher than traits for three life-history traits and lower for two. These data suggest that the evolution of life-history traits in amphibian embryos is more constrained by a species' position in the phylogeny than is the evolution of life history plasticities. The fact that an OU model of trait evolution was often a good fit to patterns of trait variation may indicate adaptive optima for traits and their plasticities.

The impact of seasonality on niche breadth, distribution range and species richness: a theoretical exploration of Janzen's hypothesis

Being invoked as one of the candidate mechanisms for the latitudinal patterns in biodiversity, Janzen's hypothesis states that the limited seasonal temperature variation in the tropics generates greater temperature stratification across elevations, which makes tropical species adapted to narrower ranges of temperatures and have lower effective dispersal across elevations than species in temperate regions. Numerous empirical studies have documented latitudinal patterns in species elevational ranges and thermal niche breadths that are consistent with the hypothesis, but the theoretical underpinnings remain unclear. This study presents the first mathematical model to examine the evolutionary processes that could back up Janzen's hypothesis and assess the effectiveness of limited seasonal temperature variation to promote speciation along elevation in the tropics. Results suggest that trade-offs in thermal tolerances provide a mechanism for Janzen's hypothesis. Limited seasonal temperature variation promotes gradient speciation not due to the reduction in gene flow that is associated with narrow thermal niche, but due to the pleiotropic effects of more stable divergent selection of thermal tolerance on the evolution of reproductive incompatibility. The proposed modelling approach also provides a potential way to test a speciation model against genetic data.

Union of phylogeography and landscape genetics

Phylogeography and landscape genetics have arisen within the past 30 y. Phylogeography is said to be the bridge between population genetics and systematics, and landscape genetics the bridge between landscape ecology and population genetics. Both fields can be considered as simply the amalgamation of classic biogeography with genetics and genomics; however, they differ in the temporal, spatial, and organismal scales addressed and the methodology used. I begin by briefly summarizing the history and purview of each field and suggest that, even though landscape genetics is a younger field (coined in 2003) than phylogeography (coined in 1987), early studies by Dobzhansky on the "microgeographic races" of Linanthus parryae in the Mojave Desert of California and Drosophila pseudoobscura across the western United States presaged the fields by over 40 y. Recent advances in theory, models, and methods have allowed researchers to better synthesize ecological and evolutionary processes in their quest to answer some of the most basic questions in biology. I highlight a few of these novel studies and emphasize three major areas ripe for investigation using spatially explicit genomic-scale data: the biogeography of speciation, lineage divergence and species delimitation, and understanding adaptation through time and space. Examples of areas in need of study are highlighted, and I end by advocating a union of phylogeography and landscape genetics under the more general field: biogeography.

Evolutionary dynamics of an expressed MHC class IIβ locus in the Ranidae (Anura) uncovered by genome walking and high-throughput amplicon sequencing

The Major Histocompatibility Complex (MHC) is a genomic region encoding immune loci that are important and frequently used markers in studies of adaptive genetic variation and disease resistance. Given the primary role of infectious diseases in contributing to global amphibian declines, we characterized the hypervariable exon 2 and flanking introns of the MHC Class IIβ chain for 17 species of frogs in the Ranidae, a speciose and cosmopolitan family facing widespread pathogen infections and declines. We find high levels of genetic variation concentrated in the Peptide Binding Region (PBR) of the exon. Ten codons are under positive selection, nine of which are located in the mammal-defined PBR. We hypothesize that the tenth codon (residue 21) is an amphibian-specific PBR site that may be important in disease resistance. Trans-species and trans-generic polymorphisms are evident from exon-based genealogies, and co-phylogenetic analyses between intron, exon and mitochondrial based reconstructions reveal incongruent topologies, likely due to different locus histories. We developed two sets of barcoded adapters that reliably amplify a single and likely functional locus in all screened species using both 454 and Illumina based sequencing methods. These primers provide a resource for multiplexing and directly sequencing hundreds of samples in a single sequencing run, avoiding the labour and chimeric sequences associated with cloning, and enabling MHC population genetic analyses. Although the primers are currently limited to the 17 species we tested, these sequences and protocols provide a useful genetic resource and can serve as a starting point for future disease, adaptation and conservation studies across a range of anuran taxa.

Elevational diversity gradients of Tibetan loaches: The relative roles of ecological and evolutionary processes

It is widely believed that species richness patterns (SRPs) are shaped by both ecological and evolutionary processes. However, the relative roles of these processes remain unclear, especially for aquatic organisms. In this study, we integrated ecological and evolutionary measures to tease apart the relative influences of these factors on the SRP of Tibetan loaches along an extensive elevational gradient. We found that the Tibetan loaches displayed a richness pattern that peaked at midelevations. The mean annual temperature (MAT), mid‐domain effect (MDE), and summed age of colonization (SAC, complex of colonization age and colonization frequency) were the main drivers, accounting for 85%, 51%, and 88% of the variations in the SRP, respectively. The three predictors had very high combined effects (MAT‐MDE‐SAC, MAT‐SAC, and MDE‐SAC were 44%, 38%, and 6%, respectively). Our analyses suggested that energy input, time‐for‐speciation, and species dispersal may directly guide the SRP or mediate it by geometric constraints. Conclusively, the SRP of the Tibetan loaches with elevation is the outcome of interactions between biogeographical processes and regional ecological conditions.

Did true frogs 'dispersify'?

The interplay between range expansion and concomitant diversification is of fundamental interest to evolutionary biologists, particularly when linked to intercontinental dispersal and/or large scale extinctions. The evolutionary history of true frogs has been characterized by circumglobal range expansion. As a lineage that survived the Eocene-Oligocene extinction event (EOEE), the group provides an ideal system to test the prediction that range expansion triggers increased net diversification. We constructed the most densely sampled, time-calibrated phylogeny to date in order to: (i) characterize tempo and patterns of diversification; (ii) assess the impact of the EOEE; and (iii) test the hypothesis that range expansion was followed by increased net diversification. We show that late Eocene colonization of novel biogeographic regions was not affected by the EOEE and surprisingly, global expansion was not followed by increased net diversification. On the contrary, the diversification rate declined or did not shift following geographical expansion. Thus, the diversification history of true frogs contradicts the prevailing expectation that amphibian net diversification accelerated towards the present or increased following range expansion. Rather, our results demonstrate that despite their dynamic biogeographic history, true frogs diversified at a relatively constantly rate, even as they colonized the major land masses of Earth.

Molecular and Morphological Study of Leaping Frogs (Anura, Ranixalidae) with Description of Two New Species

The monotypic anuran family Ranixalidae is endemic to India, with a predominant distribution in the Western Ghats, a region that is home to several unique amphibian lineages. It is also one of the three ancient anuran families that diversified on the Indian landmass long before several larger radiations of extant frogs in this region. In recent years, ranixalids have been subjected to DNA barcoding and systematic studies. Nearly half of the presently recognized species in this family have been described over the last three years, along with recognition of a new genus to accommodate three previously known members. Our surveys in the Western Ghats further suggest the presence of undescribed diversity in this group, thereby increasing former diversity estimates. Based on rapid genetic identification using a mitochondrial gene, followed by phylogenetic analyses with an additional nuclear gene and detailed morphological studies including examination of museum specimens, new collections, and available literature, here we describe two new species belonging to the genus Indirana from the Western Ghats states of Karnataka and Kerala. We also provide new genetic and morphological data along with confirmed distribution records for all the species known prior to this study. This updated systematic revision of family Ranixalidae will facilitate future studies and provide vital information for conservation assessment of these relic frogs.

Latitudinal diversity gradients in Mesozoic non-marine turtles

The latitudinal biodiversity gradient (LBG)-the pattern of increasing taxonomic richness with decreasing latitude-is prevalent in the structure of the modern biota. However, some freshwater taxa show peak richness at mid-latitudes; for example, extant Testudines (turtles, terrapins and tortoises) exhibit their greatest diversity at 25° N, a pattern sometimes attributed to recent bursts of climatically mediated species diversification. Here, we test whether this pattern also characterizes the Mesozoic distribution of turtles, to determine whether it was established during either their initial diversification or as a more modern phenomenon. Using global occurrence data for non-marine testudinate genera, we find that subsampled richness peaks at palaeolatitudes of 15-30° N in the Jurassic, 30-45° N through the Cretaceous to the Campanian, and from 30° to 60° N in the Maastrichtian. The absence of a significant diversity peak in southern latitudes is consistent with results from climatic models and turtle niche modelling that demonstrate a dearth of suitable turtle habitat in Gondwana during the Jurassic and Late Cretaceous. Our analyses confirm that the modern testudinate LBG has a deep-time origin and further demonstrate that LBGs are not always expressed as a smooth, equator-to-pole distribution.

The amphibians and reptiles of Mindanao Island, southern Philippines, II: the herpetofauna of northeast Mindanao and adjacent islands

We summarize all available amphibian and reptile species distribution data from the northeast Mindanao faunal region, including small islands associated with this subcenter of endemic vertebrate biodiversity. Together with all publicly available historical information from biodiversity repositories, we present new data from several major herpetological surveys, including recently conducted inventories on four major mountains of northeast Mindanao, and adjacent islands of Camiguin Sur, Dinagat, and Siargao. We present species accounts for all taxa, comment on unresolved taxonomic problems, and provide revisions to outdated IUCN conservation status assessments in cases where our new data significantly alter earlier classification status summaries. Together, our comprehensive analysis of this fauna suggests that the greater Mindanao faunal region possesses distinct subcenters of amphibian and reptile species diversity, and that until this area is revisited and its fauna and actually studied, with on-the-ground field work including targeted surveys of species distributions coupled to the study their natural history, our understanding of the diversity and conservation status of southern Philippine herpetological fauna will remain incomplete. Nevertheless, the northeast Mindanao geographical area (Caraga Region) appears to have the highest herpetological species diversity (at least 126 species) of any comparably-sized Philippine faunal subregion.

Phylogenetic niche conservatism explains an inverse latitudinal diversity gradient in freshwater arthropods

The underlying mechanisms responsible for the general increase in species richness from temperate regions to the tropics remain equivocal. Many hypotheses have been proposed to explain this astonishing pattern but additional empirical studies are needed to shed light on the drivers at work. Here we reconstruct the evolutionary history of the cosmopolitan diving beetle subfamily Colymbetinae, the majority of which are found in the Northern hemisphere, hence exhibiting an inversed latitudinal diversity gradient. We reconstructed a dated phylogeny using 12 genes, to investigate the biogeographical history and diversification dynamics in the Colymbetinae. We aimed to identify the role that phylogenetic niche conservatism plays in the inversed diversification pattern seen in this group. Our results suggest that Colymbetinae originated in temperate climates, which supports the hypothesis that their distribution is the result of an ancestral adaptation to temperate environmental conditions rather than tropical origins, and that temperate niche conservatism can generate and/or maintain inverse latitudinal diversity gradients.

Phylogeny of the Oniticellini and Onthophagini dung beetles (Scarabaeidae, Scarabaeinae) from morphological evidence

A phylogenetic study was conducted to hypothesize relationships of most of the genera of the Oniticellini and Onthophagini for the first time using morphological characters from a diverse array of external and internal sclerites. The monophyly and sister relationship of both tribes was found using Bayesian and parsimony analyses with heavily to moderately weighted data. An alternative hypothesis based on parsimony analyses of unweighted or slightly weighted data show a paraphyletic Oniticellini without the Onthophagini, although recognition of the subtribe Helictopleurina as a tribe would eliminate non-monophyly. Of the three Oniticellini subtribes, the Helictopleurina and Drepanocerina are monophyletic. There is no support for the monophyly of the Oniticellina or the Onthophagini subtribe Alloscelina, as currently defined. The genus Liatongus is paraphyletic, while strong support was found for monophyly of the Madagascan genus, Helictopleurus. The genus Onthophagus is never monophyletic in any analysis performed. Two new subtribes are also proposed: Liatongina subtr. n. including the genus Liatongus and Attavicina subtr. n. including the genera Attavicinus and Paroniticellus. Topological evidence shows that the ancestral oniticellines and onthophagines were all coprophagous with alternative food sources evolving relatively recently. Both myrmecophily and termitophily probably evolved only once in the onthophagines. The phylogenetic analysis supports an African origin for the two tribes, with a relatively early age for the split of the Madagascar helictopleurines from the remaining oniticellines via dispersal. Furthermore, the presence of the oniticellines in the New World is hypothesized to be due to two relatively old dispersal events via Beringia and two relatively recent trans-Atlantic invasions of the Caribbean.

Environmental constraints and call evolution in torrent-dwelling frogs

Although acoustic signals are important for communication in many taxa, signal propagation is affected by environmental properties. Strong environmental constraints should drive call evolution, favoring signals with greater transmission distance and content integrity in a given calling habitat. Yet, few empirical studies have verified this prediction, possibly due to a shortcoming in habitat characterization, which is often too broad. Here we assess the potential impact of environmental constraints on the evolution of advertisement call in four groups of torrent-dwelling frogs in the family Ranidae. We reconstruct the evolution of calling site preferences, both broadly categorized and at a finer scale, onto a phylogenetic tree for 148 species with five markers (∼3600 bp). We test models of evolution for six call traits for 79 species with regard to the reconstructed history of calling site preferences and estimate their ancestral states. We find that in spite of existing morphological constraints, vocalizations of torrent-dwelling species are most probably constrained by the acoustic specificities of torrent habitats and particularly their high level of ambient noise. We also show that a fine-scale characterization of calling sites allows a better perception of the impact of environmental constraints on call evolution.

The Australasian frog family Ceratobatrachidae in China, Myanmar and Thailand: discovery of a new Himalayan forest frog clade

In an effort to study the systematic affinities and specieslevel phylogenetic relationships of the enigmatic anurans variably assigned to the genera Ingerana or Limnonectes (family Dicroglossidae), we collected new molecular sequence data for five species including four Himalayan taxa, Limnonectes xizangensis, Lim. medogensis, Lim. alpina, [corrected] Ingerana borealis and one southeast Asian species, I. tasanae, and analyzed these together with data from previous studies involving other ostensibly related taxa. Our surprising results demonstrate unequivocally that Lim. xizangensis, Lim. medogensis and Lim. alpina [corrected] form a strongly supported clade, the sister-group of the family Australasian forest frog family Ceratobatrachidae. This discovery requires an expansion of the definition of Ceratobatrachidae and represents the first record of this family in China. These three species are distinguished from the species of Ingerana and Limnonectes by the: (1) absence of interdigital webbing of the foot, (2) absence of terminal discs on fingers and toes, (3) absence of circumarginal grooves on the fingers and toes, and (4) absence of tarsal folds. Given their phylogenetic and morphological distinctiveness, we assign them to the oldest available generic name for this clade, Liurana Dubois 1987, and transfer Liurana from Dicroglossidae to the family Ceratobatrachidae. In contrast, Ingerana tasanae was found to be clustered with strong support with the recently described genus Alcalus (Ceratobatrachidae), a small clade of otherwise Sundaic species; this constitutes a new record of the family Ceratobatrachidae for Myanmar and Thailand. Finally, Ingerana borealis clustered with the "true" Ingerana (family Dicroglossidae), for which the type species is I. tenasserimensis.