Adaptive population differentiation in phenology across a latitudinal gradient in European aspen (Populus tremula, L.): a comparison of neutral markers, candidate genes and phenotypic traits

A correct timing of growth cessation and dormancy induction represents a critical ecological and evolutionary trade-off between survival and growth in most forest trees (Rehfeldt et al. 1999; Horvath et al. 2003; Howe et al. 2003). We have studied the deciduous tree European Aspen (Populus tremula) across a latitudinal gradient and compared genetic differentiation in phenology traits with molecular markers. Trees from 12 different areas covering 10 latitudinal degrees were cloned and planted in two common gardens. Several phenology traits showed strong genetic differentiation and clinal variation across the latitudinal gradient, with Q(ST) values generally exceeding 0.5. This is in stark contrast to genetic differentiation at several classes of genetic markers (18 neutral SSRs, 7 SSRs located close to phenology candidate genes and 50 SNPs from five phenology candidate genes) that all showed F(ST) values around 0.015. We thus find strong evidence for adaptive divergence in phenology traits across the latitudinal gradient. However, the strong population structure seen at the quantitative traits is not reflected in underlying candidate genes. This result fit theoretical expectations that suggest that genetic differentiation at candidate loci is better described by F(ST) at neutral loci rather than by Q(ST) at the quantitative traits themselves.     

Disentangling complex genomic signals to understand population structure of an exploited,estuarine‐dependent flatfish

Interpreting contemporary patterns of population structure requires an understanding of the interactions among microevolutionary forces and past demographic events. Here, 4,122 SNP‐containing loci were used to assess structure in southern flounder (Paralichthys lethostigma) sampled across its range in the US Atlantic Ocean (Atlantic) and Gulf of Mexico (Gulf) and relationships among components of genomic variation and spatial and environmental variables were assessed across estuarine population samples in the Gulf. While hierarchical amova revealed significant heterogeneity within and between the Atlantic and Gulf, pairwise comparisons between samples within ocean basins demonstrated that all significant heterogeneity occurred within the Gulf. The distribution of Tajima''s D estimated at a genome‐wide scale differed significantly from equilibrium in all estuaries, with more negative values occurring in the Gulf. Components of genomic variation were significantly associated with environmental variables describing individual estuaries, and environment explained a larger component of variation than spatial proximity. Overall, results suggest that there is genetic spatial autocorrelation caused by shared larval sources for proximal nurseries (migration/drift), but that it is modified by environmentally driven differentiation (selection). This leads to conflicting signals in different parts of the genome and creates patterns of divergence that do not correspond to paradigms of strong local directional selection.     

Influence of environmental conditions at spawning sites and migration routes on adaptive variation and population connectivity in Chinook salmon

Many species that undergo long breeding migrations, such as anadromous fishes, face highly heterogeneous environments along their migration corridors and at their spawning sites. These environmental challenges encountered at different life stages may act as strong selective pressures and drive local adaptation. However, the relative influence of environmental conditions along the migration corridor compared with the conditions at spawning sites on driving selection is still unknown. In this study, we performed genome–environment associations (GEA) to understand the relationship between landscape and environmental conditions driving selection in seven populations of the anadromous Chinook salmon (Oncorhynchus tshawytscha)—a species of important economic, social, cultural, and ecological value—in the Columbia River basin. We extracted environmental variables for the shared migration corridors and at distinct spawning sites for each population, and used a Pool‐seq approach to perform whole genome resequencing. Bayesian and univariate GEA tests with migration‐specific and spawning site‐specific environmental variables indicated many more candidate SNPs associated with environmental conditions at the migration corridor compared with spawning sites. Specifically, temperature, precipitation, terrain roughness, and elevation variables of the migration corridor were the most significant drivers of environmental selection. Additional analyses of neutral loci revealed two distinct clusters representing populations from different geographic regions of the drainage that also exhibit differences in adult migration timing (summer vs. fall). Tests for genomic regions under selection revealed a strong peak on chromosome 28, corresponding to the GREB1L/ROCK1 region that has been identified previously in salmonids as a region associated with adult migration timing. Our results show that environmental variation experienced throughout migration corridors imposed a greater selective pressure on Chinook salmon than environmental conditions at spawning sites.     

Identifying Signatures of Selection in Genetic Time Series

