The sweet tabaiba or there and back again: phylogeographical history of the Macaronesian Euphorbia balsamifera.

BACKGROUND AND AIMS: Biogeographical relationships between the Canary Islands and north-west Africa are often explained by oceanic dispersal and geographical proximity. Sister-group relationships between Canarian and eastern African/Arabian taxa, the 'Rand Flora' pattern, are rare among plants and have been attributed to the extinction of north-western African populations. Euphorbia balsamifera is the only representative species of this pattern that is distributed in the Canary Islands and north-west Africa; it is also one of few species present in all seven islands. Previous studies placed African populations of E. balsamifera as sister to the Canarian populations, but this relationship was based on herbarium samples with highly degraded DNA. Here, we test the extinction hypothesis by sampling new continental populations; we also expand the Canarian sampling to examine the dynamics of island colonization and diversification. METHODS: Using target enrichment with genome skimming, we reconstructed phylogenetic relationships within E. balsamifera and between this species and its disjunct relatives. A single nucleotide polymorphism dataset obtained from the target sequences was used to infer population genetic diversity patterns. We used convolutional neural networks to discriminate among alternative Canary Islands colonization scenarios. KEY RESULTS: The results confirmed the Rand Flora sister-group relationship between western E. balsamifera and Euphorbia adenensis in the Eritreo-Arabian region and recovered an eastern-western geographical structure among E. balsamifera Canarian populations. Convolutional neural networks supported a scenario of east-to-west island colonization, followed by population extinctions in Lanzarote and Fuerteventura and recolonization from Tenerife and Gran Canaria; a signal of admixture between the eastern island and north-west African populations was recovered. CONCLUSIONS: Our findings support the Surfing Syngameon Hypothesis for the colonization of the Canary Islands by E. balsamifera, but also a recent back-colonization to the continent. Populations of E. balsamifera from northwest Africa are not the remnants of an ancestral continental stock, but originated from migration events from Lanzarote and Fuerteventura. This is further evidence that oceanic archipelagos are not a sink for biodiversity, but may be a source of new genetic variability.

Expanding the stdpopsim species catalog, and lessons learned for realistic genome simulations.

Simulation is a key tool in population genetics for both methods development and empirical research, but producing simulations that recapitulate the main features of genomic datasets remains a major obstacle. Today, more realistic simulations are possible thanks to large increases in the quantity and quality of available genetic data, and the sophistication of inference and simulation software. However, implementing these simulations still requires substantial time and specialized knowledge. These challenges are especially pronounced for simulating genomes for species that are not well-studied, since it is not always clear what information is required to produce simulations with a level of realism sufficient to confidently answer a given question. The community-developed framework stdpopsim seeks to lower this barrier by facilitating the simulation of complex population genetic models using up-to-date information. The initial version of stdpopsim focused on establishing this framework using six well-characterized model species (Adrion et al., 2020). Here, we report on major improvements made in the new release of stdpopsim (version 0.2), which includes a significant expansion of the species catalog and substantial additions to simulation capabilities. Features added to improve the realism of the simulated genomes include non-crossover recombination and provision of species-specific genomic annotations. Through community-driven efforts, we expanded the number of species in the catalog more than threefold and broadened coverage across the tree of life. During the process of expanding the catalog, we have identified common sticking points and developed the best practices for setting up genome-scale simulations. We describe the input data required for generating a realistic simulation, suggest good practices for obtaining the relevant information from the literature, and discuss common pitfalls and major considerations. These improvements to stdpopsim aim to further promote the use of realistic whole-genome population genetic simulations, especially in non-model organisms, making them available, transparent, and accessible to everyone.

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 Genetic Differentiation of Pyrrhulina (Teleostei, Characiformes) Species is Likely Influenced by Both Geographical Distribution and Chromosomal Rearrangements.

Allopatry is generally considered to be one of the main contributors to the remarkable Neotropical biodiversity. However, the role of chromosomal rearrangements including neo-sex chromosomes for genetic diversity is still poorly investigated and understood. Here, we assess the genetic divergence in five Pyrrhulina species using population genomics and combined the results with previously obtained cytogenetic data, highlighting that molecular genetic diversity is consistent with their chromosomal features. The results of a principal coordinate analysis (PCoA) indicated a clear difference among all species while showing a closer relationship of the ones located in the same geographical region. This was also observed in genetic structure analyses that only grouped P. australis and P. marilynae, which were also recovered as sister species in a species tree analysis. We observed a contradictory result for the relationships among the three species from the Amazon basin, as the phylogenetic tree suggested P. obermulleri and P. semifasciata as sister species, while the PCoA showed a high genetic difference between P. semifasciata and all other species. These results suggest a potential role of sex-related chromosomal rearrangements as reproductive barriers between these species.

