Preventing multi-resistance: New insights for managing fungal adaptation.

Sustainable crop protection is vital for food security, yet it is under threat due to the adaptation of a diverse and evolving pathogen population. Resistance can be managed by maximising the diversity of selection pressure through dose variation and the spatial and temporal combination of active ingredients. This study explores the interplay between operational drivers for maximising the sustainability of management strategies in relation to the resistance status of fungal populations. We applied an experimental evolution approach to three artificial populations of Zymoseptoria tritici, an economically significant wheat pathogen, each differing in initial resistance status. Our findings reveal that diversified selection pressure curtails the selection of resistance in naïve populations and those with low frequencies of single resistance. Increasing the number of modes of action most effectively delays resistance development, surpassing the increase in the number of fungicides, fungicide choice based on resistance risk, and temporal variation in fungicide exposure. However, this approach favours generalism in the evolved populations. The prior presence of multiple resistant isolates and their subsequent selection in populations override the effects of diversity in management strategies, thereby invalidating any universal ranking. Therefore, the initial resistance composition must be specifically considered in sustainable resistance management to address real-world field situations.

Large-scale study validates that regional fungicide applications are major determinants of resistance evolution in the wheat pathogen Zymoseptoria tritici in France.

In modern cropping systems, the near-universal use of plant protection products selects for resistance in pest populations. The emergence and evolution of this adaptive trait threaten treatment efficacy. We identified determinants of fungicide resistance evolution and quantified their effects at a large spatiotemporal scale. We focused on Zymoseptoria tritici, which causes leaf blotch in wheat. Phenotypes of qualitative or quantitative resistance to various fungicides were monitored annually, from 2004 to 2017, at about 70 sites throughout 22 regions of France (territorial units of 25 000 km2 on average). We modelled changes in resistance frequency with regional anti-Septoria fungicide use, yield losses due to the disease and the regional area under organic wheat. The major driver of resistance dynamics was fungicide use at the regional scale. We estimated its effect on the increase in resistance and relative apparent fitness of each resistance phenotype. The predictions of the model replicated the spatiotemporal patterns of resistance observed in field populations (R2 from 0.56 to 0.82). The evolution of fungicide resistance is mainly determined at the regional scale. This study therefore showed that collective management at the regional scale could effectively complete local actions.

The polyphagous plant pathogenic fungus Botrytis cinerea encompasses host-specialized and generalist populations.

The host plant is often the main variable explaining population structure in fungal plant pathogens, because specialization contributes to reduce gene flow between populations associated with different hosts. Previous population genetic analysis revealed that French populations of the grey mould pathogen Botrytis cinerea were structured by hosts tomato and grapevine, suggesting host specialization in this highly polyphagous pathogen. However, these findings raised questions about the magnitude of this specialization and the possibility of specialization to other hosts. Here we report specialization of B. cinerea populations to tomato and grapevine hosts but not to other tested plants. Population genetic analysis revealed two pathogen clusters associated with tomato and grapevine, while the other clusters co-occurred on hydrangea, strawberry and bramble. Measurements of quantitative pathogenicity were consistent with host specialization of populations found on tomato, and to a lesser extent, populations found on grapevine. Pathogen populations from hydrangea and strawberry appeared to be generalist, while populations from bramble may be weakly specialized. Our results suggest that the polyphagous B. cinerea is more accurately described as a collection of generalist and specialist individuals in populations. This work opens new perspectives for grey mould management, while suggesting spatial optimization of crop organization within agricultural landscapes.

Fashionably late partners have more fruitful encounters: Impact of the timing of co-infection and pathogenicity on sexual reproduction in Zymoseptoria tritici.

The wheat pathogen Zymoseptoria tritici is a relevant fungal model organism for investigations of the epidemiological determinants of sexual reproduction. The objective of this experimental study was to determine which intrinsic factors, including parental fitness and timing conditions of infection, affect the numbers of ascospores produced. We first performed 28 crosses on adult wheat plants in semi-controlled conditions, with 10 isolates characterized for their fitness traits. We validated the efficiency of the crossing method, opening up new perspectives for epidemiological studies. We found that the ability to reproduce sexually was determined, at least partly, by the parental genotypes. We also found that the number of ascospores released was correlated with the mean size of the sporulating lesions of the parental isolates on the one hand, and the absolute difference in the latent periods of these isolates on the other. No functional trade-off between the two modes of reproduction in Z. tritici was revealed: there was no adaptive compromise between pathogenicity (asexual multiplication on leaves) and transmission (intensity of sexual reproduction on wheat debris). Moreover, a few days' difference in the latent periods of the two parental isolates, such that one progressed more rapidly in the host tissue than the other, seemed to be slightly beneficial to ascosporogenesis. This may be because the first parental isolate breaks down host defenses, thereby facilitating infection for the other parental isolate. However, a larger difference (a few weeks), generated by leaving two to three weeks between the inoculations of the plant with the parental isolates, was clearly detrimental to ascosporogenesis. In this case, the host tissues were likely colonized by the first isolate, leaving less host resources available for the second, consistent with a competition effect during the asexual stage.

