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

Molecular Phylogeography and Intraspecific Divergences in Siberian Wildrye (Elymus sibiricus L.) Wild Populations in China, Inferred From Chloroplast DNA Sequence and cpSSR Markers

A detailed understanding of the distribution and degree of genetic variation within a species is important for determining their evolutionary potential, which in return facilitates the development of efficient conservation strategies aimed at preserving adaptive genetic variation. As an important perennial, cool-season grass in temperate Eurasia, increasing attention has been paid to Siberian wildrye (Elymus sibiricus) due to its excellent ecological utilization value and forage production potential in China, particularly in the Qinghai-Tibet Plateau (QTP) regions. In this study, we applied two chloroplast (cp) genes (matK and rbcL), three cp spacer regions (trnY-GUA∼trnD-GUC, atpH∼atpF, and rps4∼trnT-UGU), and six cpSSR markers to the genetic and phylogenetic analysis of 137 wild E. sibiricus accessions from 23 natural populations that represent the main distribution regions in China. The results show the highest genetic diversity (h = 0.913) and haplotype richness (10 haplotypes) for the QTP population, which indicates QTP as the probable diversity center and geographic origin of E. sibiricus in China. Population divergence was high, indicating a significant phylogeographic structure together with a significantly higher N_st value (N_st > G_st, P < 0.05) at the species level, QTP+XJ (combined populations from QTP and Xinjiang), QTP+NC (combined populations from QTP and North China), and XJ+NC (combined populations from Xinjiang and North China) group levels, respectively. An expansion was revealed in the distributional range of E. sibiricus in China from paleo times up to the recent past, while a dramatic range of contraction was predicted for the near future. The predicted main limiting factor for the further spread of E. sibiricus is an increasing global mean temperature. We recommend that the combination of Es-cpDNA1 and Es-cpDNA3+4+5 can be used as effective markers for phylogenetic analysis and phylogeographical history analysis of E. sibiricus. These findings shed new light on the historical population dynamics of cold-season herbs in the QTP region and the north of China and are of great significance for the future establishment of protection and collection strategies for wild E. sibiricus germplasm.

Testing Phylogeographic Hypotheses in Mepraia (Hemiptera: Reduviidae) Suggests a Complex Spatio-Temporal Colonization in the Coastal Atacama Desert

Simple Summary

Mepraia is a blood-sucking bug endemic to Chile and a vector of the parasite that causes Chagas disease. Different colonization routes have been suggested for this bug; therefore, we tested different colonization routes using DNA sequences and bioinformatics approaches to select the most probable route. Our results suggest that, after the split of Triatoma, Mepraia divided into two main groups ~2.1 Mya. The northern group would have speciated between 1.7–1.4 Mya, giving rise to M. parapatrica, M. gajardoi and to a new, still undescribed lineage (Mepraia sp.). The southern group formed M. spinolai ~1.68 Mya. We suggest that Mepraia originated from the north-central Andes due to the last Andes uplift and hyperaridity. The hyperarid cycle would have separated the southern and northern groups. Then, within the northern group, colonization would have occurred from the centre to the north and south through corridors influenced by Pleistocene climatic changes. The habitat colonized by the southern clade led to only one species (M. spinolai). Fluctuations in climatic changes probably influenced speciation strongly in this kissing bug in the Atacama Desert.

Abstract

Mepraia is a genus (Triatominae) endemic to Chile and a vector of Trypanosoma cruzi. Alternative phylogeographic hypotheses have been suggested for Mepraia. We tested different colonization routes hypothesized using mitochondrial sequences and phylogeographic approaches to select the best-supported hypothesis. Our results suggest that, after the split from the sister genus Triatoma at ~4.3 Mya, Mepraia formed two main clades at ~2.1 Mya. The northern clade diverged from Mepraia sp. ~1.7 Mya, giving rise to M. parapatrica and M. gajardoi about ~1.4 Mya. The southern clade originated M. spinolai ~1.68 Mya. We suggest that Mepraia had an origin in the north-central Andes along with orogenic processes, reinforced by hyperaridity during the Pliocene. The hyperarid cycle would have separated the southern and northern clades. Then, in the northern clade, dispersal occurred north and south from the centre through corridors during the Pleistocene Climatic Oscillations. Climate changes may have induced a major speciation process in the Atacama Desert, while the more homogeneous habitat colonized by the southern clade led to only one, but structured, species.

RAD-Seq and Ecological Niche Reveal Genetic Diversity, Phylogeny, and Geographic Distribution of Kadsura interior and Its Closely Related Species

Most plants of Kadsura have economic value and medicinal application. Among them, K. interior and its closely related species have been demonstrated to have definite efficacy. However, the taxonomy and phylogenetic relationship of Kadsura in terms of morphology and commonly used gene regions remain controversial, which adversely affects its rational application. In this study, a total of 107 individuals of K. interior, K. heteroclita, K. longipedunculata, K. oblongifolia, and K. coccinea were studied from the perspectives of genetic diversity, phylogeny, and ecology via single nucleotide polymorphisms (SNPs) developed through restriction site-associated DNA sequencing (RAD-seq). Based on these SNPs, the genetic diversity, phylogenetic reconstruction, and population genetic structure were analyzed. Subsequently, divergence time estimation and differentiation scenario simulation were performed. Meanwhile, according to the species distribution records and bioclimatic variables, the Last Glacial Maximum and current potential distributions of five species were constructed, and the main ecological factors affecting the distribution of different species were extracted. The F ST calculated showed that there was a moderate degree of differentiation among K. heteroclita, K. longipedunculata, and K. oblongifolia, and there was a high degree of genetic differentiation between K. interior and the above species. The phylogenetic tree indicated that each of the species was monophyletic. The results of population genetic structure and divergence scenario simulation and D-statistics showed that there were admixture and gene flow among K. heteroclita, K. longipedunculata, and K. oblongifolia. The results of ecological niche modeling indicated that the distribution areas and the bioclimatic variables affecting the distribution of K. interior and its related species were different. This study explored the differences in the genetic divergence and geographical distribution patterns of K. interior and its related species, clarifying the uniqueness of K. interior compared to its relatives and providing a reference for their rational application in the future.

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

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

