Search

Rivers are prominent landscape features, acting as key promoters of diversification among freshwater organisms. Albeit generally considered potential barriers to species movement, they may also facilitate gene flow and structure populations of semiaquatic species (Riverine Thruway Hypothesis, RTH). We evaluated the role of rivers on the processes responsible for current genetic variation in the semiaquatic frog Pseudis bolbodactyla, testing whether each hydrographic basin harbours distinct genetic lineages. We sequenced three markers on 166 samples from 13 localities along the Paraná (PR), Araguaia–Tocantins (AT), and São Francisco (SF) River basins in Brazil. We recovered three populations geographically matching each hydrographic basin. Our results indicate migration among basins, with the best model selected using approximate Bayesian computation, including migration between AT and SF and ancient gene flow from PR to the AT–SF ancestor. Our findings are likely related to the orogenic events in Central Brazil dating to the Late Miocene (5 Mya), when hydrographic basins and the geomorphological features of the Brazilian Shield were formed. This suggests that P. bolbodactyla probably represents a species complex, with each lineage occurring in a distinct hydrographic basin, matching the predictions of the RTH. [ABSTRACT FROM AUTHOR]

Copyright of Molecular Ecology is the property of Wiley-Blackwell and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

It has been hypothesized that the Araucaria Forest in southern Brazil underwent expansions in the past, driven either by human groups or by climate fluctuations of the Holocene and Pleistocene. Fossil pollen records of the Paraná pine Araucaria angustifolia, a dominant tree in that forest, provide some insights into when those may have occurred. Still, the timing of those expansions has never been estimated. To infer past range shifts and shed light on their main drivers, we employed next‐generation DNA sequencing (ddRADseq), machine learning, and a comprehensive database of fossil pollen records into a study of historical demographic inference and paleo‐distribution modeling of the Paraná pine. We found that A. angustifolia comprises two populations expanding at different times: one in the Mantiqueira mountain chain, the other in the southern Brazilian plateau. The southern population began to expand during the Last Glacial Period ~ 70 kya, long before human arrival in South America. Still, genetic analyses support that humans later impacted this population, resulting in lower genetic diversity, higher inbreeding, and high levels of gene flow over large distances with a weak pattern of isolation‐by‐distance. It is possible this resulted from human influence on seed dispersal and germination on the southern Brazilian plateau. The Mantiqueira population, in contrast, expanded only recently (~ 3 kya). This timing coincides with Holocene climatic changes and human settlements established further south, although, to date, there is little archeological evidence of human impact in the Mantiqueira. In addition, multitemporal species distribution models built from a combination of present‐day and pollen records infer a range expansion of the Araucaria Forest during glacial times until the cold humid HS1 event (~ 16 kya), when the forest was most widespread, with no evidence of glacial refugia. The combination of genomic and spatial analyses suggests that both human and climatic controls played a role in the dynamics of the Araucaria Forest. [ABSTRACT FROM AUTHOR]

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. Key words: acoustic communication, evolution, Pithecopus, Phyllomedusa, phylogenetic signal, trait diversification, variation. La compréhension de la variabilité des signaux acoustiques est une première étape importante pour comprendre les processus évolutifs qui ont mené à la diversité actuelle. Nous évaluons la variabilité des chants d'appel des espèces de phyllomédusidés des genres Phyllomedusa Wagler, 1830 et Pithecopus Cope, 1866 aux différentes échelles suivantes : chez un même individu, au sein de la population, entre populations, au sein de l'espèce et entre espèces. Une analyse des coefficients de variation montre un continuum de variabilité entre les paramètres acoustiques analysés, qui vont de statiques à hautement dynamiques. La majeure partie de la variation est attribuée à l'échelle interspécifique, alors que les paramètres du chant pour un même individu présentent la plus faible variation. Tous les paramètres se comportent toutefois différemment, l'intervalle entre les chants étant le paramètre le plus variable à toutes les échelles. La température ambiante influence la plupart des paramètres acoustiques temporels, alors que la taille du corps influence fortement le rythme des impulsions et la fréquence dominante. Seul le rythme des impulsions est corrélé à la distance géographique entre populations, alors que tous les paramètres présentent un signal phylogénétique significatif. À la lumière de ces résultats, nous discutons de l'importance possible de différentes forces évolutives et de l'utilisation de vocalisations à des fins taxonomiques. [Traduit par la Rédaction] Mots-clés : communication acoustique, évolution, Pithecopus, Phyllomedusa, signal phylonétique, diversification des caractères, variation., Introduction Understanding variability is fundamental for the comprehension of evolution (Hallgrimsson and Hall 2011). Darwin's observations on phenotypic variation were the basis for the development of the concept of natural [...]

