Search

Protein isoforms (PIs) play pivotal roles in regulating plant growth and development that confer adaptability to diverse environmental conditions. PIs are widely present in plants and generated through alternative splicing (AS), alternative polyadenylation (APA), alternative initiation (AI), and ribosomal frameshifting (RF) events. The widespread presence of PIs not only significantly increases the complexity of genomic information but also greatly enriches regulatory networks and enhances their flexibility. PIs may also play important roles in phenotypic diversity, ecological niche differentiation, and speciation, thereby increasing the dimensions of research in molecular ecology. However, PIs pose new challenges for the quantitative analysis, annotation, and identification of genetic regulatory mechanisms. Thus, focus on PIs make genomic and epigenomic studies both more powerful and more challenging. This review summarizes the origins, functions, regulatory patterns of isoforms, and the challenges they present for future research in molecular ecology and molecular biology., (© 2024 The Author(s). New Phytologist © 2024 New Phytologist Foundation.)

Virus spillovers from managed honey bees, Apis mellifera, are thought to contribute to the decline of wild pollinators, including bumble bees. However, data on the impact of such viruses on wild pollinators remain scarce, and the influence of landscape structure on virus dynamics is poorly understood. In this study, we deployed bumble bee colonies in an agricultural landscape and studied changes in the bumble bee virome during field placement under varying habitat composition and configuration using a multiscale analytical framework. We estimated prevalence of viruses and viral loads (i.e. number of viral genomic equivalent copies) in bumble bees before and after placing them in the field using next generation sequencing and quantitative PCR. The results show that viral loads and number of different viruses present increased during placement in the field and that the virus composition of the colonies shifted from an initial dominance of honey bee associated viruses to a higher number (in both viral loads and number of viruses present) of bumble bee associated viruses. Especially DWV-B, typical for honey bees, drastically decreased after the time in the field. Viral loads prior to placing colonies in the field showed no effect on colony development, suggesting low impacts of these viruses in field settings. Notably, we further demonstrate that increased habitat diversity results in a lower number of different viruses present in Bombus colonies, while colonies in areas with well-connected farmland patches decreased in their total viral load after field placement. Our results emphasize the importance of landscape heterogeneity and connectivity for wild pollinator health and that these influences predominate at fine spatial scales., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Protected areas are an important tool for wildlife conservation; however, research is increasingly revealing both biases and inadequacies in the global protected area network. One common criticism is that protected areas are frequently located in remote, high-elevation regions, which may face fewer threats compared to more accessible locations. To explore the conservation implications of this issue, we consider a thought experiment with seven different counterfactual scenarios for the Sunda clouded leopard's conservation on Borneo. This allows us to examine two contrasting paradigms for conservation: "proactive conservation" which prioritises areas with high biodiversity and high risk of development, and "expedient conservation" which focusses on areas with the lowest development risk. We select clouded leopards as our focal species not only because of their emerging conservation importance, but also because, as top predators, they represent both keystone species and ambassadors for wider forest biodiversity. Furthermore, a published analysis of the likely impacts of forest loss in their habitat provides a benchmark for evaluating the modelled outcomes of alternative hypothetical conservation scenarios. We find that, across all metrics, expedient reserve design offered few benefits over the business-as-usual scenario, in contrast to the much greater conservation effectiveness of proactive protected area design. This paper sheds light on the challenging trade-offs between conservation goals and the competing land uses essential for the economic development and well-being of local communities., (© 2024. The Author(s).)

The islands of Borneo and Sumatra are strongholds for biodiversity and home for many endemic species. They also have experienced amongst the highest deforestation rates globally. Both islands are undergoing massive, rapid infrastructure development, leading to further deforestation and habitat fragmentation. Here, we identify priority areas for continued functional forest connectivity across Borneo and Sumatra, using spatial models of clouded leopard (Neofelis diardi, a forest indicator species) movement, and impacts thereto from existing and future infrastructure development. We specifically measure and map the anticipated impacts on forest functional connectivity of three major infrastructure development projects (Pan Borneo Highway, Trans-Sumatran Toll Road, and the new Indonesian capital city of Nusantara). We found that core clouded leopard habitats are already highly fragmented in Sumatra, constituting only ∼13 % of the island, with potential dispersal corridors still linking some habitat fragments. In Borneo, clouded leopard core habitats cover 34 % of the island, with one large central core area and several much smaller satellite cores, which are largely unprotected (15 % protected, compared to 42 % in Sumatra). The largest negative effect on habitat connectivity was predicted for Nusantara (66 % of the total connectivity loss predicted for all three infrastructure projects), reverberating across the entirety of Borneo with the strongest effects in East Kalimantan. The Pan Borneo Highway accounted for 28 % of the total connectivity loss, affecting every province in Borneo and Brunei, with 6 % of this decrease located within protected areas. The Trans-Sumatran Toll Road had the smallest negative effect on connectivity (6 %) but only when excluding the already built segments, which, when included, produce a total negative impact similar to that of the Pan Borneo Highway., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier B.V.)

The regulation of flowering time is typically governed by transcription factors or epigenetic modifications. Transcript isoforms can play important roles in flowering regulation. Recently, transcript isoforms were discovered in the key genes, OfAP1 and OfTFL1, of the flowering regulatory network in Osmanthus fragrans. OfAP1-b generates a full-length isoform of OfAP1-b1 as well as an isoform of OfAP1-b2 that lacks the C-terminal domain. Although OfAP1-b2 does not possess an activation domain, it has a complete K domain that allows it to form heterodimers. OfAP1-b2 competes with OfAP1-b1 by binding with OfAGL24 to create non-functional and functional heterodimers. As a result, OfAP1-b1 promotes flowering while OfAP1-b2 delays flowering. OfTFL1 produces two isoforms located in different areas: OfTFL1-1 in the cytoplasm and OfTFL1-2 in the nucleus. When combined with OfFD, OfTFL1-1 does not enter the nucleus to repress AP1 expression, leading to early flowering. Conversely, when combined with OfFD, OfTFL1-2 enters the nucleus to repress AP1 expression, resulting in later flowering. Tissue-specific expression and functional conservation testing of OfAP1 and OfTFL1 support the new model's effectiveness in regulating flowering. Overall, this study provides new insights into regulating flowering time by the competition of isoforms., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Annals of Botany Company. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Large and severe wildfires, exacerbated by climate change and human behavior, are occurring more frequently in many forests across the western United States. While wildfire is a natural part of most terrestrial ecosystems, rapidly changing fire regimes have the potential to alter habitat beyond the adaptive capabilities of species. Spatial assessments of wildfire risks to species habitat may allow managers to pinpoint locations for management activities. To illustrate this, we spatially assessed wildfire risk within habitat that supports the nesting activity of the federally threatened northern spotted owl (Strix occidentalis caurina) in the California redwood coast ecoregion. To accomplish this, we built a scale-optimized ensemble nesting habitat suitability model and identified habitat with the highest wildfire hazard potential. Percent canopy cover at 100-m scale, slope at 400-m scale, and January precipitation at 800-m scale were the most influential environmental covariates for predicting northern spotted owl nesting habitat. Nearly 60% of nesting habitat was predicted to be at high or very high (>1986 index value) wildfire risks. We identified three areas in the Maple Creek Area of Humboldt County, Jackson State Demonstration Forest in Mendocino County, and Point Reyes National Seashore in Marin County, California with a high concentration of nesting habitat that are at a very high risk of experiencing high severity wildfires. We recommend these areas be targeted for future research to understand the impact of wildfire on northern spotted owl as well as management attention., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