A genetic analysis of grey squirrel (Sciurus carolinensis) populations in Ireland

The grey squirrel (Sciurus carolinensis) is an invasive rodent in Ireland that has had detrimental impacts on the native Irish red squirrel (S. vulgaris) as well as on silviculture. This invasive species spread rapidly throughout Ireland, but in recent years appears to be declining in certain areas of the country. This study analysed the genetic profile of grey squirrel populations in Ireland to gain insight into their introduction, evolutionary history in Ireland, and vulnerability to management strategies. The genetic diversity and population structure of eight grey squirrel populations in Ireland was assessed using 11 species-specific microsatellite loci, and was compared to a small population from Tennessee, U.S.A., part of the native range of the grey squirrel. This is the first time these microsatellite markers developed specifically for grey squirrels have been used to study the species in Ireland. We found low to moderate genetic diversity overall across Irish populations as well as the presence of inbreeding. One population in particular, (in Co. Kildare), was differentiated from all other populations, which could indicate genetic isolation between Irish populations or a secondary introduction of S. carolinensis to Ireland.

Microsatellite Loci Reveal High Genetic Diversity, Mutation, and Migration Rates as Invasion Drivers of Callery Pear (Pyrus calleryana) in the Southeastern United States

Pyrus calleryana Decne. (Callery pear) is a deciduous tree native to China, Japan, Korea, and Taiwan. It is a popular ornamental tree in the United States (US) with early spring blooms and vibrant fall color. There are at least 26 cultivars of P. calleryana available in the US of which "Bradford" is the most well-known. Open-pollinated P. calleryana escapees are becoming one of the most common invasive tree species in the eastern United States. Developing better management practices for invasive P. calleryana requires detailed knowledge about reproductive biology and genetic diversity of the species, however, little is currently known about genetic variability within those open-pollinated populations. We investigated genetic diversity and population structure of non-cultivated, escaped P. calleryana populations within a ∼177 km radius in the southeastern United States. Because P. calleryana exhibits a range of morphological variation with great evolutionary potential, we hypothesized that a high genetic diversity would be manifested among escaped P. calleryana. Using 15 previously developed microsatellite loci, we genotyped 180 open-pollinated P. calleryana individuals that were collected across six naturally occurring sites in Tennessee, Georgia, and South Carolina, United States. Our results demonstrated the presence of a population structure with high genetic diversity, high gene flow, and high genetic differentiation between individuals across collection sites. Our results revealed that P. calleryana populations had differentiated shortly after the introduction to the US, most likely from specimens imported from Asia, consistent with historical records and our prior findings. The high invasive potential of the species is perhaps best underscored by transformation of P. calleryana specimens introduced from Asia into escape populations at continental scale across the United States. Our data also provided novel insight into potential issues that could be problematic for the future as P. calleryana may pose a potential threat to the economy, ecology, and native biodiversity in invaded areas.

Phylogeography of a canopy-forming kelp, Eisenia bicyclis (Laminariales, Phaeophyceae), based on a genome-wide sequencing analysis

Analyses of phylogeographic patterns and genetic diversity provide fundamental information for the management and conservation of species. However, little is published about these patterns in Japanese kelp species. In this study, we conducted phylogeographic analyses of a canopy-forming kelp, Eisenia bicyclis, based on genome-wide SNPs identified by ddRAD-seq. We obtained 1,299 SNPs for 76 samples from nine localities across the distribution. STRUCTURE, NeighborNet, and discriminant analysis of principal components consistently showed high genetic differentiation among the Eastern Pacific, Central Pacific, and Sea of Japan coastal regions. Relatively strong gene flow was detected only within populations in the Eastern Pacific and in the Sea of Japan. Genetic diversity and genetic uniqueness were high in the Central Pacific and low in the Sea of Japan. These results suggest that there were at least three independent refugia corresponding to the three regions during the Last Glacial Maximum (LGM). Furthermore, relatively larger populations in the Central Pacific and smaller populations in the Sea of Japan have been maintained in the demographic history from before the LGM to the present. These phylogeographic histories were supported by an Approximate Bayesian Computation analysis. From a conservation genetics perspective, the loss of southern populations in the Central Pacific would greatly reduce the total genetic diversity of the species. Southern populations in the Sea of Japan, which have relatively low genetic diversity, may be highly vulnerable to environmental change, such as heat waves and increased feeding. Therefore, careful monitoring and conservation are needed in the two regions.

Future-proofing the koala: synergising genomic and environmental data for effective species management

Climatic and evolutionary processes are inextricably linked to conservation. Avoiding extinction in rapidly changing environments often depends upon a species' capacity to adapt in the face of extreme selective pressures. Here, we employed exon capture and high-throughput next-generation sequencing to investigate the mechanisms underlying population structure and adaptive genetic variation in the koala (Phascolarctos cinereus), an iconic Australian marsupial that represents a unique conservation challenge because it is not uniformly threatened across its range. An examination of 250 specimens representing 91 wild source locations revealed that five major genetic clusters currently exist on a continental scale. The initial divergence of these clusters appears to have been concordant with the Mid-Brunhes Transition (∼ 430-300 kya), a major climatic reorganization that increased the amplitude of Pleistocene glacial-interglacial cycles. While signatures of polygenic selection and environmental adaptation were detected, strong evidence for repeated, climate-associated range contractions and demographic bottleneck events suggests that geographically isolated refugia may have played a more significant role in the survival of the koala through the Pleistocene glaciation than in situ adaptation. Consequently, the conservation of genome-wide genetic variation must be aligned with the protection of core koala habitat to increase the resilience of threatened populations to accelerating anthropogenic threats. Finally, we propose that the five major genetic clusters identified in this study should be accounted for in future koala conservation efforts (e.g. guiding translocations), as existing management divisions in the states of Queensland and New South Wales do not reflect historic or contemporary population structure.

