Search

Large‐scale chromosome rearrangements, such as fissions and fusions, are a common feature of eukaryote evolution. They can have considerable influence on the evolution of populations, yet it remains unclear exactly how rearrangements become established and eventually fix. Rearrangements could fix by genetic drift if they are weakly deleterious or neutral, or they may instead be favoured by positive natural selection. Here, we compare genome assemblies of three closely related Brenthis butterfly species and characterize a complex history of fission and fusion rearrangements. An inferred demographic history of these species suggests that rearrangements became fixed in populations with large long‐term effective size (Ne), consistent with rearrangements being selectively neutral or only very weakly underdominant. Using a recently developed analytic framework for characterizing hard selective sweeps, we find that chromosome fusions are not enriched for evidence of past sweeps compared to other regions of the genome. Nonetheless, we do infer a strong and recent selective sweep around one chromosome fusion in the B. daphne genome. Our results suggest that rearrangements in these species likely have weak absolute fitness effects and fix by genetic drift. However, one putative selective sweep raises the possibility that natural selection may sometimes play a role in the fixation of chromosome fusions. [ABSTRACT FROM AUTHOR]

Molecular methods are routinely used to estimate the effective size of populations (Ne). However, underlying model assumptions are frequently violated to an unknown extent. Although simulations can detect sources of bias and help to adjust sampling strategies and analyses methods, additional information from empirical data can also be used to calibrate methods and improve molecular Ne estimation methods. Here, we take advantage of long‐term genetic and ecological monitoring data of the grey wolf (Canis lupus) in Germany, and detailed population genetic studies in Poland, Spain and Portugal to improve Ne estimation strategies in this species, and species with similar life history traits. We first calculated Ne from average lifetime reproductive success and detailed census data from the German population, which served as a baseline to compare to molecular estimates based on linkage disequilibrium and sibship frequency. This yielded a robust Ne/Nc estimation that we used to calibrate molecular estimates of German, Polish and Iberian wolf populations. The linkage disequilibrium method was strongly influenced by spatial genetic structure, much more than the sibship frequency method. When Ne was estimated in local neighbourhoods, both methods yielded comparable results. Estimates of the metapopulation effective size seemed to correspond generally well with the sum of the estimates of local neighbourhoods. Overall, we found that the number of packs is a good proxy of the effective population size. Using this as a rule of thumb, we evaluated for all European wolf populations the Ne 500 indicator and concluded that half of the European wolf populations do not yet fulfil this criterion. [ABSTRACT FROM AUTHOR]

High‐density genotyping methods have revolutionized the field of population and conservation genetics in the past decade. To exploit the technological and analytical advances in the field, access to high‐quality genetic material is a key component. However, access to such samples in endangered and rare animals is often challenging or even impossible. Here, we used a minimally invasive sampling method (MIS) in the endangered cave salamander Proteus anguinus, the olm, to generate thousands of genetic markers using ddRADseq for population and conservation genomic analyses. Using tail clips and MIS skin swabs taken from the same individual, we investigated genotyping data properties of the two different sampling types. We found that sufficient DNA can be extracted from swab samples to generate up to 200,000 polymorphic SNPs in divergent Proteus lineages. Swab and tissue samples were highly reproducible exhibiting low SNP genotyping error rates. We found that SNPs were most frequently (~50%) located within genic regions, while the rest mapped to mostly flanking regions of repetitive DNA. The vast majority of DNA recovered from swabbing was host DNA. However, a fraction of DNA recovered from swabs contained additional ecological information on the species, including eDNA from the surrounding environment and bacterial skin fauna. Most exogenous DNA recovered from swabs were bacteria (~80%), followed by vertebrates (~20%). Our results demonstrate that MIS can be used to (i) generate tens of thousands of ddRADseq markers for conservation and population genomic analyses and (ii) inform on the species health status and ecology from exogenous DNA. [ABSTRACT FROM AUTHOR]

Knowledge about sex‐specific difference in life‐history traits—like growth, mortality, or behavior—is of key importance for management and conservation as these parameters are essential for predictive modeling of population sustainability. We applied a newly developed molecular sex identification method, in combination with a SNP (single nucleotide polymorphism) panel for inferring the population of origin, for more than 300 large Atlantic bluefin tuna (ABFT) collected over several years from newly reclaimed feeding grounds in the Northeast Atlantic. The vast majority (95%) of individuals were genetically assigned to the eastern Atlantic population, which migrates between spawning grounds in the Mediterranean and feeding grounds in the Northeast Atlantic. We found a consistent pattern of a male bias among the eastern Atlantic individuals, with a 4‐year mean of 63% males (59%–65%). Males were most prominent within the smallest (< 230 cm) and largest (> 250 cm) length classes, while the sex ratio was close to 1:1 for intermediate sizes (230–250 cm). The results from this new, widely applicable, and noninvasive approach suggests differential occupancy or migration timing of ABFT males and females, which cannot be explained alone by sex‐specific differences in growth. Our findings are corroborated by previous traditional studies of sex ratios in dead ABFT from the Atlantic, the Mediterranean, and the Gulf of Mexico. In concert with observed differences in growth and mortality rates between the sexes, these findings should be recognized in order to sustainably manage the resource, maintain productivity, and conserve diversity within the species. [ABSTRACT FROM AUTHOR]

Adaptive evolution can facilitate species' range expansions across environmentally heterogeneous landscapes. However, serial founder effects can limit the efficacy of selection, and the evolution of increased dispersal during range expansions may result in gene flow swamping local adaptation. Here, we study how genetic drift, gene flow and selection interact during the cane toad's (Rhinella marina) invasion across the heterogeneous landscape of Australia. Following its introduction in 1935, the cane toad colonised eastern Australia and established several stable range edges. The ongoing, more rapid range expansion in north‐central Australia has occurred concomitant with an evolved increase in dispersal capacity. Using reduced representation genomic data of Australian cane toads from the expansion front and from two areas of their established range, we test the hypothesis that high gene flow constrains local adaptation at the expansion front relative to established areas. Genetic analyses indicate the three study areas are genetically distinct but show similar levels of allelic richness, heterozygosity and inbreeding. Markedly higher gene flow or recency of colonisation at the expansion front have likely hindered local adaptation at the time of sampling, as indicated by reduced slopes of genetic‐environment associations (GEAs) estimated using a novel application of geographically weighted regression that accounts for allele surfing; GEA slopes are significantly steeper in established parts of the range. Our work bolsters evidence supporting adaptation of invasive species post‐introduction and adds novel evidence for differing strengths of evolutionary forces among geographic areas with different invasion histories. [ABSTRACT FROM AUTHOR]

Human‐mediated habitat destruction has had a profound impact on increased species extinction rates and population declines worldwide. The coastal development in the United Arab Emirates (UAE) over the last two decades, serves as an example of how habitat transformation can alter the landscape of a country in just a few years. Here, we study the genomic implications of habitat transformation in the Critically Endangered Emirati Leaf‐toed Gecko (Asaccus caudivolvulus), the only endemic vertebrate of the UAE. We generate a high‐quality reference genome for this gecko, representing the first reference genome for the family Phyllodactylidae, and produce whole‐genome resequencing data for 23 specimens from 10 different species of leaf‐toed geckos. Our results show that A. caudivolvulus has consistently lower genetic diversity than any other Arabian species of Asaccus, suggesting a history of ancient population declines. However, high levels of recent inbreeding are recorded among populations in heavily developed areas, with a more than 50% increase in long runs of homozygosity within a 9‐year period. Moreover, results suggest that this species does not effectively purge deleterious mutations, hence making it more vulnerable to future stochastic threats. Overall, results show that A. caudivolvulus is in urgent need of protection, and habitat preservation must be warranted to ensure the species' survival. [ABSTRACT FROM AUTHOR]

