Great tits (Parus major) flexibly learn that herbivore‐induced plant volatiles indicate prey location: An experimental evidence with two tree species

1. When searching for food, great tits (Parus major) can use herbivore‐induced plant volatiles (HIPVs) as an indicator of arthropod presence. Their ability to detect HIPVs was shown to be learned, and not innate, yet the flexibility and generalization of learning remain unclear.
2. We studied if, and if so how, naïve and trained great tits (Parus major) discriminate between herbivore‐induced and noninduced saplings of Scotch elm (Ulmus glabra) and cattley guava (Psidium cattleyanum). We chemically analyzed the used plants and showed that their HIPVs differed significantly and overlapped only in a few compounds.
3. Birds trained to discriminate between herbivore‐induced and noninduced saplings preferred the herbivore‐induced saplings of the plant species they were trained to. Naïve birds did not show any preferences. Our results indicate that the attraction of great tits to herbivore‐induced plants is not innate, rather it is a skill that can be acquired through learning, one tree species at a time.
4. We demonstrate that the ability to learn to associate HIPVs with food reward is flexible, expressed to both tested plant species, even if the plant species has not coevolved with the bird species (i.e., guava). Our results imply that the birds are not capable of generalizing HIPVs among tree species but suggest that they either learn to detect individual compounds or associate whole bouquets with food rewards.

     

Targeted and passive environmental DNA approaches outperform established methods for detection of quagga mussels,Dreissena rostriformis bugensis in flowing water

1. The early detection of invasive non‐native species (INNS) is important for informing management actions. Established monitoring methods require the collection or observation of specimens, which is unlikely at the beginning of an invasion when densities are likely to be low. Environmental DNA (eDNA) analysis is a highly promising technique for the detection of INNS—particularly during the early stages of an invasion.
2. Here, we compared the use of traditional kick‐net sampling with two eDNA approaches (targeted detection using both conventional and quantitative PCR and passive detection via metabarcoding with conserved primers) for detection of quagga mussel, Dreissena rostriformis bugensis, a high priority INNS, along a density gradient on the River Wraysbury, UK.
3. All three molecular tools outperformed traditional sampling in terms of detection. Conventional PCR and qPCR both had 100% detection rate in all samples and outperformed metabarcoding when the target species was at low densities. Additionally, quagga mussel DNA copy number (qPCR) and relative read count (metabarcoding) were significantly influenced by both mussel density and distance from source population, with distance being the most significant predictor.
4. Synthesis and application. All three molecular approaches were more sensitive than traditional kick‐net sampling for the detection of the quagga mussel in flowing water, and both qPCR and metabarcoding enabled estimates of relative abundance. Targeted approaches were more sensitive than metabarcoding, but metabarcoding has the advantage of providing information on the wider community and consequently the impacts of INNS.

     

Agent‐based modeling of the effects of forest dynamics,selective logging,and fragment size on epiphyte communities

1. Forest canopies play a crucial role in structuring communities of vascular epiphytes by providing substrate for colonization, by locally varying microclimate, and by causing epiphyte mortality due to branch or tree fall. However, as field studies in the three‐dimensional habitat of epiphytes are generally challenging, our understanding of how forest structure and dynamics influence the structure and dynamics of epiphyte communities is scarce.
2. Mechanistic models can improve our understanding of epiphyte community dynamics. We present such a model that couples dispersal, growth, and mortality of individual epiphytes with substrate dynamics, obtained from a three‐dimensional functional–structural forest model, allowing the study of forest–epiphyte interactions. After validating the epiphyte model with independent field data, we performed several theoretical simulation experiments to assess how (a) differences in natural forest dynamics, (b) selective logging, and (c) forest fragmentation could influence the long‐term dynamics of epiphyte communities.
3. The proportion of arboreal substrate occupied by epiphytes (i.e., saturation level) was tightly linked with forest dynamics and increased with decreasing forest turnover rates. While species richness was, in general, negatively correlated with forest turnover rates, low species numbers in forests with very‐low‐turnover rates were due to competitive exclusion when epiphyte communities became saturated. Logging had a negative impact on epiphyte communities, potentially leading to a near‐complete extirpation of epiphytes when the simulated target diameters fell below a threshold. Fragment size had no effect on epiphyte abundance and saturation level but correlated positively with species numbers.
4. Synthesis: The presented model is a first step toward studying the dynamic forest–epiphyte interactions in an agent‐based modeling framework. Our study suggests forest dynamics as key factor in controlling epiphyte communities. Thus, both natural and human‐induced changes in forest dynamics, for example, increased mortality rates or the loss of large trees, pose challenges for epiphyte conservation.

     