Whole genome doubling (WGD) plays a critical role in plant evolution, yet the mechanisms underlying the maintenance of overall equilibrium following an artificial doubling event, as well as its impact on phenotype and adaptability, remain unclear. By comparing the gene expression of naturally occurring weeping forsythia diploids and colchicine-induced autotetraploids under normal growth conditions and cold stress, we identified gene expression dosage responses resulting from ploidy change. Only a small proportion of effectively expressed genes showed dosage effect, and most genes did not exhibit significant expression differences. However, the genes that showed expression dosage effect were largely random. The autotetraploids had slower overall growth rates, possibly resulting from negative gene dosage effects on zeatin synthesis and multiple metabolic delays caused by other negative dosage genes. Our comparative analysis of cold response genes in diploids and autotetraploids revealed that genes related to "response to abscisic acid" and "cold acclimation" were key factors contributing to greater cold tolerance in the autotetraploids. In particular, gene expression related to "cold acclimation" might mitigate the effects of cold stress. Taken together, our findings suggested that overall gene expression equilibrium following WGD of weeping forsythia autotetraploids was achieved through the inactivation of the majority of duplicated genes. Our research provides new insights into the mechanisms regulating expression dosage balance following polyploidization events., Competing Interests: Declaration of Competing Interest The authors declare no conflict of interest., (Copyright © 2023 Elsevier Masson SAS. All rights reserved.)

Weeping forsythia is a wide-spread shrub in China with important ornamental, medicinal and ecological values. It is widely distributed in China's warm temperate zone. In plants, WRKY transcription factors play important regulatory roles in seed germination, flower development, fruit ripening and coloring, and biotic and abiotic stress response. To date, WRKY transcription factors have not been systematically studied in weeping forsythia. In this study, we identified 79 WRKY genes in weeping forsythia and classified them according to their naming rules in Arabidopsis thaliana. Phylogenetic tree analysis showed that, except for IIe subfamily, whose clustering was inconsistent with A. thaliana clustering, other subfamily clustering groups were consistent. Cis-element analysis showed that WRKY genes related to pathogen resistance in weeping forsythia might be related to methyl jasmonate and salicylic acid-mediated signaling pathways. Combining cis-element and expression pattern analyses of WRKY genes showed that more than half of WRKY genes were involved in light-dependent development and morphogenesis in different tissues. The gene expression results showed that 13 WRKY genes were involved in drought response, most of which might be related to the abscisic acid signaling pathway, and a few of which might be regulated by MYB transcription factors. The gene expression results under cold stress showed that 17 WRKY genes were involved in low temperature response, and 9 of them had low temperature responsiveness cis-elements. Our study of WRKY family in weeping forsythia provided useful resources for molecular breeding and important clues for their functional verification., (© 2023. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

Darwin proposed that the capacity of organisms to produce more offspring that can be supported by the environment would lead to a struggle for existence, and individuals that are most fit for survival and reproduction would be selected through natural selection. Ecology is the science that studies the interaction between organisms and their environment within the context of Darwinian evolution, and an ecosystem is defined as a community of living organisms in conjunction with the nonliving components of their environment, interacting as a system. One topic that has been very much understudied and largely ignored in evolutionary biology is the overarching context of thermodynamics in controlling all biological processes and the evolution of life. Most fundamentally, organisms are self-replicating dissipative structures. Evolution is the process whereby variation in the structure of organisms have differential fitness in terms of their effectiveness at building and maintaining their structure, efficiently consuming free energy, and effectively reproducing and passing on those heritable variations, leading to change in the frequency of genetic variation and associated change in the characteristics in the population. The central process is dissipation of free energy according to the second law of thermodynamics, and evolution therefore is better conceptualized as the emergence of self-replicating dissipative structures that through natural selection become increasingly more efficient at degrading free energy. Ecosystems are linked series of dissipative structures with heat engine dynamics driven by random dissipation of energy and increasing entropy. The structure and composition of ecosystems across scales are emergent dissipative structures driven by the flow of energy and the increase in entropy. Communities and ecosystems are emergent properties of a system that has evolved to most efficiently dissipate energy and increase entropy. By focusing on the fundamental entity (energy), and the fundamental process (dissipation and disordering of energy and increasing of entropy), we are able to have a much clearer and powerful understanding of what life is, from the level of biochemistry, to evolution, to the nature of the organism itself, and to the emergent structures of ecosystems, food webs and communities.

Many protected areas worldwide increasingly resemble habitat isolates embedded in human-modified landscapes. However, establishing linkages among protected areas could significantly reduce species-loss rates. Here we present a novel method having broad applicability for assessing enhanced regional connectivity on persistence of mammal diversity. We combine theoretically-derived species relaxation rates for mammal communities with empirically-derived pathways. We assess the value of enhanced regional connectivity for two hypothetical networks of national parks in western North America: the Yellowstone-Glacier network and the Mount Rainier-North Cascades network. Linking the Yellowstone and Glacier park assemblages by eliminating barriers to movement in identified mammal dispersal pathways and by incorporating adjacent wilderness areas and known ungulate migratory routes into a protected area network would greatly enlarge available habitat. This would enhance medium to large mammal species persistence time by factor of 4.3, on average, or ~ 682 generations relative to individual parks. Similarly, linking Mount Rainier and North Cascades park assemblages would enhance mammal species persistence time by a factor of 4.3, on average, or ~305 generations relative to individual parks. Enhancing regional connectivity among western North America parks could serve as an important template for landscape-scale conservation in the 21st century., (© 2023. The Author(s).)