A major aim of evolutionary biology is to understand why patterns of genomic diversity vary within taxa and space. Large‐scale genomic studies of widespread species are useful for studying how environment and demography shape patterns of genomic divergence. Here, we describe one of the most geographically comprehensive surveys of genomic variation in a wild vertebrate to date; the great tit (Parus major) HapMap project. We screened ca 500,000 SNP markers across 647 individuals from 29 populations, spanning ~30 degrees of latitude and 40 degrees of longitude – almost the entire geographical range of the European subspecies. Genome‐wide variation was consistent with a recent colonisation across Europe from a South‐East European refugium, with bottlenecks and reduced genetic diversity in island populations. Differentiation across the genome was highly heterogeneous, with clear 'islands of differentiation', even among populations with very low levels of genome‐wide differentiation. Low local recombination rates were a strong predictor of high local genomic differentiation (FST), especially in island and peripheral mainland populations, suggesting that the interplay between genetic drift and recombination causes highly heterogeneous differentiation landscapes. We also detected genomic outlier regions that were confined to one or more peripheral great tit populations, probably as a result of recent directional selection at the species' range edges. Haplotype‐based measures of selection were related to recombination rate, albeit less strongly, and highlighted population‐specific sweeps that likely resulted from positive selection. Our study highlights how comprehensive screens of genomic variation in wild organisms can provide unique insights into spatio‐temporal evolutionary dynamics. [ABSTRACT FROM AUTHOR]

ABSTRACT Molecular methods are routinely used to estimate the effective size of populations (Ne). However, underlying model assumptions are frequently violated to an unknown extent. Although simulations can detect sources of bias and help to adjust sampling strategies and analyses methods, additional information from empirical data can also be used to calibrate methods and improve molecular Ne estimation methods. Here, we take advantage of long‐term genetic and ecological monitoring data of the grey wolf (Canis lupus) in Germany, and detailed population genetic studies in Poland, Spain and Portugal to improve Ne estimation strategies in this species, and species with similar life history traits. We first calculated Ne from average lifetime reproductive success and detailed census data from the German population, which served as a baseline to compare to molecular estimates based on linkage disequilibrium and sibship frequency. This yielded a robust Ne/Nc estimation that we used to calibrate molecular estimates of German, Polish and Iberian wolf populations. The linkage disequilibrium method was strongly influenced by spatial genetic structure, much more than the sibship frequency method. When Ne was estimated in local neighbourhoods, both methods yielded comparable results. Estimates of the metapopulation effective size seemed to correspond generally well with the sum of the estimates of local neighbourhoods. Overall, we found that the number of packs is a good proxy of the effective population size. Using this as a rule of thumb, we evaluated for all European wolf populations the Ne 500 indicator and concluded that half of the European wolf populations do not yet fulfil this criterion.

Environmental stress is a fundamental facet of life and a significant driver of natural selection in the wild. Gene expression diversity may facilitate adaptation to environmental changes, without necessary genetic change, but its role in adaptive divergence remains largely understudied in Neotropical systems. In Amazonian riparian forests, species distribution is predominantly influenced by species' waterlogging tolerance. The flooding gradient delineates distinct wetland forest types, shaping habitats and species characteristics. Here we investigated the molecular basis of environmental stress response in a tropical ground‐herb species (Ischnosiphon puberulus) to environmental variation in Amazonian riparian forests. We compared environmental variables and gene expression profiles from individuals collected in two forest types: Igapó and Terra firme in the Amazonian riparian forests. Predictable seasonal flooding poses a significant challenge in Igapó compared to Terra firme environments, with the former presenting higher water column height and longer flooding duration. Our findings suggest that contrasting environmental conditions related to flooding regimes are important drivers of population genetic differentiation and differential gene expression in I. puberulus. Enriched gene ontology terms highlight associations with environmental stresses, such as defence response, water transport, phosphorylation, root development, response to auxin, salicylic acid and oxidative stress. By uncovering key environmental stress response pathways conserved across populations, I. puberulus offers novel genetic insights into the molecular basis of plant reactions to environmental constraints found in flooded areas of this highly biodiverse neotropical ecosystem. [ABSTRACT FROM AUTHOR]

Wildlife populations are becoming increasingly fragmented by anthropogenic development. Small and isolated populations often face an elevated risk of extinction, in part due to inbreeding depression. Here, we examine the genomic consequences of urbanization in a caracal (Caracal caracal) population that has become isolated in the Cape Peninsula region of the City of Cape Town, South Africa, and is thought to number ~50 individuals. We document low levels of migration into the population over the past ~75 years, with an estimated rate of 1.3 effective migrants per generation. As a consequence of this isolation and small population size, levels of inbreeding are elevated in the contemporary Cape Peninsula population (mean FROH = 0.20). Inbreeding primarily manifests as long runs of homozygosity >10 Mb, consistent with the effects of isolation due to the rapid recent growth of Cape Town. To explore how reduced migration and elevated inbreeding may impact future population dynamics, we parameterized an eco‐evolutionary simulation model. We find that if migration rates do not change in the future, the population is expected to decline, though with a low projected risk of extinction. However, if migration rates decline or anthropogenic mortality rates increase, the potential risk of extinction is greatly elevated. To avert a population decline, we suggest that translocating migrants into the Cape Peninsula to initiate a genetic rescue may be warranted in the near future. Our analysis highlights the utility of genomic datasets coupled with computational simulation models for investigating the influence of gene flow on population viability. [ABSTRACT FROM AUTHOR]

ABSTRACT Knowledge about sex‐specific difference in life‐history traits—like growth, mortality, or behavior—is of key importance for management and conservation as these parameters are essential for predictive modeling of population sustainability. We applied a newly developed molecular sex identification method, in combination with a SNP (single nucleotide polymorphism) panel for inferring the population of origin, for more than 300 large Atlantic bluefin tuna (ABFT) collected over several years from newly reclaimed feeding grounds in the Northeast Atlantic. The vast majority (95%) of individuals were genetically assigned to the eastern Atlantic population, which migrates between spawning grounds in the Mediterranean and feeding grounds in the Northeast Atlantic. We found a consistent pattern of a male bias among the eastern Atlantic individuals, with a 4‐year mean of 63% males (59%–65%). Males were most prominent within the smallest ( 250 cm) length classes, while the sex ratio was close to 1:1 for intermediate sizes (230–250 cm). The results from this new, widely applicable, and noninvasive approach suggests differential occupancy or migration timing of ABFT males and females, which cannot be explained alone by sex‐specific differences in growth. Our findings are corroborated by previous traditional studies of sex ratios in dead ABFT from the Atlantic, the Mediterranean, and the Gulf of Mexico. In concert with observed differences in growth and mortality rates between the sexes, these findings should be recognized in order to sustainably manage the resource, maintain productivity, and conserve diversity within the species.

Anthropogenic impact is causing the decline of a large proportion of species worldwide and reduces their genetic diversity. Island species typically have smaller ranges than continental species. As a consequence, island species are particularly liable to undergo population bottlenecks, giving rise to conservation challenges such as inbreeding and unmasking of deleterious genetic load. Such challenges call for more detailed assessments of the genetic make‐up of threatened island populations. The Mascarene islands (Indian Ocean) present many prime examples, being unusual in having been pristine until first human arrival ~400 years ago, following which anthropogenic pressure was unusually intense. A threatened harrier (Circus maillardi) endemic to the westernmost island of the archipelago is a good example of the challenges faced by species that have declined to small population size following intense anthropogenic pressure. In this study, we use an extensive set of population genomic tools to quantify variation at near‐neutral and coding loci, in order to test the historical impact of human activity on this species, and evaluate the species' (mal)adaptive potential. We observed low but significant genetic differentiation between populations on the West and North‐East sides of the island, echoing observations in other endemic species. Inbreeding was significant, with a substantial fraction of samples being first or second‐degree relatives. Historical effective population sizes have declined from ~3000 to 300 individuals in the past 1000 years, with a more recent drop ~100 years ago consistent with human activity. Based on our simulations and comparisons with a close relative (Circus melanoleucos), this demographic history may have allowed purging of the most deleterious variants but is unlikely to have allowed the purging of mildly deleterious variants. Our study shows how using relatively affordable methods can reveal the massive impact that human activity may have on the genetic diversity and adaptive potential of island populations, and calls for urgent action to closely monitor the reproductive success of such endemic populations, in association with genetic studies. [ABSTRACT FROM AUTHOR]

