Systematics and biogeography of the Boana albopunctata species group (Anura, Hylidae), with the description of two new species from Amazonia

Differential Mitochondrial Genome Expression of Four Hylid Frog Species under Low-Temperature Stress and Its Relationship with Amphibian Temperature Adaptation

Extreme weather poses huge challenges for animals that must adapt to wide variations in environmental temperature and, in many cases, it can lead to the local extirpation of populations or even the extinction of an entire species. Previous studies have found that one element of amphibian adaptation to environmental stress involves changes in mitochondrial gene expression at low temperatures. However, to date, comparative studies of gene expression in organisms living at extreme temperatures have focused mainly on nuclear genes. This study sequenced the complete mitochondrial genomes of five Asian hylid frog species: Dryophytes japonicus, D. immaculata, Hyla annectans, H. chinensis and H. zhaopingensis. It compared the phylogenetic relationships within the Hylidae family and explored the association between mitochondrial gene expression and evolutionary adaptations to cold stress. The present results showed that in D. immaculata, transcript levels of 12 out of 13 mitochondria genes were significantly reduced under cold exposure (p < 0.05); hence, we put forward the conjecture that D. immaculata adapts by entering a hibernation state at low temperature. In H. annectans, the transcripts of 10 genes (ND1, ND2, ND3, ND4, ND4L, ND5, ND6, COX1, COX2 and ATP8) were significantly reduced in response to cold exposure, and five mitochondrial genes in H. chinensis (ND1, ND2, ND3, ND4L and ATP6) also showed significantly reduced expression and transcript levels under cold conditions. By contrast, transcript levels of ND2 and ATP6 in H. zhaopingensis were significantly increased at low temperatures, possibly related to the narrow distribution of this species primarily at low latitudes. Indeed, H. zhaopingensis has little ability to adapt to low temperature (4 °C), or maybe to enter into hibernation, and it shows metabolic disorder in the cold. The present study demonstrates that the regulatory trend of mitochondrial gene expression in amphibians is correlated with their ability to adapt to variable climates in extreme environments. These results can predict which species are more likely to undergo extirpation or extinction with climate change and, thereby, provide new ideas for the study of species extinction in highly variable winter climates.

Distribution and diversification of Adelphobates, emblematic poison frogs from Brazilian Amazonia

Adelphobates contains three species, and the inaccurate identification of A. quinquevittatus and the scarcity of records of A. castaneoticus complicate inference of their distributions; the latter species occurs in sympatry with A. galactonotus. Our objective was to revise the distributions of Adelphobates by compiling data and modeling habitat suitability, as range limits may be shaped by landscape features and biotic interactions. We initially analyzed the existence of operational taxonomic units within the nominal species and subsequently inferred the observed and potential distributions, taking into account the possible independent lineages for the three species, and we also generated a molecular timetree to understand the chronology of interspecific diversification events. Adelphobates quinquevittatus was found to have a more easterly distribution than previously described, and specimens with phenotypic variation were found to occur in areas inconsistent with the modeling, and A. castaneoticus was concentrated in the Tapajós-Xingu interfluve, surrounded by A. galactonotus. Models indicated that the right bank of the Xingu River is suitable for both species, indeed, both were found there. Despite Adelphobates species having their distributions delimited by major Amazonian rivers, estimated divergence times predate the formation of the modern river network, suggesting that other mechanisms were involved in their diversification.

Marine introgressions and Andean uplift have driven diversification in neotropical Monkey tree frogs (Anura, Phyllomedusinae)

The species richness in the Neotropics has been linked to environmental heterogeneity and a complex geological history. We evaluated which biogeographic processes were associated with the diversification of Monkey tree frogs, an endemic clade from the Neotropics. We tested two competing hypotheses: the diversification of Phyllomedusinae occurred either in a "south-north" or a "north-south" direction in the Neotropics. We also hypothesized that marine introgressions and Andean uplift had a crucial role in promoting their diversification. We used 13 molecular markers in a Bayesian analysis to infer phylogenetic relationships among 57 species of Phyllomedusinae and to estimate their divergence times. We estimated ancestral ranges based on 12 biogeographic units considering the landscape modifications of the Neotropical region. We found that the Phyllomedusinae hypothetical ancestor range was probably widespread throughout South America, from Western Amazon to Southern Atlantic Forest, at 29.5 Mya. The Phyllomedusines' ancestor must have initially diverged through vicariance, generally followed by jump-dispersals and sympatric speciation. Dispersal among areas occurred mostly from Western Amazonia towards Northern Andes and the South American diagonal of dry landscapes, a divergent pattern from both "south-north" and "north-south" diversification hypotheses. Our results revealed a complex diversification process of Monkey tree frogs, occurring simultaneously with the orogeny of Northern Andes and the South American marine introgressions in the last 30 million years.

A new snouted treefrog (Anura, Hylidae, Scinax) from fluvial islands of the Juruena River, southern Brazilian Amazonia

