Many populations of freshwater fishes are threatened with losses, and increasingly, the release of hatchery individuals is one strategy being implemented to support wild populations. However, stocking of hatchery individuals may pose long‐term threats to wild populations, particularly if genetic interactions occur between wild and hatchery individuals. One highly prized sport fish that has been heavily stocked throughout its range is the brook trout (Salvelinus fontinalis). In Nova Scotia, Canada, hatchery brook trout have been stocked since the early 1900s, and despite continued stocking efforts, populations have suffered declines in recent decades. Before this study, the genetic structure of brook trout populations in the province was unknown; however, given the potential negative consequences associated with hatchery stocking, it is possible that hatchery programs have adversely affected the genetic integrity of wild populations. To assess the influence of hatchery supplementation on wild populations, we genotyped wild brook trout from 12 river systems and hatchery brook trout from two major hatcheries using 100 microsatellite loci. Genetic analyses of wild trout revealed extensive population genetic structure among and within river systems and significant isolation‐by‐distance. Hatchery stocks were genetically distinct from wild populations, and most populations showed limited to no evidence of hatchery introgression (<5% hatchery ancestry). Only a single location had a substantial number of hatchery‐derived trout and was located in the only river where a local strain is used for supplementation. The amount of hatchery stocking within a watershed did not influence the level of hatchery introgression. Neutral genetic structure of wild populations was influenced by geography with some influence of climate and stocking indices. Overall, our study suggests that long‐term stocking has not significantly affected the genetic integrity of wild trout populations, highlighting the variable outcomes of stocking and the need to evaluate the consequences on a case‐by‐case basis.

The variance in phenotypic trait values is a product of environmental and genetic variation. The sensitivity of traits to environmental variation has a genetic component and is likely to be under selection. However, there are few studies investigating the evolution of this sensitivity, in part due to the challenges of estimating the environmental variance. The livestock literature provides a wealth of studies that accurately partition components of phenotypic variance, including the environmental variance, in well‐defined environments. These studies involve breeds that have been under strong selection on mean phenotype in optimal environments for many generations, and therefore represent an opportunity to study the potential evolution of trait sensitivity to environmental conditions. Here we use literature on domestic cattle to examine the evolution of micro‐environmental variance (CVR – the coefficient of residual variance) by testing for differences in expression of CVR in animals from the same breed reared in different environments. Traits that have been under strong selection did not follow a null‐expectation of an increase in CVR in heterogenous environments (e.g. grazing), a pattern that may reflect evolution of increased uniformity in heterogeneous environments. When comparing CVR across environments of different levels of optimality, here measured by trait mean, we found a reduction in CVR in the more optimal environments for both life‐history and growth traits. Selection aimed at increasing trait means in livestock breeds typically occurs in the more optimal environments, and we therefore suspect that the decreased CVR is a consequence of evolution of the expression of micro‐environmental variance in this environment. Our results highlight the heterogeneity in micro‐environmental variance across environments and point to possible connections to the intensity of selection on trait means.

The resiliency of populations and species to environmental change is dependent on the maintenance of genetic diversity, and as such quantifying diversity is central to combatting ongoing wide spread reductions in biodiversity. With the advent of next‐generation sequencing, several methods now exist for resolving fine‐scale population structure, but the comparative performance of these methods for genetic assignment has rarely been tested. Here we evaluate the performance of sequenced microsatellites and a single nucleotide polymorphism (SNP) array to resolve fine‐scale population structure in a critically important salmonid in northeastern Canada, Arctic charr (Salvelinus alpinus). We also assess the utility of sequenced microsatellites for fisheries applications by quantifying the spatial scales of movement and exploitation through genetic assignment of fishery samples to rivers of origin and comparing these results with a 29‐year tagging dataset. Self‐assignment and simulation‐based analyses of 111 genome‐wide microsatellite loci and 500 informative SNPs from 28 populations of Arctic charr in northeastern Canada identified largely river‐specific genetic structure. Despite large differences (~4X) in the number of loci surveyed between panels, mean self‐assignment accuracy was similar with the SNP panel and with the microsatellite loci (>90%). Subsequent analysis of 996 fishery‐collected samples using the microsatellite panel revealed that larger rivers contribute greater numbers of individuals to the fishery, and that coastal fisheries largely exploit individuals originating from nearby rivers, corroborating results from traditional tagging experiments. Our results demonstrate the efficacy of sequence‐based microsatellite genotyping to advance understanding of fine‐scale population structure and harvest composition in northern and understudied species.

Bacillus anthracis, the causative agent of anthrax, is a considerable global health threat affecting wildlife, livestock, and the general public. In this study, whole‐genome sequence analysis of over 350 B. anthracis isolates was used to establish a new high‐resolution global genotyping framework that is both biogeographically informative and compatible with multiple genomic assays. The data presented in this study shed new light on the diverse global dissemination of this species and indicate that many lineages may be uniquely suited to the geographic regions in which they are found. In addition, we demonstrate that plasmid genomic structure for this species is largely consistent with chromosomal population structure, suggesting vertical inheritance in this bacterium has contributed to its evolutionary persistence. This classification methodology is the first based on population genomic structure for this species and has potential use for local and broader institutions seeking to understand both disease outbreak origins and recent introductions. In addition, we provide access to a newly developed genotyping script as well as the full whole‐genome sequence analyses output for this study, allowing future studies to rapidly employ and append their data in the context of this global collection. This framework may act as a powerful tool for public health agencies, wildlife disease laboratories, and researchers seeking to utilize and expand this classification scheme for further investigations into B. anthracis evolution.

The objective of this study was to test, using a field experiment, the effect of genotypic diversity on productivity of farmed populations (Ancud and Chaica, Chile) of the domesticated red alga Agarophyton chilense (formerly known as Gracilaria chilensis), a species considered as economically important in Chile. Monoclonal and polyclonal (4 and 8 genotypes) subplots were outplanted into the mid intertidal in Metri Bay (Puerto Montt, Chile) during summer, a season in which A. chilense face higher temperatures (>18°C) and low nitrogen availability (<4.00 μmol). Ancud farm genotypes show higher growth rates in the monoclonal rather than the two polyclonal subplots. A similar tendency, yet not significant, was discernible in Chaica. In addition, whatever the population of origin of the thalli, no effect of genotypic diversity was detected neither on the agar yield and its quality, nor on the epiphyte load. Such unexpected results of a higher performance in plots with a lower genotypic diversity could be explained (a) by human‐assisted selection for dominant‐best‐performing genotypes that could counterbalance the negative effect caused by the low genotypic diversity in farms and (b) by the fact that the organisms inhabiting the algal mats do not impact the fitness of their host. Overall, the results obtained here suggest that despite farm induced selection lead to impoverished pools of genotypes, they may also have a positive effect of on the resistance of farmed populations to seasonal stressors. However, whether this may have a secondary negative effect on the longer term in a fluctuating environment remains to be determined, but may be avoided by adopting strategy of selection favoring different genotypes in space and time, as implemented in forestry.

Strong and ongoing artificial selection in domestic animals has resulted in amazing phenotypic responses that benefit humans, but often at a cost to an animal’s health, and problems related to inbreeding depression, including a higher incidence of cancer. Despite high rates of cancer in domesticated species, little attention has been devoted to exploring the hypothesis that persistent artificial selection may also favour the evolution of compensatory anticancer defences. Indeed, there is evidence for effective anti‐cancer defences found in several domesticated species associated with different cancer types. We also suggest that artificial selection can favour the “domestication” of inherited oncogenic mutations in rare instances, retaining those associated to late and/or less aggressive cancers, and that by studying these seemingly rare anti‐cancer adaptations novel cancer treatments may be found.

Effective resource management depends on our ability to partition diversity into biologically meaningful units. Recent evolutionary divergence, however, can often lead to ambiguity in morphological and genetic differentiation, complicating the delineation of valid conservation units. Such is the case with the “coregonine problem,” where recent post‐glacial radiations of coregonines into lacustrine habitats resulted in the evolution of numerous species flocks, often with ambiguous taxonomy. The application of genomics methods is beginning to shed light on this problem and the evolutionary mechanisms underlying divergence in these ecologically and economically important fishes. Here, we used restriction site‐associated DNA (RAD) sequencing to examine genetic diversity and differentiation among sympatric forms in the Coregonus artedi complex in the Apostle Islands of Lake Superior, the largest lake in the Laurentian Great Lakes. Using 29,068 SNPs, we were able to clearly distinguish among the three most common forms for the first time, as well as identify putative hybrids and potentially mis‐identified specimens. Population assignment rates for these forms using our RAD data were 93‐100% with the only mis‐assignments arising from putative hybrids, an improvement from 62‐77% using microsatellites. Estimates of pairwise differentiation (FST: 0.045‐0.056) were large given the detection of hybrids, suggesting that reduced fitness of hybrid individuals may be a potential mechanism for the maintenance of differentiation. We also used a newly built C. artedi linkage map to look for islands of genetic divergence among forms and found widespread differentiation across the genome, a pattern indicative of long‐term drift, suggesting that these forms have been reproductively isolated for a substantial amount of time. The results of this study provide valuable information that can be applied to develop well‐informed management strategies and stress the importance of re‐evaluating conservation units with genomic tools to ensure they accurately reflect species diversity.

