Interglacial microrefugia and diversification of a cactus species complex: phylogeography and palaeodistributional reconstructions forPilosocereus aurisetusand allies

Geoclimatic drivers of diversification in the largest arid and semi-arid environment of the Neotropics: Perspectives from phylogeography

The South American Dry Diagonal, also called the Diagonal of Open Formations, is a large region of seasonally dry vegetation extending from northeastern Brazil to northern Argentina, comprising the Caatinga, Cerrado, and Chaco subregions. A growing body of phylogeography literature has determined that a complex history of climatic changes coupled with more ancient geological events has produced a diverse and endemic-rich Dry Diagonal biota. However, the exact drivers are still under investigation, and their relative strengths and effects are controversial. Pleistocene climatic fluctuations structured lineages via vegetation shifts, refugium formation, and corridors between the Amazon and Atlantic forests. In some taxa, older geological events, such as the reconfiguration of the São Francisco River, uplift of the Central Brazilian Plateau, or the Miocene inundation of the Chaco by marine incursions, were more important. Here, we review the Dry Diagonal phylogeography literature, discussing each hypothesized driver of diversification and assessing degree of support. Few studies statistically test these hypotheses, with most support drawn from associating encountered phylogeographic patterns such as population structure with the timing of ancient geoclimatic events. Across statistical studies, most hypotheses are well supported, with the exception of the Pleistocene Arc Hypothesis. However, taxonomic and regional biases persist, such as a proportional overabundance of herpetofauna studies, and the under-representation of Chaco studies. Overall, both Pleistocene climate change and Neogene geological events shaped the evolution of the Dry Diagonal biota, though the precise effects are regionally and taxonomically varied. We encourage further use of model-based analyses to test evolutionary scenarios, as well as interdisciplinary collaborations to progress the field beyond its current focus on the traditional set of geoclimatic hypotheses.

Are 150 km of open sea enough? Gene flow and population differentiation in a bat-pollinated columnar cactus

Genetic differentiations and phylogeographical patterns are controlled by the interplay between spatial isolation and gene flow. To assess the extent of gene flow across an oceanic barrier, we explored the effect of the separation of the peninsula of Baja California on the evolution of mainland and peninsular populations of the long-lived columnar cactus Stenocereus thurberi. We analyzed twelve populations throughout the OPC distribution range to assess genetic diversity and structure using chloroplast DNA sequences. Genetic diversity was higher (Hd = 0.81), and genetic structure was lower (GST = 0.143) in mainland populations vs peninsular populations (Hd = 0.71, GST = 0.358 respectively). Genetic diversity was negatively associated with elevation but positively with rainfall. Two mainland and one peninsular ancestral haplotypes were reconstructed. Peninsular populations were as isolated among them as with mainland populations. Peninsular haplotypes formed a group with one mainland coastal population, and populations across the gulf shared common haplotypes giving support to regular gene flow across the Gulf. Gene flow is likely mediated by bats, the main pollinators and seed dispersers. Niche modeling suggests that during the Last Glacial Maximum (c. 130 ka), OPC populations shrank to southern locations. Currently, Stenocereus thurberi populations are expanding, and the species is under population divergence despite ongoing gene flow. Ancestral populations are located on the mainland and although vicariant peninsular populations cannot be ruled out, they are likely the result of gene flow across the seemingly formidable barrier of the Gulf of California. Still, unique haplotypes occur in the peninsula and the mainland, and peninsular populations are more structured than those on the mainland.

Phylogenomics and Biogeography of the Mammilloid Clade Revealed an Intricate Evolutionary History Arose in the Mexican Plateau

Mexico harbors ~45% of world’s cacti species richness. Their biogeography and phylogenomics were integrated to elucidate the evolutionary history of the genera Coryphantha, Escobaria, Mammillaria, Mammilloydia, Neolloydia, Ortegocactus, and Pelecyphora (Mammilloid Clade). We analyzed 52 orthologous loci from 142 complete genomes of chloroplast (103 taxa) to generate a cladogram and a chronogram; in the latter, the ancestral distribution was reconstructed with the Dispersal-Extinction-Cladogenesis model. The ancestor of these genera arose ~7 Mya on the Mexican Plateau, from which nine evolutionary lineages evolved. This region was the site of 52% of all the biogeographical processes. The lineages 2, 3 and 6 were responsible for the colonization of the arid southern territories. In the last 4 Mya, the Baja California Peninsula has been a region of prolific evolution, particularly for lineages 8 and 9. Dispersal was the most frequent process and vicariance had relevance in the isolation of cacti distributed in the south of Mexico. The 70 taxa sampled as Mammillaria were distributed in six distinct lineages; one of these presumably corresponded to this genus, which likely had its center of origin in the southern part of the Mexican Plateau. We recommend detailed studies to further determine the taxonomic circumscription of the seven genera.

Comparative phylogeography reveals the demographic patterns of neotropical ancient mountain species

Mountains are renowned for their bountiful biodiversity. Explanations on the origin of such abundant life are usually regarded to their orogenic history. However, ancient mountain systems with geological stability also exhibit astounding levels of number of species and endemism, as illustrated by the Brazilian Quartzitic Mountains (BQM) in Eastern South America. Thus, cycles of climatic changes over the last couple million years are usually assumed to play an important role in the origin of mountainous biota. These climatic oscillations potentially isolated and reconnected adjacent populations, a phenomenon known as flickering connectivity, accelerating speciation events due to range fragmentation, dispersion, secondary contact, and hybridization. To evaluate the role of the climatic fluctuations on the diversification of the BQM biota, we estimated the ancient demography of distinct endemic species of animals and plants using hierarchical approximate Bayesian computation analysis and Ecological Niche Modelling. Additionally, we evaluated if climatic oscillations have driven a genetic spatial congruence in the genetic structure of codistributed species from the Espinhaço Range, one of the main BQM areas. Our results show that the majority of plant lineages underwent a synchronous expansion over the Last Glacial Maximum (LGM, c. 21 thousand years ago), although we could not obtain a clear demographic pattern for the animal lineages. We also obtained a signal of a congruent phylogeographic break between lineages endemic to the Espinhaço Range, suggesting how ancient climatic oscillations might have driven the evolutionary history of the Espinhaço's biota.

Revisiting phylogeny, systematics, and biogeography of a Pleistocene radiation

PREMISE

Pilosocereus (Cactaceae) is an important dry forest element in all subregions and transitional zones of the neotropics, with the highest diversity in eastern Brazil. The genus is subdivided into informal taxonomic groups; however, most of these are not supported by recent molecular phylogenetic inferences. This lack of confidence is probably due to the use of an insufficient number of loci and the complexity of cactus diversification. Here, we explored the species relationships in Pilosocereus in more detail, integrating multilocus phylogenetic approaches with the assessment of the ancestral range and the effect of geography on diversification shifts.