Geographic barriers are frequently invoked to explain genetic structuring across the landscape. However, inferences on the spatial and temporal origins of population variation have been largely limited to evolutionary neutral models, ignoring the potential role of natural selection and intrinsic genomic processes known as genomic architecture in producing heterogeneity in differentiation across the genome. To test how variation in genomic characteristics (e.g. recombination rate) impacts our ability to reconstruct general patterns of differentiation between species that cooccur across geographic barriers, we sequenced the whole genomes of multiple bird populations that are distributed across rivers in southeastern Amazonia. We found that phylogenetic relationships within species and demographic parameters varied across the genome in predictable ways. Genetic diversity was positively associated with recombination rate and negatively associated with species tree support. Gene flow was less pervasive in genomic regions of low recombination, making these windows more likely to retain patterns of population structuring that matched the species tree. We further found that approximately a third of the genome showed evidence of selective sweeps and linked selection, skewing genome-wide estimates of effective population sizes and gene flow between populations toward lower values. In sum, we showed that the effects of intrinsic genomic characteristics and selection can be disentangled from neutral processes to elucidate spatial patterns of population differentiation. [ABSTRACT FROM AUTHOR]

Many phylogeographic studies on species with large ranges have found genetic–geographic structure associated with changes in habitat and physical barriers preventing or reducing gene flow. These interactions with geographic space, contemporary and historical climate, and biogeographic barriers have complex effects on contemporary population genetic structure and processes of speciation. While allopatric speciation at biogeographic barriers is considered the primary mechanism for generating species, more recently it has been shown that parapatric modes of divergence may be equally or even more common. With genomic data and better modeling capabilities, we can more clearly define causes of speciation in relation to biogeography and migration between lineages, the location of hybrid zones with respect to the ecology of parental lineages, and differential introgression of genes between taxa. Here, we examine the origins of three Nearctic milksnakes (Lampropeltis elapsoides , Lampropeltis triangulum and Lampropeltis gentilis) using genome-scale data to better understand species diversification. Results from artificial neural networks show that a mix of a strong biogeographic barrier, environmental changes, and physical space has affected genetic structure in these taxa. These results underscore conspicuous environmental changes that occur as the sister taxa L. triangulum and L. gentilis diverged near the Great Plains into the forested regions of the Eastern Nearctic. This area has been recognized as a region for turnover for many vertebrate species, but as we show here the contemporary boundary does not isolate these sister species. These two species likely formed in the mid-Pleistocene and have remained partially reproductively isolated over much of this time, showing differential introgression of loci. We also demonstrate that when L. triangulum and L. gentilis are each in contact with the much older L. elapsoides , some limited gene flow has occurred. Given the strong agreement between nuclear and mtDNA genomes, along with estimates of ecological niche, we suggest that all three lineages should continue to be recognized as unique species. Furthermore, this work emphasizes the importance of considering complex modes of divergence and differential allelic introgression over a complex landscape when testing mechanisms of speciation. [Cline; delimitation; Eastern Nearctic; Great Plains; hybrids; introgression; speciation.] [ABSTRACT FROM AUTHOR]

