Pliocene and Pleistocene events shaping the genetic diversity within the central corridor of the Brazilian Atlantic Forest

Dinoponera lucida (Formicidae; Ponerinae) is an extinction‐threatened species of ant which is endemic in the central corridor of the Atlantic Forest. We used mitochondrial sequences of the Cox1, Cox2 and Cytb genes in order to infer some aspects of the evolutionary history and phylogeography of this ant. High genetic divergence and population structure were observed for the whole species. The current pattern of D. lucida diversity seems to be shaped during different geological times: middle Pliocene, early Pleistocene and mainly late Pleistocene, when the reduction of populations generated a structure pattern of the genetic variation of this species. Our data show that this structure results from the maintenance of populations of D. lucida within very small putative refuges to the south of the central Bahia refugium. We thus argue that, for some Atlantic forest endemic species, especially those resistant to very small fragments of forest, such as D. lucida, the small putative refuges were as important as, or even more important than, larger and stable refuges for the creation and maintenance of diversity, adding another piece to the puzzle of the mechanisms underlying local endemism. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 949–960.     

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.     

ABC inference of multi‐population divergence with admixture from unphased population genomic data

Rapidly developing sequencing technologies and declining costs have made it possible to collect genome‐scale data from population‐level samples in nonmodel systems. Inferential tools for historical demography given these data sets are, at present, underdeveloped. In particular, approximate Bayesian computation (ABC) has yet to be widely embraced by researchers generating these data. Here, we demonstrate the promise of ABC for analysis of the large data sets that are now attainable from nonmodel taxa through current genomic sequencing technologies. We develop and test an ABC framework for model selection and parameter estimation, given histories of three‐population divergence with admixture. We then explore different sampling regimes to illustrate how sampling more loci, longer loci or more individuals affects the quality of model selection and parameter estimation in this ABC framework. Our results show that inferences improved substantially with increases in the number and/or length of sequenced loci, while less benefit was gained by sampling large numbers of individuals. Optimal sampling strategies given our inferential models included at least 2000 loci, each approximately 2 kb in length, sampled from five diploid individuals per population, although specific strategies are model and question dependent. We tested our ABC approach through simulation‐based cross‐validations and illustrate its application using previously analysed data from the oak gall wasp, Biorhiza pallida.     

Contrasting Demographic Structure of Short‐ and Long‐lived Pioneer Tree Species on Amazonian Forest Edges

Although tropical forests have been rapidly converted into human‐modified landscapes, tree species response to forest edges remains poorly examined. In this study, we addressed four pioneer tree species to document demographic shifts experienced by this key ecological group and make inferences about pioneer response to forest edges. All individuals with dbh ≥ 1 cm of two short‐lived (Bellucia grossularioides and Cecropia sciadophylla) and two long‐lived species (Goupia glabra and Laetia procera) were sampled in 20 1‐ha forest edge plots and 20 1‐ha forest interior plots in Oiapoque and Manaus, Northeast and Central Amazon, respectively. As expected, pioneer stem density with dbh ≥ 1 cm increased by around 10–17‐fold along forest edges regardless of species, lifespan, and study site. Edge populations of long‐lived pioneers presented 84–94 percent of their individuals in sapling/subadult size classes, whereas edge populations of short‐lived pioneers showed 56–97 percent of their individuals in adult size classes. These demographic biases were associated with negative and positive net adult recruitment of long‐ and short‐lived pioneers, respectively. Our population‐level analyses support three general statements: (1) native pioneer tree species proliferate along forest edges (i.e., increased density), at least in terms of non‐reproductive individuals; (2) pioneer response to edge establishment is not homogeneous as species differ in terms of demographic structure and net adult recruitment; and (3) some pioneer species, particularly long‐lived ones, may experience population decline due to adult sensitivity to edge‐affected habitats.     

Phylogeography and diversification of an Amazonian understorey hummingbird: paraphyly and evidence for widespread cryptic speciation in the Plio‐Pleistocene

