DIYABC v2.0: a software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism, DNA sequence and microsatellite data

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.

Coalescent models derived from birth-death processes

A coalescent model of a sample of size n is derived from a birth-death process that originates at a random time in the past from a single founder individual. Over time, the descendants of the founder evolve into a population of large (infinite) size from which a sample of size n is taken. The parameters and time of the birth-death process are scaled in N₀, the size of the present-day population, while letting N₀→∞, similarly to how the standard Kingman coalescent process arises from the Wright-Fisher model. The model is named the Limit Birth-Death (LBD) coalescent model. Simulations from the LBD coalescent model with sample size n are computationally slow compared to standard coalescent models. Therefore, we suggest different approximations to the LBD coalescent model assuming the population size is a deterministic function of time rather than a stochastic process. Furthermore, we introduce a hybrid LBD coalescent model, that combines the exactness of the LBD coalescent model model with the speed of the approximations.

Complete Chloroplast Genome Variants Reveal Discrete Long-Distance Dispersal Routes of Rhizophora in the Western Indian Ocean

Historical processes of long-distance migration and ocean-wide expansion feature the global biogeographic pattern of Rhizophora species. Throughout the Indian Ocean, Rhizophora stylosa and Rhizophora mucronata seem to be a young phylogenetic group with an expansion of R. mucronata toward the Western Indian Ocean (WIO) driven by the South Equatorial Current (SEC). Nuclear microsatellites revealed genetic patterns and breaks; however, the estimation of propagule dispersal routes requires maternally inherited cytoplasmic markers. Here, we examine the phylogeography of 21 R. mucronata provenances across a >4,200 km coastal stretch in the WIO using R. stylosa as an outgroup. Full-length chloroplast genome (164,474 bp) and nuclear ribosomal RNA cistron (8,033 bp) sequences were assembled. The boundaries, junction point, sequence orientation, and stretch between LSC/IRb/SSC/IRa/LSC showed no differences with R. stylosa chloroplast genome. A total of 58 mutations in R. mucronata encompassing transitions/transversions, insertions-deletions, and mononucleotide repeats revealed three major haplogroups. Haplonetwork, Bayesian maximum likelihood (ML), and approximate Bayesian computation (ABC) analyses supported discrete historical migration events. An ancient haplogroup A in the Seychelles and eastern Madagascar was as different from other haplogroups as from R. stylosa. A star-like haplonetwork referred as the recent range expansion of haplogroup B from northern Madagascar toward the African mainland coastline, including a single variant spanning >1,800 km across the Mozambique Channel area (MCA). Populations in the south of Delagoa Bight contained haplogroup C and was originated from a unique bottleneck dispersal event. Divergence estimates of pre- and post-Last Glacial Maximum (LGM) illustrated the recent emergence of Rhizophora mangroves in the WIO compared to other oceans. Connectivity patterns could be aligned with the directionality of major ocean currents. Madagascar and the Seychelles each harbored haplogroups A and B, albeit among spatially separated populations, explained from a different migration era. Likewise, the Aldabra Atoll harbored spatially distinct haplotypes. Nuclear ribosomal cistron (8,033 bp) variants corresponded to haplogroups and confirmed admixtures in the Seychelles and Aldabra. These findings shed new light on the origins and dispersal routes of R. mucronata lineages that have shaped their contemporary populations in large regions of the WIO, which may be the important information for defining marine conservation units both at ocean scale and at the level of small islands.

Hidden historical habitat-linked population divergence and contemporary gene flow of a deep-sea patellogastropod limpet

Hydrothermal vents and hydrocarbon seeps in the deep ocean are rare oases fueled by chemosynthesis. Biological communities inhabiting these ecosystems are often distributed in widely separated habitats, raising intriguing questions on how these organisms achieve connectivity and whether habitat types facilitate intraspecific divergence. The deep-sea patellogastropod limpet Bathyacmaea nipponica that colonizes both vents and seeps across ∼2,400 km in the Northwest Pacific is a feasible model to answer these questions. We analyzed 123 individuals from four vents and three seeps using a comprehensive method incorporating population genomics and physical ocean modeling. Genome survey sequencing and genotyping-by-sequencing resulted in 9,838 single-nucleotide polymorphisms for population genomic analyses. Genetic divergence and demographic analyses revealed four habitat-linked (i.e., three seep and one vent) genetic groups, with the vent genetic group established via the opportunistic invasion of a few limpet larvae from a nearby seep genetic group. TreeMix analysis uncovered three historical seep-to-vent migration events. ADMIXTURE and divMigrate analyses elucidated weak contemporary gene flow from a seep genetic group to the vent genetic group. Physical ocean modeling underlined the potential roles of seafloor topography and ocean currents in shaping the genetic connectivity, contemporary migration, and local hybridization of these deep-sea limpets. Our study highlighted the power of integrating genomic and oceanographic approaches in deciphering the demography and diversification of deep-sea organisms. Given the increasing anthropogenic activities (e.g., mining and gas hydrate extraction) affecting the deep ocean, our results have implications for the conservation of deep-sea biodiversity and establishment of marine protected areas.

Unraveling a genetic roadmap for improved taste in the domesticated apple

Although taste is an important aspect of fruit quality, an understanding of its genetic control remains elusive in apple and other fruit crops. In this study, we conducted genomic sequence analysis of 497 Malus accessions and revealed erosion of genetic diversity caused by apple breeding and possible independent domestication events of dessert and cider apples. Signatures of selection for fruit acidity and size, but not for fruit sugar content, were detected during the processes of both domestication and improvement. Furthermore, we found that single mutations in major genes affecting fruit taste, including Ma1, MdTDT, and MdSOT2, dramatically decrease malate, citrate, and sorbitol accumulation, respectively, and correspond to important domestication events. Interestingly, Ma1 was identified to have pleiotropic effects on both organic acid content and sugar:acid ratio, suggesting that it plays a vital role in determining fruit taste. Fruit taste is unlikely to have been negatively affected by linkage drag associated with selection for larger fruit that resulted from the pyramiding of multiple genes with minor effects on fruit size. Collectively, our study provides new insights into the genetic basis of fruit quality and its evolutionary roadmap during apple domestication, pinpointing several candidate genes for genetic manipulation of fruit taste in apple.

Genetic Differentiation and Demographic History of Three Cerris Oak Species in China Based on Nuclear Microsatellite Makers

Knowledge of interspecific divergence and population expansions/contractions of dominant forest trees in response to geological events and climatic oscillations is of major importance to understand their evolution and demography. However, the interspecific patterns of genetic differentiation and spatiotemporal population dynamics of three deciduous Cerris oak species (Q. acutissima, Q. variabilis and Q. chenii) that are widely distributed in China remain poorly understood. In this study, we genotyped 16 nuclear loci in 759 individuals sampled from 44 natural populations of these three sibling species to evaluate the plausible demographical scenarios of the closely related species. We also tested the hypothesis that macro- and microevolutionary processes of the three species had been triggered and molded by Miocene–Pliocene geological events and Quaternary climatic change. The Bayesian cluster analysis showed that Q. acutissima and Q. chenii were clustered in the same group, whereas Q. variabilis formed a different genetic cluster. An approximate Bayesian computation (ABC) analyses suggested that Q. variabilis and Q. acutissima diverged from their most common ancestor around 19.84 Ma, and subsequently Q. chenii diverged from Q. acutissima at about 9.6 Ma, which was significantly associated with the episodes of the Qinghai–Tibetan Plateau (QTP). In addition, ecological niche modeling and population history analysis showed that these three Cerris oak species repeatedly underwent considerable ‘expansion–contraction’ during the interglacial and glacial periods of the Pleistocene, although they have varying degrees of tolerance for the climatic change. Overall, these findings indicated geological and climatic changes during the Miocene–Pliocene and Pleistocene as causes of species divergence and range shifts of dominant tree species in the subtropical and warm temperature areas in China.

