SciELO - Scientific Electronic Library Online

 
vol.44 issue172Study of the magnetic semiconductor of the Fe-doped ZnO system obtained by mechanical alloyOrtalis columbiana population density in a sub-Andean forest fragment in Valle del Cauca, Colombia author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Revista de la Academia Colombiana de Ciencias Exactas, Físicas y Naturales

Print version ISSN 0370-3908

Rev. acad. colomb. cienc. exact. fis. nat. vol.44 no.172 Bogotá July/Sept. 2020  Epub July 04, 2021

https://doi.org/10.18257/raccefyn.956 

Natural sciences

The genetic structure of Pristimantis latro. (Anura: Craugastoridae) mirrors traits of their life history

La estructura genética de Pristimantis latro. (Anura: Craugastoridae) refleja rasgos de su historia de vida

Gilcilene Santana-Cornélio1 

Elciomar Araújo-de-Oliveira2 

Keila Magalhães-Xavier3 

Gabriela Wemilly Barros-da-Silva1 

Isadora França1 

Luis Reginaldo Ribeiro Rodrigues4 

Emil José Hernández-Ruz5  * 

1 Faculdade de Ciências Biológicas, Campus Universitário de Altamira Universidade Federal do Pará, Altamira, Pará, Brazil

2 Programa de Pós-Graduação em Biodiversidade e Biotecnologia da Rede BIONORTE, Universidade Federal do Amazonas, Manaus, Amazonas, Brazil

3 Faculdade de Medicina, Campus Universitário de Altamira, Universidade Federal do Pará, Altamira, Pará, Brazil

4 Programa de Pós-Graduação em Recursos Naturais da Amazônia, Universidade Federal do Oeste do Pará, Santarém, Pará, Brazil; Laboratório de Genética & Biodiversidade, Instituto de Ciências da Educação, Universidade Federal do Oeste do Pará, Santarém, Pará, Brazil

5 Programa de Pós-Graduação em Biodiversidade e Conservação, Campus Universitário de Altamira, Universidade Federal do Pará, Altamira, Pará, Brazil


Abstract

One of the main hypotheses to explain the origin of Amazonian diversity is the barrier effect of the rivers known as the riverine-barrier hypothesis, which suggests that riverine barriers isolated once continuous populations leading to differentiation and speciation. In this context, we studied the genetic structure of Pristimantis latro, a newly described species that inhabits a region under marked anthropic pressure due to expansive livestock, illegal mining, and hydroelectric dam construction. The DNA was extracted from 52 P. latro individuals and then amplified via polymerase chain reaction (PCR) using the mitochondrial 16S rRNA and COI markers. To infer the time of divergence between the P. latro localities, we built a species tree and performed an analysis of molecular variance (AMOVA) to infer the genetic differentiation between and within the P. latro populations. We found that P. latro has a marked genetic structure in the populations of the right and left margins of the Xingu River and within the Tapajós-Xingu and Xingu-Tocantins interfluvial regions and that the time of divergence between the populations of the East and West banks of the Xingu River occurred approximately 380,000 years ago. This pattern of genetic structure corresponds to that reported in recent articles for the Pristimantis genus evidencing that species without tadpoles exhibit a genetic structure explained by the hypothesis of rivers as barriers.

Keywords: Mitochondrial DNA; Isolation by distance; Rivers as barrier; Terrarana; Divergence times

Resumen

Una de las principales hipótesis para explicar el origen de la diversidad amazónica invoca el efecto de barrera de los ríos para explicar los patrones de diversidad. Esa hipótesis propone que algunos ríos pueden separar poblaciones continuas conduciendo a la diferenciación y la especiación. En ese sentido nos propusimos estudiar la estructura genética de Pristimantis latro, una especie recién descrita que habita una región bajo fuerte presión antrópica producto de la ganadería expansiva, la minería ilegal y la construcción de represas hidroeléctricas. El ADN se extrajo de 52 individuos de P. latro y se amplificó mediante reacción en cadena de la polimerasa (PCR) usando los marcadores mitocondriales 16S rRNA y COI. Para inferir el tiempo de divergencia entre las localidades de P. latro, construimos un árbol de especies e hicimos análisis de varianza molecular (AMOVA) para inferir la diferenciación genética entre y dentro de las poblaciones de P. latro. Encontramos que P. latro presenta una estructuración genética en las poblaciones de las márgenes derecha e izquierda del río Xingu y dentro de las regiones interfluviales deTapajós-Xingu y Xingu-Tocantins. Asimismo, el tiempo de divergencia entre las poblaciones de las márgenes derecha e izquierda del río Xingu aconteció hace aproximadamente 380.000 años. El patrón de estructura genética que se encontró se corresponde con el indicado en artículos recientes para el género Pristimantis, el cual revela que las especies sin renacuajos exhiben una estructura genética que responde a la hipótesis de los ríos como barreras.