In the face of a growing human footprint, understanding interactions among threatened large carnivores is fundamental to effectively mitigating anthropogenic threats and managing species. Using data from a large-scale camera trap survey, we investigated the effects of environmental and anthropogenic variables on the interspecific interaction of a carnivore guild comprising of tiger, leopard and dhole in Bhutan. We demonstrate the complex effects of human settlement density on large carnivore interactions. Specifically, we demonstrate that leopard-dhole co-occupancy probability was higher in areas with higher human settlement density. The opposite was true for tiger-leopard co-occupancy probability, but it was positively affected by large prey (gaur) abundance. These findings suggest that multi-carnivore communities across land-use gradients are spatially structured and mediated also by human presence and/or the availability of natural prey. Our findings show that space-use patterns are driven by a combination of the behavioural mechanism of each species and its interactions with competing species. The duality of the effect of settlement density on species interactions suggests that the benefits of exploiting anthropogenic environments are a trade-off between ecological opportunity (food subsidies or easy prey) and the risk of escalating conflict with humans.

As ecological data and associated analyses become more widely available, synthesizing results for effective communication with stakeholders is essential. In the case of wildlife corridors, managers in human-dominated landscapes need to identify both the locations of corridors and multiple stakeholders for effective oversight. We synthesized 5 independent studies of tiger (Panthera tigris) connectivity in central India, a global priority landscape for tiger conservation, to quantify agreement on landscape permeability for tiger movement and potential movement pathways. We used the latter analysis to identify connectivity areas on which studies agreed and stakeholders associated with these areas to determine relevant participants in corridor management. Three or more of the 5 studies' resistance layers agreed in 63% of the study area. Areas in which all studies agree on resistance were of primarily low (66%, e.g., forest) and high (24%, e.g., urban) resistance. Agreement was lower in intermediate resistance areas (e.g., agriculture). Despite these differences, the studies largely agreed on areas with high levels of potential movement: >40% of high average (top 20%) current-flow pixels were also in the top 20% of current-flow agreement pixels (measured by low variation), indicating consensus connectivity areas (CCAs) as conservation priorities. Roughly 70% of the CCAs fell within village administrative boundaries, and 100% overlapped forest department management boundaries, suggesting that people live and use forests within these priority areas. Over 16% of total CCAs' area was within 1 km of linear infrastructure (437 road, 170 railway, 179 transmission line, and 339 canal crossings; 105 mines within 1 km of CCAs). In 2019, 78% of forest land diversions for infrastructure and mining in Madhya Pradesh (which comprises most of the study region) took place in districts with CCAs. Acute competition for land in this landscape with globally important wildlife corridors calls for an effective comanagement strategy involving local communities, forest departments, Appendix 1 and infrastructure planners. This article is protected by copyright. All rights reserved., (This article is protected by copyright. All rights reserved.)

Understanding the spatial structure of genetic diversity provides insights into a populations' genetic status and enables assessment of its capacity to counteract the effects of genetic drift. Such knowledge is particularly scarce for the snow leopard, a conservation flagship species of Central Asia mountains. Focusing on a snow leopard population in the Qilian mountains of Gansu Province, China, we characterised the spatial genetic patterns by incorporating spatially explicit indices of diversity and multivariate analyses, based on different inertia levels of Principal Component Analysis (PCA). We compared two datasets differing in the number of loci and individuals. We found that genetic patterns were significantly spatially structured and were characterised by a broad geographical division coupled with a fine-scale cline of differentiation. Genetic admixture was detected in two adjoining core areas characterised by higher effective population size and allelic diversity, compared to peripheral localities. The power to detect significant spatial relationships depended primarily on the number of loci, and secondarily on the number of PCA axes. Spatial patterns and indices of diversity highlighted the cryptic structure of snow leopard genetic diversity, likely driven by its ability to disperse over large distances. In combination, the species' low allelic richness and large dispersal ability result in weak genetic differentiation related to major geographical features and isolation by distance. This study illustrates how cryptic genetic patterns can be investigated and analysed at a fine spatial scale, providing insights into the spatially variable isolation effects of both geographic distance and landscape resistance., (© 2021. The Author(s), under exclusive licence to The Genetics Society.)

Several methods have been recently proposed to calculate configurational entropy, based on Boltzmann entropy. Some of these methods appear to be fully thermodynamically consistent in their application to landscape patch mosaics, but none have been shown to be fully generalizable to all kinds of landscape patterns, such as point patterns, surfaces, and patch mosaics. The goal of this paper is to evaluate if the direct application of the Boltzmann relation is fully generalizable to surfaces, point patterns, and landscape mosaics. I simulated surfaces and point patterns with a fractal neutral model to control their degree of aggregation. I used spatial permutation analysis to produce distributions of microstates and fit functions to predict the distributions of microstates and the shape of the entropy function. The results confirmed that the direct application of the Boltzmann relation is generalizable across surfaces, point patterns, and landscape mosaics, providing a useful general approach to calculating landscape entropy.

There has been a recent surge of interest in theory and methods for calculating the entropy of landscape patterns, but relatively little is known about the thermodynamic consistency of these approaches. I posit that for any of these methods to be fully thermodynamically consistent, they must meet three conditions. First, the computed entropies must lie along the theoretical distribution of entropies as a function of total edge length, which Cushman showed was a parabolic function following from the fact that there is a normal distribution of permuted edge lengths, the entropy is the logarithm of the number of microstates in a macrostate, and the logarithm of a normal distribution is a parabolic function. Second, the entropy must increase over time through the period of the random mixing simulation, following the expectation that entropy increases in a closed system. Third, at full mixing, the entropy will fluctuate randomly around the maximum theoretical value, associated with a perfectly random arrangement of the lattice. I evaluated these criteria in a test condition involving a binary, two-class landscape using the Cushman method of directly applying the Boltzmann relation (s = klogW) to permuted landscape configurations and measuring the distribution of total edge length. The results show that the Cushman method directly applying the classical Boltzmann relation is fully consistent with these criteria and therefore fully thermodynamically consistent. I suggest that this method, which is a direct application of the classical and iconic formulation of Boltzmann, has advantages given its direct interpretability, theoretical elegance, and thermodynamic consistency.

This paper presents a multivariate textual analysis of more than 1300 papers on entropy in ecology. There are six main insights that emerged. First, there is a large body of literature that has addressed some aspect of entropy in ecology, most of which has been published in the last 5-10 years. Second, the vast majority of these papers focus on species distribution, species richness, relative abundance or trophic structure and not landscape-scale patterns or processes, pe se. Third, there have been few papers addressing landscape-level questions related to entropy. Fourth, the quantitative analysis with hierarchical clustering identified a strongly nested structure among papers that addressed entropy in ecology. Fifth, there is clear differentiation of papers focused on landscape-level applications of entropy from other papers, with landscape focused papers clustered together at each level of the hierarchy in a relatively small and closely associated group. Sixth, this group of landscape-focused papers was substructured between papers that explicitly adopted entropy measures to quantify the spatial pattern of landscape mosaics, often using variations on Boltzmann entropy, versus those that utilize Shannon entropy measures from information theory, which are not generally explicit in their assessment of spatial configuration. This review provides a comprehensive, quantitative assessment of the scope, trends and relationships among a large body of literature related to entropy in ecology and for the first time puts landscape ecological research on entropy into that context.

