Temporal and spatial diversification along the Amazonia-Cerrado transition in Neotropical treefrogs of the Boana albopunctata species group

Despite extensive research on biodiversity in Neotropical forests, biodiversity in seasonally dry, open biomes in South America has been underestimated until recently. We leverage a widespread group, Boana albopunctata, to uncover cryptic lineages and investigate the timing of diversification in Neotropical anurans with a focus on dry diagonal biomes (Cerrado, Caatinga and Chaco) and the ecotone between Amazonia and the Cerrado. We inferred a multilocus phylogeny of the B. albopunctata species group that includes 15 of 18 described species, recovered two cryptic species, and reconstructed the timing of diversification among species distributed across multiple South American biomes. One new potential species (B. aff. steinbachi), sampled in the Amazonian state of Acre, clustered within the B. calcara-fasciata species complex and is close to B. steinbachi. A second putative new species (B. aff. multifasciata), sampled in the Amazonia-Cerrado ecotone, is closely related to B. multifasciata. Lastly, we place a recently identified Cerrado lineage (B. aff. albopuncata) into the B. albopunctata species group phylogeny for the first time. Our ancestral range reconstruction showed that species in the B. albopuctata group likely dispersed from Amazonia-Cerrado into the dry-diagonal and Atlantic Forest. Intraspecies demography showed, for both B. raniceps and B. albopunctata, signs of rapid expansion across the dry diagonal. Similarly, for one clade of B. multifasciata, our analyses support an invasion of the Cerrado from Amazonia, followed by a rapid expansion across the open diagonal biomes. Thus, our study recovers several recent divergences along the Amazonia-Cerrado ecotone in northern Brazil. Tectonic uplift and erosion in the late Miocene and climate oscillations in the Pleistocene corresponded with estimated divergence times in the dry diagonal and Amazonia-Cerrado ecotone. Our study highlights the importance of these threatened open formations in the generation of biodiversity in the Neotropics.

Ecological divergence and synchronous Pleistocene diversification in the widespread South American butter frog complex

Phylogeographic studies primarily focus on the major role of landscape topography in driving lineage diversification. However, populational phylogeographic breaks may also occur as a result of either niche conservatism or divergence, in the absence of geographic barriers to gene flow. Furthermore, these two factors are not mutually exclusive and can act in concert, making it challenging to evaluate their relative importance on explaining genetic variation in nature. Herein, we use sequences of two mitochondrial and four nuclear genes to investigate the timing and diversification patterns of species pertaining to the Leptodactylus latrans complex, which harbors four morphologically cryptic species with broad distributions across environmental gradients in eastern South America. The origin of this species complex dates back to the late Miocene (ca. 5.5 Mya), but most diversification events occurred synchronically during the late Pleistocene likely as the result of ecological divergence driven by Quaternary climatic oscillations. Further, significant patterns of environmental niche divergences among species in the L. latrans complex imply that ecological isolation is the primary mode of genetic diversification, mostly because phylogenetic breaks are associated with environmental transitions rather than topographic barriers at both species and populational scale. We provided new insights about diversification patterns and processes within a species complex of broadly and continuously distributed group of frogs along South America.

Idiosyncratic responses to drivers of genetic differentiation in the complex landscapes of Isthmian Central America

Isthmian Central America (ICA) is one of the most biodiverse regions in the world, hosting an exceptionally high number of species per unit area. ICA was formed <25 million years ago and, consequently, its biotic assemblage is relatively young and derived from both colonization and in situ diversification. Despite intensive taxonomic work on the local fauna, the potential forces driving genetic divergences and ultimately speciation in ICA remain poorly studied. Here, we used a landscape genetics approach to test whether isolation by distance, topography, habitat suitability, or environment drive the genetic diversity of the regional frog assemblage. To this end, we combined data on landscape features and mitochondrial DNA sequence variation for nine codistributed amphibian species with disparate life histories. In five species, we found that at least one of the factors tested explained patterns of genetic divergence. However, rather than finding a general pattern, our results revealed idiosyncratic responses to historical and ecological processes, indicating that intrinsic life-history characteristics may determine the effect of different drivers of isolation on genetic divergence in ICA. Our work also suggests that the convergence of several factors promoting isolation among populations over a heterogeneous landscape might maximize genetic differentiation, despite short geographical distances. In conclusion, abiotic factors and geographical features have differentially affected the genetic diversity across the regional frog assemblage. Much more complex models (i.e., considering multiple drivers), beyond simple vicariance of Caribbean and Pacific lineages, are needed to better understand the evolutionary history of ICA's diverse biotas.

Variability in anuran advertisement call: a multi-level study with 15 species of monkey tree frogs (Anura, Phyllomedusidae)