Palabras clave: ADN mitocondrial; Aislamiento por distancia; Ríos como barrera; Terrarana; Tiempo de divergencia

Introduction

One of the premises about population genetics predicts that dispersion decreases with an increased geographic distance where gene flow is generally lower among geographically distant populations causing a pattern of isolation by distance (Wright, 1943). It is well known that, besides dispersion, there are vicariant events (river and mountain formation) that explain biologic diversification.

Many hypotheses have been put forward to explain the high levels of diversity in Amazonia (Da Rocha & Kaefer, 2019). One of them, known as the riverine-barrier hypothesis, invokes the barrier effect of the rivers to explain patterns of diversity (Wallace, 1852; Bates, 1874) and proposes that riverine barriers separate once continuous populations leading to differentiation and speciation. However, we know little about the genetic structure of tropical anuran populations, especially those groups with cryptic diversity in tropical regions due to the marked phenotypic conservatism present between species as compared to the marked structuring and high genetic differentiation observed between different populations (Vieites, et al., 2009; Kaefer, et al., 2013; Guayasamin, et al., 2017).

A good example of this effect on anurans is the genus Pristimantis Jiménez De la Espada 1870, the most diverse among all vertebrate groups (Fouquet, et al., 2013). This diversity may be associated with the evolution of the direct development characteristic of Terrarana, which makes them interesting for evolution studies since species with such traits have poor dispersion capacity and philopatry (Beebee & Griffiths, 2005; Fouquet, et al., 2012). As a consequence of their life history, populations tend to be highly genetically structured and retain marked signs of past evolutionary processes making them good models for testing evolutionary hypotheses (Funk, et al., 2012; Fouquet, et al., 2012).

Pristimantis latro Oliveira, Rodrigues, Kaefer, Pinto, and Hernández-Ruz 2017 is found in the eastern Amazon between the Tapajós and the Tocantins- Araguaia Rivers occupying both banks of the Xingu River (Oliveira, et al., 2017). In the region, primary forests have a high degree of fragmentation and degradation with the substitution of the land cover mainly by extensive livestock and the advance of soy crops (Almeida, et al., 2014). Given the increase in environmental threats and habitat reduction on the newly described species and the lack of knowledge about the genetic structure of P. latro, we used fragments of the mitochondrial 16s rRNA and COI genes to: 1) Evaluate patterns of genetic structure; 2) test whether there is a relationship between the genetic distance and the geographic distance of different localities; 3) to reconstruct the phylogenetic relationships among the populations, and 4) to estimate the time of divergence between the sampled localities.

Methodology

Study area and sampling

Fifty-two individuals of P. latro were collected in eight municipalities in the state of Pará, Brazil: Santarém, Medicilândia, Brasil Novo, Vitória do Xingu, Altamira, Anapu, Senador José Porfírio, and Marabá (Figure 1). The forest composition in the study area varies from dense ombrophylous forest to lowland ombrophylous forest (Instituto Brasileiro de Geografia e Estatistica - IBGE, 2012) with immense pasture areas.

Figure 1 Map of Pristimantis latro capture localities. The red circles represent the localities from the left bank of the Xingu River and the green ones the localities from the right bank. The numbers correspond to the following sampled localities: 1) Santarém; 2) Medicilândia; 3) Brasil Novo; 4) Altamira; 5) Vitória do Xingu; 6) Senador José Porfírio; 7) Anapu (Type locality), and 8) Marabá. 

The specimens collected in the field were euthanized with topical lidocaine ointment (benzotop®). The muscle or liver tissues were removed and stored in 95% ethanol. Specimens were fixed in 10% formalin, stored in 70% alcohol, and deposited in the Zoology Laboratory Collection- Adriano Giorgi of the Faculty of Biological Sciences at the Federal University of Pará - UFPA, Altamira Campus.

DNA extraction, amplification, and sequencing

The total genomic DNA was extracted with the CTAB protocol (Doyle & Doyle, 1987) and resuspended with 30 μL of ultrapure water. DNA quality was verified using 1% agarose gel in electrophoresis. Polymerase chain reaction (PCR) was performed to amplify a 600 bp fragment of the mitochondrial 16S rRNA gene and a 600 bp fragment of the mitochondrial cytochrome oxidase subunit I (COI) gene.

The amplification reaction of the mitochondrial 16S rRNA gene fragment was done under the following conditions: An initial temperature of 92°C (60 sec.), followed by 35 cycles of 92°C (60 sec), 50°C (50 sec) and 72°C (1.5 min) and a final extension at 72°C for 7 min. For the amplification of the mitochondrial COI gene fragment, we used the following conditions: 95°C (3 min), followed by 35 cycles of 95°C (of 30 sec), 53°C (30 sec), and 72°C (40 sec), and a final extension of 72°C (7 min).

