Multilocus Phylogeography of the Treefrog Scinax eurydice (Anura, Hylidae) Reveals a Plio-Pleistocene Diversification in the Atlantic Forest

We aim to evaluate the genetic structure of an Atlantic Forest amphibian species, Scinax eurydice, testing the congruence among patterns identified and proposed by the literature for Pleistocene refugia, microrefugia, and geographic barriers to gene flow such as major rivers. Furthermore, we aim to evaluate predictions of such barriers and refugia on the genetic structure of the species, such as presence/absence of dispersal, timing since separation, and population expansions/contractions. We sequenced mitochondrial and nuclear genetic markers on 94 tissue samples from 41 localities. We inferred a gene tree and estimated genetic distances using mtDNA sequences. We then ran population clustering and assignment methods, AMOVA, and estimated migration rates among populations identified through mtDNA and nDNA analyses. We used a dated species tree, skyline plots, and summary statistics to evaluate concordance between population’s distributions and geographic barriers and Pleistocene refugia. Scinax eurydice showed high mtDNA divergences and four clearly distinct mtDNA lineages. Species tree and population assignment tests supported the existence of two major clades corresponding to northeastern and southeastern Atlantic Forest in Brazil, each one composed of two other clades. Lineage splitting events occurred from late Pliocene to Pleistocene. We identified demographic expansions in two clades, and inexistent to low levels of migrations among different populations. Genetic patterns and demographic data support the existence of two northern Refuge and corroborate microrefugia south of the Doce/Jequitinhonha Rivers biogeographic divide. The results agree with a scenario of recent demographic expansion of lowland taxa. Scinax eurydice comprises a species complex, harboring undescribed taxa consistent with Pleistocene refugia. Two rivers lie at the boundaries among populations and endorse their role as secondary barriers to gene flow.


Introduction
The Atlantic Forest (hereafter AF) harbors one of the richest, most threatened and unique biotas in the world [1].A high degree of human interference (less than 12% of the original distribution remaining), high levels of endemism [567 vertebrates (2.1%) and over 8000 plants (2.7%)], and continuous rates of species descriptions characterize it as a biodiversity hotspot [1,2].Despite these impressive numbers, little was known about the region's biogeographic history, and only recently a clearer picture about the origin, distribution, and relationships of its biota has arisen [3][4][5][6][7][8][9][10][11].
Biodiversity is not evenly distributed in the AF, with centers of endemism supported by faunal inventories, phylogenetic analyses and phylogeographic studies [4,11].Phylogeographic breaks have been reported for a number of organisms, including bees [6], amphibians [12][13][14], reptiles [7,15], birds [8,16,17], bats [9], and plants [10].The causes of such breaks have been attributed to Pleistocene climatic cycles, rivers acting as barriers, and/or and geomorphological dynamics.Hence, both vicariance promoted by geographic barriers and refuge dynamics prompted by climatic fluctuations have been implicated in the AF's biota diversification [12,15,16].Notwithstanding, some phylogeographic studies involving birds pointed absence of genetic structure through these recovered breaks, as well as demographic stability during the Last Glacial Maximum [18,19].
Geographic barriers associated to neotectonic activity have been evoked to explain evolution of biota from the AF [20].Biogeographic history and evolutionary processes in the AF were deeply influenced by river systems and mountain ranges originated by tectonics dated before Pleistocene, which determined genetic patterns recognized in some taxa [7,14,21].For amphibians, speciation events have been explained by the renewal of the topography in southeastern Brazil during the Tertiary [22], primarily by the elevation that culminated in the Serra do Mar and Serra da Mantiqueira and depressions that caused the formation of rivers like the Paraíba do Sul and the Pomba-Muriaé (see [23]).However, recent phylogeographic studies focusing on the AF biota have shown that divergence events have been ongoing from the late Tertiary through the Pleistocene, evidencing the effects of more recent phenomena in shaping the current distribution of genetic diversity (eg.[12,15,20]).
Pleistocene climatic fluctuations played a central role in the stability, biodiversity distribution, and phylogeographic endemism in the Brazilian Atlantic Forest [3][4][5].Forest models predict historical stability in central and northern Atlantic Forest, corroborating paleopalynological records, areas of endemism, and mtDNA diversity for several taxa [3].Phylogeographic data for amphibians supported historically stable forest areas in the northeastern, central and southeastern [4], while phylogeographic endemism for 25 vertebrates in the Brazilian AF endorse two bioclimatic domains identified by distribution models for present and past climates [5].Phylogeographic work has now supported the subdivision between southern, central, and northern Atlantic Forest biotas [4], indicating that cryptic lineages are likely to be identified in taxa spanning these regions, especially the climatically distinct northern and southern portions [5].
The treefrog Scinax eurydice is widely distributed in the Atlantic Forest and is a species that occupies the forest edge, being able to use temporary and/or permanent pounds on the outer margin and adjacent open areas in the forest for breeding [24][25][26][27][28]. Populations exhibit phenotypic variation throughout the species' range with four distinct morphotypes recognized based on external morphology (body size and shape), color patterns, and bioacoustics (M.F. Napoli et al. "unpublished work").The wide geographic distribution allied with acoustic and morphological variation of the advertisement call provides an ideal scenario to test the congruence of genetic lineages with biogeographic domains and centers of endemism proposed for the AF.
Herein, we evaluate population structure and historical demography of S. eurydice in the AF based on one mitochondrial and three nuclear genes, using 94 tissue samples covering most of the species distribution.Specifically, we addressed the following questions: (i) how many genetically differentiated populations exist along the species distribution?(ii) how long ago did they diverge?and (iii) is there concordance between geographic distribution of genetically differentiated populations and divergence times with possible geographic barriers and/or Pleistocene climatic cycles?

