Comparative phylogeography of direct-developing frogs (Anura: Craugastoridae: Pristimantis) in the southern Andes of Colombia

Unveiling the evolutionary relationships and the high cryptic diversity in Andean rainfrogs (Craugastoridae: Pristimantis myersi group)

Background

Pristimantis is the most diverse genus of terrestrial frogs. Historically, it has been divided into several phenetic groups in order to facilitate species identification. However, in light of phylogenetic analysis, many of these groups have been shown to be non-monophyletic, denoting a high degree of morphological convergence and limited number of diagnostic traits. In this study, we focus on the Pristimantis myersi group, an assemblage of small rainfrogs distributed throughout the Andes of Ecuador and Colombia, whose external morphology is highly conserved, and its species diversity and evolutionary relationships largely unknown.

Methods

We inferred a new phylogenetic hypothesis for the frog genus Pristimantis, including all available sequences of the mtDNA 16S rRNA, as well as new DNA sequences from 175 specimens. Our sampling included 19 of the 24 species currently recognized as part of the Pristimantis myersi group.

Results

Our new evolutionary hypothesis recovered the P. myersi group as non-monophyletic and composed of 16 species. Therefore, we exclude P. albujai, P. bicantus, P. sambalan, and P. nelsongalloi in order to preserve the monophyly of the group. We discovered at least eight candidate species, most of them hidden under the names of P. leoni, P. hectus, P. festae, P. gladiator, and P. ocreatus.

Discussion

Our results reveal the occurrence of a high level of cryptic diversity to the species level within the P. myersi group and highlight the need to redefine some of its species and reassess their conservation status. We suggest that the conservation status of six species within the group need to be re-evaluated because they exhibit smaller distributions than previously thought; these species are: P. festae, P. gladiator, P. hectus, P. leoni, P. ocreatus, and P. pyrrhomerus. Finally, given that the Pristimantis myersi group, as defined in this work, is monophyletic and morphologically diagnosable, and that Trachyphrynus is an available name for the clade containing P. myersi, we implement Trachyphrynus as a formal subgenus name for the Pristimantis myersi group.

Six new species of Pristimantis (Anura: Strabomantidae) from Llanganates National Park and Sangay National Park in Amazonian cloud forests of Ecuador

We describe six new species of rainfrogs of the genus Pristimantis (Strabomantidae) from Amazonian cloud forests in Ecuador. We also present a phylogeny showing the relationships of the new species. The phylogeny is based on mitochondrial genes 16S rRNA (16S), 12 rRNA (12S), NADH-ubiquinone oxidoreductase chain 1 (ND1) and the nuclear gene recombination-activating 1 (RAG1). We also describe the osteology of two of the new species using high-resolution x-ray computed tomography. The new species belong to two clades. The first clade is sister to the subgenus Huicundomantis and includes P. tamia sp. nov., P. miktos, and P. mallii. Pristimantis tamia sp. nov. is morphologically similar to P. miktos, P. mallii, P. martiae, and P. incomptus, but differs from them by lacking vocal slits and tympanic membrane and by having light greenish blue iris. Based in our results we expand the subgenus Huicundomantis to include the P. miktos species group. The second clade is remarkable by being highly divergent and consisting exclusively of new species: P. anaiae sp. nov., P. glendae sp. nov., P. kunam sp. nov., P. resistencia sp. nov., and P. venegasi sp. nov. The new species resemble P. roni, P. yanezi, P. llanganati, P. katoptroides, P. verecundus, and P. mutabilis but can be distinguished from them by lacking vocal slits and tympanic membrane and by having large dark round areas with thin clear borders in the sacral region. All six new species occur in the eastern slopes of the Ecuadorian Andes and are known from a single locality in Llanganates or Sangay National Park. We recommend assigning all of them to the Data Deficient (DD) Red List category. Based in our high-resolution x-ray tomographies, we report the presence of structures that appear to be intercalary elements. This would be the first report of such structures in Terrarana.

The Diversity and Evolution of Sex Chromosomes in Frogs

