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

Population structure and genetic diversity of the endangered fish black shinner Pseudopungtungia nigra (Cyprinidae) in Korea: a wild and restoration population

The black shinner Pseudopungtungia nigra Mori, 1935 is an endangered fish endemic to Korea. It lives in the narrow basin of the Geumgang River, Mangyeonggang River, and Ungcheoncheon Stream, which flow into the West Sea of Korea. One population of P. nigra in Ungcheoncheon Stream has been locally exterminated once; it is now inhabiting the upper reaches of the dam through a restoration program. Efforts to identify and understand the genetic structure of these populations are important for conservation planning. Here, we analyzed genetic diversity using 21 microsatellite markers for 9 populations. The mean number of alleles ranged from 4.4 to 8.1, mean allelic richness ranged from 4.6 to 7.8, mean observed heterozygosity ranged from 0.519 to 0.702, and mean expected heterozygosity ranged from 0.540 to 0.763. All groups had recent and historical bottlenecks (P < 0.05, M-ratio < 0.68). Three groups [YD (2019), OC and UC] had significant inbreeding index values, suggesting that they were engaged in inbreeding. We observed a moderate level of genetic differentiation between MG and the rest of the population (FST = 0.135 to 0.168, P < 0.05). The genetic structure exhibited a fitting constant K = 2, along with separation between MG and the remaining populations. With respect to genetic flow, YD (2019), OC, CG, and ND shifted to the UC population (0.263 to 0.278). The genetic flow of each population was transferred only within the population; there was no gene flow among populations, except for the Ungcheoncheon Stream population. This study shows that the Ungcheoncheon Stream population needs conservation efforts to increase its genetic diversity, and the Geumgang River populations needs a conservation plan that considers the possibility of conservation and evolution through gene exchange among the populations.

Climate heterogeneity shapes phylogeographic pattern of Hippophae gyantsensis (Elaeagnaceae) in the east Himalaya-Hengduan Mountains

The interaction of recent orographic uplift and climate heterogeneity acted as a key role in the East Himalaya-Hengduan Mountains (EHHM) has been reported in many studies. However, how exactly the interaction promotes clade diversification remains poorly understood. In this study, we both used the chloroplast trnT-trnF region and 11 nuclear microsatellite loci to investigate the phylogeographic structure and population dynamics of Hippophae gyantsensis and estimate what role geological barriers or ecological factors play in the spatial genetic structure. The results showed that this species had a strong east-west phylogeographic structure, with several mixed populations identified from microsatellite data in central location. The intraspecies divergence time was estimated to be about 3.59 Ma, corresponding well with the recent uplift of the Tibetan Plateau. Between the two lineages, there was significant climatic differentiation without geographic barriers. High consistency between lineage divergence, climatic heterogeneity, and Qingzang Movement demonstrated that climatic heterogeneity but not geographic isolation drives the divergence of H. gyantsensis, and the recent regional uplift of the QTP, as the Himalayas, creates heterogeneous climates by affecting the flow of the Indian monsoon. The east group of H. gyantsensis experienced population expansion c. 0.12 Ma, closely associated with the last interglacial interval. Subsequently, a genetic admixture event between east and west groups happened at 26.90 ka, a period corresponding to the warm inter-glaciation again. These findings highlight the importance of the Quaternary climatic fluctuations in the recent evolutionary history of H. gyantsensis. Our study will improve the understanding of the history and mechanisms of biodiversity accumulation in the EHHM region.

Geographic isolation and climatic heterogeneity drive population differentiation of Rosa chinensis var. spontanea complex

Global biodiversity is contracting rapidly due to potent anthropogenic activities and severe climate change. Wild populations of Rosa chinensis var. spontanea and Rosa lucidissima are rare species endemic to China, as well as important germplasm resources for rose breeding. However, these populations are at acute risk of extinction and require urgent action to ensure their preservation. We harnessed 16 microsatellite loci to 44 populations of these species and analysed population structure and differentiation, demographic history, gene flow and barrier effect. In addition, a niche overlap test and potential distribution modelling in different time periods were also carried out. The data indicate that: (1) R. lucidissima cannot be regarded as a separate species from R. chinensis var. spontanea; (2) the Yangtze River and the Wujiang River function as barriers in population structure and differentiation, and precipitation in the coldest quarter may be the key factor for niche divergence of R. chinensis var. spontanea complex; (3) historical gene flow showed a converse tendency to current gene flow, indicating that alternate migration events of R. chinensis var. spontanea complex between south and north were a response to climate oscillations; and (4) extreme climate change will decrease the distribution range of R. chinensis var. spontanea complex, whereas the opposite will occur under a moderate scenario for the future. Our results resolve the relationship between R. chinensis var. spontanea and R. lucidissima, highlight the pivotal roles of geographic isolation and climate heterogeneity in their population differentiation, and provide an important reference for comparable conservation studies on other endangered species.

Phylogeography of the temperate grassland plant Tephroseris kirilowii (Asteraceae) inferred from multiplexed inter-simple sequence repeat genotyping by sequencing (MIG-seq) data

A group of temperate grassland plant species termed the "Mansen elements" occurs in Japan and is widely distributed in the grasslands of continental East Asia. It has been hypothesized that these species are continental grassland relicts in Japan that stretch back to a colder age, but their migration history has not been elucidated. To assess the migration history of the Mansen elements, we performed phylogeographic analyses of Tephroseris kirilowii, a member of this group, using single-nucleotide polymorphisms (SNPs) obtained from multiplexed inter-simple sequence repeat genotyping by sequencing (MIG-seq). It was estimated that the Japanese populations of T. kirilowii were divided from those of continental East Asia at 25.2 thousand years ago (ka) with 95% highest probability density interval (HPD) of 15.3-40.0 ka and that Japanese clades first diverged at 20.2 ka with 95% HPD of 10.4-30.1 ka. As the climatically suitable range during the last glacial maximum (LGM) estimated using ecological niche modeling (ENM) was limited in Japan and there was a slight genetic differentiation among Japanese populations, a post-glacial expansion of T. kirilowii in the Japanese Archipelago was indicated.

Genomic evidence for homoploid hybrid speciation in a marine mammal apex predator

