Inferring population history with DIY ABC: a user-friendly approach to approximate Bayesian computation

Intraspecific character displacement in oaks

Character displacement refers to the process by which species diverge more in sympatry due to competition for resources. This competition-driven speciation can also occur within populations, known as intraspecific character displacement (ICD). ICD can promote divergence within species by influencing intraspecific competition or encouraging the evolution of alternative phenotypes. Despite its significance, ICD remains understudied and requires further exploration. In this study, we investigate how competition influences genetic and morphological differentiation within species in sympatric and allopatric populations. We focused on Quercus serrata (in China and Japan) and Q. serrata var. brevipetiolata (found only in China), which belong to a small monophyletic group of oak species nested within Section Quercus (white oaks). Using genetic markers, we detected divergence between Chinese and Japanese populations and further diversification within China, with asymmetric historical gene flow primarily from Q. serrata (the earlier diverged species) to Q. serrata var. brevipetiolata (the later variety). Although genetic differentiation did not differ between sympatric and allopatric populations, leaf morphological variation, analyzed through the geometric morphometric method (GMM) and traditional morphological method, revealed greater trait variation in sympatry. In addition, we found an allometric growth relationship between leaf size and leaf mass of Q. serrata and Q. serrata var. brevipetiolata, with the leaf area of Q. serrata var. brevipetiolata decreasing more disproportionately to leaf mass. This suggests a resource trade-off, where Q. serrata var. brevipetiolata, the later diverged variety, adopts more resource-conservative traits in sympatry. Further analysis of trait variation with environmental factors supports these findings, while genetic variation along climate gradients showed significant responses primarily in Q. serrata, regardless of sympatric or allopatric conditions. Although neutral genetic markers are insufficient to capture selection-driven adaptive differentiation, we inferred that Q. serrata var. brevipetiolata is progressing towards ecological divergence from Q. serrata. Overall, our results highlight the role of ICD in driving morphological diversification and resource-use strategies within species in response to competitive pressures.

Next-generation sequencing-based population genetics unravels the evolutionary history of Rhodomyrtus tomentosa in China

BACKGROUND

Rhodomyrtus tomentosa (Ait.) Hassk. is useful for its ornamental, medicinal, and ecological characteristics, and is considered a "Neglected and Underutilized Crop Species". However, our understanding of the geographic structure and evolutionary history of its wild populations is limited. To address this gap, we investigated genomic data from 284 samples of R. tomentosa from 28 wild populations in southern China.

RESULTS

The genetic diversity of populations in different regions revealed the similar trends using whole-genome and RAD-seq data, and Hainan Island having a higher genetic diversity than other regions. The 28 populations clustered into three distinct groups: (a) GROUP1 on the eastern mainland within Guangdong, Fujian, and Hunan Provinces; (b) GROUP2 on the western mainland within Guangxi and Yunnan Provinces; and (c) GROUP3 on Hainan Island. Mantel tests and redundancy analyses revealed population differentiation was affected by distance and environmental factors such as annual average radiation. Demographic history and gene flow analyses indicated the mainland populations and the Hainan Island populations diverged around 0.93 MYA, with gene flow primarily occurring from Hainan Island and the coastal regions (such as Zhanjiang in Guangdong and Fangchenggang in Guangxi) towards the mainland, reflecting an expansion trend within the species. PSMC' analyses indicated that the populations of the three groups underwent a bottleneck during the Pleistocene due to glacial-interglacial cycles and geological events. Niche analysis revealed that the ice ages caused habitat contraction for the species, and populations with higher genetic diversity are generally distributed in areas with more suitable habitats.

CONCLUSIONS

This study elucidates the current genetic distribution of the species within China and suggests that drastic Pleistocene climate change and geographical events caused population divergence and fluctuations in effective population size, shaping the current genetic distribution of R. tomentosa. These findings provide a theoretical basis for the genetic conservation and improvement of R. tomentosa.

The genetic trail of the invasive mosquito species Aedes koreicus from the east to the west of Northern Italy

BACKGROUND

Aedes koreicus is native to Far East Asia and recorded in Europe since 2008. In Italy, Ae. koreicus is widespread throughout the Northern part of the peninsula, highlighting its invasive potential and spread. However, no clear clues about the dispersal patterns of the species have been collected so far.

METHODOLOGY/PRINCIPAL FINDINGS

Population genetic analyses were performed to assess the genetic structure of populations of Ae. koreicus and to make hypotheses about its dispersal patterns in Northern Italy. Ten microsatellite markers specific for Ae. koreicus were used to genotype 414 individuals from 13 populations in the pre-alpine area of Italy, and neighboring Slovenia. Basic and Bayesian population genetic analyses were performed to evaluate patterns of genetic variation, genetic structure, and demography of selected mosquito populations. While presenting a certain degree of structuring, the Italian and Slovenian populations of Ae. koreicus were poorly differentiated. Moreover, demographic analysis supports the expansion of a single population propagule of Ae. koreicus in Italy and Slovenia and provides evidence of the presence of overwintering populations in the studied area.

CONCLUSIONS/SIGNIFICANCE

Our results highlight a common origin, and stable colonization of Northern Italy and Slovenia, as a probable consequence of the expansion of a unique population. This stresses out the importance of continuous monitoring of Ae. koreicus, to finally uncover the geographic origins and entrance pathways of invasive populations and to prevent or limit further introductions.

Population genetic structure of Phaedranassa cinerea Ravenna (Amaryllidaceae) and conservation implications

BACKGROUND

Andean orography has shaped the endemism of plant species in montane forests, creating a mosaic of habitats in small and isolated areas. Understanding these endemic species' genetic diversity patterns is crucial for their conservation. Phaedranassa cinerea (Amaryllidaceae), a species restricted to the western Andes of Ecuador, is listed as "vulnerable" according to the IUCN criteria. This study seeks to determine whether there is genetic structure among and within Phaedranassa cinerea populations, estimate the timing of their genetic divergence, and recommend conservation strategies based on these genetic structure findings.

RESULTS

Using 13 microsatellites and a Bayesian approach, we analyzed the genetic differentiation of P. cinerea and possible diversification scenarios. Our results indicate that the genetic diversity of P. cinerea is lower than congeneric species. The Bayesian analysis identified two genetic groups, with no evidence of isolation by distance. Populations in the northwest of the Ecuadorean Andes have less allele richness compared to those in the southwest. Additionally, the species exhibits excess homozygosity and evidence of bottlenecks. Our Bayesian analysis suggests that the differentiation among populations was not older than 5,000 years and was as recent as 600 years ago for some of the populations. Based on the geographic distribution of the known populations, the species should be listed as endangered instead of vulnerable to extinction.

CONCLUSIONS