Approximate Bayesian computation with functional statistics.

Functional statistics are commonly used to characterize spatial patterns in general and spatial genetic structures in population genetics in particular. Such functional statistics also enable the estimation of parameters of spatially explicit (and genetic) models. Recently, Approximate Bayesian Computation (ABC) has been proposed to estimate model parameters from functional statistics. However, applying ABC with functional statistics may be cumbersome because of the high dimension of the set of statistics and the dependences among them. To tackle this difficulty, we propose an ABC procedure which relies on an optimized weighted distance between observed and simulated functional statistics. We applied this procedure to a simple step model, a spatial point process characterized by its pair correlation function and a pollen dispersal model characterized by genetic differentiation as a function of distance. These applications showed how the optimized weighted distance improved estimation accuracy. In the discussion, we consider the application of the proposed ABC procedure to functional statistics characterizing non-spatial processes.

Plant and animal protein intakes are differently associated with nutrient adequacy of the diet of French adults.

The sustainability of dietary plant and animal protein consumption is an important issue, but few data are available to ascertain the health-related impact of animal and plant protein intake. The objective was to study the relationships between animal and plant protein intake and diet quality, as reflected by an integrated index of nutrient adequacy. Using data on 1912 adults from the French Individual and National Consumption Survey 2 (2006-2007), we assessed diet quality by using the PANDiet, a unique score combining 35 probabilities of having an adequate nutrient intake, and used multiple regression analysis to study the relationship with intakes of plant, animal, and related food source proteins. After adjustments for intakes of energy, alcohol, and protein from other sources and sociodemographic factors, plant protein intake was positively associated with the PANDiet, irrespective of sex (ß = 0.50, P < 0.0001). By contrast, total and animal protein intakes were differently associated with the PANDiet according to sex, with a positive association in women (ß = 0.07 and ß = 0.08, respectively; P < 0.01) and an inverse association in men (ß = -0.07 and ß = -0.05, respectively; P < 0.01). The relationships between the PANDiet and intakes of protein from animal food sources varied: for instance, associations were inverse for red meat and poultry in men but not in women, whereas irrespective of sex, they were positive for fish, milk, and yogurt and inverse for processed meat, cheese, and eggs. These findings show that plant protein is a robust marker of a healthy diet, whereas total animal protein includes different subtypes of animal protein that largely vary in their relationship with diet quality. These data could help refine protein recommendations from a qualitative and dietary standpoint.

Antifungal alternation can be beneficial for durability but at the cost of generalist resistance.

The evolution of resistance to pesticides is a major burden in agriculture. Resistance management involves maximizing selection pressure heterogeneity, particularly by combining active ingredients with different modes of action. We tested the hypothesis that alternation may delay the build-up of resistance not only by spreading selection pressure over longer periods, but also by decreasing the rate of evolution of resistance to alternated fungicides, by applying an experimental evolution approach to the economically important crop pathogen Zymoseptoria tritici. Our results show that alternation is either neutral or slows the overall resistance evolution rate, relative to continuous fungicide use, but results in higher levels of generalism in evolved lines. We demonstrate that the nature of the fungicides, and therefore their relative intrinsic risk of resistance may underly this trade-off, more so than the number of fungicides and the rhythm of alternation. This trade-off is also dynamic over the course of resistance evolution. These findings open up new possibilities for tailoring resistance management effectively while optimizing interplay between alternation components.

Siland a R package for estimating the spatial influence of landscape.

The spatial distributions of populations are both influenced by local variables and by characteristics of surrounding landscapes. Understanding how landscape features spatially structure the frequency of a trait in a population, the abundance of a species or the species' richness remains difficult specially because the spatial scale effects of the landscape variables are unknown. Various methods have been proposed but their results are not easily comparable. Here, we introduce "siland", a general method for analyzing the effect of landscape features. Based on a sequential procedure of maximum likelihood estimation, it simultaneously estimates the spatial scales and intensities of landscape variable effects. It does not require any information about the scale of effect. It integrates two landscape effects models: one is based on focal sample site (Bsiland, b for buffer) and one is distance weighted using Spatial Influence Function (Fsiland, f for function). We implemented "siland" in the adaptable and user-friendly R eponym package. It performs landscape analysis on georeferenced point observations (described in a Geographic Information System shapefile format) and allows for effects tests, effects maps and models comparison. We illustrated its use on a real dataset by the study of a crop pest (codling moth densities).

Are Efficient-Dose Mixtures a Solution to Reduce Fungicide Load and Delay Evolution of Resistance? An Experimental Evolutionary Approach.

Pesticide resistance poses a critical threat to agriculture, human health and biodiversity. Mixtures of fungicides are recommended and widely used in resistance management strategies. However, the components of the efficiency of such mixtures remain unclear. We performed an experimental evolutionary study on the fungal pathogen Z. tritici to determine how mixtures managed resistance. We compared the effect of the continuous use of single active ingredients to that of mixtures, at the minimal dose providing full control of the disease, which we refer to as the "efficient" dose. We found that the performance of efficient-dose mixtures against an initially susceptible population depended strongly on the components of the mixture. Such mixtures were either as durable as the best mixture component used alone, or worse than all components used alone. Moreover, efficient dose mixture regimes probably select for generalist resistance profiles as a result of the combination of selection pressures exerted by the various components and their lower doses. Our results indicate that mixtures should not be considered a universal strategy. Experimental evaluations of specificities for the pathogens targeted, their interactions with fungicides and the interactions between fungicides are crucial for the design of sustainable resistance management strategies.