Evolution and conservation genetics of an insular hemiparasitic plant lineage at the limit of survival: the case of Thesium sect. Kunkeliella in the Canary Islands

PREMISE

The diversification of island flora has been widely studied. However, the role of environmental niches in insular radiation processes has been less investigated. We combined population genetic analyses with species distribution modelling to clarify the genetic relationships, diversification patterns, species niche requirements, and conservation of Thesium sect. Kunkeliella, a clade of rare hemiparasitic plants endemic to the Canaries.

METHODS

We studied the three extant Thesium species and a new taxon from La Palma Island. We developed 12 microsatellites and performed population genetic analysis and studied the demographic history of the group. To evaluate the role of niche conservatism in the diversification of the group, we performed species distribution modelling (ESM) with four algorithms.

RESULTS

All species presented moderate genetic diversity values for rare endemics. Thesium canariense (Gran Canaria) showed high differentiation, whereas T. subsucculentum, T. retamoides (Tenerife), and La Palma populations are closely related. The lineage may have undergone a recent diversification with colonization proceeding east to west, with T. canariense as sister to the others. We detected a climatic niche shift, as taxa showed different distributions across the temperature gradient. There is enough evidence to describe La Palma populations as a new species.

CONCLUSIONS

We characterized the evolutionary history of Thesium sect. Kunkeliella by integrating genetic and ecological assessments. Our results indicate that this clade has undergone a recent radiation process with niche differentiation among species. The results increase our knowledge about insular radiations and will inform the conservation management of the study species.

Molecular phylogeography and species distribution modelling evidence of 'oceanic' adaptation for Actinidia eriantha with a refugium along the oceanic-continental gradient in a biodiversity hotspot

BACKGROUND

Refugia is considered to be critical for maintaining biodiversity; while discerning the type and pattern of refugia is pivotal for our understanding of evolutionary processes in the context of conservation. Interglacial and glacial refugia have been studied throughout subtropical China. However, studies on refugia along the oceanic-continental gradient have largely been ignored. We used a liana Actinidia eriantha, which occurs across the eastern moist evergreen broad-leaved forests of subtropical China, as a case study to test hypotheses of refugia along the oceanic-continental gradient and 'oceanic' adaptation.

RESULTS

The phylogeographic pattern of A. eriantha was explored using a combination of three cpDNA markers and 38 nuclear microsatellite loci, Species distribution modelling and dispersal corridors analysis. Our data showed intermediate levels of genetic diversity [haplotype diversity (hT) = 0.498; unbiased expected heterozygosity (UHE) = 0.510] both at the species and population level. Microsatellite loci revealed five clusters largely corresponding to geographic regions. Coalescent time of cpDNA lineages was dated to the middle Pliocene (ca. 4.03 Ma). Both geographic distance and climate difference have important roles for intraspecific divergence of the species. The Zhejiang-Fujian Hilly Region was demonstrated to be a refugium along the oceanic-continental gradient of the species and fit the 'refugia in refugia' pattern. Species distribution modelling analysis indicated that Precipitation of Coldest Quarter (importance of 44%), Temperature Seasonality (29%) and Mean Temperature of Wettest Quarter (25%) contributed the most to model development. By checking the isolines in the three climate layers, we found that A. eriantha prefer higher precipitation during the coldest quarter, lower seasonal temperature difference and lower mean temperature during the wettest quarter, which correspond to 'oceanic' adaptation. Actinidia eriantha expanded to its western distribution range along the dispersal corridor repeatedly during the glacial periods.

CONCLUSIONS

Overall, our results provide integrated evidence demonstrating that the Zhejiang-Fujian Hilly Region is a refugium along the oceanic-continental gradient of Actinidia eriantha in subtropical China and that speciation is attributed to 'oceanic' adaptation. This study gives a deeper understanding of the refugia in subtropical China and will contribute to the conservation and utilization of kiwifruit wild resources in the context of climate change.

Global invasion history of the emerging plant pathogen Phytophthora multivora

BACKGROUND

global trade in living plants and plant material has significantly increased the geographic distribution of many plant pathogens. As a consequence, several pathogens have been first found and described in their introduced range where they may cause severe damage on naïve host species. Knowing the center of origin and the pathways of spread of a pathogen is of importance for several reasons, including identifying natural enemies and reducing further spread. Several Phytophthora species are well-known invasive pathogens of natural ecosystems, including Phytophthora multivora. Following the description of P. multivora from dying native vegetation in Australia in 2009, the species was subsequently found to be common in South Africa where it does not cause any remarkable disease. There are now reports of P. multivora from many other countries worldwide, but not as a commonly encountered species in natural environments.

RESULTS

a global collection of 335 isolates from North America, Europe, Africa, Australia, the Canary Islands, and New Zealand was used to unravel the worldwide invasion history of P. multivora, using 10 microsatellite markers for all isolates and sequence data from five loci from 94 representative isolates. Our population genetic analysis revealed an extremely low heterozygosity, significant non-random association of loci and substantial genotypic diversity suggesting the spread of P. multivora readily by both asexual and sexual propagules. The P. multivora populations in South Africa, Australia, and New Zealand show the most complex genetic structure, are well established and evolutionary older than those in Europe, North America and the Canary Islands.

CONCLUSIONS

according to the conducted analyses, the world invasion of P. multivora most likely commenced from South Africa, which can be considered the center of origin of the species. The pathogen was then introduced to Australia, which acted as bridgehead population for Europe and North America. Our study highlights a complex global invasion pattern of P. multivora, including both direct introductions from the native population and secondary spread/introductions from bridgehead populations.

Population Genetic Analysis of Aedes aegypti Mosquitoes From Sudan Revealed Recent Independent Colonization Events by the Two Subspecies