Straight‐billed Hermit Phaethornis bourcieri inhabits the understorey of upland terra firme forest throughout most of the Amazon basin. Currently, two allopatric taxa regarded as subspecies are recognized: P. b. bourcieri and P. b. major. However, the validity, interspecific limits and evolutionary history of these taxa are not yet fully elucidated. We use molecular characters to propose a phylogenetic hypothesis for populations and taxa grouped under Phaethornis bourcieri. Our results showed that P. bourcieri is part of the ‘Ametrornis’ clade, along with Phaethornis philippii and Phaethornis koepckeae, and that the subspecies major is more closely related to the latter two species than to populations grouped under nominate bourcieri. Our phylogenetic hypotheses recovered three main reciprocally monophyletic clades under nominate bourcieri separated by the lower Negro River and the Branco River or the Branco–Negro interfluve (clades B and C) and the upper Amazon (Solimões) or lower Marañon/Ucayali Rivers (clades C and D). Based on multi‐locus phylogeographic and population genetics approaches, we show that P. b. major is best treated as a separate species, and that P. b. bourcieri probably includes more than one evolutionary species, whose limits remain uncertain. The diversification of the ‘Ametrornis’ clade (P. bourcieri, P. philippii and P. koepckeae) is centred in the Amazon and appears to be closely linked to the formation of the modern Amazon drainage during the Plio‐Pleistocene.     

Multiple host shifts by the emerging honeybee parasite,Varroa jacobsoni

Host shifts are a key mechanism of parasite evolution and responsible for the emergence of many economically important pathogens. Varroa destructor has been a major factor in global honeybee (Apis mellifera) declines since shifting hosts from the Asian honeybee (Apis cerana) > 50 years ago. Until recently, only two haplotypes of V. destructor (Korea and Japan) had successfully host shifted to A. mellifera. In 2008, the sister species V. jacobsoni was found for the first time parasitizing A. mellifera in Papua New Guinea (PNG). This recent host shift presents a serious threat to world apiculture but also provides the opportunity to examine host shifting in this system. We used 12 microsatellites to compare genetic variation of V. jacobsoni on A. mellifera in PNG with mites on A. cerana in both PNG and surrounding regions. We identified two distinct lineages of V. jacobsoni reproducing on A. mellifera in PNG. Our analysis indicated independent host shift events have occurred through small numbers of mites shifting from local A. cerana populations. Additional lineages were found in the neighbouring Papua and Solomon Islands that had partially host shifted to A. mellifera, that is producing immature offspring on drone brood only. These mites were likely in transition to full colonization of A. mellifera. Significant population structure between mites on the different hosts suggested host shifted V. jacobsoni populations may not still reproduce on A. cerana, although limited gene flow may exist. Our studies provide further insight into parasite host shift evolution and help characterize this new Varroa mite threat to A. mellifera worldwide.     

Migration and isolation during the turbulent Ponto‐Caspian Pleistocene create high diversity in the crustacean Paramysis lacustris

The Ponto‐Caspian brackish‐water fauna inhabits estuaries and rivers of the Black, Azov and Caspian seas and is fragmented by higher salinity waters and a major interbasin watershed. The fauna is known for the high levels of endemism, complex zoogeographic histories, and as a recent source of successful invasive species. It remains debated whether the Black and Azov Sea brackish‐water populations survived unfavourable Pleistocene conditions in multiple separate refugia or whether the two seas were (repeatedly) recolonized from the Caspian. Using microsatellite and mtDNA markers, we demonstrate deep among‐ and within‐basin subdivisions in a widespread Ponto‐Caspian mysid crustacean Paramysis lacustris. Five genetic clusters were identified, but their relationships did not reflect the geography of the region. The Azov cluster was the most distinct (4–5% COI divergence), despite its geographic position in the corridor between Black and Caspian seas, and may represent a new species. In the northern Black Sea area, the Dnieper cluster was closer to the Caspian cluster than to the neighbouring Danube–Dniester–Bug populations, suggesting separate colonizations of the Black Sea. Overall, the data implied a predominant gene flow from the east to the Black Sea and highlight the importance of Caspian Sea transgressions in facilitating dispersal. Yet, the presence of distinct lineages in the Black Sea points to the persistence of isolated refugial populations that have gained diagnostic differences under presumably high mutation rates and large population sizes. The unfavourable Pleistocene periods in the Black Sea therefore appear to have promoted diversification of the brackish‐water lineages, rather than extirpated them.     

