Approximate Bayesian Computation Analysis Research Articles

Europe as a secondary distribution hub in the worldwide invasion of the potato cyst nematode Globodera rostochiensis

The potato cyst nematode Globodera rostochiensis originates from the Andean Mountain region in South America and has unintentionally been introduced to all inhabited continents. Several studies have examined the population genetic structure of this pest in various countries by using microsatellite markers. However, merging microsatellite data produced from different laboratories is challenging and can introduce uncertainty when interpreting the results. To overcome this challenge and to explore invasion routes of this pest, we have genotyped 22 G. rostochiensis populations from all continents. Within populations, the highest genetic diversity was observed in the South American populations, the European populations showed an intermediate level of genetic diversity and the remaining populations were the less diverse. This confirmed pre-existing knowledge such as a first introduction event from South America to Europe, but the less diverse populations could originate either from South America or from Europe. At the continental scale, STRUCTURE genetic clustering output indicated that North America and Asia have experienced at least two introduction events. Comparing different evolutionary scenarios, the Approximate Bayesian Computation analysis showed that Europe served as a secondary distribution centre for the invasion of G. rostochiensis into all other continents (North America, Africa, Asia and Oceania).

A wall-time minimizing parallelization strategy for approximate Bayesian computation.

Approximate Bayesian Computation (ABC) is a widely applicable and popular approach to estimating unknown parameters of mechanistic models. As ABC analyses are computationally expensive, parallelization on high-performance infrastructure is often necessary. However, the existing parallelization strategies leave computing resources unused at times and thus do not optimally leverage them yet. We present look-ahead scheduling, a wall-time minimizing parallelization strategy for ABC Sequential Monte Carlo algorithms, which avoids idle times of computing units by preemptive sampling of subsequent generations. This allows to utilize all available resources. The strategy can be integrated with e.g. adaptive distance function and summary statistic selection schemes, which is essential in practice. Our key contribution is the theoretical assessment of the strategy of preemptive sampling and the proof of unbiasedness. Complementary, we provide an implementation and evaluate the strategy on different problems and numbers of parallel cores, showing speed-ups of typically 10-20% and up to 50% compared to the best established approach, with some variability. Thus, the proposed strategy allows to improve the cost and run-time efficiency of ABC methods on high-performance infrastructure.

Population Structure and Phylogeography of Marine Gastropods Monodonta labio and M. confusa (Trochidae) along the Northwestern Pacific Coast

The genetic structure and demographic history of marine organisms are influenced by biological and ecological features, oceanic currents, and the paleo-geological effects of sea-level fluctuations. In this study, we used mitochondrial COI + 16S gene analysis to investigate the phylogeographic pattern and demography of Monodonta labio and M. confusa in Taiwan, the Ryukyu Islands, Japan, mainland China, and Korea. Our genetic analysis identified two major lineages that were not evident from the distribution patterns of different populations. The Taiwan Strait, which acted as a land bridge during Pleistocene glaciations, and the lack of strong dispersal barriers to gene flow between ocean basins after glaciations shaped the phylogeographic pattern. The genetic differentiation in the Ryukyu Islands was influenced by the specialist-generalist variation hypothesis and the Kuroshio Current. Bayesian skyline plot analyses suggested that the effective population size of M. labio and M. confusa rapidly increased approximately 0.1 and 0.075 million years ago, respectively. Our approximate Bayesian computation analysis suggested that all M. labio and M. confusa populations experienced a decline in population size following a recent population expansion and constant size, respectively. Our study provides a baseline for future investigations of the biogeographical patterns of marine gastropods in the Northwest Pacific and offers valuable insights for the management, sustainable resource utilization, and conservation of this species.

Genetic differentiation and evolution of broad-leaved evergreen shrub and tree varieties of Daphniphyllum macropodum (Daphniphyllaceae).

Tree form evolution is an important ecological specialization for woody species, but its evolutionary process with adaptation is poorly understood, especially on the microevolutionary scale. Daphniphyllum macropodum comprises two varieties: a tree variety growing in a warm temperate climate with light snowfall and a shrub variety growing in a cool temperate climate with heavy snowfall in Japan. Chloroplast DNA variations and genome-wide single-nucleotide polymorphisms across D. macropodum populations and D. teijsmannii as an outgroup were used to reveal the evolutionary process of the shrub variety. Population genetic analysis indicated that the two varieties diverged but were weakly differentiated. Approximate Bayesian computation analysis supported a scenario that assumed migration between the tree variety and the southern populations of the shrub variety. We found migration between the two varieties where the distributions of the two varieties are in contact, and it is concordant with higher tree height in the southern populations of the shrub variety than the northern populations. The genetic divergence between the two varieties was associated with snowfall. The heavy snowfall climate is considered to have developed since the middle Quaternary in this region. The estimated divergence time between the two varieties suggests that the evolution of the two varieties may be concordant with such paleoclimatic change.