Species distribution models for invasive Eurasian watermilfoil highlight the importance of data quality and limitations of discrimination accuracy metrics

AimAvailability of uniformly collected presence, absence, and abundance data remains a key challenge in species distribution modeling (SDM). For invasive species, abundance and impacts are highly variable across landscapes, and quality occurrence and abundance data are critical for predicting locations at high risk for invasion and impacts, respectively. We leverage a large aquatic vegetation dataset comprising point‐level survey data that includes information on the invasive plant Myriophyllum spicatum (Eurasian watermilfoil) to: (a) develop SDMs to predict invasion and impact from environmental variables based on presence–absence, presence‐only, and abundance data, and (b) compare evaluation metrics based on functional and discrimination accuracy for presence–absence and presence‐only SDMs.LocationMinnesota, USA.MethodsEurasian watermilfoil presence–absence and abundance information were gathered from 468 surveyed lakes, and 801 unsurveyed lakes were leveraged as pseudoabsences for presence‐only models. A Random Forest algorithm was used to model the distribution and abundance of Eurasian watermilfoil as a function of lake‐specific predictors, both with and without a spatial autocovariate. Occurrence‐based SDMs were evaluated using conventional discrimination accuracy metrics and functional accuracy metrics assessing correlation between predicted suitability and observed abundance.ResultsWater temperature degree days and maximum lake depth were two leading predictors influencing both invasion risk and abundance, but they were relatively less important for predicting abundance than other water quality measures. Road density was a strong predictor of Eurasian watermilfoil invasion risk but not abundance. Model evaluations highlighted significant differences: Presence–absence models had high functional accuracy despite low discrimination accuracy, whereas presence‐only models showed the opposite pattern.Main conclusionComplementing presence–absence data with abundance information offers a richer understanding of invasive Eurasian watermilfoil''s ecological niche and enables evaluation of the model''s functional accuracy. Conventional discrimination accuracy measures were misleading when models were developed using pseudoabsences. We thus caution against the overuse of presence‐only models and suggest directing more effort toward systematic monitoring programs that yield high‐quality data.     

Mixing tree species associated with arbuscular or ectotrophic mycorrhizae reveals dual mycorrhization and interactive effects on the fungal partners

1. Recent studies found that the majority of shrub and tree species are associated with both arbuscular mycorrhizal (AM) and ectomycorrhizal (EM) fungi. However, our knowledge on how different mycorrhizal types interact with each other is still limited. We asked whether the combination of hosts with a preferred association with either AM or EM fungi increases the host tree roots’ mycorrhization rate and affects AM and EM fungal richness and community composition.
2. We established a tree diversity experiment, where five tree species of each of the two mycorrhiza types were planted in monocultures, two‐species and four‐species mixtures. We applied morphological assessment to estimate mycorrhization rates and next‐generation molecular sequencing to quantify mycobiont richness.
3. Both the morphological and molecular assessment revealed dual‐mycorrhizal colonization in 79% and 100% of the samples, respectively. OTU community composition strongly differed between AM and EM trees. While host tree species richness did not affect mycorrhization rates, we observed significant effects of mixing AM‐ and EM‐associated hosts in AM mycorrhization rate. Glomeromycota richness was larger in monotypic AM tree combinations than in AM‐EM mixtures, pointing to a dilution or suppression effect of AM by EM trees. We found a strong match between morphological quantification of AM mycorrhization rate and Glomeromycota richness.
4. Synthesis. We provide evidence that the combination of hosts differing in their preferred mycorrhiza association affects the host''s fungal community composition, thus revealing important biotic interactions among trees and their associated fungi.

     