Aim: A long history of isolation coupled with complex landscapes makes Madagascar ideal for exploring the historical factors that have shaped patterns of population diversity and endemism. Many species‐level studies have suggested Late Quaternary climate change may have influenced population dynamics in the tropics, but unique biomes and individual species properties may have driven idiosyncratic responses to these shifts. Here, we use community‐scale population genetic data to investigate the impact of Quaternary climate cycles on shared demographic response and investigate the contributions of both biotic and abiotic factors which shape these demographic trajectories. Location: Madagascar. Taxon Reptiles and Amphibians. Methods: Using community‐scale population genetic data, we implement a hierarchical approximate Bayesian computation approach to evaluate the degree of synchronous population expansion during glacial cycles across herpetofaunal assemblages both within and across discrete biomes and taxonomic groups. We use Bayesian model averaging to identify intrinsic and extrinsic conditions predictive of individualistic demographic change. Results: We find that demographic responses are not uniform across groups, with more than 50% of all populations showing signal of recent expansion. Our explanatory models indicate species occupying narrow elevational ranges had a higher probability of expansion, while amphibian assemblages showed higher genetic diversity and greater departures from population neutrality. Expansion events were largely asynchronous, with coexpansion found in less than half of all populations. Exceptionally, 69%–74% of all humid‐restricted populations coexpanded during the start of the Last Glacial Period at around 100 kya, supporting the hypothesis of a more extensive humid forest cover for Madagascar during this time. Main conclusions: We show that differences in life history and regional biogeography have contributed to patterns of richness and endemism found across Madagascar, and historical connectivity across humid forests has been an important factor in shaping present‐day diversity and endemism on Madagascar, particularly for amphibians. [ABSTRACT FROM AUTHOR]

Phylogeography investigates historical drivers of the geographical distribution of intraspecific lineages. Particular attention has been given to ecological, climatic and geological processes in the diversification of the Neotropical biota. Several species sampled across the South American diagonal of open formations (DOF), comprising the Caatinga, Cerrado and Chaco biomes, experienced range shifts coincident with Quaternary climatic changes. However, comparative studies across different spatial, temporal and biological scales on DOF species are still meagre. Here, we combine phylogeographical model selection and machine learning predictive frameworks to investigate the influence of Pleistocene climatic changes on several plant and animal species from the DOF. We assembled mitochondrial/chloroplastic DNA sequences in public repositories and inferred the demographic responses of 44 species, comprising 70 intraspecific lineages of plants, lizards, frogs, spiders and insects. We then built a random forest model using biotic and abiotic information to identify the best predictors of demographic responses in the Pleistocene. Finally, we assessed the temporal synchrony of species demographic responses with hierarchical approximate Bayesian computation. Biotic variables related to population connectivity, gene flow and habitat preferences largely predicted how species responded to Pleistocene climatic changes, and demographic changes were synchronous primarily during the Middle Pleistocene. Although 22 (~31%) lineages underwent demographic expansion, presumably associated with the spread of aridity during the glacial Pleistocene periods, our findings suggest that nine lineages (~13%) exhibited the opposite response due to taxon‐specific attributes. [ABSTRACT FROM AUTHOR]

The article discusses the corrected taxonomic history of the North American ratsnakes (Pantherophis obsoletus Complex). It mentions that differential isolation across the genomes of geographic lineages is possible, with some chromosomal regions showing low levels of introgression and others showing considerable admixture. It states that differential level of migration among loci is expected under a process of speciation with ecological adaptation via natural selection.

Aim: Phylogeographic studies show how historical and current changes in landscapes shape the geographic distribution of genetic diversity in species of animals and plants. In particular, for the species of the Diagonal of Open Formations (DOF), the compartmentalization of the Central Brazilian Plateau (CBP) during the Tertiary and climatic oscillations during the Quaternary have often been invoked to explain the origin and current patterns of biodiversity. We investigated how landscape changes and climatic oscillations shaped the distribution and diversification history of a widespread South American treefrog. Location: South American Diagonal of Open Formations (DOF) including Caatinga, Cerrado, and Chaco biomes. Taxon: Treefrog Boana raniceps. Methods: We used a multi‐locus dataset from 288 individual frogs collected at 115 localities throughout most of the species' distribution. We used population assignment analysis, species distribution models, historical demography models, approximate Bayesian computation and landscape genetic analyses to test alternative hypotheses of diversification. Results: We found two genetic lineages that diverged during the mid‐Pleistocene with continued gene flow. Approximate Bayesian computation supported a scenario of isolation with migration until the Last Glacial Maximum, followed by more recent population expansion in north‐eastern Brazil and stability at the southwest in South America. Isolation by environment was the best predictor of genetic distance between populations, which is in accordance with their different environmental niches. As Boana raniceps is a lowland species, steep slopes in the CBP likely restrained gene flow enough to sustain population divergence. We found evidence for major range contraction during the Last Glacial Maximum, raising the possibility of synergic action of climate change and the CBP compartmentalization in regulating migration. Main conclusions: Our findings highlight how landscape and climatic changes can shape the diversification of DOF biota. Past climatic fluctuations and environmental resistance due to topography acted in concert, forming a semipermeable barrier to gene flow, promoting intraspecific differentiation in a continentally distributed species. [ABSTRACT FROM AUTHOR]