Genetic diversity is determinant for pest species' success and vector competence. Understanding the ecological and evolutionary processes that determine the genetic diversity is fundamental to help identify the spatial scale at which pest populations are best managed. In the present study, we present the first comprehensive analysis of the genetic diversity and evolution of Rhopalosiphum padi, a major pest of cereals and a main vector of the barley yellow dwarf virus (BYDV), in England. We have used a genotyping by sequencing approach to study whether i) there is any underlying population genetic structure at a national and regional scale in this pest that can disperse long distances; ii) the populations evolve as a response to environmental change and selective pressures, and; iii) the populations comprise anholocyclic lineages. Individual R. padi were collected using the Rothamsted Insect Survey's suction‐trap network at several sites across England between 2004 and 2016 as part of the RIS long‐term nationwide surveillance. Results identified two genetic clusters in England that mostly corresponded to a North – South division, although gene flow is ongoing between the two subpopulations. These genetic clusters do not correspond to different life cycles types, and cyclical parthenogenesis is predominant in England. Results also show that there is dispersal with gene flow across England, although there is a reduction between the northern and southern sites with the Southwestern population being the most genetically differentiated. There is no evidence for isolation‐by‐distance and other factors like primary host distribution, uncommon in the south and absent in the southwest, could influence the dispersal patterns. Finally, results also show no evidence for the evolution of the R. padi population, and it is demographically stable despite the ongoing environmental change. These results are discussed in view of their relevance to pest management and the transmission of BYDV.

Habitat fragmentation is a major extinction driver. Despite dramatically increasing fragmentation across the globe, its specific impacts on population connectivity across species with differing life histories remain difficult to characterize, let alone quantify. Here we investigate patterns of population connectivity in six songbird species from Singapore, a highly fragmented tropical rainforest island. Using massive panels of genome‐wide single‐nucleotide polymorphisms across dozens of samples per species, we examined population genetic diversity, inbreeding, gene flow and connectivity among species along a spectrum of ecological specificities. We found a higher resilience to habitat fragmentation in edge‐tolerant and forest‐canopy species as compared to forest‐dependent understorey insectivores. The latter exhibited levels of genetic diversity up to three times lower in Singapore than in populations from contiguous forest elsewhere. Using dense genomic and geographic sampling, we identified individual barriers such as reservoirs that effectively minimize gene flow in sensitive understorey birds, revealing that terrestrial forest species may exhibit levels of sensitivity to fragmentation far greater than previously expected. This study provides a blueprint for conservation genomics at small scales with a view to identifying preferred locations for habitat corridors, flagging candidate populations for re‐stocking with translocated individuals and improving the design of future reserves.

The giant panda is one of the most endangered mammals in the world, and many studies have revealed their evolutionary adaptation to the local environment (e.g., dietary cellulose and cyanide) on the evidences from population genetics and their gut microbiome. Here, based on the results of our analysis of the giant panda gut microbiome, we concluded that instability and resilience are the two primary characteristics of the giant panda gut microbiome. This basic information may have an impact on giant panda conservation, as well the management of other animal species.

In the open ocean without terrain boundaries, marine invertebrates with pelagic larvae can migrate long distances using ocean currents, suggesting reduced genetic diversification. Contrary to this assumption, however, genetic differentiation is often observed in marine invertebrates. In the present study, we sought to explain how population structure is established in the western Pacific Ocean, where the strong Kuroshio Current maintains high levels of gene flow from south to north, presumably promoting genetic homogeneity. We determined the population structure of the pearl oyster, Pinctada fucata, in the Indo‐Pacific Ocean using genome‐wide genotyping data from multiple sampling localities. Cluster analysis showed that the western Pacific population is distinct from that of the Indian Ocean, and that it is divided into northern (Japanese mainland) and southern (Nansei Islands, China, and Cambodia) populations. Genetic differentiation of P. fucata can be explained by geographic barriers in the Indian Ocean and a local lagoon, and by environmental gradients of sea surface temperature (SST) and oxygen concentration in the western Pacific. A genome scan showed evidence of adaptive evolution in genomic loci, possibly associated with changes in environmental factors, including SST and oxygen concentration. Furthermore, Bayesian simulation demonstrated that the past population expansion and division are congruent with ocean warming after the last glacial period. It is highly likely that the environmental gradient forms a genetic barrier that diversifies P. fucata populations in the western Pacific. This hypothesis helps to explain genetic differentiation and possible speciation of marine invertebrates.

Conservation management strategies for many highly threatened species include conservation breeding to prevent extinction and enhance recovery. Pairing decisions for these conservation breeding programmes can be informed by pedigree data to minimise relatedness between individuals in an effort to avoid inbreeding, maximise diversity, and maintain evolutionary potential. However, conservation breeding programmes struggle to use this approach when pedigrees are shallow or incomplete. While genetic data (i.e., microsatellites) can be used to estimate relatedness to inform pairing decisions, emerging evidence indicates this approach may lack precision in genetically depauperate species, and more effective estimates will likely be obtained from genomic data (i.e., thousands of genome‐wide single nucleotide polymorphisms, or SNPs). Here, we compare relatedness estimates and subsequent pairing decisions using pedigrees, microsatellites, and SNPs from whole genome resequencing approaches in two critically endangered birds endemic to New Zealand: kakī/black stilt (Himantopus novaezelandiae) and kākāriki karaka/orange‐fronted parakeet (Cyanoramphus malherbi). Our findings indicate that SNPs provide more precise estimates of relatedness than microsatellites when assessing empirical parent‐offspring and full sibling relationships. Further, our results show that relatedness estimates and subsequent pairing recommendations using PMx are most similar between pedigree‐ and SNP‐based approaches. These combined results indicate that in lieu of robust pedigrees, SNPs are an effective tool for informing pairing decisions, which has important implications for many poorly pedigreed conservation breeding programmes worldwide.

Invasive species have led to precipitous declines in biodiversity, especially in island systems. Brown (Rattus norvegicus) and black rats (R. rattus) are among the most invasive animals on the planet, with eradication being the primary tool for established island populations. The need for increased research for defining eradication units and monitoring outcomes has been highlighted as a means to maximize success. Haida Gwaii is an archipelago ~100 km off the northern coast of British Columbia, Canada, that hosts globally significant breeding populations of seabirds that are at risk due to invasive rats. Here, we paired sampling of brown (n = 287) and black (n = 291) rats across the Haida Gwaii archipelago with genotyping by sequencing (10,770–27,686 SNPs) to investigate patterns of population connectivity and infer levels/direction of gene flow among invasive rat populations in Haida Gwaii. We reconstructed three regional clusters for both species (north, central and south), with proximate populations within regions being largely more related than those that were more distant, consistent with predictions from island biogeography theory. Population assignment of recently detected individuals post‐eradication on Faraday, Murchison and the Bischof Islands revealed all were re‐invaders from Lyell Island, rather than being on‐island survivors. Based on these results, we identified six eradication units constituting single or clusters of islands that would limit the potential for reinvasion, some of which will need to be combined with biosecurity measures. Overall, our results highlight the importance of targeted research prior to conducting eradications and demonstrate a framework for applying population genomics for guiding invasive species management in island systems.

Populations may respond to environmental heterogeneity via evolutionary divergence or phenotypic plasticity. While evolutionary divergence occurs through DNA sequence differences among populations, plastic divergence among populations may be generated by changes in the epigenome. Here, we present the results of a genome‐wide comparison of DNA methylation patterns and genetic structure among four populations of Eastern oyster (Crassostrea virginica) in the northern Gulf of Mexico. We used a combination of restriction site‐associated DNA sequencing (RADseq) and reduced representation bisulfite sequencing (RRBS) to explore population structure, gene‐wide averages of FST, and DNA methylation differences between oysters inhabiting four estuaries with unique salinity profiles. This approach identified significant population structure despite a moderately low FST (0.02) across the freshwater boundary of the Mississippi river, a finding that may reflect recent efforts to restore oyster stock populations. Divergence between populations in CpG methylation was greater than for divergence in FST, likely reflecting environmental effects on DNA methylation patterns. Assessment of CpG methylation patterns across all populations identified that only 26% of methylated DNA was intergenic; and, only 17% of all differentially methylated regions (DMRs) were within these same regions. DMRs within gene bodies between sites were associated with genes known to be involved in DNA damage repair, ion transport, and reproductive timing. Finally, when assessing the correlation between genomic variation and DNA methylation between these populations, we observed population‐specific DNA methylation profiles that were not directly associated with single nucleotide polymorphisms or broader gene‐body mean FST trends. Our results suggest that C. virginica may use DNA methylation to generate environmentally responsive plastic phenotypes and that there is more divergence in methylation than divergence in allele frequencies.

