Brigitte Mangin

Mangin, B. (Auteur de Pierre Casadebaig

Eléna Cadic

Nicolas Blanchet

Marie-Claude Boniface

Sebastien Carrere

Jerome Gouzy

Ludovic Legrand

Baptiste Mayjonade

Nicolas Pouilly

Sébastien Carrère

Jérôme Gouzy

Thierry André

Marie Coque

Joël Piquemal

Marion Laporte

Patrick Vincourt

Stéphane Muños

Nicolas B Langlade
email: nicoas.lnaglade@inra.fr

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach of crop modeling and genome-wide association

de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Brief Summary

First, we developed a novel method to describe the abiotic stresses perceived by the plant using the dynamic ecophysiological model included in the SUNFLO crop model. This allowed us characterize 17 environments corresponding to a large range of cold, drought and nitrogen stresses situations.

Secondly, we characterized, in the 17 environments, the oil yield plasticity to those three abiotic stresses in a sunflower core-collection, defined as the linear response of oil yield to the modeled stresses. This allowed us to perform a genome-wide association study on these oil yield plasticity phenotypes. In this article, we describe the GWAS only for cold stress, for which we identified 9 regions bearing candidate genes. We restricted the publication of our results to cold because drought and nitrogen stress plasticities were genetically correlated and we believe the GWAS results on those phenotypes must be interpreted with caution. Nevertheless, the analysis on cold stress identified very promising candidate genes and represents a proof of concept of the general approach.

All together, our study represents a novel combination of crop modeling and quantitative genetics that addresses an important question for biology.

Introduction

Adaptation to climate change requires new crop varieties adapted to new management options. Adaptation of agriculture is a key factor to lessen the impact of climate change (Lobell et al., 2008). The crop exposition to unfavorable growing periods can be partially controlled by adapting crop management, i.e by shifting planting dates or choosing a cultivar with an adequate phenology [START_REF] Acosta-Gallegos | Phenological plasticity as an adaptation by common bean to rainfed environments[END_REF]. But such adaptations also have side-effects and while flowering can be successfully desynchronized from the period of occurrence of water deficit, crop emergence would be more exposed to cold stress. In this case, adaptation to climate change should also include the development of new crop varieties (Rosenzweig et al., 1994), with new or improved properties such as tolerance to cold or other abiotic stresses. Moreover, because the multiplicity of cultivation conditions (soils, climatic uncertainty), a single genotype can be exposed to random unfavorable growing conditions within its cultivation area, ultimately impacting the expected crop performance. To ensure a stable performance under uncertain conditions, newly developed genotypes not only need to be tolerant, i.e. adapted to a single type of environment (specialisation) but also plastic, i.e. to be able to adapt in most of growing conditions encountered in the targeted cultivation area (Sambatti and Caylor, 2007).

Phenotypic plasticity is a key process for crop productivity under climate change. One way plants will respond to these changes is through environmentally induced shifts in phenotype (phenotypic plasticity) (Nicotra et al., 2010). While the process of phenotypic plasticity is mostly studied on natural systems, its implications on crop productivity under climate change are interesting for plant breeding [START_REF] Dewitt | Integrated solutions to environmental heterogeneity: Theory of multimoment reaction norms[END_REF]Sadras et al., 2009). Empirically relationships between plant traits and environmental variables (norms of reaction) are known to vary within species in nature as well as in crop species, are considered heritable traits themselves subject to natural or artificial selection (Via and Lande, 1985;Sambatti and Caylor, 2007). However, crop growth in a fluctuating environment generates complex and dynamic interactions between plant and environment, under the control of cultural practices.

It is necessary to unravel and measure abiotic stress levels before assessing plasticity in plant traits. In these conditions, assessing plasticity in plant traits is limited by our capacity to unravel those interactions and estimate abiotic stress levels at the plant scale. Actually, each growth condition creates a unique combination of those stress levels, with possible identical combinations in different growth conditions. For example, crops growing in continental climates might be exposed to both cold (during emergence) and heat (during flowering) stresses; with high temperatures driving a strong evaporative demand and water deficit, which also limits plant nitrogen uptake from the transpired water stream (Kiani et al., 2016). Accordingly, in a given cultivation area, stresses are not independent (Vile et al., 2012) and need to be characterized and modelled prior studying their impact on plant traits.

Crop simulation and modeling can help to characterize environment from the crop point of view. Because the environment is the largest component of the phenotypic variability of most plant traits and of course crop yield, its quantitative characterization is of major importance (Lake et al., 2016). Crop simulation models are based on mathematical equations representing the crop growth and development as a function of environment (climate, soil and management). Such tools can give access to plant-level state variables, such as time-series of several abiotic stresses, in large range of growing conditions otherwise difficult to characterize with sensors. This methodology Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

was recently implemented and allowed to identify major types of water deficit patterns for rainfed wheat in the Australian target population of environments [START_REF] Chenu | Large-scale characterization of drought pattern: A continent-wide modelling approach applied to the Australian wheatbelt-spatial and temporal trends[END_REF] or for coupled thermal and water stress patterns for chickpea in Australian National Variety Trials (Lake et al., 2016). In sunflower, a crop model was developed by [START_REF] Casadebaig | SUNFLO, a model to simulate genotype-specific performance of the sunflower crop in contrasting environments[END_REF] and takes into account water, cold, heat, and nitrogen stresses to estimate their impact on grain yield and oil content.

Genetic studies of genotype-by-environment interaction (GxE) Although plasticity has long been recognized as an interesting trait in ecological and crop science since the pioneering work of Bradshaw in the 1950s, the identification of genetic variation involved in plasticity is scarce and recent. An approach to deal with phenotypic variation due to environmental effects is to develop multi-environmental QTL analysis to identify QTLxE interactions as reviewed in Van Eeuwijk et al. (2010). Although this approach was widely used for single or multiple traits, this can provide a greater sensitivity for environment-dependent QTL, inform on the stability of the QTL effect, but not on the environmental factor(s) the QTL is responding to. To overcome this issue, several studies compared two conditions, varying a single environmental factor and were able to study the genetic control of the plasticity to a specific stress for a particular trait (McKay et al., 2008;[START_REF] El-Soda | Quantitative trait loci and candidate genes underlying genotype by environment interaction in the response of arabidopsis thaliana to drought[END_REF]. et al., 2012), SUB1A for submersion in rice [START_REF] Fukao | The submergence tolerance regulator SUB1A mediates crosstalk between submergence and drought tolerance in rice[END_REF] and RAS1 for saline stress tolerance [START_REF] Ren | RAS1, a quantitative trait locus for salt tolerance and ABA sensitivity in arabidopsis[END_REF]. No gene showing a GxE interaction for drought stress nor for any yield-related trait were cloned or fine-mapped yet. This certainly lies respectively in the difficulty to characterize drought stress in natural conditions and in the complex genetic architecture of yield-related traits, which implies small effect mutations. Indeed successful examples were obtained with easy-to-phenotype traits (flowering time) and easy-to-setup environmental factors controlled one at a time.

This reductionist vision of the environment prove to be efficient for gene cloning. However, more complex and realistic approaches are needed to understand plant and crop responses to environmental factors and therefore breed for trait plasticity ultimately providing stable crops. In the context of the phenomics era [START_REF] Großkinsky | Plant phenomics and the need for physiological phenotyping across scales to narrow the genotype-to-phenotype knowledge gap[END_REF], with a greater description of the environment, and with the genomic tools accessible on many species, these approaches should flourish. However, to our knowledge, no study tackled yet the identification of the genetic control of plant plasticity to combined environmental factors in nature or field trials (Mahalingam, 2015). Thanks to prior crop modeling, this approach is now amenable, and shall take advantage of a precise description of the plant stress factors, of the statistical power of multi-environment trials and of the realistic nature of the measured traits.