The sequencing reactions were performed using the BigDye™ Terminator kit (Life Technology) following the manufacturer's protocol and using the forward primers 16S A (mtDNA) (Palumbi, 1991) and COI LCO 1590 (Folmer, et al., 1989). The products were precipitated with ethanol according to the manufacturer's recommendations (Applied Biosystems™), resuspended in 10 μL of formamide (Thermo Scientific™), and subsequently injected into the automatic ABI sequencer 3500XL (Applied Biosystems™).

Alignment of sequences and molecular dating

We aligned the sequences using the BioEdit software (Hall, 1999) using the ClustalW algorithm (Thompson, et al., 1996) and edited manually (checking the spectrograms for errors). To infer the time of divergence between the P. íatro localities, we built a species tree using the software BEAST v1.8 (Drummond, et al., 2012) using both the 16S and COI genes with the GTR + G evolutionary model chosen in the jModelTest software (Darriba et al., 2012) via corrected Akaike information criterion (AICc) and using a coalescent: constant size as tree prior. The species of Pristimantis used as outgroup were described by Oliveira, et al. (2020). The total number of sequences used for each locality in this study are shown in table 1. Given the absence of a fossil record, we used a general mutation rate for the 16S of 0.0028 per lineages per millions of years (Lymberakis, et al., 2007; Lemmon, et al., 2008) and a fixed substitution rate of 1.0 to COI mtDNA according to Fourdrilis, et al. (2016). The chain (MCMC) had 500 million generations with a tree sampled every one thousand generations. Thus, we recorded 50.000 trees and discarded the first 10.000 considering a 20% burn-in.

Table 1 Sampling localities of Pristimantis latro and respective numbers of mitochondrial genes (16S rRNA and COI) samples 

We performed such analyses assuming a lognormal relaxed molecular clock, which presumes independent rates in different branches (Drummond, et al., 2006). The previous trees were modeled according to the Yule speciation process. The burn-in value was selected by viewing the log likelihoods associated with the distribution of the tree with the Tracer v1.5 software (Drummond & Rambaut, 2007), discarding 20% of the trees with the TreeAnnotator v1.6.6 (Drummond & Rambaut, 2007).

We obtained the uncorrected pairwise genetic distances p between localities in the program MEGA 6.0 (Tamura, et al., 2013) using thep distance for the 16S rRNA and the Kimura 2 parameter (K2P) for the COI.

Population structure analysis

We only used the COI gene as it had the largest number of sequenced locations compared to the 16S. We implemented a Bayesian Analysis of Population Structure (BAPS) in the BAPS 5.0 software (Corander, et al., 2008) to find clusters formed with the COI gene sequences (eight localities). The BAPS software uses nucleotide frequencies of the samples to infer the K number of genetically different groups through Bayesian analysis allocating similar sequences in the same group. The maximum number of K chosen for the analysis was eight, which corresponded to each locality. We used the log-likelihood value to choose the most likely grouping configuration. We performed an analysis of molecular variance (AMOVA) to infer the genetic differentiation between and within the P. íatro populations using the Arlequin 3.0 program (Excoffier, et al., 2005).

To estimate the genealogical relationships among individuals, we produced a network of haplotypes based on the topology of a phylogenetic tree for the COI gene built with the TreeFinder software (Jobb, 2011) with 10,000 non-parametric bootstrap replicas (Felsentein, 1985). We built the haplotypes network using the HaploView software (Salzburger, et al., 2011) to determine the number of haplotypes shared based on likelihood calculations. We did not use 16S gene sequences in the BAPS analysis and haplotypes network.

To calculate the relative importance of isolation by distance, we used the simple Mantel test performed in the Alleles In Space (AIS) program (Miller, 2005) for the COI marker only as it presented more localities. The Mantel test involves estimating the correlation between two matrices through Pearson's correlation r coefficient. The significance of r values was estimated by 1.000 permutations.

Results

Molecular dating and genetic distance

We obtained 52 sequences of the 16S rRNA gene fragment corresponding to eight haplotypes with 521 bp, 507 preserved sites, 14 variable sites, and 11 informative sites, as well as 41 COI DNAmt sequences equivalent to 14 haplotypes with 617 bp, 579 preserved sites, 38 variable sites, and 31 informative sites. Within P latro, the first cladogenesis event occurred around 0.860 Ma (95% HPD: 0.26 - 1.87 Ma), which separated the populations of Santarém (STM) from the others of P. latro. The second cladogenesis event was between the east and west margins around 0.380 Ma (95% HPD: 0.16 - 0.67 Ma) (Figure 2).Table 2 shows all the divergence times and the HPD interval.