Hybridization is widespread and constitutes an important source of genetic variability and evolution. In animals, its role in generating novel and independent lineages (hybrid speciation) has been strongly debated, with only a few cases supported by genomic data. The South American fur seal (SAfs) Arctocephalus australis is a marine apex predator of Pacific and Atlantic waters, with a disjunct set of populations in Peru and Northern Chile [Peruvian fur seal (Pfs)] with controversial taxonomic status. We demonstrate, using complete genome and reduced representation sequencing, that the Pfs is a genetically distinct species with an admixed genome that originated from hybridization between the SAfs and the Galapagos fur seal (Arctocephalus galapagoensis) ~400,000 years ago. Our results strongly support the origin of Pfs by homoploid hybrid speciation over alternative introgression scenarios. This study highlights the role of hybridization in promoting species-level biodiversity in large vertebrates.

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes)

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide the most comprehensive phylogeographic framework for the A. obstetricans complex to date, as well as the first fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

Post-glacial colonization of the Fennoscandian coast by a plant parasitic insect with an unusual life history

Species that exhibit very peculiar ecological traits combined with limited dispersal ability pose a challenge to our understanding of ecological and evolutionary mechanisms. This is especially true when they have managed to spread over long distances, overcome physical barriers, and colonize large areas. Climate and landscape changes, trophic web relations, as well as life history all interact to shape migration routes and present-day species distributions and their population genetic structures. Here we analyzed the post-glacial colonization of northern Europe by the gall midge Contarinia vincetoxici, which is a monophagous parasite on the perennial herb White swallowwort (Vincetoxicum hirundinaria). This insect not only has a narrow feeding niche but also limited dispersal ability and an exceptionally long dormancy. Gall midge larvae (n = 329) were collected from 16 sites along its distribution range in Denmark, Sweden, and Finland. Using microsatellite loci and knowledge of the species and the regions' history, we investigated the role of landscape change, host plant distribution, insect population dynamics, and life history in shaping the population genetic structure of the insect. We devoted particular interest to the role of the insect's presumed poor dispersal capacity in combination with its exceptionally extended diapause. We found significant levels of local inbreeding (95% highest posterior density interval = 0.42-0.47), low-level within-population heterozygosity (mean H E = 0.45, range 0.20-0.61) with private alleles in all populations except two. We also found significant (p < .001) regional isolation-by-distance patterns, suggesting regularly recurring mainly short-distance dispersal. According to approximate Bayesian computations, C. vincetoxici appears to have colonized the study area via wind-aided flights from remote areas approximately 4600-700 years before present when the land has gradually risen above the sea level. Extremely long dormancy periods have allowed the species to "disperse in time", thereby aiding population persistence despite generally low census population sizes.

Climatic oscillation promoted diversification of spinous assassin bugs during Pleistocene glaciation

Insect speciation is among the most fascinating topics in evolutionary biology; however, its underlying mechanisms remain unclear. Allopatric speciation represents one of the major types of speciation and is believed to have frequently occurred during glaciation periods, when climatic oscillation may have caused suitable habitats to be fragmented repeatedly, creating geographical isolation among populations. However, supporting evidence for allopatric speciation of insects in East Asia during the Pleistocene glaciation remains lacking. We aim to investigate the effect of climatic oscillation during the Pleistocene glaciation on the diversification pattern and evolutionary history of hemipteran insects and to test the hypothesis of Pleistocene species stability using spinous assassin bugs Sclomina (Hemiptera: Reduviidae), a small genus widely distributed in southern China but was later found to have cryptic species diversity. Here, using the whole mitochondrial genome (mitogenome) and nuclear ribosomal RNA genes, we investigated both interspecific and intraspecific diversification patterns of spinous assassin bugs. Approximate Bayesian computation, ecological niche modeling, and demographic history analyses were also applied to understand the diversification process and driven factors. Our data suggest that the five species of Sclomina are highly diverged, despite three of them currently being cryptic. Speciation occurred during the Pleistocene when suitable distribution areas were possibly fragmented. Six phylogeographic groups in the type species S. erinacea were identified, among which two groups underwent expansion during the early Last Glacial Period and after Last Glacier Maximum. Our analyses suggest that this genus may have experienced climate-driven habitat fragmentation and postglacial expansion in the Pleistocene, promoting allopatric speciation and intraspecific diversification. Our results reveal underestimated species diversity in a small insect group and illustrate a remarkable example of allopatric speciation of insects in East Asia promoted by Pleistocene climatic oscillations. These findings provide important insights into the speciation processes and aid the conservation of insect species diversity.

Genetic and ecophysiological evidence that hybridization facilitated lineage diversification in yellow Camellia (Theaceae) species: a case study of natural hybridization between C. micrantha and C. flavida

BACKGROUND

Hybridization is generally considered an important creative evolutionary force, yet this evolutionary process is still poorly characterized in karst plants. In this study, we focus on natural hybridization in yellow Camellia species, a group of habitat specialists confined to karst/non-karst habitats in southwestern China.

RESULTS

Based on population genome data obtain from double digest restriction-site associated DNA (ddRAD) sequencing, we found evidence for natural hybridization and introgression between C. micrantha and C. flavida, and specifically confirmed their hybrid population, C. "ptilosperma". Ecophysiological results suggested that extreme hydraulic traits were fixed in C. "ptilosperma", these being consistent with its distinct ecological niche, which lies outside its parental ranges.

CONCLUSION

The identified hybridization event is expected to have played a role in generating novel variation during, in which the hybrid population displays different phenological characteristics and novel ecophysiological traits associated with the colonization of a new niche in limestone karst.

Phylogenetic divergence associated with climate oscillations and topology illustrates the dispersal history of Formosan sambar deer (Rusa unicolor swinhoii) in Taiwan

The island of Taiwan represents an ideal context for studying the effects of climatic oscillations and topographic variation on large herbivores due to its varied tropical to sub-tropical climate zones at different elevational ranges. We explored the phylogenetics of Formosan sambar deer (Rusa unicolor swinhoii) using the control region of the mitochondrial genome. We detected 18 haplotypes among 454 sequences across the island and grouped them into six regions based on SAMOVA, with 68.78% variance among regions. A Bayesian phylogenetic dendrogram revealed two spatially segregated genetic clades. Neutrality tests and Bayesian skyline plots uncovered different demographic expansion histories for the two clades. We further tested divergence times and chronology to propose potential phylogenetic scenarios, which were examined using approximate Bayesian computation. Finally, we present a credible hypothesis for a glacial refugium in the northern part of the Central Mountain Range. Subsequent secondary contact between the two clades during interglacial periods has led to the extant genetic structure of Formosan sambar deer.