Plotting receiver operating characteristic and precision–recall curves from presence and background data

1. The receiver operating characteristic (ROC) and precision–recall (PR) plots have been widely used to evaluate the performance of species distribution models. Plotting the ROC/PR curves requires a traditional test set with both presence and absence data (namely PA approach), but species absence data are usually not available in reality. Plotting the ROC/PR curves from presence‐only data while treating background data as pseudo absence data (namely PO approach) may provide misleading results.
2. In this study, we propose a new approach to calibrate the ROC/PR curves from presence and background data with user‐provided information on a constant c, namely PB approach. Here, c defines the probability that species occurrence is detected (labeled), and an estimate of c can also be derived from the PB‐based ROC/PR plots given that a model with good ability of discrimination is available. We used five virtual species and a real aerial photography to test the effectiveness of the proposed PB‐based ROC/PR plots. Different models (or classifiers) were trained from presence and background data with various sample sizes. The ROC/PR curves plotted by PA approach were used to benchmark the curves plotted by PO and PB approaches.
3. Experimental results show that the curves and areas under curves by PB approach are more similar to that by PA approach as compared with PO approach. The PB‐based ROC/PR plots also provide highly accurate estimations of c in our experiment.
4. We conclude that the proposed PB‐based ROC/PR plots can provide valuable complements to the existing model assessment methods, and they also provide an additional way to estimate the constant c (or species prevalence) from presence and background data.

     

How to make use of unlabeled observations in species distribution modeling using point process models

1. Species distribution modeling, which allows users to predict the spatial distribution of species with the use of environmental covariates, has become increasingly popular, with many software platforms providing tools to fit such models. However, the species observations used can have varying levels of quality and can have incomplete information, such as uncertain or unknown species identity.
2. In this paper, we develop two algorithms to classify observations with unknown species identities which simultaneously predict several species distributions using spatial point processes. Through simulations, we compare the performance of these algorithms using 7 different initializations to the performance of models fitted using only the observations with known species identity.
3. We show that performance varies with differences in correlation among species distributions, species abundance, and the proportion of observations with unknown species identities. Additionally, some of the methods developed here outperformed the models that did not use the misspecified data. We applied the best‐performing methods to a dataset of three frog species (Mixophyes).
4. These models represent a helpful and promising tool for opportunistic surveys where misidentification is possible or for the distribution of species newly separated in their taxonomy.

     

Modeling the demography of species providing extended parental care: A capture–recapture multievent model with a case study on polar bears (Ursus maritimus)

1. In species providing extended parental care, one or both parents care for altricial young over a period including more than one breeding season. We expect large parental investment and long‐term dependency within family units to cause high variability in life trajectories among individuals with complex consequences at the population level. So far, models for estimating demographic parameters in free‐ranging animal populations mostly ignore extended parental care, thereby limiting our understanding of its consequences on parents and offspring life histories.
2. We designed a capture–recapture multievent model for studying the demography of species providing extended parental care. It handles statistical multiple‐year dependency among individual demographic parameters grouped within family units, variable litter size, and uncertainty on the timing at offspring independence. It allows for the evaluation of trade‐offs among demographic parameters, the influence of past reproductive history on the caring parent''s survival status, breeding probability, and litter size probability, while accounting for imperfect detection of family units. We assess the model performance using simulated data and illustrate its use with a long‐term dataset collected on the Svalbard polar bears (Ursus maritimus).
3. Our model performed well in terms of bias and mean square error and in estimating demographic parameters in all simulated scenarios, both when offspring departure probability from the family unit occurred at a constant rate or varied during the field season depending on the date of capture. For the polar bear case study, we provide estimates of adult and dependent offspring survival rates, breeding probability, and litter size probability. Results showed that the outcome of the previous reproduction influenced breeding probability.
4. Overall, our results show the importance of accounting for i) the multiple‐year statistical dependency within family units, ii) uncertainty on the timing at offspring independence, and iii) past reproductive history of the caring parent. If ignored, estimates obtained for breeding probability, litter size, and survival can be biased. This is of interest in terms of conservation because species providing extended parental care are often long‐living mammals vulnerable or threatened with extinction.

     