Inbreeding depression, that is, the reduction of health and vigour in individuals with high inbreeding coefficients, is expected to increase with environmental, social, or physiological stress. It has therefore been predicted that sexual selection and the associated stress usually lead to higher inbreeding depression in males than in females. However, sex‐specific differences in life history may reverse that pattern during certain developmental stages. In some salmonids, for example, female juveniles start developing their gonads earlier than males who instead grow faster. We tested whether the sexes are differently affected by inbreeding during that time. To study the effects of inbreeding coefficients that may be typical for natural populations of brown trout (Salmo trutta), and also to control for potentially confounding maternal or paternal effects, we sampled males and females from the wild, used their gametes in a block‐wise full‐factorial breeding design to produce 60 full‐sib families, released the offspring as yolk‐sac larvae into the wild, sampled them 6 months later, identified their genetic sex, and used microsatellites to assign them to their parents. We used whole‐genome resequencing to calculate the kinship coefficients for each breeding pair and hence the expected average inbreeding coefficient per family. Juvenile growth could be predicted from these expected inbreeding coefficients and the genetic sex: Females reached lower body sizes with increasing inbreeding coefficient, while no such link could be found in males. This sex‐specific inbreeding depression led to the overall pattern that females were on average smaller than males by the end of their first summer. [ABSTRACT FROM AUTHOR]

Modern fisheries management strives to balance opposing goals of protection for weak stocks and opportunity for harvesting healthy stocks. Test fisheries can aid management of anadromous fishes if they can forecast the strength and timing of an annual run with adequate time to allow fisheries planning. Integration of genetic stock identification (GSI) can further maximize utility of test fisheries by resolving run forecasts into weak‐ and healthy‐stock subcomponents. Using 5 years (2017–2022) of test fishery data, our study evaluated accuracy, resolution, and lead time of predictions for stock‐specific run timing and abundance of Columbia River spring Chinook salmon (Oncorhynchus tshawytscha). We determined if this test fishery (1) could use visual stock identification (VSI) to forecast at the coarse stock resolution (i.e., classification of "lower" vs. "upriver" stocks) upon which current management is based and (2) could be enhanced with GSI to forecast at higher stock resolution. VSI accurately identified coarse stocks (83.3% GSI concordance), and estimated a proxy for abundance (catch per unit effort, CPUE) of the upriver stock in the test fishery that was correlated (R2 = 0.90) with spring Chinook salmon abundance at Bonneville dam (Rkm 235). Salmon travel rates (~8.6 Rkm/day) provided predictions with 2‐week lead time prior to dam passage. Importantly, GSI resolved this predictive ability as finely as the hatchery broodstock level. Lower river stock CPUE in the test fishery was correlated with abundance at Willamette Falls (Rkm 196, R2 = 0.62), but could not be as finely resolved as achieved for upriver stocks. We described steps to combine VSI and GSI to provide timely in‐season information and with prediction accuracy of ~12.4 mean absolute percentage error and high stock resolution to help plan Columbia River mainstem fisheries. [ABSTRACT FROM AUTHOR]

Populations of many bumblebee species are declining, with distributions shifting northwards to track suitable climates. Climate change is considered a major contributing factor. Arctic species are particularly vulnerable as they cannot shift further north, making assessment of their population viability important. Analysis of levels of whole‐genome variation is a powerful way to analyse population declines and fragmentation. Here, we use genome sequencing to analyse genetic variation in seven species of bumblebee from the Scandinavian mountains, including two classified as vulnerable. We sequenced 333 samples from across the ranges of these species in Sweden. Estimates of effective population size (NE) vary from ~55,000 for species with restricted high alpine distributions to 220,000 for more widespread species. Population fragmentation is generally very low or undetectable over large distances in the mountains, suggesting an absence of barriers to gene flow. The relatively high NE and low population structure indicate that none of the species are at immediate risk of negative genetic effects caused by high levels of genetic drift. However, reconstruction of historical fluctuations in NE indicates that the arctic specialist species Bombus hyperboreus has experienced population declines since the last ice age and we detected one highly inbred diploid male of this species close to the southern limit of its range, potentially indicating elevated genetic load. Although the levels of genetic variation in montane bumblebee populations are currently relatively high, their ranges are predicted to shrink drastically due to the effects of climate change and monitoring is essential to detect future population declines. [ABSTRACT FROM AUTHOR]

Genome‐wide single nucleotide polymorphism (SNP) genotyping platforms have become increasingly popular in characterizing livestock and wildlife populations, replacing traditional methods such as microsatellite fragment analysis. Herein, we report the development and evaluation of a novel bison SNP panel for population management and conservation. Initially, 2474 autosomal SNPs were selected from existing bison whole‐genome sequences and variable sites among bison on the GGSP bovine 50K Chip, based on minor allele frequency, data completeness, and chromosome location. Additionally, 20 mitochondrial SNPs were chosen to identify known mitochondrial haplotypes in bison according to previous research. The SNPs were further evaluated using genotyping‐by‐sequencing with 190 bison, representing the historical lineages that survived the major population crash of the late 1800s. Variants with high potential for genotyping error were filtered out, and the remaining SNPs were placed on a custom Illumina™ array. The final panel consisting of 798 autosomal and 13 mitochondrial SNPs was used to establish baseline genetic parameters, compare populations, and assign mitochondrial haplotypes in 995 bison across ten populations. These SNPs were also found to be highly informative for individual animal identification and parentage assignment. This SNP panel provides a powerful new method to establish a baseline for estimating genetic health of bison populations and a new tool for bison managers to make informed management decisions based on genetic information specific to their populations. [ABSTRACT FROM AUTHOR]

Hybrid zones among mussel species have been extensively studied in the northern hemisphere. In South America, it has only recently become possible to study the natural hybrid zones, due to the clarification of the taxonomy of native mussels of the Mytilus genus. Analysing 54 SNP markers, we show the genetic species composition and admixture in the hybrid zone between M. chilensis and M. platensis in the southern end of South America. Bayesian, non‐Bayesian clustering and re‐assignment algorithms showed that the natural hybrid zone between M. chilensis and M. platensis in the Strait of Magellan, Isla Grande de Tierra del Fuego and the Falkland Islands shows clinal architecture. The hybrid zone can be divided into three different areas: the first one is on the Atlantic coast where only pure M. platensis and hybrid were found. In the second one, inside the Strait of Magellan, pure individuals of both species and mussels with variable degrees of hybridisation coexist. In the last area at the Strait in front of Punta Arenas City, fjords on the Isla Grande de Tierra del Fuego, and at the Beagle Channel, only M. chilensis and a low number of hybrids were found. According to the proportion of hybrids, bays with protected conditions away from strong currents would give better conditions for hybridisation. We do not find evidence of any other mussel species such as M. edulis, M. galloprovincialis, M. planulatus or M. trossulus in the zone. [ABSTRACT FROM AUTHOR]

Populations of Eurasian otters Lutra lutra, one of the most widely distributed apex predators in Eurasia, have been depleted mainly since the 1950s. However, a lack of information about their genomic diversity and how they are organized geographically in East Asia severely impedes our ability to monitor and conserve them in particular management units. Here, we re‐sequenced and analyzed 20 otter genomes spanning continental East Asia, including a population at Kinmen, a small island off the Fujian coast, China. The otters form three genetic clusters (one of L. l. lutra in the north and two of L. l. chinensis in the south), which have diverged in the Holocene. These three clusters should be recognized as three conservation management units to monitor and manage independently. The heterozygosity of the East Asian otters is as low as that of the threatened carnivores sequenced. Historical effective population size trajectories inferred from genomic variations suggest that their low genomic diversity could be partially attributed to changes in the climate since the mid‐Pleistocene and anthropogenic intervention since the Holocene. However, no evidence of genetic erosion, mutation load, or high level of inbreeding was detected in the presumably isolated Kinmen Island population. Any future in situ conservation efforts should consider this information for the conservation management units. [ABSTRACT FROM AUTHOR]