METHODS

We used 28 nuclear, plastid, and mitochondrial loci from 54 plant samples of 31 Pilosocereus species for phylogenetic analyses. We used concatenated and coalescent phylogenetic trees and Bayesian models to estimate the most likely ancestral range and diversification shifts.

RESULTS

All Pilosocereus species were clustered in the same branch, except P. bohlei. The phylogenetic relationships were more associated with the geographic distribution than taxonomic affinities among taxa. The genus began diversifying during the Plio-Pleistocene transition in the Caatinga domain and experienced an increased diversification rate during the Calabrian age.

CONCLUSIONS

We recovered a well-supported multispecies coalescent phylogeny. Our results refine the pattern of rapid diversification of Pilosocereus species across neotropical drylands during the Pleistocene and highlight the need for taxonomic rearrangements in the genus. We recovered a pulse of diversification during the Pleistocene that was likely driven by multiple dispersal and vicariance events within and among the Caatinga, Cerrado, and Atlantic Forest domains.

The nuanced nature of mesic refugia in arid landscapes: a tale of two peas

BACKGROUND AND AIMS

Understanding how genetic diversity is distributed and maintained within species is a central tenet of evolutionary and conservation biology, yet is understudied in arid regions of the globe. In temperate, glaciated environments, high genetic diversity in plant species is frequently found in refugial areas, which are often associated with southern non-glaciated landscapes. In arid, unglaciated environments, landscape features providing mesic conditions are likely to be refugia, although our understanding needs more refinement in these biomes. We test whether refugia and nuclear diversity hotspots occur in high-elevation, topographically complex areas for co-distributed shrubs (Petalostylis labicheoides and Indigofera monophylla; Fabaceae) in the ancient, arid Pilbara bioregion of north-western Australia.

METHODS

We conducted extensive sampling of the Pilbara (>1400 individuals from 62 widespread populations) to detect patterns in nuclear diversity and structure based on 13-16 microsatellite loci. Evidence of historical refugia was investigated based on patterns of diversity in three non-coding chloroplast (cp) sequence regions for approx. 240 individuals per species. Haplotype relationships were defined with median-joining networks and maximum likelihood phylogenetic trees.

KEY RESULTS

We found cpDNA evidence for a high-elevation refugium in P. labicheoides but not for I. monophylla that instead exhibited extraordinary haplotype diversity and evidence for persistence across a widespread area. Nuclear diversity hotspots occurred in, but were not exclusive to, high-elevation locations and extended to adjacent, low-elevation riparian areas in both species.

CONCLUSIONS

Phylogeographic refugia in arid environments may occur in high-elevation areas for some species but not all, and may be influenced by species-specific traits: a mesic montane refugium in P. labicheoides could be related to its preference for growth in water-gaining areas, while a lack of such evidence in I. monophylla could be related to maintenance of cpDNA diversity in a large soil seed bank and dynamic evolutionary history. Mesic environments created by the intersection of topographically complex landscapes with riparian zones can be contemporary reservoirs of genetic diversity in arid landscapes.

Pleistocene climate and geomorphology drive the evolution and phylogeographic pattern of Triplophysa robusta (Kessler, 1876)

Montane systems provide excellent opportunities to study the rapid radiation influenced by geological and climatic processes. We assessed the role of Pleistocene climatic oscillations and mountain building on the evolution history of Triplophysa robusta, a cold-adapted species restricted to high elevations in China. We found seven differentiated sublineages of T. robusta, which were established during the Mid Pleistocene 0.87–0.61 Mya. The species distribution modeling (SDM) showed an expansion of T. robusta during the Last Glacial Maximum (LGM) and a considerable retraction during the Last Interglacial (LIG). The deep divergence between Clade I distributed in Qinling Mountains and Clade II in Northeastern Qinghai-Tibet Plateau (QTP) was mainly the result of a vicariance event caused by the rapid uplifting of Qinling Mountains during the Early Pleistocene. While the middling to high level of historical gene flow among different sublineages could be attributed to the dispersal events connected to the repetition of the glacial period during the Pleistocene. Our findings suggested that frequent range expansions and regressions due to Pleistocene glaciers likely have been crucial for driving the phylogeographic pattern of T. robusta. Finally, we urge a burning question in future conservation projection on the vulnerable cold-adapted species endemic to high elevations, as they would be negatively impacted by the recent rapid climate warming.

Evolutionary Genetics of Cacti: Research Biases, Advances and Prospects

Here, we present a review of the studies of evolutionary genetics (phylogenetics, population genetics, and phylogeography) using genetic data as well as genome scale assemblies in Cactaceae (Caryophyllales, Angiosperms), a major lineage of succulent plants with astonishing diversity on the American continent. To this end, we performed a literature survey (1992–2021) to obtain detailed information regarding key aspects of studies investigating cactus evolution. Specifically, we summarize the advances in the following aspects: molecular markers, species delimitation, phylogenetics, hybridization, biogeography, and genome assemblies. In brief, we observed substantial growth in the studies conducted with molecular markers in the past two decades. However, we found biases in taxonomic/geographic sampling and the use of traditional markers and statistical approaches. We discuss some methodological and social challenges for engaging the cactus community in genomic research. We also stressed the importance of integrative approaches, coalescent methods, and international collaboration to advance the understanding of cactus evolution.

Independent Evolutionary Lineages in a Globular Cactus Species Complex Reveals Hidden Diversity in a Central Chile Biodiversity Hotspot

Unraveling the processes involved in the origin of a substantial fraction of biodiversity can be a particularly difficult task in groups of similar, and often convergent, morphologies. The genus Eriosyce (Cactaceae) might present a greater specific diversity since much of its species richness might be hidden in morphological species complexes. The aim of this study was to investigate species delimitation using the molecular data of the globose cacti “E. curvispina”, which harbor several populations of unclear evolutionary relationships. We ran phylogenetic inferences on 87 taxa of Eriosyce, including nine E. curvispina populations, and by analyzing three plastid noncoding introns, one plastid and one nuclear gene. Additionally, we developed 12 new pairs of nuclear microsatellites to evaluate the population-level genetic structure. We identified four groups that originated in independent cladogenetic events occurring at different temporal depths; these groups presented high genetic diversity, and their populations were genetically structured. These results suggest a complex evolutionary history in the origin of globular cacti, with independent speciation events occurring at different time spans. This cryptic richness is underestimated in the Mediterranean flora of central Chile, and thus unique evolutionary diversity could be overlooked in conservation and management actions.

Coalescent-based species delimitation meets deep learning: Insights from a highly fragmented cactus system