Species distributions are rapidly changing as human globalization increasingly moves organisms to novel environments. In marine systems, species introductions are the result of a number of anthropogenic mechanisms, notably shipping, aquaculture/mariculture, the pet and bait trades, and the creation of canals. Marine invasions are a global threat to human and non‐human populations alike and are often listed as one of the top conservation concerns worldwide, having ecological, evolutionary, and social ramifications. Evolutionary investigations of marine invasions can provide crucial insight into an introduced species’ potential impacts in its new range, including: physiological adaptation and behavioral changes to exploit new environments; changes in resident populations, community interactions, and ecosystems; and severe reductions in genetic diversity that may limit evolutionary potential in the introduced range. This special issue focuses on current research advances in the evolutionary biology of marine invasions and can be broadly classified into a few major avenues of research: the evolutionary history of invasive populations, post‐invasion reproductive changes, and the role of evolution in parasite introductions. Together, they demonstrate the value of investigating marine invasions from an evolutionary perspective, with benefits to both fundamental and applied evolutionary biology at local and broad scales.

Understanding the mechanisms involved in pathogen adaptation to quantitative resistance in plants has a key role to play in establishing durable strategies for resistance deployment, especially in perennial crops. The erosion of quantitative resistance has been recently suspected in Cuba and the Dominican Republic for a major fungal pathogen of such a crop: Pseudocercospora fijiensis, causing black leaf streak disease on banana. This study set out to test whether such erosion has resulted from an adaptation of P. fijiensis populations, and to determine whether or not the adaptation is local. Almost 600 P. fijiensis isolates from Cuba and the Dominican Republic were sampled using a paired‐population sampling design on resistant and susceptible banana varieties. A low genetic structure of the P. fijiensis populations was detected in each country using 16 microsatellite markers. Cross‐inoculation experiments using isolates from susceptible and resistant cultivars were carried out, measuring a quantitative trait (the diseased leaf area) related to pathogen fitness on three varieties. A further analysis based on those data suggested the existence of a local pattern of adaptation to resistant cultivars in both of the study countries, due to the existence of specific (or genotype by genotype) host–pathogen interactions. However, neither cost nor benefit effects for adapted populations were found on the widely used “Cavendish” banana group. These results highlight the need to study specific host–pathogen interactions and pathogen adaptation on a wide range of quantitative resistance phenotypes in banana, in order to develop durable strategies for resistance deployment.

Incorporating species' eco‐evolutionary responses to human‐caused disturbances remains a challenge in marine management efforts. A prerequisite is knowledge of geographic structure and scale of genetic diversity and connectivity ‐ the so‐called seascape genetic patterns. The Baltic Sea is an excellent model system for studies linking seascape genetics with effects of anthropogenic stress. However, seascape genetic patterns in this area are only described for a few species and are completely unknown for invertebrate herbivores, which constitute a critical part of the ecosystem. This information is crucial for sustainable management, particularly under future scenarios of rapid environmental change. Here, we investigate the population genetic structure among 31 locations throughout the Baltic Sea, of which 45% were located in marine protected areas, in one of the most important herbivores of this region, the isopod crustacean Idotea balthica, using an array of 33,774 genome‐wide SNP markers derived from 2b‐RAD sequencing. In addition, we generate a biophysical connectivity matrix for I. balthica from a combination of oceanographic current models and estimated life history traits. We find population structure on scales of hundreds of kilometers across the Baltic Sea, where genomic patterns in most cases closely match biophysical connectivity, indicating passive transport with oceanographic currents as an important mean of dispersal in this species. We also find a reduced genetic diversity in terms of heterozygosity along the main salinity gradient of the Baltic Sea, suggesting periods of low population size. Our results provide crucial information for management of a key ecosystem species under expected changes in temperature and salinity following global climate change in a marine coastal area.

Forest ecosystems provide important ecological services and resources, from habitat for biodiversity to the production of environmentally friendly products, and play a key role in the global carbon cycle. Humanity is counting on forests to sequester and store a substantial portion of the anthropogenic carbon dioxide produced globally. However, the unprecedented rate of climate change, deforestation, and accidental importation of invasive insects and diseases are threatening the health and productivity of forests, and their capacity to provide these services. Knowledge of genetic diversity, local adaptation, and genetic control of key traits is required to predict the adaptive capacity of tree populations, inform forest management and conservation decisions, and improve breeding for productive trees that will withstand the challenges of the 21st century. Genomic approaches have well accelerated the generation of knowledge of the genetic and evolutionary underpinnings of nonmodel tree species, and advanced their applications to address these challenges. This special issue of Evolutionary Applications features 14 papers that demonstrate the value of a wide range of genomic approaches that can be used to better understand the biology of forest trees, including species that are widespread and managed for timber production, and others that are threatened or endangered, or serve important ecological roles. We highlight some of the major advances, ranging from understanding the evolution of genomes since the period when gymnosperms separated from angiosperms 300 million years ago to using genomic selection to accelerate breeding for tree health and productivity. We also discuss some of the challenges and future directions for applying genomic tools to address long‐standing questions about forest trees.

Human‐mediated transport of species outside their natural range is a rapidly growing threat to biodiversity, particularly for island ecosystems which have evolved in isolation. The genetic structure underpinning island populations will largely determine their response to increased transport and thus help to inform biosecurity management. However, this information is severely lacking for some groups, such as the soil fauna. We therefore analysed the phylogeographic structure of an indigenous and an invasive springtail species (Collembola: Poduromorpha), each distributed across multiple remote sub‐Antarctic islands, where human activity is currently intensifying. For both species, we generated a genome‐wide SNP dataset and additionally analysed all available COI barcodes. Genetic differentiation in the indigenous springtail Tullbergia bisetosa is substantial among (and, to a lesser degree, within) islands, reflecting low dispersal and historic population fragmentation, while COI patterns reveal ancestral signatures of post‐glacial recolonisation. This pronounced geographic structure demonstrates the key role of allopatric divergence in shaping the region's diversity and highlights the vulnerability of indigenous populations to genetic homogenisation via human transport. For the invasive species Hypogastrura viatica, nuclear genetic structure is much less apparent, particularly for islands linked by regular shipping, while diverged COI haplotypes indicate multiple independent introductions to each island. Thus, human transport has likely facilitated this species’ persistence since its initial colonisation, through the ongoing introduction and inter‐island spread of genetic variation. These findings highlight the different evolutionary consequences of human transport for indigenous and invasive soil species. Yet both demonstrate the need for improved intraregional biosecurity among remote island systems, where the policy focus to date has been on external introductions.

We review a recently discovered white spruce (Picea glauca) chemical defense against spruce budworm (Choristoneura fumiferana) involving hydroxyacetophenones. These defense metabolites detected in the foliage accumulate variably as the aglycons, piceol and pungenol, or the corresponding glucosides, picein and pungenin. We summarize current knowledge of the genetic, genomic, molecular, and biochemical underpinnings of this defense and its effects on C. fumiferana. We present an update with new results on the ontogenic variation and the phenological window of this defense, including analysis of transcript responses in P. glauca to C. fumiferana herbivory. We also discuss this chemical defense from an evolutionary and a breeding context.

Rising ocean temperatures and extreme temperature events have precipitated declines and local extinctions in many marine species globally, but patterns of loss are often uneven across species ranges for reasons that are poorly understood. Knowledge of the extent of local adaptation and gene flow may explain such patterns and help predict future trajectories under scenarios of climate change. We test the extent to which local differentiation in thermal tolerance is influenced by gene flow and local adaptation using a widely distributed intertidal seaweed (Hormosira banksii) from temperate Australia. Population surveys across ~2,000 km of the species range revealed strong genetic structuring at regional and local scales (global FST = 0.243) reflecting extremely limited gene flow, while common garden experiments (14‐day exposures to 15, 18, 21 °C) revealed strong site differences in early development and mortality in response to elevated temperature. Embryos from many sites spanning a longitudinal thermal gradient showed suppressed development and increased mortality to elevated water temperatures, but populations originating from warmer and more variable thermal environments tended to be less susceptible to warming. Notably, there was significant local scale variation in the thermal responses of embryos within regions which was corroborated by the finding of small‐scale genetic differences. We expect the observed genetic and phenotypic differentiation to lead to uneven responses to warming sea surface temperatures in this important marine foundation species. The study highlights the challenges of predicting species responses to thermal stress and the importance of management strategies that incorporate evolutionary potential for “climate‐proofing” marine ecosystems.

Improving our sparse knowledge of the mating and reproductive behaviour of white rhinoceros (Ceratotherium simum Burchell, 1817) is essential for the effective conservation of this iconic species. By combining morphological, physiological and habitat data with paternity assignments of 104 known mother–offspring pairs collected over a period of 13 years, we provide the most comprehensive analysis of the mating system in this species. We show that while the overall mating system was promiscuous, and both males and females produced more offspring when mating with several partners, half of all females with multiple offspring were monogamous. Additionally, we find that mating and reproductive success varied significantly among territorial males in two independent sets of males. In females, however, variation in the mating and the reproductive success was not larger than expected by random demographic fluctuations. Horn size, testosterone metabolite concentration, territory size, habitat openness and the volume of preferred food within the territory did not seem to influence male mating or reproductive success. Moreover, there was no sign of inbreeding avoidance: females tended to mate more frequently with closely related males, and one daughter produced a progeny with her father. The lack of inbreeding avoidance, in combination with the skew in male reproductive success, the partial monogamy in females and the territorial‐based mating system, jeopardizes the already low genetic variation in the species. Considering that the majority of populations are restricted to fenced reserves and private farms, we recommend taking preventive measures that aim to reduce inbreeding in white rhinoceros. A video abstract can be viewed here.