Species are being lost at an unprecedented rate during the Anthropocene. Progress has been made in clarifying how species traits influence their propensity to go extinct, but the role historical demography plays in species loss or persistence is unclear. In eastern North America, five charismatic landbirds went extinct last century, and the causes of their extinctions have been heavily debated. Although these extinctions are most often attributed to post-colonial human activity, other factors such as declining ancestral populations prior to European colonization could have made these species particularly susceptible. We used population genomic data from these extinct birds and compared them with those from four codistributed extant species. We found extinct species harboured lower genetic diversity and effective population sizes than extant species, but both extinct and non-extinct birds had similar demographic histories of population expansion. These demographic patterns are consistent with population size changes associated with glacial–interglacial cycles. The lack of support for overall population declines during the Pleistocene corroborates the view that, although species that went extinct may have been vulnerable due to low diversity or small population size, their disappearance was driven by human activities in the Anthropocene. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Genetic structure can be influenced by local adaptation to environmental heterogeneity and biogeographic barriers, resulting in discrete population clusters. Geographic distance among populations, however, can result in continuous clines of genetic divergence that appear as structured populations. Here, we evaluate the relevant importance of these three factors over a landscape characterized by environmental heterogeneity and the presence of a hypothesized biogeographic barrier in producing population genetic structure within 13 codistributed snake species using a genomic data set. We demonstrate that geographic distance and environmental heterogeneity across western North America contribute to population genomic divergence. Surprisingly, landscape features long thought to contribute to biogeographic barriers play little role in divergence community wide. Our results suggest that isolation by environment is the most important contributor to genomic divergence. Furthermore, we show that models of population clustering that incorporate spatial information consistently outperform nonspatial models, demonstrating the importance of considering geographic distances in population clustering. We argue that environmental and geographic distances as drivers of community‐wide divergence should be explored before assuming the role of biogeographic barriers. see also the Perspective by Alencar and Quental [ABSTRACT FROM AUTHOR]

For many species, climate oscillations drove cycles of population contraction during cool glacial periods followed by expansion during interglacials. Some groups, however, show evidence of uniform and synchronous expansion, while others display differences in the timing and extent of demographic change. We compared demographic histories inferred from genetic data across marine turtle species to identify responses to postglacial warming shared across taxa and to examine drivers of past demographic change at the global scale. Using coalescent simulations and approximate Bayesian computation (ABC), we estimated demographic parameters, including the likelihood of past population expansion, from a mitochondrial data set encompassing 23 previously identified lineages from all seven marine turtle species. For lineages with a high posterior probability of expansion, we conducted a hierarchical ABC analysis to estimate the proportion of lineages expanding synchronously and the timing of synchronous expansion. We used Bayesian model averaging to identify variables associated with expansion and genetic diversity. Approximately 60% of extant marine turtle lineages showed evidence of expansion, with the rest mainly exhibiting patterns of genetic diversity most consistent with population stability. For lineages showing expansion, there was a strong signal of synchronous expansion after the Last Glacial Maximum. Expansion and genetic diversity were best explained by ocean basin and the degree of endemism for a given lineage. Geographic differences in sensitivity to climate change have implications for prioritizing conservation actions in marine turtles as well as for identifying areas of past demographic stability and potential resilience to future climate change for broadly distributed taxa. [ABSTRACT FROM AUTHOR]