Spatiotemporal dynamics of fungicide resistance in the wheat pathogen Zymoseptoria tritici in France.

BACKGROUND: Management of pesticide resistance is a major issue in modern agricultural systems, particularly in the context of the broader challenge of reducing pesticide use. However, such management must be adapted to resistance dynamics, which remains challenging to predict due to its dependence on many biological traits of pests, interactions with the environment and pesticide use. We retrospectively studied the evolution of reported resistances to four modes of action (benzimidazoles, quinone outside inhibitors, sterol demethylation inhibitors and succinate dehydrogenase inhibitors), in French populations of the wheat pathogen Zymoseptoria tritici. RESULTS: We used statistical models to analyse the Performance trial network data set (2004-2017; ∼ 70 locations in France yearly). They highlighted contrasting behaviours between phenotypes, for example: (i) stable spatial distributions and colonization front structures over time, and (ii) different frequency growth rates at the national scale and between regions. CONCLUSION: We provide a quantitative description of the spatiotemporal patterns of resistance evolution for fungicides with several modes of action. Moreover, we highlight some unexpected resistance dynamics in France, with major differences between the north and south. This complex pattern of resistance evolution in French populations is consistent with previous descriptions of dynamics at the European scale. These results should make it easier to anticipate evolution locally and to improve the management of resistance. © 2019 Society of Chemical Industry.

Epidemiological trade-off between intra- and interannual scales in the evolution of aggressiveness in a local plant pathogen population.

The efficiency of plant resistance to fungal pathogen populations is expected to decrease over time, due to their evolution with an increase in the frequency of virulent or highly aggressive strains. This dynamics may differ depending on the scale investigated (annual or pluriannual), particularly for annual crop pathogens with both sexual and asexual reproduction cycles. We assessed this time-scale effect, by comparing aggressiveness changes in a local Zymoseptoria tritici population over an 8-month cropping season and a 6-year period of wheat monoculture. We collected two pairs of subpopulations to represent the annual and pluriannual scales: from leaf lesions at the beginning and end of a single annual epidemic and from crop debris at the beginning and end of a 6-year period. We assessed two aggressiveness traits-latent period and lesion size-on sympatric and allopatric host varieties. A trend toward decreased latent period concomitant with a significant loss of variability was established during the course of the annual epidemic, but not over the 6-year period. Furthermore, a significant cultivar effect (sympatric vs. allopatric) on the average aggressiveness of the isolates revealed host adaptation, arguing that the observed patterns could result from selection. We thus provide an experimental body of evidence of an epidemiological trade-off between the intra- and interannual scales in the evolution of aggressiveness in a local plant pathogen population. More aggressive isolates were collected from upper leaves, on which disease severity is usually lower than on the lower part of the plants left in the field as crop debris after harvest. We suggest that these isolates play little role in sexual reproduction, due to an Allee effect (difficulty finding mates at low pathogen densities), particularly as the upper parts of the plant are removed from the field, explaining the lack of transmission of increases in aggressiveness between epidemics.

Chikungunya virus-associated long-term arthralgia: a 36-month prospective longitudinal study.

BACKGROUND: Arthritogenic alphaviruses, including Chikungunya virus (CHIKV), are responsible for acute fever and arthralgia, but can also lead to chronic symptoms. In 2006, a Chikungunya outbreak occurred in La Réunion Island, during which we constituted a prospective cohort of viremic patients (n = 180) and defined the clinical and biological features of acute infection. Individuals were followed as part of a longitudinal study to investigate in details the long-term outcome of Chikungunya. METHODOLOGY/PRINCIPAL FINDINGS: Patients were submitted to clinical investigations 4, 6, 14 and 36 months after presentation with acute CHIKV infection. At 36 months, 22 patients with arthralgia and 20 patients without arthralgia were randomly selected from the cohort and consented for blood sampling. During the 3 years following acute infection, 60% of patients had experienced symptoms of arthralgia, with most reporting episodic relapse and recovery periods. Long-term arthralgias were typically polyarthralgia (70%), that were usually symmetrical (90%) and highly incapacitating (77%). They were often associated with local swelling (63%), asthenia (77%) or depression (56%). The age over 35 years and the presence of arthralgia 4 months after the disease onset are risk factors of long-term arthralgia. Patients with long-term arthralgia did not display biological markers typically found in autoimmune or rheumatoid diseases. These data helped define the features of CHIKV-associated chronic arthralgia and permitted an estimation of the economic burden associated with arthralgia. CONCLUSIONS/SIGNIFICANCE: This study demonstrates that chronic arthralgia is a frequent complication of acute Chikungunya disease and suggests that it results from a local rather than systemic inflammation.