A real-world implementation of a nationwide, long-term monitoring program to assess the impact of agrochemicals and agricultural practices on biodiversity

Biodiversity has undergone a major decline throughout recent decades, particularly in farmland. Agricultural practices are recognized to be an important pressure on farmland biodiversity, and pesticides are suspected to be one of the main causes of this decline in biodiversity. As part of the national plan for reduction of pesticides use (Ecophyto), the French ministry of agriculture launched the 500 ENI (nonintended effects) monitoring program in 2012 in order to assess the unintended effects of agricultural practices, including pesticide use, on biodiversity represented by several taxonomic groups of interest for farmers. This long-term program monitors the biodiversity of nontargeted species (earthworms, plants, coleoptera, and birds), together with a wide range of annual data on agricultural practices (crop rotation, soil tillage, weed control, fertilizers, chemical treatments, etc.). Other parameters (e.g., landscape and climatic characteristics) are also integrated as covariates during the analyses. This monitoring program is expected to improve our understanding of the relative contribution of the different drivers of population and community trends. Here, we present the experience of setting up the 500 ENI network for this ambitious and highly complex monitoring program, as well as the type of data it collects. The issue of data quality control and some first results are discussed. With the aim of being useful to readers who would like to set up similar monitoring schemes, we also address some questions that have arisen following the first five years of the implementation phase of the program.

Building a cluster of NLR genes conferring resistance to pests and pathogens: the story of the Vat gene cluster in cucurbits

Most molecularly characterized plant resistance genes (R genes) belong to the nucleotide-binding-site-leucine-rich-repeat (NLR) receptor family and are prone to duplication and transposition with high sequence diversity. In this family, the Vat gene in melon is one of the few R genes known for conferring resistance to insect, i.e., Aphis gossypii, but it has been misassembled and/or mispredicted in the whole genomes of Cucurbits. We examined 14 genomic regions (about 400 kb) derived from long-read assemblies spanning Vat-related genes in Cucumis melo, Cucumis sativus, Citrullus lanatus, Benincasa hispida, Cucurbita argyrosperma, and Momordica charantia. We built the phylogeny of those genes. Investigating the paleohistory of the Vat gene cluster, we revealed a step by step process beginning from a common ancestry in cucurbits older than 50 my. We highlighted Vat exclusively in the Cucumis genera, which diverged about 20 my ago. We then focused on melon, evaluating a minimum duplication rate of Vat in 80 wild and cultivated melon lines using generalist primers; our results suggested that duplication started before melon domestication. The phylogeny of 44 Vat-CDS obtained from 21 melon lines revealed gain and loss of leucine-rich-repeat domains along diversification. Altogether, we revealed the high putative recognition scale offered in melon based on a combination of SNPs, number of leucine-rich-repeat domains within each homolog and number of homologs within each cluster that might jointly confer resistance to a large pest and pathogen spectrum. Based on our findings, we propose possible avenues for breeding programs.

Spatial Distribution Patterns of Parthenolecanium corni (Hemiptera, Coccidae) and of the Ampelovirus GLRaV-1 and the Vitivirus GVA in a Commercial Vineyard

Distribution patterns of the European fruit lecanium Parthenolecanium corni (Bouché) and of grapevine leafroll-associated virus-1 (GLRaV-1) and grapevine virus A (GVA) were monitored from 2003 to 2015 in a Riesling vine plot in the northeast of France. Virus spread was compared between two periods: 2003-2008 and 2009-2014. The percentage of infected vines increased from 54 to 78% for GLRaV-1 and from 14 to 26% for GVA. The spatial distribution of viruses and of P. corni was analysed using permutation tests and revealed an aggregative pattern. Virus distribution was not associated with the density of P. corni population on grapevines. However, GLRaV-1 and GVA spread mainly from initially infected vines. New GLRaV-1 and GVA infections were more frequent on vines near primarily infected vines, first anisotropically along the row, then between neighbouring rows. Virus spread was similar to those described in literature with grapevine mealybug species. This slow vine-to-vine progression suggests that P. corni was responsible for the virus spread, in accordance with the low mobility and low transmission capacities of its local population.

Domestication rewired gene expression and nucleotide diversity patterns in tomato