Many studies propose that Quaternary climatic cycles contracted and/or expanded the ranges of species and biomes. Strong expansion-contraction dynamics of biomes presume concerted demographic changes of associated fauna. The analysis of temporal concordance of demographic changes can be used to test the influence of Quaternary climate on diversification processes. Hierarchical approximate Bayesian computation ( hABC) is a powerful and flexible approach that models genetic data from multiple species, and can be used to estimate the temporal concordance of demographic processes. Using available single-locus data, we can now perform large-scale analyses, both in terms of number of species and geographic scope. Here, we first compared the power of four alternative hABC models for a collection of single-locus data. We found that the model incorporating an a priori hypothesis about the timing of simultaneous demographic change had the best performance. Second, we applied the hABC models to a data set of seven squamate and four amphibian species occurring in the Seasonally Dry Tropical Forests (Caatinga) in northeastern Brazil, which, according to paleoclimatic evidence, experienced an increase in aridity during the Pleistocene. If this increase was important for the diversification of associated xeric-adapted species, simultaneous population expansions should be evident at the community level. We found a strong signal of synchronous population expansion in the Late Pleistocene, supporting the increase of the Caatinga during this time. This expansion likely enhanced the formation of communities adapted to high aridity and seasonality and caused regional extirpation of taxa adapted to wet forest. [ABSTRACT FROM AUTHOR]

The Atlantic Forest (AF) of Brazil has long been recognized as a biodiversity conservation hotspot. Despite decades of studies the species inventory of this biome continues to increase with the discovery of cryptic diversity and the description of new species. Different diversification mechanisms have been proposed to explain the diversity in the region, including models of forest dynamics, barriers to gene flow and dispersal. Also, sea level change is thought to have influenced coastal diversification and isolated populations on continental islands. However, the timing and mode of diversification of insular populations in the AF region were rarely investigated. Here, we analyze the phylogeography and species diversity of the small-sized direct-developing frog Ischnocnema parva . These frogs are independent from water bodies but dependent on forest cover and high humidity, and provide good models to understand forest dynamics and insular diversification. Our analysis was based on DNA sequences for one mitochondrial and four nuclear genes of 71 samples from 18 localities including two islands, São Sebastião, municipality of Ilhabela, and Mar Virado, municipality of Ubatuba, both in the state of São Paulo. We use molecular taxonomic methods to show that I. parva is composed of six independently evolving lineages, with the nominal I. parva likely endemic to the type locality. The time-calibrated species tree shows that these lineages have diverged in the Pliocene and Pleistocene, suggesting the persistence of micro-refuges of forest in the AF. For the two insular populations we used approximate Bayesian computation to test different diversification hypotheses. Our findings support isolation with migration for São Sebastião population, with ∼1 Mya divergence time, and isolation without migration for Mar Virado population, with ∼13 Kya divergence time, suggesting a combination of different processes for diversification on AF islands. [ABSTRACT FROM AUTHOR]