Both genetic drift and natural selection cause the frequencies of alleles in a population to vary over time. Discriminating between these two evolutionary forces, based on a time series of samples from a population, remains an outstanding problem with increasing relevance to modern data sets. Even in the idealized situation when the sampled locus is independent of all other loci, this problem is difficult to solve, especially when the size of the population from which the samples are drawn is unknown. A standard χ²-based likelihood-ratio test was previously proposed to address this problem. Here we show that the χ²-test of selection substantially underestimates the probability of type I error, leading to more false positives than indicated by its P-value, especially at stringent P-values. We introduce two methods to correct this bias. The empirical likelihood-ratio test (ELRT) rejects neutrality when the likelihood-ratio statistic falls in the tail of the empirical distribution obtained under the most likely neutral population size. The frequency increment test (FIT) rejects neutrality if the distribution of normalized allele-frequency increments exhibits a mean that deviates significantly from zero. We characterize the statistical power of these two tests for selection, and we apply them to three experimental data sets. We demonstrate that both ELRT and FIT have power to detect selection in practical parameter regimes, such as those encountered in microbial evolution experiments. Our analysis applies to a single diallelic locus, assumed independent of all other loci, which is most relevant to full-genome selection scans in sexual organisms, and also to evolution experiments in asexual organisms as long as clonal interference is weak. Different techniques will be required to detect selection in time series of cosegregating linked loci.     

Integrated Environmental and Genomic Analysis Reveals the Drivers of Local Adaptation in African Indigenous Chickens

Breeding for climate resilience is currently an important goal for sustainable livestock production. Local adaptations exhibited by indigenous livestock allow investigating the genetic control of this resilience. Ecological niche modeling (ENM) provides a powerful avenue to identify the main environmental drivers of selection. Here, we applied an integrative approach combining ENM with genome-wide selection signature analyses (XPEHH and Fst) and genotype−environment association (redundancy analysis), with the aim of identifying the genomic signatures of adaptation in African village chickens. By dissecting 34 agro-climatic variables from the ecosystems of 25 Ethiopian village chicken populations, ENM identified six key drivers of environmental challenges: One temperature variable—strongly correlated with elevation, three precipitation variables as proxies for water availability, and two soil/land cover variables as proxies of food availability for foraging chickens. Genome analyses based on whole-genome sequencing (n = 245), identified a few strongly supported genomic regions under selection for environmental challenges related to altitude, temperature, water scarcity, and food availability. These regions harbor several gene clusters including regulatory genes, suggesting a predominantly oligogenic control of environmental adaptation. Few candidate genes detected in relation to heat-stress, indicates likely epigenetic regulation of thermo-tolerance for a domestic species originating from a tropical Asian wild ancestor. These results provide possible explanations for the rapid past adaptation of chickens to diverse African agro-ecologies, while also representing new landmarks for sustainable breeding improvement for climate resilience. We show that the pre-identification of key environmental drivers, followed by genomic investigation, provides a powerful new approach for elucidating adaptation in domestic animals.     

Museomics Dissects the Genetic Basis for Adaptive Seasonal Coloration in the Least Weasel

Dissecting the link between genetic variation and adaptive phenotypes provides outstanding opportunities to understand fundamental evolutionary processes. Here, we use a museomics approach to investigate the genetic basis and evolution of winter coat coloration morphs in least weasels (Mustela nivalis), a repeated adaptation for camouflage in mammals with seasonal pelage color moults across regions with varying winter snow. Whole-genome sequence data were obtained from biological collections and mapped onto a newly assembled reference genome for the species. Sampling represented two replicate transition zones between nivalis and vulgaris coloration morphs in Europe, which typically develop white or brown winter coats, respectively. Population analyses showed that the morph distribution across transition zones is not a by-product of historical structure. Association scans linked a 200-kb genomic region to coloration morph, which was validated by genotyping museum specimens from intermorph experimental crosses. Genotyping the wild populations narrowed down the association to pigmentation gene MC1R and pinpointed a candidate amino acid change cosegregating with coloration morph. This polymorphism replaces an ancestral leucine residue by lysine at the start of the first extracellular loop of the protein in the vulgaris morph. A selective sweep signature overlapped the association region in vulgaris, suggesting that past adaptation favored winter-brown morphs and can anchor future adaptive responses to decreasing winter snow. Using biological collections as valuable resources to study natural adaptations, our study showed a new evolutionary route generating winter color variation in mammals and that seasonal camouflage can be modulated by changes at single key genes.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