Frogs are ideal organisms for studying sex chromosome evolution because of their diversity in sex chromosome differentiation and sex-determination systems. We review 222 anuran frogs, spanning ~220 Myr of divergence, with characterized sex chromosomes, and discuss their evolution, phylogenetic distribution and transitions between homomorphic and heteromorphic states, as well as between sex-determination systems. Most (~75%) anurans have homomorphic sex chromosomes, with XY systems being three times more common than ZW systems. Most remaining anurans (~25%) have heteromorphic sex chromosomes, with XY and ZW systems almost equally represented. There are Y-autosome fusions in 11 species, and no W-/Z-/X-autosome fusions are known. The phylogeny represents at least 19 transitions between sex-determination systems and at least 16 cases of independent evolution of heteromorphic sex chromosomes from homomorphy, the likely ancestral state. Five lineages mostly have heteromorphic sex chromosomes, which might have evolved due to demographic and sexual selection attributes of those lineages. Males do not recombine over most of their genome, regardless of which is the heterogametic sex. Nevertheless, telomere-restricted recombination between ZW chromosomes has evolved at least once. More comparative genomic studies are needed to understand the evolutionary trajectories of sex chromosomes among frog lineages, especially in the ZW systems.

A new Terrarana frog of genus Pristimantis from an unexplored cloud forest from the eastern Andes, Colombia

A new species of Pristimantis (Craugastoridae, subgenusPristimantis) is described from a relict and unexplored cloud forest in the western slope from Cordillera Oriental of the Colombian Andes. The specific name was chosen by consensus expert scientists and local people. Pristimantischamezensissp. nov. is easily distinguished from congeneric species by having a gray iris with black reticulations in life, subconical tubercles on the upper eyelid, the chin edged with irregular, dark-brown blotches, and conical heel tubercles. The phylogenetic analyses suggest that the origin and radiation of its clade may have occurred in the highlands. With the description of P.chamezensissp. nov., we identify 14 species distributed throughout the eastern slope of the Andes that are associated with the Orinoco Basin.

The enigmatic Leiosaurae clade: Phylogeography, species delimitation, phylogeny and historical biogeography of its southernmost species

The clade Leiosaurae is composed of poorly-known species endemic to the southern region of South America. The difficulties of finding these lizards in the field, and their highly conserved morphology, have limited our taxonomic knowledge and understanding of their evolutionary histories. Here, we use data collected over 9 years to study the phylogenetic history, genetic diversity, and biogeographic history of almost all the southernmost species of Leiosaurae (except P. nigroigulus), including: Leiosaurus bellii, Diplolaemus darwinii, D. bibronii, D. sexcinctus and D. leopardinus. We use a fragment of the mitochondrial cytochrome-b gene to resolve general phylogeographic patterns, and add another mitochondrial gene and eight nuclear genes to perform species delimitation and phylogenetic analyses associated with divergence times. We found evidence for three putative new species-level taxa within L. bellii and five within Diplolaemus species, indicating high levels of geographic structure. We use a time-calibrated phylogeny to estimate ranges of ancestral distributions and to generate new hypotheses about their historical biogeography.

A new species of terrestrial frog Pristimantis (Strabomantidae) from the upper basin of the Pastaza River, Ecuador

We describe a new species of Pristimantis from the montane forest of the Río Zuñag Ecological Reserve, upper basin of the Pastaza River, Ecuador. Pristimantismalliisp. n. is characterized by a snout-vent length of 11.6–21.3 mm in adult males (n = 12), 22.6–34.3 mm in adult females (n = 8), and is compared morphologically and genetically with Pristimantismiktos and with other relevant species of Pristimantis. The new species is characterized by having skin on dorsum and flanks shagreen, distinctive scapular folds, snout broadly rounded in dorsal view, upper eyelid bearing one or two subconical tubercles and some rounded tubercles, dorsum and flanks light brown to brown, with irregular dark brown marks bounded by dirty cream and groin with irregular yellowish marks.

Gene flow and Andean uplift shape the diversification of Gasteracantha cancriformis (Araneae: Araneidae) in Northern South America