Understanding the variability of acoustic signals is a first important step for the comprehension of the evolutionary processes that led to current diversity. Herein, we evaluate the variability of the advertisement call of the phyllomedusid species from the genera Phyllomedusa Wagler, 1830 and Pithecopus Cope, 1866 at different levels: intra-individual, intra-population, inter-population, intra-species, and inter-specific. An analysis of coefficients of variation showed a continuum of variability between the acoustic parameters analyzed, from static to highly dynamic. The majority of the variation was attributed to the inter-specific level, while call parameters at the intra-individual level varied the least. However, each parameter behaved differently with call interval being the most variable across all levels. Most temporal acoustic parameters were affected by environmental temperature, while pulse rate and dominant frequency were strongly influenced by body size. Only pulse rate was correlated to the geographic distance between populations, while all parameters presented a significant phylogenetic signal. Based on these results, we discuss the possible importance of different evolutionary forces and the usage of vocalizations for taxonomic purposes.

Unraveling the species diversity and relationships in the Leptodactylus mystaceus complex (Anura: Leptodactylidae), with the description of three new Brazilian species

Members of the Leptodactylus mystaceus species complex are widely distributed in forests and open formations of South America east of the Andes. Species of the complex are morphologically similar or indistinguishable among each other, but acoustic data have been the cornerstone for species discrimination across their geographic ranges. In this paper, we re-examine the monophyly, species diversity, and relationships in the L. mystaceus complex on the basis of morphology, coloration, acoustics, and DNA sequences. Morphological and color patterns originally used to the allocation of species to the L. mystaceus complex are also reassessed. Our results revealed three new species from the Brazilian Atlantic Forest and Cerrado, which are named and described herein, based mainly on acoustic and molecular data. Populations assigned to the lineage widely distributed across the South American Dry Diagonal (DD), reported in this study as L. cf. mystaceus, is likely paraphyletic with respect to the nominal species (Amazonian lineage), but additional data are still needed to address the taxonomic status of the DD lineage.

Dwarf geckos and giant rivers: the role of the São Francisco River in the evolution of Lygodactylus klugei (Squamata: Gekkonidae) in the semi-arid Caatinga of north-eastern Brazil

Species diversification can be strongly influenced by geomorphological features, such as mountains, valleys and rivers. Rivers can act as hard or soft barriers to gene flow depending on their size, speed of flow, historical dynamics and regional topographical characteristics. The São Francisco River (SFR) is the largest perennial river in the Caatinga biome in north-eastern Brazil and has been considered a barrier to gene flow and dispersal. Herein, we evaluated the role of the SFR on the evolution of Lygodactylus klugei, a small gecko from the Caatinga. Using a single-locus species delimitation method (generalized mixed Yule coalescent), we defined lineages (haploclades). Subsequently, we evaluated the role of the SFR in structuring genetic diversity in this species using a multilocus approach to quantify migration across margins. We also evaluated genetic structure based on nuclear markers, testing the number of populations found through an assignment test (STRUCTURE) across the species distribution. We recovered two mitochondrial lineages structured with respect to the SFR, but only a single population was inferred from nuclear markers. Given that we detected an influence of the SFR only on mitochondrial markers, we suggest that the current river course has acted as a relatively recent geographical barrier for L. klugei, for ~450 000 years.

The evolutionary history of Lygodactylus lizards in the South American open diagonal

The Pleistocenic Arc Hypothesis (PAH) posits that South American Seasonally Dry Tropical Forests (SDTF) were interconnected during Pleistocene glacial periods, enabling the expansion of species ranges that were subsequently fragmented in interglacial periods, promoting speciation. The lizard genus Lygodactylus occurs in Africa, Madagascar, and South America. Compared to the high diversity of African Lygodactylus, only two species are known to occur in South America, L. klugei and L. wetzeli, distributed in SDTFs and the Chaco, respectively. We use a phylogenetic approach based on mitochondrial (ND2) and nuclear (RAG-1) markers covering the known range of South American Lygodactylus to investigate (i) if they are monophyletic relative to their African congeners, (ii) if their divergence is congruent with the fragmentation of the PAH, and (iii) if cryptic diversity exists within currently recognized species. Maximum likelihood and Bayesian phylogenetic analyses recovered a well-supported monophyletic South American Lygodactylus, presumably resulting from a single trans-Atlantic dispersal event 29 Mya. Species delimitation analyses supported the existence of five putative species, three of them undescribed. Divergence times among L. klugei and the three putative undescribed species, all endemic to the SDTFs, are not congruent with the fragmentation of the PAH. However, fragmentation of the once broader and continuous SDTFs likely influenced the divergence of L. wetzeli in the Chaco and Lygodactylus sp. 3 (in a SDTF enclave in the Cerrado).

Speciation with gene flow in whiptail lizards from a Neotropical xeric biome