Weather variation affects the dispersal of grasshoppers beyond their elevational ranges

1. Understanding how abiotic conditions influence dispersal patterns of organisms is important for understanding the degree to which species can track and persist in the face of changing climate.
2. The goal of this study was to understand how weather conditions influence the dispersal pattern of multiple nonmigratory grasshopper species from lower elevation grassland habitats in which they complete their life‐cycles to higher elevations that extend beyond their range limits.
3. Using over a decade of weekly spring to late‐summer field survey data along an elevational gradient, we explored how abundance and richness of dispersing grasshoppers were influenced by temperature, precipitation, and wind speed and direction. We also examined how changes in population sizes at lower elevations might influence these patterns.
4. We observed that the abundance of dispersing grasshoppers along the gradient declined 4‐fold from the foothills to the subalpine and increased with warmer conditions and when wind flow patterns were mild or in the downslope direction. Thirty‐eight unique grasshopper species from lowland sites were detected as dispersers across the survey years, and warmer years and weak upslope wind conditions also increased the richness of these grasshoppers. The pattern of grasshoppers along the gradient was not sex biased. The positive effect of temperature on dispersal rates was likely explained by an increase in dispersal propensity rather than by an increase in the density of grasshoppers at low elevation sites.
5. The results of this study support the hypothesis that the dispersal patterns of organisms are influenced by changing climatic conditions themselves and as such, that this context‐dependent dispersal response should be considered when modeling and forecasting the ability of species to respond to climate change.

     

Mixed interactions among life history stages of two harvested related species

Climate change and harvesting can affect the ecosystems'' functioning by altering the population dynamics and interactions among species. Knowing how species interact is essential for better understanding potentially unintended consequences of harvest on multiple species in ecosystems. I analyzed how stage‐specific interactions between two harvested competitors, the haddock (Melanogrammus aeglefinus) and Atlantic cod (Gadus morhua), living in the Barents Sea affect the outcome of changes in the harvest of the two species. Using state‐space models that account for observation errors and stochasticity in the population dynamics, I run different harvesting scenarios and track population‐level responses of both species. The increasing temperature elevated the number of larvae of haddock but did not significantly influence the older age‐classes. The nature of the interactions between both species shifted from predator‐prey to competition around age‐2 to ‐3. Increased cod fishing mortality, which led to decreasing abundance of cod, was associated with an increasing overall abundance of haddock, which suggests compensatory dynamics of both species. From a stage‐specific approach, I show that a change in the abundance in one species may propagate to other species, threatening the exploited species'' recovery. Thus, this study demonstrates that considering interactions among life history stages of harvested species is essential to enhance species'' co‐existence in harvested ecosystems. The approach developed in this study steps forward the analyses of effects of harvest and climate in multi‐species systems by considering the comprehension of complex ecological processes to facilitate the sustainable use of natural resources.     

Abundance of montane salamanders over an elevational gradient

Climate change is expected to systematically alter the distribution and population dynamics of species around the world. The effects are expected to be particularly strong at high latitudes and elevations, and for ectothermic species with small ranges and limited movement potential, such as salamanders in the southern Appalachian Mountains. In this study, we sought to establish baseline abundance estimates for plethodontid salamanders (family: Plethodontidae) over an elevational gradient in Great Smoky Mountains National Park. In addition to generating these baseline data for multiple species, we describe methods for surveying salamanders that allow for meaningful comparisons over time by separating observation and ecological processes generating the data. We found that Plethodon jordani had a mid‐elevation peak (1,500 m) in abundance and Desmognathus wrighti increased in abundance with elevation up to the highest areas of the park (2025 m), whereas Eurycea wilderae increased in abundance up to 1,600 m and then plateaued with increasing uncertainty. Litter depth, herbaceous ground cover, and proximity to stream were also important predictors of abundance (dependent upon species), whereas daily temperature, precipitation, ground cover, and humidity influenced detection rates. Our data provide some of the first minimally biased information for future studies to assess changes in the abundance and distribution of salamanders in this region. Understanding abundance patterns along with detailed baseline distributions will be critical for comparisons with future surveys to understand the population and community‐level effects of climate change on montane salamanders.     