Congruent evolutionary responses of European steppe biota to late Quaternary climate change.

Quaternary climatic oscillations had a large impact on European biogeography. Alternation of cold and warm stages caused recurrent glaciations, massive vegetation shifts, and large-scale range alterations in many species. The Eurasian steppe biome and its grasslands are a noteworthy example; they underwent climate-driven, large-scale contractions during warm stages and expansions during cold stages. Here, we evaluate the impact of these range alterations on the late Quaternary demography of several phylogenetically distant plant and insect species, typical of the Eurasian steppes. We compare three explicit demographic hypotheses by applying an approach combining convolutional neural networks with approximate Bayesian computation. We identified congruent demographic responses of cold stage expansion and warm stage contraction across all species, but also species-specific effects. The demographic history of the Eurasian steppe biota reflects major paleoecological turning points in the late Quaternary and emphasizes the role of climate as a driving force underlying patterns of genetic variance on the biome level.

Integrating Cytogenetics and Population Genomics: Allopatry and Neo-Sex Chromosomes May Have Shaped the Genetic Divergence in the Erythrinus erythrinus Species Complex (Teleostei, Characiformes).

Diversity found in Neotropical freshwater fish is remarkable. It can even hinder a proper delimitation of many species, with the wolf fish Erythrinus erythrinus (Teleostei, Characiformes) being a notable example. This nominal species shows remarkable intra-specific variation, with extensive karyotype diversity found among populations in terms of different diploid chromosome numbers (2n), karyotype compositions and sex chromosome systems. Here, we analyzed three distinct populations (one of them cytogenetically investigated for the first time) that differed in terms of their chromosomal features (termed karyomorphs) and by the presence or absence of heteromorphic sex chromosomes. We combined cytogenetics with genomic approaches to investigate how the evolution of multiple sex chromosomes together with allopatry is linked to genetic diversity and speciation. The results indicated the presence of high genetic differentiation among populations both from cytogenetic and genomic aspects, with long-distance allopatry potentially being the main agent of genetic divergence. One population showed a neo-X1X2Y sexual chromosome system and we hypothesize that this system is associated with enhanced inter-population genetic differentiation which could have potentially accelerated speciation compared to the effect of allopatry alone.

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.

Multiple sex chromosomes in teleost fishes from a cytogenetic perspective: state of the art and future challenges.

Despite decades of cytogenetic and genomic research of dynamic sex chromosome evolution in teleost fishes, multiple sex chromosomes have been largely neglected. In this review, we compiled available data on teleost multiple sex chromosomes, identified major trends in their evolution and suggest further trajectories in their investigation. In a compiled dataset of 440 verified records of fish sex chromosomes, we counted 75 multiple sex chromosome systems with 60 estimated independent origins. We showed that male-heterogametic systems created by Y-autosome fusion predominate and that multiple sex chromosomes are over-represented in the order Perciformes. We documented a striking difference in patterns of differentiation of sex chromosomes between male and female heterogamety and hypothesize that faster W sex chromosome differentiation may constrain sex chromosome turnover in female-heterogametic systems. We also found no significant association between the mechanism of multiple sex chromosome formation and percentage of uni-armed chromosomes in teleost karyotypes. Last but not least, we hypothesized that interaction between fish populations, which differ in their sex chromosomes, can drive the evolution of multiple sex chromosomes in fishes. This underlines the importance of broader inter-population sampling in studies of fish sex chromosomes. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

High Genetic Diversity despite Conserved Karyotype Organization in the Giant Trahiras from Genus Hoplias (Characiformes, Erythrinidae).

In the fish genus Hoplias, two major general groups can be found, one of which is formed by the "common trahiras" (Hoplias malabaricus group) and the other by the "giant trahiras" (Hoplias lacerdae group, in addition to Hoplias aimara), which usually comprises specimens of larger body size. Previous investigations from the giant trahiras group recovered 2n = 50 meta/submetacentric chromosomes and no sex chromosome differentiation, indicating a probable conservative pattern for their karyotype organization. Here, we conducted comparative cytogenetic studies in six giant trahiras species, two of them for the first time. We employed standard and advanced molecular cytogenetics procedures, including comparative genomic hybridization (CGH), as well as genomic assessments of diversity levels and phylogenetic relationships among them. The results strongly suggest that the giant trahiras have a particular and differentiated evolutionary pathway inside the Hoplias genus. While these species share the same 2n and karyotypes, their congeneric species of the H. malabaricus group show a notable chromosomal diversity in number, morphology, and sex chromosome systems. However, at the same time, significant changes were characterized at their inner chromosomal level, as well as in their genetic diversity, highlighting their current relationships resulting from different evolutionary histories.