Figure 2 Tree of P. latro species based on the Bayesian analysis for the 16S and COI genes. The values of the branches indicate the posterior probability. The scale bar shows the number of replacements per site. Abbreviations represent the localities where P latro was studied in the Xingu region: Municipality of Santarém (STM), Municipality of Medicilandia (MED), Brazil Novo (BRN), Altamira (ATM), Vitória do Xingu (VTX), Senador José Porfirio (SEN), and Anapu (ANA), those for Pristimantis giorgii: (FUP, Fazenda Uberlândia Portel - PA; MAR, Marabá - PA, and CAX, Caxiuanã - PA), forP. pictus: (JAC, Jacareacanga; PAR, Paranaíta - MT; COT, Cotriguaçu - MT; NMU, New World - MT; CRI, Cristalino - MT), and Pristimantis sp. (TER, Terra do Meio - PA). 

Table 2 Divergence time between populations of Pristimantis latro. Terra do Meio (TER); Santarém (STM); Medicilândia (MED); Brasil Novo (BRN); Altamira (ATM); Vitória do Xingu (VTX); Senador José Porfirio (SEN), and Anapu (ANA) 

The genetic distance within P. latro presented a variation of 0.1 - 8.3% for the COI gene (Supplementary material 1,https://www.raccefyn.co/index.php/raccefyn/article/view/956/2807) and 0 - 1.7% for the 16S rRNA gene (Supplementary material 2,https://www.raccefyn.co/index.php/raccefyn/article/view/956/2808). Between the left and right banks, the distance was 2.8 - 4.6% for COI and 0 - 1.5% for 16S. The greatest genetic distance (8.3% for COI) was between P. latro from the Marabá and Senador José Porfirio municipalities with a geographic distance of 443 Km, whereas for 16S (1.7%), it was between P. latro from the Brasil Novo + Medicilândia and Santarém municipalities with a geographic distance of 240 Km. The most geographically distant localities (Santarém and Marabá, located on the opposite banks of the Xingu River, 701 Km), showed a genetic distance of 5.1% for the COI. Due to the absence of sequences for the 16S gene from Marabá, this comparison could not be made, therefore, we compared the Santarém and Anapu municipalities (opposite banks of the Xingu River) ~377 Km apart, which presented a 1% genetic distance between one another.

Population structure

Pristimantis latro presented a marked genetic structure for the COI gene (logML = -382.8427; probability = 0.99) forming five clusters (Figure 3): Cluster 1 (Santarém); Cluster 2 (Medicilândia - A); Cluster 3 (Medicilândia-B, Brasil Novo, Altamira, and Vitória do Xingu); Cluster 4 (Senador José Porfírio-A and Anapu - A); and Cluster 5 (Anapu-B, Senador José Porfírio-B, and Marabá). In the haplotypes network (Figure 4), we observed no shared of haplotype between the east and west banks of the Xingu River corroborating the result of the BAPS.

Figure 3 Population structure of P. latro from the COI mtDNA data. The vertically arranged colors represent each of the genetic groups (clusters). The numbers correspond to the sampled localities: 1) Santarém; 2) Medicilândia; 3) Brasil Novo; 4) Altamira; 5) Vitória do Xingu; 6) Senador José Porfírio; 7) Anapu, and 8) Marabá. 

Figure 4. Median-joining haplotypes network based on DNAmt COI between populations of Pristimantis latro. Each haplotype is represented by a circle where the area is proportional to its frequency. The traits indicate additional mutational steps of branches with more than one mutation. The colors indicate each location. Blue dots represent lost or unsampled haplotypes. The abbreviations correspond to the sampled localities: Santarém (STM); Medicilândia (MED); Brasil Novo (BRN); Altamira (ATM); Vitória do Xingu (VTX); Senador José Porfírio (SEN); Anapu (ANA), and Marabá (MAR). 

The result of the AMOVA for the COI gene revealed that most of the genetic differentiation is distributed among the five clusters resulting from the haplotypes network and the BAPS (79.51%) while a smaller part within them (20.49%) (ΦST = 0.79; p=0.00) (Table 3) showed a marked population structure in P. latro. The Mantel test revealed a positive and significant correlation between geographic and genetic distances (P<0.001; r=0.558) indicating a marked role of random events in the reproductive isolation between populations and genetic drifting in the process of evolutionary diversification of the studied marker.

Table 3 Analysis of molecular variance (AMOVA) estimated to evaluate the genetic differentiation between and within the five clusters of Pristimantis latro 

Discussion

Pristimantis latro showed a strong genetic structure between the populations of the left and right bank of the Xingu River, as well as in the same bank. Anurans with direct development have life strategies that promote high population structure with species dispersing only by land where habitat heterogeneity inhibits dispersal across the landscape (Voss, et al, 2001).

The genetic structure pattern found in the present study is similar to those described in other Amazon frogs by Fouquet, et al. (2015) where species without tadpoles exhibited a genetic structure associated with rivers as barriers. Also consistent with our results are those of the report by Guayasamin, et al. (2017) showing that color patterns in Pristimantis ornatissimus can be delimited by geographical barriers (elevation gradients, rivers). However, we could not make a thorough comparison with other species of the genus Pristimantis since most of the published works refer to Andean species while P. latro is found in lowlands.