Drivers and consequences of partial migration in an alpine bird species

1. Partial migration, where a portion of the population migrates between winter and summer (breeding) areas and the rest remain year‐round resident, is a common phenomenon across several taxonomic groups. Several hypotheses have been put forward to explain why some individuals migrate while others stay resident, as well as the fitness consequences of the different strategies. Yet, the drivers and consequences of the decision to migrate or not are poorly understood.
2. We used data from radio‐tagged female (n = 73) willow ptarmigan Lagopus lagopus in an alpine study area in Central Norway to test if (i) the decision to migrate was dependent on individual state variables (age and body weight), (ii) individuals repeated migratory decisions between seasons, and (iii) the choice of migratory strategy was related to reproductive success.
3. Partially supporting our prediction that migratory strategy depends on individual state, we found that juvenile birds with small body sizes were more likely to migrate, whereas large juveniles remained resident. For adult females, we found no relationship between the decision to migrate or stay resident and body weight. We found evidence for high individual repeatability of migratory decision between seasons. Migratory strategy did not explain variation in clutch size or nest fate among individuals, suggesting no direct influence of the chosen strategy on reproductive success.
4. Our results indicate that partial migration in willow ptarmigan is related to juvenile body weight, and that migratory behavior becomes a part of the individual life history as a fixed strategy. Nesting success was not affected by migratory strategy in our study population, but future studies should assess other traits to further test potential fitness consequences.

     

Remotely sensed between‐individual functional trait variation in a temperate forest

1. Trait‐based ecology holds the promise to explain how plant communities work, for example, how functional diversity may support community productivity. However, so far it has been difficult to combine field‐based approaches assessing traits at the level of plant individuals with limited spatial coverage and approaches using remote sensing (RS) with complete spatial coverage but assessing traits at the level of vegetation pixels rather than individuals. By delineating all individual‐tree crowns within a temperate forest site and then assigning RS‐derived trait measures to these trees, we combine the two approaches, allowing us to use general linear models to estimate the influence of taxonomic or environmental variation on between‐ and within‐species variation across contiguous space.
2. We used airborne imaging spectroscopy and laser scanning to collect individual‐tree RS data from a mixed conifer‐angiosperm forest on a mountain slope extending over 5.5 ha and covering large environmental gradients in elevation as well as light and soil conditions. We derived three biochemical (leaf chlorophyll, carotenoids, and water content) and three architectural traits (plant area index, foliage‐height diversity, and canopy height), which had previously been used to characterize plant function, from the RS data. We then quantified the contributions of taxonomic and environmental variation and their interaction to trait variation and partitioned the remaining within‐species trait variation into smaller‐scale spatial and residual variation. We also investigated the correlation between functional trait and phylogenetic distances at the between‐species level. The forest consisted of 13 tree species of which eight occurred in sufficient abundance for quantitative analysis.
3. On average, taxonomic variation between species accounted for more than 15% of trait variation in biochemical traits but only around 5% (still highly significant) in architectural traits. Biochemical trait distances among species also showed a stronger correlation with phylogenetic distances than did architectural trait distances. Light and soil conditions together with elevation explained slightly more variation than taxonomy across all traits, but in particular increased plant area index (light) and reduced canopy height (elevation). Except for foliage‐height diversity, all traits were affected by significant interactions between taxonomic and environmental variation, the different responses of the eight species to the within‐site environmental gradients potentially contributing to the coexistence of the eight abundant species.
4. We conclude that with high‐resolution RS data it is possible to delineate individual‐tree crowns within a forest and thus assess functional traits derived from RS data at individual level. With this precondition fulfilled, it is then possible to apply tools commonly used in field‐based trait ecology to partition trait variation among individuals into taxonomic and potentially even genetic variation, environmental variation, and interactions between the two. The method proposed here presents a promising way of assessing individual‐based trait information with complete spatial coverage and thus allowing analysis of functional diversity at different scales. This information can help to better understand processes shaping community structure, productivity, and stability of forests.

     