Delimiting species boundaries is a major goal in evolutionary biology. An increasing volume of literature has focused on the challenges of investigating cryptic diversity within complex evolutionary scenarios of speciation, including gene flow and demographic fluctuations. New methods based on model selection, such as approximate Bayesian computation, approximate likelihoods, and machine learning are promising tools arising in this field. Here, we introduce a framework for species delimitation using the multispecies coalescent model coupled with a deep learning algorithm based on convolutional neural networks (CNNs). We compared this strategy with a similar ABC approach. We applied both methods to test species boundary hypotheses based on current and previous taxonomic delimitations as well as genetic data (sequences from 41 loci) in Pilosocereus aurisetus, a cactus species complex with a sky-island distribution and taxonomic uncertainty. To validate our method, we also applied the same strategy on data from widely accepted species from the genus Drosophila. The results show that our CNN approach has a high capacity to distinguish among the simulated species delimitation scenarios, with higher accuracy than ABC. For the cactus data set, a splitter hypothesis without gene flow showed the highest probability in both CNN and ABC approaches, a result agreeing with previous taxonomic classifications and in line with the sky-island distribution and low dispersal of P. aurisetus. Our results highlight the cryptic diversity within the P. aurisetus complex and show that CNNs are a promising approach for distinguishing complex evolutionary histories, even outperforming the accuracy of other model-based approaches such as ABC.

Comparative and predictive phylogeography in the South American diagonal of open formations: Unravelling the biological and environmental influences on multitaxon demography

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.

Pollination in the campo rupestre: a test of hypothesis for an ancient tropical mountain vegetation

The campo rupestre is a Neotropical OCBIL (old, climatically buffered infertile landscape), a grassy-shrub vegetation with high species richness and endemism, characterized by rocky outcrops surrounded by grasslands distributed in South American ancient mountaintops. We tested one OCBIL prediction: the prevalence of long-distance pollinators ensuring cross-pollination across the archipelago-like landscapes of the campo rupestre. We described the pollination systems and tested whether their frequency differed across vegetation types and elevation, focusing on long-distance systems. We performed non-systematic and systematic surveys of plants and plant-pollinator interactions across the elevation gradient and vegetation types. We also reviewed the literature on campo rupestre pollination and applied an accuracy criterion to infer 11 pollination systems. The bee system was split into large bee (long-distance) and small bee (shorter distances) to test the prevalence of long-distance pollination systems. We surveyed 413 pollinator species, mostly bees (220) and flies (69). Among the 636 plant species studied, the bee pollination system was dominant (56%), followed by wind and hummingbird. Wind, small-bee and fly pollination systems increased with elevation, and small-bee and wind pollination systems prevailed in grasslands. Large-bee and hummingbird long-distance pollination systems remained unchanged with elevation and were more frequent in the highly isolated rocky outcrops corroborating the OCBIL theory.

Quaternary climatic fluctuations influence the demographic history of two species of sky-island endemic amphibians in the Neotropics

We evaluated the role of Quaternary climatic fluctuations on the demographic history and population structure of amphibian species endemic to the 'campo rupestre' in the Neotropics, evaluating their distributional shifts, demographic changes, and lineage formation from the end of Pleistocene to present. We chose two anurans endemic to the high-elevation 'campo rupestre' in the Espinhaço Range (ER) in northeastern and southeastern Brazil (Bokermannohyla alvarengai and Bokermannohyla oxente), as models to test the role of Quaternary climatic fluctuations over their distribution range in this region. We collected tissue samples throughout their distribution range and used statistical phylogeography to examine processes of divergence and population demography. We generated spatial-temporal reconstructions using Bayesian inference in a coalescent framework in combination with hind-cast projections of species distribution models (SDMs). We also used the results and literature information to test alternative diversification scenarios via approximate Bayesian computation (ABC). Our results show that Quaternary climatic fluctuations influenced the geographic ranges of both species showing population expansion during the last glacial maximum (LGM) and range contraction during interglacial periods, as inferred from selected ABC models and from past projections of SDMs. We recovered Pleistocene diversification for both species occuring in distinctly unique periods for each taxon. An older and range-restricted lineage was recovered in a geographically isolated geological massif, deserving conservation and further taxonomic study. The diversification and distribution of these amphibian species endemic to the Neotropical 'campo rupestre' were influenced by Quaternary climatic fluctuations. The expansion of cold adapted species restricted to higher elevations during glacial periods and their concomitant retraction during interglacial periods may have been crucial for producing patterns of species richness and endemism along elevation gradients in tropical and subtropical domains. Such processes may influence the evolution of the biota distributed in heterogeneous landscapes with varied topography.

An escape-to-radiate model for explaining the high plant diversity and endemism in campos rupestres†

With extraordinary levels of plant diversity and endemism, the Brazilian campos rupestres across the Espinhaço Range have a species/area ratio 40 times higher than the lowland Amazon. Although diversification drivers in campos rupestres remain a matter of debate, the Pleistocene refugium hypothesis (PRH) is often adopted as the most plausible explanation for their high diversity. The PRH has two main postulates: highland interglacial refugia and a species pump mechanism catalysed by climatic changes. We critically assessed studies on campos rupestres diversification at different evolutionary levels and conclude that most of them are affected by sampling biases, unrealistic assumptions or inaccurate results that do not support the PRH. By modelling the palaeo-range of campos rupestres based on the distribution of 1123 species of vascular plants endemic to the Espinhaço Range and using climate and edaphic variables, we projected a virtually constant suitable area for campos rupestres across the last glacial cycle. We challenge the great importance placed on Pleistocene climatic oscillations in campos rupestres plant diversification and offer an alternative explanation named escape-to-radiate model, which emphasizes niche shifts. Under this biogeographic model of diversification, the long-term fragmentation of campos rupestres combined with recurrent extinctions after genetic drift and sporadic events of adaptive radiation may provide an explanation for the current diversity and endemism in the Espinhaço Range. We conclude that long-term diversification dynamics in campos rupestres are mainly driven by selection, while most endemic diversity is ephemeral, extremely fragile and mainly driven by drift.

Iron islands in the Amazon: investigating plant beta diversity of canga outcrops

The world's largest mineral iron province, Serra dos Carajás, is home to an open vegetation known as canga, found on top of isolated outcrops rising out of the Amazon rainforest. Over one thousand vascular plants species have been recorded in these canga sites, including 38 edaphic endemics. A new survey adds to our investigation of biogeographic relationships between sixteen canga outcrops and the effect of the distance between site pairs on the number of shared species, regional species turnover and species distribution patterns. Plant collecting expeditions to the westernmost site, the Serra de Campos of São Félix do Xingu (SFX), were carried out followed by the identification of all collected specimens and the creation of a species database, built to perform biogeographical analyses. Floristic relationships among the sites were investigated regarding their similarity, using multivariate analyses. The correlation between canga areas and species richness was tested, as well as the geographical distance between pairs of outcrops and their shared species. Vascular plants at SFX total 254 species including 17 edaphic endemics. All canga sites are grouped with 25% of minimum similarity, and the SFX falls within a large subgroup of outcrops. The total species number shared between site pairs does not change significantly with geographical distance but is positively correlated with the area of each outcrop. Meanwhile, shared endemic species numbers between site pairs decline when geographical distance increases, possibly imposed by the barrier of the rainforest. Our data suggest higher shared similarity between the largest and species-richest sites as opposed to geographically nearby sites, and provide useful insight for drafting conservation and compensation measures for canga locations. The size of the canga outcrops is associated to higher floristic diversity but connectivity among islands also plays a role in their similarity.