The historical bridge between the Amazon and the Atlantic Forest of Brazil: a study of molecular phylogeography with small mammals

Abstract Aim To examine how the genetic diversity of selected taxa of forest‐dwelling small mammals is distributed between and within the major rain forest domains of Amazonia and Atlantic Forest and the intervening interior forests of Brazil, as inferred by the relationships between gene genealogies and geography. I also addressed the historical importance of the central Brazilian forests in connecting Amazon and Atlantic Forest populations of rodents and marsupials. Methods I evaluated variation in the mitochondrial cytochrome b gene to estimate the levels of sequence divergence between those taxa occurring throughout the Amazon, Atlantic Forest, and forests in the Cerrado and Caatinga regions. I inferred the hierarchical relationships between haplotypes, populations and formal taxa using the cladistic approach of maximum parsimony. I compared areas and the clades identified by superimposing cladograms on the geographical distribution of samples. The degree of concordance both in phylogeny and the depth of the nodes in these phylogenies, in addition to patterns of geographical distribution of clades, permitted me to make inferences on how, when and where the taxa differentiated. Results Sequence similarity is often greater between samples from the Atlantic Forest and either Amazon or central Brazilian forests than it is within each of the two rain forest domains. The Atlantic Forest clades are either not reciprocally monophyletic or are the sister group to all the other clades. There is some indication of northern and southern components in the Atlantic Forest. Given the geographical distribution of clades and the relatively deep levels of divergence, the central Brazilian area does not behave as a separate region but is complementary to either Amazon or Atlantic Forest. Patterns of area relationships differ across taxa, suggesting that different processes and/or historic events affected the diversification within each lineage. Main conclusions The Amazon and the Atlantic forests are not exclusive in terms of their small mammal faunas; both overlap broadly with taxa occurring in gallery forests and dry forests in central Brazil. Central Brazilian forests are an integral part of the evolutionary scenario of lowland small mammals, playing an important role as present and past habitats for rain forest species. Therefore, representatives from this area should always be included in analyses of the evolutionary history of lowland rain forest faunas. The incongruence of branching patterns among areas is in agreement with recent results presented for Neotropical passerine birds and indicates that a single hypothesis of Neotropical area relationships is unlikely. These findings reinforce the idea that speciation in the Neotropics will not be explained by any single model of vicariance or climatic changes.     

Inferring the mode of origin of polyploid species from next‐generation sequence data

Many eukaryote organisms are polyploid. However, despite their importance, evolutionary inference of polyploid origins and modes of inheritance has been limited by a need for analyses of allele segregation at multiple loci using crosses. The increasing availability of sequence data for nonmodel species now allows the application of established approaches for the analysis of genomic data in polyploids. Here, we ask whether approximate Bayesian computation (ABC), applied to realistic traditional and next‐generation sequence data, allows correct inference of the evolutionary and demographic history of polyploids. Using simulations, we evaluate the robustness of evolutionary inference by ABC for tetraploid species as a function of the number of individuals and loci sampled, and the presence or absence of an outgroup. We find that ABC adequately retrieves the recent evolutionary history of polyploid species on the basis of both old and new sequencing technologies. The application of ABC to sequence data from diploid and polyploid species of the plant genus Capsella confirms its utility. Our analysis strongly supports an allopolyploid origin of C. bursa‐pastoris about 80 000 years ago. This conclusion runs contrary to previous findings based on the same data set but using an alternative approach and is in agreement with recent findings based on whole‐genome sequencing. Our results indicate that ABC is a promising and powerful method for revealing the evolution of polyploid species, without the need to attribute alleles to a homeologous chromosome pair. The approach can readily be extended to more complex scenarios involving higher ploidy levels.     