Plant domestication has led to considerable phenotypic modifications from wild species to modern varieties. However, although changes in key traits have been well documented, less is known about the underlying molecular mechanisms, such as the reduction of molecular diversity or global gene co-expression patterns. In this study, we used a combination of gene expression and population genetics in wild and crop tomato to decipher the footprints of domestication. We found a set of 1729 differentially expressed genes (DEG) between the two genetic groups, belonging to 17 clusters of co-expressed DEG, suggesting that domestication affected not only individual genes but also regulatory networks. Five co-expression clusters were enriched in functional terms involving carbohydrate metabolism or epigenetic regulation of gene expression. We detected differences in nucleotide diversity between the crop and wild groups specific to DEG. Our study provides an extensive profiling of the rewiring of gene co-expression induced by the domestication syndrome in one of the main crop species.

Spatial and Temporal Pattern Analyses of Esca Grapevine Disease in Vineyards in France

To assess the capacity of esca to spread within vineyards of the Bordeaux region, over 8 years of annual records, containing between 1,200 and 2,300 contiguous Cabernet Sauvignon vines from 15 mature vineyards, were used for spatial statistical analyses. A group of nonparametric tests, based on join count statistics and on permutation methods, was developed to characterize the spatial structure of esca-symptomatic vines in terms of spread in any direction or within-row only. Among vineyards, a large range of spatial patterns, from random to strongly structured, associated with various prevalence rates that increased over time were observed. In four vineyards, the complex esca distribution pattern indicated different levels of clustering. By contrast, in other vineyards, only small clusters of two adjacent symptomatic vines were observed, and they were localized along rows, without enlargement over time, except in one vineyard. An analysis of spatial dependence between previously and newly symptomatic vines within k-order neighborhoods (k = 1 to 5), showed, for 5 of the 15 vineyards, that the newly symptomatic vines were located close to previously infected vines, without a favored orientation or neighbor order. All the results together suggested a limited potential for secondary local spread from neighboring symptomatic vines.

The spatial distribution of Mustelidae in France

We estimated the spatial distribution of 6 Mustelidae species in France using the data collected by the French national hunting and wildlife agency under the “small carnivorous species logbooks” program. The 1500 national wildlife protection officers working for this agency spend 80% of their working time traveling in the spatial area in which they have authority. During their travels, they occasionally detect dead or living small and medium size carnivorous animals. Between 2002 and 2005, each car operated by this agency was equipped with a logbook in which officers recorded information about the detected animals (species, location, dead or alive, date). Thus, more than 30000 dead or living animals were detected during the study period. Because a large number of detected animals in a region could have been the result of a high sampling pressure there, we modeled the number of detected animals as a function of the sampling effort to allow for unbiased estimation of the species density. For dead animals -- mostly roadkill -- we supposed that the effort in a given region was proportional to the distance traveled by the officers. For living animals, we had no way to measure the sampling effort. We demonstrated that it was possible to use the whole dataset (dead and living animals) to estimate the following: (i) the relative density -- i.e., the density multiplied by an unknown constant -- of each species of interest across the different French agricultural regions, (ii) the sampling effort for living animals for each region, and (iii) the relative detection probability for various species of interest.

Highlighting features of spatiotemporal spread of powdery mildew epidemics in the vineyard using statistical modeling on field experimental data

A greater understanding of the development of powdery mildew epidemics on vines would improve disease management by making assessments of the risk of invasion more accurate. We characterized the spatiotemporal spread of epidemics in the vineyard, quantified their variability, and identified the factors responsible for it. We described changes in the probability of infection of a leaf in a plot over time and as a function of distance from a source of disease. Logistic models were fitted to field data from artificially inoculated plots. The velocity of spread decreased along the row and increased in the direction of the prevailing winds. The rate of progression over time was plot dependent, and the velocity was dependent on the vigor of the vine (0.1 to 0.27 m day(-1) in areas of moderate vigor and 1.1 m day(-1) in areas of high vigor). When applied to a larger plot with natural primary foci, the spatiotemporal logistic model showed that the velocity and the slope of the gradient in space depended on the foci; however, the velocity remained in the same range. During the period of highest susceptibility for grape, the probability of a leaf becoming infected increased from 2.5 to 13%. Our logistic model was able to predict changes in disease over time of its extension within the plot; however, the crop heterogeneity prevented prediction of variability of disease at the vine scale.