Increases in arbovirus outbreaks in Sudan are vectored by Aedes aegypti, raising the medical importance of this mosquito. We genotyped 12 microsatellite loci in four populations of Ae. aegypti from Sudan, two from the East and two from the West, and analyzed them together with a previously published database of 31 worldwide populations to infer population structure and investigate the demographic history of this species in Sudan. Our results revealed the presence of two genetically distinct subspecies of Ae. aegypti in Sudan. These are Ae. aegypti aegypti in Eastern Sudan and Ae. aegypti formosus in Western Sudan. Clustering analysis showed that mosquitoes from East Sudan are genetically homogeneous, while we found population substructure in West Sudan. In the global context our results indicate that Eastern Sudan populations are genetically closer to Asian and American populations, while Western Sudan populations are related to East and West African populations. Approximate Bayesian Computation Analysis supports a scenario in which Ae. aegypti entered Sudan in at least two independent occasions nearly 70–80 years ago. This study provides a baseline database that can be used to determine the likely origin of new introductions for this invasive species into Sudan. The presence of the two subspecies in the country should be consider when designing interventions, since they display different behaviors regarding epidemiologically relevant parameters, such as blood feeding preferences and ability to transmit disease.

Global flyway evolution in red knots Calidris canutus and genetic evidence for a Nearctic refugium

Present‐day ecology and population structure are the legacies of past climate and habitat perturbations, and this is particularly true for species that are widely distributed at high latitudes. The red knot, Calidris canutus, is an arctic‐breeding, long‐distance migratory shorebird with six recognized subspecies defined by differences in morphology, migration behavior, and annual cycle phenology, in a global distribution thought to have arisen just since the last glacial maximum (LGM). We used nextRAD sequencing of 10,881 single‐nucleotide polymorphisms (SNPs) to assess the neutral genetic structure and phylogeographic history of 172 red knots representing all known global breeding populations. Using population genetics approaches, including model‐based scenario‐testing in an approximate Bayesian computation (ABC) framework, we infer that red knots derive from two main lineages that diverged ca. 34,000 years ago, and thus most probably persisted at the LGM in both Palearctic and Nearctic refugia, followed by at least two instances of secondary contact and admixture. Within two Beringian subspecies (C. c. roselaari and rogersi), we detected previously unknown genetic structure among sub‐populations sharing a migratory flyway, reflecting additional complexity in the phylogeographic history of the region. Conversely, we found very weak genetic differentiation between two Nearctic populations (rufa and islandica) with clearly divergent migratory phenotypes and little or no apparent contact throughout the annual cycle. Together, these results suggest that relative gene flow among migratory populations reflects a complex interplay of historical, geographical, and ecological factors.

Microsatellite Markers from Peronospora tabacina, the Cause of Blue Mold of Tobacco, Reveal Species Origin, Population Structure, and High Gene Flow

Peronospora tabacina is an obligate parasite that causes blue mold of tobacco. The pathogen reproduces primarily by sporangia, whereas the sexual oospores are rarely observed. A collection of 122 isolates of P. tabacina was genotyped using nine microsatellites to assess the population structure of individuals from subpopulations collected from central, southern, and western Europe; the Middle East; Central America; North America; and Australia. Genetic variations among the six subpopulations accounted for ∼8% of the total variation, including moderate levels of genetic differentiation, high gene flow among these subpopulations, and a positive correlation between geographic and genetic distance (r = 0.225; P < 0.001). Evidence of linkage disequilibrium (P < 0.001) showed that populations contained partially clonal subpopulations but that subpopulations from Australia and Mediterranean Europe did not. High genetic variation and population structure among samples could be explained by continuous gene flow across continents via infected transplant exchange and/or long-distance dispersal of sporangia via wind currents. This study analyzed the most numerous P. tabacina collection and allowed conclusions regarding the migration, mutation, and evolutionary history of this obligate biotrophic oomycete. The evidence pointed to the species origin in Australia and identified intracontinental and intercontinental migration patterns of this important pathogen.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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

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

More opportunities more species: Pleistocene differentiation and northward expansion of an evergreen broad-leaved tree species Machilus thunbergii (Lauraceae) in Southeast China

BACKGROUND

The broad continuum between tropical and temperate floras in Eastern Asia (EAS) are thought to be one of the main factors responsible for a prominent species diversity anomaly of temperate plants between EAS and eastern North America (ENS). However, how the broad continuum and niche evolution between tropical and temperate floras in EAS contributes to lineage divergence and species diversity remains largely unknown.

RESULTS

Population genetic structure, demography, and determinants of genetic structure [i.e., isolation-by-distance (IBD), isolation-by-resistance (IBR), and isolation-by-environment (IBE)] of Machilus thunbergii Sieb. et Zucc. (Lauraceae) were evaluated by examining sequence variation of ten low-copy nuclear genes across 43 populations in southeast China. Climatic niche difference and potential distributions across four periods (Current, mid-Holocene, the last glacial maximum, the last interglacial) of two genetic clusters were determined by niche modelling. North and south clusters of populations in M. thunbergii were revealed and their demarcation line corresponds well with the northern boundary of tropical zone in China of Zhu & Wan. The divergence time between the clusters and demographic expansion of M. thunbergii occurred after the mid-Pleistocene climate transition (MPT, 0.8-1.2 Ma). Migration rates between clusters were asymmetrical, being much greater from north to south than the reverse. Significant effects of IBE, but non-significant effects of IBD and IBR on population genetic divergence were detected. The two clusters have different ecological niches and require different temperature regimes.

CONCLUSIONS

The north-south genetic differentiation may be common across the temperate-tropical boundary in southeast China. Divergent selection under different temperature regimes (possibly above and below freezing temperature in winter) could account for this divergence pattern. The broad continuum between tropical and temperate floras in EAS may have provided ample opportunities for tropical plant lineages to acquire freezing tolerance and to colonize the temperate regions during the late-Cenozoic global cooling. Our findings shed deeper insights into the high temperate plant species diversity in EAS.

Assessing model adequacy leads to more robust phylogeographic inference