Evolutionarily significant units of the critically endangered leaf frog Pithecopus ayeaye (Anura,Phyllomedusidae) are not effectively preserved by the Brazilian protected areas network

Protected areas (PAs) are essential for biodiversity conservation, but their coverage is considered inefficient for the preservation of all species. Many species are subdivided into evolutionarily significant units (ESUs) and the effectiveness of PAs in protecting them needs to be investigated. We evaluated the usefulness of the Brazilian PAs network in protecting ESUs of the critically endangered Pithecopus ayeaye through ongoing climate change. This species occurs in a threatened mountaintop ecosystem known as campos rupestres. We used multilocus DNA sequences to delimit geographic clusters, which were further validated as ESUs with a coalescent approach. Ecological niche modeling was used to estimate spatial changes in ESUs’ potential distributions, and a gap analysis was carried out to evaluate the effectiveness of the Brazilian PAs network to protect P. ayeaye in the face of climate changes. We tested the niche overlap between ESUs to gain insights for potential management alternatives for the species. Pithecopus ayeaye contains at least three ESUs isolated in distinct mountain regions, and one of them is not protected by any PA. There are no climatic niche differences between the units, and only 4% of the suitable potential area of the species is protected in present and future projections. The current PAs are not effective in preserving the intraspecific diversity of P. ayeaye in its present and future range distributions. The genetic structure of P. ayeaye could represent a typical pattern in campos rupestres endemics, which should be considered for evaluating its conservation status.     

Unravelling the evolutionary history and future prospects of endemic species restricted to former glacial refugia

The contemporary distribution and genetic composition of biodiversity bear a signature of species’ evolutionary histories and the effects of past climatic oscillations. For many European species, the Mediterranean peninsulas of Iberia, Italy and the Balkans acted as glacial refugia and the source of range recolonization, and as a result, they contain disproportionately high levels of diversity. As these areas are particularly threatened by future climate change, it is important to understand how past climatic changes affected their biodiversity. We use an integrated approach, combining markers with different evolutionary rates and combining phylogenetic analysis with approximate Bayesian computation and species distribution modelling across temporal scales. We relate phylogeographic processes to patterns of genetic variation in Myotis escalerai, a bat species endemic to the Iberian Peninsula. We found a distinct population structure at the mitochondrial level with a strong geographic signature, indicating lineage divergence into separate glacial refugia within the Iberian refugium. However, microsatellite markers suggest higher levels of gene flow resulting in more limited structure at recent time frames. The evolutionary history of M. escalerai was shaped by the effects of climatic oscillations and changes in forest cover and composition, while its future is threatened by climatically induced range contractions and the role of ecological barriers due to competition interactions in restricting its distribution. This study warns that Mediterranean peninsulas, which provided refuge for European biodiversity during past glaciation events, may become a trap for limited dispersal and ecologically limited endemic species under future climate change, resulting in loss of entire lineages.     

Reanalysis of the biogeographical hypothesis of range expansion between robust and gracile capuchin monkeys

The data used by Lynch Alfaro et al. (2012a, Journal of Biogeography, 39 , 272–288) and the biogeographical hypothesis postulated by these authors to explain the current distribution of capuchin monkeys were reanalysed including additional cytochrome b data from Sapajus xanthosternos and Sapajus flavius. Our reconstructed phylogeny placed S. xanthosternos in a monophyletic clade representing the most basal lineage of this genus. All analyses indicated polyphyletic arrangements for several capuchin species, suggesting that incomplete lineage sorting has occurred during their evolution. These results also suggest that available molecular data lack adequate variation for accurately resolving species relationships. These results suggest that the divergence of capuchin monkey genera may have occurred in the Atlantic Forest. However, a more conclusive scenario and better resolution of the species tree requires correct identification of species, data from several unlinked nuclear loci from a higher number of individuals per species, and careful analysis of ancient DNA data from museum specimens.     

Directional selection in the evolution of elongated upper canines in clouded leopards and sabre‐toothed cats