Local adaptation and demographic history of vulnerable conifer Cephalotaxus oliveri in southern China

AbstractHaving a comprehensive understanding of genetic differentiation, responses to environmental change and demographic history is critical for genetic improvement and conservation efforts. Forest trees are an excellent resource for understanding population differentiation and adaptive genetic variation due to their ability to adapt to different climates and environments. Cephalotaxus oliveri is a relict conifer endemic to China. In this study, we generated transcriptome data and identified 17 728 high‐quality single‐nucleotide polymorphisms (SNPs) from 18 populations. We found significant negative correlations between expression diversity and nucleotide diversity within and among populations, suggesting that gene expression and nucleotide diversity have a reciprocal relationship when the species adapts to the environment. The analyses of population structure showed that C. oliveri displayed a striking genetic structure with four groups. BayeScEnv and RDA methods detected the signatures of local adaptation, and identified that 738 outlier SNPs were associated with precipitation, temperature and soil conditions across heterogeneous environmental conditions. Approximate Bayesian computation analyses showed that the first and second divergence occurred in the late Miocene (c. 10.075 million years ago [Ma]) and the middle Pleistocene transition (c. 0.815 Ma), respectively. Ecological niche modeling of C. oliveri revealed signs of westward expansion after the last glacial maximum, while it was predicted to experience significant range contractions in future climate change scenarios. Geographical factors and environmental factors in southern China have played a critical role in establishing the current genetic diversity and population structure of C. oliveri. This study provides an important reference for forest resource management and conservation for C. oliveri.

Phylogeography of 912 Cherry Accessions Insight into Independent Origins of Fruiting Cherries and Domestication Footprints of Cultivated Chinese Cherry (Prunus pseudocerasus Lindl.).

The subgenus Cerasus (Rosaceae) contain numerous fruit trees and ornamentals with high economic values. The origin and genetic divergence among various types of fruiting cherries always remain a perplexing issue. We employed three plastom fragments and ITS sequence matrices derived from 912 cherry accessions to elucidate the phylogeographic structure and genetic relationship among fruiting cherries, as well as the origin and domestication of cultivated Chinese cherry. The integration of haplotype genealogies, Approximate Bayesian computation (ABC) approach and estimation of genetic differentiation within and between different groups and lineages has facilitated the resolution of several previously unresolved questions. Firstly, distant phylogenetic relationships between Cerasus and Microcerasus accessions, as indicated by both nuclear and chloroplast data, suggested independent origins and evolution for these two taxa. Moreover, two distinct geographic origin centers (Europe and China) have been confirmed, with significant phylogeographic signals and high genetic differentiation observed between cherries from these regions. This may be attributed to long-term geographic isolation caused by Himalaya-Hengduan Mountains. Our phylogeographic analyses and ABC analysis suggested that cherries inhabiting in China may have undergone multiple hybridization events during the glacial refugia of the eastern edge and southern Himalaya-Hengduan Mountains, followed by rapid radiation throughout their current habitats during interglacial period. The discrepancy between nuclear and chloroplast data may be attributed to hybridization events and incomplete lineage sorting. Furthermore, we speculated that the domesticated Chinese cherries were derived from wild accessions in Longmenshan Fault Zones approximately 2600 years ago. We have also traced the domestication processes and dispersal routes of cultivated Chinese cherries.

Speciation in Coastal Basins Driven by Staggered Headwater Captures: Dispersal of a Species Complex, Leporinus bahiensis, as Revealed by Genome-wide SNP Data.