Phylogeographic studies base inferences on large data sets and complex demographic models, but these models are applied in ways that could mislead researchers and compromise their inference. Researchers face three challenges associated with the use of models: (i) 'model selection', or the identification of an appropriate model for analysis; (ii) 'evaluation of analytical results', or the interpretation of the biological significance of the resulting parameter estimates, delimitations, and topologies; and (iii) 'model evaluation', or the use of statistical approaches to assess the fit of the model to the data. The field collectively invests most of its energy in point (ii) without considering the other points; we argue that attention to points (i) and (iii) is essential to phylogeographic inference.

Invasion genomics uncover contrasting scenarios of genetic diversity in a widespread marine invader

Invasion rates have increased in the past 100 y irrespective of international conventions. What characterizes a successful invasion event? And how does genetic diversity translate into invasion success? Employing a whole-genome perspective using one of the most successful marine invasive species world-wide as a model, we resolve temporal invasion dynamics during independent invasion events in Eurasia. We reveal complex regionally independent invasion histories including cases of recurrent translocations, time-limited translocations, and stepping-stone range expansions with severe bottlenecks within the same species. Irrespective of these different invasion dynamics, which lead to contrasting patterns of genetic diversity, all nonindigenous populations are similarly successful. This illustrates that genetic diversity, per se, is not necessarily the driving force behind invasion success. Other factors such as propagule pressure and repeated introductions are an important contribution to facilitate successful invasions. This calls into question the dominant paradigm of the genetic paradox of invasions, i.e., the successful establishment of nonindigenous populations with low levels of genetic diversity.

Wheat speciation and adaptation: perspectives from reticulate evolution

Reticulate evolution through the interchanging of genetic components across organisms can impact significantly on the fitness and adaptation of species. Bread wheat (Triticum aestivum subsp. aestivum) is one of the most important crops in the world. Allopolyploid speciation, frequent hybridization, extensive introgression, and occasional horizontal gene transfer (HGT) have been shaping a typical paradigm of reticulate evolution in bread wheat and its wild relatives, which is likely to have a substantial influence on phenotypic traits and environmental adaptability of bread wheat. In this review, we outlined the evolutionary history of bread wheat and its wild relatives with a highlight on the interspecific hybridization events, demonstrating the reticulate relationship between species/subspecies in the genera Triticum and Aegilops. Furthermore, we discussed the genetic mechanisms and evolutionary significance underlying the introgression of bread wheat and its wild relatives. An in-depth understanding of the evolutionary process of Triticum species should be beneficial to future genetic study and breeding of bread wheat.

Pleistocene expansion, anthropogenic pressure and ocean currents: Disentangling the past and ongoing evolutionary history of Patella aspera Röding, 1798 in the archipelago of Madeira

AIMS

Rising sea-level following the Last Glacial Maximum lead to fragmentation of coastal limpet populations between islands of the Archipelago of Madeira. This fragmentation is reinforced by recent heavy exploitation reducing effective population size on Madeira Island. We use the limpet P. aspera to understand how the role of processes at different time scales (i.e. changes in the sea level and overexploitation) can influence the genetic composition of an extant species, relating these processes to reproductive phenology and seasonal shifts in ocean currents.

LOCATION

Madeira Island, Porto Santo and Desertas (Archipelago of Madeira, NE Atlantic Ocean).

TAXON

The limpet Patella aspera.

METHODS

Twelve microsatellite genetic markers were used. A power analysis was used to evaluate the power of the microsatellite markers to detect a signal of population differentiation. Long-term past migrations were assessed using a Bayesian Markov Montecarlo approach in the software MIGRATE-n to estimate mutation-scaled migration rates (M = m/μ; m, probability of a lineage immigrating per generation; μ, mutation rate). Two scenarios were evaluated using an Approximate Bayesian Computation (ABC) in the software DIYABC 2.1 (i) Scenario 1: considered a population scenario from a reduced Ne at time t₃ to a higher Ne at time t₂; and (ii) Scenario 2 considering a reduction of Ne from a time t₃ to a time t₂.

RESULTS

Colonization of the archipelago by Portuguese settlers six centuries ago probably led to an important decrease in the genetic diversity of the species (N_e). Contemporary gene flow strongly support a pattern of high asymmetric connectivity explained by the reproductive phenology of the species and spatio-temporal seasonal changes in the ocean currents. Spatio-temporal reconstructions using Bayesian methods, including coalescent and Approximate Bayesian Computation (ABC) approaches, suggest changes in the migration patterns from highly symmetric to highly asymmetric connectivity with subtle population differentiation as consequence of post-glacial maximum sea level rise during the Holocene.

MAIN CONCLUSIONS

Our results suggest that anthropogenic activity could have had serious effects on the genetic diversity of heavily exploited littoral species since the end of the Pleistocene, probably accelerating in recent years.

Phylogeographic pattern of a cryptoviviparous mangrove, Aegiceras corniculatum, in the Indo-West Pacific, provides insights for conservation actions

This study identified the historical geoclimatic factors which caused low genetic diversity and strong phylogeographic structure in a cryptoviviparous mangrove. The phylogeographic pattern was used to suggest conservation actions. Phylogeographic studies are used to understand the spatial distribution and evolution of genetic diversity, and have major conservation implications, especially for threatened taxa like the mangroves. This study aimed to assess the phylogeographic pattern of Aegiceras corniculatum, a cryptoviviparous mangrove, across its distribution range in the Indo-West Pacific (IWP) region. We genotyped 398 samples, collected from 37 populations, at four chloroplast DNA (cpDNA) loci, and identified the influence of historical processes on the contemporary population structure of the species. Low genetic diversity at the population level was observed. The evolutionary relationship between 12 cpDNA haplotypes suggested a strong phylogeographic structure, which was further validated by the clustering algorithms and proportioning of maximum variation among hierarchical population groups. The magnitude and direction of historical gene flow indicated that the species attained its wide distribution from its likely ancestral area of the Malay Archipelago. The divergence time estimates of the haplotypes indicated that the geoclimatic changes during the Pleistocene, especially the glacial sea-level changes and emergence of landmasses, hindered genetic exchange and created genetic differentiation between the phylogenetic groups. The species overwintered the last glacial maxima in multiple refugia in the IWP, as identified by the environmental niche modelling. Overall, our findings indicated that ancient glacial vicariance had influenced the present genetic composition of A. corniculatum, which was maintained by the current demographic features of this region. We discussed how these findings can be used to prioritize areas for conservation actions, restore disturbed habitats and prevent further genetic erosion.