Estimation of abomasum strongyle nematode infections in sheep at necropsy: Tentative proposals for a simplified technique

Several necropsy techniques are available for estimating the abundance of gastro-intestinal nematodes in abomasum of ruminants. Standardization of techniques is needed to allow accurate comparisons between laboratories. Here we propose a standardized technique for estimating the abundance of worms. We intend to compare the worms' number estimations in lambs and ewes based on contents and washings, to determine the uniformity of worm counts in aliquots, and to estimate the total worm number from washings. The digesta (or "contents") and the washings of the abomasum are treated separately. The worms of each subsample are diluted with water and the total number of worms is estimated on a small volume (aliquots) of these subsamples. The use of aliquots assumes that the worms are uniformly distributed in the whole volume of each subsample. We first confirmed that the use of aliquots is appropriate in most cases. We then show that the use of the washings alone allows a faster and a suitable estimation of the total worm burden for all strongyle species of the abomasum in both ewes and lambs. The evaluation of our necropsy procedure is a first step to a standardized technique which should be improved by validation in other laboratories.

Different mutations in the genome-linked protein VPg of potato virus Y confer virulence on the pvr2(3) resistance in pepper

Five different amino acid substitutions in the VPg of Potato virus Y were shown to be independently responsible for virulence toward pvr2(3) resistance gene of pepper. A consequence of these multiple mutations toward virulence involving single nucleotide substitutions is a particularly high frequency of resistance breaking (37% of inoculated plants from the first inoculation) and suggests a potentially low durability of pvr2(3) resistance. These five mutants were observed with significantly different frequencies, one of them being overrepresented. Genetic drift alone could not explain the observed distribution of virulent mutants. More plausible scenarios were obtained by taking into account either the relative substitution rates, the relative fitness of the mutants in pvr2(3) pepper plants, or both.

[An epidemiological survey of equine anaplasmosis (Anaplasma phagocytophilum) in southern France]

Anaplasmosis is caused by the bacterium Anaplasma phagocytophilum and transmitted by Ixodes spp. ticks. According to some reports the disease can be introduced into disease-free zones by migrating birds. The purpose of this study was to evaluate the seroprevalence of A. phagocytophilum in horses in the Camargue. Data concerning 424 horses were gathered and the sera were tested for A. phagocytophilum and for piroplasmoses using an enzyme-linked immunosorbent assay and a complement fixation test, respectively. The seroprevalence rates were 11.3 % for A. phagocytophilum, 64.4 % for Theileria equi and 19.7% for Babesia caballi. Stallions were less likely to produce antibodies against A. phagocytophilumthan were females or geldings (odds ratio [OR] = 0.3; p = 0.021). The presence of swallows increased the risk of infections in stables (OR = 5.18; p = 0.011). Spatial analysis showed the existence of groups of infected stables along canals and rivers (p = 0.008). These results suggest an emergence of A. phagocytophilum in the Camargue.

Prevalence of multidrug resistant (MDR) Salmonella in bovine dairy herds in western France

As a part of our effort in quantitative risk analysis of food-borne diseases, we carried out an epidemiologic study to estimate the prevalence of multidrug resistant (MDR) Salmonella in dairy herds situated in western France. The study population consisted of 489 farms in the region and manure or slurry was sampled from these operations and tested for the Salmonella spp. All strains isolated during the study were serotyped and tested for their antimicrobial susceptibility. Salmonella spp. was isolated from 8.1% (95% confidence interval (CI 95%): 4.5-13.3%) of the sampled herds. The herd prevalence of MDR Salmonella among the sampled herds was 1.9% (CI 95%: 0.5-5.4%). Spatial statistics were used to check for sampling representativeness and to determine if infected herds were clustered spatially.

Modelling the interactions between phenology and insecticide resistance genes in the codling moth Cydia pomonella