Southern Amazonia is one of the less-explored regions by anuran taxonomists. We describe a small new species of snouted treefrog, genus Scinax, from this region, from a fluvial archipelago in the Juruena River, state of Mato Grosso, Brazil. The description is based on external morphology of adults and tadpoles, advertisement call and molecular data. The species is phylogenetically related to other snouted treefrogs of the Scinax cruentomma species group and shows the most southeastern distribution in Amazonia among its close relatives. It is distinguished from congeners mainly by its larger adult body size and bilobate vocal sac that reaches the level of the pectoral fold, a reddish-brown horizontal stripe on the iris, dark melanophores or blotches on the vocal sac and the throat of females, and the uniformly brown posterior portion of the thigh. The advertisement call comprises one pulsed note emitted at regular intervals, with a duration of 189-227 ms, 30-35 pulses/note and a dominant frequency of 2,250-2,344 Hz. The type locality is suffering several environmental impacts, including illegal mining, overfishing, unsustainable agriculture, uncontrolled logging and degradation associated with the construction of new hydroelectric dams. Further study of the biology and regional distribution of the new species is required to propose mitigation measures needed for its conservation.

Characterization of Two Mitogenomes of Hyla sanchiangensis (Anura: Hylidae), with Phylogenetic Relationships and Selection Pressure Analyses of Hylidae

Hyla sanchiangensis (Anura: Hylidae) is endemic to China and is distributed across Anhui, Zhejiang, Fujian, Guangdong, Guangxi, Hunan, and Guizhou provinces. The mitogenomes of H. sanchiangensis from two different sites (Jinxiu, Guangxi, and Wencheng, Zhejiang) were sequenced. Phylogenetic analyses were conducted, including 38 mitogenomes of Hylidae from the NCBI database, and assessed the phylogenetic relationship of H. sanchiangensis within the analyzed dataset. Two mitogenomes of H. sanchiangensis showed the typical mitochondrial gene arrangement with 13 protein-coding genes (PCGs), two ribosomal RNA genes (12S rRNA and 16S rRNA), 22 transfer RNA (tRNA) genes, and one non-coding control region (D-loop). The lengths of the 12S rRNA and 16S rRNA genes from both samples (Jinxiu and Wencheng) were 933 bp and 1604 bp, respectively. The genetic distance (p-distance transformed into percent) on the basis of the mitogenomes (excluding the control region) of the two samples was calculated as 4.4%. Hyla sanchiangensis showed a close phylogenetic relationship with the clade of (H. annectans + H. tsinlingensis), which was supported by ML and BI analyses. In the branch-site model, five positive selection sites were found in the clade of Hyla and Dryophytes: Cytb protein (at position 316), ND3 protein (at position 85), and ND5 protein (at position 400) have one site, respectively, and two sites in ND4 protein (at positions 47 and 200). Based on the results, we hypothesized that the positive selection of Hyla and Dryophytes was due to their experience of cold stress in historical events, but more evidence is needed to support this conclusion.

Performance of intron 7 of the β-fibrinogen gene for phylogenetic analysis: An example using gladiator frogs, Boana Gray, 1825 (Anura, Hylidae, Cophomantinae)

Boana, the third largest genus of Hylinae, has cryptic morphological species. The potential applicability of b-fibrinogen intron 7 – FGBI7 is explored to propose a robust phylogeny of Boana. The phylogenetic potential of FGBI7 was evaluated using maximum parsimony, MrBayes, and maximum likelihood analysis. Comparison of polymorphic sites and topologies obtained with concatenated analysis of FGBI7 and other nuclear genes (CXCR4, CXCR4, RHO, SIAH1, TYR, and 28S) allowed evaluation of the phylogenetic signal of FGBI7. Mean evolutionary rates were calculated using the sequences of the mitochondrial genes ND1 and CYTB available for Boana in GenBank. Dating of Boana and some of its groups was performed using the RelTime method with secondary calibration. FGBI7 analysis revealed high values at informative sites for parsimony. The absolute values of the mean evolutionary rate were higher for mitochondrial genes than for FGBI7. Dating of congruent Boana groups for ND1, CYTB, and FGBI7 revealed closer values between mitochondrial genes and slightly different values from those of FGBI7. Divergence times of basal groups tended to be overestimated when mtDNA was used and were more accurate when nDNA was used. Although there is evidence of phylogenetic potential arising from concatenation of specific genes, FGBI7 provides well-resolved independent gene trees. These results lead to a paradigm for linking data in phylogenomics that focuses on the uniqueness of species histories and ignores the multiplicities of individual gene histories.

Historical biogeography highlights the role of Miocene landscape changes on the diversification of a clade of Amazonian tree frogs