The evolutionary heritage and ecological uniqueness of Scots pine in the Caucasus ecoregion is at risk of climate changes

Scots pine is one of the most widely occurring pines, but future projections suggest a large reduction in its range, mostly at the southern European limits. A significant part of its range is located in the Caucasus, a global hot-spot of diversity. Pine forests are an important reservoir of biodiversity and endemism in this region. We explored demographic and biogeographical processes that shaped the genetic diversity of Scots pine in the Caucasus ecoregion and its probable future distribution under different climate scenarios. We found that the high genetic variability of the Caucasian populations mirrors a complex glacial and postglacial history that had a unique evolutionary trajectory compared to the main range in Europe. Scots pine currently grows under a broad spectrum of climatic conditions in the Caucasus, which implies high adaptive potential in the past. However, the current genetic resources of Scots pine are under high pressure from climate change. From our predictions, over 90% of the current distribution of Scots pine may be lost in this century. By threatening the stability of the forest ecosystems, this would dramatically affect the biodiversity of the Caucasus hot-spot.

Genetic diversity in North American Cercis Canadensis reveals an ancient population bottleneck that originated after the last glacial maximum

Understanding of the present-day genetic diversity, population structure, and evolutionary history of tree species can inform resource management and conservation activities, including response to pressures presented by a changing climate. Cercis canadensis (Eastern Redbud) is an economically valuable understory tree species native to the United States (U.S.) that is also important for forest ecosystem and wildlife health. Here, we document and explain the population genetics and evolutionary history of this deciduous tree species across its distributed range. In this study, we used twelve microsatellite markers to investigate 691 wild-type trees sampled at 74 collection sites from 23 Eastern U.S. states. High genetic diversity and limited gene flow were revealed in wild, natural stands of C. canadensis with populations that are explained by two major genetic clusters. These findings indicate that an ancient population bottleneck occurred coinciding with the last glacial maximum (LGM) in North America. The structure in current populations likely originated from an ancient population in the eastern U.S. that survived LGM and then later diverged into two contemporary clusters. Data suggests that populations have expanded since the last glaciation event from one into several post-glacial refugia that now occupy this species’ current geographic range. Our enhanced understanding benchmarks the genetic variation preserved within this species and can direct future efforts in conservation, and resource utilization of adaptively resilient populations that present the greatest genetic and structural diversity.

Phylogeography of the East Asian grassland plant, Viola orientalis (Violaceae), inferred from plastid and nuclear restriction site-associated DNA sequencing data

To elucidate the origin and migration history of the "Mansen elements," a group of temperate grassland plants mainly distributed in northeastern Asia, phylogeographic analyses based on chloroplast DNA markers and double-digest restriction site-associated DNA sequencing (ddRAD-seq) data were performed on Viola orientalis, one of the representative species of the group. Phylogenetic analyses using ddRAD-seq data revealed that the populations of V. orientalis were clustered into five clades, among which the continental clades made of populations from Russia and Korea diverged more than 100,000 years earlier than the Japanese clades. The Japanese clade likely diverged during the last glacial period, followed by a further post-glacial divergence into the Kyushu and the Honshu subclades. Our study demonstrated that V. orientalis originated in the continental area of northeastern Asia and, during the last glacial period, has spread southward through the Korean Peninsula across the Japanese Islands. This finding supports the previously proposed evolutionary hypothesis regarding the origin and migration routes of the Mansen elements.

Sea surface temperature, rather than land mass or geographic distance, may drive genetic differentiation in a species complex of highly dispersive seabirds

Seabirds, particularly Procellariiformes, are highly mobile organisms with a great capacity for long dispersal, though simultaneously showing high philopatry, two conflicting life-history traits that may lead to contrasted patterns of genetic population structure. Landmasses were suggested to explain differentiation patterns observed in seabirds, but philopatry, isolation by distance, segregation between breeding and nonbreeding zones, and oceanographic conditions (sea surface temperatures) may also contribute to differentiation patterns. To our knowledge, no study has simultaneously contrasted the multiple factors contributing to the diversification of seabird species, especially in the gray zone of speciation. We conducted a multilocus phylogeographic study on a widespread seabird species complex, the little shearwater complex, showing highly homogeneous morphology, which led to considerable taxonomic debate. We sequenced three mitochondrial and six nuclear markers on all extant populations from the Atlantic (lherminieri) and Indian Oceans (bailloni), that is, five nominal lineages from 13 populations, along with one population from the eastern Pacific Ocean (representing the dichrous lineage). We found sharp differentiation among populations separated by the African continent with both mitochondrial and nuclear markers, while only mitochondrial markers allowed characterizing the five nominal lineages. No differentiation could be detected within these five lineages, questioning the strong level of philopatry showed by these shearwaters. Finally, we propose that Atlantic populations likely originated from the Indian Ocean. Within the Atlantic, a stepping-stone process accounts for the current distribution. Based on our divergence time estimates, we suggest that the observed pattern of differentiation mostly resulted from historical and current variation in sea surface temperatures.

Origin and dispersion pathways of guava in the Galapagos Islands inferred through genetics and historical records