A practical guide to environmental association analysis in landscape genomics

Landscape genomics is an emerging research field that aims to identify the environmental factors that shape adaptive genetic variation and the gene variants that drive local adaptation. Its development has been facilitated by next‐generation sequencing, which allows for screening thousands to millions of single nucleotide polymorphisms in many individuals and populations at reasonable costs. In parallel, data sets describing environmental factors have greatly improved and increasingly become publicly accessible. Accordingly, numerous analytical methods for environmental association studies have been developed. Environmental association analysis identifies genetic variants associated with particular environmental factors and has the potential to uncover adaptive patterns that are not discovered by traditional tests for the detection of outlier loci based on population genetic differentiation. We review methods for conducting environmental association analysis including categorical tests, logistic regressions, matrix correlations, general linear models and mixed effects models. We discuss the advantages and disadvantages of different approaches, provide a list of dedicated software packages and their specific properties, and stress the importance of incorporating neutral genetic structure in the analysis. We also touch on additional important aspects such as sampling design, environmental data preparation, pooled and reduced‐representation sequencing, candidate‐gene approaches, linearity of allele–environment associations and the combination of environmental association analyses with traditional outlier detection tests. We conclude by summarizing expected future directions in the field, such as the extension of statistical approaches, environmental association analysis for ecological gene annotation, and the need for replication and post hoc validation studies.     

Experimental demonstration of a causal relationship between heterogeneity of selection and genetic differentiation in quantitative traits

Comparisons of estimates of genetic differentiation at molecular markers (F(ST)) and at quantitative traits (Q(ST)) are a means of inferring the level and heterogeneity of selection in natural populations. However, such comparisons are questionable because they require that the influence of drift and selection on Q(ST) be detectable over possible background influences of environmental or nonadditive genetic effects on Q(ST)-values. Here we test this using an experimental evolution approach in metapopulations of Arabidopsis thaliana experiencing different levels of drift and selection heterogeneity. We estimated the intensity and heterogeneity of selection on morphological and phenological traits via selection differentials. We demonstrate that Q(ST)-values increased with increasing selection heterogeneity when genetic drift was limited. The effect of selection on Q(ST) was thus detectable despite significant genotype-by-environment interactions that most probably biased the estimates of genetic differentiation. Although they cannot be used as a direct validation of the conclusions of prior studies, our results strongly support both the relevance of Q(ST) as an estimator of genetic differentiation and the role of local selection in shaping the genetic differentiation of natural populations.     

Natural selection and the genetic differentiation of coastal and Arctic populations of the Atlantic cod in northern Norway: a test involving nucleotide sequence variation at the pantophysin (PanI) locus

To examine the role of contemporary selection in maintaining significant allele frequency differences at the pantophysin (PanI) locus among populations of the Atlantic cod, Gadus morhua, in northern Norway, we sequenced 127 PanIA alleles sampled from six coastal and two Barents Sea populations. The distributions of variable sites segregating within the PanIA allelic class were then compared among the populations. Significant differences were detected in the overall frequencies of PanIA alleles among populations within coastal and Arctic regions that was similar in magnitude to heterogeneity in the distributions of polymorphic sites segregating within the PanIA allelic class. The differentiation observed at silent sites in the PanIA allelic class contradicts the predicted effects of widescale gene flow and suggests that postsettlement selection acting on cohorts cannot be responsible for the genetic differences described between coastal and Arctic populations. Our results suggest that the marked differences observed between coastal and Arctic populations of G. morhua in northern Norway at the PanI locus reflect the action of recent diversifying selection and that populations throughout the region may be more independent than suggested by previous studies.     

The Effects of Background and Interference Selection on Patterns of Genetic Variation in Subdivided Populations