In bacteria, evolution of resistance to one antibiotic is frequently associated with increased resistance (cross‐resistance) or increased susceptibility (collateral sensitivity) to other antibiotics. Cross‐resistance and collateral sensitivity are typically evaluated at the minimum inhibitory concentration (MIC). However, these susceptibility changes are not well characterized with respect to the mutant prevention concentration (MPC), the antibiotic concentration that prevents a single‐step mutation from occurring.

Persistent insecticides sprayed onto house walls, and incorporated into insecticide‐treated bednets, provide long‐acting, cost‐effective control of vector‐borne diseases such as malaria and leishmaniasis. The high concentrations that occur immediately postdeployment may kill both resistant and susceptible insects. However, insecticide concentration, and therefore killing ability, declines in the months after deployment. As concentrations decline, resistant insects start to survive, while susceptible insects are still killed. The period of time after deployment, within which the mortality of resistant individuals is lower than that of susceptible ones, has been termed the “window of selection” in other contexts. It is recognized as driving resistance in bacteria and malaria parasites, both of which are predominantly haploid. We argue that paying more attention to these mortality differences can help understand the evolution of insecticide resistance. Because insects are diploid, resistance encoded by single genes generates heterozygotes. This gives the potential for a narrower “window of dominance,” within the window of selection, where heterozygote mortality is lower than that of susceptible homozygotes. We explore the general properties of windows of selection and dominance in driving resistance. We quantify their likely effect using data from new laboratory experiments and published data from the laboratory and field. These windows can persist months or years after insecticide deployments. Differential mortalities of resistant, susceptible and heterozygous genotypes, after public health deployments, constitute a major challenge to controlling resistance. Greater attention to mortality differences by genotype would inform strategies to reduce the evolution of resistance to existing and new insecticides.

Identifying the sources of ongoing and novel disease outbreaks is critical for understanding the diffusion of epizootic diseases. Identifying infection sources is difficult when few physical differences separate individuals with different origins. Genetic assignment procedures show great promise for assessing transmission dynamics in such situations. Here, we use genetic assignment tests to determine the source of chronic wasting disease infections in free‐ranging white‐tailed deer (Odocoileus virginianus) populations. Natural dispersal is thought to facilitate the geographic diffusion of chronic wasting disease, but egression from captive cervid populations represents an alternative source of infection that is difficult to detect due to physical similarities with wild deer. Simulated reference populations were created based on allele frequencies from 1,912 empirical microsatellite genotypes collected in four sampling subregions and five captive facilities. These reference populations were used to assess the likelihood of ancestry and assignment of 1,861 free‐ranging deer (1,834 noninfected and 27 infected) and 51 captive individuals to captive or wild populations. The ancestry (Q) and assignment scores (A) for free‐ranging deer to wild populations were high (average Qwild = 0.913 and average Awild = 0.951, respectively), but varied among subregions (Qwild = 0.800–0.947, Awild = 0.857–0.976). These findings suggest that captive egression and admixture are rare, but risk may not be spatially uniform. Ancestry and assignment scores for two free‐ranging deer with chronic wasting disease sampled in an area where chronic wasting disease was previously unobserved in free‐ranging herds indicated a higher likelihood of assignment and proportion of ancestry attributable to captive populations. While we cannot directly assign these individuals to infected facilities, these findings suggest that rare egression events may influence the epizootiology of chronic wasting disease in free‐ranging populations. Continued disease surveillance and genetic analyses may further elucidate the relative disease risk attributable to captive and wild sources.

An understanding of how genetic variability is distributed in space is fundamental for the conservation and maintenance of diversity in spatially fragmented and vulnerable populations. While fragmentation can occur from natural barriers, it can also be exacerbated by anthropogenic activities such as hydroelectric power plant development. Whatever the source, fragmentation can have significant ecological effects, including disruptions of migratory processes and gene flow among populations. In Chile, the Biobío River basin exhibits a high degree of habitat fragmentation due to the numerous hydroelectric power plants in operation, the number of which is expected to increase following new renewable energy use strategies. Here, we assessed the effects of different kinds of barriers on the genetic structure of the endemic freshwater fish Percilia irwini, knowledge that is critically needed to inform conservation strategies in light of current and anticipated further fragmentation initiatives in the system. We identified eight genetic units throughout the entire Biobío system with high effective sizes. A reduced effective size estimate was, however, observed in a single population located between two impassable barriers. Both natural waterfalls and human‐made dams were important drivers of population differentiation in this system; however, dams affect genetic diversity differentially depending on their mode of operation. Evidence of population extirpation was found in two river stretches limited by upstream and downstream dams. Significant gene flow in both directions was found among populations not separated by natural or anthropogenic barriers. Our results suggest a significant vulnerability of P. irwini populations to future dam development and demonstrate the importance of studying basin‐wide data sets with genetic metrics to understand the strength and direction of anthropogenic impacts on fish populations.

Our knowledge of the diversity of potato cyst nematodes in their native areas still remains patchy and should be improved. A previous study based on 42 Peruvian Globodera pallida populations revealed a clear south to north phylogeographic pattern, with five well‐supported clades and maximum diversity observed in the south of Peru. In order to investigate this phylogeographic pattern more closely, we genotyped a larger collection of Peruvian populations using both cathepsin L gene sequence data and a new set of 13 microsatellite loci. Using different genetic analyses (STRUCTURE, DAPC), we consistently obtained the same results that led to similar conclusions: the presence of a larger genetic diversity than previously known suggesting the presence of cryptic species in the south of Peru. These investigations also allowed us to clarify the geographic borders of the previously described G. pallida genetic clades and to update our knowledge of the genetic structure of this species in its native area, with the presence of additional clades. A distance‐based redundancy analysis (dbRDA) was also carried to understand whether there was a correlation between the population genetic differentiation and environmental conditions. This analysis showed that genetic distances observed between G. pallida populations are explained firstly by geographic distances, but also by climatic and soil conditions. This work could lead to a revision of the taxonomy that may have strong implications for risk assessment and management, especially on a quarantine species.

Research on the genetics of invasive species often focuses on patterns of genetic diversity and population structure within the introduced range. However, a growing body of literature is demonstrating the need to study how native genotypes affect both ecological and evolutionary mechanisms within the introduced range. Here, we used genotyping‐by‐sequencing to study both native and introduced ranges [based on 1,653 single nucleotide polymorphisms (SNPs)] of the amphiatlantic marine invertebrate Ciona intestinalis. A previous study using microsatellites analysed samples collected along the Swedish west coast and showed the presence of genetically distinct lineages in deep and shallow waters. Using 1,653 single nucleotide polymorphisms (SNPs) from newly collected samples (285 individuals), we first confirmed the presence of this depth‐defined genomic divergence along the Swedish coast. We then used approximate Bayesian computation to infer the historical relationship among sites from the North Sea, the English Channel and the northwest Atlantic and found evidence of ancestral divergence between individuals from deep waters off Sweden and individuals from the English Channel. This divergence was followed by a secondary contact that led to a genetic admixture between the ancestral populations (i.e., deep Sweden and English Channel), which originated the genotypes found in shallow Sweden. We then revealed that the colonization of C. intestinalis in the northwest Atlantic was as a result of an admixture between shallow Sweden and the English Channel genotypes across the introduced range. Our results showed the presence of both past and recent genetic admixture events that together may have promoted the successful colonizations of C. intestinalis. Our study suggests that secondary contacts potentially reshape the evolutionary trajectories of invasive species through the promotion of intraspecific hybridization and by altering both colonization patterns and their ecological effects in the introduced range.

Temperature plays a multidimensional role in host–pathogen interactions. As an important element of climate change, elevated world temperature resulting from global warming presents new challenges to sustainable disease management. Knowledge of pathogen adaptation to global warming is needed to predict future disease epidemiology and formulate mitigating strategies. In this study, 21 Phytophthora infestans isolates originating from seven thermal environments were acclimated for 200 days under stepwise increase or decrease of experimental temperatures and evolutionary responses of the isolates to the thermal changes were evaluated. We found temperature acclimation significantly increased the fitness and genetic adaptation of P. infestans isolates at both low and high temperatures. Low‐temperature acclimation enforced the countergradient adaptation of the pathogen to its past selection and enhanced the positive association between the pathogen's intrinsic growth rate and aggressiveness. At high temperatures, we found that pathogen growth collapsed near the maximum temperature for growth, suggesting a thermal niche boundary may exist in the evolutionary adaptation of P. infestans. These results indicate that pathogens can quickly adapt to temperature shifts in global warming. If this is associated with environmental conditions favoring pathogen spread, it will threaten future food security and human health and require the establishment of mitigating actions.