Guava (Psidium guajava) is an aggressive invasive plant in the Galapagos Islands. Determining its provenance and genetic diversity could explain its adaptability and spread, and how this relates to past human activities. With this purpose, we analyzed 11 SSR markers in guava individuals from Isabela, Santa Cruz, San Cristobal, and Floreana islands in the Galapagos, as well as from mainland Ecuador. The mainland guava population appeared genetically differentiated from the Galapagos populations, with higher genetic diversity levels found in the former. We consistently found that the Central Highlands region of mainland Ecuador is one of the most likely origins of the Galapagos populations. Moreover, the guavas from Isabela and Floreana show a potential genetic input from southern mainland Ecuador, while the population from San Cristobal would be linked to the coastal mainland regions. Interestingly, the proposed origins for the Galapagos guava coincide with the first human settlings of the archipelago. Through approximate Bayesian computation, we propose a model where San Cristobal was the first island to be colonized by guava from the mainland, and then, it would have spread to Floreana and finally to Santa Cruz; Isabela would have been seeded from Floreana. An independent trajectory could also have contributed to the invasion of Floreana and Isabela. The pathway shown in our model agrees with the human colonization history of the different islands in the Galapagos. Our model, in conjunction with the clustering patterns of the individuals (based on genetic distances), suggests that guava introduction history in the Galapagos archipelago was driven by either a single event or a series of introduction events in rapid succession. We thus show that genetic analyses supported by historical sources can be used to track the arrival and spread of invasive species in novel habitats and the potential role of human activities in such processes.

Large- and small-scale geographic structures affecting genetic patterns across populations of an Alpine butterfly

Understanding factors influencing patterns of genetic diversity and the population genetic structure of species is of particular importance in the current era of global climate change and habitat loss. These factors include the evolutionary history of a species as well as heterogeneity in the environment it occupies, which in turn can change across time. Most studies investigating spatio-temporal genetic patterns have focused on patterns across wide geographic areas rather than local variation, but the latter can nevertheless be important particularly in topographically complex areas. Here, we consider these issues in the Sooty Copper butterfly (Lycaena tityrus) from the European Alps, using genome-wide SNPs identified through RADseq. We found strong genetic differentiation within the Alps with four genetic clusters, indicating western, central, and eastern refuges, and a strong reduction of genetic diversity from west to east. This reduction in diversity may suggest that the southwestern refuge was the largest one in comparison to other refuges. Also, the high genetic diversity in the west may result from (a) admixture of different western refuges, (b) more recent demographic changes, or (c) introgression of lowland L. tityrus populations. At small spatial scales, populations were structured by several landscape features and especially by high mountain ridges and large river valleys. We detected 36 outlier loci likely under altitudinal selection, including several loci related to membranes and cellular processes. We suggest that efforts to preserve alpine L. tityrus should focus on the genetically diverse populations in the western Alps, and that the dolomite populations should be treated as genetically distinct management units, since they appear to be currently more threatened than others. This study demonstrates the usefulness of SNP-based approaches for understanding patterns of genetic diversity, gene flow, and selection in a region that is expected to be particularly vulnerable to climate change.

Phylogeography of the East Asian grassland plant, Viola orientalis (Violaceae), inferred from plastid and nuclear restriction site-associated DNA sequencing data

To elucidate the origin and migration history of the "Mansen elements," a group of temperate grassland plants mainly distributed in northeastern Asia, phylogeographic analyses based on chloroplast DNA markers and double-digest restriction site-associated DNA sequencing (ddRAD-seq) data were performed on Viola orientalis, one of the representative species of the group. Phylogenetic analyses using ddRAD-seq data revealed that the populations of V. orientalis were clustered into five clades, among which the continental clades made of populations from Russia and Korea diverged more than 100,000 years earlier than the Japanese clades. The Japanese clade likely diverged during the last glacial period, followed by a further post-glacial divergence into the Kyushu and the Honshu subclades. Our study demonstrated that V. orientalis originated in the continental area of northeastern Asia and, during the last glacial period, has spread southward through the Korean Peninsula across the Japanese Islands. This finding supports the previously proposed evolutionary hypothesis regarding the origin and migration routes of the Mansen elements.

Demographic inferences and climatic niche modelling shed light on the evolutionary history of the emblematic cold-adapted Apollo butterfly at regional scale

Cold-adapted species escape climate warming by latitudinal and/or altitudinal range shifts, and currently occur in Southern Europe in isolated mountain ranges within "sky islands". Here, we studied the genetic structure of the Apollo butterfly in five such sky islands (above 1,000 m) in France, and infer its demographic history since the last interglacial, using single nucleotide polymorphisms (ddRADseq SNPs). The Auvergne and Alps populations show strong genetic differentiation but not alpine massifs, although separated by deep valleys. Combining three complementary demographic inference methods and species distribution models (SDMs) we show that the LIG period was highly unfavourable for Apollo that probably survived in small population in the highest summits of Auvergne. The population shifted downslope and expanded eastward between LIG and LGM throughout the large climatically suitable Rhône valley between the glaciated summits of Auvergne and Alps. The Auvergne and Alps populations started diverging before the LGM but remained largely connected till the mid-Holocene. Population decline in Auvergne was more gradual but started before (~7 kya vs. 800 ya), and was much stronger with current population size ten times lower than in the Alps. In the Alps, the low genetic structure and limited evidence for isolation by distance suggest a nonequilibrium metapopulation functioning. The core Apollo population experienced cycles of contraction-expansion with climate fluctuations with largely interconnected populations overtime according to a "metapopulation-pulsar" functioning. This study demonstrates the power of combining demographic inferences and SDMs to determine past and future evolutionary trajectories of an endangered species at a regional scale.

Incipient speciation and its impact on taxonomic decision: a case study using a sky island sister-species pair of stag beetles (Lucanidae: Lucanus)