Phaedranassa cinerea shows lower genetic diversity than related species, with the most variation within populations. We identified two to four genetic groups, suggesting recent divergence along the ridges of the western Andes. The findings suggest that conservation efforts should focus on securing genetic exchange between populations to preserve the genetic diversity of P. cinerea.

Population genetics of the freshwater red alga Batrachospermum gelatinosum (Rhodophyta) II: Phylogeographic analyses reveal spatial genetic structure among and within five major drainage basins in eastern North America

The freshwater red alga Batrachospermum gelatinosum has a well-documented distribution spanning historically glaciated and unglaciated eastern North America. This alga has no known desiccation-resistant propagule; thus, long-distance dispersal events are likely rare. We predicted strong genetic structure among drainage basins and admixture among sites within basins. We predicted greater genetic diversity at lower latitude sites because they likely serve as refugia and the origin of northward, post-Pleistocene range expansion. We used 10 microsatellite loci to investigate genetic diversity from 311 gametophytes from 18 sites in five major drainage basins: South Atlantic Gulf, Mid-Atlantic, Ohio River, Great Lakes, and Northeast. Our data showed strong genetic partitioning among drainage basins and among sites within basins, yet no isolation by distance was detected. Genetic diversity varied widely among sites and was not strictly related to latitude as predicted. The results from B. gelatinosum provide strong support that each stream site contributes to the unique genetic variation within the species, potentially due to limited dispersal and the prevailing reproductive mode of intragametophytic selfing. Simulations of migration suggested post-Pleistocene dispersal from the Mid-Atlantic. Batrachospermum gelatinosum potentially persisted in refugia that were just south of the ice margins rather than in the southernmost part of its range. Research of other taxa with similar ranges could determine whether these results are generally applicable for freshwater red algae. Nevertheless, these results from B. gelatinosum add to the growing literature focused on the patterns and genetic consequences of post-Pleistocene range expansion by eastern North American biota.

Genetic structure of populations of Salvia ceratophylloides endemic to southern Calabria (southern Italy)

Assessing the degree of genetic diversity and differentiation of rare or endangered endemic species is essential to evaluate the conservation status of populations and successively implement appropriate conservation strategies. We investigated the population structure of Salvia ceratophylloides Ard., a scapose hemicryptophyte endemic to Calabria (southern Italy), both to answer questions about its genetic structure and to determine whether the actual population size has undergone significant demographic changes in the near past. The data obtained from the census showed that the populations are characterised by a greater number of adult individuals than juveniles and are on declining. The genetic analysis carried out on 99 individuals from four populations of the species under study, shows a mean expected heterozygosity value of 0.50 and an overall differentiation value of 0.083. The population structure shows that the four studied populations are distinct genetic units, genetically linked to four different ancestral gene pools. Bayesian analysis based on ABC models indicates that the present populations underwent a significant reduction in size in the past. This corresponds to the demographic decline at the end of the 19th century, which according to the literature, was due to the strong anthropic pressure (agriculture, grazing, fire and plantations) of Reggio Calabria suburbs. We can therefore conclude that populations are not affected by inbreeding and low genetic diversity and that there is no immediate danger of genetic erosion, and that the problems associated with population decline, past and present, are exclusively due to anthropogenic causes.

A Comparative Genetic Analysis of Phoenix atlantica in Cape Verde

The Cape Verde palm tree, Phoenix atlantica, holds significant ecological and cultural importance within the Cape Verde archipelago. However, its genetic distinctiveness has been questioned due to its close relationship and morphological similarity to the date palm (Phoenix dactylifera). In this study, we used an expanded sample set, 18 simple sequence repeat (SSR) markers, and a plastid minisatellite to characterize P. atlantica in Cape Verde and investigate its relationship with other Phoenix species. Our findings identify genetic markers that differentiate the P. atlantica genetic pool, including a unique fixed allele. We also provide evidence of the recent divergence of P. atlantica from Northern African date palm populations, suggesting a relatively recent colonization of Cape Verde by palm trees. Additionally, we characterized the genetic composition of palm tree populations across three Cape Verde islands, concluding that wild samples from certain populations in Boavista and Sal are best suited for establishing a seed and/or germplasm bank for replantation efforts, representing a crucial step for the conservation of Cape Verde's natural heritage. Overall, our results enhance the understanding of the historical trajectories and genetic characterization of palm trees in Africa, offering valuable insights for conservation strategies.

Across two phylogeographic breaks: Quaternary evolutionary history of a mountain aspen (Populus rotundifolia) in the Hengduan Mountains

Biogeographical barriers to gene flow are central to plant phylogeography. In East Asia, plant distribution is greatly influenced by two phylogeographic breaks, the Mekong-Salween Divide and Tanaka-Kaiyong Line, however, few studies have investigated how these barriers affect the genetic diversity of species that are distributed across both. Here we used 14 microsatellite loci and four chloroplast DNA fragments to examine genetic diversity and distribution patterns of 49 populations of Populus rotundifolia, a species that spans both the Mekong-Salween Divide and the Tanaka-Kaiyong Line in southwestern China. Demographic and migration hypotheses were tested using coalescent-based approaches. Limited historical gene flow was observed between the western and eastern groups of P. rotundifolia, but substantial flow occurred across both the Mekong-Salween Divide and Tanaka-Kaiyong Line, manifesting in clear admixture and high genetic diversity in the central group. Wind-borne pollen and seeds may have facilitated the dispersal of P. rotundifolia following prevalent northwest winds in the spring. We also found that the Hengduan Mountains, where multiple genetic barriers were detected, acted on the whole as a barrier between the western and eastern groups of P. rotundifolia. Ecological niche modeling suggested that P. rotundifolia has undergone range expansion since the last glacial maximum, and demographic reconstruction indicated an earlier population expansion around 600 Ka. The phylogeographic pattern of P. rotundifolia reflects the interplay of biological traits, wind patterns, barriers, niche differentiation, and Quaternary climate history. This study emphasizes the need for multiple lines of evidence in understanding the Quaternary evolution of plants in topographically complex areas.

Intraspecific divergence in a coastal plant, Euphorbia jolkinii, at a major biogeographic boundary in East Asia

PREMISE

Quaternary climatic fluctuations and long-distance seed dispersal across the sea are critical factors affecting the distribution of coastal plants, but the spatiotemporal nature of population expansion and distribution change of East Asian coastal plants during this period are rarely examined. To explore this process, we investigated the genome-wide phylogenetic patterns of Euphorbia jolkinii Boiss. (Euphorbiaceae), which grows widely on littoral areas of Japan, Korea, and Taiwan.

METHODS

We used plastome sequences and genome-wide single nucleotide polymorphisms in samples across the species range to reveal phylogeographic patterns and spatiotemporal distributional changes. We conducted ecological niche modeling for the present and the last glacial maximum (LGM).