Oceanographic features such as currents, waves, temperature and salinity, together with life history traits, control patterns and rates of gene flow and contribute to shaping the population genetic structure of marine organisms. Seascape genetics is an emerging discipline that adopts a spatially explicit approach to examine biotic and abiotic factors that drive gene flow in marine environments. In this study, we examined factors that contribute to genetic differentiation in two coastal Mediterranean gastropods whose geographical ranges overlap but which inhabit different environments. The two species differ in several life history traits and in their dispersal capabilities. Genetic differentiation was relatively low for the trochid species Gibbula divaricata (F ST  =0.059), and high for the vermetid species Dendropoma lebeche (F ST  =0.410). Salinity emerged as the most important variable explaining the genetic structure of both species; sea surface temperature was also important for G. divaricata. For the more sessile D. lebeche, the coastline was predicted to provide important pathways for stepping-stone connectivity and gene flow. Our results provide a greater understanding of the factors influencing marine population connectivity, which may be useful to guide marine conservation and management in the Mediterranean., (© 2021 John Wiley & Sons Ltd.)

Genetic mechanisms of species local adaptation are an emerging topic of great interest in evolutionary biology and molecular ecology. In this study, we compared the changes of physiological and phenotypic indexes and gene expression of four weeping forsythia populations under cold stress through a common garden experiment. Physiological and phenotypic results showed that there were differences in cold tolerance among populations. cold tolerance of high the latitude population (HBWZ) was the strongest, followed by the middle latitude population (SXWL), while the low latitude populations (SXHM) and (SXLJ) expressed the weakest cold tolerance. We identified significant differences in gene expression of cold tolerance related pathways and ontologies, including genes of oxylipin and isoquinoline alkaloid biosynthetic process, galactose, tyrosine and unsaturated fatty acids metabolism, among these populations under the same experimental temperature treatments. Even under the same degree of stress, there were notable differences in gene expression among natural populations. In this study, we present a working model of weeping forsythia populations which evolved in the context of different intensities of cold stress. Our study provides new insights for comprehending the genetic mechanisms of local adaptation for non-model species., (Copyright © 2021 Elsevier Masson SAS. All rights reserved.)

Background: Understanding the genetic mechanisms of local adaptation is an important emerging topic in molecular ecology and evolutionary biology., Results: Here, we identify the physiological changes and differential expression of genes among different weeping forsythia populations under drought stress in common garden experiments. Physiological results showed that HBWZ might have higher drought tolerance among four populations. RNA-seq results showed that significant differential expression in the genes responding to the synthesis of flavonoids, aromatic substances, aromatic amino acids, oxidation-reduction process, and transmembrane transport occured among four populations. By further reanalysis of results of previous studies, sequence differentiation was found in the genes related to the synthesis of aromatic substances among different weeping forsythia populations., Conclusions: Overall, our study supports the hypothesis that the dual differentiation in gene efficiency and expression increases among populations in response to heterogeneous environments and is an important evolutionary process of local adaptation. Here, we proposed a new working model of local adaptation of weeping forsythia populations under different intensities of drought stress, which provides new insights for understanding the genetic mechanisms of local adaptation for non-model species.

Southern Iran is a conservation priority area for the endangered Persian onager (Equus hemionus onager), which is threatened by habitat fragmentation and conflict with local communities. To better understand factors that influence onager conservation, we administered a questionnaire in local communities to survey their ecological knowledge, personal experience related to onager, and attitudes toward traditional solutions for reducing crop damage by onager. In addition, we used resistant kernel and factorial least-cost path analyses to identify core areas and corridors for onager movement, and spatial randomization of vehicle collisions and crossing locations to test the predictive ability of resistant kernel and factorial least-cost path predictions of movement. We found that local communities that were knowledgeable about onagers experienced less crop damage from onager compared with those who used traditional methods. Habitat connectivity models revealed that core areas of movement are highly concentrated at the center of protected areas. Some sections of core areas have been cut off by roads where most vehicle collisions with onagers occurred. We propose that effective onager conservation will require integrated landscape-level management to reduce mortality risk, protection of core areas and corridors, development of mitigation strategies to reduce vehicle collisions, and conflict mediation between local communities and onagers.

Identifying factors that create and maintain a hybrid zone is of great interest to ecology, evolution and, more recently, conservation biology. Here, we investigated the role of environmental features in shaping the spatial dynamics of a hybrid zone between the southern tigrina, Leopardus guttulus, and Geoffroy's cat, L. geoffroyi, testing for exogenous selection as the main force acting on its maintenance. These Neotropical felid species are mainly allopatric, with a restricted area of sympatry in the ecotone between the Atlantic Forest and Pampa biomes. As both biomes have experienced high rates of anthropogenic habitat alteration, we also analysed the influence of habitat conversion on the hybrid zone structure. To do this, we used 13 microsatellite loci to identify potential hybrids and generated ecological niche models for them and their parental species. We compared the influence of variables on parental species and hybrid occurrence and calculated the amount of niche overlap among them. Parental species showed different habitat requirements and predicted co-occurrence was restricted to the forest-grassland mosaic of the ecotone. However, hybrids were found beyond this area, mainly in the range of L. geoffroyi. Hybrids demonstrated higher tolerance to habitat alteration than parental types, with a probability of occurrence that was positively related with mosaics of cropland areas and remnants of natural vegetation. These results indicate that exogenous selection alone does not drive the dynamics of the hybrid zone, and that habitat conversion influences its structure, potentially favouring hybrids over parental species., (© 2021 European Society for Evolutionary Biology. Journal of Evolutionary Biology © 2021 European Society for Evolutionary Biology.)

Extant conifer species may be susceptible to rapid environmental change owing to their long generation times, but could also be resilient due to high levels of standing genetic diversity. Hybridisation between closely related species can increase genetic diversity and generate novel allelic combinations capable of fuelling adaptive evolution. Our study unravelled the genetic architecture of adaptive evolution in a conifer hybrid zone formed between Pinus strobiformis and P. flexilis. Using a multifaceted approach emphasising the spatial and environmental patterns of linkage disequilibrium and ancestry enrichment, we identified recently introgressed and background genetic variants to be driving adaptive evolution along different environmental gradients. Specifically, recently introgressed variants from P. flexilis were favoured along freeze-related environmental gradients, while background variants were favoured along water availability-related gradients. We posit that such mosaics of allelic variants within conifer hybrid zones will confer upon them greater resilience to ongoing and future environmental change and can be a key resource for conservation efforts.