Past sea level changes and geological instability along watershed boundaries have largely influenced fish distribution across coastal basins, either by dispersal via palaeodrainages now submerged or by headwater captures, respectively. Accordingly, the South American Atlantic coast encompasses several small and isolated drainages that share a similar species composition, representing a suitable model to infer historical processes. Leporinus bahiensis is a freshwater fish species widespread along adjacent coastal basins over narrow continental shelf with no evidence of palaeodrainage connections at low sea level periods. Therefore, this study aimed to reconstruct its evolutionary history to infer the role of headwater captures in the dispersal process. To accomplish this, we employed molecular-level phylogenetic and population structure analyses based on Sanger sequences (5 genes) and genome-wide SNP data. Phylogenetic trees based on Sanger data were inconclusive, but SNPs data did support the monophyletic status of L. bahiensis. Both COI and SNP data revealed structured populations according to each hydrographic basin. Species delimitation analyses revealed from 3 (COI) to 5 (multilocus approach) MOTUs, corresponding to the sampled basins. An intricate biogeographic scenario was inferred and supported by Approximate Bayesian Computation (ABC) analysis. Specifically, a staggered pattern was revealed and characterized by sequential headwater captures from basins adjacent to upland drainages into small coastal basins at different periods. These headwater captures resulted in dispersal throughout contiguous coastal basins, followed by deep genetic divergence among lineages. To decipher such recent divergences, as herein represented by L. bahiensis populations, we used genome-wide SNPs data. Indeed, the combined use of genome-wide SNPs data and ABC method allowed us to reconstruct the evolutionary history and speciation of L. bahiensis. This framework might be useful in disentangling the diversification process in other neotropical fishes subject to a reticulate geological history.

Approximate Bayesian Computation applied to time series of population genetic data disentangles rapid genetic changes and demographic variations in a pathogen population.

Adaptation can occur at remarkably short timescales in natural populations, leading to drastic changes in phenotypes and genotype frequencies over a few generations only. The inference of demographic parameters can allow understanding how evolutionary forces interact and shape the genetic trajectories of populations during rapid adaptation. Here we propose a new Approximate Bayesian Computation (ABC) framework that couples a forward and individual-based model with temporal genetic data to disentangle genetic changes and demographic variations in a case of rapid adaptation. We test the accuracy of our inferential framework and evaluate the benefit of considering a dense versus sparse sampling. Theoretical investigations demonstrate high accuracy in both model and parameter estimations, even if a strong thinning is applied to time series data. Then, we apply our ABC inferential framework to empirical data describing the population genetic changes of the poplar rust pathogen following a major event of resistance overcoming. We successfully estimate key demographic and genetic parameters, including the proportion of resistant hosts deployed in the landscape and the level of standing genetic variation from which selection occurred. Inferred values are in accordance with our empirical knowledge of this biological system. This new inferential framework, which contrasts with coalescent-based ABC analyses, is promising for a better understanding of evolutionary trajectories of populations subjected to rapid adaptation.

The role of geography, ecology, and hybridisation in the evolutionary history of Canary Island Descurainia.

Oceanic islands offer the opportunity to understand evolutionary processes underlying rapid diversification. Along with geographic isolation and ecological shifts, a growing body of genomic evidence has suggested that hybridisation can play an important role in island evolution. Here we use genotyping-by-sequencing (GBS) to understand the roles of hybridisation, ecology, and geographic isolation in the radiation of Canary Island Descurainia (Brassicaceae). We carried out GBS for multiple individuals of all Canary Island species and two outgroups. Phylogenetic analyses of the GBS data were performed using both super-matrix and gene tree approaches and hybridisation events were examined using D-statistics and Approximate Bayesian Computation (ABC). Climatic data were analysed to examine the relationship between ecology and diversification. Analysis of the supermatrix dataset resulted in a fully resolved phylogeny. Species networks suggest a hybridisation event has occurred for D. gilva, with these results being supported by ABC analysis. Strong phylogenetic signals for temperature and precipitation indicate one major ecological shift within Canary Island Descurainia. Inter-island dispersal played a significant role in the diversification of Descurainia, with evidence of only one major shift in climate preferences. Despite weak reproductive barriers and the occurrence of hybrids, hybridisation appears to have played only a limited role in the diversification of the group, with a single instance detected. The results highlight the need to use phylogenetic network approaches that can simultaneously accommodate ILS and gene flow when studying groups prone to hybridisation; patterns that might otherwise be obscured in species trees. This article is protected by copyright. All rights reserved.

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.

Genetic analyses reveal a complex introduction history of the globally invasive tree Acacia longifolia

