Atmosphere-soil carbon transfer as a function of soil depth.

The exchange of carbon between soil organic carbon (SOC) and the atmosphere affects the climate1,2 and-because of the importance of organic matter to soil fertility-agricultural productivity3. The dynamics of topsoil carbon has been relatively well quantified4, but half of the soil carbon is located in deeper soil layers (below 30 centimetres)5-7, and many questions remain regarding the exchange of this deep carbon with the atmosphere8. This knowledge gap restricts soil carbon management policies and limits global carbon models1,9,10. Here we quantify the recent incorporation of atmosphere-derived carbon atoms into whole-soil profiles, through a meta-analysis of changes in stable carbon isotope signatures at 112 grassland, forest and cropland sites, across different climatic zones, from 1965 to 2015. We find, in agreement with previous work5,6, that soil at a depth of 30-100 centimetres beneath the surface (the subsoil) contains on average 47 per cent of the topmost metre's SOC stocks. However, we show that this subsoil accounts for just 19 per cent of the SOC that has been recently incorporated (within the past 50 years) into the topmost metre. Globally, the median depth of recent carbon incorporation into mineral soil is 10 centimetres. Variations in the relative allocation of carbon to deep soil layers are better explained by the aridity index than by mean annual temperature. Land use for crops reduces the incorporation of carbon into the soil surface layer, but not into deeper layers. Our results suggest that SOC dynamics and its responses to climatic control or land use are strongly dependent on soil depth. We propose that using multilayer soil modules in global carbon models, tested with our data, could help to improve our understanding of soil-atmosphere carbon exchange.

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.

Arthropods life cycle and temperature: Beyond isomorphy hypothesis.

In poikilotherms, developmental time and adult lifespan are biological traits driven by temperature although their underlying physiological mechanisms differ. The developmental theory of ageing predicts a positive correlation between these two traits without confirming a genetic or causal relationship. The developmental rate isomorphy hypothesis established that the proportionality in the duration of each pre-imaginal stage with respect to total developmental time does not vary with temperature. This may have important evolutionary implications. We conducted an analysis with arthropods of agronomic interest to study if the hypothesis could be extended to the total lifespan of poikilotherms, including adult lifespan. We showed that isomorphy could be acceptable for a global description of this relationship. However a general model of power law type is more appropriate for characterizing the relationship between developmental time and adult lifespan. In this model, the shape of the curve is a constant characteristic, but the parameters that control it depend on the functional type of the arthropod (phytophagous, predator or parasitoid). Presumably this power law could be extended to all arthropods of agricultural interest, and more generally to all arthropod.

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.

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.

The relationship between mutation frequency and replication strategy in positive-sense single-stranded RNA viruses.

For positive-sense single-stranded RNA virus genomes, there is a trade-off between the mutually exclusive tasks of transcription, translation and encapsidation. The replication strategy that maximizes the intracellular growth rate of the virus requires iterative genome transcription from positive to negative, and back to positive sense. However, RNA viruses experience high mutation rates, and the proportion of genomes with lethal mutations increases with the number of replication cycles. Thus, intracellular mutant frequency will depend on the replication strategy. Introducing apparently realistic mutation rates into a model of viral replication demonstrates that strategies that maximize viral growth rate could result in an average of 26 mutations per genome by the time plausible numbers of positive strands have been generated, and that virus viability could be as low as 0.1 per cent. At high mutation rates or when a high proportion of mutations are deleterious, the optimal strategy shifts towards synthesizing more negative strands per positive strand, and in extremis towards a 'stamping-machine' replication mode where all the encapsidated genomes come from only two transcriptional steps. We conclude that if viral mutation rates are as high as current estimates suggest, either mutation frequency must be considerably higher than generally anticipated and the proportion of viable viruses produced extremely small, or replication strategies cannot be optimized to maximize viral growth rate. Mechanistic models linking mutation frequency to replication mechanisms coupled with data generated through new deep-sequencing technologies could play an important role in improving the estimates of viral mutation rate.