RESULTS

Genetic differentiation was observed between the northern and southern populations of E. jolkinii, separated by the major biogeographic boundary, the Tokara Gap. These two groups of populations differentiated during the glacial period and subsequently intermingled in the intermorainic areas of the central Ryukyu Islands after the LGM. Ecological niche models suggested that the potential range of E. jolkinii was restricted to southern Kyushu; however, it was widespread in the southern Ryukyu Islands and Taiwan during the LGM.

CONCLUSIONS

This study provides evidence of genetic differentiation among coastal plant populations separated by the prominent biogeographical boundary. Although coastal plants are typically expected to maintain population connectivity through sea-drifted seed dispersal, our findings suggest that genetic differences may arise because of a combination of limited gene flow and changes in climate during the glacial period.

Historical fragmentation and stepping-stone gene flow led to population genetic differentiation in a coastal seabird

Understanding the forces that shape population genetic structure is fundamental both for understanding evolutionary trajectories and for conservation. Many factors can influence the geographic distribution of genetic variation, and the extent to which local populations differ can be especially difficult to predict in highly mobile organisms. For example, many species of seabirds are essentially panmictic, but some show strong structure. Pigeon Guillemots (Cepphus columba; Charadriiformes: Alcidae) breed in small colonies scattered along the North Pacific coastline and feed in shallow nearshore waters year-round. Given their distribution, gene flow is potentially lower and population genetic structure is stronger than in most other high-latitude Northern Hemisphere seabirds. We screened variation in the mitochondrial control region, four microsatellite loci, and two nuclear introns in 202 Pigeon Guillemots representing three of five subspecies. Mitochondrial sequences and nuclear loci both showed significant population differences, although structure was weaker for the nuclear loci. Genetic differentiation was correlated with geographic distance between sampling locations for both the mitochondrial and nuclear loci. Mitochondrial gene trees and demographic modeling both provided strong evidence for two refugial populations during the Pleistocene glaciations: one in the Aleutian Islands and one farther east and south. We conclude that historical fragmentation combined with a stepping-stone model of gene flow led to the relatively strong population differentiation in Pigeon Guillemots compared to other high-latitude Northern Hemisphere seabird species. Our study adds to growing evidence that Pleistocene glaciation events affected population genetic structure not only in terrestrial species but also in coastal marine animals.

Parallel evolution of morphological and genomic selfing syndromes accompany the breakdown of heterostyly

Evolutionary transitions from outcrossing to selfing in flowering plants have convergent morphological and genomic signatures and can involve parallel evolution within related lineages. Adaptive evolution of morphological traits is often assumed to evolve faster than nonadaptive features of the genomic selfing syndrome. We investigated phenotypic and genomic changes associated with transitions from distyly to homostyly in the Primula oreodoxa complex. We determined whether the transition to selfing occurred more than once and investigated stages in the evolution of morphological and genomic selfing syndromes using 22 floral traits and both nuclear and plastid genomic data from 25 populations. Two independent transitions were detected representing an earlier and a more recently derived selfing lineage. The older lineage exhibited classic features of the morphological and genomic selfing syndrome. Although features of both selfing syndromes were less developed in the younger selfing lineage, they exhibited parallel development with the older selfing lineage. This finding contrasts with the prediction that some genomic changes should lag behind adaptive changes to morphological traits. Our findings highlight the value of comparative studies on the timing and extent of transitions from outcrossing to selfing between related lineages for investigating the tempo of morphological and molecular evolution.

Nuclear genetic diversity of head lice sheds light on human dispersal around the world

The human louse, Pediculus humanus, is an obligate blood-sucking ectoparasite that has coevolved with humans for millennia. Given the intimate relationship between this parasite and the human host, the study of human lice has the potential to shed light on aspects of human evolution that are difficult to interpret using other biological evidence. In this study, we analyzed the genetic variation in 274 human lice from 25 geographic sites around the world by using nuclear microsatellite loci and female-inherited mitochondrial DNA sequences. Nuclear genetic diversity analysis revealed the presence of two distinct genetic clusters I and II, which are subdivided into subclusters: Ia-Ib and IIa-IIb, respectively. Among these samples, we observed the presence of the two most common louse mitochondrial haplogroups: A and B that were found in both nuclear Clusters I and II. Evidence of nuclear admixture was uncommon (12%) and was predominate in the New World potentially mirroring the history of colonization in the Americas. These findings were supported by novel DIYABC simulations that were built using both host and parasite data to define parameters and models suggesting that admixture between cI and cII was very recent. This pattern could also be the result of a reproductive barrier between these two nuclear genetic clusters. In addition to providing new evolutionary knowledge about this human parasite, our study could guide the development of new analyses in other host-parasite systems.

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.

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

PREMISE

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 hybridization can play an important role in island evolution. Here we use genotyping-by-sequencing (GBS) to understand the roles of hybridization, ecology, and geographic isolation in the radiation of Canary Island Descurainia (Brassicaceae).

METHODS

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 supermatrix and gene tree approaches and hybridization events were examined using D-statistics and Approximate Bayesian Computation. Climatic data were analyzed to examine the relationship between ecology and diversification.

RESULTS

Analysis of the supermatrix data set resulted in a fully resolved phylogeny. Species networks suggest a hybridization event has occurred for D. gilva, with these results being supported by Approximate Bayesian Computation analysis. Strong phylogenetic signals for temperature and precipitation indicate one major ecological shift within Canary Island Descurainia.

CONCLUSIONS

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, hybridization 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 incomplete lineage sorting and gene flow when studying groups prone to hybridization; patterns that might otherwise be obscured in species trees.

Geographic Genetic Structure of Alectoris chukar in Türkiye: Post-LGM-Induced Hybridization and Human-Mediated Contaminations

Türkiye is considered an important evolutionary area for Chukar partridge (Alectoris chukar), since it is both a potential ancestral area and a diversification center for the species. Using 2 mitochondrial (Cty-b and D-loop) and 13 polymorphic microsatellite markers, we investigated the geographic genetic structure of A. chukar populations to determine how past climatic fluctuations and human activities have shaped the gene pool of this species in Türkiye. Our results indicate, firstly, that only A. chukar of the genus Alectoris is present in Türkiye (Anatolia and Thrace), with no natural or artificial gene flow from congenerics. Secondly, the geographic genetic structure of the species in Türkiye has been shaped by topographic heterogeneity, Pleistocene climatic fluctuations, and artificial transport by humans. Third, there appears to be three genetic clusters: Thracian, Eastern, and Western. Fourth, the post-LGM demographic expansion of the Eastern and Western populations has formed a hybrid zone in Central Anatolia (~8 kyBP). Fifth, the rate of China clade-B contamination in Türkiye is about 8% in mtDNA and about 12% in nuDNA, with the Southeastern Anatolian population having the highest contamination. Sixth, the Thracian population was the most genetically distinct, with the lowest genetic diversity and highest level of inbreeding and no China clad-B contamination. These results can contribute to the conservation regarding A. chukar populations, especially the Thracian population.