The adaptive potential of invasive species is related to the genetic diversity of the invader, which is influenced by genetic drift and natural selection. Typically, the genetic diversity of invaders is studied with neutral genetic markers; however, the expectation of reduced diversity has not been consistently supported by empirical studies. Here, we describe and interpret genetic diversity at both neutral microsatellite loci and the immune‐related MHC‐DRB locus of native and invasive populations of raccoon to better understand of how drift and selection impact patterns of genetic diversity during the invasion process. We found that despite the loss of many MHC (major histocompatibility complex) alleles in comparison with native populations, functional MHC supertypes are preserved in the invasive region. In the native raccoon population, the number of supertypes within individuals was higher than expected under a neutral model. The high level of individual functional divergence may facilitate the adaptation to local conditions in the invasive range. In the invasive populations, we also detected increased population structure at microsatellites compared to the MHC locus, further suggesting that balancing selection is acting on adaptively important regions of the raccoon genome. Finally, we found that alleles known to exhibit resistance to rabies in the native range, Prlo‐DRB*4, Prlo‐DRB*16 and Prlo‐DRB*102, were the most common alleles in the European populations, suggesting directional selection is acting on this locus. Our research shows empirical support for the importance of functional immune diversity for adaptation and survival in novel environments.

Organophosphate (OP) and carbamate (CM) insecticides are widely used in the United States and share the same mode of toxic action. Both classes are frequently documented in aquatic ecosystems, sometimes at levels that exceed aquatic life benchmarks. We previously identified a population of the nontarget amphipod, Hyalella azteca, thriving in an agricultural creek with high sediment levels of the OP chlorpyrifos, suggesting the population may have acquired genetic resistance to the pesticide. In the present study, we surveyed 17 populations of H. azteca in California to screen for phenotypic resistance to chlorpyrifos as well as genetic signatures of resistance in the acetylcholinesterase (ace‐1) gene. We found no phenotypic chlorpyrifos resistance in populations from areas with little or no pesticide use. However, there was ~3‐ to 1,000‐fold resistance in H. azteca populations from agricultural and/or urban areas, with resistance levels in agriculture being far higher than urban areas due to greater ongoing use of OP and CM pesticides. In every case of resistance in H. azteca, we identified a glycine‐to‐serine amino acid substitution (G119S) that has been shown to confer OP and CM resistance in mosquitoes and has been associated with resistance in other insects. We found that the G119S mutation was always present in a heterozygous state. Further, we provide tentative evidence of an ace‐1 gene duplication in H. azteca that may play a role in chlorpyrifos resistance in some populations. The detection of a genetically based, adaptive OP and CM resistance in some of the same populations of H. azteca previously shown to harbor a genetically based adaptive pyrethroid resistance indicates that these nontarget amphipod populations have become resistant to many of the insecticides now in common use. The terrestrial application of pesticides has provided strong selective pressures to drive evolution in a nontarget, aquatic species.

The ultimate causes of correlated evolution among sites in a genome remain difficult to teaseapart. To address this problem directly, we performed a high‐throughput search for correlated evolution among sites associated with resistance to a fluoroquinolone antibiotic using whole genome data from clinical strains of Pseudomonas aeruginosa, before validating our computational predictions experimentally. We show that for at least two sites, this correlation is underlain by epistasis. Our analysis also revealed eight additional pairs of synonymous substitutions displaying correlated evolution underlain by physical linkage, rather than selection associated with antibiotic resistance. Our results provide direct evidence that both epistasis and physical link age among sites can drive the correlated evolution identified by high throughput computational tools. In other words, the observation of correlated evolution is not by itself sufficient evidence to guarantee that the sites in question are epistatic; such a claim requires additional evidence, ideally coming from direct estimates of epistasis, based on experimental evidence.

When populations of a rare species are small, isolated and declining under climate change, some populations may become locally maladapted. Detecting this maladaptation may allow effective rapid conservation interventions, even if based on incomplete knowledge. Population maladaptation may be estimated by finding genome–environment associations (GEA) between allele frequencies and environmental variables across a local species range, and identifying populations whose allele frequencies do not fit with these trends. We can then design assisted gene flow strategies for maladapted populations, to adjust their allele frequencies, entailing lower levels of intervention than with undirected conservation action. Here, we investigate this strategy in Scottish populations of the montane plant dwarf birch (Betula nana). In genome‐wide restriction site‐associated single nucleotide polymorphism (SNP) data, we found 267 significant associations between SNP loci and environmental variables. We ranked populations by maladaptation estimated using allele frequency deviation from the general trends at these loci; this gave a different prioritization for conservation action than the Shapely Index, which seeks to preserve rare neutral variation. Populations estimated to be maladapted in their allele frequencies at loci associated with annual mean temperature were found to have reduced catkin production. Using an environmental niche modelling (ENM) approach, we found annual mean temperature (35%), and mean diurnal range (15%), to be important predictors of the dwarf birch distribution. Intriguingly, there was a significant correlation between the number of loci associated with each environmental variable in the GEA and the importance of that variable in the ENM. Together, these results suggest that the same environmental variables determine both adaptive genetic variation and species range in Scottish dwarf birch. We suggest an assisted gene flow strategy that aims to maximize the local adaptation of dwarf birch populations under climate change by matching allele frequencies to current and future environments.

Genome‐wide SNP profiling has yielded insights into the genetic structure of China indigenous pigs, but has focused on a limited number of populations. Here, we present an analysis of population structure and signals of positive selection in 42 Chinese pig populations that represent the most extensive pig phenotypic diversity in China, using genotype data of 1.1 million SNPs on customized Beadchips. This unravels the fine‐scale genetic diversity, phylogenic relationships, and population structure of these populations, which shows remarkably concordance between genetic clusters and geography with few exceptions. We also reveal the genetic contribution to North Chinese pigs from European modern pigs. Furthermore, we identify possible targets of selection in the Tibetan pig, including the well‐characterized hypoxia gene (EPAS1) and several previously unrecognized candidates. Intriguingly, the selected haplotype in the EPAS1 gene is associated with higher hemoglobin contents in Tibetan pigs, which is different from the protective role of EPAS1 in the high‐altitude adaptation in Tibetan dogs and their owners. Additionally, we present evidence for the causality between EDNRB variants and the two‐end‐black (TEB) coat color phenotype in all Chinese pig populations except the Jinhua pig. We hypothesize that distinct targets have been independently selected for the formation of the TEB phenotype in Chinese pigs of different geographic origins. This highlights the importance of characterizing population‐specific genetic determinants for heritable phenotype in diverse pig populations.

Understanding and quantifying populations' adaptive genetic variation and their response to climate change are critical to reforestation's seed source selection, forest management decisions, and gene conservation. Landscape genomics combined with geographic and environmental information provide an opportunity to interrogate forest populations' genome‐wide variation for understanding the extent to which evolutionary forces shape past and contemporary populations' genetic structure, and identify those populations that may be most at risk under future climate change. Here, we used genotyping by sequencing to generate over 11,000 high‐quality variants from Platycladus orientalis range‐wide collection to evaluate its diversity and to predict genetic offset under future climate scenarios. Platycladus orientalis is a widespread conifer in China with significant ecological, timber, and medicinal values. We found population structure and evidences of isolation by environment, indicative of adaptation to local conditions. Gradient forest modeling identified temperature‐related variables as the most important environmental factors influencing genetic variation and predicted areas with higher risk under future climate change. This study provides an important reference for forest resource management and conservation for P. orientalis.

Human‐mediated transport creates secondary contacts between genetically differentiated lineages, bringing new opportunities for gene exchange. When similar introductions occur in different places, they provide informally replicated experiments for studying hybridisation. We here examined 4,279 Mytilus mussels, sampled in Europe and genotyped with 77 ancestry‐informative markers. We identified a type of introduced mussels, called “dock mussels,” associated with port habitats and displaying a particular genetic signal of admixture between M. edulis and the Mediterranean lineage of M. galloprovincialis. These mussels exhibit similarities in their ancestry compositions, regardless of the local native genetic backgrounds and the distance separating colonised ports. We observed fine‐scale genetic shifts at the port entrance, at scales below natural dispersal distance. Such sharp clines do not fit with migration‐selection tension zone models, and instead suggest habitat choice and early‐stage adaptation to the port environment, possibly coupled with connectivity barriers. Variations in the spread and admixture patterns of dock mussels seem to be influenced by the local native genetic backgrounds encountered. We next examined departures from the average admixture rate at different loci, and compared human‐mediated admixture events, to naturally admixed populations and experimental crosses. When the same M. galloprovincialis background was involved, positive correlations in the departures of loci across locations were found; but when different backgrounds were involved, no or negative correlations were observed. While some observed positive correlations might be best explained by a shared history and saltatory colonisation, others are likely produced by parallel selective events. Altogether, genome‐wide effect of admixture seems repeatable and more dependent on genetic background than environmental context. Our results pave the way towards further genomic analyses of admixture, and monitoring of the spread of dock mussels both at large and at fine spacial scales.