The Andean uplift has played a major role in shaping the current Neotropical biodiversity. However, in arthropods other than butterflies, little is known about how this geographic barrier has impacted species historical diversification. Here, we examined the phylogeography of the widespread color polymorphic spider Gasteracantha cancriformis to evaluate the effect of the northern Andean uplift on its divergence and assess whether its diversification occurred in the presence of gene flow. We inferred phylogenetic relationships and divergence times in G. cancriformis using mitochondrial and nuclear data from 105 individuals in northern South America. Genetic diversity, divergence, and population structure were quantified. We also compared multiple demographic scenarios for this species using a model-based approach (phrapl) to determine divergence with or without gene flow. At last, we evaluated the association between genetic variation and color polymorphism. Both nuclear and mitochondrial data supported two well-differentiated clades, which correspond to populations occurring on opposite sides of the Eastern cordillera of the Colombian Andes. The final uplift of this cordillera was identified as the most likely force that shaped the diversification of G. cancriformis in northern South America, resulting in a cis- and trans-Andean phylogeographic structure for the species. We also found shared genetic variation between the cis- and trans-Andean clades, which is better explained by a scenario of historical divergence in the face of gene flow. This has been likely facilitated by the presence of low-elevation passes across the Eastern Colombian cordillera. Our work constitutes the first example in which the Andean uplift coupled with gene flow influenced the evolutionary history of an arachnid lineage.

A new species of Pristimantis from eastern Brazilian Amazonia (Anura, Craugastoridae)

In this study a new species of Pristimantis (Anura: Craugastoridae) of the P. conspicillatus species group is described. Pristimantis latrosp. n. is known only from the municipalities of Altamira, Anapu, Brasil Novo, Medicilândia, Uruará and Aveiro (Flona Tapajós, right bank of Tapajós river), in Pará state, Brazil. Morphologically, the new species distinguishes from known congeners in the group mainly by the presence of dorsal tubercles and absence of discoidal folds, smooth belly skin, as well as the presence of supernumerary tubercles on hands. The call of the new species consists of seven ascending notes, the first of which has a dominant frequency of 2635 Hz and the last 3272 Hz. Molecular analysis of the 16S mtDNA indicates a genetic distance of 8% to P. chiastonotus, its closet relative, and between 9% and 11% to populations of P. fenestratus.

Diversification of the rainfrog Pristimantis ornatissimus in the lowlands and Andean foothills of Ecuador

Geographic barriers and elevational gradients have long been recognized as important in species diversification. Here, we illustrate an example where both mechanisms have shaped the genetic structure of the Neotropical rainfrog, Pristimantis ornatissimus, which has also resulted in speciation. This species was thought to be a single evolutionary lineage distributed throughout the Ecuadorian Chocó and the adjacent foothills of the Andes. Based on recent sampling of P. ornatissimus sensu lato, we provide molecular and morphological evidence that support the validity of a new species, which we name Pristimantis ecuadorensis sp. nov. The sister species are elevational replacements of each other; the distribution of Pristimantis ornatissimus sensu stricto is limited to the Ecuadorian Chocó ecoregion (< 1100 m), whereas the new species has only been found at Andean localities between 1450-1480 m. Given the results of the Multiple Matrix Regression with Randomization analysis, the genetic difference between P. ecuadorensis and P. ornatissimus is not explained by geographic distance nor environment, although environmental variables at a finer scale need to be tested. Therefore this speciation event might be the byproduct of stochastic historic extinction of connected populations or biogeographic events caused by barriers to dispersal such as rivers. Within P. ornatissimus sensu stricto, morphological patterns and genetic structure seem to be related to geographic isolation (e.g., rivers). Finally, we provide an updated phylogeny for the genus, including the new species, as well as other Ecuadorian Pristimantis.

Comparative Phylogeography of Ethiopian anurans: impact of the Great Rift Valley and Pleistocene climate change

Background

The Ethiopian highlands are a biodiversity hotspot, split by the Great Rift Valley into two distinct systems of plateaus and mountains. The Rift Valley is currently hot and dry and acts as a barrier to gene flow for highland-adapted species. It is however unlikely that the conditions in the Rift were inhospitable to highland species during the entire Pleistocene. To assess the significance of the Ethiopian Rift as a biogeographic barrier as well as the impact Pleistocene climatic changes have had on the evolution of Ethiopian organisms, we performed phylogeographic analyses and developed present and past niche models on seven anuran species with different elevational and ecological preferences.

Results

We found that highland species on the east and the west sides of the Rift are genetically differentiated and have not experienced any detectable gene flow for at least 0.4 my. In contrast, species found at elevations lower than 2500 m do not show any population structure. We also determined that highland species have lower effective population sizes than lowland species, which have experienced a large, yet gradual, demographic expansion, starting approximately half a million year ago.