Genetic diversity and population structure of Alternaria species from tomato and potato in North Carolina and Wisconsin

Early blight (EB) caused by Alternaria linariae or Alternaria solani and leaf blight (LB) caused by A. alternata are economically important diseases of tomato and potato. Little is known about the genetic diversity and population structure of these pathogens in the United States. A total of 214 isolates of A. alternata (n = 61), A. linariae (n = 96), and A. solani (n = 57) were collected from tomato and potato in North Carolina and Wisconsin and grouped into populations based on geographic locations and tomato varieties. We exploited 220 single nucleotide polymorphisms derived from DNA sequences of 10 microsatellite loci to analyse the population genetic structure between species and between populations within species and infer the mode of reproduction. High genetic variation and genotypic diversity were observed in all the populations analysed. The null hypothesis of the clonality test based on the index of association \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {\overline{r}_{d} } \right)$$\end{document}r¯d was rejected, and equal frequencies of mating types under random mating were detected in some studied populations of Alternaria spp., suggesting that recombination can play an important role in the evolution of these pathogens. Most genetic differences were found between species, and the results showed three distinct genetic clusters corresponding to the three Alternaria spp. We found no evidence for clustering of geographic location populations or tomato variety populations. Analyses of molecular variance revealed high (> 85%) genetic variation within individuals in a population, confirming a lack of population subdivision within species. Alternaria linariae populations harboured more multilocus genotypes (MLGs) than A. alternata and A. solani populations and shared the same MLG between populations within a species, which was suggestive of gene flow and population expansion. Although both A. linariae and A. solani can cause EB on tomatoes and potatoes, these two species are genetically differentiated. Our results provide new insights into the evolution and structure of Alternaria spp. and can lead to new directions in optimizing management strategies to mitigate the impact of these pathogens on tomato and potato production in North Carolina and Wisconsin.

Genomic divergence during feralization reveals both conserved and distinct mechanisms of parallel weediness evolution

Agricultural weeds are the most important biotic constraints to global crop production, and chief among these is weedy rice. Despite increasing yield losses from weedy rice in recent years worldwide, the genetic basis of weediness evolution remains unclear. Using whole-genome sequence analyses, we examined the origins and adaptation of Japanese weedy rice. We find evidence for a weed origin from tropical japonica crop ancestry, which has not previously been documented in surveys of weedy rice worldwide. We further show that adaptation occurs largely through different genetic mechanisms between independently-evolved temperate japonica- and tropical japonica-derived strains; most genomic signatures of positive selection are unique within weed types. In addition, some weedy rice strains have evolved through hybridization between weedy and cultivated rice with adaptive introgression from the crop. Surprisingly, introgression from cultivated rice confers not only crop-like adaptive traits (such as shorter plant height, facilitating crop mimicry) but also weedy-like traits (such as seed dormancy). These findings reveal how hybridization with cultivated rice can promote persistence and proliferation of weedy rice.

Genetic divergence, admixture and subspecific boundaries in a peripheral population of the great tit, Parus major (Aves: Paridae)

Secondary contact zones have been formed between several pairs of avian species and subspecies in northern and north-eastern Iran during the post-Pleistocene and Holocene periods. Three subspecies groups out of the four in the great tit (Parus major), major, bokharensis and cinereus, are believed to have come into local or regional secondary contact in north-eastern Iran. Parus major intermedius is also known from this region and has long been suggested to have a hybrid origin from hybridization between the yellow western (major) subspecies group and the grey-coloured eastern (bokharensis or cinereus) subspecies group based on its intermediate plumage coloration. Here, we investigated the evidence of intergradation between subspecies groups and the evolutionary origin of P. m. intermedius using mitochondrial DNA sequences and microsatellites, and approximate Bayesian computation to test competing scenarios for the demographic history of the populations. Our analyses indicate a divergence origin for intermedius that resulted from expansion of the major subspecies group. Low mitochondrial diversity and high genetic differentiation in comparison with central populations suggest that intermedius is a peripheral population. Microsatellite data show no signs of nuclear admixture between the bokharensis and major subspecies groups. However, one phenotypically intermedius specimen had bokharensis mtDNA and major nuclear DNA in the assumed hybrid zone (Lotf Abad), supporting past introgression.

Addressing alpine plant phylogeography using integrative distributional, demographic and coalescent modeling

Phylogeographic studies of alpine plants have evolved considerably in the last two decades from ad hoc interpretations of genetic data to statistical model-based approaches. In this review we outline the developments in alpine plant phylogeography focusing on the recent approach of integrative distributional, demographic and coalescent (iDDC) modeling. By integrating distributional data with spatially explicit demographic modeling and subsequent coalescent simulations, the history of alpine species can be inferred and long-standing hypotheses, such as species-specific responses to climate change or survival on nunataks during the last glacial maximum, can be efficiently tested as exemplified by available case studies. We also discuss future prospects and improvements of iDDC.

Genetic Structure and Diversity of the Yellowbelly Threadfin Bream Nemipterus bathybius in the Northern South China Sea

The genetic structure and demography of the yellowbelly threadfin bream, Nemipterus bathybius, in the northern South China Sea were examined using the mitochondrial DNA cytochrome b gene (1141 bp). High levels of haplotype and nucleotide diversities (0.98 and 5.26 × 10−3, respectively) showed that all populations exhibited a high level of genetic diversity. Analysis of molecular variance (AMOVA), FST statistics, and haplotype networks suggested the absence of significant genetic differentiation along the coast of the northern South China Sea. Although the results suggested that the lack of differentiation within the population structure of N. bathybius was shaped by ocean currents, our results also showed that the Qiongzhou Strait limited their migration between Beibu Gulf and the northern South China Sea. Neutrality tests and mismatch distributions indicated population expansion, but the Bayesian skyline plots and approximate Bayesian computation approaches suggested that the population sizes recently contracted. The diversification of multiple stocks, which were induced by two ocean current systems, contributed to these discordant results. Although these analyses of demographic history revealed no evidence for recent population bottlenecks, the population demography needs to be evaluated further.

Genomic Architecture of Rapid Parallel Adaptation to Fresh Water in a Wild Fish

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. However, the genetic architecture underlying these adaptive changes is still poorly understood. Using population genomic approaches, we investigated the genomic architecture that underlies rapid parallel adaptation of Coilia nasus to fresh water by comparing four freshwater-resident populations with their ancestral anadromous population. Linkage disequilibrium network analysis and population genetic analyses revealed two putative large chromosome inversions on LG6 and LG22, which were enriched for outlier loci and exhibited parallel association with freshwater adaptation. Drastic frequency shifts and elevated genetic differentiation were observed for the two chromosome inversions among populations, suggesting that both inversions would undergo divergent selection between anadromous and resident ecotypes. Enrichment analysis of genes within chromosome inversions showed significant enrichment of genes involved in metabolic process, immunoregulation, growth, maturation, osmoregulation, and so forth, which probably underlay differences in morphology, physiology and behavior between the anadromous and freshwater-resident forms. The availability of beneficial standing genetic variation, large optimum shift between marine and freshwater habitats, and high efficiency of selection with large population size could lead to the observed rapid parallel adaptive genomic change. We propose that chromosomal inversions might have played an important role during the evolution of rapid parallel ecological divergence in the face of environmental heterogeneity in C. nasus. Our study provides insights into the genomic basis of rapid adaptation of complex traits in novel habitats and highlights the importance of structural genomic variants in analyses of ecological adaptation.