Asymmetrical over-infection as a process of plant virus emergence.

Disentangling the role of epidemiological factors in plant pathogen emergences is a prerequisite to identify the most likely future invaders. An example of emergence was recently observed in France: in 10 years, "classic" (CL) strains of Watermelon mosaic virus (WMV) were displaced at a regional scale by newly introduced "emerging" (EM) strains. Here we analyse a 3 years dataset describing the co-dynamics of CL and EM strains at field scale using state-space models estimating jointly: (i) probabilities of primary and secondary infection and (ii) probabilities of over-infecting with a CL [EM] strain a plant already infected with an EM [CL] strain. Results especially indicate that it is more than 3 times less probable for a CL strain to over-infect an EM infected plant than for an EM strain to over-infect a CL infected plant. To investigate if these asymmetric interactions can explain the CL/EM shift observed at regional scale, an exploratory model describing WMV epidemiology over several years in a landscape composed of a reservoir and a cultivated compartment is introduced. In most simulations a shift is observed and both strains do coexist in the landscape, reaching an equilibrium that depends on the probabilities of over-infection.

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.

Ecological Specialization and Rarity of Arable Weeds: Insights from a Comprehensive Survey in France.

: The definition of "arable weeds" remains contentious. Although much attention has been devoted to specialized, segetal weeds, many taxa found in arable fields also commonly occur in other habitats. The extent to which adjacent habitats are favorable to the weed flora and act as potential sources of colonizers in arable fields remains unclear. In addition, weeds form assemblages with large spatiotemporal variability, so that many taxa in weed flora are rarely observed in plot-based surveys. We thus addressed the following questions: How often do weeds occur in other habitats than arable fields? How does including field edges extend the taxonomic and ecological diversity of weeds? How does the weed flora vary across surveys at different spatial and temporal scales? We built a comprehensive dataset of weed taxa in France by compiling weed flora, lists of specialized segetal weeds, and plot-based surveys in agricultural fields, with different spatial and temporal coverages. We informed life forms, biogeographical origins and conservation status of these weeds. We also defined a broader dataset of plants occupying open habitats in France and assessed habitat specialization of weeds and of other plant species absent from arable fields. Our results show that many arable weeds are frequently recorded in both arable fields and non-cultivated open habitats and are, on average, more generalist than species absent from arable fields. Surveys encompassing field edges included species also occurring in mesic grasslands and nitrophilous fringes, suggesting spill-over from surrounding habitats. A total of 71.5% of the French weed flora was not captured in plot-based surveys at regional and national scales, and many rare and declining taxa were of Mediterranean origin. This result underlines the importance of implementing conservation measures for specialist plant species that are particularly reliant on arable fields as a habitat, while also pointing out biotic homogenization of agricultural landscapes as a factor in the declining plant diversity of farmed landscapes. Our dataset provides a reference species pool for France, with associated ecological and biogeographical information.

Climate-human interaction associated with southeast Australian megafauna extinction patterns.

The mechanisms leading to megafauna (>44 kg) extinctions in Late Pleistocene (126,000-12,000 years ago) Australia are highly contested because standard chronological analyses rely on scarce data of varying quality and ignore spatial complexity. Relevant archaeological and palaeontological records are most often also biased by differential preservation resulting in under-representated older events. Chronological analyses have attributed megafaunal extinctions to climate change, humans, or a combination of the two, but rarely consider spatial variation in extinction patterns, initial human appearance trajectories, and palaeoclimate change together. Here we develop a statistical approach to infer spatio-temporal trajectories of megafauna extirpations (local extinctions) and initial human appearance in south-eastern Australia. We identify a combined climate-human effect on regional extirpation patterns suggesting that small, mobile Aboriginal populations potentially needed access to drinkable water to survive arid ecosystems, but were simultaneously constrained by climate-dependent net landscape primary productivity. Thus, the co-drivers of megafauna extirpations were themselves constrained by the spatial distribution of climate-dependent water sources.

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.

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.

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.