Comparative cytogenetic survey of the giant bonytongue Arapaima fish (Osteoglossiformes: Arapaimidae), across different Amazonian and Tocantins/Araguaia River basins

The South American giant fishes of the genus Arapaima, commonly known as pirarucu, are one of the most iconic among Osteoglossiformes. Previously cytogenetic studies have identified their karyotype characteristics; however, characterization of cytotaxonomic differentiation across their distribution range remains unknown. In this study, we compared chromosomal characteristics using conventional and molecular cytogenetic protocols in pirarucu populations from the Amazon and Tocantins-Araguaia river basins to verify if there is differentiation among representatives of this genus. Our data revealed that individuals from all populations present the same diploid chromosome number 2n=56 and karyotype composed of 14 pairs of meta- to submetacentric and 14 pairs of subtelo- to acrocentric chromosomes. The minor and major rDNA sites are in separate chromosomal pairs, in which major rDNA sites corresponds to large heterochromatic blocks. Comparative genomic hybridizations (CGH) showed that the genome of these populations shared a great portion of repetitive elements, due to a lack of substantial specific signals. Our comparative cytogenetic data analysis of pirarucu suggested that, although significant genetic differences occur among populations, their general karyotype patterns remain conserved.(AU)

Os peixes gigantes da América do Sul do gêneroArapaima, comumente conhecidos como pirarucus, são um dos mais icônicos de Osteoglossiformes. Estudos citogenéticos prévios identificaram suas características cariotípicas, entretanto a caracterização da diferenciação citotaxonômica através de suas distribuições geográficas ainda são desconhecidas. Nesse estudo, nós comparamos características cromossômicas utilizando técnicas de citogenética clássica e molecular em populações das bacias dos rios Amazonas e Tocantins-Araguaia, a fim de verificar se há alguma diferenciação entre representantes desse gênero. Nossos dados revelaram que indivíduos de todas as populações apresentam número diploide de 2n=56 cromossomos e que seus cariótipos são compostos de 14 pares de cromossomos meta- e submetacêntricos e 14 pares de subtelo- e acrocêntricos. Os sítios maiores e menores de rDNA estão localizados em pares cromossômicos separados, onde os sítios maiores de rDNA correspondem a grandes blocos heterocromáticos. Hibridizações genômicas comparativas (CGH) mostraram que o genoma dos espécimes dessas populações é amplamente compartilhado, devido à falta de sinais substanciais específicos. Nossos dados de citogenética comparativa do pirarucu sugerem que embora diferenças genéticas significativas ocorram entre populações, os padrões cariotípicos gerais se mantêm conservados.(AU)

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 km2 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., HE = 0.390 to 0.525; HO = 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 km2 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.

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.

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.

Pleistocene niche stability and lineage diversification in the subtropical spider Araneus omnicolor (Araneidae).

The influence of Quaternary climate oscillations on the diversification of the South American fauna is being increasingly explored. However, most of these studies have focused on taxa that are endemic to tropical environments, and relatively few have treated organisms restricted to subtropical biomes. Here we used an integrative phylogeographical framework to investigate the effects of these climate events on the ecological niche and genetic patterns of the subtropical orb-weaver spider Araneus omnicolor (Araneidae). We analyzed the mitochondrial (Cytochrome Oxidase I, COI) and nuclear (Internal Transcribed Subunit II, ITS2) DNA of 130 individuals throughout the species' range, and generated distribution models in three different climate scenarios [present, Last Glacial Maximum (LGM), and Last Interglacial Maximum (LIG)]. Additionally, we used an Approximate Bayesian Computation (ABC) approach to compare possible demographic scenarios and select the hypothesis that better explains the genetic patterns of A. omnicolor. We obtained high haplotype diversity but low nucleotide variation among sequences. The population structure and demographic analyses showed discrepancies between markers, suggesting male-biased dispersal in the species. The time-calibrated COI phylogenetic inference showed a recent diversification of lineages (Middle/Late Pleistocene), while the paleoclimate modeling indicated niche stability since ~120 Kya. The ABC results agreed with the niche models, supporting a panmictic population as the most likely historical scenario for the species. These results indicate that A. omnicolor experienced no niche or population reductions during the Late Pleistocene, despite the intense landscape modifications that occurred in the subtropical region, and that other factors beside LGM and LIG climate oscillations might have contributed to the demographic history of this species. This pattern may be related to the high dispersal ability and wide environmental tolerance of A. omnicolor, highlighting the need for more phylogeographical studies with invertebrates and other generalist taxa, in order to understand the effects of Quaternary climate changes on Neotropical biodiversity.

The evolution of phylogeographic data sets.