Acacia longifolia (Sydney golden wattle) is considered one of the most problematic plant invaders in Mediterranean-type ecosystems. In this study, we investigate the species’ invasion history by comparing the genetic diversity and structure of native (Australia) and several invasive range (Brazil, Portugal, South Africa, Spain, and Uruguay) populations and by modelling different introduction scenarios using these data. We sampled 272 A. longifolia individuals – 126 from different invasive ranges and 146 from the native range – from 41 populations. We genotyped all individuals at four chloroplast and 12 nuclear microsatellite markers. From these data we calculated diversity metrics, identified chloroplast haplotypes, and estimated population genetic structure based on Bayesian assignment tests. We used Approximate Bayesian Computation (ABC) models to infer the likely introduction history into each invaded country. In Australia, population genetic structure of A. longifolia appears to be strongly shaped by the Bass Strait and we identified two genetic clusters largely corresponding to mainland Australian and Tasmanian populations. We found invasive populations to represent a mixture of these clusters. Similar levels of genetic diversity were present in native and invasive ranges, indicating that invasive populations did not go through a genetic bottleneck. Bayesian assignment tests and chloroplast haplotype frequencies further suggested a secondary introduction event between South Africa and Portugal. However, ABC analyses could not confidently identify the native source(s) of invasive populations in these two countries, probably due to the known high propagule pressure that accompanied these introductions. ABC analyses identified Tasmania as the likely source of invasive populations in Brazil and Uruguay. A definitive native source for Spanish populations could also not be identified. This study shows that tracing the introduction history of A. longifolia is difficult, most likely because of the complexity associated with the extensive movement of the species around the world. Our findings should be considered when planning management and control efforts, such as biological control, in some invaded regions.

Genetic patterns of Magnolia in the Lesser Antilles: Stepwise colonisation leading to highly inbred island ‘populations’

AbstractAimIslands are ‘nature's laboratories of evolution’. Most island biogeographical studies have focussed on remote volcanic island chains. Here, we challenge island biogeographical patterns using a slowly evolving recent colonist present on five islands in a nonlinear chronosequence island chain.LocationThe Lesser Antilles (Caribbean).TaxonMagnolia dodecapetala (Magnoliaceae).MethodsGenetic diversity was characterised using Sanger sequencing of 21 individuals amplified for 11 DNA markers, plus microsatellite data of 195 individuals genotyped with 19 simple sequence repeat (SSR) markers. Sanger sequencing data were used to construct a Bayesian phylogenetic hypothesis, while SSR markers were used to run approximate Bayesian computation (ABC) demographic analyses and calculate population statistics.ResultsBoth types of molecular data support stepwise colonisation, decoupled from known island ages. The ABC analyses support a north to south migration while the Sanger sequencing data indicate a mixture of island progression rule and stepping stone dispersal. The SSR data show strong genetic structuring per island and significant inbreeding in all populations except in Saint Lucia. The lowest genetic diversity is found in the population from Saint Vincent. A high amount of genetic linkage occurs in a subpopulation from Dominica.Main ConclusionsBiogeographical patterns in the complex geological setting of the Lesser Antilles are uncovered using a slowly evolving study species. All genetic data support treating each island as distinct Management Units for conservation and call for a re‐evaluation of the species limits. Inbreeding threatens the survival of island populations and the populations of Saint Vincent and Dominica represent conservation priorities.

Pre-domestication bottlenecks of the cultivated seaweed Gracilaria chilensis.

Gracilaria chilensis is the main cultivated seaweed in Chile. The low genetic diversity observed in the Chilean populations has been associated with the over-exploitation of natural beds and/or the founder effect that occurred during post-glacial colonization from New Zealand. How these processes have affected its evolutionary trajectory before farming and incipient domestication is poorly understood. In this study, we used 2232 single nucleotide polymorphisms (SNPs) to assess how the species' evolutionary history in New Zealand (its region of origin), the founder effect linked to transoceanic dispersion and colonization of South America, and the recent over-exploitation of natural populations have influenced the genetic architecture of G. chilensis in Chile. The contrasting patterns of genetic diversity and structure observed between the two main islands in New Zealand attest to the important effects of Quaternary glacial cycles on G. chilensis. Approximate Bayesian Computation (ABC) analyses indicated that Chatham Island and South America were colonized independently near the end of the Last Glacial Maximum and emphasized the importance of coastal and oceanic currents during that period. Furthermore, ABC analyses inferred the existence of a recent and strong genetic bottleneck in Chile, matching the period of over-exploitation of the natural beds during the 1970s, followed by rapid demographic expansion linked to active clonal propagation used in farming. Recurrent genetic bottlenecks strongly eroded the genetic diversity of G. chilensis prior to its cultivation, raising important challenges for the management of genetic resources in this incipiently domesticated species.

Contrasting genetic diversity between Planchonella obovata sensu lato (Sapotaceae) on old continental and young oceanic island populations in Japan.