From micro- to macroevolution: insights from a Neotropical bromeliad with high population genetic structure adapted to rock outcrops

Geographic isolation and reduced population sizes can lead to local extinction, low efficacy of selection and decreased speciation. However, population differentiation is an essential step of biological diversification. In allopatric speciation, geographically isolated populations differentiate and persist until the evolution of reproductive isolation and ecological divergence completes the speciation process. Pitcairnia flammea allows us to study the evolutionary consequences of habitat fragmentation on naturally disjoint rock-outcrop species from the Brazilian Atlantic Rainforest (BAF). Our main results showed low-to-moderate genetic diversity within populations, and deep population structuring caused by limited gene flow, low connectivity, genetic drift and inbreeding of long-term isolation and persistence of rock-outcrop populations throughout Quaternary climatic oscillations. Bayesian phylogenetic and model-based clustering analyses found no clear northern and southern phylogeographic structure commonly reported for many BAF organisms. Although we found two main lineages diverging by ~2 Mya during the early Pleistocene, species' delimitation analysis assigned most of the populations as independent evolving entities, suggesting an important role of disjoint rock outcrops in promoting high endemism in this rich biome. Lastly, we detected limited gene flow in sympatric populations although some hybridization and introgression were observed, suggesting a continuous speciation process in this species complex. Our data not only inform us about the extensive differentiation and limited gene flow found among Pitcairnia flammea species complex, but they also contain information about the mechanisms that shape the genetic architecture of small and fragmented populations of isolated rock outcrop of recently radiated plants.

Phylogenetic Relationships and Evolutionary Trends in the Cactus Family

Members of the cactus family are keystone species of arid and semiarid biomes in the Americas, as they provide shelter and resources to support other members of ecosystems. Extraordinary examples are the several species of flies of the genus Drosophila that lay eggs and feed in their rotting stems, which provide a model system for studying evolutionary processes. Although there is significant progress in understanding the evolution of Drosophila species, there are gaps in our knowledge about the cactus lineages hosting them. Here, we review the current knowledge about the evolution of Cactaceae, focusing on phylogenetic relationships and trends revealed by the study of DNA sequence data. During the last several decades, the availability of molecular phylogenies has considerably increased our understanding of the relationships, biogeography, and evolution of traits in the family. Remarkably, although succulent cacti have very low growth rates and long generation times, they underwent some of the fastest diversifications observed in the plant kingdom, possibly fostered by strong ecological interactions. We have a better understanding of the reproductive biology, population structure and speciation mechanisms in different clades. The recent publication of complete genomes for some species has revealed the importance of phenomena such as incomplete lineage sorting. Hybridization and polyploidization are common in the family, and have been studied using a variety of phylogenetic methods. We discuss potential future avenues for research in Cactaceae, emphasizing the need of a concerted effort among scientists in the Americas, together with the analyses of data from novel sequencing techniques.

Dispersal and local persistence shape the genetic structure of a widespread Neotropical plant species with a patchy distribution

BACKGROUND AND AIMS

Isolated populations constitute an ideal laboratory to study the consequences of intraspecific divergence, because intrinsic incompatibilities are more likely to accumulate under reduced gene flow. Here, we use a widespread bromeliad with a patchy distribution, Pitcairnia lanuginosa, as a model to infer processes driving Neotropical diversification and, thus, to improve our understanding of the origin and evolutionary dynamics of biodiversity in this highly speciose region.

METHODS

We assessed the timing of lineage divergence, genetic structural patterns and historical demography of P. lanuginosa, based on microsatellites, and plastid and nuclear sequence data sets using coalescent analyses and an Approximate Bayesian Computation framework. Additionally, we used species distribution models (SDMs) to independently estimate potential changes in habitat suitability.

KEY RESULTS

Despite morphological uniformity, plastid and nuclear DNA data revealed two distinct P. lanuginosa lineages that probably diverged through dispersal from the Cerrado to the Central Andean Yungas, following the final uplift of the Andes, and passed through long-term isolation with no evidence of migration. Microsatellite data indicate low genetic diversity and high levels of inbreeding within populations, and restricted gene flow among populations, which are likely to be a consequence of bottlenecks (or founder events), and high selfing rates promoting population persistence in isolation. SDMs showed a slight expansion of the suitable range for P. lanuginosa lineages during the Last Glacial Maximum, although molecular data revealed a signature of older divergence. Pleistocene climatic oscillations thus seem to have played only a minor role in the diversification of P. lanuginosa, which probably persisted through adverse conditions in riparian forests.

CONCLUSIONS

Our results imply drift as a major force shaping the evolution of P. lanuginosa, and suggest that dispersal events have a prominent role in connecting Neotropical open and forest biomes.

Maintaining genetic integrity with high promiscuity: Frequent hybridization with low introgression in multiple hybrid zones of Melocactus (Cactaceae)

Hybridization and introgression between species in contact/hybrid zones provide important insight into the genetic and ecological mechanisms of speciation. Cactaceae represents the most important radiation of true succulent angiosperms in the New World. This diversification continues to date, with species experiencing few intrinsic barriers to gene flow and the frequent occurrence of natural hybridization. Here, we used RAD-Seq single-nucleotide polymorphism (SNP) data to investigate the genetic architecture of hybridization in four hybrid zones hosting Melocactus concinnus and four congeneric species (M. ernestii, M. glaucescens, M. paucispinus, and M. zehntneri). Our results revealed that M. concinnus is highly promiscuous in sympatric areas and hybridizes with various species distributed in Morro do Chapéu (Diamantina Plateau, Bahia), eastern Brazil. However, the contemporary genomic introgression among the investigated species is very low (c. 2-5%), confirming that even in the face of hybridization, Melocactus species maintain their genetic integrity. The genomic cline approach showed a large fraction of loci deviating from a model of neutral introgression, where most of the loci are consistent with selection favoring parental genotypes. Our results suggest the occurrence of weak premating but strong postmating reproductive isolation in the analyzed cactus species. Furthermore, as most of the Melocactus species are restricted in distribution, hybridization might negatively affect their integrity if hybrids replace the parental species.

Response of xerophytic plants to glacial cycles in southern South America

BACKGROUND AND AIM