Historical Landscape Evolution Shaped the Phylogeography and Population History of the Cyprinid Fishes of Acrossocheilus (Cypriniformes: Cyprinidae) According to Mitochondrial DNA in Zhejiang Province, China

Geological events and landscape features, as well as changes in the climate during the Pliocene period, have shaped the distribution of genetic diversity and demographic history of freshwater fish in mainland China. In this study, we investigated the phylogeny and population genetic structure of Acrossocheilus species (A. fasciatus, A. kreyenbergii and A. wenchowensis) in the Zhejiang region by the mitochondrial cytochrome b (cyt-b) and control (D-loop) region s. Mitochondrial phylogenetic analysis revealed three major lineages (lineages A, B and C), which may represent A. fasciatus, A. wenchowensis and A. kreyenbergii, respectively. Our results revealed that A. fasciatus and A. wenchowensis diverged from A. kreyenbergii in the Zhejiang region. The uplift of the Wuyi Mountains and the Xianxia Mountains served as an important geographic barrier in the diversification of the three Acrossocheilus species. The most recent common ancestors (TMRCAs) of the three lineages and lineages A + B were dated to 1.859 and 1.614 myr, respectively. Our results indicate that the effective population sizes of A. fasciatus and A. wenchowensis remained constant from the past to the present, as supported by ABC analysis. For conservation and protection, a strategy is required because of their genetic uniqueness, and we suggest that the two regions divided by the Xianxia Mountains be regarded as different management units (Mus), conforming to the major zoological regions of the Zhejiang region.

Phylogenetic History and Phylogeographic Patterns of the European Wildcat (Felis silvestris) Populations

Disentangling phylogenetic and phylogeographic patterns is fundamental to reconstruct the evolutionary histories of taxa and assess their actual conservation status. Therefore, in this study, for the first time, the most exhaustive biogeographic history of European wildcat (Felis silvestris) populations was reconstructed by typing 430 European wildcats, 213 domestic cats, and 72 putative admixed individuals, collected across the entire species' distribution range, at a highly diagnostic portion of the mitochondrial ND5 gene. Phylogenetic and phylogeographic analyses identified two main ND5 lineages (D and W) roughly associated with domestic and wild polymorphisms. Lineage D included all domestic cats, 83.3% of putative admixed individuals, and also 41.4% of wildcats; these latter mostly showed haplotypes belonging to sub-clade Ia, that diverged about 37,700 years ago, long pre-dating any evidence for cat domestication. Lineage W included all the remaining wildcats and putative admixed individuals, spatially clustered into four main geographic groups, which started to diverge about 64,200 years ago, corresponding to (i) the isolated Scottish population, (ii) the Iberian population, (iii) a South-Eastern European cluster, and (iv) a Central European cluster. Our results suggest that the last Pleistocene glacial isolation and subsequent re-expansion from Mediterranean and extra-Mediterranean glacial refugia were pivotal drivers in shaping the extant European wildcat phylogenetic and phylogeographic patterns, which were further modeled by both historical natural gene flow among wild lineages and more recent wild x domestic anthropogenic hybridization, as confirmed by the finding of F. catus/lybica shared haplotypes. The reconstructed evolutionary histories and the wild ancestry contents detected in this study could be used to identify adequate Conservation Units within European wildcat populations and help to design appropriate long-term management actions.

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.

Wind-dispersed seeds blur phylogeographic breaks: The complex evolutionary history of Populus lasiocarpa around the Sichuan Basin

The strength of phylogeographic breaks can vary among species in the same area despite being subject to the same geological and climate history due to differences in biological traits. Several important phylogeographic breaks exist around the Sichuan Basin in Southwest China but few studies have focused on wind-dispersed plants. Here, we investigated the phylogeographic patterns and the evolutionary history of Populus lasiocarpa, a wind-pollinated and wind-dispersed tree species with a circum-Sichuan Basin distribution in southwest China. We sequenced and analyzed three plastid DNA fragments (ptDNA) and eight nuclear microsatellites (nSSRs) of 265 individuals of P. lasiocarpa from 21 populations spanning the entire distribution range. Distribution patterns based on nSSR data revealed that there are three genetic groups in P. lasiocarpa. This is consistent with the three phylogeographic breaks (Sichuan Basin, the Kaiyong Line and the 105°E line), where the Sichuan basin acts as the main barrier to gene flow between western and eastern groups. However, the distribution pattern based on ptDNA haplotypes poorly matched the phylogeographic breaks, and wind-dispersed seeds may be one of the main contributing factors. Species distribution modelling suggested a larger potential distribution in the last glacial maximum with a severe bottleneck during the last interglacial. A DIYABC model also suggested a population contraction and expansion for both western and eastern lineages. These results indicate that biological traits are likely to affect the evolutionary history of plants, and that nuclear molecular markers, which experience higher levels of gene flow, might be better indicators of phylogeographic breaks.

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

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

Exploring island syndromes: Variable matrix permeability in Phalaenopsis pulcherrima (Orchidaceae), a specialist lithophyte of tropical Asian inselbergs

INTRODUCTION

Plants confined to island-like habitats are hypothesised to possess a suite of functional traits that promote on-spot persistence and recruitment, but this may come at the cost of broad-based colonising potential. Ecological functions that define this island syndrome are expected to generate a characteristic genetic signature. Here we examine genetic structuring in the orchid Phalaenopsis pulcherrima, a specialist lithophyte of tropical Asian inselbergs, both at the scale of individual outcrops and across much of its range in Indochina and on Hainan Island, to infer patterns of gene flow in the context of an exploration of island syndrome traits.

METHODS

We sampled 323 individuals occurring in 20 populations on 15 widely scattered inselbergs, and quantified genetic diversity, isolation-by-distance and genetic structuring using 14 microsatellite markers. To incorporate a temporal dimension, we inferred historical demography and estimated direction of gene flow using Bayesian approaches.

RESULTS

We uncovered high genotypic diversity, high heterozygosity and low rates of inbreeding, as well as strong evidence for the occurrence of two genetic clusters, one comprising the populations of Hainan Island and the other those of mainland Indochina. Connectivity was greater within, rather than between the two clusters, with the former unequivocally supported as ancestral.

DISCUSSION

Despite a strong capacity for on-spot persistence conferred by clonality, incomplete self-sterility and an ability to utilize multiple magnet species for pollination, our data reveal that P. pulcherrima also possesses traits that promote landscape-scale gene flow, including deceptive pollination and wind-borne seed dispersal, generating an ecological profile that neither fully conforms to, nor fully contradicts, a putative island syndrome. A terrestrial matrix is shown to be significantly more permeable than open water, with the direction of historic gene flow indicating that island populations can serve as refugia for postglacial colonisation of continental landmasses by effective dispersers.