Two main hypotheses have been proposed to explain the diversification of the Caatinga biota. The riverine barrier hypothesis (RBH) claims that the São Francisco River (SFR) is a major biogeographic barrier to gene flow. The Pleistocene climatic fluctuation hypothesis (PCH) states that gene flow, geographic genetic structure and demographic signatures on endemic Caatinga taxa were influenced by Quaternary climate fluctuation cycles. Herein, we analyse genetic diversity and structure, phylogeographic history, and diversification of a widespread Caatinga lizard (Cnemidophorus ocellifer) based on large geographical sampling for multiple loci to test the predictions derived from the RBH and PCH. We inferred two well-delimited lineages (Northeast and Southwest) that have diverged along the Cerrado-Caatinga border during the Mid-Late Miocene (6-14 Ma) despite the presence of gene flow. We reject both major hypotheses proposed to explain diversification in the Caatinga. Surprisingly, our results revealed a striking complex diversification pattern where the Northeast lineage originated as a founder effect from a few individuals located along the edge of the Southwest lineage that eventually expanded throughout the Caatinga. The Southwest lineage is more diverse, older and associated with the Cerrado-Caatinga boundaries. Finally, we suggest that C. ocellifer from the Caatinga is composed of two distinct species. Our data support speciation in the presence of gene flow and highlight the role of environmental gradients in the diversification process.

Spermatozoa of Pseudinae (Amphibia, Anura, Hylidae), with a test of the hypothesis that sperm ultrastructure correlates with reproductive modes in anurans

We describe, for the first time, the sperm ultrastructure of the two genera of Pseudinae. Based on sperm ultrastructure, the five species herein examined can be separated into three groups: one containing Pseudis paradoxa, P. bolbodactyla, and P. tocantins, the second containing P. minuta, and the third containing Lysapsus laevis. The midpiece is similar in all species and auxiliary fibers and the undulating membrane are absent. In Pseudis a subacrosomal cone and a multilaminar structure (P. minuta) or a granular material (P. paradoxa group) are seen above the nucleus. Lysapsus laevis has only remnants of the subacrosomal cone. All species have peripheral fibers associated with the outer doublets of the axoneme. We tested the hypothesis of correlation between the presence of an undulating membrane and fertilization environments in anurans using a concentrated changes test (CCT) based on the Hay et al. (Mol Biol Evol 1995;12:928-937) hypothesis of phylogenetic relationships among anuran families. Only a subset of the resolved topologies derived from the Hay et al. (1995) cladogram, where Ranoidea is the sister-group of Sooglossidae, produced significant probabilities of the CCT. Therefore, support for the correlation between sperm ultrastructure and fertilization environments in anurans is, at best, equivocal.

Biflagellate spermatozoon of the poison-dart frogs Epipedobates femoralis and Colostethus sp. (Anura, Dendrobatidae)

This study describes the spermatozoa of the dendrobatids Epipedobates femoralis and Colostethus sp. using light and transmission electron microscopy. Both species possess a biflagellate spermatozoon, an unusual characteristic only previously reported in two anuran species belonging to the families Leptodactylidae and Racophoridae. The acrosomal complex of both species consists of a conical acrosomal vesicle and a subacrosomal cone, both of which cover the anterior portion of the nucleus, but to differing extents. In the midpiece, the centrioles are disposed parallel to each other and to the cell axis and give rise to two axonemes. Two paraxonemal rods were also seen entering the nuclear fossa. Both flagella are surrounded by a single mitochondrial collar. Each flagellum is formed by an axial fiber connected to the axoneme by an axial sheath; juxta-axonemal fibers are absent. Our data seem to support that Epipedobates femoralis should be placed in a separate clade possibly related to Colosthetus and that these two genera may not be monophyletic.

The ultrastructure of the spermatozoa of Epipedobates flavopictus (Amphibia, Anura, Dendrobatidae), with comments on its evolutionary significance

We describe, for the first time, the spermatozoon ultrastructure of a dendrobatid frog, Epipedobates flavopictus. Mature spermatozoa of E. flavopictus are filiform, with a moderately curved head and a proportionally short tail. The acrosomal vesicle is a conical structure that covers the nucleus for a considerable distance. A homogeneous subacrosomal cone lies between the acrosome vesicle and the nucleus. The nucleus contains a nuclear space at its anterior end, and electron-lucent spaces and inclusions. No perforatorium is present. In the midpiece, the proximal centriole is housed inside a deep nuclear fossa. Mitochondria are scattered around the posterior end of the nucleus and inside the undulating membrane in the anterior portion of the tail. In transverse section the tail is formed by an U-shaped axial fiber connected to the axoneme through an axial sheath, which supports the undulating membrane. The juxta-axonemal fiber is absent. The spermatozoon of E. flavopictus has several characteristics not observed before in any anurans, such as a curved axial fiber, absence of a juxta-axonemal fiber, and presence of mitochondria in the typical undulating membrane. Our results endorse the view that, in anurans, the conical perforatorium and subacrosomal cone are homologous and that Dendrobatidae should be grouped within Bufonoidea rather than Ranoidea.