How reliable is metabarcoding for pollen identification? An evaluation of different taxonomic assignment strategies by cross-validation.

Metabarcoding is a powerful tool, increasingly used in many disciplines of environmental sciences. However, to assign a taxon to a DNA sequence, bioinformaticians need to choose between different strategies or parameter values and these choices sometimes seem rather arbitrary. In this work, we present a case study on ITS2 and rbcL databases used to identify pollen collected by bees in Belgium. We blasted a random sample of sequences from the reference database against the remainder of the database using different strategies and compared the known taxonomy with the predicted one. This in silico cross-validation (CV) approach proved to be an easy yet powerful way to (1) assess the relative accuracy of taxonomic predictions, (2) define rules to discard dubious taxonomic assignments and (3) provide a more objective basis to choose the best strategy. We obtained the best results with the best blast hit (best bit score) rather than by selecting the majority taxon from the top 10 hits. The predictions were further improved by favouring the most frequent taxon among those with tied best bit scores. We obtained better results with databases containing the full sequences available on NCBI rather than restricting the sequences to the region amplified by the primers chosen in our study. Leaked CV showed that when the true sequence is present in the database, blast might still struggle to match the right taxon at the species level, particularly with rbcL. Classical 10-fold CV-where the true sequence is removed from the database-offers a different yet more realistic view of the true error rates. Taxonomic predictions with this approach worked well up to the genus level, particularly for ITS2 (5-7% of errors). Using a database containing only the local flora of Belgium did not improve the predictions up to the genus level for local species and made them worse for foreign species. At the species level, using a database containing exclusively local species improved the predictions for local species by ∼12% but the error rate remained rather high: 25% for ITS2 and 42% for rbcL. Foreign species performed worse even when using a world database (59-79% of errors). We used classification trees and GLMs to model the % of errors vs. identity and consensus scores and determine appropriate thresholds below which the taxonomic assignment should be discarded. This resulted in a significant reduction in prediction errors, but at the cost of a much higher proportion of unassigned sequences. Despite this stringent filtering, at least 1/5 sequences deemed suitable for species-level identification ultimately proved to be misidentified. An examination of the variability in prediction accuracy between plant families showed that rbcL outperformed ITS2 for only two of the 27 families examined, and that the % correct species-level assignments were much better for some families (e.g. 95% for Sapindaceae) than for others (e.g. 35% for Salicaceae).

Do pheromone traps help to reduce new attacks of Ips typographus at the local scale after a sanitary cut?

The spruce bark beetle, Ips typographus, is causing severe economic losses during epidemic phases triggered by droughts and/or windstorms. Sanitation felling and salvage logging are usually the most recommended strategies to limit the damages. However, any additional control method to limit the economic impact of an outbreak would be welcome. In this respect, the efficiency of pheromone trapping is still controversial or poorly documented. In this 2-year study (2020-2021), at the peak of a severe outbreak in Belgium, we quantified the wood volume and presence/absence of new attacks at 126 sites attacked during the previous year and within 100 m from the initial attack. Each site was randomly allocated to one of three treatments: (1) three crosstraps baited with pheromones, (2) one tree-trap baited with pheromones and treated with an insecticide and (3) control sites with no trapping device. The attacked trees of the previous year were all cut and removed before the start of the experiment and newly attacked trees were removed as they were detected. The trapping devices were only active during spring to target overwintering bark beetles that might have escaped the sanitation cuts and to limit the risk of attracting dispersing beetles from outside the patch during the summer. We found a strong decrease of the attacks relative to the previous year in all treatments, including the controls (more than 50% of the control sites had no new attacks). There was no relationship between the new attacks and the attacks of the previous year. In both years, new attacks were more frequent (presence/absence) in sites with crosstraps (95% Confidence Interval [56-84%] of the sites with new attacks) than in sites with a tree-trap (26-57% - p = 0.02) and to a lesser extent than in control sites (32-63%, p = 0.08). In 2020, the attacked volumes were slightly higher in sites with crosstraps (95% Confidence Interval [3.4-14.2 m³]) than in control sites (0.2-3.5 m³, p = 0.04) and no significant difference was found with tree-trap sites (1.1-6.2 m³, p = 0.38). In 2021, there were no significant differences between the volumes attacked in the control sites (1.8-9.4 m³), crosstraps sites (0.9-6.4 m³) and tree-trap sites (0-2.5 m³). Overall, we found no evidence in favor of the efficacy of pheromone trapping during spring to reduce economic damages at the local scale when combined with sanitation felling and during a severe outbreak. The use of baited crosstraps could even be hazardous as it seemed to increase the occurrence of new attacks probably by attracting bark beetles but failing to neutralize them.