Genomic methodologies offer unprecedented opportunities for statistically robust studies of species broadly distributed in environments conducive to high gene flow, providing valuable information for wildlife conservation and management. Here, we sequence restriction site‐associated DNA to characterize genome‐wide single nucleotide polymorphisms (SNPs) in a broadly distributed and highly migratory large pelagic fish, striped marlin (Kajikia audax). Assessment of over 4,000 SNPs resolved spatiotemporal patterns of genetic connectivity throughout the species range in the Pacific and, for the first time, Indian oceans. Individual‐based cluster analyses identified six genetically distinct populations corresponding with the western Indian, eastern Indian, western South Pacific, and eastern central Pacific oceans, as well as two populations in the North Pacific Ocean (FST = 0.0137–0.0819). FST outlier analyses identified a subset of SNPs (n = 59) putatively under the influence of natural selection and capable of resolving populations separated by comparatively high degrees of genetic differentiation. Temporal collections available for some regions demonstrated the stability of allele frequencies over three to five generations of striped marlin. Relative migration rates reflected lower levels of genetic connectivity between Indian Ocean populations (mR ≤ 0.37) compared with most populations in the Pacific Ocean (mR ≥ 0.57) and highlight the importance of the western South Pacific in facilitating gene flow between ocean basins. Collectively, our results provide novel insights into rangewide population structure for striped marlin and highlight substantial inconsistencies between genetically distinct populations and stocks currently recognized for fisheries management. More broadly, we demonstrate that species capable of long‐distance dispersal in environments lacking obvious physical barriers to movement can display substantial population subdivision that persists over multiple generations and that may be facilitated by both neutral and adaptive processes. Importantly, surveys of genome‐wide markers enable inference of population‐level relationships using sample sizes practical for large pelagic fishes of conservation concern.

Secondary contact may have important implications for ecological and evolutionary processes; however, few studies have tracked the outcomes of secondary contact from its onset in natural ecosystems. We evaluated an anadromous alewife (Alosa pseudoharengus) reintroduction project in Rogers Lake (Connecticut, USA), which contains a landlocked alewife population that was isolated as a result of colonial‐era damming. After access to the ocean was restored, adult anadromous alewife were stocked into the lake. We assessed anadromous juvenile production, the magnitude and direction of introgression, and the potential for competition between ecotypes. We obtained fin clips from all adult alewife stocked into the lake during the restoration and a sample of juveniles produced in the lake two years after the stocking began. We assessed the ancestry of juveniles using categorical assignment and pedigree reconstruction with newly developed microhaplotype genetic markers. Anadromous alewives successfully spawned in the lake and hybridized with the landlocked population. Parentage assignments revealed that male and female anadromous fish contributed equally to juvenile F1 hybrids. The presence of landlocked backcrosses shows that some hybrids were produced within the first two years of secondary contact, matured in the lake, and reproduced. Therefore, introgression appears directional, from anadromous into landlocked, in the lake environment. Differences in estimated abundance of juveniles of different ecotypes in different habitats were also detected, which may reduce competition between ecotypes as the restoration continues. Our results illustrate the utility of restoration projects to study the outcomes of secondary contact in real ecosystems.

Environmental DNA (eDNA) extracted from water samples has recently shown potential as a valuable source of population genetic information for aquatic macroorganisms. This approach offers several potential advantages compared with conventional tissue‐based methods, including the fact that eDNA sampling is noninvasive and generally more cost‐efficient. Currently, eDNA approaches have been limited to single‐marker studies of mitochondrial DNA (mtDNA), and the relationship between eDNA haplotype composition and true haplotype composition still needs to be thoroughly verified. This will require testing of bioinformatic and statistical software to correct for erroneous sequences, as well as biases and random variation in relative sequence abundances. However, eDNA‐based population genetic methods have far‐reaching potential for both basic and applied research. In this paper, we present a brief overview of the achievements of eDNA‐based population genetics to date, and outline the prospects for future developments in the field, including the estimation of nuclear DNA (nuDNA) variation and epigenetic information. We discuss the challenges associated with eDNA samples as opposed to those of individual tissue samples and assess whether eDNA might offer additional types of information unobtainable with tissue samples. Lastly, we provide recommendations for determining whether an eDNA approach would be a useful and suitable choice in different research settings. We limit our discussion largely to contemporary aquatic systems, but the advantages, challenges, and perspectives can to a large degree be generalized to eDNA studies with a different spatial and temporal focus.

Gymnosperms diverged from their sister plant clade of flowering plants 300 Mya. Morphological and functional divergence between the two major seed plant clades involved significant changes in their reproductive biology, water‐conducting systems, secondary metabolism, stress defense mechanisms, and small RNA‐mediated epigenetic silencing. The relatively recent sequencing of several gymnosperm genomes and the development of new genomic resources have enabled whole‐genome comparisons within gymnosperms, and between angiosperms and gymnosperms. In this paper, we aim to understand how genes and gene families have contributed to the major functional and morphological differences in gymnosperms, and how this information can be used for applied breeding and biotechnology. In addition, we have analyzed the angiosperm versus gymnosperm evolution of the pleiotropic drug resistance (PDR) gene family with a wide range of functionalities in plants' interaction with their environment including defense mechanisms. Some of the genes reviewed here are newly studied members of gene families that hold potential for biotechnological applications related to commercial and pharmacological value. Some members of conifer gene families can also be exploited for their potential in phytoremediation applications.

Detecting the molecular basis of local adaptation and identifying selective drivers is still challenging in nonmodel species. The use of purely population genetic approaches is limited by some characteristics of genetic systems, such as pleiotropy and polygenic control, and parallel evidence from phenotypic‐based experimental comparisons is required. In long‐lived organisms, the detection of selective pressures might also be precluded by evolutionary lag times in response to the environment. Here, we used the English yew to showcase an example of a multiscale integrative approach in a nonmodel species with limited plant and genomic resources. We combined information from two independent sources, phenotypes in a common environment and genomic data in natural populations, to investigate the signature of selection. Growth differences among populations in a common environment, and phenological patterns of both shoot elongation and male strobili maturation, were associated with climate clines, providing evidence for local adaptation and guiding us in the selection of populations for genomic analyses. We used information on over 25,000 SNPs from c. 1,200 genes to infer the demographic history and to test for molecular signatures of selection at different levels: SNP, gene, and biological pathway. Our results confirmed an overall demographic history of population decline, but we also found evidence for putative local adaptation at the molecular level. We identified or confirmed several candidate genes for positive and negative selection in forest trees, including the pseudo‐response regulator 7 (PRR7), an essential component of the circadian clock in plants. In addition, we successfully tested an approach to detect polygenic adaptation in biological pathways, allowing us to identify the flavonoid biosynthesis pathway as a candidate stress‐response pathway that deserves further attention in other plants. Finally, our study contributes to the emerging view that explaining contemporary standing genetic variation requires considering adaptation to past climates, especially for long‐lived trees.

Pinus sylvestris has a long history of basic and applied research that is relevant for both forestry and evolutionary studies. Its patterns of adaptive variation and role in forest economic and ecological systems have been studied extensively for nearly 275 years, detailed demography for a 100 years and mating system more than 50 years. However, its reference genome sequence is not yet available and genomic studies have been lagging compared to, for example, Pinus taeda and Picea abies, two other economically important conifers. Despite the lack of reference genome, many modern genomic methods are applicable for a more detailed look at its biological characteristics. For example, RNA‐seq has revealed a complex transcriptional landscape and targeted DNA sequencing displays an excess of rare variants and geographically homogenously distributed molecular genetic diversity. Current DNA and RNA resources can be used as a reference for gene expression studies, SNP discovery, and further targeted sequencing. In the future, specific consequences of the large genome size, such as functional effects of regulatory open chromatin regions and transposable elements, should be investigated more carefully. For forest breeding and long‐term management purposes, genomic data can help in assessing the genetic basis of inbreeding depression and the application of genomic tools for genomic prediction and relatedness estimates. Given the challenges of breeding (long generation time, no easy vegetative propagation) and the economic importance, application of genomic tools has a potential to have a considerable impact. Here, we explore how genomic characteristics of P. sylvestris, such as rare alleles and the low extent of linkage disequilibrium, impact the applicability and power of the tools.

Plantation‐grown trees have to cope with an increasing pressure of pest and disease in the context of climate change, and breeding approaches using genomics may offer efficient and flexible tools to face this pressure. In the present study, we targeted genetic improvement of resistance of an introduced conifer species in Canada, Norway spruce (Picea abies (L.) Karst.), to the native white pine weevil (Pissodes strobi Peck). We developed single‐ and multi‐trait genomic selection (GS) models and selection indices considering the relationships between weevil resistance, intrinsic wood quality, and growth traits. Weevil resistance, acoustic velocity as a proxy for mechanical wood stiffness, and average wood density showed moderate‐to‐high heritability and low genotype‐by‐environment interactions. Weevil resistance was genetically positively correlated with tree height, height‐to‐diameter at breast height (DBH) ratio, and acoustic velocity. The accuracy of the different GS models tested (GBLUP, threshold GBLUP, Bayesian ridge regression, BayesCπ) was high and did not differ among each other. Multi‐trait models performed similarly as single‐trait models when all trees were phenotyped. However, when weevil attack data were not available for all trees, weevil resistance was more accurately predicted by integrating genetically correlated growth traits into multi‐trait GS models. A GS index that corresponded to the breeders’ priorities achieved near maximum gains for weevil resistance, acoustic velocity, and height growth, but a small decrease for DBH. The results of this study indicate that it is possible to breed for high‐quality, weevil‐resistant Norway spruce reforestation stock with high accuracy achieved from single‐trait or multi‐trait GS.