Climate and habitat configuration limit range expansion and patterns of dispersal in a non‐native lizard

Invasive species are one of the main causes of biodiversity loss worldwide. As introduced, populations increase in abundance and geographical range, so does the potential for negative impacts on native communities. As such, there is a need to better understand the processes driving range expansion as species become established in recipient landscapes. Through an investigation into capacity for population growth and range expansion of introduced populations of a non‐native lizard (Podarcis muralis), we aimed to demonstrate how multi‐scale factors influence spatial spread, population growth, and invasion potential in introduced species. We collated location records of P. muralis presence in England, UK through data collected from field surveys and a citizen science campaign. We used these data as input for presence‐background models to predict areas of climate suitability at a national‐scale (5 km resolution), and fine‐scale habitat suitability at the local scale (2 m resolution). We then integrated local models into an individual‐based modeling platform to simulate population dynamics and forecast range expansion for 10 populations in heterogeneous landscapes. National‐scale models indicated climate suitability has restricted the species to the southern parts of the UK, primarily by a latitudinal cline in overwintering conditions. Patterns of population growth and range expansion were related to differences in local landscape configuration and heterogeneity. Growth curves suggest populations could be in the early stages of exponential growth. However, annual rates of range expansion are predicted to be low (5–16 m). We conclude that extensive nationwide range expansion through secondary introduction is likely to be restricted by currently unsuitable climate beyond southern regions of the UK. However, exponential growth of local populations in habitats providing transport pathways is likely to increase opportunities for regional expansion. The broad habitat niche of P. muralis, coupled with configuration of habitat patches in the landscape, allows populations to increase locally with minimal dispersal.     

Age‐specific habitat preference,carrying capacity,and landscape structure determine the response of population spatial variability to fishing‐driven age truncation

1. Understanding the mechanisms underlying spatial variability of exploited fish is critical for the sustainable management of fish stocks. Empirical studies suggest that size‐selective fishing can elevate fish population spatial variability (i.e., more heterogeneous distribution) through age truncation, making the population less resilient to changing environment. However, species differ in how their spatial variability responds to age truncation and the underlying mechanisms remain unclear.
2. We hypothesize that age‐specific habitat preference, together with environmental carrying capacity and landscape structure, determines the response of population spatial variability to fishing‐induced age truncation. To test these hypotheses, we design an individual‐based model of an age‐structured fish population on a two‐dimensional landscape under size‐selective fishing. Individual fish reproduces and survives, and moves between habitats according to age‐specific habitat preference and density‐dependent habitat selection.
3. Population spatial variability elevates with increasing age truncation, and the response is stronger for populations with stronger age‐specific habitat preference. On a gradient landscape, reducing carrying capacity elevates the relative importance of density dependence in habitat selection, which weakens the response of spatial variability to age truncation for populations with strong age‐specific habitat preference. On a fragmented landscape, both populations with strong and weak age‐specific habitat preferences are restricted at local optimal habitats, and reducing carrying capacity weakens the responses of spatial variability to age truncation for both populations.
4. Synthesis and applications. We demonstrate that to track and predict the changes in population spatial variability under exploitation, it is essential to consider the interactive effects of age‐specific habitat preference, carrying capacity, and landscape structure. To improve spatial management in fisheries, it is crucial to enhance empirical and theoretical developments in the methodology to quantify age‐specific habitat preference of marine fish, and to understand how climatic change influences carrying capacity and landscape continuity.

     

Camera trap reveals the co‐occurrence patterns of two sympatric muntjac species in southern Anhui Province,China: No spatial segregation