Multiple evolutionary processes influence genome‐wide allele frequencies and quantifying effects of genetic drift, and multiple forms of selection remain challenging in natural populations. Here, we investigate variation at major effect loci in contrast to patterns of neutral drift across a wide collection of steelhead (Oncorhynchus mykiss) populations that have declined in abundance due to anthropogenic impacts. Whole‐genome resequencing of 74 populations of steelhead revealed genome‐wide patterns (~8 million SNPs) consistent with expected neutral population structure. However, allelic variation at major effect loci associated with adult migration timing (chromosome 28: GREB1L/ROCK1) and age at maturity (chromosome 25: SIX6) reflected how selection has acted on phenotypic variation in contrast with neutral structure. Variation at major effect loci was influenced by evolutionary processes with differing signals between the strongly divergent Coastal and Inland lineages, while allele frequencies within and among populations within the Inland lineage have been driven by local natural selection as well as recent anthropogenic influences. Recent anthropogenic effects appeared to have influenced the frequency of major effect alleles including artificial selection for specific traits in hatchery stocks with subsequent gene flow into natural populations. Selection from environmental factors at various scales has also likely influenced variation for major effect alleles. These results reveal evolutionary mechanisms that influence allele frequencies at major effect loci that are critical for conservation of phenotypic traits and life history variation of this protected species. [ABSTRACT FROM AUTHOR]

Modern and ancient genomes are not necessarily drawn from homogeneous populations, as they may have been collected from different places and at different times. This heterogeneous sampling can be an issue for demographic inferences and results in biased demographic parameters and incorrect model choice if not properly considered. When explicitly accounted for, it can result in very complex models and high data dimensionality that are difficult to analyse. In this paper, we formally study the impact of such spatial and temporal sampling heterogeneity on demographic inference, and we introduce a way to circumvent this problem. To deal with structured samples without increasing the dimensionality of the site frequency spectrum (SFS), we introduce a new structured approach to the existing program fastsimcoal2. We assess the efficiency and relevance of this methodological update with simulated and modern human genomic data. We particularly focus on spatial and temporal heterogeneities to evidence the interest of this new SFS‐based approach, which can be especially useful when handling scattered and ancient DNA samples, as in conservation genetics or archaeogenetics. [ABSTRACT FROM AUTHOR]

Abstract Modern fisheries management strives to balance opposing goals of protection for weak stocks and opportunity for harvesting healthy stocks. Test fisheries can aid management of anadromous fishes if they can forecast the strength and timing of an annual run with adequate time to allow fisheries planning. Integration of genetic stock identification (GSI) can further maximize utility of test fisheries by resolving run forecasts into weak‐ and healthy‐stock subcomponents. Using 5 years (2017–2022) of test fishery data, our study evaluated accuracy, resolution, and lead time of predictions for stock‐specific run timing and abundance of Columbia River spring Chinook salmon (Oncorhynchus tshawytscha). We determined if this test fishery (1) could use visual stock identification (VSI) to forecast at the coarse stock resolution (i.e., classification of “lower” vs. “upriver” stocks) upon which current management is based and (2) could be enhanced with GSI to forecast at higher stock resolution. VSI accurately identified coarse stocks (83.3% GSI concordance), and estimated a proxy for abundance (catch per unit effort, CPUE) of the upriver stock in the test fishery that was correlated (R2 = 0.90) with spring Chinook salmon abundance at Bonneville dam (Rkm 235). Salmon travel rates (~8.6 Rkm/day) provided predictions with 2‐week lead time prior to dam passage. Importantly, GSI resolved this predictive ability as finely as the hatchery broodstock level. Lower river stock CPUE in the test fishery was correlated with abundance at Willamette Falls (Rkm 196, R2 = 0.62), but could not be as finely resolved as achieved for upriver stocks. We described steps to combine VSI and GSI to provide timely in‐season information and with prediction accuracy of ~12.4 mean absolute percentage error and high stock resolution to help plan Columbia River mainstem fisheries.

Abstract Genome‐wide single nucleotide polymorphism (SNP) genotyping platforms have become increasingly popular in characterizing livestock and wildlife populations, replacing traditional methods such as microsatellite fragment analysis. Herein, we report the development and evaluation of a novel bison SNP panel for population management and conservation. Initially, 2474 autosomal SNPs were selected from existing bison whole‐genome sequences and variable sites among bison on the GGSP bovine 50K Chip, based on minor allele frequency, data completeness, and chromosome location. Additionally, 20 mitochondrial SNPs were chosen to identify known mitochondrial haplotypes in bison according to previous research. The SNPs were further evaluated using genotyping‐by‐sequencing with 190 bison, representing the historical lineages that survived the major population crash of the late 1800s. Variants with high potential for genotyping error were filtered out, and the remaining SNPs were placed on a custom Illumina™ array. The final panel consisting of 798 autosomal and 13 mitochondrial SNPs was used to establish baseline genetic parameters, compare populations, and assign mitochondrial haplotypes in 995 bison across ten populations. These SNPs were also found to be highly informative for individual animal identification and parentage assignment. This SNP panel provides a powerful new method to establish a baseline for estimating genetic health of bison populations and a new tool for bison managers to make informed management decisions based on genetic information specific to their populations.

Abstract Multiple evolutionary processes influence genome‐wide allele frequencies and quantifying effects of genetic drift, and multiple forms of selection remain challenging in natural populations. Here, we investigate variation at major effect loci in contrast to patterns of neutral drift across a wide collection of steelhead (Oncorhynchus mykiss) populations that have declined in abundance due to anthropogenic impacts. Whole‐genome resequencing of 74 populations of steelhead revealed genome‐wide patterns (~8 million SNPs) consistent with expected neutral population structure. However, allelic variation at major effect loci associated with adult migration timing (chromosome 28: GREB1L/ROCK1) and age at maturity (chromosome 25: SIX6) reflected how selection has acted on phenotypic variation in contrast with neutral structure. Variation at major effect loci was influenced by evolutionary processes with differing signals between the strongly divergent Coastal and Inland lineages, while allele frequencies within and among populations within the Inland lineage have been driven by local natural selection as well as recent anthropogenic influences. Recent anthropogenic effects appeared to have influenced the frequency of major effect alleles including artificial selection for specific traits in hatchery stocks with subsequent gene flow into natural populations. Selection from environmental factors at various scales has also likely influenced variation for major effect alleles. These results reveal evolutionary mechanisms that influence allele frequencies at major effect loci that are critical for conservation of phenotypic traits and life history variation of this protected species.

Abstract Populations of Eurasian otters Lutra lutra, one of the most widely distributed apex predators in Eurasia, have been depleted mainly since the 1950s. However, a lack of information about their genomic diversity and how they are organized geographically in East Asia severely impedes our ability to monitor and conserve them in particular management units. Here, we re‐sequenced and analyzed 20 otter genomes spanning continental East Asia, including a population at Kinmen, a small island off the Fujian coast, China. The otters form three genetic clusters (one of L. l. lutra in the north and two of L. l. chinensis in the south), which have diverged in the Holocene. These three clusters should be recognized as three conservation management units to monitor and manage independently. The heterozygosity of the East Asian otters is as low as that of the threatened carnivores sequenced. Historical effective population size trajectories inferred from genomic variations suggest that their low genomic diversity could be partially attributed to changes in the climate since the mid‐Pleistocene and anthropogenic intervention since the Holocene. However, no evidence of genetic erosion, mutation load, or high level of inbreeding was detected in the presumably isolated Kinmen Island population. Any future in situ conservation efforts should consider this information for the conservation management units.

Understanding the rate and extent to which populations can adapt to novel environments at their ecological margins is fundamental to predicting the persistence of biological communities during ongoing and rapid global change. Recent range expansion in response to climate change in the UK butterfly Aricia agestis is associated with the evolution of novel interactions with a larval food plant, and the loss of its ability to use an ancestral host species. Using ddRAD analysis of 61,210 variable SNPs from 261 females from throughout the UK range of this species, we identify genomic regions at multiple chromosomes that are associated with evolutionary responses, and their association with demographic history and ecological variation. Gene flow appears widespread throughout the range, despite the apparently fragmented nature of the habitats used by this species. Patterns of haplotype variation between selected and neutral genomic regions suggest that evolution associated with climate adaptation is polygenic, resulting from the independent spread of alleles throughout the established range of this species, rather than the colonization of pre‐adapted genotypes from coastal populations. These data suggest that rapid responses to climate change do not depend on the availability of pre‐adapted genotypes. Instead, the evolution of novel forms of biotic interaction in A. agestis has occurred during range expansion, through the assembly of novel genotypes from alleles from multiple localities. [ABSTRACT FROM AUTHOR]