Species delimitation can be difficult when the divergence between focal taxa is in the incipient stage of speciation, because conflicting results are expected among different data sets, and the species limits can differ depending on the species concept applied. We studied speciation history and investigated the impact on taxonomic decision-making when using different types of data in a Taiwanese endemic sister-species pair of stag beetles, Lucanus miwai and Lucanus yulaoensis, from sky island habitats. We showed that the two geographical taxa can be diagnosed by male mandibular shape. We found two mitochondrial co1 lineages with pairwise sequence divergence &gt; 3%; however, L. miwai might not be monophyletic. The result of our multispecies coalescent-based species delimitation using five nuclear loci supported the evolutionary independence of the two sister species, but the calculated values of the genealogical divergence index (gdi) corresponded to the ambiguous zone of species delimitation. We also showed that post-divergence gene flow is unlikely. Our study demonstrates challenges in the delineation of incipient species, but shows the importance of understanding the speciation history and adopting integrative approaches to reconcile seemingly conflicting results before making evolutionarily relevant taxonomic decisions.

Spatial dynamics of Chinese Muntjac related to past and future climate fluctuations

Climate fluctuations in the past and in the future are likely to result in population expansions, shifts, or the contraction of the ecological niche of many species, and potentially leading to the changes in their geographical distributions. Prediction of suitable habitats has been developed as a useful tool for the assessment of habitat suitability and resource conservation to protect wildlife. Here, we model the ancestral demographic history of the extant modern Chinese Muntjac Muntiacus reevesi populations using approximate Bayesian computation (ABC) and used the maximum entropy model to simulate the past and predict the future spatial dynamics of the species under climate oscillations. Our results indicated that the suitable habitats for the M. reevesi shifted to the Southeast and contracted during the Last Glacial Maximum, whereas they covered a broader and more northern position in the Middle Holocene. The ABC analyses revealed that the modern M. reevesi populations diverged in the Middle Holocene coinciding with the significant contraction of the highly suitable habitat areas. Furthermore, our predictions suggest that the potentially suitable environment distribution for the species will expand under all future climate scenarios. These results indicated that the M. reevesi diverged in the recent time after the glacial period and simultaneously as its habitat’s expanded in the Middle Holocene. Furthermore, the past and future climate fluctuation triggered the change of Chinese muntjac spatial distribution, which has great influence on the Chinese muntjac’s population demographic history.

Evidence for independent domestication of sheep mtDNA lineage A in India and introduction of lineage B through Arabian sea route

India ranks the second in the world in terms of its sheep population with approximately 74.26 million represented by 44 well-described breeds in addition to several non-descript populations. Genetic diversity and phylogeography of Indian sheep breeds remain poorly understood, particularly for south Indian breeds. To have a comprehensive view of the domestication history of Indian sheep, we sequenced the mitochondrial DNA (mtDNA) control region (D-loop) and cytochrome b gene (CYTB) of 16 Indian domestic sheep breeds, most of them (13) from the south India. We analysed these sequences along with published data of domestic and wild sheep from different countries, including India. The haplotype diversity was relatively high in Indian sheep, which were classified into the three known mtDNA lineages, namely A, B and C. Lineage A was predominant among Indian sheep whereas lineages B and C were observed at low frequencies but C was restricted to the breeds of north and east India. The median joining network showed five major expanding haplogroups of lineage A (A1–A5). Out of which, A2, A4 and A5 were more frequent in Indian sheep in contrast to breeds from other parts of the world. Among the 27 Indian sheep breeds analysed, Mandya and Sonadi breeds were significantly different from other Indian breeds in the MDS analyses. This was explained by a very high contribution of lineage B into these two breeds. The Approximate Bayesian Computation (ABC) provided evidence for the domestication of lineage A sheep in the Indian subcontinent. Contrary to the current knowledge, we also found strong support for the introduction of lineage B into Indian subcontinent through sea route rather than from the Mongolian Plateau. The neighbour-joining tree of domestic and wild sheep revealed the close genetic relationship of Indian domestic sheep with Pakistani wild sheep O. vignei blanfordi. Based on our analyses and archaeological evidences, we suggest the Indian subcontinent as one of the domestication centres of the lineage A sheep, while lineage B sheep might have arrived into India from elsewhere via Arabian sea route. To the best of our knowledge, this is the first comprehensive study on Indian sheep where we have analysed more than 740 animals belonging to 27 sheep breeds raised in various regions of India. Our study provides insight into the understanding of the origin and migratory history of Indian sheep.

A recovery network leads to the natural recolonization of an archipelago and a potential trailing edge refuge

Rapid environmental change is reshaping ecosystems and driving species loss globally. Carnivore populations have declined and retracted rapidly and have been the target of numerous translocation projects. Success, however, is complicated when these efforts occur in novel ecosystems. Identifying refuges, locations that are resistant to environmental change, within a translocation framework should improve population recovery and persistence. American martens (Martes americana) are the most frequently translocated carnivore in North America. As elsewhere, martens were extirpated across much of the Great Lakes region by the 1930s and, despite multiple translocations beginning in the 1950s, martens remain of regional conservation concern. Surprisingly, martens were rediscovered in 2014 on the Apostle Islands of Lake Superior after a putative absence of >40 yr. To identify the source of martens to the islands and understand connectivity of the reintroduction network, we collected genetic data on martens from the archipelago and from all regional reintroduction sites. In total, we genotyped 483 individual martens, 43 of which inhabited the Apostle Islands (densities 0.42-1.46 km^-2 ). Coalescent analyses supported the contemporary recolonization of the Apostle Islands with progenitors likely originating from Michigan, which were sourced from Ontario. We also identified movements by a first-order relative between the Apostle Islands and the recovery network. We detected some regional gene flow, but in an unexpected direction: individuals moving from the islands to the mainland. Our findings suggest that the Apostle Islands were naturally recolonized by progeny of translocated individuals and now act as a source back to the reintroduction sites on the mainland. We suggest that the Apostle Islands, given its protection from disturbance, complex forest structure, and reduced carnivore competition, will act as a potential refuge for marten along their trailing range boundary and a central node for regional recovery. Our work reveals that translocations, even those occurring along southern range boundaries, can create recovery networks that function like natural metapopulations. Identifying refuges, locations that are resistant to environmental change, within these recovery networks can further improve species recovery, even within novel environments. Future translocation planning should a priori identify potential refuges and sources to improve short-term recovery and long-term persistence.