ISOLATION BY DISTANCE

Modeling the velocity of evolving lineages and predicting dispersal patterns

Accurate estimation of the dispersal velocity or speed of evolving organisms is no mean feat. In fact, existing probabilistic models in phylogeography or spatial population genetics generally do not provide an adequate framework to define velocity in a relevant manner. For instance, the very concept of instantaneous speed simply does not exist under one of the most popular approaches that models the evolution of spatial coordinates as Brownian trajectories running along a phylogeny. Here, we introduce a family of models-the so-called Phylogenetic Integrated Velocity (PIV) models-that use Gaussian processes to explicitly model the velocity of evolving lineages instead of focusing on the fluctuation of spatial coordinates over time. We describe the properties of these models and show an increased accuracy of velocity estimates compared to previous approaches. Analyses of West Nile virus data in the United States indicate that PIV models provide sensible predictions of the dispersal of evolving pathogens at a one-year time horizon. These results demonstrate the feasibility and relevance of predictive phylogeography in monitoring epidemics in time and space.

Population genomics informs the management of harvested snappers across north-western Australia

Failure to consider population structure when managing harvested fishes increases the risk of stock depletion, yet empirical estimates of population structure are often lacking for important fishery species. In this study, we characterise genetic variation in single nucleotide polymorphisms (SNPs) to assess population structure for three harvested species of tropical snappers across the broad (up to 300 km wide) and extensive (~ 4000 km) continental shelf of north-western Australia. Comparisons across ~ 300 individuals per species, showed remarkably similar patterns of genetic structure among Lutjanus sebae (red emperor), L. malabaricus (saddletail snapper) and Pristipomoides multidens (goldband snapper) despite subtle differences in biological and ecological traits. Low levels of genetic subdivision were reflected in an isolation by distance relationship where genetic connectivity increased with geographic proximity. This indicates extensive but not unlimited dispersal across the north-western Australian shelf. Our findings provide evidence of connectivity between current management areas, violating the assumption of multiple independent stocks. Spatial stock assessment models may be more suitable for the management of these species however demographic connectivity rates cannot be accurately estimated from the conventional population genetic approaches applied in this study. We recommend that managers aim to maintain adequate spawning biomass across current management areas, and assess stocks at finer scales, where practical.

Evolutionary and Ecological Processes Determining the Properties of the Matrix

AbstractThe G matrix is the matrix of additive genetic variances and covariances for a vector of phenotypes. Here we apply the classical theory for the balance among selection, genetic drift, and mutations to find the contributions to G from each locus for populations at stasis. The fitness is approximated by a linear function of phenotypes, with coefficients affected by environmental fluctuations. We show that the G matrix can be decomposed into four additive components generated by selection, drift, mutations, and environmental fluctuations. Selection is on average counteracted by the other three processes included in Fisher's concept of deterioration of the environment, generating considerable changes in mean phenotypes. The theory illustrates that neither Fisher's fundamental theorem nor Lande's classical gradient formula is sufficient for assessing adaptive changes through time unless the deteriorations are corrected for. This applies for populations at stasis, but also for populations that are subject to long-term evolutionary changes. The theory also indicates several possible comparative studies for investigations of deteriorating effects. Our analyses also suggest that the factor loadings to the eigenvector of the G matrix with the lowest eigenvalue will rather accurately indicate the relative contributions from different phenotype components to fitness. This is information notoriously difficult to obtain in natural populations.

Influence of geographic isolation and the environment on gene flow among phenotypically diverse lizards

Lizards in the genus Anolis comprise hundreds of species that display a wide range of phenotypic variation closely related to their environment. One example is the Guadeloupean anole (Anolis marmoratus ssp.) that display extreme phenotypic variation, primarily in adult male color and pattern, with twelve described subspecies on the archipelago. Here we examine the relationship between phenotypic and genetic divergence among five subspecies on the two main islands and test the role of geographic isolation and the environment in reducing gene flow. We also examined two offshore island populations to assess the impact of complete geographic isolation on gene flow. We analyzed color phenotypes by measuring spectral reflectance and genomic diversity using SNPs. Genetic divergence was correlated with dorsolateral head and body color phenotypes, and slope and geographic distance were nearly equivalent at explaining this divergence. There was minimal genome-wide divergence at neutral loci among phenotypically disparate subspecies on the two main islands and their differentiation is consistent with a model of divergence with gene flow. Our spatial visualization of gene flow showed an impact of environmental features consistent with a hypothesis of ecologically driven divergence. Nonetheless, subspecies on the two main islands remain interconnected by substantial gene flow and their phenotypic variation is likely maintained at selection-gene flow equilibrium by divergent selection at loci associated with their color phenotypes. Greater isolation, such as inhabiting a remote island, may be required for reducing gene flow. Our findings highlight the role of the environment, adaptation, and geographic isolation on gene flow.

Population Genomic Scans for Natural Selection and Demography

Uncovering the fundamental processes that shape genomic variation in natural populations is a primary objective of population genetics. These processes include demographic effects such as past changes in effective population size or gene flow between structured populations. Furthermore, genomic variation is affected by selection on nonneutral genetic variants, for example, through the adaptation of beneficial alleles or balancing selection that maintains genetic variation. In this article, we discuss the characterization of these processes using population genetic models, and we review methods developed on the basis of these models to unravel the underlying processes from modern population genomic data sets. We briefly discuss the conditions in which these approaches can be used to infer demography or identify specific nonneutral genetic variants and cases in which caution is warranted. Moreover, we summarize the challenges of jointly inferring demography and selective processes that affect neutral variation genome-wide.

Modeling the velocity of evolving lineages and predicting dispersal patterns

Accurate estimation of the dispersal velocity or speed of evolving organisms is no mean feat. In fact, existing probabilistic models in phylogeography or spatial population genetics generally do not provide an adequate framework to define velocity in a relevant manner. For instance, the very concept of instantaneous speed simply does not exist under one of the most popular approaches that models the evolution of spatial coordinates as Brownian trajectories running along a phylogeny (Lemey et al., 2010). Here, we introduce a new family of models - the so-called "Phylogenetic Integrated Velocity" (PIV) models - that use Gaussian processes to explicitly model the velocity of evolving lineages instead of focusing on the fluctuation of spatial coordinates over time. We describe the properties of these models and show an increased accuracy of velocity estimates compared to previous approaches. Analyses of West Nile virus data in the U.S.A. indicate that PIV models provide sensible predictions of the dispersal of evolving pathogens at a one-year time horizon. These results demonstrate the feasibility and relevance of predictive phylogeography in monitoring epidemics in time and space.

Population connectivity and size reductions in the Anthropocene: the consequence of landscapes and historical bottlenecks in white forsythia fragmented habitats

BACKGROUND

White forsythia (Abeliophyllum distichum) is an endangered Korean Peninsula endemic that has been subjected to recent population genomics studies using SNPs via RAD sequencing. Here, we primarily employed the often underutilized haplotype information from RAD loci to further describe the species' previously uninvestigated haplotype-based genomic variation and structure, and genetic-geographic characteristics and gene flow patterns among its five earlier identified genetic groups. We also inferred the time of past events that may have impacted the effective population size of these groups, as well as the species' potential future distribution amidst the warming climate and anthropogenic threats.

RESULTS

Our findings emphasized the recognition of the species' regional patterns of genetic structure, and the role of topography and its associated gene flow patterns as some of the possible factors that may have influenced the species' present-day fragmented population distribution. The inferred bottleneck events during the Anthropocene, some of which aligned with the time of historical catastrophic events on the Peninsula (e.g., the Korean War), were revealed to have contributed to the generally low effective population size of its five lineages, particularly those with marginal distributional range. Future distribution under both optimistic and pessimistic climatic scenarios suggests unlikely suitable habitats for these populations to expand from their current range limits, at least in the next 80 years.

CONCLUSIONS