Abstract The dengue mosquito, Aedes aegypti (Linnaeus, 1762), is a highly invasive and medically significant vector of dengue, yellow fever, chikungunya and Zika viruses, whose global spread can be attributed to increased globalization in the 15th through 20th century. Records of the invasion history of Ae. aegypti across Southeast Asia are sparse and there is little knowledge regarding the invasion routes that the species exploited to gain a foothold in the Indo‐Pacific. Likewise, a broad and geographically thorough investigation of Ae. aegypti population genetics in the Indo‐Pacific is lacking, despite this region being highly impacted by diseases transmitted by this species. We assess 11 nuclear microsatellites and mitochondrial COI sequences, coupled with widespread sampling through the Indo‐Pacific region to characterise population structure at a broad geographic scale. We also perform a comprehensive literature search to collate documentation of the first known records of Ae. aegypti at various locations in the Indo‐Pacific. We revealed additional spatial population genetic structure of Ae. aegypti in Southeast Asia, the Indo‐Pacific and Australasia compared with previous studies and find differentiation between multiple Queensland and Torres Strait Islands populations. We also detected additional genetic breaks within Australia, Indonesia and Malaysia. Characterising the structure of previously unexplored populations through this region enhances the understanding of the population structure of Ae. aegypti in Australasia and Southeast Asia and may assist predictions of future mosquito movement, informing control strategies as well as assessing the risk of new invasion pathways.

Abstract The ability of local populations to adapt to future climate conditions is facilitated by a balance between short range dispersal allowing local buildup of adaptively beneficial alleles, and longer dispersal moving these alleles throughout the species range. Reef building corals have relatively low dispersal larvae, but most population genetic studies show differentiation only over 100s of km. Here, we report full mitochondrial genome sequences from 284 tabletop corals (Acropora hyacinthus) from 39 patch reefs in Palau, and show two signals of genetic structure across reef scales from 1 to 55 km. First, divergent mitochondrial DNA haplotypes exist in different proportions from reef to reef, causing PhiST values of 0.02 (p = 0.02). Second, closely related sequences of mitochondrial Haplogroups are more likely to be co‐located on the same reefs than expected by chance alone. We also compared these sequences to prior data on 155 colonies from American Samoa. In these comparisons, many Haplogroups in Palau were disproportionately represented or absent in American Samoa, and inter‐regional PhiST = 0.259. However, we saw three instances of identical mitochondrial genomes between locations. Together, these data sets suggest two features of coral dispersal revealed by occurrence patterns in highly similar mitochondrial genomes. First, the Palau‐American Samoa data suggest that long distance dispersal in corals is rare, as expected, but that it is common enough to deliver identical mitochondrial genomes across the Pacific. Second, higher than expected co‐occurrence of Haplogroups on the same Palau reefs suggests greater retention of coral larvae on local reefs than predicted by many current oceanographic models of larval movement. Increased attention to local scales of coral genetic structure, dispersal, and selection may help increase the accuracy of models of future adaptation of corals and of assisted migration as a reef resilience intervention.

Abstract Early blight of potato is caused by the fungal pathogen Alternaria solani and is an increasing problem worldwide. The primary strategy to control the disease is applying fungicides such as succinate dehydrogenase inhibitors (SDHI). SDHI‐resistant strains, showing reduced sensitivity to treatments, appeared in Germany in 2013, shortly after the introduction of SDHIs. Two primary mutations in the SDH complex (SdhB‐H278Y and SdhC‐H134R) have been frequently found throughout Europe. How these resistances arose and spread, and whether they are linked to other genomic features, remains unknown. For this project, we performed whole‐genome sequencing for 48 A. solani isolates from potato fields across Europe to better characterize the pathogen's genetic diversity in general and understand the development and spread of the genetic mutations that lead to SDHI resistance. The isolates can be grouped into seven genotypes. These genotypes do not show a geographical pattern but appear spread throughout Europe. We found clear evidence for recombination on the genome, and the observed admixtures might indicate a higher adaptive potential of the fungus than previously thought. Yet, we cannot link the observed recombination events to different Sdh mutations. The same Sdh mutations appear in different, non‐admixed genetic backgrounds; therefore, we conclude they arose independently. Our research gives insights into the genetic diversity of A. solani on a genome level. The mixed occurrence of different genotypes, apparent admixture in the populations, and evidence for recombination indicate higher genomic complexity than anticipated. The conclusion that SDHI tolerance arose multiple times independently has important implications for future fungicide resistance management strategies. These should not solely focus on preventing the spread of isolates between locations but also on limiting population size and the selective pressure posed by fungicides in a given field to avoid the rise of new mutations in other genetic backgrounds.

The genetic consequences of the subdivision of populations are regarded as significant to long‐term evolution, and research has shown that the scale and speed at which this is now occurring is critically reducing the adaptive potential of most species which inhabit human‐impacted landscapes. Here, we provide a rare and, to our knowledge, the first analysis of this process while it is happening and demonstrate a method of evaluating the effect of mitigation measures such as fauna crossings. We did this by using an extensive genetic data set collected from a koala population which was intensely monitored during the construction of linear transport infrastructure which resulted in the subdivision of their population. First, we found that both allelic richness and effective population size decreased through the process of population subdivision. Second, we predicted the extent to which genetic drift could impact genetic diversity over time and showed that after only 10 generations the resulting two subdivided populations could experience between 12% and 69% loss in genetic diversity. Lastly, using forward simulations we estimated that a minimum of eight koalas would need to disperse from each side of the subdivision per generation to maintain genetic connectivity close to zero but that 16 koalas would ensure that both genetic connectivity and diversity remained unchanged. These results have important consequences for the genetic management of species in human‐impacted landscapes by showing which genetic metrics are best to identify immediate loss in genetic diversity and how to evaluate the effectiveness of any mitigation measures. [ABSTRACT FROM AUTHOR]

Accelerating climate change and habitat loss make it imperative that plans to conserve biodiversity consider species' ability to adapt to changing environments. However, in biomes where biodiversity is highest, the evolutionary mechanisms responsible for generating adaptative variation and, ultimately, new species are frequently poorly understood. African rainforests represent one such biome, as decadal debates continue concerning the mechanisms generating African rainforest biodiversity. These debates hinge on the relative importance of geographic isolation versus divergent natural selection across environmental gradients. Hindering progress is a lack of robust tests of these competing hypotheses. Because African rainforests are severely at‐risk due to climate change and other anthropogenic activities, addressing this long‐standing debate is critical for making informed conservation decisions. We use demographic inference and allele frequency‐environment relationships to investigate mechanisms of diversification in an African rainforest skink, Trachylepis affinis, a species inhabiting the gradient between rainforest and rainforest‐savanna mosaic (ecotone). We provide compelling evidence of ecotone speciation, in which gene flow has all but ceased between rainforest and ecotone populations, at a level consistent with infrequent hybridization between sister species. Parallel patterns of genomic, morphological, and physiological divergence across this environmental gradient and pronounced allele frequency‐environment correlation indicate speciation is mostly probably driven by ecological divergence, supporting a central role for divergent natural selection. Our results provide strong evidence for the importance of ecological gradients in African rainforest speciation and inform conservation strategies that preserve the processes that produce and maintain biodiversity. [ABSTRACT FROM AUTHOR]

Current genetic methods of population assessment in conservation biology have been challenged by genome‐scale analyses due to their quantitatively novel insights. These analyses include assessments of runs‐of‐homozygosity (ROH), genomic evolutionary rate profiling (GERP), and mutational load. Here, we aim to elucidate the relationships between these measures using three divergent ungulates: white‐tailed deer, caribou, and mountain goat. The white‐tailed deer is currently expanding, while caribou are in the midst of a significant decline. Mountain goats remain stable, having suffered a large historical bottleneck. We assessed genome‐wide signatures of inbreeding using the inbreeding coefficient F and %ROH (FROH) and identified evolutionarily constrained regions with GERP. Mutational load was estimated by identifying mutations in highly constrained elements (CEs) and sorting intolerant from tolerant (SIFT) mutations. Our results showed that F and FROH are higher in mountain goats than in caribou and white‐tailed deer. Given the extended bottleneck and low Ne of the mountain goat, this supports the idea that the genome‐wide effects of demographic change take time to accrue. Similarly, we found that mountain goats possess more highly constrained CEs and the lowest dN/dS values, both of which are indicative of greater purifying selection; this is also reflected by fewer mutations in CEs and deleterious mutations identified by SIFT. In contrast, white‐tailed deer presented the highest mutational load with both metrics, in addition to dN/dS, while caribou were intermediate. Our results demonstrate that extended bottlenecks may lead to reduced diversity and increased FROH in ungulates, but not necessarily an increase in mutational load, probably due to the purging of deleterious alleles in small populations. [ABSTRACT FROM AUTHOR]