A historic religious sanctuary may have preserved ancestral genetics of Japanese sika deer (Cervus nippon)

Deer have been a major resource for human populations for thousands of years. Anthropogenic activities, such as hunting, have influenced the genetic structure and distribution of deer populations. In Japan, wild Japanese sika deer (Cervus nippon) have been hunted since ancient times but have also been historically protected as sacred animals in several sanctuaries. Sika deer have been protected for over a thousand years in the religious sanctuary around the Kasuga Taisha Shrine on the Kii Peninsula, located in the center of Japan. Here, we used short sequence repeats (SSR) and mitochondrial DNA (mtDNA) to investigate the genetic diversity, population structure, and demography of Japanese sika deer inhabiting the Kii Peninsula, Japan, and discuss possible anthropogenic influences. Using SSR, three distinct genetic groups were distinguished on the Kii Peninsula: an Eastern genetic group, a Western genetic group, and an isolated genetic group with individuals in the religious sanctuary of Kasuga Taisha Shrine in Nara city. The isolated genetic sanctuary group had only the mtDNA haplotype S4. The SSR genotype data suggested a newer divergence time of the genetic groups of the religious sanctuary than would have occurred as a result of Late Quaternary climate change. This time scale coincided with the establishment of the sanctuary with Kasuga Taisha Shrine. Thus, the religious protection conserved genetic variation over a thousand years.

Suspicions of two bridgehead invasions of Xylella fastidiosa subsp. multiplex in France

Of American origin, a wide diversity of Xylella fastidiosa strains belonging to different subspecies have been reported in Europe since 2013 and its discovery in Italian olive groves. Strains from the subspecies multiplex (ST6 and ST7) were first identified in France in 2015 in urban and natural areas. To trace back the most probable scenario of introduction in France, the molecular evolution rate of this subspecies was estimated at 3.2165 × 10^-7 substitutions per site per year, based on heterochronous genome sequences collected worldwide. This rate allowed the dating of the divergence between French and American strains in 1987 for ST6 and in 1971 for ST7. The development of a new VNTR-13 scheme allowed tracing the spread of the bacterium in France, hypothesizing an American origin. Our results suggest that both sequence types were initially introduced and spread in Provence-Alpes-Côte d'Azur (PACA); then they were introduced in Corsica in two waves from the PACA bridgehead populations.

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.

Population Genetics of Oncomelania hupensis Snails from New-Emerging Snail Habitats in a Currently Schistosoma japonicum Non-Endemic Area

Schistosomiasis is still one of the most significant neglected tropical diseases worldwide, and China is endemic for Schistosoma japonicum. With its great achievement in schistosomiasis control, the government of China has set the goal to eliminate the parasitic disease at the country level by 2030. However, one major challenge is the remaining huge areas of habitats for the intermediate host Oncomelania hupensis. This is further exacerbated by an increasing number of new emerging snail habitats reported each year. Therefore, population genetics on snails in such areas will be useful in evaluation of snail control effect and/or dispersal. We then sampled snails from new emerging habitats in Taicang of Jiangsu, China, a currently S. japonicum non-endemic area from 2014 to 2017, and performed population genetic analyses based on nine microsatellites. Results showed that all snail populations had low genetic diversity, and most genetic variations originated from within snail populations. The estimated effective population size for the 2015 population was infinitive. All snails could be separated into two clusters, and further DIYABC analysis revealed that both the 2016 and the 2017 populations may derive from the 2015, indicating that the 2017 population must have been missed in the field survey performed in 2016. These findings may have implications in development of more practical guidelines for snail monitoring and control.

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.

GADMA2: more efficient and flexible demographic inference from genetic data

BACKGROUND

Inference of complex demographic histories is a source of information about events that happened in the past of studied populations. Existing methods for demographic inference typically require input from the researcher in the form of a parameterized model. With an increased variety of methods and tools, each with its own interface, the model specification becomes tedious and error-prone. Moreover, optimization algorithms used to find model parameters sometimes turn out to be inefficient, for instance, by being not properly tuned or highly dependent on a user-provided initialization. The open-source software GADMA addresses these problems, providing automatic demographic inference. It proposes a common interface for several likelihood engines and provides global parameters optimization based on a genetic algorithm.

RESULTS

Here, we introduce the new GADMA2 software and provide a detailed description of the added and expanded features. It has a renovated core code base, new likelihood engines, an updated optimization algorithm, and a flexible setup for automatic model construction. We provide a full overview of GADMA2 enhancements, compare the performance of supported likelihood engines on simulated data, and demonstrate an example of GADMA2 usage on 2 empirical datasets.

CONCLUSIONS

We demonstrate the better performance of a genetic algorithm in GADMA2 by comparing it to the initial version and other existing optimization approaches. Our experiments on simulated data indicate that GADMA2's likelihood engines are able to provide accurate estimations of demographic parameters even for misspecified models. We improve model parameters for 2 empirical datasets of inbred species.

Genetic Variability and Admixture Zones in the Italian Populations of Turkey Oak (Quercus cerris L.)

The Turkey oak (Quercus cerris L.) is widely distributed in Italy, where it is the ecologically dominant oak on sandy and acidic soil. In this work, we analysed 23 natural populations by means of eight SSR (microsatellite) markers, to obtain the first synthetic map of genetic variability for this species and to study its dispersion during the Holocene, due to the possibility that at least one refugium during the Last Glacial Maximum was in Italy. The analyses showed a good amount of genetic variability together with fair differentiation between populations, as indicated by F_ST = 0.059. A Bayesian analysis of the amount of admixture among populations revealed the presence of four putative gene pools of origin and a rough subdivision of the populations according to their geographic location, as confirmed by the spatial analysis. No evidence for the existence of putative refugial populations was found; however, this study paves the way for the planning of conservation strategies also with regard to the relationship between Turkey oak and other oak species in Italy.

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.

New insights on patterns of genetic admixture and phylogeographic history in Iberian high mountain populations of midwife toads