The small effective population size and landscape-driven limited gene flow among white forsythia populations will remain a big challenge for the conservation management of the species' already fragmented population distribution. To help mitigate these impacts, the merging of various research approaches and the use of genomic data to their full potential is recommended to provide the optimized knowledge-based tools for the conservation of this endangered species, and other similar plants under pressure.

Revealing the evolutionary history and contemporary population structure of Pacific salmon in the Fraser River through genome resequencing

The Fraser River once supported massive salmon returns. However, over the last century, the largest returns have consistently been less than half of the recorded historical maximum. There is substantial interest from surrounding communities and governments to increase salmon returns for both human use and functional ecosystems. To generate resources for this endeavor, we resequenced genomes of Chinook (Oncorhynchus tshawytscha), coho (Oncorhynchus kisutch), and sockeye salmon (Oncorhynchus nerka) from the Fraser River at moderate coverage (∼16×). A total of 954 resequenced genomes were analyzed, with 681 collected specifically for this study from tissues sampled between 1997 and 2021. An additional 273 were collected from previous studies. At the species level, Chinook salmon appeared to have 1.6-2.1× more SNPs than coho or sockeye salmon, respectively. This difference may be attributable to large historical declines of coho and sockeye salmon. At the population level, 3 Fraser River genetic groups were identified for each species using principal component and admixture analyses. These were consistent with previous research and supports the continued use of these groups in conservation and management efforts. Environmental factors and a migration barrier were identified as major factors influencing the boundaries of these genetic groups. Additionally, 20 potentially adaptive loci were identified among the genetic groups. This information may be valuable in new management and conservation efforts. Furthermore, the resequenced genomes are an important resource for contemporary genomics research on Fraser River salmon and have been made publicly available.

Notes from the taxonomic disaster zone: Evolutionary drivers of intractable species boundaries in an Australian lizard clade (Scincidae: Ctenotus)

Genomic-scale datasets, sophisticated analytical techniques, and conceptual advances have disproportionately failed to resolve species boundaries in some groups relative to others. To understand the processes that underlie taxonomic intractability, we dissect the speciation history of an Australian lizard clade that arguably represents a "worst-case" scenario for species delimitation within vertebrates: the Ctenotus inornatus species group, a clade beset with decoupled genetic and phenotypic breaks, uncertain geographic ranges, and parallelism in purportedly diagnostic morphological characters. We sampled hundreds of localities to generate a genomic perspective on population divergence, structure, and admixture. Our results revealed rampant paraphyly of nominate taxa in the group, with lineages that are either morphologically cryptic or polytypic. Isolation-by-distance patterns reflect spatially continuous differentiation among certain pairs of putative species, yet genetic and geographic distances are decoupled in other pairs. Comparisons of mitochondrial and nuclear gene trees, tests of nuclear introgression, and historical demographic modelling identified gene flow between divergent candidate species. Levels of admixture are decoupled from phylogenetic relatedness; gene flow is often higher between sympatric species than between parapatric populations of the same species. Such idiosyncratic patterns of introgression contribute to species boundaries that are fuzzy while also varying in fuzziness. Our results suggest that "taxonomic disaster zones" like the C. inornatus species group result from spatial variation in the porosity of species boundaries and the resulting patterns of genetic and phenotypic variation. This study raises questions about the origin and persistence of hybridizing species and highlights the unique insights provided by taxa that have long eluded straightforward taxonomic categorization.

Peruvian origin and global invasions of five continents by the highly damaging agricultural pest Liriomyza huidobrensis (Diptera: Agromyzidae)

Identification of the geographic origin of invasive species can be critical to effective management and amelioration of negative impacts in the introduced range. Liriomyza huidobrensis is a polyphagous leafmining fly that is a devastating pest of many vegetable and floriculture crops around the world. Considered native to South and possibly Central America, L. huidobrensis became invasive in the 1980s and has since spread to at least 30 countries on five continents. We used phylogeographic analysis of over 2 kb of mitochondrial cytochrome oxidase I and II sequence data from 403 field-collected specimens from both native and introduced populations to investigate the geographic origins of invasive L. huidobrensis worldwide. Within South America, there was substantial genetic variation, as well as the strong phylogeographic structure typical of a native range. In contrast, leafminers from the introduced range and Central America all contained little genetic variation and shared the same small set of haplotypes. These haplotypes trace to Peru as the ultimate geographic origin of invasive populations. Central America is rejected as part of the original geographic range of L. huidobrensis. Within Peru, the primary export region of Lima shared an extremely similar pattern of reduced haplotype variation to the invasive populations. An additional 18 specimens collected at US ports of entry did not share the same haplotype profile as contemporary invasive populations, raising perplexing questions on global pathways and establishment success in this species.

Understanding species limits through the formation of phylogeographic lineages

The outcomes of speciation across organismal dimensions (e.g., ecological, genetic, phenotypic) are often assessed using phylogeographic methods. At one extreme, reproductively isolated lineages represent easily delimitable species differing in many or all dimensions, and at the other, geographically distinct genetic segments introgress across broad environmental gradients with limited phenotypic disparity. In the ambiguous gray zone of speciation, where lineages are genetically delimitable but still interacting ecologically, it is expected that these lineages represent species in the context of ontology and the evolutionary species concept when they are maintained over time with geographically well-defined hybrid zones, particularly at the intersection of distinct environments. As a result, genetic structure is correlated with environmental differences and not space alone, and a subset of genes fail to introgress across these zones as underlying genomic differences accumulate. We present a set of tests that synthesize species delimitation with the speciation process. We can thereby assess historical demographics and diversification processes while understanding how lineages are maintained through space and time by exploring spatial and genome clines, genotype-environment interactions, and genome scans for selected loci. Employing these tests in eight lineage-pairs of snakes in North America, we show that six pairs represent 12 "good" species and that two pairs represent local adaptation and regional population structure. The distinct species pairs all have the signature of divergence before or near the mid-Pleistocene, often with low migration, stable hybrid zones of varying size, and a subset of loci showing selection on alleles at the hybrid zone corresponding to transitions between distinct ecoregions. Locally adapted populations are younger, exhibit higher migration, and less ecological differentiation. Our results demonstrate that interacting lineages can be delimited using phylogeographic and population genetic methods that properly integrate spatial, temporal, and environmental data.

Habitat fragmentation strongly restricts gene flow in endangered ectomycorrhizal fungal populations: Evidence from Rhizopogon togasawarius, specific to Pseudotsuga japonica, across the entire distribution range

Habitat fragmentation reduces gene flow, causing genetic differentiation and diversity loss in endangered species through genetic drift and inbreeding. However, the impact of habitat fragmentation on ectomycorrhizal (ECM) fungi remains unexplored, despite their critical roles in forest ecosystems. Here, we investigated the population genetic structure and the demographic history of Rhizopogon togasawarius, the ECM fungus specifically colonizing the host tree Pseudotsuga japonica, across its entire distribution range (>200 km). These two species are designated as endangered species on the IUCN Red List since they are found only in small, fragmented forests in Japan. We analysed 236 R. togasawarius individuals from five remaining populations across the Kii Peninsula and the Shikoku Island, separated by a sea channel. Simple sequence repeat analyses using 20 loci revealed strong genetic differentiation among populations (FST = 0.255), even significant in the nearest population pair separated by a distance of only 8 km (FST = 0.075), indicating extremely limited gene flow between populations. DIYABC-RF analyses implied that population divergence occurred approximately 6000 generations ago between the two regions, and nearly 1500 generations ago between the nearest populations within Shikoku Island, related to past climate events. Because of prolonged genetic isolation, significant inbreeding was confirmed in four of five populations, where effective population sizes became very small (Ne = 9.0-58.0). Although evaluation of extinction risks for microorganisms is challenging, our conservation genetic results indicated that habitat fragmentation increases extinction risk through population genetic mechanisms, and therefore should not be overlooked in biodiversity conservation efforts.

Demographic inference for spatially heterogeneous populations using long shared haplotypes