Conclusions

The pattern we report here is consistent with the increasingly warmer and drier conditions of the Pleistocene in East Africa, which resulted in the expansion of savanna, the fragmentation of forests and the shrinking of highland habitats. Climatic niche models indicated that the Rift is currently non suitable for most of the studied species, but it could have been a more permeable barrier during the Last Glacial Maximum. However, considering the strong genetic structure of highland species, we hypothesize that the barrier mechanisms at the Rift are not only climatic but also topographical.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0774-1) contains supplementary material, which is available to authorized users.

Comparative Phylogeography Reveals Cryptic Diversity and Repeated Patterns of Cladogenesis for Amphibians and Reptiles in Northwestern Ecuador

Comparative phylogeography allow us to understand how shared historical circumstances have shaped the formation of lineages, by examining a broad spectrum of co-distributed populations of different taxa. However, these types of studies are scarce in the Neotropics, a region that is characterized by high diversity, complex geology, and poorly understood biogeography. Here, we investigate the diversification patterns of five lineages of amphibians and reptiles, co-distributed across the Choco and Andes ecoregions in northwestern Ecuador. Mitochondrial DNA and occurrence records were used to determine the degree of geographic genetic divergence within species. Our results highlight congruent patterns of parapatric speciation and common geographical barriers for distantly related taxa. These comparisons indicate similar biological and demographic characteristics for the included clades, and reveal the existence of two new species of Pristimantis previously subsumed under P. walkeri, which we describe herein. Our data supports the hypothesis that widely distributed Chocoan taxa may generally experience their greatest opportunities for isolation and parapatric speciation across thermal elevational gradients. Finally, our study provides critical information to predict which unstudied lineages may harbor cryptic diversity, and how geology and climate are likely to have shaped their evolutionary history.

Spatial diversity patterns of Pristimantis frogs in the Tropical Andes

Although biodiversity gradients have been widely documented, the factors governing broad-scale patterns in species richness are still a source of intense debate and interest in ecology, evolution, and conservation biology. Here, we tested whether spatial hypotheses (species-area effect, topographic heterogeneity, mid-domain null model, and latitudinal effect) explain the pattern of diversity observed along the altitudinal gradient of Andean rain frogs of the genus Pristimantis. We compiled a gamma-diversity database of 378 species of Pristimantis from the tropical Andes, specifically from Colombia to Bolivia, using records collected above 500 m.a.s.l. Analyses were performed at three spatial levels: Tropical Andes as a whole, split in its two main domains (Northern and Central Andes), and split in its 11 main mountain ranges. Species richness, area, and topographic heterogeneity were calculated for each 500-m-width elevational band. Spatial hypotheses were tested using linear regression models. We examined the fit of the observed diversity to the mid-domain hypothesis using randomizations. The species richness of Pristimantis showed a hump-shaped pattern across most of the altitudinal gradients of the Tropical Andes. There was high variability in the relationship between area and species richness along the Tropical Andes. Correcting for area effects had little impact in the shape of the empirical pattern of biodiversity curves. Mid-domain models produced similar gradients in species richness relative to empirical gradients, but the fit varied among mountain ranges. The effect of topographic heterogeneity on species richness varied among mountain ranges. There was a significant negative relationship between latitude and species richness. Our findings suggest that spatial processes partially explain the richness patterns of Pristimantis frogs along the Tropical Andes. Explaining the current patterns of biodiversity in this hot spot may require further studies on other possible underlying mechanisms (e.g., historical, biotic, or climatic hypotheses) to elucidate the factors that limit the ranges of species along this elevational gradient.

Co-Infection by Chytrid Fungus and Ranaviruses in Wild and Harvested Frogs in the Tropical Andes