Multiple Quaternary glacial refugia in the Iberian Peninsula, commonly known as "refugia within refugia", allowed diverging populations to come into contact and admix, potentially boosting substantial mito-nuclear discordances. In this study, we employ a comprehensive set of mitochondrial and nuclear markers to shed light onto the drivers of geographical differentiation in Iberian high mountain populations of the midwife toads Alytes obstetricans and A. almogavarii from the Pyrenees, Picos de Europa and Guadarrama Mountains. In the three analysed mountain regions, we detected evidence of extensive mito-nuclear discordances and/or admixture between taxa. Clustering analyses identified three major divergent lineages in the Pyrenees (corresponding to the eastern, central and central-western Pyrenees), which possibly recurrently expanded and admixed during the succession of glacial-interglacial periods that characterised the Late Pleistocene, and that currently follow a ring-shaped diversification pattern. On the other hand, populations from the Picos de Europa mountains (NW Iberian Peninsula) showed a mitochondrial affinity to central-western Pyrenean populations and a nuclear affinity to populations from the central Iberian Peninsula, suggesting a likely admixed origin for Picos de Europa populations. Finally, populations from the Guadarrama Mountain Range (central Iberian Peninsula) were depleted of genetic diversity, possibly as a consequence of a recent epidemic of chytridiomycosis. This work highlights the complex evolutionary history that shaped the current genetic composition of high mountain populations, and underscores the importance of using a multilocus approach to better infer the dynamics of population divergence.

Clear phylogeographical structures shed light on the origin and dispersal of the aquatic boreal plant Hippuris vulgaris

Aquatic plants are an important ecological group in the arctic flora; however, their evolutionary histories remain largely unknown. In order to deepen our understanding of the evolution of these plants, we explored the phylogeographical structure of an aquatic boreal plant Hippuris vulgaris in a broad geographical sampling from Eurasia and North America using the chloroplast intergenic spacer psbA-trnH and seven nuclear microsatellite loci. Two closely-related species H. lanceolata and H. tetraphylla were also included because of their taxonomic controversy. Both chloroplast DNA sequences and nuclear microsatellite data revealed three genetic lineages with distinct distribution ranges. Incongruence between nuclear and chloroplast DNA lineages occurred in 14 samples from Russian Far East and Europe caused by inter-lineage hybridization. No private haplotypes or independent genetic clusters were evident in H. lanceolata or H. tetraphylla, suggesting that these two species should be considered conspecific ecotypes of H. vulgaris. Analysis using Approximate Bayesian Computation-Random Forest approach suggests that Hippuris vulgaris originated in China, followed by dispersal into Russia plus Northeast China, then successively westwards into Europe and North America, and finally into the Russian Far East from both North America and Russia plus Northeast China. This study is the first to elucidate the historical dispersal processes of a circumarctic aquatic plant across the entirety of its range.

Genome-wide SNPs in the spiny lobster Panulirus homarus reveal a hybrid origin for its subspecies

Background

Evolutionary divergence and speciation often occur at a slower rate in the marine realm due to the higher potential for long-distance reproductive interaction through larval dispersal. One common evolutionary pattern in the Indo-Pacific, is divergence of populations and species at the peripheries of widely-distributed organisms. However, the evolutionary and demographic histories of such divergence are yet to be well understood. Here we address these issues by coupling genome-wide SNP data with mitochondrial DNA sequences to test the patterns of genetic divergence and possible secondary contact among geographically distant populations of the highly valuable spiny lobster Panulirus homarus species complex, distributed widely through the Indo-Pacific, from South Africa to the Marquesas Islands.

Result

After stringent filtering, 2020 SNPs were used for population genetic and demographic analyses, revealing strong regional structure (FST = 0.148, P < 0001), superficially in accordance with previous analyses. However, detailed demographic analyses supported a much more complex evolutionary history of these populations, including a hybrid origin of a North-West Indian Ocean (NWIO) population, which has previously been discriminated morphologically, but not genetically. The best-supported demographic models suggested that the current genetic relationships among populations were due to a complex series of past divergences followed by asymmetric migration in more recent times.

Conclusion

Overall, this study suggests that alternating periods of marine divergence and gene flow have driven the current genetic patterns observed in this lobster and may help explain the observed wider patterns of marine species diversity in the Indo-Pacific.

Diversity, migration routes, and worldwide population genetic structure of Lecanosticta acicola, the causal agent of brown spot needle blight

Lecanosticta acicola is a pine needle pathogen causing brown spot needle blight that results in premature needle shedding with considerable damage described in North America, Europe, and Asia. Microsatellite and mating type markers were used to study the population genetics, migration history, and reproduction mode of the pathogen, based on a collection of 650 isolates from 27 countries and 26 hosts across the range of L. acicola. The presence of L. acicola in Georgia was confirmed in this study. Migration analyses indicate there have been several introduction events from North America into Europe. However, some of the source populations still appear to remain unknown. The populations in Croatia and western Asia appear to originate from genetically similar populations in North America. Intercontinental movement of the pathogen was reflected in an identical haplotype occurring on two continents, in North America (Canada) and Europe (Germany). Several shared haplotypes between European populations further suggests more local pathogen movement between countries. Moreover, migration analyses indicate that the populations in northern Europe originate from more established populations in central Europe. Overall, the highest genetic diversity was observed in south-eastern USA. In Europe, the highest diversity was observed in France, where the presence of both known pathogen lineages was recorded. Less than half of the observed populations contained mating types in equal proportions. Although there is evidence of some sexual reproduction taking place, the pathogen spreads predominantly asexually and through anthropogenic activity.

Multilocus Phylogeography and Population Genetic Analyses of Opsariichthys hainanensis Reveal Pleistocene Isolation Followed by High Gene Flow around the Gulf of Tonkin

The ichthyofauna of continental islands is characterized by immigration through a land bridge due to fluctuating sea levels. Hainan Island is adjacent to the southern margin of mainland China and provides opportunities for understanding the origin and diversification of freshwater fishes. The aim of our study was to evaluate the level of genetic variation and phylogeographic structure of Opsariichthys hainanensis on Hainan Island and mainland China, using mtDNA cyt b gene (1140 bp) and D-loop (926 bp), nuclear RAG1 gene (1506 bp), and 12 microsatellite loci. Mitochondrial phylogenetic analysis identified five major lineages according to the geographical distribution from different populations. We suggested that two dispersal events occurred: the population in the Changhua River migrated to the Red River (Lineage B), and the populations in the South Hainan region moved northwards to the North Hainan region. However, populations in Northwest Hainan Island dispersed to the populations around the Gulf of Tonkin (Lineage A1) and populations in Northeast Hainan Island dispersed to the populations in mainland China (Lineage A2). Our results indicated that the populations of O. hainanensis suffered a bottleneck event followed by a recent population expansion supported by the ABC analysis. We suggest that O. hainanensis populations were found mostly in the lowlands and a lack of suitable freshwater habitat in southern mainland China and Hainan during the Last Interglacial period, and then expansion occurred during the Last Glacial Maximum.