The competitive relationship and coexistence pattern among close related species have long been one of the hot issues in ecological research. Interspecies interactions can exert important influences on the local distribution of rare species. Black muntjac Muntiacus crinifrons is an endemic species to eastern China, currently restricted to limited regions. In contrast, Chinese muntjac Muntiacus reevesi is the most common and widespread deer in southern China. Both species co‐occur in southern Anhui and western Zhejiang Province. Little is known about the interaction of these two sympatric‐related species. In this study, to investigate the site use determinants and co‐occurrence pattern of the two sympatric muntjac species, we conducted a camera trap survey across about 250 km² in mountainous area of southern Anhui Province, China. We adopted a multistep approach to incorporate habitat preferences while modeling occupancy and detection. We found that the two species did not separate along elevation gradient (range from 400 m to 1,400 m) as described in previous studies. Results of single‐species occupancy models indicated that elevation had positive effects on the site use of both species, while slope had an opposite influence on their site use. Positive effects of elevation on the site use implied that both species try to avoid human interference at low elevations. Significant negative effect of slope on the site use of black muntjac suggested that the species prefer habitat with gentle slope and avoided steep. Co‐occurrence models and species interaction factors provided evidence that the two muntjac species had an independent occupancy (ψ ^{BM CM} = ψ ^{BM cm}, SIF = 1) and exhibited a positive species interaction in detection probability (p ^BM < r ^{BM CM}). Combined with the results of previous studies, we suggested that it was fine differentiation in microhabitats and food resources utilization rather spatial or temporal segregation that allowed the two species co‐occurrence. The site use determinants revealed in our study would be useful for the habitat conservation and restoration for the rare black muntjac, and the co‐occurrence pattern of the two sympatric muntjac species could provide useful information for deep understanding of the coexistence mechanism among forest‐dwelling ungulates.     

Effects of consumer surface sterilization on diet DNA metabarcoding data of terrestrial invertebrates in natural environments and feeding trials

1. DNA metabarcoding is an emerging tool used to quantify diet in environments and consumer groups where traditional approaches are unviable, including small‐bodied invertebrate taxa. However, metabarcoding of small taxa often requires DNA extraction from full body parts (without dissection), and it is unclear whether surface contamination from body parts alters presumed diet presence or diversity.
2. We examined four different measures of diet (presence, rarefied read abundance, richness, and species composition) for a terrestrial invertebrate consumer (the spider Heteropoda venatoria) both collected in its natural environment and fed an offered diet item in contained feeding trials using DNA metabarcoding of full body parts (opisthosomas). We compared diet from consumer individuals surface sterilized to remove contaminants in 10% commercial bleach solution followed by deionized water with a set of unsterilized individuals.
3. We found that surface sterilization did not significantly alter any measure of diet for consumers in either a natural environment or feeding trials. The best‐fitting model predicting diet detection in feeding trial consumers included surface sterilization, but this term was not statistically significant (β = −2.3, p‐value = .07).
4. Our results suggest that surface contamination does not seem to be a significant concern in this DNA diet metabarcoding study for consumers in either a natural terrestrial environment or feeding trials. As the field of diet DNA metabarcoding continues to progress into new environmental contexts with various molecular approaches, we suggest ongoing context‐specific consideration of the possibility of surface contamination.

     

Population dynamics of little brown bats (Myotis lucifugus) at summer roosts: Apparent survival,fidelity, abundance,and the influence of winter conditions

1. White‐nose syndrome (WNS) has caused the death of millions of bats, but the impacts have been more difficult to identify in western North America. Understanding how WNS, or other threats, impacts western bats may require monitoring other roosts, such as maternity roosts and night roosts, where bats aggregate in large numbers.
2. Little brown bats (Myotis lucifugus) are experiencing some of the greatest declines from WNS. Estimating survival and understanding population dynamics can provide valuable data for assessing population declines and informing conservation efforts.
3. We conducted a 5‐year mark–recapture study of two M. lucifugus roosts in Colorado. We used the robust design model to estimate apparent survival, fidelity, and abundance to understand population dynamics, and environmental covariates to understand how summer and winter weather conditions impact adult female survival. We compared the fidelity and capture probability of M. lucifugus between colonies to understand how bats use such roosts.
4. Overwinter survival increased with the number of days with temperatures below freezing (β > 0.100, SE = 0.003) and decreased with the number of days with snow cover (β < −0.40, SE < 0.13). Adult female fidelity was higher at one maternity roost than the other. Overwinter and oversummer adult female survival was high (>0.90), and based on survival estimates and fungal‐swabbing results, we believe these populations have yet to experience WNS.
5. Recapture of M. lucifugus using antennas that continuously read passive integrated transponder tags allows rigorous estimation of bat population parameters that can elucidate trends in abundance and changes in survival. Monitoring populations at summer roosts can provide unique population ecology data that monitoring hibernacula alone may not. Because few adult males are captured at maternity colonies, and juvenile males have low fidelity, additional effort should focus on understanding male M. lucifugus population dynamics.

     