"Jumping Jack": Genomic Microsatellites Underscore the Distinctiveness of Closely Related Pseudoperonospora cubensis and Pseudoperonospora humuli and Provide New Insights Into Their Evolutionary Past

Downy mildews caused by obligate biotrophic oomycetes result in severe crop losses worldwide. Among these pathogens, Pseudoperonospora cubensis and P. humuli, two closely related oomycetes, adversely affect cucurbits and hop, respectively. Discordant hypotheses concerning their taxonomic relationships have been proposed based on host-pathogen interactions and specificity evidence and gene sequences of a few individuals, but population genetics evidence supporting these scenarios is missing. Furthermore, nuclear and mitochondrial regions of both pathogens have been analyzed using microsatellites and phylogenetically informative molecular markers, but extensive comparative population genetics research has not been done. Here, we genotyped 138 current and historical herbarium specimens of those two taxa using microsatellites (SSRs). Our goals were to assess genetic diversity and spatial distribution, to infer the evolutionary history of P. cubensis and P. humuli, and to visualize genome-scale organizational relationship between both pathogens. High genetic diversity, modest gene flow, and presence of population structure, particularly in P. cubensis, were observed. When tested for cross-amplification, 20 out of 27 P. cubensis-derived gSSRs cross-amplified DNA of P. humuli individuals, but few amplified DNA of downy mildew pathogens from related genera. Collectively, our analyses provided a definite argument for the hypothesis that both pathogens are distinct species, and suggested further speciation in the P. cubensis complex.

Moderate Genetic Diversity and Demographic Reduction in the Threatened Giant Anteater, Myrmecophaga tridactyla

In general, large mammal species with highly specialized feeding behavior and solitary habits are expected to suffer genetic consequences from habitat loss and fragmentation. To test this hypothesis, we analyzed the genetic diversity distribution of the threatened giant anteater inhabiting a human-modified landscape. We used 10 microsatellite loci to assess the genetic diversity and population structure of 107 giant anteaters sampled in the Brazilian Central-Western region. No genetic population structuring was observed in this region suggesting no gene flow restriction within the studied area. On the other hand, the moderate level of genetic diversity (Ho = 0.54), recent bottleneck detected and inbreeding (F_is, 0.13; p ≤ 0.001) signatures suggest potential impacts on the genetic variation of this Xenarthra. Additionally, a previous demographic reduction was suggested. Thus, considering the increased human-promoted impacts across the entire area of distribution of the giant anteater, our results can illustrate the potential effects of these disturbances on the genetic variation, allowing us to request the long-term conservation of this emblematic species.

Extensive human-mediated jump dispersal within and across the native and introduced ranges of the invasive termite Reticulitermes flavipes

As native ranges are often geographically structured, invasive species originating from a single source population only carry a fraction of the genetic diversity present in their native range. The invasion process is thus often associated with a drastic loss of genetic diversity resulting from a founder event. However, the fraction of diversity brought to the invasive range may vary under different invasion histories, increasing with the size of the propagule, the number of reintroduction events, and/or the total genetic diversity represented by the various source populations in a multiple-introduction scenario. In this study, we generated a SNP data set for the invasive termite Reticulitermes flavipes from 23 native populations in the eastern United States and six introduced populations throughout the world. Using population genetic analyses and approximate Bayesian computation random forest, we investigated its worldwide invasion history. We found a complex invasion pathway with multiple events out of the native range and bridgehead introductions from the introduced population in France. Our data suggest that extensive long-distance jump dispersal appears common in both the native and introduced ranges of this species, probably through human transportation. Overall, our results show that similar to multiple introduction events into the invasive range, admixture in the native range prior to invasion can potentially favour invasion success by increasing the genetic diversity that is later transferred to the introduced range.

Climate-induced range shifts shaped the present and threaten the future genetic variability of a marine brown alga in the Northwest Pacific

Glaciation-induced environmental changes during the last glacial maximum (LGM) have strongly influenced species' distributions and genetic diversity patterns in the northern high latitudes. However, these effects have seldom been assessed on sessile species in the Northwest Pacific. Herein, we chose the brown alga Sargassum thunbergii to test this hypothesis, by comparing present population genetic variability with inferred geographical range shifts from the LGM to the present, estimated with species distribution modelling (SDM). Projections for contrasting scenarios of future climate change were also developed to anticipate genetic diversity losses at regional scales. Results showed that S. thunbergii harbours strikingly rich genetic diversity and multiple divergent lineages in the centre-northern range of its distribution, in contrast with a poorer genetically distinct lineage in the southern range. SDM hindcasted refugial persistence in the southern range during the LGM as well as post-LGM expansion of 18 degrees of latitude northward. Approximate Bayesian computation (ABC) analysis further suggested that the multiple divergent lineages in the centre-northern range limit stem from post-LGM colonization from the southern survived lineage. This suggests divergence due to demographic bottlenecks during range expansion and massive genetic diversity loss during post-LGM contraction in the south. The projected future range of S. thunbergii highlights the threat to unique gene pools that might be lost under global changes.

Contrasting Phylogeographic Patterns in Lumnitzera Mangroves Across the Indo-West Pacific

Mangroves are ecologically important forest communities in tropical and subtropical coasts, the effective management of which requires understanding of their phylogeographic patterns. However, these patterns often vary among different species, even among ecologically similar taxa or congeneric species. Here, we investigated the levels and patterns of genetic variation within Lumnitzera consisting of two species (L. racemosa and L. littorea) with nearly sympatric ranges across the Indo-West Pacific (IWP) region by sequencing three chloroplast DNA regions (for both species) and genotyping 11 nuclear microsatellite loci (for L. littorea). Consistent with findings in studies on other mangrove species, we found that both L. racemosa and L. littorea showed relatively high genetic variation among populations but low genetic variation within populations. Haplotype network and genetic clustering analyses indicated two well-differentiated clades in both L. racemosa and L. littorea. The relationship between geographic and genetic distances and divergence time estimates of the haplotypes indicated that limited dispersal ability of the propagules, emergence of land barriers during ancient sea-level changes, and contemporary oceanic circulation pattern in the IWP influenced the current population structure of the two species. However, the position of genetic break was found to vary between the two species: in L. racemosa, strong divergence was observed between populations from the Indian Ocean and the Pacific Ocean possibly due to land barrier effect of the Malay Peninsula; in L. littorea, the phylogeographic pattern was created by a more eastward genetic break along the biogeographic barrier identified as the Huxley's line. Overall, our findings strongly supported previous hypothesis of mangrove species divergence and revealed that the two Lumnitzera species have different phylogeographic patterns despite their close genetic relationship and similar current geographic distribution. The findings also provided references for the management of Lumnitzera mangroves, especially for the threatened L. littorea.

The demographic history of Castanea sativa Mill. in southwest Europe: A natural population structure modified by translocations