We found a spatial pattern of distribution of genetic variability similar to that proposed in the rivers-as-barrier hypothesis by Wallace (1852). Such a pattern has been found in several studies of different terrestrial vertebrates (Antonelli, et al., 2010). Finally, the sample design used in the present study did not evaluate if there was "permeability" of the rivers to gene flow in the respective extensions. Thus, our results add to a growing body of knowledge indicating that genetic variability at intra and interspecific levels in amphibians is related to the transposition of large Amazonian rivers (Funk, et al., 2007; Fouquet, et al., 2012; Kaefer, et al., 2012; 2013; Simões, et al., 2014).

The municipalities of Medicilândia (left bank of the Xingu River) and Senador José Porfírio and Anapu (right bank) presented more than one cluster per locality supporting local genetic structuring. Genetic structuring in amphibians at distances less than 5 km has already been documented (Lampert, et al., 2003; Andersen, et al., 2004; Burns, et al., 2004) showing that limited gene flow can promote differentiation in local populations mainly due to habitat fragmentation. The Mantel test explains that 31% (r2) of the variation is caused by geographic distance, but is not responsible for all structuring, which indicates that other factors are associated with this pronounced structuring (rivers).

The separation between the left and right banks of the Xingu River (middle of the Pleistocene) is relatively recent, however, it played an important role in population structuring (BAPS in figure 3) and the absence of haplotypes shared between the banks (haplotypes network in figure 4).Ribas, et al. (2012) suggested that the Xingu River acted a geographic barrier separating Psophia interjecta and P. dextralis around 0.8 -0.3 Ma and was responsible for the separation of younger lineages. Recent rivers may represent recent or semipermeable barriers, as found for Amphibia and Squamata in the Tapajós River (Moraes, et al., 2016).

Conclusions

Pristimantis latro presents genetic structure between the left and right banks of the Xingu River, as well as within the Tapajós-Xingu and Xingu-Tocantins interfluvial regions. The pattern of genetic structure found here corresponds to that described in recent articles for the genus Pristimantis, which reveals that species without tadpoles exhibit a genetic structure explained by the hypothesis of rivers as barriers.

Acknowledgments

We thank Faculdade de Ciências Biológicas Zoology Laboratory at the Universidade Federal do Pará, Campus Universitário de Altamira, the three anonymous reviewers for their contributions, and Dr. Jesus Rivas for reviewing the English of the abstract.

References

Almeida, A.S. de, Vieira, I.C.G., Barros, M.N.R., Rocha, D.P.N. (2014). Área de endemismo Belém e Xingu: configuração e espacialização do uso da terra e da cobertura vegetal. In: Thaise Emilio; Flávio Luizão. (Org.). Cenários para Amazônia: clima, biodiversidade e uso da terra. 1ed.Manaus: INPA, p. 57-66. [ Links ]

Andersen, L.W., Fog, K., Damgaard, C. (2004). Habitat fragmentation causes bottlenecks and inbreeding in the European tree frog (Hyla arborea). Proceedings of the Royal Society B: Biological Sciences. 271: 1293-1302. Doi: 10.1098/rspb.2004.2720 [ Links ]

Antonelli, A., Quijada-Mascareñas, A., Crawford, A.J., Bates, J.M., Velazco, P.M., Wuster, W. (2010). Molecular studies and phylogeography of Amazonian tetrapods and their relation to geological and climatic models. Molecular studies of Amazonian tetrapods. En C. Hoorn and F. P. Wesselingh (Ed), Amazonia: Landscape and Species Evolution: A look into the past (pp. 385-404). Wiley-Blackwell [ Links ]

Bates, H. (1874). The naturalist on the River Amazon. John Murray, London. 506 pp. [ Links ]

Beebee, T.J.C., Griffiths, R.A. (2005). The amphibian decline crisis: A watershed for conservation biology? Biological Conservation. 125: 271-285. Doi: 10.1016/j.biocon.2005.04.009 [ Links ]

Bogart, J. P. (1991). The influence of life history on karyotypic evolution in frogs. Amphibian Cytogenetics and Evolution. Academic Press, San Diego. p. 233-258. [ Links ]

Burns, E.L., Eldridge, M.D.B., Houlden, B.A. (2004). Microsatellite variation and population structure in a declining Australian Hylid Litoria aurea. Molecular Ecology. 13: 1745-1757. Doi: 10.1111/j.1365-294X.2004.02190.x [ Links ]

Corander, J., Marttinen, P., Sirén, J., Tang, J. (2008). Enhanced Bayesian modelling in BAPS software for learning genetic structures of populations. BMC bioinformatics. 9: 1-14. Doi: 10.1186/1471-2105-9-539 [ Links ]