Interacting forces of predation and fishing affect species’ maturation size

1. Fishing is a strong selective force and is supposed to select for earlier maturation at smaller body size. However, the extent to which fishing‐induced evolution is shaping ecosystems remains debated. This is in part because it is challenging to disentangle fishing from other selective forces (e.g., size‐structured predation and cannibalism) in complex ecosystems undergoing rapid change.
2. Changes in maturation size from fishing and predation have previously been explored with multi‐species physiologically structured models but assumed separation of ecological and evolutionary timescales. To assess the eco‐evolutionary impact of fishing and predation at the same timescale, we developed a stochastic physiologically size‐structured food‐web model, where new phenotypes are introduced randomly through time enabling dynamic simulation of species'' relative maturation sizes under different types of selection pressures.
3. Using the model, we carried out a fully factorial in silico experiment to assess how maturation size would change in the absence and presence of both fishing and predation (including cannibalism). We carried out ten replicate stochastic simulations exposed to all combinations of fishing and predation in a model community of nine interacting fish species ranging in their maximum sizes from 10 g to 100 kg. We visualized and statistically analyzed the results using linear models.
4. The effects of fishing on maturation size depended on whether or not predation was enabled and differed substantially across species. Fishing consistently reduced the maturation sizes of two largest species whether or not predation was enabled and this decrease was seen even at low fishing intensities (F = 0.2 per year). In contrast, the maturation sizes of the three smallest species evolved to become smaller through time but this happened regardless of the levels of predation or fishing. For the four medium‐size species, the effect of fishing was highly variable with more species showing significant and larger fishing effects in the presence of predation.
5. Ultimately our results suggest that the interactive effects of predation and fishing can have marked effects on species'' maturation sizes, but that, at least for the largest species, predation does not counterbalance the evolutionary effect of fishing. Our model also produced relative maturation sizes that are broadly consistent with empirical estimates for many fish species.

     

Context‐dependent outcomes in a reproductive mutualism between two freshwater fish species

1. The development of encompassing general models of ecology is precluded by underrepresentation of certain taxa and systems. Models predicting context‐dependent outcomes of biotic interactions have been tested using plants and bacteria, but their applicability to higher taxa is largely unknown.
2. We examined context dependency in a reproductive mutualism between two stream fish species: mound nest‐building bluehead chub Nocomis leptocephalus and mountain redbelly dace Chrosomus oreas, which often uses N. leptocephalus nests for spawning. We hypothesized that increased predator density and decreased substrate availability would increase the propensity of C. oreas to associate with N. leptocephalus and decrease reproductive success of both species.
3. In a large‐scale in situ experiment, we manipulated egg predator density and presence of both symbionts (biotic context), and replicated the experiment in habitats containing high‐ and low‐quality spawning substrate (abiotic context).
4. Contradictory to our first hypothesis, we observed that C. oreas did not spawn without its host. The interaction outcome switched from commensalistic to mutualistic with changing abiotic and biotic contexts, although the net outcome was mutualistic.
5. The results of this study yielded novel insight into how context dependency operates in vertebrate mutualisms. Although the dilution effect provided by C. oreas positively influenced reproductive success of N. leptocephalus, it was not enough to overcome both egg predation and poor spawning habitat quality. Outcomes of the interaction may be ultimately determined by associate density. Studies of context dependency in vertebrate systems require detailed knowledge of species life‐history traits.