In the codling moth Cydia pomonella (L), insecticide resistance genes have been associated with pleiotropic effects affecting phenology. In this paper, we investigated whether an increase in the frequency of insecticide resistance in field populations of C pomonella was likely to entail significant divergences in the temporal occurrence of both susceptible and insecticide-resistant individuals. For this purpose, we built a phenological model that provided suitable predictions of the distinct and diverging seasonal evolutions of populations of a susceptible and two insecticide-resistant (at two and three loci) homozygous genotypes of C pomonella. Model simulations for each genotype were further compared with pheromone trap catches recorded in a field insecticide-treated population over an 8-year period (from 1992 to 2000), which reflected the progressive annual increase in the frequency of resistance in southeastern France. We found a significant delay in field adult emergence relative to those predicted by the homozygous susceptible model, and the magnitude of such a delay was positively correlated with increasing frequencies of insecticide resistance in the sampled field population of C pomonella. Adult emergence predicted in the theoretical population that was homozygous for resistance at two loci converged with those recorded in the field during the investigated 8-year period. This suggested that the pleiotropic effects of resistance were likely to result in a significant phenological segregation of insecticide-resistant alleles in the field. The results of this study emphasized the potential for pest populations exposed to chemical selection to evolve qualitatively with respect to phenology. This may raise critical questions regarding the use of phenological modelling as a forecasting tool for appropriate resistance management strategies that would take into account the diverging seasonal evolutions of both insecticide resistance and susceptibility.

Constraints on adaptive mutations in the codling moth Cydia pomonella (L.): measuring fitness trade-offs and natural selection

Adaptive changes in populations encountering a new environment are often constrained by deleterious pleiotropic interactions with ancestral physiological functions. Evolutionary responses of populations can thus be limited by natural selection under fluctuating environmental conditions, if the adaptive mutations are associated with pleiotropic fitness costs. In this context, we have followed the evolution of the frequencies of insecticide-resistant mutants of Cydia pomonella when reintroduced into an untreated environment. The novel set of selective forces after removal of insecticide pressure led to the decline of the frequencies of resistant phenotypes over time, suggesting that the insecticide-adapted genetic variants were selected against the absence of insecticide (with a selective coefficient estimated at 0.11). The selective coefficients were also estimated for both the major cytochrome P450-dependent monooxygenase (MFO) and the minor glutathione S-transferase (GST) systems (0.17 and negligible, respectively), which have been previously shown to be involved in resistance. The involvement of metabolic systems acting both through xenobiotic detoxification and biosynthetic pathways of endogenous compounds may be central to explaining the deleterious physiological consequences resulting from pleiotropy of adaptive changes. The estimation of the magnitude of the fitness cost associated with insecticide resistance in C. pomonella suggests that resistance management strategies exclusively based on insecticide alternations would be unlikely to delay such a selection process.

A point pattern model of the spread of foot-and-mouth disease

The spatial spread of foot-and-mouth disease (FMD) is influenced by several sources of spatial heterogeneity: heterogeneity of the exposure to the virus, heterogeneity of the animal density and heterogeneity of the networks formed by the contacts between farms. A discrete space model assuming that farms can be reduced to points is proposed to handle these different factors. The farm-to-farm process of transmission of the infection is studied using point-pattern methodology. Farm management, commercial exchanges, possible airborne transmission, etc. cannot be explicitly taken into account because of lack of data. These latter factors are introduced via surrogate variables such as herd size and distance between farms. The model is built on the calculation of an infectious potential for each farm. This method has been applied to the study of the 1967-1968 FMD epidemic in UK and allowed us to evaluate the spatial variation of the probability of infection during this epidemic. Maximum likelihood estimation has been conducted conditional on the absence of data concerning the farms which were not infected during the epidemic. Model parameters have then been tested using an approximated conditional-likelihood ratio test. In this case study, results and validation are limited by the lack of data, but this model can easily be extended to include other information such as the effect of wind direction and velocity on airborne spread of the virus or the complex interactions between the locations of farms and the herd size. It can also be applied to other diseases where point approximation is convenient. In the context of an increase of animal density in some areas, the model explicitly incorporates the density and known epidemiological characteristics (e.g. incubation period) in the calculation of the probability of FMD infection. Control measures such as vaccination or slaughter can be simply introduced, respectively, as a reduction of the susceptible population or as a reduction of the source of infection.

A stochastic flowering model describing an asynchronically flowering set of trees