In domesticated species, translocation of materials can alter natural demographic patterns; this may have occurred in Castanea sativa (European chestnut), a species conserved in several refuges, with a long domestication history for nut production. Bayesian analysis of population genetic structure in marginal areas and in the centre of C. sativa range, considering spatial information and making corrections for unbalanced size, allowed visualization of a genetic structure that was subsequently analysed by approximate Bayesian computation to assess its natural demographic history and test the origin of some hypothetical translocated germplasm. We obtained evidence of C. sativa population contraction during the earliest Pleistocene, resulting in a split into eastern (Greek) and western (Italian and Iberian) populations. The northern Iberian population, in the Eurosiberian area, is one of the vestiges that remained after the contraction that split the global Tertiary population. A secondary encounter, occurred from Middle to Upper Pleistocene, which explains the hybrid origin of the Western Mediterranean population present in Italy and in the centre and south of the Iberian Peninsula. It has been demonstrated that a germplasm translocation from Italy to the Central Iberian Range may have occurred. Recent translocations have hybridized with the local northern Iberian population, producing naturalized populations with high diversity. The populations of C. sativa in southwestern Europe have a genetic structure compatible with a natural origin, in which signs of population contractions and expansions caused by climatic oscillations since the Late Miocene have been imprinted.

Mitochondrial introgression by ancient admixture between two distant lacustrine fishes in Sulawesi Island

Sulawesi, an island located in a biogeographical transition zone between Indomalaya and Australasia, is famous for its high levels of endemism. Ricefishes (family Adrianichthyidae) are an example of taxa that have uniquely diversified on this island. It was demonstrated that habitat fragmentation due to the Pliocene juxtaposition among tectonic subdivisions of this island was the primary factor that promoted their divergence; however, it is also equally probable that habitat fusions and resultant admixtures between phylogenetically distant species may have frequently occurred. Previous studies revealed that some individuals of Oryzias sarasinorum endemic to a tectonic lake in central Sulawesi have mitochondrial haplotypes that are similar to the haplotypes of O. eversi, which is a phylogenetically related but geologically distant (ca. 190 km apart) adrianichthyid endemic to a small fountain. In this study, we tested if this reflects ancient admixture of O. eversi and O. sarasinorum. Population genomic analyses of genome-wide single-nucleotide polymorphisms revealed that O. eversi and O. sarasinorum are substantially reproductively isolated from each other. Comparison of demographic models revealed that the models assuming ancient admixture from O. eversi to O. sarasinorum was more supported than the models assuming no admixture; this supported the idea that the O. eversi-like mitochondrial haplotype in O. sarasinorum was introgressed from O. eversi. This study is the first to demonstrate ancient admixture of lacustrine or pond organisms in Sulawesi beyond 100 km. The complex geological history of this island enabled such island-wide admixture of lacustrine organisms, which usually experience limited migration.

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.

Genomics-informed models reveal extensive stretches of coastline under threat by an ecologically dominant invasive species

Explaining why some species are widespread, while others are not, is fundamental to biogeography, ecology, and evolutionary biology. A unique way to study evolutionary and ecological mechanisms that either limit species' spread or facilitate range expansions is to conduct research on species that have restricted distributions. Nonindigenous species, particularly those that are highly invasive but have not yet spread beyond the introduced site, represent ideal systems to study range size changes. Here, we used species distribution modeling and genomic data to study the restricted range of a highly invasive Australian marine species, the ascidian Pyura praeputialis This species is an aggressive space occupier in its introduced range (Chile), where it has fundamentally altered the coastal community. We found high genomic diversity in Chile, indicating high adaptive potential. In addition, genomic data clearly showed that a single region from Australia was the only donor of genotypes to the introduced range. We identified over 3,500 km of suitable habitat adjacent to its current introduced range that has so far not been occupied, and importantly species distribution models were only accurate when genomic data were considered. Our results suggest that a slight change in currents, or a change in shipping routes, may lead to an expansion of the species' introduced range that will encompass a vast portion of the South American coast. Our study shows how the use of population genomics and species distribution modeling in combination can unravel mechanisms shaping range sizes and forecast future range shifts of invasive species.

Low levels of genetic differentiation with isolation by geography and environment in populations of Drosophila melanogaster from across China

The fruit fly, Drosophila melanogaster, is a model species in evolutionary studies. However, population processes of this species in East Asia are poorly studied. Here we examined the population genetic structure of D. melanogaster across China. There were 14 mitochondrial haplotypes with 10 unique ones out of 23 known from around the globe. Pairwise F_ST values estimated from 15 novel microsatellites ranged from 0 to 0.11, with geographically isolated populations showing the highest level of genetic uniqueness. STRUCTURE analysis identified high levels of admixture at both the individual and population levels. Mantel tests indicated a strong association between genetic distance and geographical distance as well as environmental distance. Full redundancy analysis (RDA) showed that independent effects of environmental conditions and geography accounted for 62.10% and 31.58% of the total explained genetic variance, respectively. When geographic variables were constrained in a partial RDA analysis, the environmental variables bio2 (mean diurnal air temperature range), bio13 (precipitation of the wettest month), and bio15 (precipitation seasonality) were correlated with genetic distance. Our study suggests that demographic history, geographical isolation, and environmental factors have together shaped the population genetic structure of D. melanogaster after its introduction into China.

The DNA history of a lonely oak: Quercus humboldtii phylogeography in the Colombian Andes

The climatic and geological changes that occurred during the Quaternary, particularly the fluctuations during the glacial and interglacial periods of the Pleistocene, shaped the population demography and geographic distribution of many species. These processes have been studied in several groups of organisms in the Northern Hemisphere, but their influence on the evolution of Neotropical montane species and ecosystems remains unclear. This study contributes to the understanding of the effect of climatic fluctuations during the late Pleistocene on the evolution of Andean mountain forests. First, we describe the nuclear and plastidic DNA patterns of genetic diversity, structure, historical demography, and landscape connectivity of Quercus humboldtii, which is a typical species in northern Andean montane forests. Then, these patterns were compared with the palynological and evolutionary hypotheses postulated for montane forests of the Colombian Andes under climatic fluctuation scenarios during the Quaternary. Our results indicated that populations of Q. humboldtii have high genetic diversity and a lack of genetic structure and that they have experienced a historical increase in connectivity from the last glacial maximum (LGM) to the present. Furthermore, our results showed a dramatic reduction in the effective population size followed by an expansion before the LGM, which is consistent with the results found by palynological studies, suggesting a change in dominance in Andean forests that may be related to ecological factors rather than climate change.

A dynamic history of admixture from Mediterranean and Carpathian glacial refugia drives genomic diversity in the bank vole

Understanding the historical contributions of differing glacial refugia is key to evaluating the roles of microevolutionary forces, such as isolation, introgression, and selection in shaping genomic diversity in present-day populations. In Europe, where both Mediterranean and extra-Mediterranean (e.g., Carpathian) refugia of the bank vole (Clethrionomys glareolus) have been identified, mtDNA indicates that extra-Mediterranean refugia were the main source of colonization across the species range, while Mediterranean peninsulas harbor isolated, endemic lineages. Here, we critically evaluate this hypothesis using previously generated genomic data (>6,000 SNPs) for over 800 voles, focusing on genomic contributions to bank voles in central Europe, a key geographic area in considering range-wide colonization. The results provide clear evidence that both extra-Mediterranean (Carpathian) and Mediterranean (Spanish, Calabrian, and Balkan) refugia contributed to the ancestry and genomic diversity of bank vole populations across Europe. Few strong barriers to dispersal and frequent admixture events in central Europe have led to a prominent mid-latitude peak in genomic diversity. Although the genomic contribution of the centrally located Carpathian refugium predominates, populations in different parts of Europe have admixed origins from Mediterranean (28%-47%) and the Carpathian (53%-72%) sources. We suggest that the admixture from Mediterranean refugia may have provisioned adaptive southern alleles to more northern populations, facilitating the end-glacial spread of the admixed populations and contributing to increased bank vole diversity in central Europe. This study adds critical details to the complex end-glacial colonization history of this well-studied organism and underscores the importance of genomic data in phylogeographic interpretation.