While global amphibian declines are associated with the spread of Batrachochytrium dendrobatidis (Bd), undetected concurrent co-infection by other pathogens may be little recognized threats to amphibians. Emerging viruses in the genus Ranavirus (Rv) also cause die-offs of amphibians and other ectotherms, but the extent of their distribution globally, or how co-infections with Bd impact amphibians are poorly understood. We provide the first report of Bd and Rv co-infection in South America, and the first report of Rv infections in the amphibian biodiversity hotspot of the Peruvian Andes, where Bd is associated with extinctions. Using these data, we tested the hypothesis that Bd or Rv parasites facilitate co-infection, as assessed by parasite abundance or infection intensity within individual adult frogs. Co-infection occurred in 30% of stream-dwelling frogs; 65% were infected by Bd and 40% by Rv. Among terrestrial, direct-developing Pristimantis frogs 40% were infected by Bd, 35% by Rv, and 20% co-infected. In Telmatobius frogs harvested for the live-trade 49% were co-infected, 92% were infected by Bd, and 53% by Rv. Median Bd and Rv loads were similar in both wild (Bd = 10^1.2 Ze, Rv = 10^2.3 viral copies) and harvested frogs (Bd = 10^3.1 Ze, Rv = 10^2.7 viral copies). While neither parasite abundance nor infection intensity were associated with co-infection patterns in adults, these data did not include the most susceptible larval and metamorphic life stages. These findings suggest Rv distribution is global and that co-infection among these parasites may be common. These results raise conservation concerns, but greater testing is necessary to determine if parasite interactions increase amphibian vulnerability to secondary infections across differing life stages, and constitute a previously undetected threat to declining populations. Greater surveillance of parasite interactions may increase our capacity to contain and mitigate the impacts of these and other wildlife diseases.

Insights from Integrative Systematics Reveal Cryptic Diversity in Pristimantis Frogs (Anura: Craugastoridae) from the Upper Amazon Basin

Pluralistic approaches to taxonomy facilitate a more complete appraisal of biodiversity, especially the diversification of cryptic species. Although species delimitation has traditionally been based primarily on morphological differences, the integration of new methods allows diverse lines of evidence to solve the problem. Robber frogs (Pristimantis) are exemplary, as many of the species show high morphological variation within populations, but few traits that are diagnostic of species. We used a combination of DNA sequences from three mitochondrial genes, morphometric data, and comparisons of ecological niche models (ENMs) to infer a phylogenetic hypothesis for the Pristimantis acuminatus complex. Molecular phylogenetic analyses revealed a close relationship between three new species—Pristimantis enigmaticus sp. nov., P. limoncochensis sp. nov. and P. omeviridis sp. nov.—originally confused with Pristimantis acuminatus. In combination with morphometric data and geographic distributions, several morphological characters such as degree of tympanum exposure, skin texture, ulnar/tarsal tubercles and sexual secondary characters (vocal slits and nuptial pads in males) were found to be useful for diagnosing species in the complex. Multivariate discriminant analyses provided a successful classification rate for 83–100% of specimens. Discriminant analysis of localities in environmental niche space showed a successful classification rate of 75–98%. Identity tests of ENMs rejected hypotheses of niche equivalency, although not strongly because the high values on niche overlap. Pristimantis acuminatus and P. enigmaticus sp. nov. are distributed along the lowlands of central–southern Ecuador and northern Peru, in contrast with P. limoncochensis sp. nov. and P. omeviridis sp. nov., which are found in northern Ecuador and southern Colombia, up to 1200 m in the upper Amazon Basin. The methods used herein provide an integrated framework for inventorying the greatly underestimated biodiversity in Amazonia.

A likelihood inference of historical biogeography in the world's most diverse terrestrial vertebrate genus: diversification of direct-developing frogs (Craugastoridae: Pristimantis) across the Neotropics

The geology of the northern Andean region has driven the evolutionary history of Neotropical fauna through the creation of barriers and connections that have resulted in speciation and dispersal events, respectively. One of the most conspicuous groups of anuran fauna in the Andes and surrounding areas is the direct-developing species of the genus Pristimantis. We investigated the molecular phylogenetic placement of 12 species from the montane Andes of Colombia in a broader geographical context with a new genus-level phylogeny in order to identify the role of Andean orogeny over the last 40million years and the effect of elevational differences in diversification of Pristimantis. We examined the biogeographic history of the genus using ancestral range reconstruction by biogeographic regions and elevational ranges. We recognized the middle elevational band (between 1000 and 3000m) in the Northwestern Andes region of Colombia and Ecuador as a focal point for the origin and radiation of Pristimantis species. Additionally, we found several Andean migrations toward new habitats in Central Andes and Merida Andes for some species groups. We suggest that the paleogeological changes in the Northwestern Andes were the main promoter of speciation in Pristimantis, and may have served as a corridor for the dispersion of lowland species.