The dengue mosquito, Aedes aegypti (Linnaeus, 1762), is a highly invasive and medically significant vector of dengue, yellow fever, chikungunya and Zika viruses, whose global spread can be attributed to increased globalization in the 15th through 20th century. Records of the invasion history of Ae. aegypti across Southeast Asia are sparse and there is little knowledge regarding the invasion routes that the species exploited to gain a foothold in the Indo‐Pacific. Likewise, a broad and geographically thorough investigation of Ae. aegypti population genetics in the Indo‐Pacific is lacking, despite this region being highly impacted by diseases transmitted by this species. We assess 11 nuclear microsatellites and mitochondrial COI sequences, coupled with widespread sampling through the Indo‐Pacific region to characterise population structure at a broad geographic scale. We also perform a comprehensive literature search to collate documentation of the first known records of Ae. aegypti at various locations in the Indo‐Pacific. We revealed additional spatial population genetic structure of Ae. aegypti in Southeast Asia, the Indo‐Pacific and Australasia compared with previous studies and find differentiation between multiple Queensland and Torres Strait Islands populations. We also detected additional genetic breaks within Australia, Indonesia and Malaysia. Characterising the structure of previously unexplored populations through this region enhances the understanding of the population structure of Ae. aegypti in Australasia and Southeast Asia and may assist predictions of future mosquito movement, informing control strategies as well as assessing the risk of new invasion pathways. [ABSTRACT FROM AUTHOR]

Abstract Many coral reef fishes display remarkable genetic and phenotypic variation across their geographic ranges. Understanding how historical and contemporary processes have shaped these patterns remains a focal question in evolutionary biology since they reveal how diversity is generated and how it may respond to future environmental change. Here, we compare the population genomics and demographic histories of a commercially and ecologically important coral reef fish, the common coral grouper (Plectropomus leopardus [Lacépède 1802]), across two adjoining regions (the Great Barrier Reef; GBR, and the Coral Sea, Australia) spanning approximately 14 degrees of latitude and 9 degrees of longitude. We analysed 4548 single nucleotide polymorphism (SNP) markers across 11 sites and show that genetic connectivity between regions is low, despite their relative proximity (~100 km) and an absence of any obvious geographic barrier. Inferred demographic histories using 10,479 markers suggest that the Coral Sea population was founded by a small number of GBR individuals and that divergence occurred ~190 kya under a model of isolation with asymmetric migration. We detected population expansions in both regions, but estimates of contemporary effective population sizes were approximately 50% smaller in Coral Sea sites, which also had lower genetic diversity. Our results suggest that P. leopardus in the Coral Sea have experienced a long period of isolation that precedes the recent glacial period (~10–120 kya) and may be vulnerable to localized disturbances due to their relative reliance on local larval replenishment. While it is difficult to determine the underlying events that led to the divergence of the Coral Sea and GBR lineages, we show that even geographically proximate populations of a widely dispersed coral reef fish can have vastly different evolutionary histories.

Processes governing genetic diversity and adaptive potential in reef‐building corals are of interest both for fundamental evolutionary biology and for reef conservation. Here, we investigated the possibility of "sweepstakes reproductive success" (SRS) in a broadcast spawning coral, Acropora hyacinthus, at Yap Island, Micronesia. SRS is an extreme yearly variation in the number of surviving offspring among parents. It is predicted to generate genetically differentiated, low‐genetic‐diversity recruit cohorts, containing close kin individuals. We have tested these predictions by comparing genetic composition of size classes (adults and juveniles) at several sites on the island of Yap. We did see the genome‐wide dip in genetic diversity in juveniles compared to adults at two of the four sites; however, both adults and juveniles varied in genetic diversity across sites, and there was no detectable genetic structure among juveniles, which does not conform to the classical SRS scenario. Yet, we have identified a pair of juvenile siblings at the site where juveniles had the lowest genetic diversity compared to adults, an observation that is hard to explain without invoking SRS. While further support for SRS is needed to fully settle the issue, we show that incorporating SRS into the Indo‐West Pacific coral metapopulation adaptation model had surprisingly little effect on mean rates of coral cover decline during warming. Still, SRS notably increases year‐to‐year variation in coral cover throughout the simulation. [ABSTRACT FROM AUTHOR]

The ability of local populations to adapt to future climate conditions is facilitated by a balance between short range dispersal allowing local buildup of adaptively beneficial alleles, and longer dispersal moving these alleles throughout the species range. Reef building corals have relatively low dispersal larvae, but most population genetic studies show differentiation only over 100s of km. Here, we report full mitochondrial genome sequences from 284 tabletop corals (Acropora hyacinthus) from 39 patch reefs in Palau, and show two signals of genetic structure across reef scales from 1 to 55 km. First, divergent mitochondrial DNA haplotypes exist in different proportions from reef to reef, causing PhiST values of 0.02 (p = 0.02). Second, closely related sequences of mitochondrial Haplogroups are more likely to be co‐located on the same reefs than expected by chance alone. We also compared these sequences to prior data on 155 colonies from American Samoa. In these comparisons, many Haplogroups in Palau were disproportionately represented or absent in American Samoa, and inter‐regional PhiST = 0.259. However, we saw three instances of identical mitochondrial genomes between locations. Together, these data sets suggest two features of coral dispersal revealed by occurrence patterns in highly similar mitochondrial genomes. First, the Palau‐American Samoa data suggest that long distance dispersal in corals is rare, as expected, but that it is common enough to deliver identical mitochondrial genomes across the Pacific. Second, higher than expected co‐occurrence of Haplogroups on the same Palau reefs suggests greater retention of coral larvae on local reefs than predicted by many current oceanographic models of larval movement. Increased attention to local scales of coral genetic structure, dispersal, and selection may help increase the accuracy of models of future adaptation of corals and of assisted migration as a reef resilience intervention. [ABSTRACT FROM AUTHOR]

Widespread introduced species can be leveraged to investigate the genetic, ecological and adaptive processes underlying rapid evolution and range expansion, particularly the contributions of genetic diversity to adaptation. Rhinella marina, the cane toad, has been a focus of invasion biology for decades in Australia. However, their introduction history in North America is less clear. Here, we investigated the roles of introduction history and genetic diversity in establishment success of cane toads across their introduced range. We used reduced representation sequencing (ddRAD) to obtain 34,000 SNPs from 247 toads in native (French Guiana, Guyana, Ecuador, Panama, Texas) and introduced (Bermuda, southern Florida, northern Florida, Hawaiʻi, Puerto Rico) populations. Unlike all other cane toad introductions, we found that Florida populations were more closely related to native Central American lineages (R. horribilis), than to native Southern American lineages (R. marina). Furthermore, we found high levels of diversity and population structure in the native range, corroborating suggestions that R. marina is a species complex. We also found that introduced populations exhibit only slightly lower genetic diversity than native populations. Together with demographic analyses, this indicates founding populations of toads in Florida were larger than previously reported. Lastly, within R. marina, only one of 245 putatively adaptive SNPs showed fixed differences between native and introduced ranges, suggesting that putative selection in these introduced populations is based upon existing genetic variation. Our findings highlight the importance of genetic sequencing in understanding biological introductions and hint at the role of standing genetic variation in range expansion. [ABSTRACT FROM AUTHOR]