Dermacentor reticulatus - a tick on its way from glacial refugia to a panmictic Eurasian population

The ornate dog tick (Dermacentor reticulatus) shows a recently expanding geographic distribution. Knowledge on its intraspecific variability, population structure, rate of genetic diversity and divergence, including its evolution and geographic distribution, is crucial to understand its dispersal capacity. All such information would help to evaluate the potential risk of future spread of associated pathogens of medical and veterinary concern. A set of 865 D. reticulatus ticks was collected from 65 localities across 21 countries, from Portugal in the west to Kazakhstan and southern Russia in the east. Cluster analyses of 16 microsatellite loci were combined with nuclear (ITS2, 18S) and mitochondrial (12S, 16S, COI) sequence data to uncover the ticks' population structures and geographical patterns. Approximate Bayesian computation was applied to model evolutionary relationships among the found clusters. Low variability and a weak phylogenetic signal showing an east-west cline were detected both for mitochondrial and nuclear sequence markers. Microsatellite analyses revealed three genetic clusters, where the eastern and western cluster gradient was supplemented by a third, northern cluster. Alternative scenarios could explain such a tripartite population structure by independent formation of clusters in separate refugia, limited gene flow connected with isolation by distance causing a "bipolar pattern", and the northern cluster deriving from admixture between the eastern and western populations. The best supported demographic scenario of this tick species indicates that the northern cluster derived from admixture between the eastern and western populations 441 (median) to 224 (mode) generations ago, suggesting a possible link with the end of the Little Ice Age in Europe.

A wide hybrid zone mediated by precipitation contributed to confused geographical structure of Scutiger boulengeri

Confused geographical structure of a population and mitonuclear discordance are shaped by a combination of rapid changes in population demographics and shifts in ecology. In this study, we generated a time-calibrated phylogeny of Scutiger boulengeri, an endemic Xizang alpine toad occurring in mountain streams on the Qinghai-Xizang (Tibet) Plateau (QTP). Based on three mitochondrial DNA (mtDNA) genes, eight clades were assigned to three deeply divergent lineages. Analysis of nuclear DNA (nuDNA) genes revealed three distinct clusters without geographic structure, indicating significantly high rates of gene flow. Coalescent theory framework analysis (approximate Bayesian computation model DIYABC and Migrate-N) suggested that divergence of the main intraspecific clusters was the result of hybridization after secondary contact in the Holocene around 0.59 million years ago (Ma). The ratio of mtDNA F _ST (fixation index) to nuDNA F _ST was 2.3, thus failing to show male-biased dispersal. Geographic cline analysis showed that a wide hybrid zone was initially established in southwestern China, without significant reproductive isolation but with strong introgression in S. boulengeri, suggesting high hybrid fitness. Furthermore, mtDNA genes exhibited isolation by distance (IBD) while nuDNA genes exhibited significant isolation by environment (IBE). Results suggested that mitonuclear discordance may have initially been caused by geographic isolation, followed by precipitation-mediated hybridization, producing a wide hybrid zone and geographic structure confusion of nuDNA genes in S. boulengeri. This study indicated that complicated historical processes may have led to specific genetic patterns, with a specific climate factor facilitating gene flow in the system.

Genetic variation of endangered Jankowski’s Bunting (Emberiza jankowskii): High connectivity and a moderate history of demographic decline

Introduction

Continued discovery of “mismatch” patterns between population size and genetic diversity, involving wild species such as insects, amphibians, birds, mammals, and others, has raised issues about how population history, especially recent dynamics under human disturbance, affects currently standing genetic variation. Previous studies have revealed high genetic diversity in endangered Jankowski’s Bunting. However, it is unclear how the demographic history and recent habitat changes shape the genetic variation of Jankowski’s Bunting.

Methods

To explore the formation and maintenance of high genetic diversity in endangered Jankowski’s Bunting, we used a mitochondrial control region (partial mtDNA CR) and 15 nuclear microsatellite markers to explore the recent demographic history of Jankowski’s Bunting, and we compared the historical and contemporary gene flows between populations to reveal the impact of habitat change on population connectivity. Specifically, we aimed to test the following hypotheses: (1) Jankowski’s Bunting has a large historical Ne and a moderate demographic history; and (2) recent habitat change might have no significant impact on the species’ population connectivity.

Results

The results suggested that large historical effective population size, as well as severe but slow population decline, may partially explain the high observable genetic diversity. Comparison of historical (over the past 4Ne generations) and contemporary (1–3 generations) gene flow indicated that the connectivity between five local populations was only marginally affected by landscape changes.

Discussion

Our results suggest that high population connectivity and a moderate history of demographic decline are powerful explanations for the rich genetic variation in Jankowski’s Bunting. Although there is no evidence that the genetic health of Jankowski’s Bunting is threatened, the time-lag effects on the genetic response to recent environmental changes is a reminder to be cautious about the current genetic characteristics of this species. Where possible, factors influencing genetic variation should be integrated into a systematic framework for conducting robust population health assessments. Given the small contemporary population size, inbreeding, and ecological specialization, we recommend that habitat protection be maintained to maximize the genetic diversity and population connectivity of Jankowski’s Bunting.

Transcriptome-based analyses of adaptive divergence between two closely related spruce species on the Qinghai-Tibet plateau and adjacent regions

Speciation among populations connected by gene flow is driven by adaptation to different environments, but underlying gene-environment associations remain largely unknown. Here, 162 individuals from 32 populations were sampled to obtain 191,648 independent single nucleotide polymorphisms (SNPs) across the genomes of two closely related spruce species, Picea asperata and Picea crassifolia, which occur on the Qinghai-Tibet Plateau and in surrounding regions. Using the SNP data set, genotype-environment associations and demographic modelling were used to examine local adaptation and genetic divergence between these two species. While morphologically similar, the two Picea species were genetically differentiated in multiple analyses. These species diverged despite continuous gene flow, and their initial divergence was dated back to the late Quaternary. The effective population sizes of both species have expanded since their divergence, as confirmed by niche distribution simulations. A total of 6365 genes were associated with the tested environmental variables; of these, 41 were positively selected in P. asperata and were mainly associated with temperature, while 83 were positively selected in P. crassifolia and were primarily associated with precipitation. These results deepen our understanding of the adaptive divergence and demographic histories of these two spruce species and highlight the importance of genomic data in deciphering the environmental selection underlying Quaternary interspecific divergence.