It is well known that most new mutations that affect fitness exert deleterious effects and that natural populations are often composed of subpopulations (demes) connected by gene flow. To gain a better understanding of the joint effects of purifying selection and population structure, we focus on a scenario where an ancestral population splits into multiple demes and study neutral diversity patterns in regions linked to selected sites. In the background selection regime of strong selection, we first derive analytic equations for pairwise coalescent times and F_ST as a function of time after the ancestral population splits into two demes and then construct a flexible coalescent simulator that can generate samples under complex models such as those involving multiple demes or nonconservative migration. We have carried out extensive forward simulations to show that the new methods can accurately predict diversity patterns both in the nonequilibrium phase following the split of the ancestral population and in the equilibrium between mutation, migration, drift, and selection. In the interference selection regime of many tightly linked selected sites, forward simulations provide evidence that neutral diversity patterns obtained from both the nonequilibrium and equilibrium phases may be virtually indistinguishable for models that have identical variance in fitness, but are nonetheless different with respect to the number of selected sites and the strength of purifying selection. This equivalence in neutral diversity patterns suggests that data collected from subdivided populations may have limited power for differentiating among the selective pressures to which closely linked selected sites are subject.     

Adaptive phenotypic and genomic divergence in the common chaffinch (Fringilla coelebs) following niche expansion within a small oceanic island

According to models of ecological speciation, adaptation to adjacent, contrasting habitat types can lead to population divergence given strong enough environment-driven selection to counteract the homogenizing effect of gene flow. We tested this hypothesis in the common chaffinch (Fringilla coelebs) on the small island of La Palma, Canary Islands, where it occupies two markedly different habitats. Isotopic (δ¹³C, δ¹⁵N) analysis of feathers indicated that birds in the two habitats differed in ecosystem and/or diet, and analysis of phenotypic traits revealed significant differences in morphology and plumage colouration that are consistent with ecomorphological and ecogeographical predictions respectively. A genome-wide survey of single-nucleotide polymorphism revealed marked neutral structure that was consistent with geography and isolation by distance, suggesting low dispersal. In contrast, loci putatively under selection identified through genome-wide association and genotype-environment association analyses, revealed a marked adaptive divergence between birds in both habitats. Loci associated with phenotypic and environmental differences among habitats were distributed across the genome, as expected for polygenic traits involved in local adaptation. Our results suggest a strong role for habitat-driven local adaptation in population divergence in the chaffinches of La Palma, a process that appears to be facilitated by a strong reduction in effective dispersal distances despite the birds' high dispersal capacity.     

Effects of forest disturbance on the fitness of an endemic rodent in a biodiversity hotspot

Praomys delectorum occurs abundantly in both disturbed and intact forests in the Ukaguru Mountains within the Eastern Arc Mountains (EAM), Morogoro, Tanzania. While previous studies have reported that anthropogenic disturbances such as grazing, wood cutting, and harvesting have a positive effect on the population density of P. delectorum, the impact of habitat disturbance on its demographic traits is still unknown. We performed a capture–mark–recapture study in both disturbed and intact forests from June 2018 to February 2020 in order to investigate the effects of habitat disturbance on abundance and two demographic traits: survival and maturation of P. delectorum in the Ukaguru Mountains. We found no variation in abundance or maturation between intact and disturbed forests, but habitat type did affect survival. However, this effect was sex‐dependent since female survival was higher in disturbed forests, while male survival remained similar across the two forest types potentially due to differences in predation pressure or food availability between the two habitats. Continuous demographic monitoring of P. delectorum in EAM is necessary given that the increasing human population surrounding the landscape is leading to higher deforestation rates and expansion of the pine plantation in the forest reserve.     

Divergent selection and heterogeneous migration rates across the range of Sitka spruce (Picea sitchensis)

Gene flow and effective population size (N(e)) should depend on a population's position within its range: those near the edges are expected to have smaller N(e) and lower relative emigration rates, whereas those nearer the centre should have larger N(e) and higher relative emigration rates. In species with continuous ranges, this phenomenon may limit the ability of peripheral populations to respond to divergent selection. Here, we employ Sitka spruce as a model to test these predictions. We previously genotyped 339 single nucleotide polymorphisms (SNPs) in 410 individuals from 13 populations, and used these data to identify putative targets of divergent selection, as well as to explore the extent to which central-peripheral structure may impede adaptation. Fourteen SNPs had outlier F(ST) estimates suggestive of divergent selection, of which nine were previously associated with phenotypic variation in adaptive traits (timing of autumn budset and cold hardiness). Using coalescent simulations, we show that populations from near the centre of the range have higher effective populations sizes than those from the edges, and that central populations contribute more migrants to marginal populations than the reverse. Our results suggest that while divergent selection appears to have shaped allele frequencies among populations, asymmetrical movement of alleles from the centre to the edges of the species range may affect the adaptive capacity of peripheral populations. In southern peripheral populations, the movement of cold-adapted alleles from the north represents a significant impediment to adaptation under climate change, while in the north, movement of warm-adapted alleles from the south may enhance adaptation.     