The South American sea lion (Otaria flavescens) is widely distributed along the southern Atlantic and Pacific coasts of South America with a history of significant commercial exploitation. We aimed to evaluate the population genetic structure and the evolutionary history of South American sea lion along its distribution by analyses of mitochondrial DNA (mtDNA) and 10 nuclear microsatellites loci. We analyzed 147 sequences of mtDNA control region and genotyped 111 individuals of South American sea lion for 10 microsatellite loci, representing six populations (Peru, Northern Chile, Southern Chile, Uruguay (Brazil), Argentina and Falkland (Malvinas) Islands) and covering the entire distribution of the species. The mtDNA phylogeny shows that haplotypes from the two oceans comprise two very divergent clades as observed in previous studies, suggesting a long period (>1 million years) of low inter-oceanic female gene flow. Bayesian analysis of bi-parental genetic diversity supports significant (but less pronounced than mitochondrial) genetic structure between Pacific and Atlantic populations, although also suggested some inter-oceanic gene flow mediated by males. Higher male migration rates were found in the intra-oceanic population comparisons, supporting very high female philopatry in the species. Demographic analyses showed that populations from both oceans went through a large population expansion ~10,000 years ago, suggesting a very similar influence of historical environmental factors, such as the last glacial cycle, on both regions. Our results support the proposition that the Pacific and Atlantic populations of the South American sea lion should be considered distinct evolutionarily significant units, with at least two managements units in each ocean. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Scincine lizards in Madagascar form an endemic clade of about 60 species exhibiting a variety of ecomorphological adaptations. Several subclades have adapted to burrowing and convergently regressed their limbs and eyes, resulting in a variety of partial and completely limbless morphologies among extant taxa. However, patterns of limb regression in these taxa have not been studied in detail. Here we fill this gap in knowledge by providing a phylogenetic analysis of DNA sequences of three mitochondrial and four nuclear gene fragments in an extended sampling of Malagasy skinks, and microtomographic analyses of osteology of various burrowing taxa adapted to sand substrate. Based on our data we propose to (i) consider Sirenoscincus Sakata & Hikida, 2003, as junior synonym of Voeltzkowia Boettger, 1893; (ii) resurrect the genus name Grandidierina Mocquard, 1894, for four species previously included in Voeltzkowia; and (iii) consider Androngo Brygoo, 1982, as junior synonym of Pygomeles Grandidier, 1867. By supporting the clade consisting of the limbless Voeltzkowia mira and the forelimb-only taxa V. mobydick and V. yamagishii, our data indicate that full regression of limbs and eyes occurred in parallel twice in the genus Voeltzkowia (as hitherto defined) that we consider as a sand-swimming ecomorph: in the Voeltzkowia clade sensu stricto the regression first affected the hindlimbs and subsequently the forelimbs, whereas the Grandidierina clade first regressed the forelimbs and subsequently the hindlimbs following the pattern prevalent in squamates. Timetree reconstructions for the Malagasy Scincidae contain a substantial amount of uncertainty due to the absence of suitable primary fossil calibrations. However, our preliminary reconstructions suggest rapid limb regression in Malagasy scincids with an estimated maximal duration of 6 MYr for a complete regression in Paracontias, and 4 and 8 MYr respectively for complete regression of forelimbs in Grandidierina and hindlimbs in Voeltzkowia. [ABSTRACT FROM AUTHOR]

Species distributed across vast continental areas and across major biomes provide unique model systems for studies of biotic diversification, yet also constitute daunting financial, logistic and political challenges for data collection across such regions. The tree frog Dendropsophus minutus (Anura: Hylidae) is a nominal species, continentally distributed in South America, that may represent a complex of multiple species, each with a more limited distribution. To understand the spatial pattern of molecular diversity throughout the range of this species complex, we obtained DNA sequence data from two mitochondrial genes, cytochrome oxidase I (COI) and the 16S rhibosomal gene (16S) for 407 samples of D. minutus and closely related species distributed across eleven countries, effectively comprising the entire range of the group. We performed phylogenetic and spatially explicit phylogeographic analyses to assess the genetic structure of lineages and infer ancestral areas. We found 43 statistically supported, deep mitochondrial lineages, several of which may represent currently unrecognized distinct species. One major clade, containing 25 divergent lineages, includes samples from the type locality of D. minutus. We defined that clade as the D. minutus complex. The remaining lineages together with the D. minutus complex constitute the D. minutus species group. Historical analyses support an Amazonian origin for the D. minutus species group with a subsequent dispersal to eastern Brazil where the D. minutus complex originated. According to our dataset, a total of eight mtDNA lineages have ranges >100,000 km2. One of them occupies an area of almost one million km2 encompassing multiple biomes. Our results, at a spatial scale and resolution unprecedented for a Neotropical vertebrate, confirm that widespread amphibian species occur in lowland South America, yet at the same time a large proportion of cryptic diversity still remains to be discovered. [ABSTRACT FROM AUTHOR]