Genetic and epigenetic variation separately contribute to range expansion and local metalliferous habitat adaptation during invasions of Chenopodium ambrosioides into China

BACKGROUND AND AIMS

Invasive plants often colonize wide-ranging geographical areas with various local microenvironments. The specific roles of epigenetic and genetic variation during such expansion are still unclear. Chenopodium ambrosioides is a well-known invasive alien species in China that can thrive in metalliferous habitats. This study aims to comprehensively understand the effects of genetic and epigenetic variation on the successful invasion of C. ambrosioides.

METHODS

We sampled 367 individuals from 21 heavy metal-contaminated and uncontaminated sites with a wide geographical distribution in regions of China. We obtained environmental factors of these sampling sites, including 13 meteorological factors and the contents of four heavy metals in soils. Microsatellite markers were used to investigate the demographic history of C. ambrosioides populations in China. We also analysed the effect of epigenetic variation on metalliferous microhabitat adaptation using methylation-sensitive amplified polymorphism (MSAP) markers. A common garden experiment was conducted to compare heritable phenotypic variations among populations.

KEY RESULTS

Two distinct genetic clusters that diverged thousands of years ago were identified, suggesting that the eastern and south-western C. ambrosioides populations in China may have originated from independent introduction events without recombination. Genetic variation was shown to be a dominant determinant of phenotypic differentiation relative to epigenetic variation, and further affected the geographical distribution pattern of invasive C. ambrosioides. The global DNA unmethylation level was reduced in metalliferous habitats. Dozens of methylated loci were significantly associated with the heavy metal accumulation trait of C. ambrosioides and may contribute to coping with metalliferous microenvironments.

CONCLUSIONS

Our study of C. ambrosioides highlighted the dominant roles of genetic variation in large geographical range expansion and epigenetic variation in local metalliferous habitat adaptation.

Possible northern persistence of Siebold's beech, Fagus crenata, at its northernmost distribution limit on an island in Japan Sea: Okushiri Island, Hokkaido

Siebold's beech, Fagus crenata, is widely distributed across the Japanese Archipelago and islands in Japan Sea. Similar to the northern limit of the geographical distribution of F. crenata on the mainland of Hokkaido, the northern limit of the distribution of F. crenata on islands in the Japan Sea is observed on Okushiri Island (ca 42°N). To understand the genetic relationships of F. crenata on Okushiri Island, we examined chloroplast (cp) DNA haplotypes and 11 nuclear microsatellite (SSR) loci among 1,838 individuals from 44 populations from Okushiri Island, mainland Hokkaido, and the northern part of the Tohoku region on Honshu Island. We identified 2 cpDNA haplotypes, which represent not only populations on the Japan Sea coast but also those on the Pacific coast and this suggested the Okushiri Island populations might not be formed by single colonization. Genetic diversity of the Okushiri Island populations of nuclear SSR was not lower than the mainland and the STRUCTURE analysis revealed the Okushiri Island individuals were admixed between Hokkaido and Tohoku clusters. Approximate Bayesian computation inferred that divergence between Tohoku and Hokkaido, and admixture between two populations which generated Okushiri populations occurred before the last glacial maximum (LGM), that is, 7,890 (95% hyper probability density (HPD): 3,420 - 9,910) and 3,870 (95% HPD: 431- 8,540) generations ago, respectively. These inferences were well supported by a geological history which suggested an isolation of Okushiri Island from Hokkaido started prior to the Middle Pleistocene. We discuss the possible persistence of F. crenata during the last glacial maximum on northern islands in the Japan Sea such as Okushiri Island.

Survivors and colonizers: Contrasting biogeographic histories reconciled in the Antarctic freshwater copepod Boeckella poppei

Two main hypotheses have been proposed to explain the contemporary distribution of Antarctic terrestrial biota. We assess whether the current distribution of maritime Antarctic populations of the freshwater copepod Boeckella poppei is the result of (1) a post-Last Glacial Maximum (LGM) colonization, or whether (2) the species survived in regional glacial refugia throughout the LGM and earlier glaciations. Using 438 specimens from 34 different sampling sites across Southern South America, South Georgia, South Orkney Islands, South Shetland Islands, and the Antarctic Peninsula, we analyzed mitochondrial and nuclear sequences to uncover patterns of genetic diversity and population structure. We also performed median-joining haplotype network, phylogenetic reconstruction, and divergence time analyses. Finally, we evaluated past demographic changes and historical scenarios using the Approximate Bayesian Computation (ABC) method. Our data support the existence of two clades with different and contrasting biogeographic histories. The first clade has been present in maritime Antarctica since at least the mid-Pleistocene, with the South Orkney Islands the most likely refugial area. The second clade has a broader distribution including southern South America, South Georgia, South Shetland Islands, and the Antarctic Peninsula. The ABC method identified long-distance dispersal (LDD) colonization event(s) from southern South America to South Georgia and the maritime Antarctic after the LGM deglaciation, supporting more recent colonization of Antarctic locations. The current Antarctic and sub-Antarctic distribution of B. poppei is likely derived from two independent biogeographic events. The combination of both (1) post-LGM colonization from southern South America and (2) longer-term persistence in in situ regional refugia throughout glacial periods challenges current understanding of the biogeographic history of Antarctic freshwater biota. Re-colonization of ice-impacted Antarctic areas would have occurred following a LDD and Establishment model, pointing to the existence of possible post-dispersal barriers, despite widely assumed high passive dispersal capacity in freshwater invertebrates.

Convolutional Neural Networks as Summary Statistics for Approximate Bayesian Computation

Approximate Bayesian Computation is widely used in systems biology for inferring parameters in stochastic gene regulatory network models. Its performance hinges critically on the ability to summarize high-dimensional system responses such as time series into a few informative, low-dimensional summary statistics. The quality of those statistics acutely impacts the accuracy of the inference task. Existing methods to select the best subset out of a pool of candidate statistics do not scale well with large pools of several tens to hundreds of candidate statistics. Since high quality statistics are imperative for good performance, this becomes a serious bottleneck when performing inference on complex and high-dimensional problems. This paper proposes a convolutional neural network architecture for automatically learning informative summary statistics of temporal responses. We show that the proposed network can effectively circumvent the statistics selection problem of the preprocessing step for ABC inference. The proposed approach is demonstrated on two benchmark problem and one challenging inference problem learning parameters in a high-dimensional stochastic genetic oscillator. We also study the impact of experimental design on network performance by comparing different data richness and data acquisition strategies.

Population Genetics of the Blueberry Gall Midge, Dasineura oxycoccana (Diptera: Cecidomyiidae), on Blueberry and Cranberry and Testing Invasion Scenarios