A lack of optimal gene combinations, as well as low levels of genetic diversity, is often associated with the formation of species range margins. Conservation efforts rely on predictive modelling using abiotic variables and assessments of genetic diversity to determine target species and populations for controlled breeding, germplasm conservation and assisted migration. Biotic factors such as interspecific competition and hybridization, however, are largely ignored, despite their prevalence across diverse taxa and their role as key evolutionary forces. Hybridization between species with well‐developed barriers to reproductive isolation often results in the production of offspring with lower fitness. Generation of novel allelic combinations through hybridization, however, can also generate positive fitness consequences. Despite this possibility, hybridization‐mediated introgression is often considered a threat to biodiversity as it can blur species boundaries. The contribution of hybridization towards increasing genetic diversity of populations at range margins has only recently gathered attention in conservation studies. We assessed the extent to which hybridization contributes towards range dynamics by tracking spatio‐temporal changes in the central location of a hybrid zone between two recently diverged species of pines: Pinus strobiformis and P. flexilis. By comparing geographic cline centre estimates for global admixture coefficient with morphological traits associated with reproductive output, we demonstrate a northward shift in the hybrid zone. Using a combination of spatially explicit, individual‐based simulations and linkage disequilibrium variance partitioning, we note a significant contribution of adaptive introgression towards this northward movement, despite the potential for differences in regional population size to aid hybrid zone movement. Overall, our study demonstrates that hybridization between recently diverged species can increase genetic diversity and generate novel allelic combinations. These novel combinations may allow range margin populations to track favourable climatic conditions or facilitate adaptive evolution to ongoing and future climate change.

Future climate change has been predicted to disrupt local adaptation in many perennial plants, such as forest trees, but the magnitude and location of these effects are thus far poorly understood. Here, we assess local adaptation to current climate in European aspen (Populus tremula) by using environmental association analyses to identify genetic variants associated with two representative climate variables describing current day variation in temperature and precipitation. We also analysed patterns of genetic differentiation between southern and northern populations and observe that regions of high genetic differentiation are enriched for SNPs that are significantly associated with climate. Using variants associated with climate, we examined patterns of isolation by distance and environment and used spatial modelling to predict the geographic distribution of genomic variation in response to two scenarios of future climate change. We show that climate conditions at a northern reference site will correspond to climate conditions experienced by current day populations located 4–8 latitude degrees further south. By assessing the relationship between phenotypic traits and vegetative fitness, we also demonstrate that southern populations harbour genetic variation that likely would be adaptive further north under both climate change scenarios. Current day populations at the lagging edge of the distribution in Sweden can therefore serve as sources for introducing adaptive alleles onto northern populations, but the likelihood of this largely depends on naturally occurring levels of gene flow.

Identifying genetic variants responsible for phenotypic variation under selective pressure has the potential to enable productive gains in natural resource conservation and management. Despite this potential, identifying adaptive candidate loci is not trivial, and linking genotype to phenotype is a major challenge in contemporary genetics. Many of the population genetic approaches commonly used to identify adaptive candidates will simultaneously detect false positives, particularly in nonmodel species, where experimental evidence is seldom provided for putative roles of the adaptive candidates identified by outlier approaches. In this study, we use outcomes from population genetics, phenotype association, and gene expression analyses as multiple lines of evidence to validate candidate genes. Using lodgepole and jack pine as our nonmodel study species, we analyzed 17 adaptive candidate loci together with 78 putatively neutral loci at 58 locations across Canada (N > 800) to determine whether relationships could be established between these candidate loci and phenotype related to mountain pine beetle susceptibility. We identified two candidate loci that were significant across all population genetic tests, and demonstrated significant changes in transcript abundance in trees subjected to wounding or inoculation with the mountain pine beetle fungal associate Grosmannia clavigera. Both candidates are involved in central physiological processes that are likely to be invoked in a trees response to stress. One of these two candidate loci showed a significant association with mountain pine beetle attack status in lodgepole pine. The spatial distribution of the attack‐associated allele further coincides with other indicators of susceptibility in lodgepole pine. These analyses, in which population genetics was combined with laboratory and field experimental validation approaches, represent first steps toward linking genetic variation to the phenotype of mountain pine beetle susceptibility in lodgepole and jack pine, and provide a roadmap for more comprehensive analyses.

Hybridization dynamics between co‐occurring species in environments where human‐mediated changes take place are important to quantify for furthering our understanding of human impacts on species evolution and for informing management. The allis shad Alosa alosa (Linnaeus, 1758) and twaite shad Alosa fallax (Lacépède, 1803), two clupeids sister species, have been severely impacted by human activities across Europe. The shrinkage of A. alosa distribution range along with the decline of the remaining populations' abundance threatens its persistence. The main objective was to evaluate the extent of hybridization and introgression between those interacting species. We developed a set of 77 species‐specific SNP loci that allowed a better resolution than morphological traits as they enabled the detection of hybrids up to the third generation. Variable rates of contemporary hybridization and introgression patterns were detected in 12 studied sites across the French Atlantic coast. Mitochondrial markers revealed a cyto‐nuclear discordance almost invariably involving A. alosa individuals with an A. fallax mitochondrial DNA and provided evidence of historical asymmetric introgression. Overall, contemporary and historical introgression revealed by nuclear and mitochondrial markers strongly suggests that a transfer of genes occurs from A. fallax toward A. alosa genome since at least four generations. Moreover, the outcomes of introgression greatly depend on the catchments where local processes are thought to occur. Undoubtedly, interspecific interaction and gene flow should not be overlooked when considering the management of those species.

Hybridization is a natural process at species range boundaries, but increasing numbers of species are hybridizing due to direct or indirect human activities. In such cases of anthropogenic hybridization, subsequent introgression can threaten the survival of native species. To date, many such systems have been studied with too few genetic markers to assess the level of threat resulting from advanced backcrossing. Here, we use 44,999 single nucleotide polymorphisms (SNPs) and the ADMIXTURE program to study two areas of Scotland where a panel of 22 diagnostic microsatellites previously identified introgression between native red deer (Cervus elaphus) and introduced Japanese sika (Cervus nippon). In Kintyre, we reclassify 26% of deer from the pure species categories to the hybrid category whereas in the NW Highlands we only reclassify 2%. As expected, the reclassified individuals are mostly advanced backcrosses. We also investigate the ability of marker panels selected on different posterior allele frequency criteria to find hybrids assigned by the full marker set and show that in our data, ancestry informative markers (i.e. those that are highly differentiated between the species, but not fixed) are better than diagnostic markers (those markers that are fixed between the species) because they are more evenly distributed in the genome. Diagnostic loci are concentrated on the X chromosome to the detriment of autosomal coverage.

The root‐knot nematodes are the most devastating worms to worldwide agriculture with Meloidogyne incognita being the most widely distributed and damaging species. This parasitic and ecological success seems surprising given its supposed obligatory clonal reproduction. Clonal reproduction has been suspected based on cytological observations but, so far, never confirmed by population genomics data. As a species, M. incognita is highly polyphagous with thousands of host plants. However, different M. incognita isolates present distinct and overlapping patterns of host compatibilities. Historically, four “host races” had been defined as a function of ranges of compatible and incompatible plants. In this study, we used population genomics to assess whether (a) reproduction is actually clonal in this species, (b) the host races follow an underlying phylogenetic signal or, rather represent multiple independent transitions, and (c) how genome variations associate with other important biological traits such as the affected crops and geographical distribution. We sequenced the genomes of 11 M. incognita isolates across Brazil that covered the four host races in replicates. By aligning the genomic reads of these isolates to the M. incognita reference genome assembly, we identified point variations. Analysis of linkage disequilibrium and 4‐gametes test showed no evidence for recombination, corroborating the clonal reproduction of M. incognita. The few point variations between the isolates showed no significant association with the host races, the geographical origin of the samples, or the crop on which they have been collected. Addition of isolates from other locations around the world confirmed this lack of underlying phylogenetic signal. This suggests multiple gains and losses of parasitic abilities and adaptations to different environments account for the broad host spectrum and wide geographical distribution of M. incognita and thus to its high economic impact. This surprising adaptability without sex poses both evolutionary and agro‐economic challenges.

Sequencing technologies and bioinformatic approaches are now available to resolve the challenges associated with complex and heterozygous genomes. Increased access to less expensive and more effective instrumentation will contribute to a wealth of high‐quality plant genomes in the next few years. In the meantime, more than 370 tree species are associated with public projects in primary repositories that are interrogating expression profiles, identifying variants, or analyzing targeted capture without a high‐quality reference genome. Genomic data from these projects generates sequences that represent intermediate assemblies for transcriptomes and genomes. These data contribute to forest tree biology, but the associated sequence remains trapped in supplemental files that are poorly integrated in plant community databases and comparative genomic platforms. Successful implementation of life science cyberinfrastructure is improving data standards, ontologies, analytic workflows, and integrated database platforms for both model and non‐model plant species. Unique to forest trees with large populations that are long‐lived, outcrossing, and genetically diverse, the phenotypic and environmental metrics associated with georeferenced populations are just as important as the genomic data sampled for each individual. To address questions related to forest health and productivity, cyberinfrastructure must keep pace with the magnitude of genomic and phenomic sampling of larger populations. This review examines the current landscape of cyberinfrastructure, with an emphasis on best practices and resources to align community data with the Findable, Accessible, Interoperable, and Reusable (FAIR) guidelines.