Day geckos of the Phelsuma lineata group are widespread in Madagascar and have been historically split into numerous species and subspecies based almost exclusively on differences in coloration and body size. We apply phylogenetic and phylogeographic methods to examine the biogeography and taxonomy of these lizards, including explicit tests of various biogeographic predictions and based on a molecular data set covering much of the distribution ranges of all species and subspecies of P. lineata, P. dorsivittata, P. comorensis, P. hoeschi, P. kely, and P. pusilla in Madagascar (and the Comoros archipelago for P. comorensis). Sequences of the mitochondrial 16S rRNA and the nuclear RAG-1 gene fragment were determined from 376 samples, and a multigene mtDNA phylogeny of the species group was constructed for the main phylogroups identified in the 16S haplotype network. We used the 16S sequences to estimate the geographic location of the ancestor of each major mtDNA clade and to infer their demographic history using a variety of statistical tools. Our phylogeny separates the taxa analyzed into two well-supported major subclades mainly occurring in the north respectively east of the island. Mismatch distribution of samples together with rejection of neutrality, the results of Bayesian Skyline Plots analysis, and a star-like network suggests a recent demographic expansion for the P. l. lineata lineage into the eastern lowlands, while the highland ( P. l. elanthana) and northern clades ( P. dorsivittata and P. l. punctulata) show signatures of rather stable populations. A major genetic discontinuity observed coincided with a northern lowland stretch that separates mid-altitude rainforests in the north from those in the center and south. Our analysis points to numerous unsolved taxonomic problems in this group of geckos, especially in the small-sized taxa ( P. hoeschi, P. kely, P. pusilla), and provides a basis for a future comprehensive taxonomic revision, which will require integrative analysis of molecular, morphological and chromatic data as well as careful examination of type specimens. [ABSTRACT FROM AUTHOR]

Based on concordant differences in male advertisement call, tadpole morphology, and absence of haplotype sharing in mitochondrial and nuclear DNA markers, we describe a new species of treefrog from Ranomafana National Park in the southern central east of Madagascar. In its adult stage Boophis narinsi sp. n. is highly similar to its sister species, Boophis majori, but appears to differ in having longer hindlimbs. The genetic divergences between these two species (2.5-3.3% in a fragment of the 16S rRNA gene, depending on fragment length and individual haplotype analyzed) are below the threshold typically characterizing distinct species of anurans. Together with their relatively small and largely overlapping ranges and their sympatric occurrence in Ranomafana National Park, this indicates that they potentially could have originated rather recently by adaptive speciation under parapatric or sympatric conditions. Most studies on amphibian speciation have so far by default assumed vicariant speciation. We suggest that alternative speciation scenarios should be considered in future works and characterize settings in which more reliable assessments of adaptive parapatric or sympatric speciation could be carried out. [ABSTRACT FROM AUTHOR]

Graphical abstract Highlights • Extensive variation among gene trees jeopardizes the phylogenetic inference. • Different species tree inference methods yielded contradictory topologies. • The subspecies Neurergus strauchii barani is elevated to species level. Abstract We reconstruct the molecular phylogeny of Near Eastern mountain brook newts of the genus Neurergus (family Salamandridae) based on newly determined RADseq data, and compare the outcomes of concatenation-based phylogenetic reconstruction with species-tree inference. Furthermore, we test the current taxonomy of Neurergus (with four species: Neurergus strauchii, N. crocatus, N. kaiseri, and N. derjugini) against coalescent-based species-delimitation approaches of our genome-wide genetic data set. While the position of N. strauchii as sister species to all other Neurergus species was consistent in all of our analyses, the phylogenetic relationships between the three remaining species changed depending on the applied method. The concatenation approach, as well as quartet-based species-tree inference, supported a topology with N. kaiseri as the closest relative to N. derjugini , while full-coalescent species-tree inference approaches supported N. crocatus as sister species of N. derjugini. Investigating the individual signal of gene trees highlighted an extensive variation among gene histories, most likely resulting from incomplete lineage sorting. Coalescent-based species-delimitation models suggest that the current taxonomy might underestimate the species richness within Neurergus and supports seven species. Based on the current sampling, our analysis suggests that N. strauchii , N. derjugini and N. kaiseri might each be subdivided into further species. However, as amphibian species are known to be composed of deep conspecific lineages that do not always warrant species status, these results need to be cautiously interpreted in an integrative taxonomic framework. We hypothesize that the rather shallow divergences detected within N. kaiseri and N. derjugini likely reflect an ongoing speciation process and thus require further investigation. On the contrary, the much deeper genetic divergence found between the two morphologically and geographically differentiated subspecies of N. strauchii leads us to propose that N. s. barani should be considered a distinct species, Neurergus barani Öz, 1994. [ABSTRACT FROM AUTHOR]

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). [ABSTRACT FROM AUTHOR]