We introduce a modified spatial Λ-Fleming-Viot process to model the ancestry of individuals in a population occupying a continuous spatial habitat divided into two areas by a sharp discontinuity of the dispersal rate and effective population density. We derive an analytical formula for the expected number of shared haplotype segments between two individuals depending on their sampling locations. This formula involves the transition density of a skew diffusion which appears as a scaling limit of the ancestral lineages of individuals in this model. We then show that this formula can be used to infer the dispersal parameters and the effective population density of both regions, using a composite likelihood approach, and we demonstrate the efficiency of this method on a range of simulated data sets.

Of Clams and Clades: Genetic Diversity and Connectivity of Small Giant Clams (Tridacna maxima) in the Southern Pacific Ocean

Giant clams (Tridacna and Hippopus) are large marine bivalves occupying tropical and subtropical reefs in the Indo-Pacific. Giant clam populations have declined in many areas of the Indo-Pacific and continue to be threatened by harvesting and environmental change. The small giant clam (Tridacna maxima) occurs throughout the Indo-Pacific and has been subject to several phylogeographic studies across its range. However, given its broad range, there are several areas where the genetic diversity and connectivity of T. maxima populations has not been characterised. Here, we analyse the mitochondrial marker cytochrome oxidase 1 (CO1) to examine the genetic diversity and connectivity of T. maxima in two regions: Australia's Coral Sea Marine Park and the Cook Islands. Samples were collected from 13 reefs within the Coral Sea Marine Park and ten islands within the Cook Islands archipelago. Tridacna maxima across the sampled region of the Coral Sea did not display any population structure, whereas significant population structure was detected for T. maxima within the Cook Islands. For the Cook Islands, most pairwise comparisons involving an island in the northern group (Manihiki) were significant, as were comparisons for Palmerston (a more centrally located island) and the southern islands, Rarotonga and Mangaia. Both regions displayed high haplotype diversities (> 0.90), indicating that they are important repositories of genetic diversity. Additional CO1 data from throughout T. maxima's distribution showed that the Coral Sea clams belonged to the clade occurring in the South-Western Pacific Ocean, whilst those from the Cook Islands belonged to a unique clade found in the Central Pacific Ocean. This clade extended from Fiji in the west to French Polynesia in the east and the atolls of Palmyra and Tarawa (Kiribati) in the north. Our assessment of genetic diversity and population structure in these regions will assist with management decisions for the species.

Comparing Genetic Ne Reconstructions Over Time With Long-Time Wolf Monitoring Data in Two Populations

Many methods are now available to calculate N e , but their performance varies depending on assumptions. Although simulated data are useful to discover certain types of bias, real empirical data supported by detailed known population histories allow us to discern how well methods perform with actual messy and complex data. Here, we focus on two genomic data sets of grey wolf populations for which population size changes of the past 40-120 years are well documented. We use this background to explore in what detail we can retrieve the known population history from these populations, in the light of pitfalls relating to population history, sampling design and the change in the spatial scale at which N e is estimated as we go further back in time. The Scandinavian wolf population was founded in the early 1980s from a few individuals and has gradually expanded up to 510 wolves. Although the founder event of the Scandinavian population was detected by GONE, the founding effective population size was strongly overestimated when the most recent samples were used, but less so when older samples were considered. Nevertheless, the present-day N e corresponds to theoretical expectations. The western Great Lakes wolf population of Minnesota is the only population in the contiguous United States that persisted throughout the 20th century, surviving intense persecution. We found a good concordance between the estimated N e and trends in census size data, but the reconstruction of N e clearly highlights the difficulty of interpreting results in spatially structured populations that underwent demographic fluctuations.

Identification of non-model mammal species using the MinION DNA sequencer from Oxford Nanopore

Background

The Neotropics harbors the largest species richness of the planet; however, even in well-studied groups, there are potentially hundreds of species that lack a formal description, and likewise, many already described taxa are difficult to identify using morphology. Specifically in small mammals, complex morphological diagnoses have been facilitated by the use of molecular data, particularly from mitochondrial sequences, to obtain accurate species identifications. Obtaining mitochondrial markers implies the use of PCR and specific primers, which are largely absent for non-model organisms. Oxford Nanopore Technologies (ONT) is a new alternative for sequencing the entire mitochondrial genome without the need for specific primers. Only a limited number of studies have employed exclusively ONT long-reads to assemble mitochondrial genomes, and few studies have yet evaluated the usefulness of such reads in multiple non-model organisms.

Methods

We implemented fieldwork to collect small mammals, including rodents, bats, and marsupials, in five localities in the northern extreme of the Cordillera Central of Colombia. DNA samples were sequenced using the MinION device and Flongle flow cells. Shotgun-sequenced data was used to reconstruct the mitochondrial genome of all the samples. In parallel, using a customized computational pipeline, species-level identifications were obtained based on sequencing raw reads (Whole Genome Sequencing). ONT-based identifications were corroborated using traditional morphological characters and phylogenetic analyses.

Results

A total of 24 individuals from 18 species were collected, morphologically identified, and deposited in the biological collection of Universidad EAFIT. Our different computational pipelines were able to reconstruct mitochondrial genomes from exclusively ONT reads. We obtained three new mitochondrial genomes and eight new molecular mitochondrial sequences for six species. Our species identification pipeline was able to obtain accurate species identifications for up to 75% of the individuals in as little as 5 s. Finally, our phylogenetic analyses corroborated the identifications from our automated species identification pipeline and revealed important contributions to the knowledge of the diversity of Neotropical small mammals.

Discussion

This study was able to evaluate different pipelines to reconstruct mitochondrial genomes from non-model organisms, using exclusively ONT reads, benchmarking these protocols on a multi-species dataset. The proposed methodology can be applied by non-expert taxonomists and has the potential to be implemented in real-time, without the need to euthanize the organisms and under field conditions. Therefore, it stands as a relevant tool to help increase the available data for non-model organisms, and the rate at which researchers can characterize life specially in highly biodiverse places as the Neotropics.

Habitat connectivity drives panda recovery

Globally, the majority of habitat loss is irreversible, and most species will never recover their former ranges. We have learned a great deal about what leads to population decline and extinction, but less about recovery. The recently downlisted giant panda provides a unique opportunity to understand the mechanisms of species recovery. In our study, we estimate giant panda suitable habitats, population density, and gene flow across landscapes to fully investigate the direct and indirect ecological mechanisms underlying bold conservation strategies. We found that the Giant Panda National Survey has modestly but systematically underestimated population size. China's effort to mitigate anthropogenic disturbances was associated with increased panda population density through improving habitat quality and reducing habitat fragmentation. Enhanced landscape connectivity reduced inbreeding via gene flow but indirectly increased inbreeding temporarily due to high local panda density. Although the panda's recovery has been geographically uneven, we provide evidence for improving connectivity and gene flow resulting from conservation efforts. If these processes can be sustained and improved, the panda's path to recovery will be less encumbered by loss of genetic diversity, fostering hope that the present rate of recovery will not be stalled. Findings from this study will not only help guide future giant panda conservation management but also provide a model for how a more mechanistic examination of the genetic processes underlying species recovery can foster the development of more effective strategies for endangered species recovery.

Genetic investigation of population structure in Atlantic chub mackerel, Scomber colias Gmelin, 1789 along the West African coast

Sustainable management of transboundary fish stocks hinges on accurate delineation of population structure. Genetic analysis offers a powerful tool to identify potential subpopulations within a seemingly homogenous stock, facilitating the development of effective, coordinated management strategies across international borders. Along the West African coast, the Atlantic chub mackerel (Scomber colias) is a commercially important and ecologically significant species, yet little is known about its genetic population structure and connectivity. Currently, the stock is managed as a single unit in West African waters despite new research suggesting morphological and adaptive differences. Here, eight microsatellite loci were genotyped on 1,169 individuals distributed across 33 sampling sites from Morocco (27.39°N) to Namibia (22.21°S). Bayesian clustering analysis depicts one homogeneous population across the studied area with null overall differentiation (F ST = 0.0001ns), which suggests panmixia and aligns with the migratory potential of this species. This finding has significant implications for the effective conservation and management of S. colias within a wide scope of its distribution across West African waters from the South of Morocco to the North-Centre of Namibia and underscores the need for increased regional cooperation in fisheries management and conservation.