Quaternary glaciations strongly affected the distribution of species from arid and semi-arid environments, as temperature drops were accompanied by strong fluctuations in rainfall. In this study, we examined the response of xerophytic species to glacial cycles, determining the genetic patterns and climatic niche of Echinopsis chiloensis var. chiloensis, an endemic columnar cactus of arid and semi-arid regions of Chile.

METHODS

We analysed 11 polymorphic microsatellites for 130 individuals from 13 populations distributed across the entire distribution of the species. We examined genetic diversity and structure, identified possible patterns of isolation by distance (IBD) and tested two competing population history scenarios using Approximate Bayesian Computation. The first scenario assumes a constant population size while the second includes a bottleneck in the southern population. The latter scenario assumed that the southernmost populations experienced a strong contraction during glaciation, followed by a postglacial expansion; by contrast, the area of the northernmost populations remained as a stable refugium. We also used ecological niche modelling (ENM) to evaluate the location and extension of suitable areas during the Last Glacial Maximum (LGM) and the mid-Holocene.

KEY RESULTS

We found a decline in genetic diversity towards high latitudes and a significant IBD pattern that together with ENM predictions suggest that E. chiloensis var. chiloensis experienced range contraction northwards during wet-cold conditions of the LGM, followed by expansion during aridification of the mid-Holocene. In addition to IBD, we detected the presence of a strong barrier to gene flow at 32°30'S, which according to coalescence analysis occurred 44 kyr BP. The resulting genetic clusters differed in realized climatic niche, particularly in the variables related to precipitation.

CONCLUSIONS

Our results suggest that the cactus E. chiloensis var. chiloensis experienced range contraction and fragmentation during the wet-cold conditions of the LGM, which may have facilitated ecological differentiation between northern and southern populations, promoting incipient speciation.

Pleistocene diversification in an ancient lineage: a role for glacial cycles in the evolutionary history of Dioon Lindl. (Zamiaceae)

PREMISE OF THE STUDY

Recent estimates of crown ages for cycad genera (Late Miocene) challenge us to consider what processes have produced the extant diversity of this ancient group in such relatively little time. Pleistocene climate change has driven major shifts in species distributions in Mexico and may have led to speciation in the genus Dioon by forcing populations to migrate up in elevation, thereby becoming separated by topography.

METHODS

We inferred orthologs from transcriptomes of five species and sequenced these in 42 individuals representing all Dioon species. From these data and published plastid sequences, we inferred dated species trees and lineage-specific diversification rates.

KEY RESULTS

Analyses of 84 newly sequenced nuclear orthologs and published plastid data confirm four major clades within Dioon, all of Pleistocene age. Gene tree analysis, divergence dates, and an increase in diversification rate support very recent and rapid divergence of extant taxa.

CONCLUSIONS

This study confirms the Pleistocene age of Dioon species and implicates Pleistocene climate change and established topography in lineage spitting. These results add to our understanding of the cycads as evolutionarily dynamic lineages, not relicts or evolutionary dead ends. We also find that well-supported secondary calibration points can be reliable in the absence of fossils. Our hypothesis of lineage splitting mediated by habitat shifts may be applicable to other taxa that are restricted to elevation specific ecotones.

Genetic and Historical Colonization Analyses of an Endemic Savanna Tree, Qualea grandiflora, Reveal Ancient Connections Between Amazonian Savannas and Cerrado Core

The evolutionary processes underlying the high diversity and endemism in the Cerrado, the most extensive Neotropical savanna, remain unclear, including the factors promoting the presence and evolution of savanna enclaves in the Amazon forest. In this study, we investigated the effects of past climate changes on genetic diversity, dynamics of species range and the historical connections between the savanna enclaves and Cerrado core for Qualea grandiflora, a tree species widely distributed in the biome. Totally, 40 populations distributed in the Cerrado core and Amazon savannas were analyzed using chloroplast and nuclear DNA sequences. We used phylogeographic, coalescent and ecological niche modeling approaches. Genetic data revealed a phylogeographic structure shaped by Pleistocene climatic oscillations. An eastern-western split in the Cerrado core was observed. The central portion of the Cerrado core harbored most of the sampled diversity for cpDNA. Ecological niche models predicted the presence of a large historical refuge in this region and multiple small refuges in peripheral areas. Relaxed Random Walk (RRW) models indicated the ancestral population in the north-western border of the central portion of the Cerrado core and cyclical dynamics of colonization related to Pleistocene climatic oscillations. Central and western ancient connections between Cerrado core and Amazonian savannas were observed. No evidence of connections among the Amazonian savannas was detected. Our study highlights the importance of Pleistocene climatic oscillations for structuring the genetic diversity of Q. grandiflora and complex evolutionary history of ecotonal areas in the Cerrado. Our results do not support the recent replacement of a large area in the Amazon forest by savanna vegetation. The Amazonian savannas appear to be fragmented and isolated from each other, evolving independently a long ago.

Highly structured genetic diversity of Bixa orellana var. urucurana, the wild ancestor of annatto, in Brazilian Amazonia

Annatto (Bixa orellana L.) is a tropical American crop, commercially valuable due to its application in the food and cosmetics industries as a natural dye. The wild ancestor of cultivated annatto is B. orellana var. urucurana. Although never cultivated, this variety occurs in open forests and anthropogenic landscapes, and is always associated with riparian environments. In this study, we evaluated the genetic diversity and structure of B. orellana var. urucurana populations in Brazilian Amazonia using 16 microsatellite loci. We used Ecological Niche Modeling (ENM) to characterize the potential geographical range of this variety in northern South America. We analyzed 170 samples from 10 municipalities in the states of Rondônia, Pará and Roraima. A total of 194 alleles was observed, with an average of 12.1 alleles per locus. Higher levels of expected (H_E) than observed (H_O) heterozygosities were found for all populations. Bayesian analysis, Neighbor-Joining dendrograms and PCAs suggest the existence of three strongly structured groups of populations. A strong and positive correlation between genetic and geographic distances was found, suggesting that genetic differentiation might be caused by geographic isolation. From species distribution modelling, we detected that South Rondônia, Madre di Dios River basin, Llanos de Mojos, Llanos de Orinoco and eastern Ecuador are highly suitable areas for wild annatto to occur, providing additional targets for future exploration and conservation. Climatic adaptation analyses revealed strong differentiation among populations, suggesting that precipitation plays a key role in wild annatto’s current and potential distribution patterns.

Weak population structure and no genetic erosion in Pilosocereus aureispinus: A microendemic and threatened cactus species from eastern Brazil