Empirical phylogeographic studies have progressively sampled greater numbers of loci over time, in part motivated by theoretical papers showing that estimates of key demographic parameters improve as the number of loci increases. Recently, next-generation sequencing has been applied to questions about organismal history, with the promise of revolutionizing the field. However, no systematic assessment of how phylogeographic data sets have changed over time with respect to overall size and information content has been performed. Here, we quantify the changing nature of these genetic data sets over the past 20 years, focusing on papers published in Molecular Ecology. We found that the number of independent loci, the total number of alleles sampled and the total number of single nucleotide polymorphisms (SNPs) per data set has improved over time, with particularly dramatic increases within the past 5 years. Interestingly, uniparentally inherited organellar markers (e.g. animal mitochondrial and plant chloroplast DNA) continue to represent an important component of phylogeographic data. Single-species studies (cf. comparative studies) that focus on vertebrates (particularly fish and to some extent, birds) represent the gold standard of phylogeographic data collection. Based on the current trajectory seen in our survey data, forecast modelling indicates that the median number of SNPs per data set for studies published by the end of the year 2016 may approach ~20,000. This survey provides baseline information for understanding the evolution of phylogeographic data sets and underscores the fact that development of analytical methods for handling very large genetic data sets will be critical for facilitating growth of the field.

Interglacial microrefugia and diversification of a cactus species complex: phylogeography and palaeodistributional reconstructions for Pilosocereus aurisetus and allies.

The role of Pleistocene climate changes in promoting evolutionary diversification in global biota is well documented, but the great majority of data regarding this subject come from North America and Europe, which were greatly affected by glaciation. The effects of Pleistocene changes on cold- and/or dry-adapted species in tropical areas where glaciers were not present remain sparsely investigated. Many such species are restricted to small areas surrounded by unfavourable habitats, which may represent potential interglacial microrefugia. Here, we analysed the phylogeographic structure and diversification history of seven cactus species in the Pilosocereus aurisetus complex that are restricted to rocky areas with high diversity and endemism within the Neotropical savannas of eastern South America. We combined palaeodistributional estimates with standard phylogeographic approaches based on two chloroplast DNA regions (trnT-trnL and trnS-trnG), exon 1 of the nuclear gene PhyC and 10 nuclear microsatellite loci. Our analyses revealed a phylogeographic history marked by multiple levels of distributional fragmentation, isolation leading to allopatric differentiation and secondary contact among divergent lineages within the complex. Diversification and demographic events appear to have been affected by the Quaternary climatic cycles as a result of isolation in multiple patches of xerophytic vegetation. These small patches presently harbouring P. aurisetus populations seem to operate as microrefugia, both at present and during Pleistocene interglacial periods; the role of such microrefugia should be explored and analysed in greater detail.

Cross-species amplification of microsatellites reveals incongruence in the molecular variation and taxonomic limits of the Pilosocereus aurisetus group (Cactaceae).

The Pilosocereus aurisetus group contains eight cactus species restricted to xeric habitats in eastern and central Brazil that have an archipelago-like distribution. In this study, 5-11 microsatellite markers previously designed for Pilosocereus machrisii were evaluated for cross-amplification and polymorphisms in ten populations from six species of the P. aurisetus group. The genotypic information was subsequently used to investigate the genetic relationships between the individuals, populations, and species analyzed. Only the Pmac101 locus failed to amplify in all of the six analyzed species, resulting in an 88 % success rate. The number of alleles per polymorphic locus ranged from 2 to 12, and the most successfully amplified loci showed at least one population with a larger number of alleles than were reported in the source species. The population relationships revealed clear genetic clustering in a neighbor-joining tree that was partially incongruent with the taxonomic limits between the P. aurisetus and P. machrisii species, a fact which parallels the problematic taxonomy of the P. aurisetus group. A Bayesian clustering analysis of the individual genotypes confirmed the observed taxonomic incongruence. These microsatellite markers provide a valuable resource for facilitating large-scale genetic studies on population structures, systematics and evolutionary history in this group.

Isolation, characterization, and cross-species amplification of polymorphic microsatellite markers for Pilosocereus machrisii (Cactaceae).

PREMISE OF THE STUDY: Microsatellite primers were developed in Pilosocereus machrisii, a columnar cactus with a patchy distribution in eastern tropical South America, to assess its genetic diversity and population structure. METHODS AND RESULTS: Eleven microsatellite markers were developed, of which one was monomorphic among 51 individuals from two populations. The number of alleles per polymorphic locus ranged from two to eight, and the total number of alleles was 57. From the 11 isolated loci, nine were successfully amplified in the other four Pilosocereus species. CONCLUSIONS: The results showed that these markers will be useful for studies of genetic diversity, patterns of gene flow, and population genetic structure in P. machrisii, as well as across other congeneric species.