Da Rocha, D.G., Kaefer, I.L. (2019). What has become of the refugia hypothesis to explain biological diversity in Amazonia? Ecology and Evolution. 9: 4302-4309. Doi:10.1002/ece3.5051 [ Links ]

Darriba, D., Taboada, G.L., Doallo, R., Posada, D. (2012). JModelTest 2: More models, new heuristics and parallel computing. Nature Methods. 9: 772. Doi: 10.1038/nmeth.2109 [ Links ]

Doyle, J.J., Doyle, J.L. (1987). Isolation of plant DNA from fresh tissue. Focus. 12 (1): 13-15. [ Links ]

Drummond, A.J., Ho, S.Y.W., Phillips, M.J., Rambaut, A. (2006). Relaxed phylogenetics and dating with confidence. PLoS Biology. 4: 699-710. Doi: 10.1371/journal.pbio.0040088 [ Links ]

Drummond, A.J., Rambaut, A. (2007). BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evolutionary Biology. 7: 1-8. Doi: 10.1186/1471-2148-7-214 [ Links ]

Drummond, A, J., Suchard, M.A., Xie, D., Rambaut, A. (2012). Bayesian Phylogenetics with BEAUti and the BEAST 1.7. Molecular Biology and Evolution. 29: 1969-1973. Doi: 10.1093/molbev/mss075 [ Links ]

Excoffier, L., Laval, G., Schneider, S. (2005). Arlequin 3.0: An integrated software package for population genetics data analysis. Evolution Bioinformatics Online. 1: 47-50. [ Links ]

Felsenstein, J. (1985). Confidence Limits on Phylogenies: An Approach Using the Bootstrap. Evolution. 39: 783-791. Doi: 10.2307/2408678 [ Links ]

Folmer, O., Black, M., Hoeh, W., Lutz, R., Vrijenhoek, R. (1994). DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Molecular Marine Biology Biotechnology. 3: 294-299. [ Links ]

Fouquet, A., Ledoux, J.B., Dubut, V., Noonan, B.P., Scotti, I. (2012). The interplay of dispersal limitation, rivers, and historical events shapes the genetic structure of an Amazonian frog. Biological Journal of the Linnean Society. 106: 356-373. Doi: 10.1111/j.1095-8312.2012.01871.x [ Links ]

Fouquet, A., Martínez, Q., Courtois, E.A., Dewynter, M., Pineau, K., Gaucher, P., Blanc, M., Marty, C., Kok, P.J.R. (2013). A new species of the genus Pristimantis (Amphibia, Craugastoridae) associated with the moderately evelated massifs of French Guiana. Zootaxa. 3750: 569-586. Doi: 10.11646/zootaxa.3750.5.8 [ Links ]

Fouquet, A., Courtois, E.A., Baudain, D., Lima, J.D., Souza, S.M., Noonan, B.P., Rodrigues, M.T. (2015). The trans-riverine genetic structure of 28 Amazonian frog species is dependent on life history. Journal of Tropical Ecology. 31: 361-373. Doi: 10.1017/S0266467415000206 [ Links ]

Fourdrilis, S., Mardulyn, P., Hardy, O.J., Jordaens, K., Freitas Martins, A.M. de, Backeljau, T. (2016). Mitochondrial DNA hyperdiversity and its potential causes in the marine periwinkle Melarhaphe neritoides (Mollusca: Gastropoda). PeerJ. 4: e2549. Doi: 10.7717/peerj.2549 [ Links ]

Funk, W.C., Caldwell, J.P., Peden, C.E., Padial, J.M., De la Riva, I., Cannatella, D.C. (2007). Tests of biogeographic hypotheses for diversification in the Amazonian forest frog, Physalaemuspetersi. Molecular Phylogenetics and Evolution. 44: 825-837. Doi: 10.1016/j.ympev.2007.01.012 [ Links ]

Funk, W.C., Caminer, M., Ron, S.R. (2012). High levels of cryptic species diversity uncovered in Amazonian frogs. Proceedings of the Royal Society of London. Series B, Biological sciences. 279: 1806-1814. Doi: 10.1098/rspb.2011.1653 [ Links ]

Guayasamin, J.M., Hutter, C.R., Tapia, E.E., Culebras, J., Pyron, R.A., Morochz, C., Funk, W.C., Arteaga, A. (2017). Diversification of the rainfrog Pristimantis ornatissimus in the lowlands and Andean foothills of Ecuador. PLoS ONE. 12: 1-21. [ Links ]

Hall, T.A. (1999). BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium: 95-98. [ Links ]

Heinicke, M.P., Duellman, W.E., Trueb, L., Means, D.B., Macculloch, R.D., Hedges, S.B. (2009). A new frog family (Anura: Terrarana) from South America and an expanded direct-developing clade revealed by molecular phylogeny. Zootaxa. 2111: 1-35. Doi: 10.11646/zootaxa.2211.1.1 [ Links ]