Succulent cacti (Cactaceae) are among the most threatened taxonomic groups assessed to date. Here we evaluated the genetic diversity and population structure of a narrow endemic columnar cactus Pilosocereus aureispinus. This species is only found in a small area of c. 300 km² of rocky savanna from eastern Brazil and it is currently classified as vulnerable (VU) on the International Union for Conservation of Nature (IUCN) red list. Eight microsatellite loci were genotyped for 91 individuals from four localities of the known P. aureispinus range. In contrast with expectations for narrow endemic species, we found relatively high levels of genetic diversity (e.g., H_E = 0.390 to 0.525; H_O = 0.394 to 0.572) and very low population structure based on the variation of six loci. All the analyzed individuals were clustered in one unique genetic group in assignment tests. We also generated the sequences of two plastid markers (trnT-trnL and psbD-trnT) and found no variation on a subsample of 39 individuals. We used Landsat 8 images and Normalized Difference Vegetation Index to estimate a potential extent of occurrence of c. 750 km² for this species. Our results showed that P. aureispinus is not suffering from erosion of nuclear genetic variability due to its narrow distribution. However, we advocate that because of the extremely limited extent of occurrence, the ongoing anthropogenic disturbances in its habitat, and phylogenetic distinctiveness of P. aureispinus, this species should be classified as endangered (EN) on the IUCN Red List.

Species delimitation in the Stenocereus griseus (Cactaceae) species complex reveals a new species, S. huastecorum

The Stenocereus griseus species complex (SGSC) has long been considered taxonomically challenging because the number of taxa belonging to the complex and their geographical boundaries remain poorly understood. Bayesian clustering and genetic distance-based methods were used based on nine microsatellite loci in 377 individuals of three main putative species of the complex. The resulting genetic clusters were assessed for ecological niche divergence and areolar morphology, particularly spination patterns. We based our species boundaries on concordance between genetic, ecological, and morphological data, and were able to resolve four species, three of them corresponding to S. pruinosus from central Mexico, S. laevigatus from southern Mexico, and S. griseus from northern South America. A fourth species, previously considered to be S. griseus and commonly misidentified as S. pruinosus in northern Mexico showed significant genetic, ecological, and morphological differentiation suggesting that it should be considered a new species, S. huastecorum, which we describe here. We show that population genetic analyses, ecological niche modeling, and morphological studies are complementary approaches for delimiting species in taxonomically challenging plant groups such as the SGSC.

The xeric side of the Brazilian Atlantic Forest: The forces shaping phylogeographic structure of cacti

In order to investigate biogeographic influences on xeric biota in the Brazilian Atlantic Forest (BAF), a biodiversity hotspot, we used a monophyletic group including three cactus taxa as a model to perform a phylogeographic study: Cereus fernambucensis subsp. fernambucensis, C. fernambucensis subsp. sericifer, and C. insularis. These cacti are allopatric and grow in xeric habitats along BAF, including isolated granite and gneiss rock outcrops (Inselbergs), sand dune vegetation (Restinga forest), and the rocky shore of an oceanic archipelago (islands of Fernando de Noronha). The nucleotide information from nuclear gene phytochrome C and plastid intergenic spacer trnS-trnG was used to perform different approaches and statistical analyses, comprising population structure, demographic changes, phylogenetic relationships, and biogeographic reconstruction in both spatial and temporal scales. We recovered four allopatric population groups with highly supported branches in the phylogenetic tree with divergence initiated in the middle Pleistocene: southern distribution of C. fernambucensis subsp. fernambucensis, northern distribution of C. fernambucensis subsp. fernambucensis together with C. insularis, southern distribution of C. fernambucensis subsp. sericifer, and northern distribution of C. fernambucensis subsp. sericifer. Further, the results suggest that genetic diversity of population groups was strongly shaped by an initial colonization event from south to north followed by fragmentation. The phylogenetic pattern found for C. insularis is plausible with peripatric speciation in the archipelago of Fernando de Noronha. To explain the phylogeographic patterns, the putative effects of both climatic and sea level changes as well as neotectonic activity during the Pleistocene are discussed.

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.

Genetic diversity and differentiation patterns in Micromeria from the Canary Islands are congruent with multiple colonization dynamics and the establishment of species syngameons

Background

Especially on islands closer to the mainland, such as the Canary Islands, different lineages that originated by multiple colonization events could have merged by hybridization, which then could have promoted radiation events (Herben et al., J Ecol 93: 572–575, 2005; Saunders and Gibson, J Ecol 93: 649–652, 2005; Caujapé-Castells, Jesters, red queens, boomerangs and surfers: a molecular outlook on the diversity of the Canarian endemic flora, 2011). This is an alternative to the scenario where evolution is mostly driven by drift (Silvertown, J Ecol 92: 168–173, 2004; Silvertown et al., J Ecol 93: 653–657, 2005). In the former case hybridization should be reflected in the genetic structure and diversity patterns of island species. In the present work we investigate Micromeria from the Canary Islands by extensively studying their phylogeographic pattern based on 15 microsatellite loci and 945 samples. These results are interpreted according to the hypotheses outlined above.

Results

Genetic structure assessment allowed us to genetically differentiate most Micromeria species and supported their current classification. We found that populations on younger islands were significantly more genetically diverse and less differentiated than those on older islands. Moreover, we found that genetic distance on younger islands was in accordance with an isolation-by-distance pattern, while on the older islands this was not the case. We also found evidence of introgression among species and islands.

Conclusions

These results are congruent with a scenario of multiple colonizations during the expansion onto new islands. Hybridization contributes to the grouping of multiple lineages into highly diverse populations. Thus, in our case, islands receive several colonization events from different sources, which are combined into sink populations. This mechanism is in accordance with the surfing syngameon hypothesis. Contrary to the surfing syngameon current form, our results may reflect a slightly different effect: hybridization might always be related to colonization within the archipelago as well, making initial genetic diversity to be high to begin with. Thus the emergence of new islands promotes multiple colonization events, contributing to the establishment of hybrid swarms that may enhance adaptive ability and radiation events. With time, population sizes grow and niches start to fill. Consequently, gene-flow is not as effective at maintaining the species syngameon, which allows genetic differentiation and reproductive isolation to be established between species. This process contributes to an even further decrease in gene-flow between species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1031-y) contains supplementary material, which is available to authorized users.

Alternative glacial-interglacial refugia demographic hypotheses tested on Cephalocereus columna-trajani (Cactaceae) in the intertropical Mexican drylands