The salmonid fish Brown trout is iconic as a model for the application of conservation genetics to understand and manage local interspecific variation. However, there is still scant information about relationships between local and large‐scale population structure, and to what extent geographical and environmental variables are associated with barriers to gene flow. We used information from 3,782 mapped SNPs developed for the present study and conducted outlier tests and gene–environment association (GEA) analyses in order to examine drivers of population structure. Analyses comprised >2,600 fish from 72 riverine populations spanning a central part of the species' distribution in northern Europe. We report hitherto unidentified genetic breaks in population structure, indicating strong barriers to gene flow. GEA loci were widely spread across genomic regions and showed correlations with climatic, abiotic and geographical parameters. In some cases, individual loci showed consistent GEA across the geographical regions Britain, Europe and Scandinavia. In other cases, correlations were observed only within a sub‐set of regions, suggesting that locus‐specific variation was associated with local processes. A paired‐population sampling design allowed us to evaluate sampling effects on detection of outlier loci and GEA. Two widely applied methods for outlier detection (pcadapt and bayescan) showed low overlap in loci identified as statistical outliers across sub‐sets of data. Two GEA analytical approaches (LFMM and RDA) showed good correspondence concerning loci associated with specific variables, but LFMM identified five times more statistically significant associations than RDA. Our results emphasize the importance of carefully considering the statistical methods applied for the hypotheses being tested in outlier analysis. Sampling design may have lower impact on results if the objective is to identify GEA loci and their population distribution. Our study provides new insights into trout populations, and results have direct management implications in serving as a tool for identification of conservation units.

Sandeels are an ecologically important group of fishes; they are a key part of the food chain serving as food for marine mammals, seabirds and fish. Sandeels are further targeted by a large industrial fishery, which has led to concern about ecosystem effects. In the North Sea, the lesser sandeel Ammodytes marinus is by far the most prevalent species of sandeel in the fishery. Management of sandeel in the North Sea plus the Kattegat is currently divided into seven geographical areas, based on subtle differences in demography, population dynamics and results from simulations of larval dispersal. However, little is known about the underlying genetic population structure. In this study, we used 2,522 SNPs derived from restriction site‐associated DNA sequencing (RADseq) typed in 429 fish representing four main sandeel management areas. Our main results showed (a) a lack of a clear spatially defined genetic structure across the majority of genetic markers and (b) the existence of a group of at least 13 SNPs under strong linkage disequilibrium which together separate North Sea sandeel into three haplotype clusters, suggestive of one or more structural variants in the genome. Analyses of the spatial distribution of these putative structural variants suggest at least partial reproductive isolation of sandeel in the western management area along the Scottish coast, supporting a separate management. Our results highlight the importance of the application of a large number of markers to be able to detect weak patterns of differentiation. This study contributes to increasing the genetic knowledge of this important exploited species, and results can be used to improve our understanding of population dynamics and stock structure.

The world is changing at a rapid rate, threatening extinction for a large part of the world's biota. One potential response to those altered conditions is to evolve so as to be able to persist in place. Such evolution includes not just traits themselves, but also the phenotypic plasticity of those traits. We used individual‐based simulations to explore the potential of an evolving phenotypic plasticity to increase the probability of persistence in the response to either a step change or continual, directional change in the environment accompanied by within‐generation random environmental fluctuations. Populations could evolve by altering both their nonplastic and plastic genetic components. We found that phenotypic plasticity enhanced survival and adaptation if that plasticity was not costly. If plasticity was costly, for it to be beneficial the phenotypic magnitude of plasticity had to be great enough in the initial generations to overcome those costs. These results were not sensitive to either the magnitude of the within‐generation correlation between the environment of development and the environment of selection or the magnitude of the environmental fluctuations, except for very small phenotypic magnitudes of plasticity. So, phenotypic plasticity has the potential to enhance survival; however, more data are needed on the ubiquity of trait plasticity, the extent of costs of plasticity, and the rate of mutational input of genetic variation for plasticity.

[This corrects the article DOI: 10.1111/eva.12607.].

Ongoing speciation in the most important African malaria vectors gives rise to cryptic populations, which differ remarkably in their behavior, ecology, and capacity to vector malaria parasites. Understanding the population structure and the drivers of genetic differentiation among mosquitoes is crucial for effective disease control because heterogeneity within vector species contributes to variability in malaria cases and allow fractions of populations to escape control efforts. To examine population structure and the potential impacts of recent large-scale control interventions, we have investigated the genomic patterns of differentiation in mosquitoes belonging to the Anopheles nili group-a large taxonomic group that diverged ~3 Myr ago. Using 4,343 single nucleotide polymorphisms (SNPs), we detected strong population structure characterized by high-FST values between multiple divergent populations adapted to different habitats within the Central African rainforest. Delineating the cryptic species within the Anopheles nili group is challenging due to incongruence between morphology, ribosomal DNA, and SNP markers consistent with incomplete lineage sorting and/or interspecific gene flow. A very high proportion of loci are fixed (FST = 1) within the genome of putative species, which suggests that ecological and/or reproductive barriers are maintained by strong selection on a substantial number of genes.

Technological and methodological advances have facilitated the use of genetic data to infer census population size (Nc) in natural populations, particularly where traditional mark-and-recapture is challenging. The effective number of breeders (Nb) describes how many adults effectively contribute to a cohort and is often correlated with Nc. Predicting Nc from Nb or vice versa in species with overlapping generations has important implications for conservation by permitting (i) estimation of the more difficult to quantify variable and (ii) inferences of Nb/Nc relationships in related species lacking data. We quantitatively synthesized Nb/Nc relationships in three salmonid fishes where sufficient data have recently accumulated. Mixed-effects models were analysed in which each variable was included as a dependent variable or predictor term (Nb from Nc and vice versa). Species-dependent Nb/Nc slope estimates were significantly positive in two of three species. Variation in species slopes was likely due to varying life histories and reinforce caution when inferring Nb/Nc from taxonomically related species. Models provided maximum probable estimates for Nb and Nc for two species. However, study, population and year effects explained substantial amounts of variation (39%-57%). Consequently, prediction intervals were wide and included or were close to zero for all population sizes and species; model predictive utility was limited. Cost-benefit trade-offs when estimating Nb and/or Nc were also discussed using a real-world system example. Our findings based on salmonids suggest that no short cuts currently exist when estimating population size and researchers should focus on quantifying the variable of interest or be aware of caveats when inferring the desired variable because of cost or logistics. We caution that the salmonid species examined share life-history traits that may obscure relationships between Nb and Nc. Sufficient data on other taxa were unavailable; additional research examining Nb/Nc relationships in species with potentially relevant life-history trait differences (e.g., differing survival curves) is needed.

Limestone karsts in southern China are characterized by high edaphic and topographic heterogeneity and host high levels of species richness and endemism. However, the evolutionary mechanisms for generating such biodiversity remain poorly understood. Here, we performed species delimitation, population genetic analyses, simulations of gene flow and analyses of floral morphological traits to infer the geographic history of speciation in a species complex of Primulina eburnea from limestone karsts of southern China. Using Bayesian species delimitation, we determined that there are seven distinct species that correspond well to the putative morphological species. Species tree reconstruction, Structure and Neighbor-Net analyses all recovered four lineages in agreement with currently species geographic boundaries. High levels of genetic differentiation were observed both within and among species. Isolation-migration coalescent analysis provides evidence for significant but low gene flow among species. Approximate Bayesian computation (ABC) analysis supports a scenario of historical gene flow rather than recent contemporary gene flow for most species divergences. Finally, we found no evidence of divergent selection contributing to population differentiation of a suite of flower traits. These results support the prevalence of allopatric speciation and highlight the role of geographic isolation in the diversification process. At small geographic scales, limited hybridization occurred in the past between proximate populations but did not eliminate species boundaries. We conclude that limited gene flow might have been the predominant evolutionary force in promoting population differentiation and speciation.

Insecticide resistance is an escalating global issue for a wide variety of agriculturally important pests. The genetic basis and biochemical mechanisms of resistance are well characterized in some systems, but little is known about the ecological aspects of insecticide resistance. We therefore designed a laboratory experiment to quantify the progression of phosphine resistance in Tribolium castaneum populations subject to different immigration regimes and selection pressures. Mated resistant females were added to originally susceptible populations under two distinct migration rates, and in addition, half of the populations in each migration treatment were exposed to selection pressures from phosphine fumigation. The progression of phosphine resistance was assessed by screening beetles for the resistance allele at rph2. Phosphine resistance increased slowly in the low migration treatment and in the absence of selection, as expected. But at the higher migration rate, the increase in frequency of the resistance allele was lower than predicted. These outcomes result from the high levels of polyandry known in T. castaneum females in the laboratory, because most of the Generation 1 offspring (86%) were heterozygous for the rph2 allele, probably because resistant immigrant females mated again on arrival. Phosphine resistance was not fixed by fumigation as predicted, perhaps because susceptible gametes and eggs survived fumigation within resistant females. In terms of phosphine resistance progression in populations exposed to selection, the effect of fumigation negated the difference in migration rates. These results demonstrate how species-specific traits relating to the mating system may shape the progression of insecticide resistance within populations, and they have broad implications for the management of phosphine resistance in T. castaneum in the field. We specify and discuss how these mating system attributes need to be accounted for when developing guidelines for resistance management.