The emergence of a cryptic lineage and cytonuclear discordance through past hybridization in the Japanese fire-bellied newt, Cynops pyrrhogaster (Amphibia: Urodela)

Discrepancies in geographic variation patterns between nuclear DNA and mitochondrial DNA (mtDNA) are the result of the complicated differentiation processes in organisms and are the key to understanding their true evolutionary processes. The genetic differentiation of the northern and Southern-Izu lineages of the Japanese newt, Cynops pyrrhogaster, was investigated through their single nucleotide polymorphism variations obtained via multiplexed ISSR genotyping by sequencing (MIG-seq). We found three genetic groups (Tohoku, N-Kanto and S-Kanto), that were not detected by mtDNA variations, in the northern lineage. N-Kanto has intermediate genetic characteristics between Tohoku and S-Kanto. The genetic groups are now moderately isolated from each other and have unique genetic characteristics. An estimation of the evolutionary history using the approximate Bayesian computation (ABC) approach suggested that Tohoku diverged from the common ancestor of S-Kanto and S-Izu. Then, S-Kanto and S-Izu split, and the recent hybridization between Tohoku and S-Kanto gave rise to N-Kanto. The origin of N-Kanto through the hybridization is relatively young and seems to be related to changes in the distributions of Tohoku and S-Kanto as a result of climatic oscillation in the Pleistocene. We conclude that the mitochondrial genome of S-Kanto was captured in Tohoku and that the original mitochondrial genome of Tohoku was entirely removed through hybridization.

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.

Phylogeographic structure of common sage (Salvia officinalis L.) reveals microrefugia throughout the Balkans and colonizations of the Apennines

Studying the population-genetic and phylogeographic structures of a representative species of a particular geographical region can not only provide us with information regarding its evolutionary history, but also improve our understanding of the evolutionary processes underlying the patterns of species diversity in that area. By analysing eight highly polymorphic microsatellite loci and two chloroplast DNA regions, we have investigated the influence of Pleistocene climate fluctuations on the evolutionary history of Salvia officinalis L. (common sage). The populations with the highest genetic diversity were located in the central parts of the Balkan distribution range. A large group of closely related haplotypes was distributed throughout the Balkans and the central Apennines, while the private lineage occupied the southern Apennines. In addition, two highly differentiated lineages were scattered only over the Balkans. The results suggest that a single refugium of the studied species from the last glacial period was located in the central part of the range in the Balkans. Numerous microrefugia, probably spanning several glaciation cycles, were scattered across the Balkans, while colonisation of the Apennines from the Balkans occurred at least on two occasions.

Population history, genetic variation, and conservation status of European white elm (Ulmus laevis Pall.) in Poland

Key message

The core populations of the European white elm (Ulmus laevis Pall.) located in Poland maintained slightly higher level of genetic diversity compared to the peripheral populations of this species.

Context

The most severe threat to elms is the loss of natural habitat under the pressures of agriculture and forestry as well as urbanization. The reductions in European white elm populations as well as populations of other elm species have also been caused by Dutch elm disease (DED). Previous studies have indicated a low level of genetic variation in Ulmus leavis Pall. However, in Poland, the genetic resources and demographic history of U. laevis populations remain poorly documented.

Aims

The genetic resources of U. laevis in Poland were identified and characterized. Additionally, tests were performed to identify potential bottleneck signatures and effective population sizes of the examined populations.

Methods

Polymorphism was analyzed using a set of six nuclear microsatellite markers (nSSRs) for 1672 individuals from 41 populations throughout the species range in Poland.

Results

(1) A moderate level of genetic variation was found. (2) A low genetic differentiation and lack of population structuring were identified. (3) Evidence of reduction in population size was found as a consequence of severe, past bottlenecks.

Conclusion

The loss of genetic diversity of U. laevis probably occurred in their refugia or shortly after the postglacial recolonization. This loss may have been affected by past DED pandemics similar to those seen at present.

Genetic pattern and demographic history of cutlassfish (Trichiurus nanhaiensis) in South China Sea by the influence of Pleistocene climatic oscillations

Trichiurus nanhaiensis is one of the most important commercial fish species in the South China Sea. This study aimed to investigate the level of genetic variation and population genetic structure of T. nanhaiensis in the South China Sea for the first time, using 281 individuals collected from seven locations along the coast of mainland China, Taiwan, and Hainan Island. A high level of haplotype diversity and low nucleotide diversity were detected in the mitochondrial DNA cyt b gene and nuDNA RYR 3 gene. The overall expected heterozygosity (He = 0.693) among the seven populations ranged from 0.681 to 0.706 in microsatellite DNA data, which revealed high levels of genetic diversity. Significant genetic differentiation was found in Taidong populations in Taiwan, revealing the prevention of gene flow caused by the Kuroshio Current. Two major lineages based on the cyt b gene suggested that the Taiwan Strait acted as a geographic barrier for T. nanhaiensis during the glacier periods in the late Pleistocene. The Bayesian skyline plot also revealed that population demographic expansion of T. nanhaiensis was estimated to have occurred in 0.1 Mya. Our results indicated that all populations of T. nanhaiensis had experienced a recent genetic bottleneck following recent expansion based on ABC analysis.

Global phylogeography of ridley sea turtles (Lepidochelys spp.): evolution, demography, connectivity, and conservation