In this study, we developed a novel approach that combined crop modeling and quantitative genetics to identify the genetic basis of oil yield plasticity in sunflower. Our methodology consisted in estimating four abiotic stress indicators in a range of environments using a crop model and in selecting the indicators explaining the best the grain yield plasticity of commercial varieties. The plant-level stress estimations on these varieties were used to characterize each experimental sites for the three selected abiotic stresses. We could then estimate the plasticity of oil yield for each stress and for every hybrid of a sunflower diversity panel cultivated in every site. Therefore, as a demonstration on cold stress, we could successfully perform a genome-wide association study to identify genomic regions putatively involved in oil yield plasticity to cold.

Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

Material and methods

Plant material

Association mapping was carried out on a panel of 317 inbred lines from INRA and sunflower breeding companies. This panel was a subset of the core collection of 384 inbred lines of [START_REF] Cadic | Combined linkage and association mapping of flowering time in sunflower (helianthus annuus l.)[END_REF] chosen for its diversity from an initial set of 752 inbred lines [START_REF] Coque | Sunflower line core collections for association studies and phenomics[END_REF]. It was comprised of both elite lines, parents of commercial hybrids, and lines with introgressions from several wild Helianthus accessions, including H. annuus, H. argophyllus and H. petiolaris.

Oil yield was observed on testcross progeny obtained by crossing panel lines with testers according to their status (maintainers of cytoplasmic male sterility [B-lines] or fertility restorers [R-lines]), as described in [START_REF] Cadic | Combined linkage and association mapping of flowering time in sunflower (helianthus annuus l.)[END_REF]. R-lines were crossed with the two CMS PET1 counterparts of B-line proprietary testers (FS71501 or AT0521) while the B-lines were crossed with two R-line testers: 83HR4gms and SOLR001M. 83HR4gms is derived from the 83HR4 line and was converted to female by the introduction of a genetic male sterility, SOLR001M is a proprietary line carrying PEF1 cytoplasmic male sterility (Serieys, 1984;[START_REF] Crouzillat | Cytoplasmic male sterility in sunflower: Comparison of molecular biology and genetic studies[END_REF] which it maintains, although it is a restorer for classical PET1 cytoplasm (Table 1).

Description of the multi-environment trial (MET)

From 2008 to 2010, eight locations located in the center and Southwest of France were planted with the testcross progeny. In six location, trials were conducted with and without irrigation, providing a total of 17 location x treatment x year combinations, (designated as environments). The panel lines were evaluated on the same tester in each environment (Table 1). Each experiment was an Augmented-Design (Federer, 1961) formed of blocks, with 24 or 30 entries replicated in two sub-blocks. Each sub-block was randomized separately and contained two to four control hybrids.

The climatic variability on experimental locations was summarized by computing the mean air temperature, the sum of water inputs (rainfall and irrigation) and the climatic water deficit (difference between precipitations, P and potential evapotranspiration, PET) on the cropping period (see supplementary Figure S1 and Table S1). No environment had a climatic water deficit (ranging from -177 to -458 mm), meaning that the climatic evaporative demand was always above the water supply (even accounting for irrigation). Rainfall on the cropping period ranged from a low 36 mm to 368.5 mm and the average amount of irrigation was 74 mm. Trials were performed on various soils depth leading to a soil water capacity (SWC) from 112 mm to 240 mm. Mean nitrogen fertilization was 60 kg/ha (eq. mineral nitrogen).

Among the 17 environments, we discarded 3 environments (AI09_I, AI09_NI, CO09_NI). The first two were outliers for the observed yield and the SUNFLO model failed at simulating yield phenotypes close to the observed one for the controls. The last one did not exhibit genotypic effect for the panel phenotypes, the estimated genotypic variance was judged non-significantly different from zero by a Z-ratio test in the linear mixed model used to correct for micro-environment effects and to predict the genotypic value of the yield panel lines.

Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

Intra environment phenotypic data analysis

Within each environment, the oil yield was adjusted for micro-environment effects using ASReml-R [START_REF] Butler | ASReml-r reference manual[END_REF], as described in [START_REF] Cadic | Combined linkage and association mapping of flowering time in sunflower (helianthus annuus l.)[END_REF]. A linear mixed model (named naïve in Cadic et al., 2012) with a random effect for the genotypic value of the panel lines, including blocks and sub-blocks as fixed effects was compared to two spatial models. The spatial models included (a) random effects of row and column or (b) a first-order autoregressive process in the residuals to take into account autocorrelation between neighbour plots. The three models were compared using the Akaike criterion (AIC) and the best linear unbiased predictors (BLUP) of the genotypic values were extracted from the best model according to AIC, for the next step of analysis.

Estimation and choice of abiotic stresses in each environment using simulation

SUNFLO is a process-based model for the sunflower crop which was developed to simulate the grain yield and oil concentration as a function of time, environment (soil and climate), management practice and genetic diversity [START_REF] Debaeke | Simulation de la réponse variétale du tournesol à l'environnement à l'aide du modèle SUNFLO[END_REF][START_REF] Casadebaig | SUNFLO, a model to simulate genotype-specific performance of the sunflower crop in contrasting environments[END_REF][START_REF] Casadebaig | SUNFLO, a model to simulate genotype-specific performance of the sunflower crop in contrasting environments[END_REF]. The model simulates the main soil and plant processes: root growth, soil water and nitrogen content, plant transpiration and nitrogen uptake, leaf expansion and senescence and biomass accumulation, as a function of main environmental constraints (temperature, radiation, water and nitrogen deficit).

This model is based on a conceptual framework initially proposed by Monteith (1977) and now shared by a large family of crop models (e.g. [START_REF] Jones | The dSSAT cropping system model[END_REF]Keating et al., 2003;[START_REF] Brisson | An overview of the crop model STICS[END_REF]. In this framework, the daily crop dry biomass (DM t ) is calculated as an ordinary difference equation (Equation 1) function of incident photosynthetically active radiation (PAR, MJ m -2 ), light interception efficiency (1-exp -kLAI t ) and radiation use efficiency (RUE t , g MJ -1 , Monteith (1994)). The light interception efficiency is based on Beer-Lambert's law as a function of leaf area index (LAI t ) and light extinction coefficient (k).

DM t =DM t-1 +RUE t (1-exp -kLAI t )PAR t (1)
Broad scale processes of this framework, the dynamics of LAI, photosynthesis (RUE) and biomass allocation to grains were split into finer processes (e.g plant phenologic development, leaf expansion and senescence, response functions to environmental stresses) to reveal genotypic specificity and to allow the emergence of GxE interactions. Globally, the SUNFLO crop model has about 50 equations and 64 parameters (43 plant-related traits and 21 environment-related). When evaluated on the presented MET dataset, the SUNFLO model predicted accurately the performance of control hybrids across environments: the root of mean square error (RMSE) was 0.3 t ha -1 , relative RMSE was 9.6 %, bias was -0.14 t ha -1 .

Using the SUNFLO model, we computed two indicators (continuous and discrete) per type of considered abiotic stresses to characterize the different environments (Table 2). Each indicator was integrated over three key periods:

vegetative stage (veg), flowering period (flo) and grain filling period (fil). We also considered the sum over two periods and during the whole cropping period, for a total of seven time periods per indicator. All these indicators corresponded to the mean stress felt by the control hybrids during the above periods.

Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

For water stress, the fraction of transpirable soil water (FTSW) represents yield limitation through water deficit (integration of 1 minus FTSW); ETR is the conditional sum of days, if the ratio of the real evapotranspiration (ET)

to potential evapotranspiration (PET) was less than 0.6 (threshold for photosynthesis limitation). For cold stress, LTs is the conditional sum of days if mean air temperature was below 20°C and LTi represent low temperatures impact on photosynthesis (integration of 1 minus equation 2). Heat stress indicators were computed following the same logic, albeit representing high temperatures impact on photosynthesis. Equations ( 2) and ( 3) are used in the crop model to define the radiation use efficiency (RUE) response to temperature (Villalobos et al., 1996). For nitrogen deficit, NAB is the amount of absorbed nitrogen in the considered cropping period and NNI is the sum of 1 minus nitrogen nutrition index, which indicates crop nitrogen deficit (Lemaire and Meynard, 1997;[START_REF] Debaeke | A species-specific critical nitrogen dilution curve for sunflower (helianthus annuus l.)[END_REF]. Three sunflower varieties (Melody, Pacific, Pegasol) which were used as controls within the environments had their phenotypic characteristics previously included in the SUNFLO model (Supplementary Table S2). The soil characteristics, the crop management and climate data were collected to allow simulation for each environment.

LTRUE=      0 ifT m <T b T m T ol -T b - T b T ol -T b if T b <T m <T ol 1 else (2) HTRUE=      1 ifT m <T ou T m T ou -T c - T c T ou -T c if T ou <T m <T c 0 else (3)
One soil characteristic, the soil depth, is difficult to observe and has a strong impact on the simulated yield. Instead of using the approximated value given by the breeders and the farmers, we adjusted this parameter by minimizing the empirical mean square error between the observed and the simulated yields.

Model choice to select the stress indicators was made with the AIC using the native R function lm. AIC of each model was computed for all three control varieties and model choice was made on the mean over the three Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

Genetic study

Estimation of plasticities of oil yield to water, nitrogen and cold stresses in the diversity panel

In order to get plasticity phenotypes that reflect the responses of the panel lines to the different abiotic factors, we adjusted the BLUP phenotype of each panel line with the following linear model:

Y ij =a i +b i DS j +c i CS j +d i NS j +ε ij (4)
where Y ij is the BLUP of the phenotype for the ith genotype, a i the potential phenotype in an environment with average stresses, DS j the water stress indicator in the jth environment calculated as the mean over the 3 controls, b i the slope linked to the water stress, CS j the cold stress indicator in the jth environment calculated as the mean over the 3 controls, c i the slope linked to the cold stress, NS j the nitrogen stress indicator in the jth environment calculated as the mean over the 3 controls, d i the slope linked to the nitrogen stress and ε ij the residual variance.

The covariates DS, CS, and NS were centered before computing the regression model. The three estimated covariate coefficients b i , ĉ i , and d i of this model are the plasticity phenotypes of interest, i.e. the genotype slopes in response to water, cold and nitrogen stresses respectively.

To compare model (Error! Reference source not found.) to a more simple regression model, we also fitted the BLUP phenotype of each panel line using the cold stress indicator as a single regressor.

Before computing the above linear models, the oil yield missing data were imputed using the missMDA R package (Josse et al., 2012;[START_REF] Josse | missMDA: A package for handling missing values in multivariate data analysis[END_REF]. All recorded traits (24 traits for a total of 198 traits x environments) were used to impute missing oil yield values. However, panel lines with less than six observations per oil yield trait over the environments were discarded.

Combined stresses in a single multi-stress index of plasticity

In order to compare the panel lines for the stability against multiple stresses together, we defined a multi-stress plasticity index accounting for the variance-covariance of the stress slopes, equals to:

( b i ,ĉ i , d i )V -1 ( b i ,ĉ i , d i ) t (5)
where V denotes the (3 by 3) variance-covariance matrix.

To visualize the strategy of the panel lines versus the three stresses, we drew a star representation using the star Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

Genotyping and map building.

A set of 197,914 SNP were used to produce an AXIOM® genotyping 96-array (Affymetrix, Santa Clara, CA, USA). These SNP were selected from either genomic re-sequencing or transcriptomic experiments. An additional set of 6,800 non-polymorphic sequences were added as controls. Combined with internal technical controls, the AXIOM® genotyping 96-array was designed with a total of 445,876 probesets. proposed by Mangin et al. (2012). They were mapped to the same position of the mapped SNP in maximum LD.

For this, we used as LD statistics the maximum of the r 2 V and r 2 VS measurements that correct for relatedness and for both structure and relatedness, respectively.

Association tests

Association mapping was based on a set of 65,534 SNP with MAF > 0.05. Similarly to our previous work [START_REF] Cadic | Combined linkage and association mapping of flowering time in sunflower (helianthus annuus l.)[END_REF], two association models were performed using EMMA (Kang et al., 2008)

Sl i =  c X ic α c +M il Θ l +u i +e i ( 6
) Sl i
is the slope for ith line, X ic is the line status, α c is the effect of the line status c, M il is the genotype of the ith line at locus l, Θ l is the effect of locus l. α c and Θ l are considered to be fixed effects. u i is the random polygenic effect modelling genetic relatedness with Var(u)=σ Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

tests (M eff

) based on the eigen-values of the SNP correlation matrix as proposed by Li and Ji (2005). We computed M eff by using blocks of 250 markers and assuming that the blocks were independent.

In order to conduct a multi-loci analysis, the genotypic data was imputed using Beagle v4.0 [START_REF] Browning | Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering[END_REF] and MLMM (Segura et al., 2012) was used for this purpose. We stopped the forward approach of MLMM when the variance of the polygenic term was non significant, which is a little different compared to the initial forward procedure of MLMM that stops when the estimator of this polygenic variance is equal to 0. The significance of the polygenic variance component was judged with a log-likelihood ratio test and a risk of 1%.

This log-likelihood ratio test compared the model with and without the polygenic effect. Two times the difference between the log-likelihood of the models is known to follow asymptotically a mixture between a Dirac at 0 and a Chi-square distribution with 1 degree of freedom (Self and Liang, 1987). When the forward approach stopped, a SNP was judged associated if its Bonferonni corrected p-value, using M eff as the number of independent tests, was inferior to the chosen type I error. The associated SNP were put all together with the polygenic term to create the final multi-loci linear model. The associated SNP effect, their reevaluated p-values were computed with the base R function lm in this final model [START_REF] Core | R: A language and environment for statistical computing[END_REF].

Functional annotation of associated SNP

Context sequences of associated SNP were compared to the sunflower genome (line XRQ) sequenced using the Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

Results

Estimation of abiotic stresses in field environment.

The characterization of abiotic stresses at the plant level through modelling and simulation explained observed yield variability better than climate-based indicators (Figure 1). For example, the correlation between yield observed in field experiments and water deficit computed from (1) climate data only (P -PET), (2) climate, soil, and management data (P -PET + SWC + irrigation) and (3) simulated plant data (ETR, defined in Table 2) was gradually stronger (respectively -0.01, -0.65 and -0.86). The correlation between yield and water deficit index increased in strength with the "proximity" of the regressor to the plant, revealing the expected negative impact of water deficit on crop yield.

The correlation between different type of abiotic stresses (Figure 2), indicated that cold (LT) and heat stresses Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

Model selection to estimate the best combination of abiotic stresses.

The model selected as having the best AIC on average over the control genotypes was a three indicator model, including a cold stress indicator during the vegetative stage (LTi veg ), a water stress indicator during the vegetative and the flowering period (FTSW veg+flo ), and a nitrogen indicator during the whole growth period (NAB). Table 3 presents the p-value of the Fisher test for each control in this best model and the proportion of variance explained by this best model for each control.

Each stress indicator had a significant effect on the grain yield for each control genotype. The highest impact was due to cold stress for Melody (p-value = 2.39 10 -4 ) and the smallest was due to water stress for Pacific (p-value = 2.91 10 -2 ) as indicated in Table 3. All together, these indicators explained very well the grain yield variability of each control since the percentages of explained variance were 90%, 93% and 94% for Pegasol, Pacific and Melody respectively.

Although cold stress has a strong impact on yield (-0.35 q h -1 d -1 ), the relatively low number of day of cold stress in the MET (8.7 on average, variation in the MET shown on Figure 3) reduces its impact, whereas water stress impacts the most yield due to the high number of days of stress (33.9 on average) even if sunflower is relatively tolerant (-0.14 q h -1 d -1 ).

Plasticity in the diversity panel.

As we described the actual cold, nitrogen and water stresses felt by sunflower in 14 environments, we also measured oil yield in a diversity panel in those environments. This allowed us to compute plasticities of oil yield to abiotic stresses for each line in the panel as the slopes of regression of oil yield to individual stress indicators.

These three slopes represented how plastic was the response of a line faced to a given stress and are therefore referred as plasticity phenotypes later. Minimum, mean, maximum, and variance of these plasticities as well as their correlation are presented in Table 4 and their distribution histograms are in supplementary Figure S2. Cold stress plasticity appears genetically independent to water and nitrogen stress plasticities. On the contrary, nitrogen and water stress plasticity are highly correlated (Pearson correlation of 0.62) suggesting common genetic control of these traits in our panel. Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

A multi-stress plasticity index.

Following the estimation of the three single stress plasticities, we were interested in calculating a multi-stress plasticity index to describe the general abiotic tolerance of every line. This index is a weighted sum of the three abiotic stress plasticities taking into account the correlation between them, specifically the one between drought and nitrogen. This index allowed us to rank the panel lines and to describe the different strategies observed in the panel to tolerate combined stresses. This is illustrated on Figure 5 that shows the multi-stress plasticity index of every panel lines against its mean oil yield in the MET. Stable panel lines, with a small multi-stress plasticity index, as well as panel lines belonging to the border of the point cloud were highlighted by a star representation (a triangle in our case) illustrating their abiotic stress tolerance strategy. First, we observed that unstable panel lines were generally more sensitive to cold stress compared to stable lines. This was confirmed by a significant correlation between the multi-stress plasticity index and the percentage of plasticity due to the cold stress (p-value of 6 10 -4 ) although this was not observed for nitrogen or water stress (p-value of 0.28 and 0.10 respectively).

Second, we could not identify a consensus strategy for stable panel lines. As examples: (i) the three similarly stable panel lines (multi-stress index around 45) having high oil yield (around 15) developed tolerance to all stresses but showed different plasticity patterns, (ii) the two most stable panel lines had opposite profile of stability with the most stable being sensitive to cold and the second to water.

Genetic map.

In order to position the genomic regions controlling the abiotic stress plasticity, we constructed a genetic map with the markers genotyped on the diversity panel using three RIL populations. The final genetic map was composed of 89,979 markers positioned in 4,782 genetic positions for a total distance of 1,398.5 cM. On this map, 4,094 markers were used to build the consensus map between the INEDI (XRQxPSC8) and FUxPAZ2 RIL populations.

Among the other markers, 27,663 were mapped using the INEDI population alone, 13,807 were mapped using the FUxPAZ2 population alone, 29,586 were located thanks to a genotyping by sequencing map on the RIL population RHA801xRHA280 (Kane et al., 2011) and the remaining 14,829 markers were placed by linkage disequilibrium. Details of these maps can be found in the supplementary Table S3.

Association study of oil yield plasticity to cold stress.

For the genetic analysis, we focused, as a proof of principle of the approach, on the plasticity to cold stress as it is the most impacting stress and appeared genetically independent from the others. Among the 65,534 association tests, the effective number of tests (Li and Ji, 2005) was estimated around 14,000. Using this effective number of tests, we kept SNP associated with a Bonferroni corrected p-value of 7.1 10 -6 for a family wise type I error of 10%.

Only 2 SNP were detected by a marker by marker association analysis using EMMA. The most significant SNP Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

or restorer status as a fixed structure effect. The second is located at the center of LG 17 and was detected without the line status effect.

We completed this study by a forward approach of the multi-loci association analysis (MLMM) for both models (with or without the maintainer/restorer status). Both models stopped with 6 SNP among which four were judged associated. None of them was in common between the two models but one SNP corresponded to the previously identified SNP in LG05.64. The most associated SNP was located at the end of the LG13.

In total 9 QTL could be identified using MLMM and EMMA procedures with or without population structures: two located on LG 5 and 10 and one on LG 9, 13, 14, 16 and 17. Their phenotypic effects on cold stress plasticity varies from 10 to 21% of the average plasticity in the panel (Table 5).

Genes located in QTL controlling oil yield plasticity to cold stress.

We were interested in genes containing associated SNP to link the genetic identification to molecular and physiological processes putatively involved in cold tolerance. All associated SNP were located within coding sequences as expected from the AXIOM® genotyping array design. Functional annotation of the corresponding genes pointed out homologues of: NPF3.1, LTP, CYS6, NPF5.3, GMII, RPD1, PPX1, HAOX2 and IAR4 (from the most to the least significantly associated, as shown in Table 5). Strikingly, two close homologues of oligopeptide transporters (NPF3.1 and NPF5.3) are present in associated QTL on LG 5 and 14 and two homologues of genes involved in root development (RPD1 and IAR4) on LG 9 and 17. In addition, homologues of a lipid transfer protein, a cystatin, an alpha-mannosidase, a protein phosphatase and an aldolase are also in QTL Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

Discussion

In this work, we developed a novel method to characterize the abiotic stress levels on different environments by using crop modelling and simulation and subsequently exploited it to identify genetic control of stress plasticity.

We implemented this environment characterisation method on 17 locations from a multi-environment trial for sunflower and four abiotic stresses (water, nitrogen, cold and heat). The SUNFLO model [START_REF] Casadebaig | SUNFLO, a model to simulate genotype-specific performance of the sunflower crop in contrasting environments[END_REF] was used to simulate stress patterns dynamics for three varieties used as controls in each location and integrated indicators were computed from these data considering different crop phenological stages and physiological processes (eight stress indicators over seven periods). We used a model selection approach to select the best linear model among combinations of stress indicators used as regressors for yield. Water, nitrogen and cold stresses were retained as the most explicative abiotic stresses for yield variability in this MET. Using reaction norms as conceptual reference, we computed abiotic stress plasticities as the slopes of the linear regression of oil yield on selected stress indicators. We then conducted an association study with a panel of 317 lines genotyped for nearly 65,000 markers on oil yield plasticity for cold stress as its was the most impacting stress (per time unit) and was not correlated to other abiotic stress plasticities.

Crop modelling helped to analyze abiotic stress patterns and to explain their impact on yield.

In a location, stress indicators can be climate-based (precipitations minus evapotranspiration), crop-based (simple water balance including soil water capacity and irrigation) or plant-based (simulated dynamic water balance). We observed that the observed grain yield was best explained by plant-based stress indicators because the interactions between climate, leaf area dynamics, plant stomatal conductance (isohydric vs anisohydric behaviours, e.g. [START_REF] Casadebaig | Thresholds for leaf expansion and transpiration response to soil water deficit in a range of sunflower genotypes[END_REF] for sunflower) and management practices can be partly reproduced by the crop model algorithm. This is for example illustrated by environment rankings, where some irrigated locations (CO08_I, CO09_I) still display a high level of water stress while a rainfed one (CA10_NI, VE09_NI) show reduced water deficit.

Abiotic stresses also do not have the same impact on crop physiology according to their timing of occurrence during the crop cycle (Table 3). Among the seven possible combinations between main crop phases (i.e. vegetative, flowering, grain filling), we indeed observed that the relevance of these timings was specific to the type of abiotic stress. For cold stress, the detection of early crop growth (vegetative period) was expected because this stress occurrence is strongly determined by the climate (low temperature during crop installation).

However, in continental climates, where sunflower is mainly grown, we can also observe cold temperatures at the end of the crop cycle. For water stress, where interactions between crop growth and climate variability are more important, vegetative and flowering periods were identified, which is consistent with numerous previous reports on sunflower [START_REF] Blanchet | Articulation des divers processus d'adaptation à la sécheresse et comportements globaux du tournesol[END_REF][START_REF] Cabelguenne | EPICphase, a version of the EPIC model simulating the effects of water and nitrogen stress on biomass and yield, taking account of developmental stages: Validation on maize, sunflower, sorghum, soybean and winter wheat[END_REF]. Regarding nitrogen stress, the importance of this process over the whole crop cycle was highlighted. Indeed, recent reports indicated that post-flowering nitrogen absorption could also be significant [START_REF] Andrianasolo | Source and sink indicators for determining nitrogen, plant density and genotype effects on oil and protein contents in sunflower achenes[END_REF]. Remarkably, heat stress was not identified as a major contribution to yield variability in the MET. Actually, depreciative effect of high temperatures on photosynthesis were caused by temperatures that were almost never reached in our experimental conditions. Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

source not found.) was only significant in one location (CO09) and null (9/17 locations) or weak in the others (Figures 2 andS1).

Genetic control of plasticity of oil yield to cold stress

The most and fourth most associated SNP pointed to two genes highly homologous to NPF3.1 and NPF5.3 that are both oligopeptide transporters. These two independent association signals on chromosomes 5 and 14 strongly suggest a role of oligopeptide transport in tolerance to cold stress observed in our experimental conditions i.e. when young plants are exposed to chilling. In plants, these transporters are key players in nitrogen nutrition and therefore plantlet growth. The importance of oligopeptide transport to tolerate cold is corroborated by the demonstrated molecular adaptation of this transporter family in antarctic icefish (Chionodraco hamatus) adapted to sub-zero temperatures (Maffia et al., 2003;[START_REF] Rizzello | Protein cold adaptation strategy via a unique seven-amino acid domain in the icefish (chionodraco hamatus) PEPFS71501 transporter[END_REF]. The role in N nutrition of these transporters in animal and plants indicates that nutrient transport can be a limiting factor at low temperature that likely limits remobilization of seed stocks and/or absorption and transport of N from roots to aerial organs.

The second most associated SNP is located in a putative lipid transfer protein (LTP) (QTL LG16.48). Many LTP have been reported to be transcriptionally induced by freezing (reviewed in Liu et al., 2015) and over-expression of LTP3 provided freezing tolerance in A. thaliana [START_REF] Guo | Lipid transfer protein 3 as a target of MYB96 mediates freezing and drought stress in arabidopsis[END_REF]. This action could be due to membrane stabilization as demonstrated by Hincha et al. (2001) in preventing chloroplastic damages induced by freezing, or through its role in seed lipid mobilization during germination and seedling growth [START_REF] Pagnussat | On the role of a lipid-transfer protein. arabidopsis ltp3 mutant is compromised in germination and seedling growth[END_REF] as shown in sunflower for another LTP [START_REF] Pagnussat | Unexpected localization of a lipid transfer protein in germinating sunflower seeds[END_REF]. Another candidate genes is the cystatin CYS6 homologue located in QTL LG13.72. Several homologues of cystatin were shown to be induced during cold exposure in barley [START_REF] Gaddour | A constitutive cystatin-encoding gene from barley (icy) responds differentially to abiotic stimuli[END_REF], maize (Massonneau et al., 2005), wheat (Talanova et al., 2012), and increase, when over-expressed, cold tolerance in Arabidopsis (Zhang et al., 2008). Based on [START_REF] Prins | Cysteine proteinases regulate chloroplast protein content and composition in tobacco leaves: A model for dynamic interactions with ribulose-1,5-bisphosphate carboxylase/oxygenase (rubisco) vesicular bodies[END_REF], the sunflower cystatin could provide a better regulation of Rubisco turnover in chloroplasts in cold conditions.

The alternative oxidase HAOX2 homologue found in QTL LG10.34 constitutes another good candidate gene for cold tolerance. The Alternative Oxidase Pathway (AOP) has been described in many plants to be involved in cold stress response as a biochemical protection against overproduction of Reactive Oxygen Species due to the cold inhibition of the electron transport chain in mitochondria [START_REF] Feng | Chilling tolerance of wheat seedlings is related to an enhanced alternative respiratory pathway[END_REF]. Furthermore, the AOP was shown to participate in differential cold-sensitivity between two maize genotypes [START_REF] Ribas-Carbo | The electron partitioning between the cytochrome and alternative respiratory pathways during chilling recovery in two cultivars of maize differing in chilling sensitivity[END_REF].

Interestingly, two genes (homologues to RPD1, IAR4) were found (QTL LG09.27, LG17.49) and their Arabidopsis counterparts share similar features: both are involved in root development and both mutants show temperature-sensitive phenotypes (at 20°C and 28°C) (Konishi and Sugiyama, 2006;[START_REF] Quint | Arabidopsis IAR4 modulates auxin response by regulating auxin homeostasis[END_REF]. This suggests that root setting could also be temperature-and genotype-dependent in sunflower. All together, the functional annotation of QTL associated to cold stress plasticity of oil yield identified several candidate genes and physiological processes. Most of them were already described to be involved in cold stress tolerance in other plants which supports our study and indicates a probably short genetic distance between associated SNP and causal mutations.

To complete our understanding of how these processes that act during the early growth phase of sunflower, impact the final seed yield, further studies putting in relation dynamic measurements of plant growth rate and cell physiology would be enlightening.

Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

Phenotypic plasticity and tradeoff for potential yield

In the studied multi-environment network, mean oil yield was positively correlated with a high sensitivity to environmental stresses, indicating that a global gain in performance was generally associated to a higher yield instability (Figure 5). This is globally in accordance with previous claims that an increase in phenotypic plasticity allowed to achieve better yield stability across environments but at the expense of greater performance in low stress conditions (Sadras et al., 2009). However the presence of both high-yielding and stress-tolerant genotypes suggests this general observation can be genetically by-passed and leaves room for more efficient sunflower varieties.

Stability of complex traits such as yield or fitness depends on plasticity of numerous intermediate traits (likely physiological and developmental processes) that are yet unknown. In our approach, we studied directly the plasticity of the complex trait with the idea of stabilizing it. The molecular and physiological processes pointed out by the genetic analysis allow us to identify some candidate processes: oligopeptide transport, root development, ROS scavenging, chloroplast and mitochondrial physiology (i.e. intermediate traits). In this kind of approach, we can wonder whether we could detect key regulators such as transcription factors (none in our case).

Indeed, genetic variation in those would likely have trade-off effects on various physiological processes. Then, they would impact intermediate traits significantly but with possibly opposite effects and at the end, no significant impact on the resulting complex trait.

On the breeding strategy point of view, the considered environmental stresses (temperature, nitrogen and water)

do not necessarily coexist in the French target population of environments: e.g. south-western production regions are exposed to drought and heat stress while the temperatures in northern regions are low enough and necessitate short crop cycle cultivars. This lack of spatial superposition of environmental constraints allow to exploit the differential sensitivities in the studied genetic material and adapt cultivar choice or breeding according to local growing environment. On the other hand, breeding for phenotypic plasticity of intermediate traits would potentially result in resilience to increasingly unpredictable environments (Nicotra et al., 2010).

Conclusion

Improving crop performance in low-input cropping systems requires a coordinated improvement of genotypes and agronomical practices (Sadras and Denison, 2016). In these growth-limiting conditions, abiotic stresses occur in combined and dynamic patterns. Therefore, disentangling those using models allows to understand their specific impacts on complex traits (such as yield) and the genetic factors potentially reducing those.

In our study, precise characterization of water, cold, heat and nitrogen stresses allowed accurate identification of nine QTL and underlying genes controlling stress plasticity. This joint approach between crop modelling and quantitative genetics also permitted estimations of allelic variation in natural conditions.

Such inter-disciplinary approach should be useful to conduct different breeding strategies and adapt crop to Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved. Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved. Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved. Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

Table 4 Variation and correlation of oil yield plasticity for different abiotic stresses in the diversity panel. Minimum, mean, maximum, and variance of the plasticity phenotypes and their correlations. The plasticity phenotypes are calculated as the slopes in a linear model including the three stress indicators that best characterized the environments. a cold stress plasticity (q h -1 d -1 ), b nitrogen stress plasticity (q h -1 kg -1 ), c water stress plasticity (q h -1 d -1 ). Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved. For water stress, the fraction of transpirable soil water (FTSW) represent yield limitation through water deficit (integration of 1 minus FTSW); ETR is the conditional sum of days, if the ratio of the real evapotranspiration (ET) to potential evapotranspiration (PET) was less than 0.6 (threshold for photosynthesis limitation). For cold stress, LTs is the conditional sum of days if mean air temperature was below 20°C and LTi represent low temperatures impact on photosynthesis (integration of 1 minus equation 2). Heat stress indicators were computed following the same logic, albeit representing high temperatures impact on photosynthesis. Equations ( 2) and (3) are used in the crop model to define the radiation use efficiency (RUE) response to temperature (Villalobos et al., 1996). For nitrogen deficit, NAB is the amount of absorbed nitrogen in the considered cropping period and NNI is the sum of 1 minus nitrogen nutrition index, which indicates crop nitrogen deficit (Lemaire and Meynard, 1997;[START_REF] Debaeke | A species-specific critical nitrogen dilution curve for sunflower (helianthus annuus l.)[END_REF].

Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved. Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

  controls. Compared linear models, which fitted the grain yield with the indicators, were limited to combinations of only one indicator per type of stresses. All models integrating one, two, three or four stresses were compared; in total 50,624 models were computed and compared. The R function lm was also used to compute the p-value of the Fisher test for each indicator including in the best model. Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  function of the R package graphics. The three plasticity phenotypes were first taken in absolute value and scaled to 0-1 in order to have three comparable values with values close to 0 noted the stability and values close to 1 noted the instability. Then they were normalized by panel line to sum to 1 in order to get comparable values for all panel lines. The star representation was applied on these scaled and normalized values that can be interpreted as the percentage of stability dedicated to each stress. Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  where K ais is an AIS matrix and Var(e)=σ 2 e I where I denotes the identity matrix. Multiple testing correction was achieved using an approximate effective number of Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  //www.heliagene.org/HanXRQ-SUNRISE/) by blast analysis. In case of ambiguous positioning on the genome, we retained the chromosomal position in accordance to genetic map location of the marker. Similarly, we positioned associated SNP on the reference transcriptome (https://www.heliagene.org/HaT13l/). Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

(

  HT) were naturally the most highly negatively correlated combination. The nitrogen indicator NAB was not correlated to other indicators except the negative correlation with the water stress indicator (ETR, particularly during grain filling). The other nitrogen indicator NNI, which was not correlated to NAB showed both a positive correlation with heat stress indicators and a negative correlation with cold indicators. The water stress indicator FTSW during grain filling was also positively correlated to heat stress indicators (HT) and nitrogen stress (NNI), and negatively correlated to cold stress indicators. Within a single type of stress, correlations between indicators computed during different cropping periods were generally positive and high (Figure 2), ranging 0.40 for water deficit (FTSW) to 0.98 for heat stress (HTi), with nitrogen and cold indicators in-between. Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  To illustrate the importance of multi-stress indicator modelling, we compared regression of oil yield in a single stress indicator model (including only cold stress indicator) and in the multi-stress indicator model in Figure 4. Oil yields of the most sensitive and tolerant lines were plotted against the cold stress indicator. Point clouds are closer to the regression lines in the multi-stress indicator model indicating a better characterization of cold stress impact in this modelling approach. Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  was located at the end of LG 5 in a QTL named LG05.64, and was detected with a model including the maintainer Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  controlling oil yield plasticity to cold stress. Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  According to the current parameterization of the simulation model, heat stress indicator (Error! Reference Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  climate variability through local adaptation: maximizing performance in a given environment type, or through global adaptation: maximizing yield stability over different environment types. Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017). Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach This article is protected by copyright. All rights reserved. Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  dt Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  . (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  . (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

  . (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approachThis article is protected by copyright. All rights reserved.

Figure 1

 1 Figure 1 Relation between observed grain yield and several abiotic stress indicators. Panels A, B, and C display the regression line between grain yield and water stress indicators, computed at different levels: using climatic data only (panel A, precipitation -potential evapotranspiration, Pearson correlation (r) of -0.01), using both climatic and crop data (panel B, precipitation -potential evapotranspiration + irrigation + soil water capacity, r = -0.65, p-value = 9.6 10 -6 ) and using simulated plant data (panel C, evapotranspiration ratio, r = -0.86, p-value = 3.3 10 -12). Panel D displays grain yield predicted as a function of combined abiotic stresses (linear model of water, nitrogen, and cold stress indicators, r = 0.91, root mean square error of 0.21 t ha -1 ).

Figure 2

 2 Figure2Heat map of the correlations between stress indicators. The three selected stress indicators are indicated in blue, yellow and green for water, cold and nitrogen respectively. For water stress, the fraction of transpirable soil water (FTSW) represent yield limitation through water deficit (integration of 1 minus FTSW); ETR is the conditional sum of days, if the ratio of the real evapotranspiration (ET) to potential evapotranspiration (PET) was less than 0.6 (threshold for photosynthesis limitation). For cold stress, LTs is the conditional sum of days if mean air temperature was below 20°C and LTi represent low temperatures impact on photosynthesis (integration of 1 minus equation 2). Heat stress indicators were computed following the same logic, albeit representing high temperatures impact on photosynthesis. Equations (2) and (3) are used in the crop model to define the radiation use efficiency (RUE) response to temperature(Villalobos et al., 1996). For nitrogen deficit, NAB is the amount of absorbed nitrogen in the considered cropping period and NNI is the sum of 1 minus nitrogen nutrition index, which indicates crop nitrogen deficit(Lemaire and Meynard, 1997;[START_REF] Debaeke | A species-specific critical nitrogen dilution curve for sunflower (helianthus annuus l.)[END_REF].

Figure 3

 3 Figure 3 Range of variation in abiotic stresses in the multi-environment trial. Water (A), cold (B), and nitrogen stress (C) indicators, computed by the SUNFLO crop model for each environment, and averaged over the three control genotypes. The SUNFLO crop model was used to compute stress indicators for three control commercial hybrids (Melody, Pacific, Pegasol) according to the developmental stage of the crop: vegetative growth (green), flowering (yellow), and seed filling (brown).

Figure 4

 4 Figure 4 Regression of oil yield against the cold stress indicator for the most tolerant (green triangles) and sensitive (red squares) panel lines showing the ability of multi-stress modelling to better characterize the environment. A) Regression in a single (cold) stress indicator model. B) Regression in a multi-stress (cold, nitrogen, water) indicator model.

Figure 5

 5 Figure 5 Relation between the multi-stress plasticity index and the average oil yield showing differences in stress tolerance strategies in the diversity panel of sunflower lines. The three-branch star represents the strategy of the most extreme panel lines in response to combined abiotic stresses. The length of star branches represents the relative plasticity against the corresponding stress: longer equals more sensitive.

  Genomic DNA from the 317 panel lines and two recombinant inbred lines (RIL) populations INEDI and FUxPAZ2 obtained from the cross between XRQ and PSC8 lines (180 RIL) and from the cross between FU and PAZ2 lines (87 RIL) respectively, were genotyped with the AXIOM® array. All hybridization experiments were performed by Affymetrix and the genotypic data were obtained with the GTC software (Affymetrix).From the 197,914 SNP, 35,562 were polymorphic between XRQ and PSC8 and 28,529 between FU and PAZ2.We used CarthaGène v1.3[START_REF] Givry S De | CAR(H)(T)AGene: Multipopulation integrated genetic and radiation hybrid mapping[END_REF] to build the genetic maps of the INEDI population and the FUxPAZ2 population separately. For the INEDI population, we added to the set of AXIOM® SNP the markers

previously mapped by

[START_REF] Cadic | Combined linkage and association mapping of flowering time in sunflower (helianthus annuus l.)[END_REF]

which allowed completing this map and assigning AXIOM® markers to appropriate linkage group (LG). We built a consensus map with common markers of both previous maps using Biomercator v4.0 (Sosnowski et al., 2012) and we projected the specific markers of the two previous maps on this consensus map. Unmapped AXIOM® SNP were placed by BLAST analysis on the RHA280xRHA801 genetic map based on 454 sequencing (Kane et al., 2011) and finally projected on the INEDI and FUxPAZ2 consensus map. The remaining AXIOM® SNP were located by computing the linkage disequilibrium (LD) measurements

Table 1 Details on location, treatment, year, testers and observed hybrids for the 17 environments

 1 . a The locations are designated as follows: AI: Aigrefeuille (Center West), CA: Castelnaudary (South West), CO : Cornebarrieu (South West), GA: Gaillac (South West), VE: Verdun (South West), LO: Loudun (Center West), SE: Segoufielle (South West), CHA: Chateauroux (Center). b I for irrigation, NI for non-irrigation. Environment Location a Treatment b Year Tester for B-line Tester for R-line Number of hybrids

	AI08_I	AI	I	2008	83HR4gms	FS71501	192
	AI08_NI	AI	NI	2008	83HR4gms	FS71501	193
	CO09_I	CO	I	2009	83HR4gms	FS71501	278
	CO09_NI	CO	NI	2009	83HR4gms	FS71501	278
	GA09_I	GA	I	2009	83HR4gms	FS71501	275
	GA09_NI	GA	NI	2009	83HR4gms	FS71501	274
	LO10_NI	LO	NI	2010	83HR4gms	FS71501	284
	VE10_I	VE	I	2010	83HR4gms	FS71501	289
	AI09_I	AI	I	2009	SOLR001M	AT0521	280
	AI09_NI	AI	NI	2009	SOLR001M	AT0521	280
	VE09_I	VE	I	2009	SOLR001M	AT0521	273
	VE09_NI	VE	NI	2009	SOLR001M	AT0521	273
	CA10_NI	C1	NI	2010	SOLR001M	AT0521	306
	CO08_I	CO	I	2008	SOLR001M	AT0521	249
	CO08_NI	CO	NI	2008	SOLR001M	AT0521	249
	SE10_NI	SE	NI	2010	SOLR001M	AT0521	285
	CHA10_I	CHA	NI	2010	SOLR001M	AT0521	306

Table 2 Description of abiotic stress indicators simulated by the crop model.

 2 The SUNFLO crop model was used to simulate the interactions between plant growth and available environmental resources. The evolution of resource level and abiotic constraints during the crop cycle was summarized via computing eight stress indicators during seven cropping periods: vegetative, flowering, grain filling, and their combination. 1 [x] equals 1 if x is true and 0 else.

	Stress	Symbol	Description	Unit	Formula
	temperature	HTi	high temperature (continuous)	-	  1-HTRUE dt
	temperature	HTs	high temperature (discrete)	d	 1 [T m >28]
	temperature	LTi	low temperature (continuous)	-	  1-LTRUE dt
	temperature	LTs	low temperature (discrete)	d	 1 [T m <20]
	water	FTSW	Edaphic water deficit (continuous)	-	  1-FTSW dt
	water	ETR	Edaphic water deficit (discrete)	d	 1 [ET/PET<0.6]
	nitrogen	NAB	Absorbed nitrogen	kg ha -1	  NAB dt

Table 3 Results of the regression between grain yield and abiotic environmental indicators: cold, nitrogen, and water stresses.

 3 The plasticity of the three control genotypes (i.e., commercial varieties) are presented as their slopes and p-values of the Fisher test in the best-fit model, with the proportion of explained variance for each control genotype. a LTi_veg: integration of low temperatures during the vegetative stage, b NAB: absorbed nitrogen during the whole growth period, c FTSW_veg+flo: fraction of the transpirable soil water during the vegetative and the flowering periods. The best explanatory linear model was determined with AIC and was selected from all regression models with one to four linear regressors and one indicator per abiotic stress.

						Unit	Melody	Pacific	Pegasol
	Yield					q ha	-1	32.41	29.18	31.30
	Cold stress	a		slope	q ha	-1 d -1	-0.36	-0.36	-0.32
	Cold stress			p-value					2.39 10 -4	1.97 10	-3	1.95 10 -3
	Nitrogen stress	b	slope	q ha	-1 kg	-1	0.27	0.22	0.25
	Nitrogen stress	p-value					1.16 10 -3	7.66 10	-3	6.6 10 -4
	Water stress	c		slope	q ha	-1 d -1	-0.16	-0.12	-0.13
	Water stress		p-value					3.27 10

Table 5 Positions and estimated effects of SNP and genes associated with oil yield plasticity to cold stress.

 5 Genetic distance determined on a consensus map built from three RIL populations (see details in Materials and Methods). Physical position was determined via BLASTing against the genome of line XRQ (early access to version HanXRQv1.1). Association tests were performed using the MLMM procedure unless otherwise stated and either with or without the maintainer/restorer status structure (noted as B/R status). Effects of SNP were estimated using a linear model including all associated SNP. Their relative effect compared with the average plasticity of the diversity panel is indicated. Sunflower gene names correspond to the genome HannXRQ v1.1 (early access). The closest Arabidopsis thaliana homologue is indicated with its TAIR code, name, and a brief functional description.

	SNP name	QTL	LG	Genetic	Physical	Association	Association model	Effect	Effect
				distance	position (bp)	test p-value		on plasticity	percentage
				(cM)
	AX-84511295	LG05.17	5	16.6	14,284,297	1.54 10 -8	Without B/R status	2.80 10-2	12%
	AX-84248033	LG05.64	5	64.4	206,893,338	5.70 10 -6	EMMA with B/R status	3.15 10 -2	13%
						3.64 10 -6	MLMM with B/R status	3.34 10 -2	14%
	AX-84488969	LG09.27	9	27.4	203,648,151	6.76 10 -7	With B/R status	-2.35 10 -2	-10%
	AX-84358846	LG10.34	10	34.1	220,548,129	3.87 10 -6	Without B/R status	-5.06 10 -2	-21%
	AX-84316040	LG10.45	10	44.8	93,526,595	5.11 10 -7	With B/R status	-3.14 10 -2	-13%
					or 217,685,071
	AX-84586324	LG13.72	13	72.4	190,712,390	1.17 10 -7	Without B/R status	-2.94 10 -2	-12%
	AX-84436021	LG14.27	14	27.0	168,814,054	2.70 10 -7	With B/R status	3.61 10 -2	15%
	AX-84337313	LG16.48	16	47.7	129,668,002	6.65 10 -8	Without B/R status	-2.72 10 -2	-12%
	AX-84507515	LG17.49	17	49.1	153,030,442	6.81 10 -6	EMMA without B/R status	3.55 10 -2	15%
	(Table 5 continued...)
	SNP name	Sunflower gene		TAIR code		Description	Arabidopsis name
	AX-84511295	HannXRQ_Chr05g0131961	AT1G68570		Oligopeptide transporter		NPF3.1
	AX-84248033	HannXRQ_Chr05g0159711	AT4G26720		Protein phosphatase		PPX1
	AX-84488969	HannXRQ_Chr09g0274811	AT4G33495	Root development, Temperature sensitive	RPD1
	AX-84358846	HannXRQ_Chr10g0311881	AT3G14150		ROS protection		HAOX2
	AX-84316040	HannXRQ_Chr10g0292461	AT5G14950	N-glycan mannose hyper-osmotic salinity stress	GMII
		or HannXRQ_Chr10g0311021
	AX-84586324	HannXRQ_Chr13g0424331	AT3G12490	Cystatin, protease inhibitor, Tolerance to abiotic stresses	CYS6
							incl. cold
	AX-84436021	HannXRQ_Chr14g0460351	AT5G46040		Oligopeptide transporter		NPF5.3
	AX-84337313	HannXRQ_Chr16g0519611	AT3G07450		Lipid Transfer Protein

Acknowledgement

This work benefited from the GENOPLANTE program ''HP1' ' (2001-2004), the ''SUNYFUEL'' project, financially supported by the French National Research Agency (ANR-07-GPLA-0022, 2008-2011), the ''OLEOSOL '' project (2009-2012) with the financial support from the Midi Pyrénées Region, the European Fund for Regional Development (EFRD), and the French Fund for Competitiveness Clusters (FUI), and the ''SUNRISE'' project of the French National Research Agency (ANR-11-BTBR-0005, 2012-2019). Version postprint Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.

Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved. Josse J, Chavent M, Liquet B, Husson F. 2012. Handling missing values with regularized iterative multiple correspondence analysis. Journal of Classification 29, 91-116. Kane NC, Gill N, Lindhauer MG et al. 2011. Progress towards a reference genome for sunflower. Botany 89, 429-437. Kang HM, Zaitlen NA, Wade CM, Kirby A, Heckerman D, Daly MJ, Eskin E. 2008. Efficient control of population structure in model organism association mapping. Genetics 178, 1709-1723. Keating BA, Carberry PS, Hammer GL, Probert ME, Robertson MJ, Holzworth D, Huth NI, Hargreaves JNG, Meinke H, Hochman Z. 2003. An overview of APSIM, a model designed for farming systems simulation. European Journal of Agronomy 18, 267-288. Kesari R, Lasky JR, Villamor JG, Marais DLD, Chen Y-JC, Liu T-W, Lin W, Juenger TE, Verslues PE. 2012. Intron-mediated alternative splicing of arabidopsis p5CS1 and its association with natural variation in proline and climate adaptation. Proceedings of the National Academy of Sciences 109, 9197-9202. Kiani M, Gheysari M, Mostafazadeh-Fard B, Majidi MM, Karchani K, Hoogenboom G. 2016. Effect of the interaction of water and nitrogen on sunflower under drip irrigation in an arid region. Agricultural Water Management 171, 162-172. Konishi M, Sugiyama M. 2006. A novel plant-specific family gene, ROOT PRIMORDIUM DEFECTIVE 1, is required for the maintenance of active cell proliferation. Plant Physiology 140, 591-602. Lake L, Chenu K, Sadras V. 2016. Patterns of water stress and temperature for australian chickpea production. Crop and Pasture Science. Lecoeur J, Poiré-Lassus R, Christophe A, Pallas B, Casadebaig P, Debaeke P, Vear F, Guilioni L. 2011. Quantifying physiological determinants of genetic variation for yield potential in sunflower. SUNFLO: a model-based analysis. Functional Plant Biology 38, 246-259. Lemaire G, Meynard J. 1997. Use of the nitrogen nutrition index for the analysis of agronomical data. In: Diagnosis of the Nitrogen Status in Crops. Ed. G. Lemaire. Springer-Verlag, Berlin, 45-55. Li J, Ji L. 2005. Adjusting multiple testing in multilocus analyses using the eigenvalues of a correlation matrix. Heredity 95, 221-227. Liu F, Zhang X, Lu C, Zeng X, Li Y, Fu D, Wu G. 2015. Non-specific lipid transfer proteins in plants: Presenting new advances and an integrated functional analysis. Journal of Experimental Botany 66, 5663-5681. Lobell DB, Burke MB, Tebaldi C, Mastrandrea MD, Falcon WP, Naylor RL. 2008. Prioritizing climate change adaptation needs for food security in 2030. Science 319, 607-610. Maffia M, Rizzello A, Acierno R, Verri T, Rollo M, Danieli A, Döring F, Daniel H, Storelli C. 2003. Characterisation of intestinal peptide transporter of the antarctic haemoglobinless teleost chionodraco hamatus. Journal of Experimental Biology 206, 705-714. Mahalingam R. 2015. Consideration of combined stress: A crucial paradigm for improving multiple stress tolerance in plants. In: Mahalingam R, ed. Combined stresses in plants. Springer International Publishing, 1-25. Version postprint Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved. Mangin B, Siberchicot A, Nicolas S, Doligez A, This P, Cierco-Ayrolles C. 2012. Novel measures of linkage disequilibrium that correct the bias due to population structure and relatedness. Heredity 108, 285-291. Massonneau A, Condamine P, Wisniewski JP, Zivy M, Rogowsky PM. 2005. Maize cystatins respond to developmental cues, cold stress and drought. Biochimica Et Biophysica Acta-Gene Structure and Expression 1729, 186-199. McKay J, Richards J, Nemali K, Sen S, Mitchell-Olds T, Boles S, Stahl E, Wayne T, Juenger T. 2008. Genetics of drought adaptation in arabidopsis thaliana II. QTL analysis of a new mapping population, kAS-1 x tSU-1. Evolution 62, 3014. Monteith JL. 1977. Climate and the Efficiency of Crop Production in Britain. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 281, 277-294. Monteith JL. 1994. Validity of the correlation between intercepted radiation and biomass. Agricultural and Forest Meteorology 68, 213-220. Nicotra AB, Atkin OK, Bonser SP et al. 2010. Plant phenotypic plasticity in a changing climate. Trends in Plant Science 15, 684-692.

Comment citer ce document : Mangin, B. (Auteur de correspondance), Casadebaig, P., Cadic, E., Blanchet, N., Boniface, M.-C., Carrere, S., Gouzy, J., Legrand, L., Mayjonade, B., Pouilly, N., André, T., Coque, M., Piquemal, J., Laporte, M., Vincourt, P., Munos, S., Langlade, N. (Auteur de correspondance) (2017).

Genetic control of plasticity of oil yield for combined abiotic stresses using a joint approach

This article is protected by copyright. All rights reserved.