Evolutionary dynamics of quantitative variation in an adaptive trait at the regional scale: The case of zinc hyperaccumulation in Arabidopsis halleri

Metal hyperaccumulation in plants is an ecological trait whose biological significance remains debated, in particular because the selective pressures that govern its evolutionary dynamics are complex. One of the possible causes of quantitative variation in hyperaccumulation may be local adaptation to metalliferous soils. Here, we explored the population genetic structure of Arabidopsis halleri at fourteen metalliferous and nonmetalliferous sampling sites in southern Poland. The results were integrated with a quantitative assessment of variation in zinc hyperaccumulation to trace local adaptation. We identified a clear hierarchical structure with two distinct genetic groups at the upper level of clustering. Interestingly, these groups corresponded to different geographic subregions, rather than to ecological types (i.e., metallicolous vs. nonmetallicolous). Also, approximate Bayesian computation analyses suggested that the current distribution of A. halleri in southern Poland could be relictual as a result of habitat fragmentation caused by climatic shifts during the Holocene, rather than due to recent colonization of industrially polluted sites. In addition, we find evidence that some nonmetallicolous lowland populations may have actually derived from metallicolous populations. Meanwhile, the distribution of quantitative variation in zinc hyperaccumulation did separate metallicolous and nonmetallicolous accessions, indicating more recent adaptive evolution and diversifying selection between metalliferous and nonmetalliferous habitats. This suggests that zinc hyperaccumulation evolves both ways—towards higher levels at nonmetalliferous sites and lower levels at metalliferous sites. Our results open a new perspective on possible evolutionary relationships between A. halleri edaphic types that may inspire future genetic studies of quantitative variation in metal hyperaccumulation.     

3种高山植物的物候和种群分布格局在融雪梯度上的变化

在青藏高原东部的一个高山雪床,沿着融雪梯度分别设置早融、中间和晚融3个融雪部位,然后测定川西小黄菊（Pyrethrum tatsienense）、长叶火绒草（Leontopodium longifolium）和圆穗蓼（Polygonum macrophyllum）在3个融雪部位上的物候差异以及种群分布格局的变化。结果表明：从早融到晚融的梯度上,3个物种的物候期都不同程度地有所推迟。其中,开始生长的时间推迟12~14 d,始花期推迟6~8 d,盛花期推迟6 d左右,但同一种植物在不同的融雪部位上的衰老枯黄期趋于一致,这标志着在晚融部位同一植物的生长期要缩短。在种群层次上,长叶火绒草和圆穗蓼的分布格局随着融雪的推迟都发生了一定的变化,基本上表现为从早融部位的集群分布到中间或晚融部位的随机分布。川西小黄菊在各个融雪部位上都表现为集群分布,但集群的强度随融雪的推迟逐渐减弱。     

Population Structure and Growth Patterns of Opuntia echios var. gigantea along an Elevational Gradient in the Galípagos Islands1

Cacti growing in forests potentially experience growth limitation due to reduced light availability. To test this hypothesis, we studied the population structure of Opuntia echios var. gigantea at 15 sites on the south side of Isla Santa Cruz, Galípagos Islands, Ecuador. Populations were located in communities ranging from arid scrub at low elevations to closed‐canopy tropical dry forest at higher elevations. Ordination confirmed the existence of a strong elevation‐vegetation gradient. Opuntia abundance peaked at lower elevations (ca 30 m), with lower densities in closed‐canopy sites. For populations in scrub vegetation, density declined fairly regularly with plant height. Populations in forested sites had few plants of intermediate height, suggesting periodic recruitment. Scrub populations had random dispersion, while those in forests were aggregated. The change in spatial pattern may be related to a change in primary reproductive mode from asexual propagation via fallen fruits to propagation via fallen cladodes. Height was significantly correlated with stem diameter. Intercepts of these relationships increased toward higher elevations, probably in response to the increasing height of the surrounding canopy.     