Globally distributed marine taxa are well suited for investigations of biogeographic impacts on genetic diversity, connectivity, and population demography. The sea turtle genus Lepidochelys includes the wide-ranging and abundant olive ridley (L. olivacea), and the geographically restricted and ‘Critically Endangered’ Kemp’s ridley (L. kempii). To investigate their historical biogeography, we analyzed a large dataset of mitochondrial DNA (mtDNA) sequences from olive (n = 943) and Kemp’s (n = 287) ridleys, and genotyped 15 nuclear microsatellite loci in a global sample of olive ridleys (n = 285). We found that the ridley species split ~ 7.5 million years ago, before the Panama Isthmus closure. The most ancient mitochondrial olive ridley lineage, located in the Indian Ocean, was dated to ~ 2.2 Mya. Both mitochondrial and nuclear markers revealed significant structure for olive ridleys between Atlantic (ATL), East Pacific (EP), and Indo-West Pacific (IWP) areas. However, the divergence of mtDNA clades was very recent (< 1 Mya) with low within- clade diversity, supporting a recurrent extinction-recolonization model for these ocean regions. All data showed that ATL and IWP groups were more closely related than those in the EP, with mtDNA data supporting recent recolonization of the ATL from the IWP. Individual olive ridley dispersal between the ATL, EP, and IN/IWP could be interpreted as more male- than female-biased, and genetic diversity was lowest in the Atlantic Ocean. All populations showed signs of recent expansion, and estimated time frames were concordant with their recent colonization history. Investigating species abundance and distribution changes over time is central to evolutionary biology, and this study provides a historical biogeographic context for marine vertebrate conservation and management.

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

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

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

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

Geography vs. past climate: the drivers of population genetic structure of the Himalayan langur

Background

Contemporary species distribution, genetic diversity and evolutionary history in many taxa are shaped by both historical and current climate as well as topography. The Himalayas show a huge variation in topography and climatic conditions across its entire range, and have experienced major climatic fluctuations in the past. However, very little is known regarding how this heterogenous landscape has moulded the distribution of Himalayan fauna. A recent study examined the effect of these historical events on the genetic diversity of the Himalayan langurs in Nepal Himalaya. However, this study did not include the samples from the Indian Himalayan region (IHR). Therefore, here we revisit the questions addressed in the previous study with a near complete sampling from the IHR, along with the samples from the Nepal Himalaya. We used the mitochondrial Cytochrome-b (Cyt-b, 746 bp) region combined with multiple phylogeographic analyses and palaeodistribution modelling.

Results

Our dataset contained 144 sequences from the IHR as well as the Nepal Himalaya. Phylogenetic analysis showed a low divergent western clade nested within high divergent group of eastern lineages and in the network analysis we identified 22 haplotypes over the entire distribution range of the Himalayan langurs. Samples from the Nepal Himalaya showed geographically structured haplotypes corresponding to different river barriers, whereas samples from IHR showed star-like topology with no structure. Our statistical phylogeography analysis using diyABC supported the model of east to west colonisation of these langurs with founder event during colonisation. Analysis of demographic history showed that the effective population size of the Himalayan langurs decreased at the onset of last glacial maximum (LGM) and started increasing post LGM. The palaeodistribution modelling showed that the extent of suitable habitat shifted from low elevation central Nepal, and adjoining parts of north India, during LGM to the western Himalaya at present.

Conclusion

The current genetic diversity and distribution of Himalayan langurs in the Nepal Himalaya has been shaped by river barriers, whereas the rivers in the IHR had relatively less time to act as a strong genetic barrier after the recent colonisation event. Further, the post LGM expansion could have had confounding effect on Himalayan langur population structure in both Nepal Himalaya and IHR.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02054-1.

Demographic history and population genetic structure of Anisakis pegreffii in the cutlassfish Trichiurus japonicus along the coast of mainland China and Taiwan

Studying the genetic diversity of nematode parasite populations is crucial to gaining insight into parasite infection dynamics and informing parasite phylogeography. Anisakiasis is a zoonotic disease caused by the consumption of infectious third-stage larvae (L3) of Anisakis spp. carried by marine fish. In the present study, a total of 206 mitochondrial DNA sequences (cytochrome c oxidase 2, cox2) were used to study the genetic diversity, genetic structure, and historical demography of twelve A. pegreffii populations from Trichiurus japonicas along the coast of mainland China and Taiwan. Two distinct evolutionary lineages of A. pegreffii and no significant genealogical structures corresponding to sampling localities suggested that isolation in the marginal seas shaped their patterns of phylogeographic distribution along the coast of mainland China and Taiwan during glaciation with lower sea levels. Furthermore, pairwise F_ST values and AMOVA did not indicate any significant genetic differentiation among groups with no relation to the geographic area, which might be attributed to fewer barriers to gene flow as well as large population sizes. The results of the neutrality test, mismatch distribution, and Bayesian skyline plot analyses showed that entire population underwent population expansion during the late Pleistocene. Analysis of the demographic history revealed that A. pegreffii underwent historical lineage diversification and admixture due to secondary contact based on ABC analysis. The present research represents the first definitive population structure and demographic history across sampling locations of A. pegreffii along the coast of mainland China and Taiwan.

Population genetic diversity and structure of Rhinogobius candidianus (Gobiidae) in Taiwan: Translocation and release

Rhinogobius candidianus is a freshwater goby distributed in north, northwest, west, and south Taiwan, but this species has been introduced to east Taiwan and became dominant. To investigate its native population genetic diversity and structure and evaluate the sources and diversity of translocated populations, the mitochondrial DNA control region and cytochrome b gene (1981 bp) from 220 specimens were analyzed. These results indicated that (1) the east populations originated from two sources in west Taiwan; (2) translocated populations exist in east Taiwan and south Taiwan; (3) many populations have likely been moved secondarily by human intervention; (4) the effective size of the populations had declined greatly; (5) within the native populations, the ancestral populations colonized Taiwan during the land bridge phase in the Pleistocene through north Taiwan; (6) the landform changes in Taiwan shaped the population structure; and (7) the landforms of the coastline during glaciation also shaped the native range. The low-level genetic diversity, high population differentiation, and population decline greatly suggest the need for resource management and conservation interventions. Four clades (α-δ) should be managed as four distinct evolutionarily significant units, while the translocated populations should be managed as separate management units. Moreover, the translocated populations in east Taiwan should be evaluated and monitored carefully.

Palaeoclimatic changes resulted in range expansion and subsequent divergence in brown fur seals, Arctocephalus pusillus

Past climatic change as a driving force of marine diversification is still largely unclear, particularly for Southern Hemisphere species. Here, we present a case using the brown fur seal, Arctocephalus pusillus, assessing the geographical structure and demographic history using mitochondrial and nuclear data. Results show the two previously defined subspecies (one from Australia and the other from southern Africa) are phylogeographically distinct. Migration analyses based on nuclear data suggest the absence of migrants among the two genetically close assemblages. The demographic history of A. pusillus is characterized by a glacial population expansion (approx. 18 kya) in the southern African lineage, which coincides with time estimates of population expansion of prey species of seals. Approximate Bayesian calculations support an eastward dispersal event during the Last Glacial Maximum when sea levels were lower, followed by a postglacial divergence event, approximately 13 kya. The demographic history of the brown fur seal in the Southern Oceans provides support that recent palaeoclimatic changes could have facilitated expansions in some marine species and that postglacial sea-level rise may have acted as a dispersal barrier for species mostly confined to continental shelves.