The phylogeographic history of a range disjunction in eastern North America: the role of post-glacial expansion into newly suitable habitat

PREMISE

A disjunct distribution, where a species' geographic range is discontinuous, can occur through vicariance or long-distance dispersal. Approximately 75 North American plant species exhibit a ~650 km disjunction between the Ozark and Appalachian regions. This disjunction is attributed to biogeographic forces including: (1) Eocene-Oligocene vicariance by the formation of the Mississippi embayment; (2) Pleistocene vicariance from interglacial flooding; (3) post-Pleistocene northward colonization from separate glacial refugia; (4) Hypsithermal vicariance due to climate fluctuations; and (5) recent long-distance dispersal. We investigated which of these pathways most likely gave rise to the Appalachian-Ozark disjunction in Delphinium exaltatum.

METHODS

We genotyped populations of D. exaltatum from five Ozark and seven Appalachian localities, analyzed genetic structure, tested the order and timing of divergences using DIYABC, and conducted niche reconstructions up to 21,000 years before present (YBP).

RESULTS

Populations fell into five main genetic clusters, i.e., a group in the central Appalachians, and four "lowland" groups. DIYABC analyses showed the central Appalachian and lowland lineages diverging 11,300 to17,000 YBP, and the lowland groups diverging 6800 to 10,900 YBP. Niche reconstructions showed that suitable climate for the central Appalachian lineage experienced large spatial discontinuity starting 14,000 YBP, such that divergence and persistence before this period is less plausible than divergence thereafter.

CONCLUSIONS

Our results did not fully support any of the original hypotheses. Rather, the oldest divergence likely occurred after 13,500 YBP through expansion into newly opened habitat in the Appalachians. The Appalachian-Ozark disjunction likely resulted from northward dispersal of the lowland lineage as climate warmed during the Holocene.

Speciation along a latitudinal gradient: The origin of the Neotropical cycad sister pair Dioon sonorense-D. vovidesii (Zamiaceae)

Latitude is correlated with environmental components that determine the distribution of biodiversity. In combination with geographic factors, latitude-associated environmental variables are expected to influence speciation, but empirical evidence on how those factors interplay is scarce. We evaluated the genetic and environmental variation among populations in the pair of sister species Dioon sonorense-D. vovidesii, two cycads distributed along a latitudinal environmental gradient in northwestern Mexico, to reveal their demographic histories and the environmental factors involved in their divergence. Using genome-wide loci data, we determined the species delimitation, estimated the gene flow, and compared multiple demographic scenarios of divergence. Also, we estimated the variation of climatic variables among populations and used ecological niche models to test niche overlap between species. The effect of geographic and environmental variables on the genetic variation among populations was evaluated using linear models. Our results suggest the existence of a widespread ancestral population that split into the two species ~829 ky ago. The geographic delimitation along the environmental gradient occurs in the absence of major geographic barriers, near the 28th parallel north, where a zonation of environmental seasonality exists. The northern species, D. vovidesii, occurs in more seasonal environments but retains the same niche of the southern species, D. sonorense. The genetic variation throughout populations cannot be solely explained by stochastic processes; the latitudinal-associated seasonality has been an additive factor that strengthened the species divergence. This study represents an example of how speciation can be achieved by the effect of the latitude-associated factors on the genetic divergence among populations.

Demographic history and local adaptation of Myripnois dioica (Asteraceae) provide insight on plant evolution in northern China flora

The flora of northern China forms the main part of the Sino-Japanese floristic region and is located in a south-north vegetative transect in East Asia. Phylogeographic studies have demonstrated that an arid belt in this region has promoted divergence of plants in East Asia. However, little is known about how plants that are restricted to the arid belt of flora in northern China respond to climatic oscillation and environmental change. Here, we used genomic-level data of Myripnois dioica across its distribution as a representative of northern China flora to reconstruct plant demographic history, examine local adaptation related to environmental disequilibrium, and investigate the factors related to effective population size change. Our results indicate M. dioica originated from the northern area and expanded to the southern area, with the Taihang Mountains serving as a physical barrier promoting population divergence. Genome-wide evidence found strong correlation between genomic variation and environmental factors, specifically signatures associated with local adaptation to drought stress in heterogeneous environments. Multiple linear regression analyses revealed joint effects of population age, mean temperature of coldest quarter, and precipitation of wettest month on effective population size (Ne). Our current study uses M. dioica as a case for providing new insights into the evolutionary history and local adaptation of northern China flora and provides qualitative strategies for plant conservation.

Genetic divergence between isolated populations of the North Island New Zealand Rifleman (Acanthisitta chloris granti) implicates ancient biogeographic impacts rather than recent habitat fragmentation

This research investigates the extent and causal mechanisms of genetic population divergence in a poorly flighted passerine, the North Island Rifleman or Titipounamu (Acanthisitta chloris granti). While this species has a historically widespread distribution, anthropogenic forest clearance has resulted in a highly fragmented current distribution. We conducted analyses of mitochondrial DNA (COI and Control Region) and 12 nuclear DNA microsatellites to test for population divergence and estimate times of divergence. diyabc and biogeobears were then used to assess likely past dispersal scenarios based on both mtDNA and nDNA. The results reveal several significantly divergent lineages across the North Island of New Zealand and indicate that some populations have been isolated for extensive periods of time (0.7-4.9 mya). Modeling indicated a dynamic history of population connectivity, with a drastic restriction in gene flow between three geographic regions, followed by a more recent re-establishment of connectivity. Our analyses indicate the dynamic influence of key geological and climatological events on the distribution of genetic diversity in this species, including support for the genetic impact of old biogeographic boundaries such as the Taupo Line and Cockayne's Line, rather than recent anthropogenic habitat fragmentation. These findings present a rare example of an avian species with a genetic history more like that of flightless taxa and so provide new general insights into vicariant processes affecting populations of passerines with limited dispersal.

Hierarchical genetic structure in an evolving species complex: Insights from genome wide ddRAD data in Sebastes mentella

The diverse biology and ecology of marine organisms may lead to complex patterns of intraspecific diversity for both neutral and adaptive genetic variation. Sebastes mentella displays a particular life-history as livebearers, for which existence of multiple ecotypes has been suspected to complicate the genetic population structure of the species. Double digest restriction-site associated DNA was used to investigate genetic population structure in S. mentella and to scan for evidence of selection. In total, 42,288 SNPs were detected in 277 fish, and 1,943 neutral and 97 tentatively adaptive loci were selected following stringent filtration. Unprecedented levels of genetic differentiation were found among the previously defined ‘shallow pelagic’, ‘deep pelagic’ and ‘demersal slope’ ecotypes, with overall mean F_ST = 0.05 and 0.24 in neutral and outlier SNPs, respectively. Bayesian computation estimated a concurrent and historical divergence among these three ecotypes and evidence of local adaptation was found in the S. mentella genome. Overall, these findings imply that the depth-defined habitat divergence of S. mentella has led to reproductive isolation and possibly adaptive radiation among these ecotypes. Additional sub-structuring was detected within the ‘shallow’ and ‘deep’ pelagic ecotypes. Population assignment of individual fish showed more than 94% agreement between results based on SNP and previously generated microsatellite data, but the SNP data provided a lower estimate of hybridization among the ecotypes than that by microsatellite data. We identified a SNP panel with only 21 loci to discriminate populations in mixed samples based on a machine-learning algorithm. This first SNP based investigation clarifies the population structure of S. mentella, and provides novel and high-resolution genomic tools for future investigations. The insights and tools provided here can readily be incorporated into the management of S. mentella and serve as a template for other exploited marine species exhibiting similar complex life history traits.