Where the "ruber" Meets the Road: Using the Genome of the Red Diamond Rattlesnake to Unravel the Evolutionary Processes Driving Venom Evolution

Understanding the proximate and ultimate causes of phenotypic variation is fundamental in evolutionary research, as such variation provides the substrate for selection to act upon. Although trait variation can arise due to selection, the importance of neutral processes is sometimes understudied. We presented the first reference-quality genome of the Red Diamond Rattlesnake (Crotalus ruber) and used range-wide 'omic data to estimate the degree to which neutral and adaptive evolutionary processes shaped venom evolution. We characterized population structure and found substantial genetic differentiation across two populations, each with distinct demographic histories. We identified significant differentiation in venom expression across age classes with substantially reduced but discernible differentiation across populations. We then used conditional redundancy analysis to test whether venom expression variation was best predicted by neutral divergence patterns or geographically variable (a)biotic factors. Snake size was the most significant predictor of venom variation, with environment, prey availability, and neutral sequence variation also identified as significant factors, though to a lesser degree. By directly including neutrality in the model, our results confidently highlight the predominant, yet not singular, role of life history in shaping venom evolution.

Dealing With Assumptions and Sampling Bias in the Estimation of Effective Population Size: A Case Study in an Amphibian Population

Accurately estimating effective population size (N e) is essential for understanding evolutionary processes and guiding conservation efforts. This study investigates N e estimation methods in spatially structured populations using a population of moor frog (Rana arvalis) as a case study. We assessed the behaviour of N e estimates derived from the linkage disequilibrium (LD) method as we changed the spatial configuration of samples. Moor frog eggs were sampled from 25 breeding patches (i.e., separate vernal ponds, ditches or parts of larger fens) within a single population, revealing an isolation-by-distance pattern and a local spatial genetic structure. Varying buffer sizes around each patch were used to examine the impact of sampling window size on the estimation of effective number of breeders (N b). Our results indicate a downward bias in LD N b estimates with increasing buffer size, suggesting an underestimation of N b. The observed bias is attributed to LD resulting from including genetically divergent individuals (mixture-LD) confounding LD due to drift. This emphasises the significance of considering even subtle spatial genetic patterns. The implications of these findings are discussed, emphasising the need to account for spatial genetic structure to accurately assess population viability and inform conservation efforts. This study contributes to our understanding of the challenges associated with N e estimation in spatially structured populations and underscores the importance of refining methodologies to address population-specific spatial dynamics for effective conservation planning and management.

Artificial intelligence enables unified analysis of historical and landscape influences on genetic diversity

While genetic variation in any species is potentially shaped by a range of processes, phylogeography and landscape genetics are largely concerned with inferring how environmental conditions and landscape features impact neutral intraspecific diversity. However, even as both disciplines have come to utilize SNP data over the last decades, analytical approaches have remained for the most part focused on either broad-scale inferences of historical processes (phylogeography) or on more localized inferences about environmental and/or landscape features (landscape genetics). Here we demonstrate that an artificial intelligence model-based analytical framework can consider both deeper historical factors and landscape-level processes in an integrated analysis. We implement this framework using data collected from two Brazilian anurans, the Brazilian sibilator frog (Leptodactylus troglodytes) and granular toad (Rhinella granulosa). Our results indicate that historical demographic processes shape most the genetic variation in the sibulator frog, while landscape processes primarily influence variation in the granular toad. The machine learning framework used here allows both historical and landscape processes to be considered equally, rather than requiring researchers to make an a priori decision about which factors are important.

Disentangling effects of dispersal, environment and anthropogenic barriers on functional connectivity in aquatic systems

Disentangling the roles of structural landscape factors and animal movement behaviour can present challenges for practitioners managing landscapes to maintain functional connectivity and achieve conservation goals. We used a landscape genetics approach to combine robust demographic, behavioural and genetic datasets with spatially explicit simulations to evaluate the effects of anthropogenic barriers (dams, culverts) and natural landscape resistance (gradient, elevation) affecting dispersal behaviour, genetic connectivity and genetic structure in a resident population of Westslope Cutthroat Trout (Oncorhynchus clarkii lewisi). Analyses based on 10 years of sampling effort revealed a pattern of restricted dispersal, and population genetics identified discrete population clusters between distal tributaries and the mainstem stream and no structure within the mainstem stream. Demogenetic simulations demonstrated that, for this population, the effects of existing anthropogenic barriers on population structure are redundant with effects of restricted dispersal associated with the underlying environmental resistance. Our approach provides an example of how extensive field sampling combined with landscape genetics can be incorporated into spatially explicit simulation modelling to explore how, together, movement ecology and landscape resistance can be used to inform decisions around restoration and connectivity.

Ecological filtering drives rapid spatiotemporal dynamics in fish skin microbiomes

Skin microbiomes provide vital functions, yet knowledge about the drivers and processes structuring their species assemblages is limited-especially for non-model organisms. In this study, fish skin microbiome was assessed by high throughput sequencing of amplicon sequence variants from metabarcoding of V3-V4 regions in the 16S rRNA gene on fish hosts subjected to the following experimental manipulations: (i) translocation between fresh and brackish water habitats to investigate the role of environment; (ii) treatment with an antibacterial disinfectant to reboot the microbiome and investigate community assembly and priority effects; and (iii) maintained alone or in pairs to study the role of social environment and inter-host dispersal of microbes. The results revealed that fish skin microbiomes harbour a highly dynamic microbial composition that was distinct from bacterioplankton communities in the ambient water. Microbiome composition first diverged as an effect of translocation to either the brackish or freshwater habitat. When the freshwater individuals were translocated back to brackish water, their microbiome composition converged towards the fish microbiomes in the brackish habitat. In summary, external environmental conditions and individual-specific factors jointly determined the community composition dynamics, whereas inter-host dispersal had negligible effects. The dynamics of the microbiome composition was seemingly non-affected by reboot treatment, pointing towards high resilience to disturbance. The results emphasised the role of inter-individual variability for the unexplained variation found in many host-microbiome systems, although the mechanistic underpinnings remain to be identified.

Ecological, genetic and geographical divergence explain differences in colouration among sunbird species (Nectariniidae)

How extravagant ornamental traits evolve is a key question in evolutionary biology. Bird plumages are among the most elaborate ornaments, displaying almost all colours of the rainbow. Why and how birds evolved to be so colourful remains an open question with multiple and sometimes competing hypotheses. Different colours in different patches (i.e. body parts) might have different functions and thus result from different forms of selection (e.g. natural vs. sexual selection). Here we test the influence of three factors on colour diversity in sunbirds: (1) geographical distance, (2) differences in light environment and (3) phylogenetic distances. We show that both natural and sexual selection affect the evolution of sunbird colouration, but that their extent and direction differs between sexes, and varies with the extent of species overlap and across different patches on the body. Even though overlap in light environment partially explains colour differences among species, no colour metric (brightness, hue or chroma) covaries with light environment. Our results suggest that multiple forms of selection influence the colouration of different colour patches in different ways across an organism's body, highlighting the need to investigate colouration as a network of individual but inter-connected colour patches. These results are likely to be generalizable across the multitude of colourful animals.

A spatial approach to jointly estimate Wright's neighborhood size and long-term effective population size