Geographic variation in adaptation at the molecular level: a case study of plant immunity genes

Natural selection imposed by interacting species frequently varies among geographic locations and can lead to local adaptation, where alternative phenotypes are found in different populations. Little is known, however, about whether geographically variable selection acting on traits that mediate species interactions is consistent or strong enough to influence patterns of nucleotide variation at individual loci. To investigate this question, we examined patterns of nucleotide diversity and population structure at 16 plant innate immunity genes, with putative functions in defending plants against pathogens or herbivores, from six populations of teosinte (Zea mays ssp. parviglumis). Specifically, we tested whether patterns of population structure and within-population diversity at immunity genes differed from patterns found at nonimmunity (reference) loci and from neutral expectations derived from coalescent simulations of structured populations. For the majority of genes, we detected no strong evidence of geographically variable selection. However, in the wound-induced serine protease inhibitor (wip1), which inhibits the hydrolysis of dietary proteins in insect herbivores, one population showed unusually high levels of genetic differentiation, very low levels of nucleotide polymorphism, and was fixed for a novel replacement substitution in the active site of the protein. Taken together, these data suggest that wip1 experienced a recent selective sweep in one geographic region; this pattern may reflect local adaptation or an ongoing species-wide sweep. Overall, our results indicate that a signature of local adaptation at the molecular level may be uncommon-particularly for traits that are under complex genetic control.     

Toward a conceptual framework for managing and conserving marine habitats: A case study of kelp forests in the Salish Sea

Kelp forests are in decline across much of their range due to place‐specific combinations of local and global stressors. Declines in kelp abundance can lead to cascading losses of biodiversity and productivity with far‐reaching ecological and socioeconomic consequences. The Salish Sea is a hotspot of kelp diversity where many species of kelp provide critical habitat and food for commercially, ecologically, and culturally important fish and invertebrate species. However, like other regions, kelp forests in much of the Salish Sea are in rapid decline. Data gaps and limited long‐term monitoring have hampered attempts to identify and manage for specific drivers of decline, despite the documented urgency to protect these important habitats. To address these knowledge gaps, we gathered a focus group of experts on kelp in the Salish Sea to identify perceived direct and indirect stressors facing kelp forests. We then conducted a comprehensive literature review of peer‐reviewed studies from the Salish Sea and temperate coastal ecosystems worldwide to assess the level of support for the pathways identified by the experts, and we identified knowledge gaps to prioritize future research. Our results revealed major research gaps within the Salish Sea and highlighted the potential to use expert knowledge for making informed decisions in the region. We found high support for the pathways in the global literature, with variable consensus on the relationship between stressors and responses across studies, confirming the influence of local ecological, oceanographic, and anthropogenic contexts and threshold effects on stressor–response relationships. Finally, we prioritized areas for future research in the Salish Sea. This study demonstrates the value expert opinion has to inform management decisions. These methods are readily adaptable to other ecosystem management contexts, and the results of this case study can be immediately applied to kelp management.     

Response of Deer Mice (Peromyscus maniculatus) to shrubs in shortgrass prairie: linking small-scale movements and the spatial distribution of individuals

1. Organisms respond to the abundance and spatial distribution of resources at both individual and population scales but there have been relatively few attempts to link insights from studies of these different phenomena, especially for wide-ranging vertebrates.
2. Deer Mice ( Peromyscus maniculatus ) were live-trapped and tracked across a gradient of shrub cover in shortgrass prairie to determine patterns of abundance, microhabitat use and movements.
3. In areas with few shrubs, mice preferred shrub microhabitats and their movement trails were relatively straight. Both trail tortuosity and population density increased with increasing shrub cover, so that mice tended to accumulate in areas where their trails were most convoluted. However, movements were also linear where shrubs were dense, presumably because mice could achieve the benefits of association with shrubs without travelling directly beneath them. Areas with dense shrubs also had high but variable population densities, suggesting that other factors such as intraspecific interactions may have affected movements on a larger, home-range scale.
4. Apparent thresholds in the selective vs random use of shrubs, movement patterns and abundance occurred over a narrow range of shrub cover where shrubs were most aggregated, underscoring the importance of both shrub density and dispersion. Non-linear relationships in the response to resources at both behavioural and population scales thus may complicate our ability to predict abundance from individual movements across a broad range of resource distributions.