Discovery of cryptic plant diversity on the rooftops of the Alps

High elevation temperate mountains have long been considered species poor owing to high extinction or low speciation rates during the Pleistocene. We performed a phylogenetic and population genomic investigation of an emblematic high-elevation plant clade (Androsace sect. Aretia, 31 currently recognized species), based on plant surveys conducted during alpinism expeditions. We inferred that this clade originated in the Miocene and continued diversifying through Pleistocene glaciations, and discovered three novel species of Androsace dwelling on different bedrock types on the rooftops of the Alps. This highlights that temperate high mountains have been cradles of plant diversity even during the Pleistocene, with in-situ speciation driven by the combined action of geography and geology. Our findings have an unexpected historical relevance: H.-B. de Saussure likely observed one of these species during his 1788 expedition to the Mont Blanc and we describe it here, over two hundred years after its first sighting.

Genome-wide SNPs redefines species boundaries and conservation units in the freshwater mussel genus Cyprogenia of North America

Detailed information on species delineation and population genetic structure is a prerequisite for designing effective restoration and conservation strategies for imperiled organisms. Phylogenomic and population genomic analyses based on genome-wide double digest restriction-site associated DNA sequencing (ddRAD-Seq) data has identified three allopatric lineages in the North American freshwater mussel genus Cyprogenia. Cyprogenia stegaria is restricted to the Eastern Highlands and displays little genetic structuring within this region. However, two allopatric lineages of C. aberti in the Ozark and Ouachita highlands exhibit substantial levels (mean uncorrected F_ST = 0.368) of genetic differentiation and each warrants recognition as a distinct evolutionary lineage. Lineages of Cyprogenia in the Ouachita and Ozark highlands are further subdivided reflecting structuring at the level of river systems. Species tree inference and species delimitation in a Bayesian framework using single nucleotide polymorphisms (SNP) data supported results from phylogenetic analyses, and supports three species of Cyprogenia over the currently recognized two species. A comparison of SNPs generated from both destructively and non-destructively collected samples revealed no significant difference in the SNP error rate, quality and amount of ddRAD sequence reads, indicating that nondestructive or trace samples can be effectively utilized to generate SNP data for organisms for which destructive sampling is not permitted.

Recent demographic histories of temperate deciduous trees inferred from microsatellite markers

Background

Accurate inference of demographic histories for temperate tree species can aid our understanding of current climate change as a driver of evolution. Microsatellites are more suitable for inferring recent historical events due to their high mutation rates. However, most programs analyzing microsatellite data assume a strict stepwise mutation model (SMM), which could cause false detection of population shrinkage when microsatellite mutation does not follow SMM.

Results

This study aims to reconstruct the recent demographic histories of five cool-temperate tree species in Eastern Asia, Quercus mongolica, Q. liaotungensis, Juglans cathayensis, J. mandshurica and J. ailantifolia, by using 19 microsatellite markers with two methods considering generalized stepwise mutation model (GSM) (MIGRAINE and VarEff). Both programs revealed that all the five species experienced expansions after the Last Glacial Maximum (LGM). Within butternuts, J. cathayensis experienced a more serious bottleneck than the other species, and within oaks, Q. mongolica showed a moderate increase in population size and remained stable after the expansion. In addition, the point estimates of the multistep mutation proportion in the GSM model (p_GSM) for all five species were between 0.50 and 0.65, indicating that when inferring population demographic history of the cool-temperate forest species using microsatellite markers, it is better to assume a GSM rather than a SMM.

Conclusions

This study provides the first direct evidence that five cool-temperate tree species in East Asia have experienced expansions after the LGM with microsatellite data. Considering the mutation model of microsatellite has a vital influence on demographic inference, combining multiple programs such as MIGRAINE and VarEff can effectively reduce errors caused by inappropriate model selection and prior setting.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01805-w.

Extending approximate Bayesian computation with supervised machine learning to infer demographic history from genetic polymorphisms using DIYABC Random Forest

Simulation-based methods such as approximate Bayesian computation (ABC) are well-adapted to the analysis of complex scenarios of populations and species genetic history. In this context, supervised machine learning (SML) methods provide attractive statistical solutions to conduct efficient inferences about scenario choice and parameter estimation. The Random Forest methodology (RF) is a powerful ensemble of SML algorithms used for classification or regression problems. Random Forest allows conducting inferences at a low computational cost, without preliminary selection of the relevant components of the ABC summary statistics, and bypassing the derivation of ABC tolerance levels. We have implemented a set of RF algorithms to process inferences using simulated data sets generated from an extended version of the population genetic simulator implemented in DIYABC v2.1.0. The resulting computer package, named DIYABC Random Forest v1.0, integrates two functionalities into a user-friendly interface: the simulation under custom evolutionary scenarios of different types of molecular data (microsatellites, DNA sequences or SNPs) and RF treatments including statistical tools to evaluate the power and accuracy of inferences. We illustrate the functionalities of DIYABC Random Forest v1.0 for both scenario choice and parameter estimation through the analysis of pseudo-observed and real data sets corresponding to pool-sequencing and individual-sequencing SNP data sets. Because of the properties inherent to the implemented RF methods and the large feature vector (including various summary statistics and their linear combinations) available for SNP data, DIYABC Random Forest v1.0 can efficiently contribute to the analysis of large SNP data sets to make inferences about complex population genetic histories.

Global population genetic structure and demographic trajectories of the black soldier fly, Hermetia illucens

BACKGROUND

The black soldier fly (Hermetia illucens) is the most promising insect candidate for nutrient-recycling through bioconversion of organic waste into biomass, thereby improving sustainability of protein supplies for animal feed and facilitating transition to a circular economy. Contrary to conventional livestock, genetic resources of farmed insects remain poorly characterised. We present the first comprehensive population genetic characterisation of H. illucens. Based on 15 novel microsatellite markers, we genotyped and analysed 2862 individuals from 150 wild and captive populations originating from 57 countries on seven subcontinents.

RESULTS

We identified 16 well-distinguished genetic clusters indicating substantial global population structure. The data revealed genetic hotspots in central South America and successive northwards range expansions within the indigenous ranges of the Americas. Colonisations and naturalisations of largely unique genetic profiles occurred on all non-native continents, either preceded by demographically independent founder events from various single sources or involving admixture scenarios. A decisive primarily admixed Polynesian bridgehead population serially colonised the entire Australasian region and its secondarily admixed descendants successively mediated invasions into Africa and Europe. Conversely, captive populations from several continents traced back to a single North American origin and exhibit considerably reduced genetic diversity, although some farmed strains carry distinct genetic signatures. We highlight genetic footprints characteristic of progressing domestication due to increasing socio-economic importance of H. illucens, and ongoing introgression between domesticated strains globally traded for large-scale farming and wild populations in some regions.

CONCLUSIONS