Early blight of potato is caused by the fungal pathogen Alternaria solani and is an increasing problem worldwide. The primary strategy to control the disease is applying fungicides such as succinate dehydrogenase inhibitors (SDHI). SDHI‐resistant strains, showing reduced sensitivity to treatments, appeared in Germany in 2013, shortly after the introduction of SDHIs. Two primary mutations in the SDH complex (SdhB‐H278Y and SdhC‐H134R) have been frequently found throughout Europe. How these resistances arose and spread, and whether they are linked to other genomic features, remains unknown. For this project, we performed whole‐genome sequencing for 48 A. solani isolates from potato fields across Europe to better characterize the pathogen's genetic diversity in general and understand the development and spread of the genetic mutations that lead to SDHI resistance. The isolates can be grouped into seven genotypes. These genotypes do not show a geographical pattern but appear spread throughout Europe. We found clear evidence for recombination on the genome, and the observed admixtures might indicate a higher adaptive potential of the fungus than previously thought. Yet, we cannot link the observed recombination events to different Sdh mutations. The same Sdh mutations appear in different, non‐admixed genetic backgrounds; therefore, we conclude they arose independently. Our research gives insights into the genetic diversity of A. solani on a genome level. The mixed occurrence of different genotypes, apparent admixture in the populations, and evidence for recombination indicate higher genomic complexity than anticipated. The conclusion that SDHI tolerance arose multiple times independently has important implications for future fungicide resistance management strategies. These should not solely focus on preventing the spread of isolates between locations but also on limiting population size and the selective pressure posed by fungicides in a given field to avoid the rise of new mutations in other genetic backgrounds. [ABSTRACT FROM AUTHOR]

Many coral reef fishes display remarkable genetic and phenotypic variation across their geographic ranges. Understanding how historical and contemporary processes have shaped these patterns remains a focal question in evolutionary biology since they reveal how diversity is generated and how it may respond to future environmental change. Here, we compare the population genomics and demographic histories of a commercially and ecologically important coral reef fish, the common coral grouper (Plectropomus leopardus [Lacépède 1802]), across two adjoining regions (the Great Barrier Reef; GBR, and the Coral Sea, Australia) spanning approximately 14 degrees of latitude and 9 degrees of longitude. We analysed 4548 single nucleotide polymorphism (SNP) markers across 11 sites and show that genetic connectivity between regions is low, despite their relative proximity (~100 km) and an absence of any obvious geographic barrier. Inferred demographic histories using 10,479 markers suggest that the Coral Sea population was founded by a small number of GBR individuals and that divergence occurred ~190 kya under a model of isolation with asymmetric migration. We detected population expansions in both regions, but estimates of contemporary effective population sizes were approximately 50% smaller in Coral Sea sites, which also had lower genetic diversity. Our results suggest that P. leopardus in the Coral Sea have experienced a long period of isolation that precedes the recent glacial period (~10–120 kya) and may be vulnerable to localized disturbances due to their relative reliance on local larval replenishment. While it is difficult to determine the underlying events that led to the divergence of the Coral Sea and GBR lineages, we show that even geographically proximate populations of a widely dispersed coral reef fish can have vastly different evolutionary histories. [ABSTRACT FROM AUTHOR]

Polymorphism for immune functions can explain significant variation in health and reproductive success within species. Drastic loss in genetic diversity at such loci constitutes an extinction risk and should be monitored in species of conservation concern. However, effective implementations of genome‐wide immune polymorphism sets into high‐throughput genotyping assays are scarce. Here, we report the design and validation of a microfluidics‐based amplicon sequencing assay to comprehensively capture genetic variation in Alpine ibex (Capra ibex). This species represents one of the most successful large mammal restorations recovering from a severely depressed census size and a massive loss in diversity at the major histocompatibility complex (MHC). We analysed 65 whole‐genome sequencing sets of the Alpine ibex and related species to select the most representative markers and to prevent primer binding failures. In total, we designed ~1,000 amplicons densely covering the MHC, further immunity‐related genes as well as randomly selected genome‐wide markers for the assessment of neutral population structure. Our analysis of 158 individuals shows that the genome‐wide markers perform equally well at resolving population structure as RAD‐sequencing or low‐coverage genome sequencing data sets. Immunity‐related loci show unexpectedly high degrees of genetic differentiation within the species. Such information can now be used to define highly targeted individual translocations. Our design strategy can be realistically implemented into genetic surveys of a large range of species. In conclusion, leveraging whole‐genome sequencing data sets to design targeted amplicon assays allows the simultaneous monitoring of multiple genetic risk factors and can be translated into species conservation recommendations. [ABSTRACT FROM AUTHOR]

Understanding the consequences of local adaptation at the genomic diversity is a central goal in evolutionary genetics of natural populations. In species with large continuous geographical distributions the phenotypic signal of local adaptation is frequently clear, but the genetic basis often remains elusive. We examined the patterns of genetic diversity in Pinus sylvestris, a keystone species in many Eurasian ecosystems with a huge distribution range and decades of forestry research showing that it is locally adapted to the vast range of environmental conditions. Making P. sylvestris an even more attractive subject of local adaptation study, population structure has been shown to be weak previously and in this study. However, little is known about the molecular genetic basis of adaptation, as the massive size of gymnosperm genomes has prevented large scale genomic surveys. We generated a both geographically and genomically extensive dataset using a targeted sequencing approach. By applying divergence-based and landscape genomics methods we identified several loci contributing to local adaptation, but only few with large allele frequency changes across latitude. We also discovered a very large (ca. 300 Mbp) putative inversion potentially under selection, which to our knowledge is the first such discovery in conifers. Our results call for more detailed analysis of structural variation in relation to genomic basis of local adaptation, emphasize the lack of large effect loci contributing to local adaptation in the coding regions and thus point out the need for more attention toward multi-locus analysis of polygenic adaptation.

Abstract Poleward species range shifts have been predicted to result from climate change, and many observations have confirmed such movement. Poleward shifts may represent a homogeneous shift in distribution, seasonal northward movement of specific populations, or colonization processes at the poleward edge of the distribution. The ecosystem of the Bering Sea has been changing along with the climate, moving from an arctic to a subarctic system. Several fish species have been observed farther north than previously reported and in increasing abundances. We examined one of these fish species, Pacific cod, in the northern Bering Sea (NBS) to assess whether they migrated from another stock in the eastern Bering Sea (EBS), Gulf of Alaska, or Aleutian Islands, or whether they represent a separate population. Genetic analyses using 3,599 single nucleotide polymorphism markers indicated that nonspawning cod collected in August 2017 in the NBS were similar to spawning stocks of cod in the EBS. This result suggests escalating northward movement of the large EBS stock during summer months. Whether the cod observed in the NBS migrate south during winter to spawn or remain in the NBS as a sink population is unknown.