Instituto Brasileiro de Geografia e Estatistica - IBGE (2012). Manual técnico da vegetação brasileira. Rio de Janeiro. 271p. [ Links ]

Jobb, G., Haeseler, A. V., Strimmer, K. (2011). TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics. BMC Evolutionary Biology, 4:18. doi: 10.1186/1471-2148-4-18. [ Links ]

Kaefer, I.L., Kaefer, I.L., Tsuji-Nishikido, B.M., Lima, A.P. (2012). Beyond the river: Underlying determinants of population acoustic signal variability in Amazonian direct-developing Allobates (Anura: Dendrobatoidea). Acta Ethologica. 15: 187-194. Doi: 10.1007/s10211-012-0126-0 [ Links ]

Kaefer, I.L., Tsuji-Nishikido, B.M., Mota, E.P., Farias, I.P., Lima, A.P. (2013). The Early Stages of Speciation in Amazonian Forest Frogs: Phenotypic Conservatism despite Strong Genetic Structure. Evolutionary Biology. 40: 228-245. Doi: 10.1007/s11692-012-9205-4. [ Links ]

Lampert, K.P., Rand, A.S., Mueller, U.G., Ryan, M.J. (2003). Fine-scale genetic pattern and evidence for sex-biased dispersal in the tungara frog, Physalaemus pustulosus. Molecular Ecology. 12: 3325-3334. Doi: 10.1046/j.1365-294X.2003.02016.x [ Links ]

Lemmon, E.M., Lemmon, A.R., Collins, J.T., Cannatella, D.C. (2008). A new North American chorus frog species (Amphibia: Hylidae: Pseudacris) from the south-central United States. Zootaxa. 1675: 1-30. [ Links ]

Lymberakis, P., Poulakakis, N., Manthalou, G., Tsigenopoulos, C.S., Magoulas, A., Mylonas, M. (2007). Mitochondrial phylogeography of Rana (Pelophylax) populations in the Eastern Mediterranean region. Molecular Phylogenetics and Evolution. 44: 115-125. Doi: 10.1016/j.ympev.2007.03.009 [ Links ]

Palumbi, S.R., Martin, A., Romano, S., McMillan, W.O., Stice, L., Grabowski, G. (1991) 'The Simple Fool's Guide to PCR, Version 2.0.'' Privately published document compiled by S. Palumbi, Dept. Zoology, Univ. Hawaii. [ Links ]

Miller, M.P. (2005). Alleles In Space (AIS): Computer Software for the Joint Analysis of Interindividual Spatial and Genetic Information. Journal of Heredity. 96: 722-724. Doi: 10.1093/jhered/esi119 [ Links ]

Moraes, L.J.C.L., Pavan, D., Barros, M.C., Ribas, C.C. (2016). The combined influence of riverine barriers and flooding gradients on biogeographical patterns for amphibians and squamates in south-eastern Amazonia. Journal of Biogeography. 43: 2113-2124. Doi: 10.1111/jbi.12756 [ Links ]

Oliveira, E.A., Rodrigues, L.R., Kaefer, I.L., Pinto, K.C., Hernández-Ruz, E.J. (2017). A new species of Pristimantis from eastern Brazilian Amazonia (Anura: Craugastoridae). ZooKeys. 2017: 101-129. Doi: 10.3897/zookeys.687.13221 [ Links ]

Oliveira, E.A., Silva, L.A., Silva, E.A.P., Guimaraes, K.L.A., Penhacek, M.P., Martínez, J.G., Rodrigues, L.R.R., Santana, D.J., Hernández-Ruz, E.J. (2020). Four new species of Pristimantis Jiménez de la Espada, 1870 (Anura: Craugastoridae) in the eastern Amazon. PLoS One. 15 (3): e0229971. [ Links ]

Ribas, C.C., Aleixo, A., Nogueira, A.C.R., Miyaki, C.Y., Cracraft, J., Andre, A. (2012). A palaeobiogeographic model for biotic diversification within Amazonia over the past three million years. Proceedings of the Royal Society of London. Series B, Biological sciences. 279: 681-689. Doi: 10.1098/rspb.2011.1120 [ Links ]

Rodríguez, A., Bomer, M., Pabijan, M., Gehara, M., Haddad, C.F.B. (2015). Genetic divergence in tropical anurans: Deeper phylogeographic structure in forest specialists and in topographically complex regions. Evolutionary Biology. 29: 765-785. [ Links ]

Salzburger, W., Ewing, G.B., von Haeseler, A. (2011). The performance of phylogenetic algorithms in estimating haplotype genealogies with migration. Molecular Ecology. 20: 1952-1963. Doi: 10.1111/j.1365-294X.2011.05066.x [ Links ]