An assumption of correlative landscape genetic methods is that genetic differentiation at neutral markers arises solely from the degree to which the intervening landscape between individuals or populations resists gene flow. However, this assumption is violated when gene flow occurs into the sampled population from an unsampled, differentiated deme. This may happen when sampling within only a portion of a population's extent or when closely related species hybridize with the sampled population. In both cases, violation of the modelling assumptions has the potential to reduce landscape genetic model selection accuracy and result in poor inferences. We used individual-based population genetic simulations in complex landscapes within a model selection framework to explore the potential confounding effect of gene flow from unsampled demes. We hypothesized that as gene flow from outside the sampling extent increased, model selection accuracy would decrease due to the formation of a hybrid zone where allele frequencies were perturbed in a way that was not correlated with effective distances between sampled individuals. Surprisingly, we found this expectation was unfounded, because the reduced accuracy due to admixture was counteracted by an increase in allelic diversity as alleles spread from the unsampled deme into the sampled population. These new alleles increased the power to detect landscape genetic relationships and even slightly improving model selection accuracy overall. This is a reassuring result, suggesting that sampling the full extent of a population or related species that may hybridize may be unnecessary, as long as other well-established sampling requirements are met., (© 2020 John Wiley & Sons Ltd.)

Iran lies at the southernmost range limit of brown bears globally. Therefore, understanding the habitat associations and patterns of population connectivity for brown bears in Iran is relevant for the species' conservation. We applied species distribution modeling to predict habitat suitability and connectivity modeling to identify population core areas and corridors. Our results showed that forest density, topographical roughness, NDVI and human footprint were the most influential variables in predicting brown bear distribution. The most crucial core areas and corridor networks for brown bear are concentrated in the Alborz and Zagros Mountains. These two core areas were predicted to be fragmented into a total of fifteen isolated patches if dispersal of brown bear across the landscape is limited to 50,000 cost units, and aggregates into two isolated habitat patches if the species is capable of dispersing 400,000 cost units. We found low overlap between corridors, and core habitats with protected areas, suggesting that the existing protected area network may not be adequate for the conservation of brown bear in Iran. Our results suggest that effective conservation of brown bears in Iran requires protection of both core habitats and the corridors between them, especially outside Iran's network of protected areas.

The phenotype of trees is determined by the relationships and interactions among genetic and environmental influences. Understanding the patterns and processes that are responsible for phenotypic variation is facilitated by studying the relationships between phenotype and the environment among many individuals across broad ecological and climatic gradients. We used Pinus strobiformis , which has a wide latitudinal distribution, as a model species to: (a) estimate the relative importance of different environmental factors in predicting these morphological traits and (b) characterize the spatial patterns of standing phenotypic variation of cone and seed traits across the species' range. A large portion of the total variation in morphological characteristics was explained by ecological, climatic and geographical variables (54.7% collectively). The three climate, vegetation and geographical variable groups, each had similar total ability to explain morphological variation (43.4%, 43.8%, 51.5%, respectively), while the topographical variable group had somewhat lower total explanatory power (36.9%). The largest component of explained variance (33.6%) was the four-way interaction of all variable sets, suggesting that there is strong covariation in environmental, climate and geographical variables in their relationship to morphological traits of southwest white pine across its range. The regression results showed that populations in more humid and warmer climates expressed greater cone length and seed size. This may in part facilitate populations of P. strobiformis in warmer and wetter portions of its range growing in dense, shady forest stands, because larger seeds provide greater resources to germinants at the time of germination. Our models provide accurate predictions of morphological traits and important insights regarding the factors that contribute to their expression. Our results indicate that managers should be conservative during reforestation efforts to ensure match between ecotypic variation in seed source populations. However, we also note that given projected large range shift due to climate change, managers will have to balance the match between current ecotypic variation and expected range shift and changes in local adaptive optima under future climate conditions., (Copyright © 2020 Leal-Sáenz, Waring, Menon, Cushman, Eckert, Flores-Rentería, Hernández-Díaz, López-Sánchez, Martínez-Guerrero and Wehenkel.)

Understanding the genetic basis underlying the local adaptation of nonmodel species is a fundamental goal in evolutionary biology. In this study, we explored the genetic mechanisms of the local adaptation of Forsythia suspensa using genome sequence and population genomics data obtained from specific-locus amplified fragment sequencing. We assembled a high-quality reference genome of weeping forsythia (Scaffold N50 = 7.3 Mb) using ultralong Nanopore reads. Then, genome-wide comparative analysis was performed for 15 natural populations of weeping forsythia across its current distribution range. Our results revealed that candidate genes associated with local adaptation are functionally correlated with solar radiation, temperature and water variables across heterogeneous environmental scenarios. In particular, solar radiation during the period of fruit development and seed drying after ripening, cold, and drought significantly contributed to the adaptive differentiation of F. suspensa . Natural selection exerted by environmental factors contributed substantially to the population genetic structure of F. suspensa . Our results supported the hypothesis that adaptive differentiation should be highly pronounced in the genes involved in signal crosstalk between different environmental variables. Our population genomics study of F. suspensa provides insights into the fundamental genetic mechanisms of the local adaptation of plant species to climatic gradients., Competing Interests: Conflict of interestThe authors declare that they have no conflict of interest., (© The Author(s) 2020.)

Replicated multiple scale species distribution models (SDMs) have become increasingly important to identify the correct variables determining species distribution and their influences on ecological responses. This study explores multi-scale habitat relationships of the snow leopard ( Panthera uncia ) in two study areas on the Qinghai-Tibetan Plateau of western China. Our primary objectives were to evaluate the degree to which snow leopard habitat relationships, expressed by predictors, scales of response, and magnitude of effects, were consistent across study areas or locally landcape-specific. We coupled univariate scale optimization and the maximum entropy algorithm to produce multivariate SDMs, inferring the relative suitability for the species by ensembling top performing models. We optimized the SDMs based on average omission rate across the top models and ensembles' overlap with a simulated reference model. Comparison of SDMs in the two study areas highlighted landscape-specific responses to limiting factors. These were dependent on the effects of the hydrological network, anthropogenic features, topographic complexity, and the heterogeneity of the landcover patch mosaic. Overall, even accounting for specific local differences, we found general landscape attributes associated with snow leopard ecological requirements, consisting of a positive association with uplands and ridges, aggregated low-contrast landscapes, and large extents of grassy and herbaceous vegetation. As a means to evaluate the performance of two bias correction methods, we explored their effects on three datasets showing a range of bias intensities. The performance of corrections depends on the bias intensity; however, density kernels offered a reliable correction strategy under all circumstances. This study reveals the multi-scale response of snow leopards to environmental attributes and confirms the role of meta-replicated study designs for the identification of spatially varying limiting factors. Furthermore, this study makes important contributions to the ongoing discussion about the best approaches for sampling bias correction., Competing Interests: The Authors declare no competing interests., (© 2020 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