Spatially continuous patterns of genetic differentiation, which are common in nature, are often poorly described by existing population genetic theory or methods that assume either panmixia or discrete, clearly definable populations. There is therefore a need for statistical approaches in population genetics that can accommodate continuous geographic structure, and that ideally use georeferenced individuals as the unit of analysis, rather than populations or subpopulations. In addition, researchers are often interested in describing the diversity of a population distributed continuously in space; this diversity is intimately linked to both the dispersal potential and the population density of the organism. A statistical model that leverages information from patterns of isolation by distance to jointly infer parameters that control local demography (such as Wright's neighborhood size), and the long-term effective size (Ne) of a population would be useful. Here, we introduce such a model that uses individual-level pairwise genetic and geographic distances to infer Wright's neighborhood size and long-term Ne. We demonstrate the utility of our model by applying it to complex, forward-time demographic simulations as well as an empirical dataset of the two-form bumblebee (Bombus bifarius). The model performed well on simulated data relative to alternative approaches and produced reasonable empirical results given the natural history of bumblebees. The resulting inferences provide important insights into the population genetic dynamics of spatially structured populations.

Scots pine - panmixia and the elusive signal of genetic adaptation

Scots pine is the foundation species of diverse forested ecosystems across Eurasia and displays remarkable ecological breadth, occurring in environments ranging from temperate rainforests to arid tundra margins. Such expansive distributions can be favored by various demographic and adaptive processes and the interactions between them. To understand the impact of neutral and selective forces on genetic structure in Scots pine, we conducted range-wide population genetic analyses on 2321 trees from 202 populations using genotyping-by-sequencing, reconstructed the recent demography of the species and examined signals of genetic adaptation. We found a high and uniform genetic diversity across the entire range (global FST 0.048), no increased genetic load in expanding populations and minor impact of the last glacial maximum on historical population sizes. Genetic-environmental associations identified only a handful of single-nucleotide polymorphisms significantly linked to environmental gradients. The results suggest that extensive gene flow is predominantly responsible for the observed genetic patterns in Scots pine. The apparent missing signal of genetic adaptation is likely attributed to the intricate genetic architecture controlling adaptation to multi-dimensional environments. The panmixia metapopulation of Scots pine offers a good study system for further exploration into how genetic adaptation and plasticity evolve under gene flow and changing environment.

Isolation by environment and its consequences for range shifts with global change: Landscape genomics of the invasive plant common tansy

Invasive species are a growing global economic and ecological problem. However, it is not well understood how environmental factors mediate invasive range expansion. In this study, we investigated the recent and rapid range expansion of common tansy across environmental gradients in Minnesota, USA. We densely sampled individuals across the expanding range and performed reduced representation sequencing to generate a dataset of 3071 polymorphic loci for 176 individuals. We used non-spatial and spatially explicit analyses to determine the relative influences of geographic distance and environmental variation on patterns of genomic variation. We found no evidence for isolation by distance but strong evidence for isolation by environment, indicating that environmental factors may have modulated patterns of range expansion. Land use classification and soils were particularly important variables related to population structure although they operated on different spatial scales; land use classification was related to broad-scale patterns and soils were related to fine-scale patterns. All analyses indicated a distinctive genetic cluster in the most recently invaded portion of the range. Individuals from the far northwestern range margin were separated from the remainder of the range by reduced migration, which was associated with environmental resistance. This portion of the range was invaded primarily in the last 15 years. Ecological niche models also indicated that this cluster was associated with the expansion of the niche. While invasion is often assumed to be primarily influenced by dispersal limitation, our results suggest that ongoing invasion and range shifts with climate change may be strongly affected by environmental heterogeneity.

Environmental and anthropogenic factors mediating the functional connectivity of the mesquite lizard along the eastern Trans-Mexican Volcanic Belt

Functional connectivity, the extent to which a landscape facilitates or impedes the dispersal of individuals across the landscape, is a key factor for the survival of species. Anthropogenic activities, such as urbanization, agriculture and roads, negatively impact functional connectivity of most species, particularly low-vagility species like lizards. Here, we examine how a landscape modified by anthropogenic activities affects the functional connectivity, at both broad and fine scales, of a widely distributed generalist lizard Sceloporus grammicus in the eastern Trans-Mexican Volcanic Belt, Mexico. We estimated for the first time the species' genetic structure, gene flow and functional connectivity in agricultural and forest zones using genomic data, a comprehensive landscape characterization and novel methods including gravity models. Our results showed not only marked genetic differentiation across the study region but also that functional connectivity is maintained for tens of kilometres despite S. grammicus low vagility. Specifically, we found that substrate and air temperature facilitated connectivity over broad and fine scales, respectively, while agricultural cover, relative humidity and slope were important for connectivity and gene flow. Contrastingly, forest cover and roads favoured (broad-scale) and limited (fine-scale) connectivity, likely associated with movement facilitated by small forest patches and with thermoregulation. Altogether, these results support that S. grammicus alternates its thermoregulatory behaviour depending on the distance travelled and the habitat environmental conditions, and that it can disperse through relatively modified landscapes, mainly using agricultural zones. The information obtained is crucial to understanding the response of lizards to current anthropogenic pressures and their potential to adapt.

Genome resequencing reveals population divergence and local adaptation of blacklegged ticks in the United States

Tick vectors and tick-borne disease are increasingly impacting human populations globally. An important challenge is to understand tick movement patterns, as this information can be used to improve management and predictive modelling of tick population dynamics. Evolutionary analysis of genetic divergence, gene flow and local adaptation provides insight on movement patterns at large spatiotemporal scales. We develop low coverage, whole genome resequencing data for 92 blacklegged ticks, Ixodes scapularis, representing range-wide variation across the United States. Through analysis of population genomic data, we find that tick populations are structured geographically, with gradual isolation by distance separating three population clusters in the northern United States, southeastern United States and a unique cluster represented by a sample from Tennessee. Populations in the northern United States underwent population contractions during the last glacial period and diverged from southern populations at least 50 thousand years ago. Genome scans of selection provide strong evidence of local adaptation at genes responding to host defences, blood-feeding and environmental variation. In addition, we explore the potential of low coverage genome sequencing of whole-tick samples for documenting the diversity of microbial pathogens and recover important tick-borne pathogens such as Borrelia burgdorferi. The combination of isolation by distance and local adaptation in blacklegged ticks demonstrates that gene flow, including recent expansion, is limited to geographical scales of a few hundred kilometres.

The coalescent in finite populations with arbitrary, fixed structure

The coalescent is a stochastic process representing ancestral lineages in a population undergoing neutral genetic drift. Originally defined for a well-mixed population, the coalescent has been adapted in various ways to accommodate spatial, age, and class structure, along with other features of real-world populations. To further extend the range of population structures to which coalescent theory applies, we formulate a coalescent process for a broad class of neutral drift models with arbitrary - but fixed - spatial, age, sex, and class structure, haploid or diploid genetics, and any fixed mating pattern. Here, the coalescent is represented as a random sequence of mappings [Formula: see text] from a finite set G to itself. The set G represents the "sites" (in individuals, in particular locations and/or classes) at which these alleles can live. The state of the coalescent, Ct:G→G, maps each site g∈G to the site containing g's ancestor, t time-steps into the past. Using this representation, we define and analyze coalescence time, coalescence branch length, mutations prior to coalescence, and stationary probabilities of identity-by-descent and identity-by-state. For low mutation, we provide a recipe for computing identity-by-descent and identity-by-state probabilities via the coalescent. Applying our results to a diploid population with arbitrary sex ratio r, we find that measures of genetic dissimilarity, among any set of sites, are scaled by 4r(1-r) relative to the even sex ratio case.

Some reef-building corals only disperse metres per generation

Understanding the dispersal potential of different species is essential for predicting recovery trajectories following local disturbances and the potential for adaptive loci to spread to populations facing extreme environmental changes. However, dispersal distances have been notoriously difficult to estimate for scleractinian corals, where sexually (as gametes or larvae) or asexually (as fragments or larvae) derived propagules disperse through vast oceans. Here, we demonstrate that generational dispersal distances for sexually produced propagules can be indirectly inferred for corals using individual-based isolation-by-distance (IbD) analyses by combining reduced-representation genomic sequencing with photogrammetric spatial mapping. Colonies from the genus Agaricia were densely sampled across plots at four locations and three depths in Curaçao. Seven cryptic taxa were found among the three nominal species (Agaricia agaricites, Agaricia humilis and Agaricia lamarcki), with four taxa showing generational dispersal distances within metres (two taxa within A. agaricites and two within A. humilis). However, no signals of IbD were found in A. lamarcki taxa and thus these taxa probably disperse relatively longer distances. The short distances estimated here imply that A. agaricites and A. humilis populations are reliant on highly localized replenishment and demonstrate the need to estimate dispersal distances quantitatively for more coral species.

On the connections between the spatial Lambda-Fleming-Viot model and other processes for analysing geo-referenced genetic data

The introduction of the spatial Lambda-Fleming-Viot model (ΛV) in population genetics was mainly driven by the pioneering work of Alison Etheridge, in collaboration with Nick Barton and Amandine Véber about ten years ago (Barton et al., 2010; Barton et al., 2013). The ΛV model provides a sound mathematical framework for describing the evolution of a population of related individuals along a spatial continuum. It alleviates the "pain in the torus" issue with Wright and Malécot's isolation by distance model and is sampling consistent, making it a tool of choice for statistical inference. Yet, little is known about the potential connections between the ΛV and other stochastic processes generating trees and the spatial coordinates along the corresponding lineages. This work focuses on a version of the ΛV whereby lineages move rapidly over small distances. Using simulations, we show that the induced ΛV tree-generating process is well approximated by a birth-death model. Our results also indicate that Brownian motions modelling the movements of lines of descent along birth-death trees do not generally provide a good approximation of the ΛV due to habitat boundaries effects that play an increasingly important role in the long run. Accounting for habitat boundaries through reflected Brownian motions considerably increases the similarity to the ΛV model however. Finally, we describe efficient algorithms for fast simulation of the backward and forward in time versions of the ΛV model.

Substrate specialisation drives an unexpectedly diverse radiation in barking geckos (Ptenopus: Gekkonidae)

Barking geckos (genus Ptenopus) are terrestrial, burrowing lizards endemic to southern Africa, currently with three recognised species. Two species are range-restricted (P. kochi and P. carpi) and display clear differences in substrate preference (soft sand vs. hard gravel). The third and most widespread species, P. garrulus, occurs on a variety of substrates of differing hardness, across potential geographic barriers, and over a steep climatic gradient. Variations in morphology and advertisement calls indicates that P. garrulus may be a species complex. Two subspecies of P. garrulus are currently recognised: P. g. maculatus and P. g. garrulus. To investigate species boundaries, we produced the first comprehensive phylogeny for the genus. We used a novel application of multiple regression on matrices models to assess multiple environmental drivers of diversification, as contrasted to isolation by distance. We show that P. kochi, P. carpi, and P. g. garrulus are valid species, but that P. g. maculatus is a paraphyletic complex of five previously unrecognised taxa. Specialisation onto different substrates was likely the main driver of divergence, with parapatric occurrence of two to four clades occurring at each of the three substrate transition zones identified a priori. The region encompasses diverse bioclimatic regions and potential geographic barriers, and these likely played a role in some divergence events.

Modelling the demographic history of human North African genomes points to a recent soft split divergence between populations

BACKGROUND

North African human populations present a complex demographic scenario due to the presence of an autochthonous genetic component and population substructure, plus extensive gene flow from the Middle East, Europe, and sub-Saharan Africa.

RESULTS

We conducted a comprehensive analysis of 364 genomes to construct detailed demographic models for the North African region, encompassing its two primary ethnic groups, the Arab and Amazigh populations. This was achieved through an Approximate Bayesian Computation with Deep Learning (ABC-DL) framework and a novel algorithm called Genetic Programming for Population Genetics (GP4PG). This innovative approach enabled us to effectively model intricate demographic scenarios, utilizing a subset of 16 whole genomes at > 30X coverage. The demographic model suggested by GP4PG exhibited a closer alignment with the observed data compared to the ABC-DL model. Both point to a back-to-Africa origin of North African individuals and a close relationship with Eurasian populations. Results support different origins for Amazigh and Arab populations, with Amazigh populations originating back in Epipaleolithic times, while GP4PG supports Arabization as the main source of Middle Eastern ancestry. The GP4PG model includes population substructure in surrounding populations (sub-Saharan Africa and Middle East) with continuous decaying gene flow after population split. Contrary to ABC-DL, the best GP4PG model does not require pulses of admixture from surrounding populations into North Africa pointing to soft splits as drivers of divergence in North Africa.

CONCLUSIONS

We have built a demographic model on North Africa that points to a back-to-Africa expansion and a differential origin between Arab and Amazigh populations.

Signatures of selective sweeps in continuous-space populations

Selective sweeps describe the process by which an adaptive mutation arises and rapidly fixes in the population, thereby removing genetic variation in its genomic vicinity. The expected signatures of selective sweeps are relatively well understood in panmictic population models, yet natural populations often extend across larger geographic ranges where individuals are more likely to mate with those born nearby. To investigate how such spatial population structure can affect sweep dynamics and signatures, we simulated selective sweeps in populations inhabiting a two-dimensional continuous landscape. The maximum dispersal distance of offspring from their parents can be varied in our simulations from an essentially panmictic population to scenarios with increasingly limited dispersal. We find that in low-dispersal populations, adaptive mutations spread more slowly than in panmictic ones, while recombination becomes less effective at breaking up genetic linkage around the sweep locus. Together, these factors result in a trough of reduced genetic diversity around the sweep locus that looks very similar across dispersal rates. We also find that the site frequency spectrum around hard sweeps in low-dispersal populations becomes enriched for intermediate-frequency variants, making these sweeps appear softer than they are. Furthermore, haplotype heterozygosity at the sweep locus tends to be elevated in low-dispersal scenarios as compared to panmixia, contrary to what we observe in neutral scenarios without sweeps. The haplotype patterns generated by these hard sweeps in low-dispersal populations can resemble soft sweeps from standing genetic variation that arose from substantially older alleles. Our results highlight the need for better accounting for spatial population structure when making inferences about selective sweeps.

Pathogenicity is associated with population structure in a fungal pathogen of humans

Aspergillus flavus is a clinically and agriculturally important saprotrophic fungus responsible for severe human infections and extensive crop losses. We analyzed genomic data from 250 (95 clinical and 155 environmental) A. flavus isolates from 9 countries, including 70 newly sequenced clinical isolates, to examine population and pan-genome structure and their relationship to pathogenicity. We identified five A. flavus populations, including a new population, D, corresponding to distinct clades in the genome-wide phylogeny. Strikingly, > 75% of clinical isolates were from population D. Accessory genes, including genes within biosynthetic gene clusters, were significantly more common in some populations but rare in others. Population D was enriched for genes associated with zinc ion binding, lipid metabolism, and certain types of hydrolase activity. In contrast to the major human pathogen Aspergillus fumigatus, A. flavus pathogenicity in humans is strongly associated with population structure, making it a great system for investigating how population-specific genes contribute to pathogenicity.

Estimating scale-specific and localized spatial patterns in allele frequency

Characterizing spatial patterns in allele frequencies is fundamental to evolutionary biology because these patterns contain evidence of underlying processes. However, the spatial scales at which gene flow, changing selection, and drift act are often unknown. Many of these processes can operate inconsistently across space, causing nonstationary patterns. We present a wavelet approach to characterize spatial pattern in allele frequency that helps solve these problems. We show how our approach can characterize spatial patterns in relatedness at multiple spatial scales, i.e. a multilocus wavelet genetic dissimilarity. We also develop wavelet tests of spatial differentiation in allele frequency and quantitative trait loci (QTL). With simulation, we illustrate these methods under different scenarios. We also apply our approach to natural populations of Arabidopsis thaliana to characterize population structure and identify locally adapted loci across scales. We find, for example, that Arabidopsis flowering time QTL show significantly elevated genetic differentiation at 300-1,300 km scales. Wavelet transforms of allele frequencies offer a flexible way to reveal geographic patterns and underlying evolutionary processes.

Unraveling the genomic landscape of Campylorhynchus wrens along western Ecuador's precipitation gradient: Insights into hybridization, isolation by distance, and isolation by the environment

Environmental gradients have the potential to influence genetic differentiation among populations ultimately leading to allopatric speciation. However, environmental gradients can also facilitate hybridization between closely related taxa. We investigated a putative hybrid zone in western Ecuador, involving two polytypic wren species (Aves: Troglodytidae), Campylorhynchus zonatus and C. fasciatus. Our study addressed two primary questions: (1) Is there evidence of population structure and genetic admixture between these taxa in western Ecuador? and (2) What are the relative contributions of isolation by distance and isolation by the environment to the observed genetic differentiation along the environmental gradient in this region? We analyzed 4409 single-nucleotide polymorphisms (SNPs) from 112 blood samples sequenced using ddRadSeq and a de novo assembly. The optimum number of genetic clusters ranged from 2 to 4, aligning with geographic origins, known phylogenetics, and physical or ecological constraints. We observed notable transitions in admixture proportions along the environmental gradient in western Ecuador between C. z. brevirostris and the northern and southern genetic clusters of C. f. pallescens. Genetic differentiation between the two C. f. pallescens populations could be attributed to an unreported potential physical barrier in central western Ecuador, where the proximity of the Andes to the coastline restricts lowland habitats, limiting dispersal and gene flow, especially among dry-habitat specialists. The observed admixture in C. f. pallescens suggests that this subspecies may be a hybrid between C. z. brevirostris and C. fasciatus, with varying degrees of admixture in western Ecuador and northwestern Peru. We found evidence of isolation by distance, while isolation by the environment was less pronounced but still significant for annual mean precipitation and precipitation seasonality. This study enhances our understanding of avian population genomics in tropical regions.

Neutral processes and taxonomic scale drive beta species-genetic diversity correlations in a submesophotic tropical reef fish

If similar evolutionary forces maintain intra- and interspecific diversity, patterns of diversity at both levels of biological organization can be expected to covary across space. Although this prediction of a positive species-genetic diversity correlation (SGDC) has been tested for several taxa in natural landscapes, no study has yet evaluated the influence of the community delineation on these SGDCs. In this study, we focused on tropical fishes of the Indo-Pacific Ocean, using range-wide single nucleotide polymorphism data for a deep-sea fish (Etelis coruscans) and species presence data of 4878 Teleostei species. We investigated whether a diversity continuum occurred, for different community delineations (subfamily, family, order and class) and spatial extents, and which processes explained these diversity patterns. We found no association between genetic diversity and species richness (α-SGDC), regardless of the community and spatial extent. In contrast, we evidenced a positive relationship between genetic and species dissimilarities (β-SGDC) when the community was defined at the subfamily or family level of the species of interest, and when the Western Indian Ocean was excluded. This relationship was related to the imprint of dispersal processes across levels of biological organization in Lutjanidae. However, this positive β-SGDC was lost when considering higher taxonomic communities and at the scale of the entire Indo-Pacific, suggesting different responses of populations and communities to evolutionary processes at these scales. This study provides evidence that the taxonomic scale at which communities are defined and the spatial extent are pivotal to better understand the processes shaping diversity across levels of biological organization.

City divided: Unveiling family ties and genetic structuring of coyotes in Seattle

Linear barriers pose significant challenges for wildlife gene flow, impacting species persistence, adaptation, and evolution. While numerous studies have examined the effects of linear barriers (e.g., fences and roadways) on partitioning urban and non-urban areas, understanding their influence on gene flow within cities remains limited. Here, we investigated the impact of linear barriers on coyote (Canis latrans) population structure in Seattle, Washington, where major barriers (i.e., interstate highways and bodies of water) divide the city into distinct quadrants. Just under 1000 scats were collected to obtain genetic data between January 2021 and December 2022, allowing us to identify 73 individual coyotes. Notably, private allele analysis underscored limited interbreeding among quadrants. When comparing one quadrant to each other, there were up to 16 private alleles within a single quadrant, representing nearly 22% of the population allelic diversity. Our analysis revealed weak isolation by distance, and despite being a highly mobile species, genetic structuring was apparent between quadrants even with extremely short geographic distance between individual coyotes, implying that Interstate 5 and the Ship Canal act as major barriers. This study uses coyotes as a model species for understanding urban gene flow and its consequences in cities, a crucial component for bolstering conservation of rarer species and developing wildlife friendly cities.

Neutral markers reveal complex population structure across the range of a widespread songbird

Understanding how both contemporary and historical physical barriers influence gene flow is key to reconstructing evolutionary histories and can allow us to predict species' resilience to changing environmental conditions. During the last glacial maximum (LGM), many high latitude North American bird species were forced into glacial refugia, including mountain bluebirds (Silia currucoides). Within their current breeding range, mountain bluebirds still experience a wide variety of environmental conditions and barriers that may disrupt gene flow and isolate populations. Using single nucleotide polymorphisms (SNPs) obtained through restriction site-associated DNA sequencing, we detected at least four genetically distinct mountain bluebird populations. Based on this structure, we determined that isolation-by-distance, the northern Rocky Mountains, and discontinuous habitat are responsible for the low connectivity and the overall history of each population going back to the last glacial maximum. Finally, we identified five candidate genes under balancing selection and three loci under diversifying selection. This study provides the first look at connectivity and gene flow across the range of these high-altitude and high latitude songbirds.

Centrality to the metapopulation is more important for population genetic diversity than habitat area or fragmentation

Drift and gene flow affect genetic diversity. Given that the strength of genetic drift increases as population size decreases, management activities have focused on increasing population size through preserving habitats to preserve genetic diversity. Few studies have empirically evaluated the impacts of drift and gene flow on genetic diversity. Kryptolebias marmoratus, henceforth 'rivulus', is a small killifish restricted to fragmented New World mangrove forests with gene flow primarily associated with ocean currents. Rivulus form distinct populations across patches, making them a well-suited system to test the extent to which habitat area, fragmentation and connectivity are associated with genetic diversity. Using over 1000 individuals genotyped at 32 microsatellite loci, high-resolution landcover data and oceanographic simulations with graph theory, we demonstrate that centrality (connectivity) to the metapopulation is more strongly associated with genetic diversity than habitat area or fragmentation. By comparing models with and without centrality standardized by the source population's genetic diversity, our results suggest that metapopulation centrality is critical to genetic diversity regardless of the diversity of adjacent populations. While we find evidence that habitat area and fragmentation are related to genetic diversity, centrality is always a significant predictor with a larger effect than any measure of habitat configuration.

Geoclimatic drivers of diversification in the largest arid and semi-arid environment of the Neotropics: Perspectives from phylogeography

The South American Dry Diagonal, also called the Diagonal of Open Formations, is a large region of seasonally dry vegetation extending from northeastern Brazil to northern Argentina, comprising the Caatinga, Cerrado, and Chaco subregions. A growing body of phylogeography literature has determined that a complex history of climatic changes coupled with more ancient geological events has produced a diverse and endemic-rich Dry Diagonal biota. However, the exact drivers are still under investigation, and their relative strengths and effects are controversial. Pleistocene climatic fluctuations structured lineages via vegetation shifts, refugium formation, and corridors between the Amazon and Atlantic forests. In some taxa, older geological events, such as the reconfiguration of the São Francisco River, uplift of the Central Brazilian Plateau, or the Miocene inundation of the Chaco by marine incursions, were more important. Here, we review the Dry Diagonal phylogeography literature, discussing each hypothesized driver of diversification and assessing degree of support. Few studies statistically test these hypotheses, with most support drawn from associating encountered phylogeographic patterns such as population structure with the timing of ancient geoclimatic events. Across statistical studies, most hypotheses are well supported, with the exception of the Pleistocene Arc Hypothesis. However, taxonomic and regional biases persist, such as a proportional overabundance of herpetofauna studies, and the under-representation of Chaco studies. Overall, both Pleistocene climate change and Neogene geological events shaped the evolution of the Dry Diagonal biota, though the precise effects are regionally and taxonomically varied. We encourage further use of model-based analyses to test evolutionary scenarios, as well as interdisciplinary collaborations to progress the field beyond its current focus on the traditional set of geoclimatic hypotheses.

Niche modelling and landscape genetics of the yellow-legged hornet (Vespa velutina): An integrative approach for evaluating central-marginal population dynamics in Europe

Genetic diversity is an important biological trait for a successful invasion. During the expansion across a new territory, an invasive species may face unprecedented ecological conditions that will determine its demography and genetic diversity. The first record of the yellow-legged hornet (Vespa velutina) in Europe dates back to 2004 in France, from where it has successfully spread through a large territory in the continent, including Italy, Spain and Portugal. Integrative approaches offer a powerful strategy to detect and understand patterns of genetic variation in central and marginal populations. Here, we have analysed the relationship between genetic diversity parameters inferred from 15 V. velutina nuclear DNA microsatellite loci, and geographical and environmental drivers, such as the distance to the introduction focus, environmental suitability and distance to native and invasive niche centroids. Our results revealed a central-marginal dynamic, where allelic richness decreased towards the edge of the expansion range. The low environmental suitability of the territories invaded by marginal populations could prevent a diverse population from establishing and reducing the genetic diversity in populations at the expansion edge. Moreover, Markov chain Monte Carlo analysis showed both geographical and environmental distances were influencing population genetic differentiation. This study highlights the importance of combining genetic analysis with geographical and environmental drivers to understand genetic trends of invasive species to new environment.

Limits to species' range: the tension between local and global adaptation

We know that heritable variation is abundant, and that selection causes all but the smallest populations to rapidly shift beyond their original trait distribution. So then, what limits the range of a species? There are physical constraints and also population genetic limits to the effectiveness of selection, ultimately set by population size. Global adaptation, where the same genotype is favoured over the whole range, is most efficient when based on a multitude of weakly selected alleles and is effective even when local demes are small, provided that there is some gene flow. In contrast, local adaptation is sensitive to gene flow and may require alleles with substantial effect. How can populations combine the advantages of large effective size with the ability to specialise into local niches? To what extent does reproductive isolation help resolve this tension? I address these questions using eco-evolutionary models of polygenic adaptation, contrasting discrete demes with continuousspace.

Environment influences the genetic structure and genetic differentiation of Sassafras tzumu (Lauraceae)

BACKGROUND

Sassafras tzumu, an elegant deciduous arboreal species, belongs to the esteemed genus Sassafras within the distinguished family Lauraceae. With its immense commercial value, escalating market demands and unforeseen human activities within its natural habitat have emerged as new threats to S. tzumu in recent decades, so it is necessary to study its genetic diversity and influencing factors, to propose correlative conservation strategies.

RESULTS

By utilizing genotyping-by-sequence (GBS) technology, we acquired a comprehensive database of single nucleotide polymorphisms (SNPs) from a cohort of 106 individuals sourced from 13 diverse Sassafras tzumu natural populations, scattered across various Chinese mountainous regions. Through our meticulous analysis, we aimed to unravel the intricate genetic diversity and structure within these S. tzumu populations, while simultaneously investigating the various factors that potentially shape genetic distance. Our preliminary findings unveiled a moderate level of genetic differentiation (FST = 0.103, p < 0.01), accompanied by a reasonably high genetic diversity among the S. tzumu populations. Encouragingly, our principal component analysis painted a vivid picture of two distinct genetic and geographical regions across China, where gene flow appeared to be somewhat restricted. Furthermore, employing the sophisticated multiple matrix regression with randomization (MMRR) analysis method, we successfully ascertained that environmental distance exerted a more pronounced impact on genetic distance when compared to geographical distance (βE = 0.46, p < 0.01; βD = 0.16, p < 0.01). This intriguing discovery underscores the potential significance of environmental factors in shaping the genetic landscape of S. tzumu populations.

CONCLUSIONS

The genetic variance among populations of S. tzumu in our investigation exhibited a moderate degree of differentiation, alongside a heightened level of genetic diversity. The environmental distance of S. tzumu had a greater impact on its genetic diversity than geographical distance. It is of utmost significance to formulate and implement meticulous management and conservation strategies to safeguard the invaluable genetic resources of S. tzumu.

Evolutionary graph theory beyond single mutation dynamics: on how network-structured populations cross fitness landscapes

Spatially resolved datasets are revolutionizing knowledge in molecular biology, yet are under-utilized for questions in evolutionary biology. To gain insight from these large-scale datasets of spatial organization, we need mathematical representations and modeling techniques that can both capture their complexity, but also allow for mathematical tractability. Evolutionary graph theory utilizes the mathematical representation of networks as a proxy for heterogeneous population structure and has started to reshape our understanding of how spatial structure can direct evolutionary dynamics. However, previous results are derived for the case of a single new mutation appearing in the population and the role of network structure in shaping fitness landscape crossing is still poorly understood. Here we study how network-structured populations cross fitness landscapes and show that even a simple extension to a two-mutational landscape can exhibit complex evolutionary dynamics that cannot be predicted using previous single-mutation results. We show how our results can be intuitively understood through the lens of how the two main evolutionary properties of a network, the amplification and acceleration factors, change the expected fate of the intermediate mutant in the population and further discuss how to link these models to spatially resolved datasets of cellular organization.

Synthesizing environmental, epidemiological and vector and parasite genetic data to assist decision making for disease elimination

We present a framework for identifying when conditions are favourable for transmission of vector-borne diseases between communities by incorporating predicted disease prevalence mapping with landscape analysis of sociological, environmental and host/parasite genetic data. We explored the relationship between environmental features and gene flow of a filarial parasite of humans, Onchocerca volvulus, and its vector, blackflies in the genus Simulium. We generated a baseline microfilarial prevalence map from point estimates from 47 locations in the ecological transition separating the savannah and forest in Ghana, where transmission of O. volvulus persists despite onchocerciasis control efforts. We generated movement suitability maps based on environmental correlates with mitochondrial population structure of 164 parasites from 15 communities and 93 vectors from only four sampling sites, and compared these to the baseline prevalence map. Parasite genetic distance between sampling locations was significantly associated with elevation (r = .793, p = .005) and soil moisture (r = .507, p = .002), while vector genetic distance was associated with soil moisture (r = .788, p = .0417) and precipitation (r = .835, p = .0417). The correlation between baseline prevalence and parasite resistance surface maps was stronger than that between prevalence and vector resistance surface maps. The centre of the study area had high prevalence and suitability for parasite and vector gene flow, potentially contributing to persistent transmission and suggesting the importance of re-evaluating transmission zone boundaries. With suitably dense sampling, this framework can help delineate transmission zones for onchocerciasis and would be translatable to other vector-borne diseases.

Effective dispersal patterns in prairie plant species across human-modified landscapes

Effective dispersal among plant populations is dependent on vector behaviour, landscape features and availability of adequate habitats. To capture landscape feature effects on dispersal, studies must be conducted at scales reflecting single-generation dispersal events (mesoscale). Many studies are conducted at large scales where genetic differentiation is due to dispersal occurring over multiple generations, making it difficult to interpret the effects of specific landscape features on vector behaviour. Genetic structure at the mesoscale may be determined by ecological and evolutionary processes, such as the consequences of vector behaviour on patterns of gene flow. We used chloroplast haplotypes and nuclear genome SNP surveys to identify landscape features influencing seed and pollen dispersal at a mesoscale within the Rogue River Valley in southern Oregon. We evaluated biotic and abiotic vector behaviour by contrasting two annual species with differing dispersal mechanisms; Achyrachaena mollis (Asteraceae) is a self-pollinating and anemochoric species, and Plectritis congesta (Caprifoliaceae) is biotically pollinated with barochoric seeds. Using landscape genetics methods, we identified features of the study region that conduct or restrict dispersal. We found chloroplast haplotypes were indicative of historic patterns of gene flow prior to human modification of landscapes. Seed dispersal of A. mollis was best supported by models of isolation by distance, while seed-driven gene flow of P. congesta was determined by the distribution of preserved natural spaces and quality habitat. Nuclear genetic structure was driven by both pollen and seed dispersal, and both species responded to contemporary landscape changes, such as urban and agricultural conversion, and habitat availability.