Mosaic Evolution of Molecular Pathways for Sex Pheromone Communication in a Butterfly.

Unraveling the origin of molecular pathways underlying the evolution of adaptive traits is essential for understanding how new lineages emerge, including the relative contribution of conserved ancestral traits and newly evolved derived traits. Here, we investigated the evolutionary divergence of sex pheromone communication from moths (mostly nocturnal) to butterflies (mostly diurnal) that occurred ~119 million years ago. In moths, it is the females that typically emit pheromones to attract male mates, but in butterflies males emit pheromones that are used by females for mate choice. The molecular bases of sex pheromone communication are well understood in moths, but they have remained relatively unexplored in butterflies. We used a combination of transcriptomics, real time qPCR, and phylogenetics to identify genes involved in the different steps (i.e., production, regulation, and reception) of sex pheromone communication of the butterfly Bicyclus anynana. Our results show that the biosynthesis and reception of sex pheromones relies both on moth-specific gene families (reductases) and on more ancestral insect gene families (desaturases, olfactory receptors, odorant binding proteins). Interestingly, B. anynana appears to use what was believed to be the moth-specific neuropeptide Pheromone Biosynthesis Activating Neuropeptide (PBAN) for regulating sex pheromone production. Altogether, our results suggest that a mosaic pattern best explains how sex pheromone communication evolved in butterflies, with some molecular components derived from moths, and others conserved from more ancient insect ancestors. This is the first large-scale investigation of the genetic pathways underlying sex pheromone communication in a butterfly.

Comparison of qPCR and Metabarcoding Methods as Tools for the Detection of Airborne Inoculum of Forest Fungal Pathogens.

Forest diseases caused by invasive fungal pathogens are becoming more common, sometimes with dramatic consequences to forest ecosystems. The development of early detection systems is necessary for efficient surveillance and to mitigate the impact of invasive pathogens. Windborne spores are an important pathway for introduction of fungal pathogens into new areas; the design of spore trapping devices adapted to forests, capable of collecting different types of spores, and aligned with development of efficient molecular methods for detection of the pathogen, should help forest managers anticipate new disease outbreaks. Two types of Rotorod samplers were evaluated for the collection of airborne inoculum of forest fungal pathogens with a range of spore sizes in five forest types. Detection was by specific quantitative PCR (qPCR) and by high-throughput sequencing (HTS) of amplified internal transcribed spacer sequences using a new bioinformatic pipeline, FungiSearch, developed for diagnostic purposes. Validation of the pipeline was conducted on mock communities of 10 fungal species belonging to different taxa. Although the sensitivity of the new HTS pipeline was lower than the specific qPCR, it was able to detect a wide variety of fungal pathogens. FungiSearch is easy to use, and the reference database is updatable, making the tool suitable for rapid identification of new pathogens. This new approach combining spore trapping and HTS detection is promising as a diagnostic tool for invasive fungal pathogens.

Sexual selection contributes to partial restoration of phenotypic robustness in a butterfly.

Phenotypic variation is the raw material for selection that is ubiquitous for most traits in natural populations, yet the processes underlying phenotypic evolution or stasis often remain unclear. Here, we report phenotypic evolution in a mutant line of the butterfly Bicyclus anynana after outcrossing with the genetically polymorphic wild type population. The comet mutation modifies two phenotypic traits known to be under sexual selection in this butterfly: the dorsal forewing eyespots and the pheromone-producing structures. The original comet mutant line was inbred and remained phenotypically stable for at least seven years, but when outcrossed to the wild type population the outcrossed comet line surprisingly recovered the wild type phenotype within 8 generations at high (27 °C), but not at low (20 °C), developmental temperatures. Male mating success experiments then revealed that outcrossed comet males with the typical comet phenotype suffered from lower mating success, while mating success of outcrossed comet males resembling wild types was partially restored. We document a fortuitous case where the addition of genetic polymorphism around a spontaneous mutation could have allowed partial restoration of phenotypic robustness. We further argue that sexual selection through mate choice is likely the driving force leading to phenotypic robustness in our system.