We compared the population genetic structure between populations of the blueberry gall midge-Dasineura oxycoccana (Johnson) (Diptera: Cecidomyiidae)-from blueberry and cranberry and determined the genetic relationships among geographical subgroups by genotyping 632 individuals from 31 different populations from their native USA regions (New Jersey, Michigan, and Georgia) and from invaded Korean regions using 12 microsatellite loci. Our population genetic analyses showed a clear separation between the two host-associated D. oxycoccana populations from blueberry and cranberry. Using data from only the blueberry-associated D. oxycoccana populations, we identified five genetically isolated subgroups. An analysis of the approximate Bayesian computation suggests that the invasive D. oxycoccana population from Korea appears to have been introduced from an unsampled source population rather than directly from its native range. Our findings will allow for an easier identification of the source of D. oxycoccana into newly invaded regions, as well as to determine their association with blueberry and cranberry, which based on our results can be considered as two distinct species.

Response to: A paleoecological context to assess the development of oak forest in Colombia: A comment on Zorrilla-Azcué, S., González-Rodríguez, A., Oyama, K., González, M.A. & Rodríguez-Correa, H., The DNA history of a lonely oak: Quercus humboldtii phylogeography in the Colombian Andes. Ecology and Evolution 2021, doi:10.1002/ece3.7529

In this response, we address comments and clarify the rationale behind the choice of hypotheses aimed to describe the Quercus humboldtii phylogeography in the Colombian Andes. Finally, we explain our disagreement with the conclusions of a previous critique, since these are not necessarily adequate under the implemented population genetics approach.

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.

ddRAD-Seq reveals evolutionary insights into population differentiation and the cryptic phylogeography of Hyporhamphus intermedius in Mainland China

Species differentiation and local adaptation in heterogeneous environments have attracted much attention, although little is known about the mechanisms involved. Hyporhamphus intermedius is an anadromous, brackish-water halfbeak that is widely distributed in coastal areas and hyperdiverse freshwater systems in China, making it an interesting model for research on phylogeography and local adaptation. Here, 156 individuals were sampled at eight sites from heterogeneous aquatic habitats to examine environmental and genetic contributions to phenotypic divergence. Using double-digest restriction-site-associated DNA sequencing (ddRAD-Seq) in the specimens from the different watersheds, 5498 single nucleotide polymorphisms (SNPs) were found among populations, with obvious population differentiation. We find that present-day Mainland China populations are structured into distinct genetic clusters stretching from southern and northern ancestries, mirroring geography. Following a transplant event in Plateau Lakes, there were virtually no variations of genetic diversity occurred in two populations, despite the fact two main splits were unveiled in the demographic history. Additionally, dorsal, and anal fin traits varied widely between the southern group and the others, which highlighted previously unrecognized lineages. We then explore genotype-phenotype-environment associations and predict candidate loci. Subgroup ranges appeared to correspond to geographic regions with heterogeneous hydrological factors, indicating that these features are likely important drivers of diversification. Accordingly, we conclude that genetic and phenotypic polymorphism and a moderate amount of genetic differentiation occurred, which might be ascribed to population subdivision, and the impact of abiotic factors.

Assessing the Genetic Identity of Tuscan Sweet Chestnut (Castanea sativa Mill.)

The sweet chestnut (Castanea sativa Mill.) is an important species of European trees, studied for both ecological and economic reasons. Its cultivation in the Italian peninsula can be linked to the Roman period and has been documented, especially in the Tuscan region, for centuries. We sampled 131 grafted trees from three separate areas to determine the genetic variability between populations and assess genetic identity for different varieties of trees, which is useful for future breeding programs and propagation efforts. Molecular analyses were performed using eight microsatellite loci. A total of 98 alleles was detected with an average of 12.3 alleles per locus. We found high levels of genetic diversity within the varieties of the same area, ranging between He = 0.682–0.745. Of the eight loci, seven were found to be at Hardy-Weinberg equilibrium. (FST values Differentiation between cultivation areas was significant between 0.052–0.147) with the two Southern Tuscan areas showing the closest relationship as also indicated by Bayesian inference of the population structure, which revealed the existence of three ancestral gene pools of origin. Demographic events were detected by a coalescent-based approximate Bayesian computation in two of the homogeneous clusters. This work is a step forward for the conservation of this iconic species, albeit at a regional level, as chestnut varieties have never received the full attention of breeders.

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.

riceExplorer: Uncovering the Hidden Potential of a National Genomic Resource Against a Global Database

Agricultural crop breeding programs, particularly at the national level, typically consist of a core panel of elite breeding cultivars alongside a number of local landrace varieties (or other endemic cultivars) that provide additional sources of phenotypic and genomic variation or contribute as experimental materials (e.g., in GWAS studies). Three issues commonly arise. First, focusing primarily on core development accessions may mean that the potential contributions of landraces or other secondary accessions may be overlooked. Second, elite cultivars may accumulate deleterious alleles away from nontarget loci due to the strong effects of artificial selection. Finally, a tendency to focus solely on SNP-based methods may cause incomplete or erroneous identification of functional variants. In practice, integration of local breeding programs with findings from global database projects may be challenging. First, local GWAS experiments may only indicate useful functional variants according to the diversity of the experimental panel, while other potentially useful loci-identifiable at a global level-may remain undiscovered. Second, large-scale experiments such as GWAS may prove prohibitively costly or logistically challenging for some agencies. Here, we present a fully automated bioinformatics pipeline (riceExplorer) that can easily integrate local breeding program sequence data with international database resources, without relying on any phenotypic experimental procedure. It identifies associated functional haplotypes that may prove more robust in determining the genotypic determinants of desirable crop phenotypes. In brief, riceExplorer evaluates a global crop database (IRRI 3000 Rice Genomes) to identify haplotypes that are associated with extreme phenotypic variation at the global level and recorded in the database. It then examines which potentially useful variants are present in the local crop panel, before distinguishing between those that are already incorporated into the elite breeding accessions and those only found among secondary varieties (e.g., landraces). Results highlight the effectiveness of our pipeline, identifying potentially useful functional haplotypes across the genome that are absent from elite cultivars and found among landraces and other secondary varieties in our breeding program. riceExplorer can automatically conduct a full genome analysis and produces annotated graphical output of chromosomal maps, potential global diversity sources, and summary tables.

Revisiting the comparative phylogeography of unglaciated eastern North America: 15 years of patterns and progress