We document the dynamic population genetic history of a cosmopolitan dipteran of South American origin shaped by striking geographic patterns. These reflect both ancient dispersal routes, and stochastic and heterogeneous anthropogenic introductions during the last century leading to pronounced diversification of worldwide structure of H. illucens. Upon the recent advent of its agronomic commercialisation, however, current human-mediated translocations of the black soldier fly largely involve genetically highly uniform domesticated strains, which meanwhile threaten the genetic integrity of differentiated unique local resources through introgression. Our in-depth reconstruction of the contemporary and historical demographic trajectories of H. illucens emphasises benchmarking potential for applied future research on this emerging model of the prospering insect-livestock sector.

Isolation by geographical distance after release from Pleistocene refugia explains genetic and phenotypic variation in Xylotrupes siamensis (Coleoptera: Scarabaeidae)

Consistent and objective species delimitation is crucial to biodiversity studies, but challenges remain when conflicting taxonomic decisions have been made because different data sets and analytical methods were used to delineate species. In the rhinoceros beetle, Xylotrupes siamensis, the use of different morphological characters has resulted in taxonomic disagreement between studies. We used three molecular loci (mitochondrial CO1 and nuclear ITS2 and H3) to investigate the genetic divergence between populations exhibiting different male horn phenotypes. We also applied an approximate Bayesian computation approach to test alternative historical hypotheses that might explain the present genetic diversity among geographical populations. Furthermore, we used species distribution models to estimate the temporal variation in the geographical distribution of suitable habitats. The results show that the two phenotypic taxa within X. siamensis are not genetically structured and that their genetic structure can be explained using isolation by geographical distance. The emergence of the two phenotypic taxa might have been associated with historical isolation in separate refugia. However, spatial expansion and genetic interchange between populations might have gradually eroded the spatial genetic structure. We demonstrate that understanding the historical processes responsible for phenotypic divergence and genetic diversity among current populations could help with making evolutionarily coherent taxonomic decisions.

The Asian plethodontid salamander preserves historical genetic imprints of recent northern expansion

The Korean Peninsula, located at the southern tip of Northeast Asia, has never been covered by ice sheets and was a temperate refugium during the Pleistocene. Karsenia koreana, the sole Asian plethodontid salamander species, occurs only on the southern half of the Korean Peninsula and is thought to have found various climatic refugia. Despite its phylogenetic and biogeographic importance, no population-level genetic analysis has been performed on this species. Here we study the population genetic structure of K. koreana using mitochondrial and microsatellite loci to understand the recent historical dispersion process that shaped its current distribution. Overall, the genetic distance between populations correlated well with the spatial distance, and the genetic structure among populations showed signs of a unilateral northward expansion from a southernmost refugium population. Given the distinct genetic structure formed among the populations, the level of historical gene flow among populations appears to have been very low. As the estimated effective population size of K. koreana was also small, these results suggest that the small, restricted populations of K. koreana are extremely vulnerable to environmental changes that may require high levels of genetic diversity to cope with. Thus, special management strategies are needed to preserve these remnant populations.

Late Pleistocene Expansion of Small Murid Rodents across the Palearctic in Relation to the Past Environmental Changes

We investigated the evolutionary history of the striped field mouse to identify factors that initiated its past demographic changes and to shed light on the causes of its current genetic structure and trans-Eurasian distribution. We sequenced mitochondrial cyt b from 184 individuals, obtained from 35 sites in central Europe and eastern Mongolia. We compared genetic analyses with previously published historical distribution models and data on environmental and climatic changes. The past demographic changes displayed similar population trends in the case of recently expanded clades C1 and C3, with the glacial (MIS 3–4) expansion and postglacial bottleneck preceding the recent expansion initiated in the late Holocene and were related to environmental changes during the upper Pleistocene and Holocene. The past demographic trends of the eastern Asian clade C3 were correlated with changes in sea level and the formation of new land bridges formed by the exposed sea shelf during the glaciations. These data were supported by reconstructed historical distribution models. The results of our genetic analyses, supported by the reconstruction of the historical spatial distributions of the distinct clades, confirm that over time the local populations mixed as a consequence of environmental and climatic changes resulting from cyclical glaciation and the interglacial period during the Pleistocene.

Ring distribution patterns-diversification or speciation? Comparative phylogeography of two small mammals in the mountains surrounding the Sichuan Basin

Studying the genetic differentiation in a unique geographical area contributes to understanding the process of speciation. Here, we explore the spatial genetic structure and underlying formation mechanism of two congeneric small mammal species (Apodemus draco and A. chevrieri), which are mainly distributed in the mountains surrounding the lowland Sichuan Basin, southwest China. We applied a set of comparative phylogeographical analyses to determine their genetic diversification patterns, combining mitochondrial (Cytb and COI) and nuclear (microsatellite loci) markers, with dense sampling throughout the range (411 A. draco from 21 sites and 191 A. chevrieri from 22 sites). Moreover, we performed three complementary statistical methods to investigate the correlation between genotype and geographical and environmental components, and predicted the potential suitable distributional range under the present and historical climate conditions. Our results suggest that both species have experienced allopatric differentiation and admixture in historical periods, resulting in a ring-shape diversification, under the barrier effect of the Sichuan Basin. We infer that the tectonic events of the Qinghai-Tibetan Plateau and climatic oscillations during the Quaternary played an important role on the genetic divergence of the two species by providing environmental heterogeneity and geographical variation. Our study reveals a case of two sympatric small mammals following a ring-shaped diversification pattern and provides insight into the process of differentiation.

Estimating divergence times from DNA sequences

The patterns of genetic variation within and among individuals and populations can be used to make inferences about the evolutionary forces that generated those patterns. Numerous population genetic approaches have been developed in order to infer evolutionary history. Here, we present the "Two-Two (TT)" and the "Two-Two-outgroup (TTo)" methods; two closely related approaches for estimating divergence time based in coalescent theory. They rely on sequence data from two haploid genomes (or a single diploid individual) from each of two populations. Under a simple population-divergence model, we derive the probabilities of the possible sample configurations. These probabilities form a set of equations that can be solved to obtain estimates of the model parameters, including population split times, directly from the sequence data. This transparent and computationally efficient approach to infer population divergence time makes it possible to estimate time scaled in generations (assuming a mutation rate), and not as a compound parameter of genetic drift. Using simulations under a range of demographic scenarios, we show that the method is relatively robust to migration and that the TTo method can alleviate biases that can appear from drastic ancestral population size changes. We illustrate the utility of the approaches with some examples, including estimating split times for pairs of human populations as well as providing further evidence for the complex relationship among Neandertals and Denisovans and their ancestors.

Localized environmental heterogeneity drives the population differentiation of two endangered and endemic Opisthopappus Shih species

BACKGROUND

Climate heterogeneity not only indirectly shapes the genetic structures of plant populations, but also drives adaptive divergence by impacting demographic dynamics. The variable localized climates and topographic complexity of the Taihang Mountains make them a major natural boundary in Northern China that influences the divergence of organisms distributed across this region. Opisthopappus is an endemic genus of the Taihang Mountains that includes only two spatially partitioned species Opisthopappus longilobus and Opisthopappus taihangensis. For this study, the mechanisms behind the genetic variations in Opisthopappus populations were investigated.

RESULTS