A general stochastic model is presented that simulates the time course of flowering of individual trees and populations, integrating the synchronization of flowering both between and within trees. Making some hypotheses, a simplified expression of the model, called the 'shoot' model, is proposed, in which the synchronization of flowering both between and within trees is characterized by specific parameters. Two derived models, the 'tree' model and the 'population' model, are presented. They neglect the asynchrony of flowering, respectively, within trees, and between and within trees. Models were fitted and tested using data on flowering of Psidium cattleianum observed at study sites at elevations of 200, 520 and 890 m in Reunion Island. The 'shoot' model fitted the data best and reproduced the strong irregularities in flowering shown by empirical data. The asynchrony of flowering in P. cattleianum was more pronounced within than between trees. Simulations showed that various flowering patterns can be reproduced by the 'shoot' model. The use of different levels of organization of the general model is discussed.

A reference map of Cucumis melo based on two recombinant inbred line populations

A composite genetic melon map was generated based on two recombinant inbred line (RI) populations. By analyzing the segregation of 346 AFLPs, 113 IMAs and phenotypic characters on a RI population of 163 individuals derived from the cross Védrantais x PI 161375, a first map was constructed. About 20% of the molecular markers were skewed, and the residual heterozygosity was estimated at 4.43% which was not significantly different from the theoretical value of 4.2%. The genome distribution of molecular markers among the 12 linkage groups was not different from a random distribution with the exception of linkage group XII which was found significantly less populated. The genome distributions of IMAs and AFLPs were complementary. AFLPs were found mainly in the middle of each linkage group and sometimes clustered, whereas IMAs were found mainly at the end. A total of 318 molecular markers, mainly AFLP and IMA markers, were mapped on 63 RIs of the second population, Védrantais x PI 414723. Comparison of the maps enables one to conclude that AFLPs and IMAs of like molecular size, amplified with the same primer combination, correspond to the same genetic locus. Both maps were joined through 116 common markers comprising 106 comigrating AFLPs/IMAs, plus five SSRs and five phenotypic markers. The integrated melon map contained 668 loci issuing from the segregation of 1,093 molecular markers in the two RI populations. The composite map spanned 1,654 cM on 12 linkage groups which is the haploid number of chromosomes in melon. Thirty two known-function probes, i.e. known-function genes (9) and morphological traits (23), were included in this map. In addition, the composite map was anchored to previously published maps through SSRs, RFLPs and phenotypic characters.

Testing local independence between two point processes

Dependencies between two types of points in a spatial point process can be due either to a real dependence between the two types or to the dependence on common underlying variables. We propose a global test for dependence between two point processes that is valid for a wide range of models. In contrast with previously proposed methods, it is based on a number of local test statistics, which makes it possible to map the local association between the two processes. The behavior of the test is evaluated by a simulation study. It is then applied to a vegetation pattern data set from Burkina Faso.

Testing local dependence of spatial structures on images

Associations between two spatial processes can be due to a real dependence between the two processes or to the dependence on common underlying variables. We propose to test the existence of a real dependence by use of local tests, leading to a global test of real dependence and a map of local interactions. We present first how classical interaction tests based on random rotations between completely observed processes such as those developed by Berman (Berman. Appl. Statist. (1986) 35, 54-62), can be integrated in local analyses. For this purpose, tests are first performed locally, and the distribution of their p-values is then compared to the corresponding value under the null hypothesis. A similar approach is proposed to test non-stationarity of a point pattern by using distance statistics popularized by Diggle (Diggle. Statistical Analysis of Spatial Point Patterns. (1983) Academic Press, New York). The problem of testing the interaction between a random field and a censoring area pattern process is discussed and an approach similar to the preceding ones is then proposed. The methods are mainly applied to agricultural examples but they can be applied to any microscopical images for which one wishes to analyse the spatial structure.

Effect of three ectomycorrhizal fungi on growth and phosphorus uptake of Pinus silvestris seedlings at increasing phosphorus levels

The objective of this study was to determine the different effects of three ectomycorrhizal fungi (Laccaria laccata, Hebeloma crustuliniforme, and Thelephora terrestris) on the growth and mineral nutrition of Pinus silvestris at different levels of soluble phosphorus. Even a low intensity of infection by Laccaria laccata stimulated Pinus silvestris growth greatly. The ability of this ectomycorrhizal fungus to increase Pinus silvestris growth seemed to be more related to its capacity to produce growth substances than to its capacity to stimulate phosphorus uptake. The poor efficiency of Hebeloma crustuliniforme compared with Laccaria laccata at any level of phosphorus could result from differences in diversion of carbohydrates from the host to fungal structures.