Population sampling
We obtained 94 tissue samples from 41 localities of S. eurydice from several Brazilian zoological collections, which cover most of the range of the species (Fig 1 ; S1 Table, see in Supporting Information).We used S. x-signatus, S. auratus, S. nebulosus and S. berthae as outgroups.Tissue samples for outgroups were obtained from Brazilian and Paraguayan zoological collections (see S1 Table ).

Gene trees and mtDNA genetic distances
We estimated a gene tree for ND2 using Maximum Likelihood (ML) and Bayesian Inference (BI), as implemented inRAxMLGUI v.1.3[33] and MrBayes 3 [34], respectively.Here and in all subsequent analyses substitution models were selected using JmodelTest2 [35] according to the AIC (Akaike Information Criterion).The ML analysis was performed with 100 runs and 1000 bootstrap replications (ML + 'thorough Bootstrap', 'n threads.' 4).For the BI analysis we performed 2 independent MCMC runs, each with four chains for 20 million steps sampling every 1000 generations and discarding the first two million as burn-in.Parameter convergence was checked in Tracer v1.6 (http://beast.bio.ed.ac.uk/Tracer).We calculated the genetic distance between recognized ND2 clades (net average distances-Da), using the Tamura-Nei

Population assignment and basic population genetics statistics
We estimated the most likely number of groups (k) using the Bayesian clustering method implemented in Structure 2.3.3 [37].The mitochondrial fragment was not used in this analysis.We assumed linkage model and independent allele frequencies.The analyses were repeated 10 times for each value of k from 1 to 10.Each k run consisted of six million MCMC generations with a burn-in of one million.To obtain the most likely number of k we used the Evannos's K [38] method implemented in Structure Harvester [39] and the Bayes rule Pr(K) following [40].Finally, ten runs with the most likely k were combined using Clumpp 1.1.2[41] with the Greedy algorithm and 1000 permutations.
To determine the level of population structure among the defined groups we implemented an analysis of molecular variance (AMOVA) for mtDNA and nDNA separately, using the software ARLEQUIN 3.11 [42].Localities with at least two individuals were assigned to populations.The significance was obtained using 1000 permutations.We calculated the following summary statistics for each defined group in DnaSP: haplotype diversity (Hd), nucleotide diversity (π), and their standard deviation (SD); Tajima's D and Fu's F s neutrality tests, and the population size change test Ramos-Onsins & Rozas R 2 for each locus.The significances were calculated using 1000 coalescent simulations.

Divergence time estimates
We estimated divergence times using the species tree method implemented in Ã BEAST 1.8.0 [43].Furthermore, the Isolation with Migration method implemented in IMa2 [44] produces divergence time estimates, that were compared to Ã BEAST estimates.
We estimated the species tree by using Ã BEAST at the CIPRES Science Gateway [45].This method estimates divergence time, population size change through time and species tree topology.Species assignment followed results of population assignment test (see section above and results for details).We used the strict clock model for all gene partition.We fixed ND2 substitution rate at 0.00957 mutation/site/million years [46].For the nuclear loci substitution priors we used a default lognormal distribution (mean = 0 and sdtv = 1) and allowed the program to co-estimate rates.We implemented a piecewise linear and constant root and Yule Process as priors for the species tree model.We performed three independent MCMC runs of 200 million generations, sampling every 10,000 generations with a 10% burn-in.We checked for chain mixing, effective sample sizes (ESS), and convergence between runs using Tracer v1.6, making sure ESS values were higher than 200 for all parameters.The resulting trees were combined and summarized using TreeAnotator.We visualized the entire population of sampled trees in DensiTree 2.0.1 [47].

Migration
IMa2 implements an Isolation-with-Migration model (IM) [44] and estimates effective population size, migration rate and divergence time.With few loci at hand it becomes hard to estimate all parameters of a four-population model [48].Therefore, we ran a series of analyses of IM two-population model for all six pairwise combinations of the populations found [48].We compare the two-population model results with a four-population model analysis with reduced number of parameters (migration only among current populations, ancestral populations had zero migration; exponential prior distribution with mean 0.6 for the migration parameter).
For each two-population analysis we ran IMa2 in M-mode with 20 chains and a geometric heating model (first parameter: 0.96; second parameter: 0.9), sampling 200,000 genealogies with a burn-in of 500,000 steps.We assumed a HKY substitution model for ND2, RAG1 and β-fibint7; and an Infinite Site model (IS) for 28S.Upper prior limits for population parameters were defined following the IMa2 manual (q = 16, t = 6.2, m = 0.6).For the four population model we ran IMa2 in M-mode with 40 chains and geometric heating model (first parameter: 0.975; second parameter: 0.75) and tree topology according to Species tree results.We sampled 100,000 genealogies after a burn-in of 1,000,000 steps.A HKY model was assumed for all genes.In all IMa2 analysis we accessed appropriate chain mixing by ensuring convergent results among 4 replicates with different seed numbers, and by checking the burntrend plot.Migration among populations was verified according to the likelihood ratio test (LRT), which accesses if the posterior distribution of the migration parameter significantly departs from zero.Estimated parameters were transformed to demographic quantities using the geometric mean of the mutation rates estimated by the Species tree analysis (see above).Posterior distribution graphs were generated using R programing environment [49] and ggplot2 R-package [50].

Historical demography and population sizes
To evaluate the population size changes through time we used the Extended Bayesian Skyline Plot method (EBSP) in Ã BEAST v1.8.0 applied separately for each population identified by gene tree analysis.Two independent runs of EBSP were obtained with the following parameters: an initial tree generated by UPGMA, linear model, 200 million steps, parameters sampled every 10,000 steps and 10% burn-in (if necessary, the values of burn-in were changed to suit the ESS values > 200).We used an uncorrelated lognormal distribution for the mitochondrial fragment and a strict clock for nuclear fragments.To calibrate the EBSP we used the same strategy implemented in the Ã BEAST analysis (see section above).Parameter convergence between runs and performance analysis (considered ESS > 200) were conducted in Tracer v1.6 (http://beast.bio.ed.ac.uk/Tracer).We also obtained population size estimates from the IMa2 analysis (see migration section above).

Gene trees and mtDNA genetic distance
We obtained 69 ND2 sequences (969 bp, 35 localities) comprising 49 haplotypes.ML and BI mitochondrial gene trees (GTR+G model) confirmed the monophyly of S. eurydice (Fig 1).Three clades were recovered, NE1, NE2 and SE, corresponding to northeastern and southeastern Brazil separated by the Jequitinhonha River (NE1+ NE2 and SE) and the São Francisco River (NE1 and NE2) (Fig 1).The NE1 clade includes samples from north of the São Francisco River and the NE2 clade comprises samples between the rivers São Francisco and Jequitinhonha.The SE present subclades: SE1 encompassing frogs from Rio de Janeiro State and SE2 including samples above the Jequitinhonha River.Genetic distances between mitochondrial clades ranged from 8% to 16% (Table 1).For nuclear genes we sequenced 91 samples for RAG1 (429 bp), 94 samples for 28S (768bp), and 58 samples for β-fibint7.Phased nuclear sequences resulted in 158 alleles for RAG1, 98 for β-fibint7, and 188 for 28S with P > 0.7.Twenty-four alleles for RAG1 and eighteen for β-fibint7 had low probability with PHASE and were removed from further analysis.The PHI test excluded recombination for β -fibint7 and 28S (P = 0.2539 and P = 1.0, respectively), but not for RAG1 (p = 0.0432).We applied the Four Gamete test to RAG1 to identify the largest non-recombinant block.We identified 46 pairs of sites with four gametic types and a minimum of five recombinant events.For RAG1 three non-recombinant blocks with 64 bp, 113 bp and 136 bp were identified, of which the longest was used in further analyses.The best-fitted substitution models of nuclear genes were as follow: HKY for RAG1 and GTR for β-fibint7 and 28S.
Haplotype networks for nuclear loci showed variable levels of geographic structure (Fig 2).RAG1 showed lower geographic structure among populations than mitochondrial gene trees (ND2) and β-fibint7 showed a northeast-southeast distribution of haplotypes.In contrast, 28S haplotypes showed geographical structure concordant with major mitochondrial clades (NE, SE).Given the mitochondrial phylogeny, the k populations recovered by Evanno's method, and the number of populations selected through Bayes Rule, we chose to assign individuals to four populations.Because admixture in nuclear genes precludes us from unequivocally assigning individuals to four different populations, we used the mitochondrial phylogeny to assign individuals to four different populations.

Migration
Mean values of migration rates for two-population and four-population models are largely in congruent.However, peak estimates for migration from SE2 to SE1 and from NE2 to SE1 are much lower in the four-population model (S4 Table, S3 and S4 Figs).One potential explanation is that for the two-population model, when there is migration with a third not included population, the IMa2 model is violated and migration rates can be biased [51].Nevertheless, the exponential prior used for the migration parameter in the four-population model can also cause this difference in peak estimates.IMa two-population model analyses estimated significant bidirectional migration between populations NE2 and SE2; and populations SE1 and SE2, with higher population migration rates in the southward direction from NE2 to SE2, and from SE2 to SE1.Unidirectional migration was also significant from population NE2 to population SE1 (S4 Table ,

Historical demography and population sizes
Neutrality and population size change tests for ND2 indicated significant signals of demographic expansion in clades NE2 and SE2 (R2 for clade NE2, D and R2 for SE2).For RAG1 D and Fs tests were significant for samples comprising clade NE2 and for β-fibint7 only the Fs test was significant for samples that correspond to clade SE2.We did not find any evidence for demographic expansion in 28S (S3 Table ).These results suggest that clades NE1 and SE1 did not experience events of demographic expansion.
The EBSP evidenced that mtDNA clades of S. eurydice (NE1, NE2, SE1 and SE2) experienced distinct demographic histories.Changes in effective population sizes (Ne) through time were observed only for SE1 and NE2.We evidenced a reduction of Ne in clade SE1, which started about 130,000 years ago (Fig 5), in the late Pleistocene.On other hand, clade NE2 presented a demographic expansion starting about 400,000 years ago.In addition, we observed population growth in clade SE2, beginning 500,000 years ago.
IMa2 four-population model estimates of current population size suggest SE2 being the largest population with around 1.3 million individuals.NE1 represents the smallest population

Discussion
Two different climatic domains based on congruent patterns of phylogeographic endemism for vertebrates have been recently proposed for the AF [5].These domains are associated with broadly different climatic regimes, which divide the AF in two large northern and southern  blocks, roughly at the Doce River.Our major break was found near the Jequitinhonha River, which reaches the Atlantic Ocean ~400 km north of the Doce River.Still, our sampling in the region between these two rivers is relatively poor, and because migrants from NE and SE populations across these rivers have been identified in other species studied [14], it is not possible to reject that the geographic separation between our clades lies actually near the Doce River.Nevertheless, other studies have also shown major breaks at the Jequitinhonha River valley for lizards [7], birds [8,16,52], and amphibians [4,14,53].At last, the Doce and Jequitinhonha rivers are relatively close (considering the 3200 km North-South distribution of the AF, over 29°in latitude), and because they are not the barrier per se, representing in fact a climatic turnover, it is likely that the position of phylogeographic breaks will vary among clades because of the different, clade-specific responses to the presence of rivers.
Although not necessarily being the sole cause of such breaks, rivers can act as secondary barriers, partially restricting gene flow among populations [54].This seems the case of the Doce River, which coincidently lies within a region of significant climatic turnover [5].Still, examples of phylogeographic breaks for amphibians and reptiles along the AF that match the position of current rivers abound (e.g., [12,14,15,53,55]).This partial restriction is likely to be clade-specific and may be the source of differences among the exact position of phylogeographic breaks between SE and NE vertebrate clades identified in the literature.The geographic structure observed between S. eurydice mitochondrial clades may therefore be increased by the action of rivers as secondary barriers, restricting gene flow and promoting allopatry between the two northern clades (São Francisco River) and between northern and southern clades (Jequitinhonha River).Nevertheless, the signal for these barriers is not pervasive for all markers, and no river can be pointed as a barrier between the two southern mitochondrial clades.
Indeed, nuclear and mitochondrial markers frequently show different phylogeographic resolutions for the same biogeographic history due to differences in mutation rates, recombination, and effective population sizes [12,56].Such differences are expected to occur more often in recently diverged populations or species.Accordingly, our nuclear markers supported the separation of NE and SE clades.Still, the two SE clades, although clearly distinct mitochondrially (average net distance of 10%, Table 1), were not fully recovered with nuclear genes.This may result from a secondary contact (with male biased migrations), incomplete lineage sorting, or both.The divergence observed for mitochondrial clades suggest that populations of a formally widespread species may have persisted in isolation for some time, likely in forest refugia and microrefugia (sensu [57]).
Microrefugia have been originally proposed to help explain the persistence in time of highmountain species and communities in the tepuis of Venezuela, and their biogeographic significance for the establishment of present day biotas is undisputable [57].In the AF, paleo modeling indicated two stable areas north of the Doce River, the Bahia and Pernambuco refuges, while areas south of the river have been predicted to be less stable [3].Although such models match phylogeographic data for several taxa [5], they also contradict studies that have found high phylogeographic endemism and diversity on the southern portions of the AF [12,13,21,58,59].Microrefugia may have contributed to the preservation of diversity at higher latitudes of the AF [5,13].In fact, the same region where clade SE1 was found (Rio de Janeiro) has been shown to harbor endemic lineages consistent with the existence of a microrefugium during the Pleistocene [13,60].
Pleistocene climatic fluctuations have been implicated in the diversification of several animal taxa in the AF [4,8,[16][17][18][19]61].Cladogenetic events for S. eurydice exhibit spatial and temporal concordance with previous works, which recovered a range of temporal divergences between the Pliocene and Pleistocene in the AF [15,52].Phylogeographic structure and divergence dates show that the separation between the mitochondrial clades of S. eurydice is consistent with allopatric diversification associated with Pliocene-Pleistocene climatic cycles.The two divergent clades in the Bahia and Pernambuco refuges that do not exchange migrants, the significant migration from NE to SE, and a highly divergent mitochondrial clade in high altitude areas at the southern portion of the AF support both refuges predicted north of the Doce River [3][4][5] and the existence of microrefugia at the southern portion of the AF [3,5,13].Hence, a once widely distributed S. eurydice may have had its distribution fragmented and restricted to microrefuges in the southern AF and to large refuges in Pernambuco and Bahia, with a recent re-establishment of contact by a directional migration from climatically-stable northeastern populations to the southeast and from lowland southeastern populations (SE2) to upland populations (SE1).
During the Last Glacial Maximum rainforests were restricted to mountainous areas and high altitude species were less affected by cold and dry climates of Pleistocene glaciations [55].Lowland species, in contrast, should have benefited from warmer climates and present signs of demographic expansions following the end of glaciations.Because southern populations of S. eurydice occupy less stable portions of the AF [3], they should present more prominent signals of demographic alterations.Indeed, SE1 and SE2 were the only populations showing significant contractions/expansions.However, SE1 presented a significant contraction while SE2 showed signals of demographic expansion.This contrasting result may reflect different adaptive niches for each population (and therefore different responses to the end of glaciations).Alternatively, the demographic expansion detected for SE2 may result from the high levels of migration from NE populations (which violate EBSP assumptions), the expansion of suitable lowland forested areas in warmer periods, or both.Dates recovered for these demographic fluctuations suggest a recent effect of temperature increasing for SE1 (~150-200k years ago), but a much older expansion for SE2 (~500-600k years ago).Nevertheless, demographic expansions within the same timeframe have been identified in this region of the AF for several other taxa, including frogs [14] and birds [18] suggesting the pervasive effects of historical events affecting several taxa in the region.Such scenario of different demographic responses along with high mitochondrial differentiation (lowest Da value of 8%), wide geographic distribution, and lack of migration among NE populations suggest that S. eurydice may actually represent more than one species.Future researches with S. eurydice populations using an integrative approach and inferences of hybridization in secondary contact zones could be used to describe and delimit the spatial patterns of variation in their geographic range in Atlantic Forest and provide a more accurate hypothesis for its taxonomy.

Conclusions
Scinax eurydice is genetically structured across the Atlantic Forest, showing high mtDNA divergence and low to inexistent levels of migration among populations.Genetic patterns and demographic data match predictions for the existence of large Refugia in the Northern AF (Bahia and Pernambuco), as well as microrefugia south of the Doce/Jequitinhonha Rivers biogeographic division.Patterns of gene flow, demography, and divergence dates agree with a scenario of recent demographic expansion of lowland taxa, likely favored by the end of major glaciation events and the expansion of suitable habitats.Species tree, demographic parameters, temporal and spatial population structures recovered endorse S. eurydice as a complex of species.Two rivers lie at the boundaries among populations (Jequitinhonha and São Francisco), a pattern repeatedly found in the AF, which suggests that the role of rivers as secondary barriers must be re-evaluated.The current distribution of these populations and their origin are linked to the recent history of the AF, underscoring the importance of taxonomic reappraisals based on multiple lines of evidences to ensure the preservation not only of species, but of the evolutionary processes responsible for the origin of current biodiversity.

Fig 1 .
Fig 1. Localities sampled for mtDNA of Scinax eurydice and mitochondrial gene tree.The gene tree was obtained by Maximum Likelihood and Bayesian inference for the ND2 locus with support values (Bootstrap/Bayesian Posterior Probabilities) shown by major nodes.The blue circles, green triangles, red hexagons and yellow squares on the map indicate the locations sampled (mtDNA ND2) for each geographic group identified.The area shaded in dark gray on the Brazilian map represents states sampled for present work.The area in light gray in the inset represents the Atlantic Rainforest (based on SOS Mata Atlântica 2013).Haplotypes (H) represented on the tree topology are available in S1 Table.NE (Northeastern) and SE (Southeastern) correspond to populations identified.Numbers inside circles represent rivers listed as possible secondary barrier.doi:10.1371/journal.pone.0154626.g001

For
nuclear genes alone, the most likely number of populations according to Evanno's method was 2 (see S1 Fig in Supporting Information), with a second peak for k = 4. Structure tends to represent the most basal hierarchical structure in the dataset, and the second peak may represent a more refined structure of the dataset [38].Indeed, Bayes rule [Pr(k)] recovered 4 as the best number of k (S3 Table).k = 2 shows a strong separation between northeastern (NE1 + NE2) and southeastern (SE1 + SE2) clades consistent with the basal separation recovered in the mitochondrial phylogeny (Fig 1).Q graphs for k = 4 do not match the mitochondrial phylogeny (Fig 1) and also exhibits admixture possibly due to gene flow (Fig 3, see Migration below).

Tamura-
Nei distance (Da) within and between clades recognized by genes trees (ND2).doi:10.1371/journal.pone.0154626.t003with around 170 thousand individuals.NE2 and SE1 have intermediate sizes with around 700 thousand individuals (S5 Fig).Two-population model estimates of current population sizes overlap four-population model estimates, but peak estimates suggest larger populations.This is expected if migration from a third population is not taken into account [48,51].Ancestral population size estimates for two-population models suggest an ancestral population of 250 thousand individuals for NE1 and NE2 populations, and around 2 million individuals for SE1 and SE2 (S6 Fig).The same ancestral populations are poorly estimated in the four-population model, posterior distributions are rather flat with the ancestral of NE1 and NE2 having the highest probability at zero.The population size estimate of the ancestral of all for populations shows a peak at around 2.1 million individuals (S6 Fig).

Fig 4 .Table 4 . 3 )
Fig 4. Chronogram of species tree from the combined mitochondrial and nuclear DNA data generated on *BEAST.Numbers by nodes represent clade posterior probabilities.Areas where many species trees agree on topology and/or branch lengths are densely colored in blue.The most frequent topologies are represented by shaded consensus tree.White circles with letters at nodes are divergence times and 95% (HPD) of high posterior density of divergence times between clades according to Table 4.doi:10.1371/journal.pone.0154626.g004

Fig 5 .
Fig 5. Extended Bayesian Skyline Plot (EBSP) results for populations of Scinax eurydice.The dotted horizontal line shows the median estimate of the EBSP and the gray area shows the upper and lower 95% highest posterior density (HPD).Each graph represents the EBSP analysis for: (A) clade NE1; (B) clade NE2; (C) SE1; (D) clade SE2.The dotted vertical line illustrates the approximate time when events of demographic expansion started in each clade of S. eurydice.doi:10.1371/journal.pone.0154626.g005

Table 3 .
Diversity indexes for each fragment and clade.
H, number of haplotypes; Hd, percentage of haplotype diversity and respective standard deviations (SD); π, percentage of nucleotide diversity and respective standard deviations; N, number of sequences; D, Tajima test; R2, Ramos-Onsins test; Fs, Fu test; np, no polymorphism; ns, no significant value; *, P < 0.05.