In a landmark comparative phylogeographic study, “Comparative phylogeography of unglaciated eastern North America,” Soltis et al. (Molecular Ecology, 2006, 15, 4261) identified geographic discontinuities in genetic variation shared across taxa occupying unglaciated eastern North America and proposed several common biogeographical discontinuities related to past climate fluctuations and geographic barriers. Since 2006, researchers have published many phylogeographical studies and achieved many advances in genotyping and analytical techniques; however, it is unknown how this work has changed our understanding of the factors shaping the phylogeography of eastern North American taxa. We analyzed 184 phylogeographical studies of eastern North American taxa published between 2007 and 2019 to evaluate: (1) the taxonomic focus of studies and whether a previously detected taxonomic bias towards studies focused on vertebrates has changed over time, (2) the extent to which studies have adopted genotyping technologies that improve the resolution of genetic groups (i.e., NGS DNA sequencing) and analytical approaches that facilitate hypothesis‐testing (i.e., divergence time estimation and niche modeling), and (3) whether new studies support the hypothesized biogeographic discontinuities proposed by Soltis et al. (Molecular Ecology, 2006, 15, 4261) or instead support new, previously undetected discontinuities. We observed little change in taxonomic focus over time, with studies still biased toward vertebrates. Although many technological and analytical advances became available during the period, uptake was slow and they were employed in only a small proportion of studies. We found variable support for previously identified discontinuities and identified one new recurrent discontinuity. However, the limited resolution and taxonomic breadth of many studies hindered our ability to clarify the most important climatological or geographical factors affecting taxa in the region. Broadening the taxonomic focus to include more non‐vertebrate taxa, employing technologies that improve genetic resolution, and using analytical approaches that improve hypothesis testing are necessary to strengthen our inference of the forces shaping the phylogeography of eastern North America.

Efficient ancestry and mutation simulation with msprime 1.0

Stochastic simulation is a key tool in population genetics, since the models involved are often analytically intractable and simulation is usually the only way of obtaining ground-truth data to evaluate inferences. Because of this, a large number of specialized simulation programs have been developed, each filling a particular niche, but with largely overlapping functionality and a substantial duplication of effort. Here, we introduce msprime version 1.0, which efficiently implements ancestry and mutation simulations based on the succinct tree sequence data structure and the tskit library. We summarize msprime's many features, and show that its performance is excellent, often many times faster and more memory efficient than specialized alternatives. These high-performance features have been thoroughly tested and validated, and built using a collaborative, open source development model, which reduces duplication of effort and promotes software quality via community engagement.

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.

Origin, Persistence, and Vulnerability to Climate Changes of Podocarpus Populations in Central African Mountains

Background and objectives—Podocarpus latifolius (synonym of P. milanjianus) is a key tree representative of Afromontane forests where it is highly threatened by climate and land-use changes. While large populations occur in East Africa, only a few isolated and usually small populations remain in western Central Africa (Cameroon to Angola). Studying the evolutionary history of such relictual populations can thus be relevant to understand their resilience under changing environments. Materials and Methods—we developed nine polymorphic nuclear microsatellites (nSSRs) to estimate genetic variability, (historical) gene flow, and demographic changes among natural populations from Central to East Africa. Results—despite the extended distribution range of P. latifolius, a strong isolation-by-distance pattern emerges at the intra-population scale, indicating low seed and pollen dispersal capacities. Central African populations display a lower genetic diversity (He = 0.34 to 0.61) and are more differentiated from each other (FST = 0.28) than are East African populations (He = 0.65 to 0.71; FST = 0.10), suggesting high genetic drift in the Central African populations. Spatial genetic structure reveals past connections between East and West Africa but also a gene flow barrier across the equator in western Central Africa. Demographic modelling anchors the history of current lineages in the Pleistocene and supports a strong demographic decline in most western populations during the last glacial period. By contrast, no signature of demographic change was detected in East African populations. Conclusions—in Cameroon, our results exclude a recent (re)colonization from one source population of all mountain ranges, but rather indicate long-term persistence of populations in each mountain with fluctuating sizes. A higher impact of genetic drift and further loss of diversity can be expected by survival through climatically unfavorable periods in such small refugial populations. Tracking the Quaternary legacy of podocarp populations is thus essential for their conservation since there is a temporal gap between environment crises and an ecological/genetic answer at the population level.

A brief history and popularity of methods and tools used to estimate micro-evolutionary forces

Population genetics is a field of research that predates the current generations of sequencing technology. Those approaches, that were established before massively parallel sequencing methods, have been adapted to these new marker systems (in some cases involving the development of new methods) that allow genome-wide estimates of the four major micro-evolutionary forces-mutation, gene flow, genetic drift, and selection. Nevertheless, classic population genetic markers are still commonly used and a plethora of analysis methods and programs is available for these and high-throughput sequencing (HTS) data. These methods employ various and diverse theoretical and statistical frameworks, to varying degrees of success, to estimate similar evolutionary parameters making it difficult to get a concise overview across the available approaches. Presently, reviews on this topic generally focus on a particular class of methods to estimate one or two evolutionary parameters. Here, we provide a brief history of methods and a comprehensive list of available programs for estimating micro-evolutionary forces. We furthermore analyzed their usage within the research community based on popularity (citation bias) and discuss the implications of this bias for the software community. We found that a few programs received the majority of citations, with program success being independent of both the parameters estimated and the computing platform. The only deviation from a model of exponential growth in the number of citations was found for the presence of a graphical user interface (GUI). Interestingly, no relationship was found for the impact factor of the journals, when the tools were published, suggesting accessibility might be more important than visibility.

Influence of Pleistocene climatic oscillations on the phylogeography and demographic history of endemic vulnerable trees (section Magnolia) of the Tropical Montane Cloud Forest in Mexico

The Tropical Montane Cloud Forest (TMCF) is a highly dynamic ecosystem that has undergone frequent spatial changes in response to the interglacial-glacial cycles of the Pleistocene. These climatic fluctuations between cold and warm cycles have led to species range shifts and contractions-expansions, resulting in complex patterns of genetic structure and lineage divergence in forest tree species. In this study, we sequenced four regions of the chloroplast DNA (trnT-trnL, trnK5-matk, rpl32-trnL, trnS-trnG) for 20 populations and 96 individuals to evaluate the phylogeography, historical demography, and paleodistributions of vulnerable endemic TMCF trees in Mexico: Magnolia pedrazae (north-region), M. schiedeana (central-region), and M. schiedeana population Oaxaca (south-region). Our data recovered 49 haplotypes that showed a significant phylogeographic structure in three regions: north, central, and south. Bayesian Phylogeographic and Ecological Clustering (BPEC) analysis also supported the divergence in three lineages and highlighted the role of environmental factors (temperature and precipitation) in genetic differentiation. Our historical demography analyses revealed demographic expansions predating the Last Interglacial (LIG, ~125,000 years ago), while Approximate Bayesian Computation (ABC) simulations equally supported two contrasting demographic scenarios. The BPEC and haplotype network analyses suggested that ancestral haplotypes were geographically found in central Veracruz. Our paleodistributions modeling showed evidence of range shifts and expansions-contractions from the LIG to the present, which suggested the complex evolutionary dynamics associated to the climatic oscillations of the Pleistocene. Habitat management of remnant forest fragments where large and genetically diverse populations occur in the three TMCF regions analyzed would be key for the conservation of these magnolia populations.

