Divergence at the edges: peripatric isolation in the montane spiny throated reed frog complex
© Lawson et al. 2015
Received: 12 February 2015
Accepted: 20 May 2015
Published: 1 July 2015
Peripatric speciation and peripheral isolation have uncertain importance in species accumulation, and are largely overshadowed by assumed dominance of allopatric modes of speciation. Understanding the role of different speciation mechanisms within biodiversity hotspots is central to understanding the generation of biological diversity. Here, we use a phylogeographic analysis of the spiny-throated reed frogs and examine sister pairings with unbalanced current distributional ranges for characteristics of peripatric speciation. We further investigate whether forest/grassland mosaic adapted species are more likely created through peripatric speciation due to instability of this habitat type.
We reconstructed a multi-locus molecular phylogeny of spiny-throated reed frogs which we then combined with comparative morphologic data to delimit species and analyze historical demographic change; identifying three new species. Three potential peripatric speciation events were identified along with one case of allopatric speciation. Peripatric speciation is supported through uneven potential and realized distributions and uneven population size estimates based on field collections. An associated climate shift was observed in most potentially peripatric splits. Morphological variation was highest in sexually dimorphic traits such as body size and gular shape, but this variation was not limited to peripatric species pairs as hypothesized. The potentially allopatric species pair showed no niche shifts and equivalent effective population sizes, ruling out peripatry in that speciation event. Two major ecological niche shifts were recovered within this radiation, possibly as adaptations to occupy areas of grassland that became more prevalent in the last 5 million years. Restricted and fluctuating grassland mosaics within forests might promote peripatric speciation in the Eastern Arc Biodiversity Hotspot (EABH).
In our case study, peripatric speciation appears to be an important driver of diversity within the EABH biodiversity hotspot, implying it could be a significant speciation mechanism in highly fragmented ecosystems. Extensive peripatric speciation in this montane archipelago may explain the abundance of discrete lineages within the limited area of the EABH, as inferred in remote island archipelagos. Future phylogenetic studies incorporating demographic and spatial analyses will clarify the role of peripatric speciation in creating biodiversity hotspots.
Isolation has been historically considered the most frequent mechanism to limit gene flow, promoting divergence of lineages leading to speciation [1, 2]. The two dominant proposed forms of isolation have been vicariant allopatric (dichopatric) and peripatric, while sympatric and parapatric speciation appear less common  (see below). The classical model of vicariant allopatry starts with an ancestral distribution, which is split into large, disjointed populations ; hereafter referred to as “vicariant allopatry” or simply “allopatry”. Numerous examples have been cited demonstrating the process of allopatry in driving speciation (summarized in [3, 4]), and the signal of vicariant barriers limiting gene flow has been shown in many comparative phylogeographic studies [5–8]. In peripatric speciation, a small population is isolated at the edge of a larger population and rapidly diversifies, often crossing a previous barrier to dispersal and/or occupying a newly available ecological niche/habitat type. Over short periods of evolutionary time, peripatric events leave different evolutionary footprints from allopatry in terms of range, historical demography, and degree of ecological and morphological divergence . Parapatric speciation events also involve a shift in habitat requirements or ecological niche, though unlike in peripatric speciation, there are no physical barriers separating divergent lineages. Sympatric speciation, where species diverge within the same environment without spatial or ecological barriers, has rarely if ever been observed . New analytical approaches using multiple types of data enable testing for these processes in greater detail [10–13].
Different speciation processes produce predictable outcomes in the molecular, morphological, and geographic relationships among lineages. In contrast to the stable population sizes associated with allopatric speciation, a peripatric speciation event yields genetic variation consistent with bottlenecks and low founder population numbers . Peripatric species are also more likely to be ecologically divergent, as they have crossed a barrier to dispersal that previously limited occupation of new habitat at the edge of the ancestral range [3, 15]. Allopatric species, in contrast, more often retain their ecological niche due to the similarity of habitats isolated by vicariant barriers [16–18]. Morphologically, divergence in physical traits, especially those shaped by sexual selection, are predicted to be especially affected in peripatric populations both due to drift in small populations and as reinforcement between sister lineages in close proximity [19–22]. Quantifying molecular, ecological, and morphological metrics among sister lineages can therefore distinguish between these two modes (allopatry and peripatry) in species formation.
Across the African continent, speciation patterns are largely consistent with fragmentation of formally continuous forests through fluctuating climate [23–28]. Such patterns occur in the Eastern Afromontane Biodiversity Hotspot (EABH). Here, isolation of forest endemics on refugial mountain tops appears to be a common driver of exceptional levels of endemism and species diversity (e.g., [26, 29–31]). However, in part due to incomplete phylogenetic and phylogeographic sampling and imprecise distribution models, little attention has been paid to whether regional-scale fragmentation (vicariance/allopatric) or local-scale expansions or niche shifts (peripatric) are the dominant speciation process within this sky-island system.
We evaluate an insular montane radiation of spiny-throated reed frogs within the EABH and assess regional-scale and local-scale speciation mechanisms. Deciphering these two processes for individual radiations is imperative to understand diversification in tropical ecosystems, as both niche conservatism and niche divergence appear to play a major role in amphibian speciation [32, 33]. Establishing the speciation processes involved in the spiny-throated reed frog radiation allows for direct comparisons with other EABH endemic radiations as well as linking these processes to montane and island archipelagos around the world.
Within the EABH, the arrangement of montane areas of varying size and distance from each other provides an ideal setting in which to compare mechanisms of divergence. The mountains of the EABH system are dominated today by two forest vegetation types: rainforest found approximately between 700-1800 m or lower depending upon the mountain block and aspect (referred to as “rainforest” hereafter), and a montane forest/grassland mosaic (referred to as “mosaic” hereafter) above ca. 1800 m in the southern EABH . The extent and distribution of these biomes throughout the EABH has varied through time, with a gradual shift towards increased grasslands from 4 mya onward . The areas of high montane forest/grassland mosaic have historically expanded and contracted and at times likely vanished or become restricted to extremely narrow refugia as forests expanded . Lower elevation rainforest habitats fluctuated to a lesser extent, and were likely persistently stable on montane slopes despite habitat cycling in the adjacent lowland savannahs (classically outlined by ; see [38–41]). The larger fluctuations, smaller potential distributions, and the relative newness of mosaic habitat compared to mid-elevation rainforests, may make the mosaic habitat biome more susceptible to peripatry mechanisms than the more persistent rainforest biome. Though other biome-specific patterns of diversification are known [42, 43], whether peripatric speciation could also be biome-specific is less well documented .
Summary statistics of molecular loci used in this study
Summary statistics for each locus
Combined nuclear loci
Number of samples
Theta S sd
Theta Pi sd
P Taj D
Number of samples
Theta S sd
Theta Pi sd
P Taj D
Species Delimitations from BFD *BEAST stepping stone and pathsampling and BPP
*BEAST stepping stone and pathsampling results from BFD
Species groups by locality
All tanzanian spinigularis
Eastern Arc spinigularis
BPP species Delimitation from 10,000 iterations
# of species
Species groups by locality
Effective population size, population bottlenecks and expansions
Species, altitudinal range, habitat and available area of occurrence from IUCN estimates
East Usambara: 900–1100 m
Nguru: 900–1000 m
Uluguru: 980 m
Livingstone: 2010 m
Montane forest edge
Njombe: 2010 m
Montane forest edge and grassland
Udzungwa: 1680–1970 m
Malawi: 690 m
Submontane forest and forest edge
Mozambique: 1250 m
West Usambara: 1310–1650 m
Submontane forest and forest edge
Rubeho: 1660 m
Montane forest edge
BPP estimates of Theta (θ) show essentially stable mean and 95 % confidence intervals for all splits within the tree (e.g., θ from median gamma prior across nodes: ~ 0.003 (~0.0007, 0.007)). However, estimates from the limited-distribution species are based on 1–3 individuals, and thus are not informative. Stable θ estimates of the species with larger distributions (and corresponding divergence nodes) agree with EBSP population estimates that these populations have been stable through time without marked reductions or expansions (Additional file 2: Figure S2).
ANOVA and Kruskal-Wallis (KW) analyses of morphology by species
Anova sum Sq
Anova mean Sq
Anova sum Sq
Anova mean Sq
Environmental niche divergence, timing of shifts, and range sizes
In testing potentially peripatric sister species pairs for niche shifts, H. davenporti falls outside of the PCA confidence ellipse probability threshold for H. burgessi into different niche space (Fig. 5). The presence of H. davenporti and H. burgessi in two divergent areas of niche space (mosaic vs. rainforest) implies a major niche shift at the time of speciation. Hyperolius ukwiva localities were outside of the PCA confidence ellipse probability threshold for H. minutissimus (Fig. 5), but only separated by PC2. The appropriate comparison for H. tanneri in a phylogenetic context is between H. tanneri and the ancestor of H. burgessi, H. spinigularis, and H. davenporti. To approximate their ancestral niche, we combined distributional data from these three species into a lumped distribution, which we then compared to H. tanneri in a PCA analysis as above. No separation was seen between H. tanneri and the combined niche space of the “spinigularis” group (H. burgessi, H. spinigularis, and H. davenporti) (Additional file 2: Figure S3).
The two cases of adaptation to grasslands (basal split and H. burgessi/H. davenporti) occurred at differing times (Fig. 5, ~5 MYA and ~ 1 MYA) and therefore a single event cannot have caused these shifts. Both events coincide with periods of increased climatic cycling and a gradual increase in coverage from the forest/grassland mosaic from the previously extensively distributed rainforests of this region. Geographic range estimates (Table 3), underline range asymmetries between sister species .
New lineages within the spiny-throated reed frog clade
Prior to this study and our complementary taxonomic revision , only three species were recognized within the spiny-throated reed frog radiation within the highlands of East Africa: H. spinigularis, H. tanneri, and H. minutissimus. Our field surveys in isolated mountain ranges and associated molecular analyses and morphological assessments of members within this group have documented three new lineages within this radiation: H. burgessi, H. davenporti, and H. ukwiva, each corresponding to a different mountain block or blocks within the EABH. In addition to morphological diagnoses, species delimitation methods based on molecular data were in general agreement in the recognition of each of these species. Deeper population-level sampling should help clarify the levels of divergence within and between species. Nodal support was strong in most methods, though *BEAST showed weak support at basal nodes. Increased sampling of additional molecular markers may help clarify relationships in these deeply divergent lineages.
Peripatric speciation or not
Three species within this radiation are known from single localities in the mountains of the EABH (H. tanneri, H. davenporti, H. ukwiva), raising the question of whether these lineages are the result of peripatric speciation events when isolated from their broadly distributed sister species. Though evidence is limited, many cases of speciation within the spiny-throated reed frogs are consistent with peripatric speciation. In support of peripatric origins, the three species with limited distributions all have smaller realized and potential distributions than their broadly distributed sister lineages (Table 3). Uneven population sizes and population bottlenecks could not be explicitly tested in peripatric species as field observations showed very low population densities and genetic material was limited, yet the rarity of populations and low population sizes of the limited-distribution species were apparent from repeated surveys across localities . In contrast, their sister lineages showed much larger distributions with large and stable population sizes through time. Minimally, the lack of bottleneck patterns from the well sampled, broadly distributed species and limited distribution of sister lineages fits with our a priori predictions for peripatric speciation.
We expected peripatric lineages would differ in climate/habitat, and this was found in 2/3 pairs. Though a habitat shift is not necessary for peripatric speciation, these shifts are less likely in vicariant speciation. The expectation of accentuated divergence in sexually-dimorphic traits in peripatric lineages compared to allopatric was not altogether supported, as all species differed in SUL, gular shape, and the arrangement of gular spines which we assume to be involved in some aspect of sexual selection. There is, however, a complete loss of gular spines in H. tanneri differing substantially from its sister group. This rapid morphological change in H. tanneri, potentially the result of sexual selection, is what would be predicted for peripatric speciation.
Within the allopatric populations of both H. burgessi and H. spinigularis, we find conserved ecological niche/environmental tolerance despite population distances greater than most peripatric species pairs, and no evidence of population bottlenecks. As most species within this radiation are known from single mountain blocks, we were concerned that morphological variation would be high even within species on different mountain blocks (H. spinigularis and H. burgessi have distributions spanning multiple mountain blocks). However, this was not the case in H. burgessi, where morphology was conserved despite substantial divergence times (ca. 0.3-1 mya) between different blocks. Sufficient morphological data are not yet available for the newly discovered Mozambique population of H. spinigularis to determine if they are similarly conserved in relation to Malawian specimens.
Distinguishing between allopatric and peripatric speciation can be difficult, as these processes anchor a continuum of progressively asymmetric isolation. The precise point at which a process ceases to be allopatric and becomes peripatric is therefore blurred. The cases of mixed support in the indicators we examined might therefore be a consequence of the extent to which some of the events were truly “peripatric” or “allopatric” events. Likewise, the indicators themselves (current range sizes, potential ecological niche shifts, and accelerated morphological diversification) may yield ambiguous results when evaluating species pairs that are deeply divergent. With the current data available on the spiny-throated reed frogs, it seems possible that peripatric speciation plays a role in generating biodiversity within the forests and grasslands of the EABH. In systems such as the EABH, which contain a high concentration of old endemics with limited and fragmented distributions in isolated highlands , a community approach may be necessary to clarify the prevalence of these speciation mechanisms in order to more fully understand the formation of its rich biota.
Biome specific speciation mechanisms
Based on the comparisons outlined here from spiny-throated hyperoliids found across these habitats – it seems that peripatric processes may be more pronounced in higher elevation species compared to those in more persistent rainforests (e.g., ). If there were a correlation between the less persistent montane forest/grassland mosaic sites and peripatric speciation, we would expect a common biogeographic pattern repeated across other species. Little data from other amphibians of peripatry has been reported worldwide (e.g., [41, 54, 55]), however, and information is particularly lacking for amphibians within the EABH. Similarly little is known about the phylogenetic relationships and demographic history among EABH high-altitude species and whether patterns exist that might be consistent with peripatry. Evidence in analogous systems, however, suggests that such a distinct boundary in biomes – as in rainforest and grasslands - might be important for generating specific speciation mechanisms. For example, ecotones across Savanna-Forest mosaics of Central Africa been suggested to be important in generating rainforest diversity through peri – and para – patric processes ; however, the generality of such patterns has not be quantified. The role of distinct habitat boundaries promoting speciation processes clearly deserves further investigation, potentially such situations might favor fundamentally different speciation processes compared to other areas (such as vicariance within lower elevation forest habitats). Likewise, there are also clear conservation implications that arise from these considerations. If high elevation species consistently exhibit limited population numbers and historical population genetic bottlenecks, they may be at increased risk if exposed to global anthropogenic pressures such as habitat destruction, fragmentation, and climate change.
Molecular, morphological, range, and niche analyses provide mixed support for cases of peripatric speciation in spiny-throated reed frogs. The highest support is in species pairs exhibiting extreme potential and realized distributional asymmetry, niche shifts, and morphological divergence driven by sexually selected traits: H. burgessi vs. H. davenporti; H. minutissimus vs. H. ukwiva. The older break between H. tanneri and the “spinigularis” clade of H. spinigularis/burgessi/davenporti is less clear, though the limited distribution and unique loss of gular spines (in this clade) in H. tanneri both support the potential for a peripatric split across the 10 km gap separating their current distributions. Though low survey numbers and extremely limited range size in potentially peripatric populations limited our ability to employ demographic reconstructions to confirm peripatric speciation from a signature of founding bottlenecks within the narrowly-distributed species, these conditions themselves support that these are small and divergent populations on the periphery of sister lineages with larger geographic distributions. The emergence of the highland forest/grassland mosaic ecosystems during the same historical interval of this radiation may have facilitated a peripatric mode of speciation as these novel habitats created dispersal opportunities within the larger ranges of rainforest adapted species. Whether the high elevation species that occur along the edge of forest and grassland areas are more amenable to peripatric speciation than rainforest-restricted species awaits additional assessment across more taxonomic groups. However, biome specific peripatric speciation process appears to be supported in some EABH hyperoliids.
We sampled 40 individuals (Hyperolius spinigularis = 8, H. tanneri = 3, H. burgessi = 21, H. ukwiva = 1, H. minutissimus = 6, H. davenporti = 1) from nine mountain blocks, covering the entire range of the spiny-throated reed frog complex (Fig. 1, Sample sizes, Museum numbers and GPS coordinates in Additional file 1: Table S1). Outgroup members of the genus Hyperolius and the family Hyperoliidae were also collected from the field. We collected full voucher specimens, tissue samples, photographs, and microhabitat information. All specimens for this study were collected in accordance with animal ethics guidelines established in the institutions of authors (including Field Museum of Natural History, Science Museo of Trento, and University of Basel).
Primers use in this study
Crawford 2003 
Crawford 2003 
Cmyc H int
Lawson 2010 
Cmyc L int
Lawson 2010 
Wiens et al. 2005 
Wiens et al. 2005 
Lawson 2010 
Lawson 2010 
This study – designed internal to Hoegg et al. 2004 
This study – designed internal to Hoegg et al. 2004 
Alignments were performed using MUSCLE  as described in . Alignments were unambiguous in ND2, POMC, and Rag-1, and easily aligned in C-myc despite the indel-prone intron region. All protein-coding regions were translated into amino acids to check for errors. Substitution rates for each locus are described in Lawson (2010) : ND2 0.00957/lineage/my ; C-myc 0.0006334/lineage/my ; POMC 0.000721/lineage/my ; and Rag-1 0.00042/lineage/my , and used for approximate estimates of diversification events. One nuclear allele from each individual (either cloned or phased as outlined above) was randomly discarded to create a 1:1 mtDNA:nuDNA dataset. Summary information for each locus shown in Table 1.
Phylogenetic reconstruction and divergence time estimation
Sequences from  and 10 additional East African Hyperolius and Afrixalus were included in initial analyses to assess monophyly of the spiny-throated reed frogs (Additional files) including an additional proposed close relative, H. parkeri , using phylogenetic methods described below. A number of datasets were used for initial assessment of phylogenetic relationships: species tree (single representative per species and all outgroups), population/species tree (all spiny-throated specimens (ingroup) and all outgroups), and population-level analysis of only spiny-throated specimens. Monophyly of the ingroup individuals was strongly supported in all cases. Final Bayesian inference (BI) and maximum likelihood (ML) phylogenetic analyses were then completed with the full ingroup dataset (all spiny-throated specimens). RAxML (ML) and BPP (BI) analysis included a single arbitrary close outgroup (H. mitchelli), while the BEAST (BI) and *BEAST (BI) rate-calibrated trees were completed without an outgroup for improved precision of branch length estimates. RAxML (version 7.0.4)  analyses used the rapid hill climbing algorithm and the GTRGAMMA substitution model  partitioned by gene and codon. BEAST (versions 1.8 and 2.1.3)  and *BEAST (starBEAST) , were partitioned by locus and codon position (SRD06 model) . BPP3 [69, 70] parameters listed below in species delimitation. For both ML and BI analyses, model parameters were independently optimized for each partition. All analyses were run for 10 million generations, and run twice to ensure stability (results not combined). The first 10 % of total generations were discarded as burnin for both convergence and tree estimates. Convergence was investigated using Tracer (version 1.6)  through a visual inspection of adequate mixing and ESS estimates >200. The maximum clade credibility tree was calculated for BEAST and *BEAST trees using TreeAnnotator in BEAST. ML node support was evaluated by non-parametric bootstrapping  with 1000 replicates performed in RAxML. The Majority Rule Consensus tree from BPP3 is reported from Evolver in PAML 4 .
BEAST analysis was run with a coalescent, constant size tree-prior and a strict molecular clock (as recommended for recent population-level analyses). *BEAST and BPP3 analysis used individuals from the same mountain block as discrete units (after confirmation of monophyly from individual-based tree constructions) with the exception of a single unit containing individuals from the East Usambara Mountains and Nguru Mountains which were not mutually monophyletic. A Birth-Death tree prior was used in *BEAST due to a priori hypotheses that this lineage has likely undergone local population extinctions through its history.
The program BPP3 [69, 70] was used to jointly infer the species tree and species limits. Default priors and settings were used: gamma prior = G(2, 1000) with mean 2/2000 = 0.001 for population size parameters (s, θ); gamma prior G(2, 1000) for the species tree root age (τ0); all other divergence time parameters are assigned the Dirichlet prior  equation two). Analyses were also run with a very large (1,10) and a very small mean (1, 10000) for θ and τ0 to ensure stability of results. A heredity file (G(4,5)) and locus rate file were incorporated to account for the combined mitochondrial and nuclear datasets and associated mutation rates. Each analysis was run twice to confirm consistency between runs. Species delimitation was assessed with rjMCMC algorithm 0, e = 2.5.
To establish species boundaries within the spiny-throated reed frogs, a Bayesian General Mixed Yule-Coalescent (bGMYC) model was implemented in the package bGMYC  in R v. 3.0.3  using 100 random trees from the concatenated BEAST. Simulations were set at 50,000 generations with 40,000 burn-in, sampling every 100th generation. The upper threshold for the number of species was set at 40, to match the number of tips on the tree.
Bayes Factor species Delimitation (BFD; ) was used to compare three alternative species scenarios using stepping stone and pathsampling analysis in *BEAST summarized in Table 2. Based on the depth of nodes and morphological divergence, H. tanneri, H. minutissimus, and H. ukwiva are all considered separate species in all analyses. All *BEAST parameters are as above. These approaches were used to define units for analyzing molecular, morphological and ecological parameters among lineages.
Effective population size
BPP3 was used to estimate coalescent-scaled population sizes (θ = 4Neμ) and time of divergence (τ = μ t) throughout this radiation on the fixed species tree recovered from all phylogenetic analyses. The parameters used are the same as in the joint species delimitation and species tree analysis.
To test for signs of demographic expansion or contraction we implemented multilocus extended Bayesian skyline plots (EBSP) in BEAST v. 1.8 on populations with n > 5 (Mulanje Massif, Malawi; East Usambara Mountains, Tanzania; Nguru Mountains, Tanzania; Udzungwa Mountains, Tanzania) and all individuals within the “spinigularis” clade to calculate population size through time and the probability of bottlenecks or expansion. All EBSP analyses used a strict clock. Operators were modified according to author recommendations and analyses were run for 10 million generations to obtain adequate ESS values.
To determine morphological distinctiveness among members of the species cluster, we took commonly-used limb and cranial measurements for 17 traits , including leg bone (e.g. femur, tibia) and foot lengths and head width (see below). Many of these traits are often assumed to be adaptive, as they are correlated with food source (e.g., head width), perching (e.g., foot length), and other life history traits. If species are morphologically divergent, it is possible that competition or adaptation to different environmental conditions reinforced separation within this radiation. Morphological measurements made in this analysis include: Snout-Urostyle Length (SUL), Head Width (HW), Head Length Diagonal from corner of mouth (HLD), Head Length Diagonal from jawbone end (HLDJ), Nostril-Snout (NS), Inter-narial (IN), Eye to Nostril (EN), Eye Distance (EE), Inter-orbital (IO), Tibiafibula Length (TL), Thigh Length (THL), Tibiale Fibulare Length (TFL), Foot Length (FL), Forelimb Length (FLL), Hand Length (HL), Width of Gular Flap (WGF) Height of Gular Flap (HGF). We also scored specimens for traits associated with sexual selection including presence and distribution of gular spines (which are only present during the breeding season for species that possess them) and color patterns. Morphological measurements were taken of 114 mature specimens using Mitutoyo Absolute Digimatic Calipers (CD-6"C) (Sample sizes and Measurements in supplementary materials). In order to assess the morphological distinctness of these species, we conducted Principal Component analyses on log-transformed data of males and females separately with all variables centered and scaled in R  (95 % confidence ellipse probability threshold from ggbiplot package). To further evaluate whether species were significantly different, we did a Permutational Multivariate Analysis of Variance Using Distance Matrices using the adonis function from the vegan  R package.
Kruskal–Wallis and ANOVA were used to assess variable divergence (standardized by body size) between species for use as diagnostic characters in R.
Niche similarity of all species was assessed by Principle Component Analysis (PCA) using bioclim variables associated with GPS coordinates (Additinal file 1: Table S1) in the MASS and ggbiplot packages in R (95 % confidence ellipse probability threshold). This method allows at least a preliminary view of habitat similarity for severely range-restricted taxa such as the members of this complex. We evaluated both a full bioclim dataset (all 19 standard bioclim variables) and a reduced dataset with Pearson’s correlation coefficients below 0.7: Mean Diurnal Range, Temperature Seasonality, Temperature Annual Range, Mean Temperature of Coldest Quarter, Precipitation of Wettest Month, Precipitation Seasonality, Precipitation of Driest Quarter, Precipitation of Warmest Quarter, Precipitation of Coldest Quarter (30 arc-second resolution, WorldClim database ). In reconstructing ecological shifts across the phylogenetic tree, we rely on the assumption that a species niche remains stable through time after divergence. Though this appears to be widely true [16, 81, 82], there is increased uncertainty associated with assessing extremely locally adapted lineages were spatial autocorrelation might have a more pronounced effect.
Range estimates and field observations of population densities
Species’ range estimates were calculated using standardized approaches conducted on species being evaluated for the IUCN Red List of Threatened Species™ . Estimates of population size were made using molecular techniques when sufficient numbers of individuals were available (above). For those species with very few individuals ever found (H. tanneri, H. ukwiva, H. davenporti), estimates were qualitatively assessed from sampling effort over multiple field seasons (Lawson, Menegon, Loader, personal observations, ) in comparison to all known localities from other spiny-throated reed frogs (H. spinigularis, H. burgessi, H. minutissimus).
We would like to thank J. Miller, H. Müller, D. Blackburn, D. Portik, E. Tonelli, C. Barratt, T. Price, J. Luedtke, & L. Coals for tissues, data and comments on early versions of this manuscript. Also thanks to M. Munissi, E. Mulungu, Nyson, L. Mahler, and B. Zimkus for help in the field. K. Howell, D. Moyer, L. Mazibuko, and J. Poynton provided extensive knowledge of African Amphibians and collecting logistics in Tanzania and Malawi. P. Kaliba was instrumental for fieldwork in Malawi. Special thanks to TAWIRI, COSTEC, Tanzanian Wildlife and Forestry departments for research permits (RCA 2006, 2007-72-Na-2006-19, RCA 2001–272; RCA 2007–153, RCA 2009-306-NA-2009-201, 2011-239-NA-2011-82), and Museums of Malawi for Malawian permits. Thanks to K. Feldheim and the Pritzker lab for help on molecular work and for creating the taq used in this study. Lawson’s PhD work was funded through a GAANN educational grant to the Committee on Evolutionary Biology, University of Chicago. This project was funded through National Science Foundation DDEP grant (NSF 0841641), Hinds Fund of the University of Chicago, and grants through the Council on Africa (Field Museum of Natural History). SPL was funded to conduct surveys and lab work by the following institutes: University of Basel, Swiss National Science Foundation (grant number 31003A-133067), SCNAT, Freiwillige Akademische Gesellschaft Basel, and the Council on Africa (Field Museum of Natural History). Sponsors were not involved in collection, analysis, or interpretation of data, or writing and submitting this manuscript. Finally, we thank the editor (John Wiens) and reviewers for helpful and constructive criticisms that improved the manuscript.
- Mayr E. Animal Species and Evolution. Cambridge, MA: Harvard University Press; 1963.View ArticleGoogle Scholar
- Lomolino MV, Riddle BR, Whittaker RJ, Brown JH. Biogeography. 4th ed. Sunderland, MA: Sinauer Associates; 2010.Google Scholar
- Coyne JA, Orr HA. Speciation. Sinauer Associates Sunderland, MA: Sunderland, MA; 2004.Google Scholar
- Zink RM. The Geography of Speciation: Case Studies from Birds. Evo Edu Outreach. 2012;5:541–6.View ArticleGoogle Scholar
- Hickerson MJ, Stahl EA, Lessios HA. Test for Simultaneous Divergence Using Approximate Bayesian Computation. Evolution. 2006;60:2435–53.PubMedView ArticleGoogle Scholar
- Soltis DE, Morris AB, McLachlan JS, Manos PS, Soltis PS. Comparative phylogeography of unglaciated eastern North America. Mol Ecol. 2006;15:4261–93.PubMedView ArticleGoogle Scholar
- Carstens BC, Richards CL. Integrating Coalescent and Ecological Niche Modeling in Comparative Phylogeography. Evolution. 2007;61:1439–54.PubMedView ArticleGoogle Scholar
- Victoriano PF, Ortiz JC, Benavides E, Adams BJ, Sites Jr JW. Comparative phylogeography of codistributed species of Chilean Liolaemus (Squamata: Tropiduridae) from the central-southern Andean range. Mol Ecol. 2008;17:2397–416.PubMedView ArticleGoogle Scholar
- Hickerson MJ, Meyer CP. Testing comparative phylogeographic models of marine vicariance and dispersal using a hierarchical Bayesian approach. BMC Evol Biol. 2008;8:322.PubMed CentralPubMedView ArticleGoogle Scholar
- Kozak KH, Graham CH, Wiens JJ. Integrating GIS-based environmental data into evolutionary biology. Trends Ecol Evol. 2008;23:141–8.PubMedView ArticleGoogle Scholar
- Knowles LL. Statistical phylogeography. Annu Rev Ecol Evol Syst. 2009;40:593–612.View ArticleGoogle Scholar
- Freckleton RP, Jetz W. Space versus phylogeny: disentangling phylogenetic and spatial signals in comparative data. Proc R Soc B. 2009;276:21–30.PubMed CentralPubMedView ArticleGoogle Scholar
- Pigot AL, Phillimore AB, Owens IPF, Orme CDL. The Shape and Temporal Dynamics of Phylogenetic Trees Arising from Geographic Speciation. Syst Biol. 2010;59:660–73.PubMedView ArticleGoogle Scholar
- Losos J, Glor R. Phylogenetic comparative methods and the geography of speciation. Trends Ecol Evol. 2003;18:220–7.View ArticleGoogle Scholar
- Blair C, Heckman KL, Russell AL, Yoder AD. Multilocus coalescent analyses reveal the demographic history and speciation patterns of mouse lemur sister species. BMC Evol Biol. 2014;14:57.PubMed CentralPubMedView ArticleGoogle Scholar
- Peterson AT, Soberón J, Sánchez-Cordero V. Conservatism of Ecological Niches in Evolutionary Time. Science. 1999;285:1265–7.PubMedView ArticleGoogle Scholar
- Schluter D. The Ecology of Adaptive Radiation. New York: Oxford University Press; 2000.Google Scholar
- Dieckmann U, Doebeli M, Metz J, Tautz D. Adaptive Speciation. Cambridge, UK: Cambridge Univ. Press; 2004.View ArticleGoogle Scholar
- Kaneshiro KY. The dynamics of sexual selection and founder effects in speciation formation. In: Kaneshiro KY, Giddings L, Anderson W, editors. Genetics, speciation, and the founder principle. Oxford, UK: Oxford University Press; 1989. p. 279–96.Google Scholar
- Panhuis TM, Butlin R, Zuk M, Tregenza T. Sexual selection and speciation. Trends Ecol Evol. 2001;16:364–71.PubMedView ArticleGoogle Scholar
- Ödeen A, Florin A-B. Sexual selection and peripatric speciation: the Kaneshiro model revisited. J Evol Biol. 2002;15:301–6.View ArticleGoogle Scholar
- Carson HL, Templeton AR. Genetic Revolutions in Relation to Speciation Phenomena: The Founding of New Populations. Annu Rev Ecol Evol Syst. 1984;15:97–131.View ArticleGoogle Scholar
- Roy MS. Recent diversification in African greenbuls (Pycnonotidae: Andropadus) supports a montane speciation model. Proc R Soc Lond B. 1997;264:1337–44.View ArticleGoogle Scholar
- Linder HP. The radiation of the Cape flora, southern Africa. Biol Rev. 2003;78:597–638.PubMedView ArticleGoogle Scholar
- Bowie RCK, Fjeldså J, Hackett SJ, Bates JM, Crowe TM. Coalescent models reveal the relative roles of ancestral polymorphism, vicariance, and dispersal in shaping phylogeographical structure of an African montane forest robin. Mol Phylo Evol. 2006;38:171–88.View ArticleGoogle Scholar
- Lawson LP. The discordance of diversification: evolution in the tropical-montane frogs of the Eastern Arc Mountains of Tanzania. Mol Ecol. 2010;19:4046–60.PubMedView ArticleGoogle Scholar
- Linder HP, de Klerk HM, Born J, Burgess ND, Fjeldså J, Rahbek C. The partitioning of Africa: statistically defined biogeographical regions in sub-Saharan Africa. J Biogeogr. 2012;39:1189–205.View ArticleGoogle Scholar
- Demos TC, Kerbis Peterhans JC, Agwanda B, Hickerson MJ. Uncovering cryptic diversity and refugial persistence among small mammal lineages across the Eastern Afromontane biodiversity hotspot. Mol Phylo Evol. 2014;71:41–54.View ArticleGoogle Scholar
- Ceccarelli FS, Menegon M, Tolley KA, Tilbury CR, Gower DJ, Laserna MH, Kasahun R, Rodriguez-Prieto A, Hagmann R, Loader SP. Evolutionary relationships, species delimitation and biogeography of Eastern Afromontane horned chameleons (Chamaeleonidae: Trioceros). Mol Phylo Evol. 2014;80:125–36.View ArticleGoogle Scholar
- Poynton JC, Loader SP, Sherratt E, Clarke BT. Amphibian Diversity in East African Biodiversity Hotspots: Altitudinal and latitudinal Patterns. Biodivers Conserv. 2007;16:1103–18.View ArticleGoogle Scholar
- Fjeldså J, Bowie RCK, Rahbek C. The Role of Mountain Ranges in the Diversification of Birds. Annu Rev Ecol Evol Syst. 2012;43:249–65.View ArticleGoogle Scholar
- Hua X, Wiens JJ. Latitudinal Variation in Speciation Mechanisms in Frogs. Evolution. 2010;64:429–43.PubMedView ArticleGoogle Scholar
- Kozak KH, Wiens JJ. Climatic zonation drives latitudinal variation in speciation mechanisms. Proc R Soc B. 2007;274:2995–3003.PubMed CentralPubMedView ArticleGoogle Scholar
- Lovett JC. Eastern Arc moist forest flora. In: Biogeography and ecology of the rain forests of eastern Africa. Cambridge, UK: Cambridge University Press; 1993. p. 33–55.View ArticleGoogle Scholar
- De Menocal PB. Plio-Pleistocene African Climate. Science. 1995;270:53–9.View ArticleGoogle Scholar
- Meadows ME, Linder HP. A Palaeoecological Perspective on the Origin of Afromontane Grasslands. J Biogeogr. 1993;20:345–55.View ArticleGoogle Scholar
- Hamilton AC. The quaternary history of African forests: its relevance to conservation. Afr J Ecol. 1981;19:1–6.View ArticleGoogle Scholar
- Lovett JC, Marchant R, Taplin J, Küper W. The oldest rainforests in Africa: stability or resilience for survival and diversity. In: Purvis A, Gittleman JL, Brooks TM, editors. Phylogeny and conservation. Cambridge, UK: Cambridge University Press; 2005. p. 115–25.Google Scholar
- Marchant R, Mumbi C, Behera S, Yamagata T. The Indian Ocean dipole – the unsung driver of climatic variability in East Africa. Afr J Ecol. 2007;45:4–16.View ArticleGoogle Scholar
- Mumbi CT, Marchant R, Hooghiemstra H, Wooller MJ. Late Quaternary vegetation reconstruction from the Eastern Arc Mountains, Tanzania. Quat Res. 2008;69:326–41.View ArticleGoogle Scholar
- Loader SP, Sara Ceccarelli F, Menegon M, Howell KM, Kassahun R, Mengistu AA, Saber SA, Gebresenbet F, de Sá R, Davenport TRB, Larson JG, Müller H, Wilkinson M, Gower DJ. Persistence and stability of Eastern Afromontane forests: evidence from brevicipitid frogs. J Biogeogr. 2014;41:1781–92.View ArticleGoogle Scholar
- Crisp M, Cook L, Steane D. Radiation of the Australian flora: what can comparisons of molecular phylogenies across multiple taxa tell us about the evolution of diversity in present–day communities? Phil Trans R Soc Lond B. 2004;359:1551–71.View ArticleGoogle Scholar
- Cantalapiedra JL, Hernández Fernández M, Morales J. Biomic Specialization and Speciation Rates in Ruminants (Cetartiodactyla, Mammalia): A Test of the Resource-Use Hypothesis at the Global Scale. PLoS One. 2011;6, e28749.PubMed CentralPubMedView ArticleGoogle Scholar
- Johnson SD, Linder HP. Speciation and extinction in the Greater Cape Floristic Region. In: Allsopp N, Colville JF, Verboom GA, editors. Fynbos: Ecology, Evolution, and Conservation of a Megadiverse Region. Oxford, UK: Oxford University Press; 2014. p. 119.Google Scholar
- Loader SP, Lawson LP, Portik DM, Menegon M. Three new species of spiny throated reed frogs (Anura: Hyperoliidae) from Evergreen forests of Tanzania. in press. BMC Res Notes. 8(1):167Google Scholar
- Schiøtz A. Two New Hyperolius (Anura) from Tanzania. Steenstrupia. 1982;8:269–76.Google Scholar
- Schiøtz A. Treefrogs of Africa. Chimairath ed. Germany: Frankfurt am Main; 1999.Google Scholar
- Schiøtz A, Westergaard MM. Notes on some Hyperolius (Anura: Hyperoliidae) from Tanzania, with supplementary information on two recently described species. Steenstrupia. 2000;25:1–9.Google Scholar
- Portik DM, Mulungu EA, Sequeira D, McEntee JP. Herpetological surveys of the Serra Jeci and Namuli massifs, Mozambique, and an annotated checklist of the southern Afromontane archipelago. Herpetol Rev. 2013;44:394–406.Google Scholar
- Cerling TE, Bowman JR, O’Neil JR. An isotopic study of a fluvial-lacustrine sequence: The Plio-Pleistocene koobi fora sequence, East Africa. Palaeogeogr Palaeoclimatol Palaeoecol. 1988;63:335–56.View ArticleGoogle Scholar
- Trauth MH, Larrasoaña JC, Mudelsee M. Trends, rhythms and events in Plio-Pleistocene African climate. Quat Sci Rev. 2009;28:399–411.View ArticleGoogle Scholar
- IUCN SSC Amphibian Specialist Group 2013. Hyperolius spinigularis/Hyperolius minutissimus. The IUCN Red List of Threatened Species. Version 2014.3. <www.iucnredlist.org> Downloaded on 01 June 2015.
- Rovero F, Menegon M, Fjeldså J, Collett L, Doggart N, Leonard C, Norton G, Owen N, Perkin A, Spitale D, Ahrends A, Burgess ND. Targeted vertebrate surveys enhance the faunal importance and improve explanatory models within the Eastern Arc Mountains of Kenya and Tanzania. Diversity Distrib. 2014;20:1438–49.View ArticleGoogle Scholar
- Rovito SM. Lineage divergence and speciation in the Web-toed Salamanders (Plethodontidae: Hydromantes) of the Sierra Nevada, California. Mol Ecol. 2010;19:4554–71.PubMedView ArticleGoogle Scholar
- Rovito SM, Wake DB, Papenfuss TJ, Parra-Olea G, Muñoz-Alonso A, Vásquez-Almazán CR. Species formation and geographical range evolution in a genus of Central American cloud forest salamanders (Dendrotriton). J Biogeogr. 2012;39:1251–65.View ArticleGoogle Scholar
- Smith TB, Wayne RK, Girman DJ, Bruford MW. A Role for Ecotones in Generating Rainforest Biodiversity. Science. 1997;276:1855–7.View ArticleGoogle Scholar
- Seutin G, White BN, Boag PT. Preservation of avian blood and tissue samples for DNA analyses. Can J Zool. 1991;69:82–90.View ArticleGoogle Scholar
- Hoegg S, Vences M, Brinkmann H, Meyer A. Phylogeny and Comparative Substitution Rates of Frogs Inferred from Sequences of Three Nuclear Genes. Mol Biol Evol. 2004;21:1188–200.PubMedView ArticleGoogle Scholar
- Wiens JJ, Fetzner JW, Parkinson CL, Reeder TW. Hylid frog phylogeny and sampling strategies for speciose clades. Sys Biol. 2005;54:719–748.View ArticleGoogle Scholar
- Edgar RC. MUSCLE: Multiple Sequence Alignment with High Accuracy and High Throughput. Nucl Acids Res. 2004;32:1792–7.PubMed CentralPubMedView ArticleGoogle Scholar
- Crawford AJ. Relative Rates of Nucleotide Substitution in Frogs. J Mol Evol. 2003;57:636–41.PubMedView ArticleGoogle Scholar
- Evans BJ, McGuire JA, Brown RM, Andayani N, Supriatna J. A coalescent framework for comparing alternative models of population structure with genetic data: evolution of Celebes toads. Biol Lett. 2008;4:430–3.PubMed CentralPubMedView ArticleGoogle Scholar
- Pickersgill M. Frog Search, Results of Expeditions to Southern and Eastern Africa. Chimairath ed. Germany: Frankfurt am Main; 2007.Google Scholar
- Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.PubMedView ArticleGoogle Scholar
- Stamatakis A, Blagojevic F, Nikolopoulos DS, Antonopoulos CD. Exploring New Search Algorithms and Hardware for Phylogenetics: RAxML Meets the IBM Cell. J VLSI Sign Process Syst Sign Im. 2007;48:271–86.View ArticleGoogle Scholar
- Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:1–8.View ArticleGoogle Scholar
- Heled J, Drummond AJ. Bayesian inference of species trees from multilocus data. Mol Biol Evol. 2010;27:570–80.PubMed CentralPubMedView ArticleGoogle Scholar
- Shapiro B, Rambaut A, Drummond AJ. Choosing Appropriate Substitution Models for the Phylogenetic Analysis of Protein-Coding Sequences. Mol Biol Evol. 2006;23:7–9.PubMedView ArticleGoogle Scholar
- Rannala B, Yang Z. Bayes estimation of species divergence times and ancestral population sizes using DNA sequences from multiple loci. Genetics. 2003;164:1645–56.PubMed CentralPubMedGoogle Scholar
- Yang Z, Rannala B. Bayesian species delimitation using multilocus sequence data. PNAS. 2010;107:9264–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Rambaut A, Drummond AJ. Tracer ver. 1.5. Program available at 2007; http://tree.bio.ed.ac.uk/software/tracer/.
- Felsenstein J. Phylogenies and the comparative method. Am Nat. 1985;125:1–15.View ArticleGoogle Scholar
- Yang Z. PAML 4: Phylogenetic Analysis by Maximum Likelihood. Mol Biol Evol. 2007;24:1586–91.PubMedView ArticleGoogle Scholar
- Reid NM, Carstens BC. Phylogenetic estimation error can decrease the accuracy of species delimitation: a Bayesian implementation of the general mixed Yule-coalescent model. BMC Evol Biol. 2012;12:196.PubMed CentralPubMedView ArticleGoogle Scholar
- R Development Core Team. R: A language and environment for statistical computing. Austria: R Foundation for Statistical Computing Vienna; 2011. 3.Google Scholar
- Grummer JA, Bryson RW, Reeder TW. Species Delimitation Using Bayes Factors: Simulations and Application to the Sceloporus scalaris Species Group (Squamata: Phrynosomatidae). Syst Biol. 2014;63:119–33.PubMedView ArticleGoogle Scholar
- Moen DS, Wiens JJ. Phylogenetic Evidence for Competitively Driven Divergence: Body-Size Evolution in Caribbean Treefrogs (hylidae: Osteopilus). Evolution. 2009;63:195–214.PubMedView ArticleGoogle Scholar
- Wickham H. ggplot2: Elegant Graphics for Data Analysis. New York: Springer; 2009.View ArticleGoogle Scholar
- Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H. Vegan: Community Ecology Package. R Package Version 2.0-10. 2013.Google Scholar
- Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005;25:1965–78.View ArticleGoogle Scholar
- Wiens JJ, Ackerly DD, Allen AP, Anacker BL, Buckley LB, Cornell HV, Damschen EI, Jonathan Davies T, Grytnes J-A, Harrison SP, Hawkins BA, Holt RD, McCain CM, Stephens PR. Niche conservatism as an emerging principle in ecology and conservation biology. Ecol Lett. 2010;13:1310–24.PubMedView ArticleGoogle Scholar
- Martínez-Meyer E, Townsend Peterson A, Hargrove WW. Ecological niches as stable distributional constraints on mammal species, with implications for Pleistocene extinctions and climate change projections for biodiversity. Global Ecol Biogeogr. 2004;13:305–14.View ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.