Genomics reveals the history of a complex plant invasion and improves the management of a biological invasion from the South African-Australian biotic exchange

Many plants exchanged in the global redistribution of species in the last 200 years, particularly between South Africa and Australia, have become threatening invasive species in their introduced range. Refining our understanding of the genetic diversity and population structure of native and alien populations, introduction pathways, propagule pressure, naturalization, and initial spread, can transform the effectiveness of management and prevention of further introductions. We used 20,221 single nucleotide polymorphisms to reconstruct the invasion of a coastal shrub, Chrysanthemoides monilifera ssp. rotundata (bitou bush) from South Africa, into eastern Australia (EAU), and Western Australia (WAU). We determined genetic diversity and population structure across the native and introduced ranges and compared hypothesized invasion scenarios using Bayesian modeling. We detected considerable genetic structure in the native range, as well as differentiation between populations in the native and introduced range. Phylogenetic analysis showed the introduced samples to be most closely related to the southern-most native populations, although Bayesian analysis inferred introduction from a ghost population. We detected strong genetic bottlenecks during the founding of both the EAU and WAU populations. It is likely that the WAU population was introduced from EAU, possibly involving an unsampled ghost population. The number of private alleles and polymorphic SNPs successively decreased from South Africa to EAU to WAU, although heterozygosity remained high. That bitou bush remains an invasion threat in EAU, despite reduced genetic diversity, provides a cautionary biosecurity message regarding the risk of introduction of potentially invasive species via shipping routes.

The genetic consequences of captive breeding, environmental change and human exploitation in the endangered peninsular pronghorn

Endangered species with small population sizes are susceptible to genetic erosion, which can be detrimental to long-term persistence. Consequently, monitoring and mitigating the loss of genetic diversity are essential for conservation. The Peninsular pronghorn (Antilocapra americana peninsularis) is an endangered pronghorn subspecies that is almost entirely held in captivity. Captive breeding has increased the number of pronghorns from 25 founders in 1997 to around 700 individuals today, but it is unclear how the genetic diversity of the captive herd may have changed over time. We therefore generated and analysed data for 16 microsatellites spanning 2009-2021. We detected a decline in heterozygosity and an increase in the proportion of inbred individuals over time. However, these trends appear to have been partially mitigated by a genetically informed breeding management attempt that was implemented in 2018. We also reconstructed the recent demographic history of the Peninsular pronghorn, revealing two sequential population declines putatively linked to the desertification of the Baja California peninsula around 6000 years ago, and hunting and habitat loss around 500 years ago, respectively. Our results provide insights into the genetic diversity of an endangered antelope and indicate the potential for genetically informed management to have positive conservation outcomes.

Phylogeography of Tridentiger bifasciatus (Gobiidae) in the Northwestern Pacific

The shimofuri goby (Tridentiger bifasciatus) is native to marine, brackish, and fresh waters along the coasts of the northwest Pacific. Our study examined the population genetic structure, diversity, and demography of T. bifasciatus in the China Seas, including the Yellow Sea, East China Sea, and South China Sea, using the sequences of mitochondrial DNA [mtDNA cytochrome b (cyt b) gene and d-loop region] and nuclear DNA [nuDNA ryanodine receptor 3 (Ryr3) gene]. The mtDNA dataset revealed a significant population differentiation, but the nuDNA dataset displayed the absence of genetic differentiation. The discordance between these two datasets was accounted for by population admixture, selection, and incomplete lineage sorting. Although the mtDNA and nuDNA displayed a discordant population structure, these genetic markers revealed the same population history: (1) the populations retreated into two refugia during glaciations and (2) the populations declined recently. Our study revealed that after glaciations, the re-flooding in Taiwan Strait did not shape the migrations of the southern lineage from the South China Sea to the East China Sea, and displayed that two mtDNA lineages have diverged before they migrated southward during glaciations. These results offer important resources for the further study of conservation genetics.

EggLib 3: a Python package for population genetics and genomics

Rapid and repeatable polymorphism analyses have become a necessity with the current amount of genomic data that can be collected in many organisms. Traditionally such analyses are conducted using a variety of tools in combination, often requiring numerous format translation and manipulation. Here we present a massively updated version of our previous software package EggLib, intended to alleviate such costly and error-prone tinkering with the data. EggLib has been streamlined into a Python package and thoroughly updated and optimized to accommodate modern-day sized dataset. We show the main characteristics of the package making it a tool of choice to perform population genetics analyses. Once the data are imported (whatever their encoding), they can be filtered, edited, analyzed and compared to coalescent simulations very easily and efficiently. Furthermore the list of diversity and polymorphism statistics that can now be calculated has been greatly expanded. The software and its full documentation are available at https://egglib.org/.

Glacial vicariance and oceanic circulation shape population structure of the coastal legume Derris trifoliata in the Indo-West Pacific

PREMISE

The phylogeography of coastal plant species is shaped by contemporary and historical biogeographic processes. In this study, we aim to decipher the phylogeography of Derris trifoliata, a woody legume of relatively recent origin and wide distribution, in coastal areas in the Indo-West Pacific (IWP) region.

METHODS

Genetic diversity and population structure were assessed by analyzing six nuclear and three chloroplast DNA sequences from 30 populations across the species' range. Phylogeography was inferred by estimating gene flow, divergence time, historical population size changes, and historical habitat suitability using paleoclimatic niche modeling.

RESULTS

High genetic diversity was observed at the species level. The populations of three oceanic regions included in this study (i.e., Indian Ocean, South China Sea, and Pacific Ocean) formed distinct clades and likely diverged during the late Pleistocene. Potential barriers to gene flow were identified, including the Sunda and Sahul shelves, geographic distance, and current patterns of oceanic circulation. Analysis of changes in population size supported the bottleneck model, which was strengthened by estimates of habitat suitability across paleoclimatic conditions.

CONCLUSIONS

The once widespread distribution of D. trifoliata was fragmented by changes in climatic suitability and biogeographic barriers that arose following sea-level changes during the Pleistocene. In addition, contemporary patterns of oceanic circulation and geographic distance between populations appear to maintain genetic differentiation across its distribution in the IWP.