Agricultural intensification, with its associated habitat loss and fragmentation, is among the most important drivers of the ongoing pollination crisis. In this quasi-experimental study, conducted in intensively managed vineyards in southwestern Switzerland, we tested the separate and interdependent effects of habitat amount and fragmentation on the foraging activity and reproductive performance of bumblebee Bombus t. terrestris colonies. Based on a factorial design, we selected a series of spatially replicated study sites across a dual gradient of habitat amount (area of ground-vegetated vineyards) and fragmentation (density of ground-vegetated vineyard fields) in a landscape predominantly consisting of vineyards with bare grounds. The foraging activity of individual bumblebees was measured using the radio frequency identification (RFID) technology, and we assessed final colony size to measure reproductive performance. We found an interactive effect of habitat amount and fragmentation on colony size. More specifically, the degree of fragmentation had a negative effect on bumblebee colony size when the amount of habitat was low, while it had a weak positive effect on colony size in landscapes with high amounts of habitat. At the level of individual vineyard fields, ground vegetation cover exerted a positive effect on bumblebee colony size. Fragmentation, but not habitat amount, significantly influenced foraging activity, with more foraging trips in sites with lower degrees of fragmentation. Our results emphasise the importance of studying the separate and interdependent effects of habitat amount and fragmentation to understand their influence on pollinators, providing guidance for optimising the spatial configuration of agricultural landscapes from a biodiversity viewpoint.

We implemented multilocus selection in a spatially-explicit, individual-based framework that enables multivariate environmental gradients to drive selection in many loci as a new module for the landscape genetics programs, CDPOP and CDMetaPOP. Our module simulates multilocus selection using a linear additive model, providing a flexible platform to evaluate a wide range of genotype-environment associations. Importantly, the module allows simulation of selection in any number of loci under the influence of any number of environmental variables. We validated the module with individual-based selection simulations under Wright-Fisher assumptions. We then evaluated results for simulations under a simple landscape selection model. Next, we simulated individual-based multilocus selection across a complex selection landscape with three loci linked to three different environmental variables. Finally, we demonstrated how the program can be used to simulate multilocus selection under varying selection strengths across different levels of gene flow in a landscape genetics framework. This new module provides a valuable addition to the study of landscape genetics, allowing for explicit evaluation of the contributions and interactions between gene flow and selection-driven processes across complex, multivariate environmental and landscape conditions., (© 2019 John Wiley & Sons Ltd.)