Time-to-death approach to reveal chronic and cumulative toxicity of a fungicide for honeybees not revealed with the standard ten-day test.

Synthetic fungicides are pesticides widely used in agriculture to control phytopathogenic fungi. The systemicity, persistency and intense application of some of these fungicides, such as boscalid, leads to long periods of exposure for honeybees via contaminated water, pollen and nectar. We exposed adult honeybees in the lab to food contaminated with boscalid for 33 days instead of the standard 10-day test. Most of the toxic effects were observed after 10 days. The median time to death (LT50) ranged from 24.9 days (lowest concentration) to 7.1 days (highest concentration) and was significantly shorter in all cases than with the control (32.0 days). The concentration and dietary doses of boscalid inducing 50% mortality (LC50 and LDD50, respectively) decreased strongly with the time of exposure: LC50 = 14,729 and 1,174 mg/l and LDD50 = 0.318 and 0.0301 mg bee-1 day-1 at days 8 and 25, respectively. We found evidence of reinforced toxicity when exposure is prolonged, but with an unusual pattern: no cumulative toxicity is observed until 17-18 days, when a point of inflexion appears that suggests a reduced capacity of bees to deal with the toxicant. Our results show the importance of time-to-death experiments rather than fixed-duration studies for evaluating chronic toxicity.

The challenges of predicting pesticide exposure of honey bees at landscape level.

To evaluate the risks of pesticides for pollinators, we must not only evaluate their toxicity but also understand how pollinators are exposed to these xenobiotics in the field. We focused on this last point and modeled honey bee exposure to pesticides at the landscape level. Pollen pellet samples (n = 60) from 40 Belgian apiaries were collected from late July to October 2011 and underwent palynological and pesticide residue analyses. Areas of various crops around each apiary were measured at 4 spatial scales. The most frequently detected pesticides were the fungicides boscalid (n = 19, 31.7%) and pyrimethanil (n = 10, 16.7%) and the insecticide dimethoate (n = 10, 16.7%). We were able to predict exposure probability for boscalid and dimethoate by using broad indicators of cropping intensity, but it remained difficult to identify the precise source of contamination (e.g. specific crops in which the use of the pesticide is authorized). For pyrimethanil, we were not able to build any convincing landscape model that could explain the contamination. Our results, combined with the late sampling period, strongly suggest that pesticides applied to crops unattractive to pollinators, and therefore considered of no risk for them, may be sources of exposure through weeds, drift to neighboring plants, or succeeding crops.

Nest Predation Deviates from Nest Predator Abundance in an Ecologically Trapped Bird.

In human-modified environments, ecological traps may result from a preference for low-quality habitat where survival or reproductive success is lower than in high-quality habitat. It has often been shown that low reproductive success for birds in preferred habitat types was due to higher nest predator abundance. However, between-habitat differences in nest predation may only weakly correlate with differences in nest predator abundance. An ecological trap is at work in a farmland bird (Lanius collurio) that recently expanded its breeding habitat into open areas in plantation forests. This passerine bird shows a strong preference for forest habitat, but it has a higher nest success in farmland. We tested whether higher abundance of nest predators in the preferred habitat or, alternatively, a decoupling of nest predator abundance and nest predation explained this observed pattern of maladaptive habitat selection. More than 90% of brood failures were attributed to nest predation. Nest predator abundance was more than 50% higher in farmland, but nest predation was 17% higher in forest. Differences between nest predation on actual shrike nests and on artificial nests suggested that parent shrikes may facilitate nest disclosure for predators in forest more than they do in farmland. The level of caution by parent shrikes when visiting their nest during a simulated nest predator intrusion was the same in the two habitats, but nest concealment was considerably lower in forest, which contributes to explaining the higher nest predation in this habitat. We conclude that a decoupling of nest predator abundance and nest predation may create ecological traps in human-modified environments.

Honeybee colony disorder in crop areas: the role of pesticides and viruses.