Genetic diversity of plant populations on islands is likely to be influenced by characteristics such as island origin (oceanic or continental) and their age, size, and distance to continental landmasses. In Japan, Planchonella obovata sensu lato which is found on both continental and oceanic islands of varying age, size, and distance to East Asian continental areas-is an ideal system in which to investigate the factors influencing genetic diversity of island plant species. In this study, we examined the genetic diversity of P. obovata s.l. populations, in the context of the species population genetic structure, demography, and between island migration, from 668 individuals, 28 populations and 14 islands including both continental (the Yaeyama Islands) and oceanic islands (the Daito, Bonin, and Volcano Islands) using 11 microsatellite markers. The Yaeyama and Volcano Islands respectively had the highest and lowest genetic diversity, and island origin and age significantly affected genetic diversity. Clustering analysis revealed that populations were grouped into Bonin, Volcano, and Yaeyama + Daito groups. However, Bonin and Volcano groups were distinct despite the relatively short geographical distance between them. Approximate Bayesian Computation analysis suggested that the population size was stable in Bonin and Yaeyama + Daito groups, whereas population reduction occurred in Volcano group, and migration between groups were very limited. Younger oceanic islands showed lower genetic diversity, probably due to limited gene flow and a lack of time to accumulate unique alleles. Genetic structure was generally consistent with the geographic pattern of the islands, but in Volcano, a limited number of founders and limited gene flow among islands are likely to have caused the large genetic divergence observed.

The antique genetic plight of the Mediterranean monk seal (Monachus monachus).

Disentangling the impact of Late Quaternary climate change from human activities can have crucial implications on the conservation of endangered species. We investigated the population genetics and demography of the Mediterranean monk seal (Monachus monachus), one of the world's most endangered marine mammals, through an unprecedented dataset encompassing historical (extinct) and extant populations from the eastern North Atlantic to the entire Mediterranean Basin. We show that Cabo Blanco (Western Sahara/Mauritania), Madeira, Western Mediterranean (historical range) and Eastern Mediterranean regions segregate into four populations. This structure is probably the consequence of recent drift, combined with long-term isolation by distance (R2 = 0.7), resulting from prevailing short-distance (less than 500 km) and infrequent long-distance dispersal (less than 1500 km). All populations (Madeira especially), show high levels of inbreeding and low levels of genetic diversity, seemingly declining since historical time, but surprisingly not being impacted by the 1997 massive die-off in Cabo Blanco. Approximate Bayesian Computation analyses support scenarios combining local extinctions and a major effective population size decline in all populations during Antiquity. Our results suggest that the early densification of human populations around the Mediterranean Basin coupled with the development of seafaring techniques were the main drivers of the decline of Mediterranean monk seals.

Molecular insights into the invasion dynamics of Carcinus crabs in South Africa

Knowledge of the introduction history and spread dynamics of invasive species can provide important insights for management (Ens et al. in Environ Rev (in press), 2022), however such information is often unavailable for accidental introductions. Here we infer how the European shore crab, Carcinus maenas, and its congener, the Mediterranean shore crab, C. aestuarii, were introduced to and spread within South Africa. We do this using nuclear microsatellite data and Bayesian assignment tests and Approximate Bayesian Computation (ABC) modelling that included samples from the native and other invasive ranges of these two species. We also compared the genetic diversity and structure of one of the South African populations during and after intensive management, with that of another, unmanaged, population. South African populations had higher genetic diversity than invasive Carcinus populations from elsewhere in the world. Moreover, the ABC analyses suggest that South African populations originated from an admixture event between individuals of C. maenas from a population in the native range and an invasive population from Canada. We also identified instances of hybridisation between Carcinus maenas and C. aestuarii in South Africa. South African populations showed no genetic structure, suggesting either extensive migration between them or that populations arose from the same initial introduction. Management of Carcinus did not affect genetic diversity or structure, and we suspect that the management duration was insufficient to target a full generation of crabs. Together these results suggest multiple introductions and/or high propagule pressure to South Africa, crab (larval or adult) movement between existing populations, and some hybridisation. For eradication from South Africa to be achieved, management would need to concurrently target all known invasive populations and clearly establish that new introductions could be prevented.

The fast invasion of Europe by the box tree moth: an additional example coupling multiple introduction events, bridgehead effects and admixture events