Highlights: [•] We present an integrative approach of phylogeny and populations genetics. [•] We analyzed 900 individuals of 37 populations for nuclear and mitochondrial markers. [•] We use a multilocus coalescent species tree approach to estimate species divergence. [•] We estimated a much more recent divergence time compared with previous studies. [•] We found differences in population history in the sister species of Salamandrina. [Copyright &y& Elsevier]

• Demographic history of bird species occurring in the Atlantic Forest sky islands. • Current and ancestral genetic diversity is associated with climate and elevation. • Supervised machine learning for testing the demographic history of multiple species. Montane organisms responded to Quaternary climate change by tracking suitable habitat along elevational gradients. However, it is unclear whether these past climatic dynamics generated predictable patterns of genetic diversity in co-occurring montane taxa. To test if the genetic variation is associated with historical changes in the elevational distribution of montane habitats, we integrated paleoclimatic data and a model selection approach for testing the demographic history of five co-distributed bird species occurring in the southern Atlantic Forest sky islands. We found that changes in historical population sizes and current genetic diversity are attributable to habitat dynamics among time periods and the current elevational distribution of populations. Taxa with populations restricted to the more climatically dynamic southern mountain block (SMB) had, on average, a six-fold demographic expansion, whereas the populations from the northern mountain block (NMB) remained constant. In the current configuration of the southern Atlantic Forest montane habitats, populations in the SMB have more widespread elevational distributions, occur at lower elevations, and harbor higher levels of genetic diversity than NMB populations. Despite the apparent coupling of demographic and climatic oscillations, our data rejected simultaneous population structuring due to historical habitat fragmentation. Demographic modeling indicated that the species had different modes of differentiation, and varied in the timing of divergence and the degree of gene flow across mountain blocks. Our results suggest that the heterogeneous distribution of genetic variation in birds of the Atlantic Forest sky islands is associated with the interplay between topography and climate of distinct mountains, leading to predictable patterns of genetic diversity. [ABSTRACT FROM AUTHOR]

Highlights: [•] Haplotype sharing in nuclear genes but no evidence for hybridization. [•] Simultaneous split of two species pairs possibly due to post-Messinian vicariance. [•] Steroidal toxins occur in high concentrations in all Salamandra species. [•] One alkaloid (samandarone) found also in other salamandrid genera. [Copyright &y& Elsevier]

[Display omitted] • Species in the Leptodactylus latrans complex showed old origin during late Miocene. • Species limits and genetic breaks coincides with ecologically transitional areas rather than barriers. • Intraspecific lineages exhibited synchronous divergence within late Pleistocene and population expansion since Last Interglacial. • Pleistocene climatic oscillations prompted ecologically adapted lineages along eastern South America. • Rivers might not be important in the diversification of Leptodactylus latrans 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 importances 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 scales. We provided new insights about diversification patterns and processes within a species complex of broadly and continuously distributed group of frogs along South America. [ABSTRACT FROM AUTHOR]

Highlights: [•] Phylogeography was reconstructed from >1000 samples of two widespread and large co-distributed European frog species. [•] Rana temporaria included five main mtDNA lineages of which four occur in the Iberian Peninsula. [•] Rana dalmatina is genetically uniform across its range except for a divergent mitochondrial lineage in southern Italy. [•] Rana dalmatina had consistently lower intrapopulational genetic variation (GV) than R. temporaria. [•] GV differences were found in mitochondrial and nuclear sequences, microsatellites and transcriptome SNPs. [Copyright &y& Elsevier]

Highlights: [•] A comprehensive analysis of genetic variation using mitochondrial and nuclear markers. [•] Three species could be reliable delimited but the taxonomic status of E. iberia, E. jaumei, and E. limbatus is questioned. [•] The group originated in eastern Cuba and only one of the main lineages expanded to the west in the Pleistocene. [•] Skin alkaloids were detected in all lineages of the group. [•] Alkaloid sequestration evolved first and was followed by miniaturization, diurnality and contrasting coloration. [ABSTRACT FROM AUTHOR]