Population Genetics for Inferring Introduction Sources of the Oriental Fruit Fly, Bactrocera dorsalis: A Test for Quarantine Use in Korea

To infer the introduction sources of the oriental fruit fly, Bactrocera dorsalis, we used a mitochondrial marker to reconstruct the haplotype network and 15 microsatellite loci to reveal genetic structure and relationships between the geographically or temporally different collections from Asia. We performed Approximate Bayesian computations to infer a global origin and a source of the quarantine collections found in Korea. As a result, the 40 populations were divided into three groups, of which genetic similarity is not related to the geographic vicinity. Korean samples had a similar genetic structure to Taiwan and Thailand ones. Our results suggest that the place of origin of the B. dorsalis specimens found in Korea's border quarantine is likely to be Taiwan or Thailand. As the global origin of B. dorsalis, we estimated that Taiwan and Thailand were most likely the global origins of Southeast Asian populations by testing hypothetical scenarios by the approximate Bayesian computation analyses. Our results will allow easier identification of the source region of the forthcoming invasion of quarantined B. dorsalis specimens.

Signature of mid-Pleistocene lineages in the European silver fir (Abies alba Mill.) at its geographic distribution margin

In a conservation and sustainable management perspective, we identify the ecological, climatic, and demographic factors responsible for the genetic diversity patterns of the European silver fir (Abies alba Mill.) at its southwestern range margin (Pyrenees Mountains, France, Europe). We sampled 45 populations throughout the French Pyrenees and eight neighboring reference populations in the Massif Central, Alps, and Corsica. We genotyped 1,620 individuals at three chloroplast and ten nuclear microsatellite loci. We analyzed within- and among-population genetic diversity using phylogeographic reconstructions, tests of isolation-by-distance, Bayesian population structure inference, modeling of demographic scenarios, and regression analyses of genetic variables with current and past environmental variables. Genetic diversity decreased from east to west suggesting isolation-by-distance from the Alps to the Pyrenees and from the Eastern to the Western Pyrenees. We identified two Pyrenean lineages that diverged from a third Alpine-Corsica-Massif Central lineage 0.8 to 1.1 M years ago and subsequently formed a secondary contact zone in the Central Pyrenees. Population sizes underwent contrasted changes, with a contraction in the west and an expansion in the east. Glacial climate affected the genetic composition of the populations, with the western genetic cluster only observed in locations corresponding to the coldest past climate and highest elevations. The eastern cluster was observed over a larger range of temperatures and elevations. All demographic events shaping the current spatial structure of genetic diversity took place during the Mid-Pleistocene Transition, long before the onset of the Holocene. The Western Pyrenees lineage may require additional conservation efforts, whereas the eastern lineage is well protected in in situ gene conservation units. Due to past climate oscillations and the likely emergence of independent refugia, east-west oriented mountain ranges may be important reservoir of genetic diversity in a context of past and ongoing climate change in Europe.

Four times out of Europe: Serial invasions of the winter moth, Operophtera brumata, to North America

Reconstructing the geographic origins of non-native species is important for studying the factors that influence invasion success, however; these analyses can be constrained by the amount of diversity present in the native and invaded regions, and by changes in the genetic background of the invading population following bottlenecks and/or hybridization events. Here we explore the geographical origins of the invasive winter moth (Operopthera brumata L.) that has caused widespread defoliation to forests, orchards, and crops in Nova Scotia, British Columbia, Oregon, and the northeastern United States. It is not known whether these represent independent introductions to North America, or a "stepping stone" spread among regions. Using a combination of Bayesian assignment and approximate Bayesian computation methods, we analysed a population genetic data set of 24 microsatellite loci. We estimate that winter moth was introduced to North America on at least four occasions, with the Nova Scotian and British Columbian populations probably being introduced from France and Sweden, respectively; the Oregonian population probably being introduced from either the British Isles or northern Fennoscandia; and the population in the northeastern United States probably being introduced from somewhere in Central Europe. We discuss the impact of genetic bottlenecks on analyses meant to determine region of origin.

Tracing the Origin of Korean Invasive Populations of the Spotted Lanternfly, Lycorma delicatula (Hemiptera: Fulgoridae)

Lycorma delicatula (White) suddenly arrived in Korea where it rapidly spread out in the central region of Korea and caused serious damage to grape vineyards. To trace the source region of its invasiveness, population genetic structures were compared between the native region, China, and the introduced regions, Korea and Japan. We examined 762 individuals from 38 different population collections using 15 microsatellite loci. Both principal coordinate and structure analyses displayed that the Chinese populations were separated into three subgroups which were located significantly far apart from each other. Among them, the Shanghai population was located closest to most Korean populations. Based on the genetic relationships and structures, it was revealed that the multiple introductions into Korea occurred at least three times. In addition, the Shanghai population was strongly estimated to be a source of initial invasive populations of Korea. In addition, analysis of the approximate Bayesian computation suggested simultaneous spread from two distant locations early in the invasion by artificial transportation of the host plants bearing egg masses. Our population genetics study can provide a precedent case with regards to identifying spreads by anthropogenic outcomes in other invasive regions.

Dispersal and mating patterns determine the fate of naturally dispersed populations: evidence from Bombina orientalis

Background

In contrast to the explosive increase of a population following biological invasion, natural dispersal, i.e., when a population disperses from its original range into a new range, is a passive process that is affected by resources, the environment, and other factors. Natural dispersal is also negatively impacted by genetic drift and the founder effect. Although the fates of naturally dispersed populations are unknown, they can adapt evolutionarily over time to the new environment. Can naturally dispersed populations evolve beneficial adaptive strategies to offset these negative effects to maintain their population in a stable state?

Results

The current study addressed this question by focusing on the toad Bombina orientalis, the population of which underwent natural dispersal following the Last Glacial Maximum in Northeast Asia. Population genetic approaches were used to determine the genetic structure, dispersal pattern, and mating system of the population of B. orientalis in northeast China (Northern population). The results showed that this northern population of B. orientalis is a typical naturally dispersed population, in which the stable genetic structure and high level of genetic diversity of the population have been maintained through the long-distance biased dispersal behavior of males and the pattern of promiscuity within the population.

Conclusions

Our findings suggest that naturally dispersed populations can evolve effective adaptive strategies to maintain a stable population. Different species may have different strategies. The relevance of these maintenance mechanisms for naturally dispersed populations provide a new perspective for further understanding the processes of speciation and evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01844-3.