Historic demography changes of plant species adapted to New World arid environments could be consistent with either the Glacial Refugium Hypothesis (GRH), which posits that populations contracted to refuges during the cold-dry glacial and expanded in warm-humid interglacial periods, or with the Interglacial Refugium Hypothesis (IRH), which suggests that populations contracted during interglacials and expanded in glacial times. These contrasting hypotheses are developed in the present study for the giant columnar cactus Cephalocereus columna-trajani in the intertropical Mexican drylands where the effects of Late Quaternary climatic changes on phylogeography of cacti remain largely unknown. In order to determine if the historic demography and phylogeographic structure of the species are consistent with either hypothesis, sequences of the chloroplast regions psbA-trnH and trnT-trnL from 110 individuals from 10 populations comprising the full distribution range of this species were analysed. Standard estimators of genetic diversity and structure were calculated. The historic demography was analysed using a Bayesian approach and the palaeodistribution was derived from ecological niche modelling to determine if, in the arid environments of south-central Mexico, glacial-interglacial cycles drove the genetic divergence and diversification of this species. Results reveal low but statistically significant population differentiation (F_ST = 0.124, P < 0.001), although very clear geographic clusters are not formed. Genetic diversity, haplotype network and Approximate Bayesian Computation (ABC) demographic analyses suggest a population expansion estimated to have taken place in the Last Interglacial (123.04 kya, 95% CI 115.3–130.03). The species palaeodistribution is consistent with the ABC analyses and indicates that the potential area of palaedistribution and climatic suitability were larger during the Last Interglacial and Holocene than in the Last Glacial Maximum. Overall, these results suggest that C. columna-trajani experienced an expansion following the warm conditions of interglacials, in accordance with the GRH.

Demographic stability and high historical connectivity explain the diversity of a savanna tree species in the Quaternary

BACKGROUND AND AIMS

Cyclic glaciations were frequent throughout the Quaternary and this affected species distribution and population differentiation worldwide. The present study reconstructed the demographic history and dispersal routes of Eugenia dysenterica lineages and investigated the effects of Quaternary climate change on its spatial pattern of genetic diversity.

METHODS

A total of 333 individuals were sampled from 23 populations and analysed by sequencing four regions of the chloroplast DNA and the internal transcribed spacer of the nuclear DNA. The analyses were performed using a multi-model inference approach based on ecological niche modelling and statistical phylogeography.

KEY RESULTS

Coalescent simulation showed that population stability through time is the most likely scenario. The palaeodistribution dynamics predicted by the ecological niche models revealed that the species was potentially distributed across a large area, extending over Central-Western Brazil through the last glaciation. The lineages of E. dysenterica dispersed from Central Brazil towards populations at the northern, western and south-eastern regions. A historical refugium through time may have favoured lineage dispersal and the maintenance of genetic diversity.

CONCLUSIONS

The results suggest that the central region of the Cerrado biome is probably the centre of distribution of E. dysenterica and that the spatial pattern of its genetic diversity may be the outcome of population stability throughout the Quaternary. The lower genetic diversity in populations in the south-eastern Cerrado biome is probably due to local climatic instability during the Quaternary.

Phylogeny strongly drives seed dormancy and quality in a climatically buffered hotspot for plant endemism

BACKGROUND AND AIMS

Models of costs and benefits of dormancy (D) predict that the evolutionarily stable strategy in long-term stable environments is for non-dormancy (ND), but this prediction remains to be tested empirically. We reviewed seed traits of species in the climatically buffered, geologically stable and nutrient-impoverished campo rupestre grasslands in Brazil to test the hypothesis that ND is favoured over D. We examined the relative importance of life-history traits and phylogeny in driving the evolution of D and assessed seed viability at the community level.

METHODS

Germination and viability data were retrieved from 67 publications and ND/D was determined for 168 species in 25 angiosperm families. We also obtained the percentage of embryoless, viable and dormant seeds for 74 species. Frequencies of species with dormant and non-dormant seeds were compared with global databases of dormancy distribution.

KEY RESULTS

The majority of campo rupestre taxa (62·5 %) had non-dormant seeds, and the ND/D ratio was the highest for any vegetation type on Earth. Dormancy was unrelated to other species life-history traits, suggesting that contemporary factors are poor predictors of D. We found a significant phylogenetic structure in the dormancy categorical trait. Dormancy diversity was highly skewed towards the root of the phylogenetic tree and there was a strong phylogenetic signal in the data, suggesting a major role of phylogeny in determining the evolution of D versus ND and seed viability. Quantitative analysis of the data revealed that at least half of the seeds produced by 46 % of the surveyed populations were embryoless and/or otherwise non-viable.

CONCLUSIONS

Our results support the view that long-term climatic and geological stability favour ND. Seed viability data show that campo rupestre species have a markedly low investment in regeneration from seeds, highlighting the need for specific in situ and ex situ conservation strategies to avoid loss of biodiversity.

Climatic drivers of leaf traits and genetic divergence in the tree Annona crassiflora: a broad spatial survey in the Brazilian savannas

The Cerrado is the largest South American savanna and encompasses substantial species diversity and environmental variation. Nevertheless, little is known regarding the influence of the environment on population divergence of Cerrado species. Here, we searched for climatic drivers of genetic (nuclear microsatellites) and leaf trait divergence in Annona crassiflora, a widespread tree in the Cerrado. The sampling encompassed all phytogeographic provinces of the continuous area of the Cerrado and included 397 individuals belonging to 21 populations. Populations showed substantial genetic and leaf trait divergence across the species' range. Our data revealed three spatially defined genetic groups (eastern, western and southern) and two morphologically distinct groups (eastern and western only). The east-west split in both the morphological and genetic data closely mirrors previously described phylogeographic patterns of Cerrado species. Generalized linear mixed effects models and multiple regression analyses revealed several climatic factors associated with both genetic and leaf trait divergence among populations of A. crassiflora. Isolation by environment (IBE) was mainly due to temperature seasonality and precipitation of the warmest quarter. Populations that experienced lower precipitation summers and hotter winters had heavier leaves and lower specific leaf area. The southwestern area of the Cerrado had the highest genetic diversity of A. crassiflora, suggesting that this region may have been climatically stable. Overall, we demonstrate that a combination of current climate and past climatic changes have shaped the population divergence and spatial structure of A. crassiflora. However, the genetic structure of A. crassiflora reflects the biogeographic history of the species more strongly than leaf traits, which are more related to current climate.

Genetic structure of Micromeria (Lamiaceae) in Tenerife, the imprint of geological history and hybridization on within-island diversification

Geological history of oceanic islands can have a profound effect on the evolutionary history of insular flora, especially in complex islands such as Tenerife in the Canary Islands. Tenerife results from the secondary connection of three paleo‐islands by a central volcano, and other geological events that further shaped it. This geological history has been shown to influence the phylogenetic history of several taxa, including genus Micromeria (Lamiaceae). Screening 15 microsatellite markers in 289 individuals representing the eight species of Micromeria present in Tenerife, this study aims to assess the genetic diversity and structure of these species and its relation with the geological events on the island. In addition, we evaluate the extent of hybridization among species and discuss its influence on the speciation process. We found that the species restricted to the paleo‐islands present lower levels of genetic diversity but the highest levels of genetic differentiation suggesting that their ranges might have contracted over time. The two most widespread species in the island, M. hyssopifolia and M. varia, present the highest genetic diversity levels and a genetic structure that seems correlated with the geological composition of the island. Samples from M. hyssopifolia from the oldest paleo‐island, Adeje, appear as distinct while samples from M. varia segregate into two main clusters corresponding to the paleo‐islands of Anaga and Teno. Evidence of hybridization and intraspecific migration between species was found. We argue that species boundaries would be retained despite hybridization in response to the habitat's specific conditions causing postzygotic isolation and preserving morphological differentiation.