As in many other locations in the world, honeybee colony losses and disorders have increased in Belgium. Some of the symptoms observed rest unspecific and their causes remain unknown. The present study aims to determine the role of both pesticide exposure and virus load on the appraisal of unexplained honeybee colony disorders in field conditions. From July 2011 to May 2012, 330 colonies were monitored. Honeybees, wax, beebread and honey samples were collected. Morbidity and mortality information provided by beekeepers, colony clinical visits and availability of analytical matrix were used to form 2 groups: healthy colonies and colonies with disorders (n = 29, n = 25, respectively). Disorders included: (1) dead colonies or colonies in which part of the colony appeared dead, or had disappeared; (2) weak colonies; (3) queen loss; (4) problems linked to brood and not related to any known disease. Five common viruses and 99 pesticides (41 fungicides, 39 insecticides and synergist, 14 herbicides, 5 acaricides and metabolites) were quantified in the samples.The main symptoms observed in the group with disorders are linked to brood and queens. The viruses most frequently found are Black Queen Cell Virus, Sac Brood Virus, Deformed Wing Virus. No significant difference in virus load was observed between the two groups. Three acaricides, 5 insecticides and 13 fungicides were detected in the analysed samples. A significant correlation was found between the presence of fungicide residues and honeybee colony disorders. A significant positive link could also be established between the observation of disorder and the abundance of crop surface around the beehive. According to our results, the role of fungicides as a potential stressor for honeybee colonies should be further studied, either by their direct and/or indirect impacts on bees and bee colonies.

Heritability and artificial selection on ambulatory dispersal distance in Tetranychus urticae: effects of density and maternal effects.

Dispersal distance is understudied although the evolution of dispersal distance affects the distribution of genetic diversity through space. Using the two-spotted spider mite, Tetranychus urticae, we tested the conditions under which dispersal distance could evolve. To this aim, we performed artificial selection based on dispersal distance by choosing 40 individuals (out of 150) that settled furthest from the home patch (high dispersal, HDIS) and 40 individuals that remained close to the home patch (low dispersal, LDIS) with three replicates per treatment. We did not observe a response to selection nor a difference between treatments in life-history traits (fecundity, survival, longevity, and sex-ratio) after ten generations of selection. However, we show that heritability for dispersal distance depends on density. Heritability for dispersal distance was low and non-significant when using the same density as the artificial selection experiments while heritability becomes significant at a lower density. Furthermore, we show that maternal effects may have influenced the dispersal behaviour of the mites. Our results suggest primarily that selection did not work because high density and maternal effects induced phenotypic plasticity for dispersal distance. Density and maternal effects may affect the evolution of dispersal distance and should be incorporated into future theoretical and empirical studies.

Maladaptive habitat selection of a migratory passerine bird in a human-modified landscape.

In human-altered environments, organisms may preferentially settle in poor-quality habitats where fitness returns are lower relative to available higher-quality habitats. Such ecological trapping is due to a mismatch between the cues used during habitat selection and the habitat quality. Maladaptive settlement decisions may occur when organisms are time-constrained and have to rapidly evaluate habitat quality based on incomplete knowledge of the resources and conditions that will be available later in the season. During a three-year study, we examined settlement decision-making in the long-distance migratory, open-habitat bird, the Red-backed shrike (Lanius collurio), as a response to recent land-use changes. In Northwest Europe, the shrikes typically breed in open areas under a management regime of extensive farming. In recent decades, Spruce forests have been increasingly managed with large-size cutblocks in even-aged plantations, thereby producing early-successional vegetation areas that are also colonised by the species. Farmland and open areas in forests create mosaics of two different types of habitats that are now occupied by the shrikes. We examined redundant measures of habitat preference (order of settlement after migration and distribution of dominant individuals) and several reproductive performance parameters in both habitat types to investigate whether habitat preference is in line with habitat quality. Territorial males exhibited a clear preference for the recently created open areas in forests with higher-quality males settling in this habitat type earlier. Reproductive performance was, however, higher in farmland, with higher nest success, offspring quantity, and quality compared to open areas in forests. The results showed strong among-year consistency and we can therefore exclude a transient situation. This study demonstrates a case of maladaptive habitat selection in a farmland bird expanding its breeding range to human-created open habitats in plantations. We discuss the reasons that could explain this decision-making and the possible consequences for the population dynamics and persistence.