Using SNP and InDel data coupled with geographic and climatic information, significant genetic differentiation was found to exist either between Opisthopappus populations or two species. All studied populations were divided into two genetic groups with the differentiation of haplotypes between the groups. At approximately 17.44 Ma of the early Miocene, O. taihangensis differentiated from O. longilobus under differing precipitation regimes due to the intensification of the Asian monsoon. Subsequently, intraspecific divergence might be induced by the dramatic climatic transformation from the mid- to late Miocene. During the Pleistocene period, the rapid uplift of the Taihang Mountains coupled with violent climatic oscillations would further promote the diversity of the two species. Following the development of the Taihang Mountains, its complex topography created geographical and ecological heterogeneity, which could lead to spatiotemporal isolation between the Opisthopappus populations. Thus the adaptive divergence might occur within these intraspecific populations in the localized heterogeneous environment of the Taihang Mountains.

CONCLUSIONS

The localized environmental events through the integration of small-scale spatial effects impacted the demographic history and differentiation mechanism of Opisthopappus species in the Taihang Mountains. The results provide useful information for us to understand the ecology and evolution of organisms in the mountainous environment from population and species perspective.

Allelic Diversity at Abiotic Stress Responsive Genes in Relationship to Ecological Drought Indices for Cultivated Tepary Bean, Phaseolus acutifolius A. Gray, and Its Wild Relatives

Some of the major impacts of climate change are expected in regions where drought stress is already an issue. Grain legumes are generally drought susceptible. However, tepary bean and its wild relatives within Phaseolus acutifolius or P. parvifolius are from arid areas between Mexico and the United States. Therefore, we hypothesize that these bean accessions have diversity signals indicative of adaptation to drought at key candidate genes such as: Asr2, Dreb2B, and ERECTA. By sequencing alleles of these genes and comparing to estimates of drought tolerance indices from climate data for the collection site of geo-referenced, tepary bean accessions, we determined the genotype x environmental association (GEA) of each gene. Diversity analysis found that cultivated and wild P. acutifolius were intermingled with var. tenuifolius and P. parvifolius, signifying that allele diversity was ample in the wild and cultivated clade over a broad sense (sensu lato) evaluation. Genes Dreb2B and ERECTA harbored signatures of directional selection, represented by six SNPs correlated with the environmental drought indices. This suggests that wild tepary bean is a reservoir of novel alleles at genes for drought tolerance, as expected for a species that originated in arid environments. Our study corroborated that candidate gene approach was effective for marker validation across a broad genetic base of wild tepary accessions.

Genetic diversity and population structure in Vitis species illustrate phylogeographic patterns in eastern North America

Geographical distribution and diversity of current plant species have been strongly shaped by climatic oscillations during the Quaternary. Analysing the resulting divergence among species and differentiation within species is crucial to understand the evolution of taxa like the Vitis genus, which provides very useful genetic resources for grapevine improvement and might reveal original recolonization patterns due to growth habit and dispersal mode. Here, we studied the genetic structure in natural populations of three species from eastern North America: Vitis aestivalis, V. cinerea and V. riparia using different marker types. Vitis aestivalis and V. cinerea showed higher diversity than V. riparia. The two former species are less differentiated, confirming an earlier divergence of V. riparia. V. aestivalis and V. riparia exhibited different genetic groups on both sides of the Appalachian Mountains that could mirror different recolonization routes from southern refugia. Genetic structure was stronger in V. cinerea, for which two varieties (var. berlandieri and var. cinerea) are morphologically recognized. Our results confirm this distinction and suggest the existence of three other lineages within var. cinerea. These discontinuities appear linked to adaptation of var. berlandieri to dry and limy areas of Texas and partially to the Mississippi River Valley. Rapid range expansions from refugia upon climate warming are also suggested by the low linkage disequilibrium values observed. Furthermore, large variation for downy mildew resistance was observed in the three species. Our findings appeared consistent with the vegetation history of eastern North America.

A Whole-Genome Scan for Association with Invasion Success in the Fruit Fly Drosophila suzukii Using Contrasts of Allele Frequencies Corrected for Population Structure

Evidence is accumulating that evolutionary changes are not only common during biological invasions but may also contribute directly to invasion success. The genomic basis of such changes is still largely unexplored. Yet, understanding the genomic response to invasion may help to predict the conditions under which invasiveness can be enhanced or suppressed. Here, we characterized the genome response of the spotted wing drosophila Drosophila suzukii during the worldwide invasion of this pest insect species, by conducting a genome-wide association study to identify genes involved in adaptive processes during invasion. Genomic data from 22 population samples were analyzed to detect genetic variants associated with the status (invasive versus native) of the sampled populations based on a newly developed statistic, we called C2, that contrasts allele frequencies corrected for population structure. We evaluated this new statistical framework using simulated data sets and implemented it in an upgraded version of the program BayPass. We identified a relatively small set of single-nucleotide polymorphisms that show a highly significant association with the invasive status of D. suzukii populations. In particular, two genes, RhoGEF64C and cpo, contained single-nucleotide polymorphisms significantly associated with the invasive status in the two separate main invasion routes of D. suzukii. Our methodological approaches can be applied to any other invasive species, and more generally to any evolutionary model for species characterized by nonequilibrium demographic conditions for which binary covariables of interest can be defined at the population level.

Algerian Olive Germplasm and Its Relationships with the Central-Western Mediterranean Varieties Contributes to Clarify Cultivated Olive Diversification

Olive tree with its main final product, olive oil, is an important element of Mediterranean history, considered the emblematic fruit of a civilization. Despite its wide diffusion and economic and cultural importance, its evolutionary and phylogenetic history is still difficult to clarify. As part of the Mediterranean basin, Algeria was indicated as a secondary diversification center. However, genetic characterization studies from Maghreb area, are currently underrepresented. In this context, we characterized 119 endemic Algerian accessions by using 12 microsatellite markers with the main goal to evaluate the genetic diversity and population structure. In order to provide new insights about the history of olive diversification events in the Central-Western Mediterranean basin, we included and analyzed a sample of 103 Italian accessions from Sicily and, a set of molecular profiles of cultivars from the Central-Western Mediterranean area. The phylogenetic investigation let us to evaluate genetic relationships among Central-Mediterranean basin olive germplasm, highlight new synonymy cases to support the importance of vegetative propagation in the cultivated olive diffusion and consolidate the hypothesis of more recent admixture events occurrence. This work provided new information about Algerian germplasm biodiversity and contributed to clarify olive diversification process.

Europe as a bridgehead in the worldwide invasion history of grapevine downy mildew, Plasmopara viticola

Europe is the historical cradle of viticulture, but grapevines (Vitis vinifera) have been increasingly threatened by pathogens of American origin. The invasive oomycete Plasmopara viticola causes downy mildew, one of the most devastating grapevine diseases worldwide. Despite major economic consequences, its invasion history remains poorly understood. We analyzed a comprehensive dataset of ∼2,000 samples, collected from the most important wine-producing countries, using nuclear and mitochondrial gene sequences and microsatellite markers. Population genetic analyses revealed very low genetic diversity in invasive downy mildew populations worldwide and little evidence of admixture. All the invasive populations originated from only one of the five native North American lineages, the one parasitizing wild summer grape (V. aestivalis). An approximate Bayesian computation-random forest approach allowed inferring the worldwide invasion scenario of P. viticola. After an initial introduction into Europe, invasive European populations served as a secondary source of introduction into vineyards worldwide, including China, South Africa, and twice independently, Australia. Only the invasion of Argentina probably represents a tertiary introduction, from Australia. Our findings provide a striking example of a global pathogen invasion resulting from secondary dispersal of a successful invasive population. Our study will also help designing quarantine regulations and efficient breeding for resistance against grapevine downy mildew.