Model-based analysis supports interglacial refugia over long-dispersal events in the diversification of two South American cactus species

Pilosocereus machrisii and P. aurisetus are cactus species within the P. aurisetus complex, a group of eight cacti that are restricted to rocky habitats within the Neotropical savannas of eastern South America. Previous studies have suggested that diversification within this complex was driven by distributional fragmentation, isolation leading to allopatric differentiation, and secondary contact among divergent lineages. These events have been associated with Quaternary climatic cycles, leading to the hypothesis that the xerophytic vegetation patches which presently harbor these populations operate as refugia during the current interglacial. However, owing to limitations of the standard phylogeographic approaches used in these studies, this hypothesis was not explicitly tested. Here we use Approximate Bayesian Computation to refine the previous inferences and test the role of different events in the diversification of two species within P. aurisetus group. We used molecular data from chloroplast DNA and simple sequence repeats loci of P. machrisii and P. aurisetus, the two species with broadest distribution in the complex, in order to test if the diversification in each species was driven mostly by vicariance or by long-dispersal events. We found that both species were affected primarily by vicariance, with a refuge model as the most likely scenario for P. aurisetus and a soft vicariance scenario most probable for P. machrisii. These results emphasize the importance of distributional fragmentation in these species, and add support to the hypothesis of long-term isolation in interglacial refugia previously proposed for the P. aurisetus species complex diversification.

Anonymous nuclear markers data supporting species tree phylogeny and divergence time estimates in a cactus species complex in South America

Supportive data related to the article “Anonymous nuclear markers reveal taxonomic incongruence and long-term disjunction in a cactus species complex with continental-island distribution in South America” (Perez et al., 2016) [1]. Here, we present pyrosequencing results, primer sequences, a cpDNA phylogeny, and a species tree phylogeny.

Lineage-specific evolutionary rate in plants: Contributions of a screening for Cereus (Cactaceae)

PREMISE OF THE STUDY

Predictable chloroplast DNA (cpDNA) sequences have been listed for the shallowest taxonomic studies in plants. We investigated whether plastid regions that vary between closely allied species could be applied for intraspecific studies and compared the variation of these plastid segments with two nuclear regions.

METHODS

We screened 16 plastid and two nuclear intronic regions for species of the genus Cereus (Cactaceae) at three hierarchical levels (species from different clades, species of the same clade, and allopatric populations).

RESULTS

Ten plastid regions presented interspecific variation, and six of them showed variation at the intraspecific level. The two nuclear regions showed both inter- and intraspecific variation, and in general they showed higher levels of variability in almost all hierarchical levels than the plastid segments.

DISCUSSION

Our data suggest no correspondence between variation of plastid regions at the interspecific and intraspecific level, probably due to lineage-specific variation in cpDNA, which appears to have less effect in nuclear data. Despite the heterogeneity in evolutionary rates of cpDNA, we highlight three plastid segments that may be considered in initial screenings in plant phylogeographic studies.

Anonymous nuclear markers reveal taxonomic incongruence and long-term disjunction in a cactus species complex with continental-island distribution in South America

The Pilosocereus aurisetus complex consists of eight cactus species with a fragmented distribution associated to xeric enclaves within the Cerrado biome in eastern South America. The phylogeny of these species is incompletely resolved, and this instability complicates evolutionary analyses. Previous analyses based on both plastid and microsatellite markers suggested that this complex contained species with inherent phylogeographic structure, which was attributed to recent diversification and recurring range shifts. However, limitations of the molecular markers used in these analyses prevented some questions from being properly addressed. In order to better understand the relationship among these species and make a preliminary assessment of the genetic structure within them, we developed anonymous nuclear loci from pyrosequencing data of 40 individuals from four species in the P. aurisetus complex. The data obtained from these loci were used to identify genetic clusters within species, and to investigate the phylogenetic relationship among these inferred clusters using a species tree methodology. Coupled with a palaeodistributional modelling, our results reveal a deep phylogenetic and climatic disjunction between two geographic lineages. Our results highlight the importance of sampling more regions from the genome to gain better insights on the evolution of species with an intricate evolutionary history. The methodology used here provides a feasible approach to develop numerous genealogical molecular markers throughout the genome for non-model species. These data provide a more robust hypothesis for the relationship among the lineages of the P. aurisetus complex.

Using Next Generation RAD Sequencing to Isolate Multispecies Microsatellites for Pilosocereus (Cactaceae)

Microsatellite markers (also known as SSRs, Simple Sequence Repeats) are widely used in plant science and are among the most informative molecular markers for population genetic investigations, but the development of such markers presents substantial challenges. In this report, we discuss how next generation sequencing can replace the cloning, Sanger sequencing, identification of polymorphic loci, and testing cross-amplification that were previously required to develop microsatellites. We report the development of a large set of microsatellite markers for five species of the Neotropical cactus genus Pilosocereus using a restriction-site-associated DNA sequencing (RAD-seq) on a Roche 454 platform. We identified an average of 165 microsatellites per individual, with the absolute numbers across individuals proportional to the sequence reads obtained per individual. Frequency distribution of the repeat units was similar in the five species, with shorter motifs such as di- and trinucleotide being the most abundant repeats. In addition, we provide 72 microsatellites that could be potentially amplified in the sampled species and 22 polymorphic microsatellites validated in two populations of the species Pilosocereus machrisii. Although low coverage sequencing among individuals was observed for most of the loci, which we suggest to be more related to the nature of the microsatellite markers and the possible bias inserted by the restriction enzymes than to the genome size, our work demonstrates that an NGS approach is an efficient method to isolate multispecies microsatellites even in non-model organisms.

Trans-Amazon dispersal potential for Crotalus durissus during Pleistocene climate events

Two disjunct distributional areas of Crotalus durissus(Neotropical rattlesnake) are in open habitats north and south of the Amazon Basin and are presently separated by humid rainforest habitats. We used ecological niche modeling to identify and investigate potential dispersal pathways for this species between the two areas during the late Pleistocene. Niches estimated for the two populations did not differ significantly. Our analyses indicated two possible, but a single most likely, potential routes of dispersal during the last glacial cycle. These results are important to understanding the history of Amazon Basin humid forest biotas, as they suggest agents of isolation among putative humid forest refugia in the form of dry forest and scrub, and associated biotas.