Extremely developed or specialized traits such as the elongated upper canines of extinct sabre‐toothed cats are often not analogous to those of any extant species, which limits our understanding of their evolutionary cause. However, an extant species may have undergone directional selection for a similar extreme phenotype. Among living felids, the clouded leopard, Neofelis nebulosa, has exceptionally long upper canines for its body size. We hypothesized that directional selection generated the elongated upper canines of clouded leopards in a manner similar to the process in extinct sabre‐toothed cats. To test this, we developed an approach that compared the effect of directional selection among lineages in a phylogeny using a simulation of trait evolution and approximate Bayesian computation. This approach was applied to analyse the evolution of upper canine length in the Felidae phylogeny. Our analyses consistently showed directional selection favouring longer upper canines in the clouded leopard lineage and a lineage leading to the sabre‐toothed cat with the longest upper canines, Smilodon. Most of our analyses detected an effect of directional selection for longer upper canines in the lineage leading to another sabre‐toothed cat, Homotherium, although this selection may have occurred exclusively in the primitive species. In all the analyses, the clouded leopard and Smilodon lineages showed comparable directional selection. This implies that clouded leopards share a selection advantage with sabre‐toothed cats in having elongated upper canines.     

Colonization of the Mediterranean basin by the vector biting midge species Culicoides imicola: an old story

Understanding the demographic history and genetic make‐up of colonizing species is critical for inferring population sources and colonization routes. This is of main interest for designing accurate control measures in areas newly colonized by vector species of economically important pathogens. The biting midge Culicoides imicola is a major vector of orbiviruses to livestock. Historically, the distribution of this species was limited to the Afrotropical region. Entomological surveys first revealed the presence of C. imicola in the south of the Mediterranean basin by the 1970s. Following recurrent reports of massive bluetongue outbreaks since the 1990s, the presence of the species was confirmed in northern areas. In this study, we addressed the chronology and processes of C. imicola colonization in the Mediterranean basin. We characterized the genetic structure of its populations across Mediterranean and African regions using both mitochondrial and nuclear markers, and combined phylogeographical analyses with population genetics and approximate Bayesian computation. We found a west/east genetic differentiation between populations, occurring both within Africa and within the Mediterranean basin. We demonstrated that three of these groups had experienced demographic expansions in the Pleistocene, probably because of climate changes during this period. Finally, we showed that C. imicola could have colonized the Mediterranean basin in the Late Pleistocene or Early Holocene through a single event of introduction; however, we cannot exclude the hypothesis involving two routes of colonization. Thus, the recent bluetongue outbreaks are not linked to C. imicola colonization event, but rather to biological changes in the vector or the virus.     

Reconstructing the evolutionary history of an endangered subspecies across the changing landscape of the Great Central Valley of California

Identifying historic patterns of population genetic diversity and connectivity is a primary challenge in efforts to re‐establish the processes that have generated and maintained genetic variation across natural landscapes. The challenge of reconstructing pattern and process is even greater in highly altered landscapes where population extinctions and dramatic demographic fluctuations in remnant populations may have substantially altered, if not eliminated, historic patterns. Here, we seek to reconstruct historic patterns of diversity and connectivity in an endangered subspecies of woodrat that now occupies only 1–2 remnant locations within the highly altered landscape of the Great Central Valley of California. We examine patterns of diversity and connectivity using 14 microsatellite loci and sequence data from a mitochondrial locus and a nuclear intron. We reconstruct temporal change in habitat availability to establish several historical scenarios that could have led to contemporary patterns of diversity, and use an approximate Bayesian computation approach to test which of these scenarios is most consistent with our observed data. We find that the Central Valley populations harbour unique genetic variation coupled with a history of admixture between two well‐differentiated species of woodrats that are currently restricted to the woodlands flanking the Valley. Our simulations also show that certain commonly used analytical approaches may fail to recover a history of admixture when populations experience severe bottlenecks subsequent to hybridization. Overall our study shows the strength of combining empirical and simulation analyses to recover the history of populations occupying highly altered landscapes.