The origins of geographic races in wide‐ranging species are poorly understood. In Texas, the texanus subspecies of Helianthusannuus has long been thought to have acquired its defining phenotypic traits via introgression from a local congener, H. debilis, but previous tests of this hypothesis were inconclusive. Here, we explore the origins of H. a. texanus using whole genome sequencing data from across the entire range of H. annuus and possible donor species, as well as phenotypic data from a common garden study. We found that although it is morphologically convergent with H. debilis, H. a. texanus has conflicting signals of introgression. Genome wide tests (Patterson's D and TreeMix) only found evidence of introgression from H. argophyllus (sister species to H. annuus and also sympatric), but not H. debilis, with the exception of one individual of 109 analysed. We further scanned the genome for localized signals of introgression using PCAdmix and found minimal but nonzero introgression from H. debilis and significant introgression from H. argophyllus in some populations. Given the paucity of introgression from H. debilis, we argue that the morphological convergence observed in Texas is probably from standing genetic variation. We also found that genomic differentiation in H. a. texanus is mostly driven by large segregating inversions, several of which have signatures of natural selection based on haplotype frequencies. [ABSTRACT FROM AUTHOR]

Various hypotheses have been proposed regarding the origin of the plant pathogen Phytophthora cinnamomi. P. cinnamomi is a devastating, highly invasive soilborne pathogen associated with epidemics of agricultural, horticultural and forest plantations and native ecosystems worldwide. We conducted a phylogeographic analysis of populations of this pathogen sampled in Asia, Australia, Europe, southern and northern Africa, South America, and North America. Based on genotyping‐by‐sequencing, we observed the highest genotypic diversity in Taiwan and Vietnam, followed by Australia and South Africa. Mating type ratios were in equal proportions in Asia as expected for a sexual population. Simulations based on the index of association suggest a partially sexual, semi‐clonal mode of reproduction for the Taiwanese and Vietnamese populations while populations outside of Asia are clonal. Ancestral area reconstruction provides new evidence supporting Taiwan as the ancestral area, given our sample, indicating that this region might be near or at the centre of origin for this pathogen as speculated previously. The Australian and South African populations appear to be a secondary centre of diversity following migration from Taiwan or Vietnam. Our work also identified two panglobal, clonal lineages PcG1‐A2 and PcG2‐A2 of A2 mating type found on all continents. Further surveys of natural forests across Southeast Asia are needed to definitively locate the actual centre of origin of this important plant pathogen. [ABSTRACT FROM AUTHOR]

Since allozymes were first used to assess genetic diversity in the 1960s and 1970s, biologists have attempted to characterize gene pools and conserve the diversity observed in domestic crops, livestock, zoos and (more recently) natural populations. Recently, some authors have claimed that the importance of genetic diversity in conservation biology has been greatly overstated. Here, we argue that a voluminous literature indicates otherwise. We address four main points made by detractors of genetic diversity's role in conservation by using published literature to firmly establish that genetic diversity is intimately tied to evolutionary fitness, and that the associated demographic consequences are of paramount importance to many conservation efforts. We think that responsible management in the Anthropocene should, whenever possible, include the conservation of ecosystems, communities, populations and individuals, and their underlying genetic diversity. [ABSTRACT FROM AUTHOR]

After centuries of human hunting, the Eurasian beaver Castor fiber had disappeared from most of its original range by the end of the 19th century. The surviving relict populations are characterized by both low genetic diversity and strong phylogeographical structure. However, it remains unclear whether these attributes are the result of a human-induced, late Holocene bottleneck or already existed prior to this reduction in range. To investigate genetic diversity in Eurasian beaver populations during the Holocene, we obtained mitochondrial control region DNA sequences from 48 ancient beaver samples and added 152 modern sequences from GenBank. Phylogeographical analyses of the data indicate a differentiation of European beaver populations into three mitochondrial clades. The two main clades occur in western and eastern Europe, respectively, with an early Holocene contact zone in eastern Europe near a present-day contact zone. A divergent and previously unknown clade of beavers from the Danube Basin survived until at least 6000 years ago, but went extinct during the transition to modern times. Finally, we identify a recent decline in effective population size of Eurasian beavers, with a stronger bottleneck signal in the western than in the eastern clade. Our results suggest that the low genetic diversity and the strong phylogeographical structure in recent beavers are artefacts of human hunting-associated population reductions. While beaver populations have been growing rapidly since the late 19th century, genetic diversity within modern beaver populations remains considerably reduced compared to what was present prior to the period of human hunting and habitat reduction.

The M and S forms of Anopheles gambiae have been the subject of intense study, but are morphologically indistinguishable and can only be identified using molecular techniques. PCR-based assays to distinguish the two forms have been designed and applied widely. However, the application of these assays towards identifying hybrids between the two forms, and backcrossed hybrids in particular, has been problematic as the currently available diagnostic assays are based on single locus and/or are located within a multicopy gene. Here, we present an alternative genotyping method for detecting hybridization and introgression between M and S molecular forms based on a multilocus panel of single-nucleotide polymorphisms (SNPs) fixed between the M and S forms. The panel of SNPs employed is located in so-called islands of divergence leading us to describe this method as the 'Divergence Island SNP' (DIS) assay. We show this multilocus SNP genotyping approach can robustly and accurately detect F1 hybrids as well as backcrossed individuals.

Population genetic analyses often use summary statistics to describe patterns of genetic variation and provide insight into evolutionary processes. Among the most fundamental of these summary statistics are π and dXY, which are used to describe genetic diversity within and between populations, respectively. Here, we address a widespread issue in π and dXY calculation: systematic bias generated by missing data of various types. Many popular methods for calculating π and dXY operate on data encoded in the variant call format (VCF), which condenses genetic data by omitting invariant sites. When calculating π and dXY using a VCF, it is often implicitly assumed that missing genotypes (including those at sites not represented in the VCF) are homozygous for the reference allele. Here, we show how this assumption can result in substantial downward bias in estimates of π and dXY that is directly proportional to the amount of missing data. We discuss the pervasive nature and importance of this problem in population genetics, and introduce a user‐friendly UNIX command line utility, pixy, that solves this problem via an algorithm that generates unbiased estimates of π and dXY in the face of missing data. We compare pixy to existing methods using both simulated and empirical data, and show that pixy alone produces unbiased estimates of π and dXY regardless of the form or amount of missing data. In summary, our software solves a long‐standing problem in applied population genetics and highlights the importance of properly accounting for missing data in population genetic analyses. [ABSTRACT FROM AUTHOR]

Understanding the consequences of local adaptation at the genomic diversity is a central goal in evolutionary genetics of natural populations. In species with large continuous geographical distributions the phenotypic signal of local adaptation is frequently clear, but the genetic basis often remains elusive. We examined the patterns of genetic diversity in Pinus sylvestris, a keystone species in many Eurasian ecosystems with a huge distribution range and decades of forestry research showing that it is locally adapted to the vast range of environmental conditions. Making P. sylvestris an even more attractive subject of local adaptation study, population structure has been shown to be weak previously and in this study. However, little is known about the molecular genetic basis of adaptation, as the massive size of gymnosperm genomes has prevented large scale genomic surveys. We generated a both geographically and genomically extensive dataset using a targeted sequencing approach. By applying divergence-based and landscape genomics methods we identified several loci contributing to local adaptation, but only few with large allele frequency changes across latitude. We also discovered a very large (ca. 300 Mbp) putative inversion potentially under selection, which to our knowledge is the first such discovery in conifers. Our results call for more detailed analysis of structural variation in relation to genomic basis of local adaptation, emphasize the lack of large effect loci contributing to local adaptation in the coding regions and thus point out the need for more attention toward multi-locus analysis of polygenic adaptation. [ABSTRACT FROM AUTHOR]

It is hard to imagine a world without food‐associated microbes. The production of bread, wine, beer, salami, coffee, chocolate, cheese and many other foods and beverages all rely on specific microbes. In cheese, myriad microbial species collaborate to yield the complex organoleptic properties that are appreciated by millions of people worldwide. In the early days of cheese making, these complex communities emerged spontaneously from the natural flora associated with the raw materials, the equipment, the production environment or craftsmen involved in the production process. However, in some cases, the microbes shifted their natural habitat to the new cheese‐associated environment. The most obvious cause of this is backslopping, where part of a fermented product is used to inoculate the next batch. In addition, some microbes may simply adhere to the tools used in the production process. These microbial communities gradually adapted to the novel man‐made niches, a process referred to as "domestication." Domestication is associated with specific genomic and phenotypic changes and ultimately leads to lineages that are genetically and phenotypically distinct from their wild ancestors. In this issue of Molecular Ecology, Dumas et al. have investigated a prime example of cheese‐associated microbes, the fungus Penicillium roqueforti. The authors identified several hallmarks of domestication in the genome and phenome of this species, allowing them to hypothesize about the origin of blue‐veined cheese fungi domestication, and the specific evolutionary processes involved in adaptation to the cheese matrix. [ABSTRACT FROM AUTHOR]

Since the early Holocene, fish population genetics in the Laurentian Great Lakes have been shaped by the dual influences of habitat structure and post‐glacial dispersal. Riverscape genetics theory predicts that longitudinal habitat corridors and unidirectional downstream water‐flow drive the downstream accumulation of genetic diversity, whereas post‐glacial dispersal theory predicts that fish genetic diversity should decrease with increasing distance from glacial refugia. This study examines populations of seven native fish species codistributed above and below the 58 m high Niagara Falls – a hypothesized barrier to gene flow in aquatic species. A better understanding of Niagara Falls' role as a barrier to gene flow and dispersal is needed to identify drivers of Great Lakes genetic diversity and guide strategies to limit exotic species invasions. We used genome‐wide SNPs and coalescent models to test whether populations are: (a) genetically distinct, consistent with the Niagara Falls barrier hypothesis; (b) more genetically diverse upstream, consistent with post‐glacial expansion theory, or downstream, consistent with the riverscape habitat theory; and (c) have migrated either upstream or downstream past Niagara Falls. We found that genetic diversity is consistently greater below Niagara Falls and the falls are an effective barrier to migration, but two species have probably dispersed upstream past the falls after glacial retreat yet before opening of the Welland Canal. Models restricting migration to after opening of the Welland Canal were generally rejected. These results help explain how river habitat features affect aquatic species' genetic diversity and highlight the need to better understand post‐glacial dispersal pathways. [ABSTRACT FROM AUTHOR]