The diversification processes underlying why Amazonia hosts the most species-rich vertebrate fauna on earth remain poorly understood. We studied the spatio-temporal diversification of a tree frog clade distributed throughout Amazonia (Anura: Hylidae: Osteocephalus, Tepuihyla, and Dryaderces) and tested the hypothesis that Miocene mega wetlands located in western and central Amazonia impacted connectivity among major biogeographic areas during extensive periods. We assessed the group’s diversity through DNA-based (16S rRNA) species delimitation to identify Operational Taxonomic Units (OTUs) from 557 individuals. We then selected one terminal for each OTU (n = 50) and assembled a mitogenomic matrix (~14,100 bp; complete for 17 terminals) to reconstruct a Bayesian, time-calibrated phylogeny encompassing nearly all described species. Ancestral area reconstruction indicates that each genus was restricted to one of the major Amazonian biogeographic areas (western Amazonia, Guiana Shield and Brazilian Shield, respectively) between ~10 and 20 Mya, suggesting that they diverged and diversified in isolation during this period around the Pebas mega wetland. After 10 Mya and the transition to the modern configuration of the Amazon River watershed, most speciation within each genus continued to occur within each area. In Osteocephalus, only three species expanded widely across Amazonia (< 6 Mya), and all were pond-breeders. Species with other breeding modes remained mostly restricted to narrow ranges. The spectacular radiation of Osteocephalus was probably driven by climatic stability, habitat diversity and the acquisition of new reproductive modes along the Andean foothills and western Amazonia. Our findings add evidence to the importance of major hydrological changes during the Miocene on biotic diversification in Amazonia.

A new nurse frog of the Allobates tapajos species complex (Anura: Aromobatidae) from the upper Madeira River, Brazilian Amazonia

Cryptic diversity is extremely common in widespread Amazonian anurans, but especially in nurse frogs of the genus Allobates. There is an urgent need to formally describe the many distinct but unnamed species, both to enable studies of their basic biology but especially to facilitate conservation of threatened environments in which many are found. Here, we describe through integrative taxonomy a new species of the Allobates tapajos species complex from the upper Madeira River, southwestern Amazonia. Species delimitation analyses based on molecular data are congruent and delimit five candidate species in addition to A. tapajos sensu stricto. The new species is recovered as sister to A. tapajos clade F, a candidate species from Teles-Pires River, southeastern Amazonia. The new species differs from nominal congeners in adult and larval morphology and in male advertisement call. Egg deposition sites differ between east and west banks of the upper Madeira River, but there is no evidence of corresponding morphologic or bioacoustic differentiation. The new species appears to be restricted to riparian forests; its known geographic range falls entirely within the influence zone of reservoirs of two large dams, which underscores the urgent need of a conservation assessment through long-term monitoring. This region harbors the richest assemblage of Allobates reported for Brazilian Amazonia, with six nominal species and four additional candidate species awaiting formal description.

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.

Prepollex diversity and evolution in Cophomantini (Anura: Hylidae: Hylinae)

Several species of Cophomantini are known to have an enlarged prepollex, commonly modified as an osseous spine. We surveyed the osteology and myology of the prepollex and associated elements of 94 of the 190 species of Cophomantini, sampling all genera, except Nesorohyla. Two distinct prepollex morphologies were found: a blade-shaped and a spine-shaped morphology. We described the observed variation in 17 discrete characters to study their evolution in the most inclusive phylogenetic hypothesis for Cophomantini. Both morphologies evolved multiple times during the evolutionary history of this clade, but the origin of the spine-shaped distal prepollex in Boana and Bokermannohyla is ambiguous. The articulation of metacarpal II with the prepollex through a medial expansion of the metacarpal proximal epiphysis is a synapomorphy for Boana. The shape of the curve of the spine, and a large post-articular process of the distal prepollex, are synapomorphies for the Boana pulchella group, the latter being homoplastic in the Bokermannohyla martinsi group. Muscle character states associated with the spine-shaped prepollex are plesiomorphic for Cophomantini. We discuss evolution, function, behaviour and sexual dimorphism related to the prepollical elements. A bony spine is associated with fights between males, but forearm hypertrophy could be more related with habitat than with territorial combat.

Subtle environmental variation affects phenotypic differentiation of shallow divergent treefrog lineages in Amazonia

Amazonia harbours a vast biotic and ecological diversity, enabling investigation of the effects of microevolutionary processes and environmental variation on species diversification. Integrative approaches combining phenotypic and genetic variation can improve our knowledge on diversification processes in megadiverse regions. Here, we investigate the influence of environmental and geographic variation on the genetic and morphological differentiation in the Amazonian Boana calcarata-fasciata (Anura: Hylidae) species complex. We analysed the variation of one mtDNA gene from individuals of different forest environments, and assessed their phylogenetic relationships and species limits to define the lineages to perform a phenotypic-environmental approach. We collected morphological data (head shape and size) using 3D models and investigated the phylogenetic signal, evolutionary model and influence of environmental variables on morphology. We verified associations between environmental and geographical distances with morphological and genetic variation using distance-based redundancy analyses and Mantel tests. We found an even higher cryptic diversity than already recognized within the species complex. Body size and head shape varied among specimens, but did not present phylogenetic signal, diverging under a selective evolutionary model. Our results show that diverse factors have influenced morphological and genetic variation, but environmental conditions such as vegetation cover, precipitation and climate change velocity influenced morphological diversification. Possible population-level mechanisms such as parallel morphological evolution or plastic responses to similar environments could account for such patterns in these typical Amazonian treefrogs.