Simões, P.I., Stow, A., Hõdl, W., Amézquita, A., Farias, I.P., Lima, A.P. (2014). The value of including intraspecific measures of biodiversity in environmental impact surveys is highlighted by the Amazonian brilliant-thighed frog (Allobates femoralis). Tropical Conservation Science. 7: 811-828. Doi: 10.1177/194008291400700416 [ Links ]

Tamura, K., Stecher, G., Peterson, D., Filipski, A., Kumar, S. (2013). MEGA6: Molecular evolutionary genetics analysis version 6.0. Molecular Biology and Evolution. 30: 2725-2729. Doi: 10.1093/molbev/mst197 [ Links ]

Thompson, J.D., Higgins, D.G., Gibson, T.J. (1996). CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Research. 22: 4673-4680. Doi: 10.1093/nar/22.22.4673 [ Links ]

Vieites, D.R., Wollenberg, K.C., Andreone, F., Kohler, J., Glaw, F., Vences, M. (2009) Vast underestimation of Madagascar's biodiversity evidenced by an integrative amphibian inventory. Proceedings of the National Academy of Sciences. 106: 8267-8272. Doi: 10.1073/pnas.0810821106 [ Links ]

Vos, C.C., Antonisse-de Jong, A.G., Goedhart, P.W., Smulders, M.J.M. (2001). Genetic similarity as a measure for connectivity between fragmented populations of the moor frog (Rana arvalis). Heredity. 86: 598-608. PMID: 11554976 [ Links ]

Wallace, A.R. (1852). On the Monkeys of the Amazon. Zoological Society of London. 20: 107-110. [ Links ]

Wells, K.D. (2007). The Ecology and Behavior of Amphibians. Chicago, USA: The University of Chicago Press. Pp. 1148 [ Links ]

Wright, S. (1943). Isolation by distance. Genetics. 28: 114-138. Doi: 10.5194/isprs-Archives-XLII-5-W1-419-2017 [ Links ]

Citation: Santana-Cornélio G, Araújo-de -Oliveira E, Magalhães-Xavier K, et al. The genetic structure of Pristimantis íatro (Anura: Craugastoridae) mirrors traits of their life history. Rev. Acad. Colomb. Cienc. Ex. Fis. Nat. 44(172): 729-739, julio-septiembre de 2020. doi: https://doi.org/10.18257/raccefyn.956

Editor: Martha Patricia Ramírez Pinilla

Supplementary informationSupplementary material 1: Genetic distance p uncorrected (in %) COI between sampling locations of P. latro; 1) P. latro (Altamira); 2) P. latro (Brasil Novo); 3) P. latro (Medicilândia); 4) P. latro (Santarém); 5) P. latro (Vitória do Xingu); 6) P. latro (Anapu); 7) P. latro (Senador José Porfirio); 8) P. fenestratus (Borba); 9) P. fenestratus (Manaus); 10) P. latro (Marabá); 11) Pristimantis sp.; (Maranhão); 12) P. zeuctotylus (Monte Alegre); 13) P. bipunctatus (Peru) and 14) P. vilarsi (Colombia). See the Supplementary material 1 in https://www.raccefyn.co/index.php/racceiyn/article/view/956/2807. Supplementary material 2: Genetic distance p uncorrected (in %) 16S rRNA between sampling locations of P. latro 1) P. latro (Santarém, PA); 2) P. latro (Anapu, PA); 3) P. latro (Altamira, PA); 4) P. latro (Senador José Porfírio, PA); 5) P. latro (Altamira, PA); 6); P. latro (Medicilândia, PA); 7) P. latro (Brasil Novo, PA); 8) P. fenestratus (Manaus, AM); 9) P. psamaipatae (Bolivia); 10) P. koehleri (Bolívia); 11) P. chiastonotus (French Guyana) 12) P. fenestratus (Borba, AM); 13) P. fenestratus (Bolivia); 14) P. fenestratus (Borba, AM); 15) P. vilarsi (Colombia); 16) P. conspicillatus (Ecuador); 17) P. zeeuctotylus (Monte Alegre, PA); 18) P. achatinus (Colombia); 19) P. zeeuctotylus (Obidos, PA); 20) P. zeeuctotylus (Alenquer, PA) and 21) P. skydmainos (Peru). See the Supplementary material 2 in https://www.raccefyn.co/index.php/raccefyn/article/view/956/2808

Authors contribution GSC laboratory work and writing of the manuscript; EAO laboratory work, statistical analysis, writing and revision of the manuscript; KMX statistical analysis (AMOVA) and review of the manuscript; GWBS revision of the manuscript; IF manuscript writing; LRRR writing and revision of the manuscript; EJHR writing and revision of the manuscript.

Conflicts of interest None

Received: July 31, 2019; Accepted: March 27, 2020

*Corresponding autor: Emil José Hernández-Ruz; emilhjh@yahoo.com

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License