Identifying the invasion routes of non-native species is crucial to understanding invasions and customizing management strategies. The box tree moth, Cydalima perspectalis, is native to Asia and was recently accidentally introduced into Europe as a result of the ornamental plant trade. Over the last 15 years, it has spread across the continent and has reached the Caucasus and Iran. It is threatening Buxus trees in both urban areas and forests. To investigate the species’ invasion routes, native and invasive box tree moth populations were sampled, and moth’s genetic diversity and structure were compared using microsatellite markers. Our approximate Bayesian computation analyses strongly suggest that invasion pathways were complex. Primary introductions originating from eastern China probably occurred independently twice in Germany and once in the Netherlands. There were also possibly bridgehead effects, where at least three invasive populations may have served as sources for other invasive populations within Europe, with indication of admixture between the two primary invasive populations. The bridgehead populations were likely those in the countries that play a major role in the ornamental plant trade in Europe, notably Germany, the Netherlands, and Italy. All these invasion processes likely facilitated its fast expansion across Europe and illustrate the role played by the ornamental plant trade not only in the moth’s introduction from China but also in the species’ spread across Europe, leading to an invasion with a complex pattern.

Hormonal regulation of ovarian follicle growth in humans: Model-based exploration of cycle variability and parameter sensitivities

We present a modelling and simulation framework for the dynamics of ovarian follicles and key hormones along the hypothalamic-pituitary–gonadal axis throughout consecutive human menstrual cycles. All simulation results (hormone concentrations and ovarian follicle sizes) are in biological units and can easily be compared to clinical data. The model takes into account variability in follicles’ response to stimulating hormones, which introduces variability between cycles. The growth of ovarian follicles in waves is an emergent property in our model simulations and further supports the hypothesis that follicular waves are also present in humans. We use Approximate Bayesian Computation and cluster analysis to construct a population of virtual subjects and to study parameter distributions and sensitivities. The model can be used to compare and optimize treatment protocols for ovarian hyperstimulation, thus potentially forming the integral part of a clinical decision support system in reproductive endocrinology.

Biogeography pattern of the marine angiosperm Cymodocea nodosa in the eastern Mediterranean Sea related to the quaternary climatic changes.

We investigated the population dynamics of a highly clonal marine angiosperm, Cymodocea nodosa, in the eastern Mediterranean Sea, to identify the historical dynamics, demography, and connectivity of the species in the area. Eighteen microsatellite loci were used in conjunction with coalescent methods to investigate the genetic structure and demographic history of C. nodosa meadows. Approximate Bayesian computation (ABC) modeling was used to examine the pattern of divergence over time in the context of environmental change over the course of the Quaternary period. ABC analysis revealed an initial split of the C. nodosa populations between the north‐western, northern, and north‐eastern Aegean Sea during the Pleistocene epoch, followed by a more recent divergence of the north‐western population and the central‐western part of the Aegean Sea. According to the results, the most parsimonious historical scenario is that of a pervasive genetic signature of the effects of the drop in sea level during the Pleistocene epoch. This scenario supports the isolation of the north‐western, north, and north‐eastern area, and the subsequent recolonization after post‐glaciation sea level rise that may explain the north‐western differentiation as well present‐day detected dispersion of C. nodosa.

Evolutionary history and colonization patterns of the wing dimorphic grasshopper Dichroplus vittatus in two Argentinean biomes

Quaternary climate oscillations and modification of the environment by humans have played an important role in shaping species distribution and genetic structure of modern species. Here, population genetic parameters were inferred from the analysis of 168 individuals belonging to 11 populations of the South American grasshopper, Dichroplus vittatus, distributed in two Argentinean Biomes (Grassland and Savanna), by sequencing a 543 bp of the mitochondrial COI gene. Overall, we detected considerable haplotype diversity and low nucleotide diversity. AMOVA analyses showed a significant degree of differentiation among Biomes and between populations. Two major mitochondrial lineages can be distinguished. The haplogroup containing the most common haplotype split 17,000 years BP while the haplogroup including the second most common haplotype has a divergence date of about 11,700 years. Approximate Bayesian Computation (ABC) analyses showed that the palaeodemographic scenario that best fitted our data is consistent with a hypothesis of divergence from an ancestral population and subsequent admixture with Grassland-Savanna (South–North) direction. Our results suggest that populations located in both Biomes would derive from a single ancestral population that colonized the region after the Last Glacial Maximum and Grassland would have a more ancestral origin than Savanna. Further, our results emphasize the importance of human-mediated dispersal in the reconfiguration of genetic diversity of species with potential pest capacity.