Background: In forest trees, genetic markers have been used to understand the genetic architecture of natural populations, identify quantitative trait loci, infer gene function, and enhance tree breeding. Recently, new, efficient technologies for genotyping thousands to millions of single nucleotide polymorphisms (SNPs) have finally made large-scale use of genetic markers widely available. These methods will be exceedingly valuable for improving tree breeding and understanding the ecological genetics of Douglas-fir, one of the most economically and ecologically important trees in the world., Results: We designed SNP assays for 55,766 potential SNPs that were discovered from previous transcriptome sequencing projects. We tested the array on ~ 2300 related and unrelated coastal Douglas-fir trees (Pseudotsuga menziesii var. menziesii) from Oregon and Washington, and 13 trees of interior Douglas-fir (P. menziesii var. glauca). As many as ~ 28 K SNPs were reliably genotyped and polymorphic, depending on the selected SNP call rate. To increase the number of SNPs and improve genome coverage, we developed protocols to 'rescue' SNPs that did not pass the default Affymetrix quality control criteria (e.g., 97% SNP call rate). Lowering the SNP call rate threshold from 97 to 60% increased the number of successful SNPs from 20,669 to 28,094. We used a subset of 395 unrelated trees to calculate SNP population genetic statistics for coastal Douglas-fir. Over a range of call rate thresholds (97 to 60%), the median call rate for SNPs in Hardy-Weinberg equilibrium ranged from 99.2 to 99.7%, and the median minor allele frequency ranged from 0.198 to 0.233. The successful SNPs also worked well on interior Douglas-fir., Conclusions: Based on the original transcriptome assemblies and comparisons to version 1.0 of the Douglas-fir reference genome, we conclude that these SNPs can be used to genotype about 10 K to 15 K loci. The Axiom genotyping array will serve as an excellent foundation for studying the population genomics of Douglas-fir and for implementing genomic selection. We are currently using the array to construct a linkage map and test genomic selection in a three-generation breeding program for coastal Douglas-fir.

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., Competing Interests: None declared., (© 2019 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Habitat loss is the greatest threat to biodiversity in Borneo, and to anticipate and combat its effects it is important to predict the pattern of loss and its consequences. Borneo is a region of extremely high biodiversity from which forest is being lost faster than in any other. The little-known Sunda clouded leopard (Neofelis diardi) is the top predator in Borneo and is likely to depend critically on habitat connectivity that is currently being rapidly lost to deforestation. We modeled the effects of landscape fragmentation on population size, genetic diversity and population connectivity for the Sunda clouded leopard across the entirety of Borneo. We modelled the impacts of land use change between the years 2000, 2010 and projected forwards to 2020. We found substantial reductions across all metrics between 2000 and 2010: the proportion of landscape connected by dispersal fell by approximately 12.5% and the largest patch size declined by around 15.1%, leading to a predicted 11.4% decline in clouded leopard numbers. We also predict that these trends will accelerate greatly towards 2020, with the percentage of the landscape being connected by dispersal falling by about 57.8%, the largest patch size falling by around 62.8% and the predicted clouded leopard population falling by 62.5% between 2010 and 2020. We predicted that these large declines in clouded leopard population size and connectivity will also substantially reduce the genetic diversity of the remaining clouded leopard population., Competing Interests: The authors have declared that no competing interests exist.

We evaluated how differences between two empirical resistance models for the same geographic area affected predictions of gene flow processes and genetic diversity for the Mexican spotted owl ( Strix occidentalis lucida ). The two resistance models represented the landscape under low- and high-fragmentation parameters. Under low fragmentation, the landscape had larger but highly concentrated habitat patches, whereas under high fragmentation, the landscape had smaller habitat patches that scattered across a broader area. Overall habitat amount differed little between resistance models. We tested eight scenarios reflecting a factorial design of three factors: resistance model (low vs. high fragmentation), isolation hypothesis (isolation-by-distance, IBD, vs. isolation-by-resistance, IBR), and dispersal limit of species (200 km vs. 300 km). Higher dispersal limit generally had a positive but small influence on genetic diversity. Genetic distance increased with both geographic distance and landscape resistance, but landscape resistance displayed a stronger influence. Connectivity was positively related to genetic diversity under IBR but was less important under IBD. Fragmentation had a strong negative influence on the spatial patterns of genetic diversity and effective population size (Ns). Despite habitats being more concentrated and less widely distributed, the low-fragmentation landscape had greater genetic diversity than the high-fragmentation landscape, suggesting that highly concentrated but larger habitat patches may provide a genetic refuge for the Mexican spotted owl.

Niche differentiation, the partitioning of resources along one or more axes of a species' niche hyper-volume, is widely recognised as an important mechanism for sympatric species to reduce interspecific competition and predation risk, and thus facilitate co-existence. Resource partitioning may be facilitated by behavioural differentiation along three main niche dimensions: habitat, food and time. In this study, we investigate the extent to which these mechanisms can explain the coexistence of an assemblage of five sympatric felids in Borneo. Using multi-scale logistic regression, we show that Bornean felids exhibit differences in both their broad and fine-scale habitat use. We calculate temporal activity patterns and overlap between these species, and present evidence for temporal separation within this felid guild. Lastly, we conducted an all-subsets logistic regression to predict the occurrence of each felid species as a function of the co-occurrence of a large number of other species and showed that Bornean felids co-occurred with a range of other species, some of which could be candidate prey. Our study reveals apparent resource partitioning within the Bornean felid assemblage, operating along all three niche dimension axes. These results provide new insights into the ecology of these species and the broader community in which they live and also provide important information for conservation planning for this guild of predators., Competing Interests: The authors have declared that no competing interests exist.

Conservation of large carnivores, such as the African lion, requires preservation of extensive core habitat areas, linkages between them, and mitigation of human-wildlife conflict. However, there are few rigorous examples of efforts that prioritized conservation actions for all three of these critical components. We used an empirically optimized resistance surface to calculate resistant kernel and factorial least cost path predictions of population connectivity and conflict risk for lions across the Kavango-Zambezi Transfrontier Conservation Area (KAZA) and surrounding landscape. We mapped and ranked the relative importance of (1) lion dispersal areas outside National Parks, (2) corridors between the key areas, and (3) areas of highest human-lion conflict risk. Spatial prioritization of conservation actions is critical given extensive land use redesignations that are reducing the extent and increasing the fragmentation of lion populations. While our example focuses on lions in southern Africa, it provides a general approach for rigorous, empirically based comprehensive conservation planning based on spatial prioritization., Competing Interests: The authors have declared that no competing interests exist.

Entropy and the second law of thermodynamics are fundamental concepts that underlie all natural processes and patterns. Recent research has shown how the entropy of a landscape mosaic can be calculated using the Boltzmann equation, with the entropy of a lattice mosaic equal to the logarithm of the number of ways a lattice with a given dimensionality and number of classes can be arranged to produce the same total amount of edge between cells of different classes. However, that work seemed to also suggest that the feasibility of applying this method to real landscapes was limited due to intractably large numbers of possible arrangements of raster cells in large landscapes. Here I extend that work by showing that: (1) the proportion of arrangements rather than the number with a given amount of edge length provides a means to calculate unbiased relative configurational entropy, obviating the need to compute all possible configurations of a landscape lattice; (2) the edge lengths of randomized landscape mosaics are normally distributed, following the central limit theorem; and (3) given this normal distribution it is possible to fit parametric probability density functions to estimate the expected proportion of randomized configurations that have any given edge length, enabling the calculation of configurational entropy on any landscape regardless of size or number of classes. I evaluate the boundary limits (4) for this normal approximation for small landscapes with a small proportion of a minority class and show it holds under all realistic landscape conditions. I further (5) demonstrate that this relationship holds for a sample of real landscapes that vary in size, patch richness, and evenness of area in each cover type, and (6) I show that the mean and standard deviation of the normally distributed edge lengths can be predicted nearly perfectly as a function of the size, patch richness and diversity of a landscape. Finally, (7) I show that the configurational entropy of a landscape is highly related to the dimensionality of the landscape, the number of cover classes, the evenness of landscape composition across classes, and landscape heterogeneity. These advances provide a means for researchers to directly estimate the frequency distribution of all possible macrostates of any observed landscape, and then directly calculate the relative configurational entropy of the observed macrostate, and to understand the ecological meaning of different amounts of configurational entropy. These advances enable scientists to take configurational entropy from a concept to an applied tool to measure and compare the disorder of real landscapes with an objective and unbiased measure based on entropy and the second law.

Interactions between extrinsic factors, such as disruptive selection and intrinsic factors, such as genetic incompatibilities among loci, often contribute to the maintenance of species boundaries. The relative roles of these factors in the establishment of reproductive isolation can be examined using species pairs characterized by gene flow throughout their divergence history. We investigated the process of speciation and the maintenance of species boundaries between Pinus strobiformis and Pinus flexilis. Utilizing ecological niche modelling, demographic modelling and genomic cline analyses, we illustrated a divergence history with continuous gene flow. Our results supported an abundance of advanced generation hybrids and a lack of loci exhibiting steep transition in allele frequency across the hybrid zone. Additionally, we found evidence for climate-associated variation in the hybrid index and niche divergence between parental species and the hybrid zone. These results are consistent with extrinsic factors, such as climate, being an important isolating mechanism. A build-up of intrinsic incompatibilities and of coadapted gene complexes is also apparent, although these appear to be in the earliest stages of development. This supports previous work in coniferous species demonstrating the importance of extrinsic factors in facilitating speciation. Overall, our findings lend support to the hypothesis that varying strength and direction of selection pressures across the long lifespans of conifers, in combination with their other life history traits, delays the evolution of strong intrinsic incompatibilities., (© 2018 John Wiley & Sons Ltd.)

Anthropogenic migration barriers fragment many populations and limit the ability of species to respond to climate-induced biome shifts. Conservation actions designed to conserve habitat connectivity and mitigate barriers are needed to unite fragmented populations into larger, more viable metapopulations, and to allow species to track their climate envelope over time. Landscape genetic analysis provides an empirical means to infer landscape factors influencing gene flow and thereby inform such conservation actions. However, there are currently many methods available for model selection in landscape genetics, and considerable uncertainty as to which provide the greatest accuracy in identifying the true landscape model influencing gene flow among competing alternative hypotheses. In this study, we used population genetic simulations to evaluate the performance of seven regression-based model selection methods on a broad array of landscapes that varied by the number and type of variables contributing to resistance, the magnitude and cohesion of resistance, as well as the functional relationship between variables and resistance. We also assessed the effect of transformations designed to linearize the relationship between genetic and landscape distances. We found that linear mixed effects models had the highest accuracy in every way we evaluated model performance; however, other methods also performed well in many circumstances, particularly when landscape resistance was high and the correlation among competing hypotheses was limited. Our results provide guidance for which regression-based model selection methods provide the most accurate inferences in landscape genetic analysis and thereby best inform connectivity conservation actions., (Published 2017. This article is a U.S. Government work and is in the public domain in the USA.)

Human-carnivore conflict is a primary driver of carnivore declines worldwide and resolving these conflicts is a conservation priority. However, resources to mitigate conflicts are limited and should be focused on areas of highest priority. We conducted 820 semistructured interviews with community members living within Kenya's Maasai Mara ecosystem. A multiscale analysis was used to determine the influence of husbandry and environmental factors on livestock depredation inside livestock enclosures (bomas). Areas with a high proportion of closed habitat and protected areas had the highest risk of depredation. Depredation was most likely to occur at weak bomas and at households where there were fewer dogs. We used the results to identify potential conflict hotspots by mapping the probability of livestock depredation across the landscape. 21.4% of the landscape was classified as high risk, and within these areas, 53.4% of the households that were interviewed had weak bomas. Synthesis and applications . With limited resources available to mitigate human-carnivore conflicts, it is imperative that areas are identified where livestock is most at risk of depredation. Focusing mitigation measures on high-risk areas may reduce conflict and lead to a decrease in retaliatory killings of predators.

A major aim of landscape genetics is to understand how landscapes resist gene flow and thereby influence population genetic structure. An empirical understanding of this process provides a wealth of information that can be used to guide conservation and management of species in fragmented landscapes and also to predict how landscape change may affect population viability. Statistical approaches to infer the true model among competing alternatives are based on the strength of the relationship between pairwise genetic distances and landscape distances among sampled individuals in a population. A variety of methods have been devised to quantify individual genetic distances, but no study has yet compared their relative performance when used for model selection in landscape genetics. In this study, we used population genetic simulations to assess the accuracy of 16 individual-based genetic distance metrics under varying sample sizes and degree of population genetic structure. We found most metrics performed well when sample size and genetic structure was high. However, it was much more challenging to infer the true model when sample size and genetic structure was low. Under these conditions, we found genetic distance metrics based on principal components analysis were the most accurate (although several other metrics performed similarly), but only when they were derived from multiple principal components axes (the optimal number varied depending on the degree of population genetic structure). Our results provide guidance for which genetic distance metrics maximize model selection accuracy and thereby better inform conservation and management decisions based upon landscape genetic analysis., (Published 2017. This article is a U.S. Government work and is in the public domain in the USA.)

Gene flow is an evolutionary process that supports genetic connectivity and contributes to the capacity of species to adapt to environmental change. Yet, for most species, little is known about the specific environmental factors that influence genetic connectivity, or their effects on genetic diversity and differentiation. We used a landscape genetic approach to understand how geography and climate influence genetic connectivity in a foundation riparian tree (Populus angustifolia), and their relationships with specieswide patterns of genetic diversity and differentiation. Using multivariate restricted optimization in a reciprocal causal modelling framework, we quantified the relative contributions of riparian network connectivity, terrestrial upland resistance and climate gradients on genetic connectivity. We found that (i) all riparian corridors, regardless of river order, equally facilitated connectivity, while terrestrial uplands provided 2.5× more resistance to gene flow than riparian corridors. (ii) Cumulative differences in precipitation seasonality and precipitation of the warmest quarter were the primary climatic factors driving genetic differentiation; furthermore, maximum climate resistance was 45× greater than riparian resistance. (iii) Genetic diversity was positively correlated with connectivity (R 2  = 0.3744, p = .0019), illustrating the utility of resistance models for identifying landscape conditions that can support a species' ability to adapt to environmental change. From these results, we present a map highlighting key genetic connectivity corridors across P. angustifolia's range that if disrupted could have long-term ecological and evolutionary consequences. Our findings provide recommendations for conservation and restoration management of threatened riparian ecosystems throughout the western USA and the high biodiversity they support., (© 2017 John Wiley & Sons Ltd.)

Winter bait stations are becoming a commonly used technique for multispecies inventory and monitoring but a technical evaluation of their effectiveness is lacking. Bait stations have three components: carcass attractant, remote camera, and hair snare. Our 22,975 km 2 mountainous study area was stratified with a 5 × 5 km sampling grid centered on northern Idaho and including portions of Washington, Montana, and British Columbia. From 2010-14, we conducted 563 sampling sessions at 497 bait stations in 453 5 × 5 km cells. We evaluated the effectiveness of cameras and hair snare DNA collection at stations to detect species and individual animals, factors affecting DNA viability, the effectiveness of re-visiting stations, and the influence of elevation, seasonality, and latency on species detections. Cameras were more effective at detecting multiple species than DNA hair snaring. Length of deployment time and elevation increased genetic species ID success but individual ID success rates were increased only by collecting hairs earlier in the season. Re-visiting stations did not change camera or genetic species detection results but did increase the number of individual genotypes identified. Marten and fisher were detected quickly while bobcat and coyote showed longer latency to detection. Seasonality significantly affected coyote and bobcat detections but not marten, fisher, or weasel. Multispecies bait station study design should incorporate mixed elevation sites with stratified seasonality. Priority should be given to including cameras as components of bait stations over hair snares, unless there is a specific genetic goal to the study. A hair snare component should be added, however, if individual ID or genetic data are necessary. Winter stations should be deployed a minimum of 45-60 days to allow for detection of low density species and species with long latency to detection times. Hair samples should be collected prior to DNA-degrading late season rain events. Re-visiting stations does not change which species are detected at stations; therefore, studies with objectives to delineate species presence or distribution will be more effective if they focus on deploying more stations across a broader landscape in lieu of surveying the same site multiple times.

Recent population declines to the high elevation western North America foundation species whitebark pine, have been driven by the synergistic effects of the invasive blister rust pathogen, mountain pine beetle (MPB), fire exclusion, and climate change. This has led to consideration for listing whitebark pine (WBP) as a threatened or endangered species under the Endangered Species Act, which has intensified interest in developing management strategies for maintaining and restoring the species. An important, but poorly studied, aspect of WBP restoration is the spatial variation in adaptive genetic variation and the potential of blister rust resistant strains to maintain viable populations in the future. Here, we present a simulation modeling framework to improve understanding of the long-term genetic consequences of the blister rust pathogen, the evolution of rust resistance, and scenarios of planting rust resistant genotypes of whitebark pine. We combine climate niche modeling and eco-evolutionary landscape genetics modeling to evaluate the effects of different scenarios of planting rust-resistant genotypes and impacts of wind field direction on patterns of gene flow. Planting scenarios showed different levels for local extirpation of WBP and increased population-wide blister rust resistance, suggesting that the spatial arrangement and choice of planting locations can greatly affect survival rates of whitebark pine. This study presents a preliminary, but potentially important, framework for facilitating the conservation of whitebark pine.