Marion Prudent

Mathilde M Causse

Michel Génard

Pasquale Tripodi

Silvana Grandillo

Nadia Bertin
email: nadia.bertin@avignon.inra.fr

Causse Mathilde

Génard Michel

Tripodi Pasquale

Grandillo Silvana

Genetic and physiological analysis of tomato fruit weight and composition: influence of carbon availability on QTL detection

Keywords: Carbon allocation, physiological processes, fruit load, genetic variability, genotype x environment interaction, quality, QTL, Solanum lycopersicum, Solanum chmielewskii, sugar content

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

Introduction

Plant domestication is the genetic modification of wild species to meet human needs [START_REF] Doebley | The molecular genetics of crop domestication[END_REF]. In tomato, wild species originated in the Andes region where their tiny sweet fruits propagate the species. Tomato domestication dramatically increased fruit yield, and changes mainly occurred in fruit morphology or plant development. The most obvious evolution has been the massive increase in fruit size [START_REF] Tanksley | The genetic, developmental and molecular bases of fruit size and shape variation in tomato[END_REF][START_REF] Bai | Domestication and breeding of tomatoes: What have we gained and what can we gain in the future?[END_REF][START_REF] Cong | Regulatory change in YABBY-like transcription factor led to evolution of extreme fruit size during tomato domestication[END_REF] which was associated with a reduction in sugar content and subsequently in sweet flavour. Solublesolids content, mainly constituted by sugars, is high in wild tomato species such as Solanum pimpinellifolium or Solanum chmielewskii with more than 6% and 10% of the fruit fresh weight, respectively [START_REF] Rick | High soluble-solids content in large fruited tomato lines derived from a wild green-fruited species[END_REF], whereas most of fresh-market tomatoes contain less than 4% soluble solids. A negative relationship between fresh weight and soluble solids content in the fruit [START_REF] Goldenberg | Genetic and phenotypic correlation between weight and dry matter content of tomato fruits and their heritabilities[END_REF][START_REF] Ibarbia | Tomato fruit size and quality interrelationships[END_REF] hampers the transfer of high soluble-solids content from wild species into cultivated varieties. Several reviews summarize the numerous experiments of quantitative trait loci (QTL) mapping for tomato fresh weight and soluble solids content, and their co-localizations [START_REF] Grandillo | Identifying the loci responsible for natural variation in fruit size and shape in tomato[END_REF][START_REF] Causse | Traditional and enhanced breeding for quality traits in tomato[END_REF][START_REF] Foolad | Genome mapping and molecular breeding of tomato[END_REF]. QTL for these two traits are frequently co-localized with opposite allele effects.

Such QTL co-localizations can be due either to genetic linkage, and recombination could modify the relationships, or to fruit physiology and these relationships cannot be genetically modified. [START_REF] Yousef | Evaluation of breeding utility of a chromosomal segment from Lycopersicon chmielewskii that enhances cultivated tomato soluble solids[END_REF] identified chromosomal segments from S. chmielewskii that had a positive influence on fruit soluble solids while maintaining fruit size unaltered. Similarly, the alleles of S. pennellii at the QTL Lin5 for soluble solids content increases sugar content without reducing total yield [START_REF] Fridman | Zooming in on a quantitative trait for tomato yield using interspecific introgressions[END_REF].

In breeding programs, interactions between genotype and environment limit the possibility to increase sugar content. These interactions were studied in wheat through comparisons of various trial locations [START_REF] Robert | Structuring genotype x environment interaction for quality traits in bread wheat, in two multi-location series of trials[END_REF][START_REF] Matus-Cadiz | Genotype x environment interaction for grain color in hard white spring wheat[END_REF] or through impact of stresses (drought, salt) [START_REF] Snape | Dissecting gene x environmental effects on wheat yields via QTL and physiological analysis[END_REF][START_REF] Williams | The influences of genotype, environment, and genotype x environment interaction on wheat quality[END_REF]. In tomato, QTL detection under different saline conditions revealed a QTL specific to saline conditions [START_REF] Villalta | Comparative QTL analysis of salinity tolerance in terms of fruit yield using two solanum populations of F 7 lines[END_REF], indicating strong interactions between QTL and environment.

Understanding the relation between fruit fresh weight and sugar content in tomato fruits requires identifying the key processes underlying this relationship, both from the genetic and physiological points of view. The integration of physiological and genetic approaches allows a better dissection of the genetic basis of complex traits as well as their interactions with the environment. For instance, several co-localizations of QTL for fruit weight, fruit dry matter content or fruit sugar content and QTL for ecophysiological parameters were identified in peach [START_REF] Quilot | Analysing the genetic control of peach fruit quality through an ecophysiological model combined with a QTL approach[END_REF], suggesting which physiological processes were involved in the genetic variation.

Fruit weight and composition depend on the balance between inward and outward fluxes to/from fruit (mostly water and carbon), which involve many different processes. Transpiration leads to a water loss [START_REF] Wu | Relationship between skin speckle, soluble solids content and transpiration rate in nectarines[END_REF], and may decrease the fruit fresh weight and concentrate the soluble compounds. Cell division and cell expansion determine final fruit size and carbohydrate dilution within cells [START_REF] Bohner | Effects of fruit set sequence and defoliation on cell number, cell size and hormone levels of tomato fruits (Lycopersicon esculentum Mill.) within a truss[END_REF][START_REF] Ho | Tomato[END_REF]. The number of cells influences the structural dry matter through the amount of cell wall. The size of cells mainly affects the capacity to store soluble dry matter. Finally, the number of seeds may interfere with cell division and cell expansion through the production of hormones [START_REF] Rylski | Fruit set and development of seeded and seedless tomato fruits under diverse regimes of temperature and pollination[END_REF][START_REF] Gillaspy | Fruits : a developmental perspective[END_REF]. At the plant level, characteristics of the leaves (carbon source) could affect the sugar production via photosynthesis.

Carbon supply can be modified by environmental stresses or cultural practices. For instance, fruit thinning reduces the competition for carbon, and thus promotes fruit size and sugar content in several species, like peach fruits [START_REF] Morandi | Carbohydrate availability affects growth and metabolism in peach fruit[END_REF], apple fruits [START_REF] Link | Significance of flower and fruit thinning on fruit quality[END_REF], mandarin fruits [START_REF] Kubo | Sucrose accumulation and its related enzyme activities in the juice sacs of satsuma mandarin fruit from trees with different crop loads[END_REF], and papaya fruits [START_REF] Zhou | Defoliation and fruit removal effects on papaya fruit production, sugar accumulation, and sucrose metabolism[END_REF]. In tomato, this treatment has previously been carried out, increasing in similar proportions fresh and dry weight [START_REF] Heuvelink | Effect of fruit load on dry matter partitioning in tomato[END_REF].

The aims of this study were first to determine and dissect the relative influence of genotypes and fruit load on fruit weight and sugar content, then to identify which processes were underlying these traits and their relationships, at both genetic and physiological levels. For this purpose, a set of physiological, biochemical and morphological parameters related to the main processes of fruit growth and plant development were measured in a population consisting of 20 introgression lines derived from a cross between a wild tomato species (Solanum chmielewskii) with small green fruits and a medium sized cultivar (Solanum lycopersicum L.). QTL mapping and interactions between QTL and carbohydrate availability were examined over two years, by comparing plants under two contrasted fruit loads. One fruit per truss was the condition with no limitation for carbohydrate supply and corresponded to maximum genotypic potential. The second condition was similar to what is currently applied for genetic studies: a free load condition (without fruit pruning), with competition for assimilates among fruits.

Materials and methods

Plant material

The study was performed using the Solanum lycopersicum line 'Moneyberg' (hereafter referred to as M) and 20 indeterminate lines carrying single or multiple introgressions of the Solanum chmielewskii LA1840 in the background of Moneyberg, kindly provided by Keygene (The Netherlands). Each line was named by the chromosomal number and the location of the largest introgression. For instance, genotype C3a was the line that contained an introgressed fragment at the top of chromosome 3, while genotype C3d possessed an introgression at the bottom of chromosome 3.

Growth conditions and experimental treatments

Seeds were sown at the end of February, and a total of 400 plants were grown at a density of 3.6 plants m -2 in a ground bed greenhouse in Avignon (Southern France) at day-night temperature set points of 24/16°C during spring 2006 (March-July) and 25/15°C during spring 2007 (March-July).

Plants were randomly distributed in two blocks containing each 200 plants and facing respectively North and South. Plant nutrition and chemical pest and disease control followed commercial practices and plants were conducted on single vine. Starting from anthesis of the first truss, flowers were pollinated with an electrical shaker every two to three days.

For each genotype, ten plants were randomly selected in the first block while nine plants were randomly selected in the second block. On 12 plants of each genotype, trusses were pruned to one fruit (low fruit load, LL) while on seven other plants trusses were not pruned (high fruit load, HL).

Under HL condition, the average number of fruit sets per truss within the population was 5.3. On each inflorescence of the LL plants, all the flowers except the second one were removed just after fruit set. The fruit removal experiment concerned the first nine trusses. All the plants were stopped two leaves above the ninth truss.

Observations and measurements

Plant development

Anthesis time, achieved as flower fully opened, was recorded three times a week in order to determine fruit age and fruit development duration (Dura expressed in days) considered as time between anthesis and red ripe stage. Plant development traits, used as indicators of plant vigour and carbohydrate supply, were measured on nine contrasted genotypes (C11b, C12d, C3a, C3c, C4c, C4d, C7a, C8e and C9d) and M. For each genotype and treatment, four randomly selected plants (two per block) were measured for the number of leaves (LfN) until the ninth truss and the height of the fourth truss (H4t expressed in cm), at the end of the growing season. The area of five representative leaves were measured using a planimeter, and then the total leaf area (LfA expressed in cm²/plant) was calculated by multiplying the mean leaf area by the number of leaves.

Dry weight of the five leaves was assessed after five days in a ventilated oven at 80°C and the specific leaf weight (SLW expressed in g.cm -²) was calculated.

Fruit cuticular conductance and cracking

Fruit transpiration is an important process involved in fruit growth and water content. Fruit cracking and fruit cuticular conductance were thus measured. Fruit surface conductance to water vapor diffusion (CutC expressed in cm.h -1 ) was measured as described in [START_REF] Gibert | Modelling the effect of fruit growth on surface conductance to water vapour diffusion[END_REF], only on fruits harvested in 2007, and at 21 days post anthesis, as this stage corresponds to the visual absence of cuticular macro-cracks. Five fruits per genotype and per fruit load, located on the fifth truss (flower 2), were harvested on five plants randomly selected within the two blocks. Then, fruit cracking was estimated on five red ripe fruits, harvested for each genotype and treatment on five different plants randomly selected within the two blocks. In order to eliminate bias due to competition within and among trusses, only the second fruits of the fourth trusses were harvested.

On each fruit, cheek, suture and height diameters were measured using a calliper and the total cuticle area was calculated assuming spherical form. The cuticular macro-crack area (FCr expressed in percentage of cracked area relative to the total cuticle area) was calculated by drawing their outlines on tracing paper, cutting the area and measuring it with a planimeter.

Fruit composition, cell number and size

The five red ripe fruits analyzed for FCr were used for fruit composition, and cell measurements.

Fruit fresh weight was measured (FW expressed in g / fruit), jelly and seeds were removed from the fruits and seeds were counted (SdN). Then, pericarp tissue including external, internal and transverse parts was weighted. One half of the fruit pericarps was ground in liquid nitrogen using a blender. Powders were freeze-dried and stored at -20°C prior to sugar extraction. Pericarp dry weight (DW expressed in g) and dry matter content (DMC expressed in g / 100 g FW) were measured after lyophilization of tomato powders. Sugars were extracted from the pericarp according to the method described in [START_REF] Gomez | A new procedure for extraction and measurement of soluble sugars in ligneous plants[END_REF]. The main sugar contents (glucose, fructose and sucrose) were quantified by enzymatic assay in 96-well microplates, as detailed in [START_REF] Gomez | The microplate reader: an efficient tool for the separate enzymatic analysis of sugars in plant tissuesvalidation of a micro-method[END_REF]. Sugar contents were expressed relative to the pericarp fresh weight (SUGfw in g / 100 g FW) or to the pericarp dry weight (SUGdw in g / 100 g DW). The pericarp structural carbon content (StrCfw) was estimated as follow:

StrCfw = DMC -SUGfw [1]
where StrCfw, DMC and SUGfw are expressed in g / 100 g FW. This relationship was experimentally checked on 30 red ripe fruits, by the extraction of pericarp insoluble material after sugar eliminations from lyophilized powders.

On the second half of the fruit pericarp, cell division and expansion were evaluated by measuring the number (ClN) and mean size (ClS expressed in nL) of pericarp cells, according to the method described in [START_REF] Bertin | Number of cells in tomato fruit depending on fruit position and source-sink balance during plant development[END_REF].

DNA markers and assays

A set of PCR-based markers consisting of 130 Conserved Ortholog Set II (COSII) markers [START_REF] Wu | Combining bioinformatics and phylogenetics to identify large sets of single-copy orthologous genes (COSII) for comparative, evolutionary and systematic studies: A test case in the euasterid plant clade[END_REF]; http://www.sgn.cornell.edu/markers/cosii_markers.pl) and 3 Simple Sequence Repeats (SSR) [START_REF] Frary | Development of a set of PCR-based anchor markers encompassing the tomato genome and evaluation of their usefulness for genetics and breeding experiments[END_REF], covering all 12 tomato chromosomes, were used to genotype the population of 20 S. chmielewskii LA1840 introgression lines (IL).

Genomic DNA of the two parent lines and of the 20 IL was extracted from leaf tissue of three weeks old plants according to the protocol of [START_REF] Fulton | Microprep protocol for extraction of DNA from tomato and other herbaceous plants[END_REF].

The PCR of each marker was performed on both parents (Moneyberg and LA1840) in 25 l reactions containing 50 ng of template DNA, 2.5 pmol of each forward and reverse primer, 1X

Colorless GoTaq ® Flexi Buffer (Promega) and, 0.2 mM dNTPs, and 0.2 U GoTaq ® DNA Polymerase (Promega). The reactions were amplified using a DNA Engine (PTC-200) Peltier

Thermal Cycler (Biorad). Amplification consisted of an initial denaturation for 5 min at 94° C, followed by 35 cycles of amplification with denaturation at 94° C for 30 s, annealing at 55° C for 45 s, and extension at 72° C for 60 s, with a final cycle of 72° C for 5 min. Following amplification, PCR products were analyzed by electrophoresis on 1% agarose gels in 1X TAE buffer for 1-2 h at 100 V and room temperature. Direct fragment length polymorphism was detected for 12% of the COSII markers. In the other cases, the PCR products of the parent lines were digested with different frequent cutter restriction enzymes including TaqI, HinfI, AluI, DraI, RsaI, MseI and electrophoresed through 2% agarose to identify polymorphisms. If no polymorphism was detected with these enzymes, then the PCR products of the two parent lines were sequenced. For this purpose, amplicons showing a single band on agarose gels were cleaned up from nucleotides and residual primers using 10 μl of PCR product and 2 μl of ExoSAP-IT ® , and the mixture was incubated at 37°C for 1h and then at 80°C for 15 minutes. A 10 μl sequencing reaction volume was prepared using 2 μl of the cleaned PCR product, 1 μl of either the forward or reverse primer, 1 μl of 5X BigDye Buffer 3.1 and 1 μl BigDye terminator v3.1 (Applied Biosystems). The sequencing PCR consisted of initial denaturation at 96°C for 60s and 25 cycles of 10 seconds at 96°C, 5 sec at 55°C, and 4 min at 60°C. The sequences were obtained through ABI Prism 3100 Genetic Analyzer automated sequencer (Applied Biosystems). Restriction maps were predicted using CAPS designer (http://www.sgn.cornell.edu/tools/caps_designer).

For some of the regions where polymorphic COSII markers were not found, SSR were surveyed in order to increase the map saturation. PCR reactions and thermocycling were as described for the COSII markers; PCR products were separated on 2.5% agarose gels in 1X TAE buffer for 1-2 h at 100 V and room temperature.

Sequences of the primers used for the COSII and SSR assays are available on the SGN website (http://www.sgn.cornell.edu). The polymorphic COSII and SSR markers were then assayed on the 20 IL.

Statistical analyses

Altogether, data obtained in this study contained very few missing values, and all variables were normally distributed. After checking that block effect was not significant, the effects of year, genotype, fruit load and their interactions on each trait were analyzed using a three-way ANOVA.

The linear model used was:

ijkl jk ik ij k j i ijkl Y                 ) ( ) ( ) ( [2]
where Y ijkl is the trait value for plant l on genotype i, fruit load j and year k; is the general mean ;

 i is the effect of genotype i,  j is the effect of fruit load j,  k is the effect of year k; is the effect of the interaction between fruit load j and year k and  ijkl is the error term.

The percentage of variation due to each factor was then calculated as:

   i i i i SS SS V 100 [3]
where SS is the sum of squares of the factor i. For each trait t, the percentage of fruit load variation ( FL ) from HL to LL was calculated following:

        100     HL t HL t LL t t FL    [4]
where   HL t  and   LL t  are the general means of the trait t under HL and LL, respectively.

For each genotype, the effect of fruit load on fruit fresh weight (FW), pericarp dry weight (DW), pericarp dry matter content (DMC), and pericarp sugar contents (SUGfw and SUGdw), was analyzed by a two-way analysis of variance following this model:

jkl jk k j jkl Y           ) ( [5]
where Y jkl is the trait value for plant l, fruit load j and year k; is the general mean ;  j is the effect of fruit load j,  k is the effect of year k, jk ) (  is the effect of the interaction between fruit load j and year k and  jkl is the error term.

Multivariate analyses were conducted using the MANOVA function in R (http://www.rproject.org). Genotypic and residual covariances were estimated for each trait pair under both fruit loads by pooling the two years and all genotypes as explained in [START_REF] Holland | Estimating genotypic correlations and their standard errors using multivariate restricted maximum likelihood estimation with SAS Proc MIXED[END_REF]. Covariance matrices were then used to calculate genotypic and residual correlation coefficients (hereafter called inter and intra genotypic correlation coefficients).

QTL analysis was performed under each fruit load condition, first separately on both years by using a one-way ANOVA and then by pooling the two years in a two-way ANOVA. The withingenotype mean squares from these ANOVA were used to carry out Dunnett multiple comparison test [START_REF] Dunnett | Pairwise multiple comparisons in the homogeneous variance, unequal sample size case[END_REF] in SAS version 9.1.3 (SAS Institute, Inc.), to determine which genotypes were significantly different (at the 0.05 probability level) from the parent conferring the genetic background (Moneyberg), meaning which genotypes carried a QTL. The QTL effects are presented as percentages of difference from Moneyberg.

Results

Genotype and fruit load effect on fruit and plant trait variation

As individual sugar contents (glucose, fructose and sucrose contents) were highly correlated together and to total sugar content (r≥ 0.8, data not shown), only the total pericarp sugar content was considered in the following analysis. Analysis of variance showed that the genotype significantly influenced all traits (Table 1), and that fruit load significantly influenced most of the traits, except seed number, cell size, leaf number and specific leaf weight. Most of the traits showed higher average values under LL than under HL, except the sugar content (SUGdw), the cuticular conductance (CutC), and the fruit development duration (Dura). Fruit load mostly affected fruit cracking, pericarp dry weight and fruit weight which increased by about 1000%, 87% and 55%, respectively, under LL condition. Genotype x fruit load interactions were also significant for all traits except the seed number, the height of the 4 th truss and the leaf number. Year effects and interactions between genotype, fruit load and year were also found significant. Globally the percentages of variations due to interactions between year and genotype or between year and fruit load were lower than those of genotype or fruit load.

Fruit load effect was analyzed for each genotype on fruit fresh weight, pericarp dry weight, dry matter content and sugar contents (Fig. 1). Linear positive relationships were found between HL and LL conditions for fruit fresh weight, pericarp dry weight and dry matter content. For sugar contents, the correlation between HL and LL values was lower (SUGfw) or non significant (SUGdw). The fruit fresh weight and pericarp dry weight of all genotypes significantly increased from HL to LL, enlarging twice their range of variation (Fig. 1B and Fig. 1C) and the dry matter content (DMC) significantly increased for all the genotypes except for C3a. Five groups of genotypes were identified according to their response to low fruit load. The first group contained two genotypes and was characterized by significant increase in sugar contents (SUGdw and SUGfw). The second group was the most represented within this population as it contained 13 genotypes (including M): it was distinguished by increase in SUGfw but no change in SUGdw.

Groups 3, 4 and 5 (three, two and one genotypes, respectively) showed no change in SUGfw, associated with a decrease (Gr 3) or no change (Gr 4 and Gr 5) in SUGdw. The belonging to these groups was not dependent on particular values of fruit fresh weight or sugar contents.

Inter-and intra-genotypic correlations among traits

Genotypic and residual correlations among traits were evaluated by analyzing inter-and intragenotypic correlation coefficients among variables, respectively (Table 2). Some correlations were common to inter-and intra-genotypic levels and to both fruit loads. Fruit fresh weight was positively correlated to pericarp dry weight, seed number and cell number. Pericarp dry weight was positively correlated to cell number. Dry matter content and sugar content (SUGdw) were positively correlated to SUGfw and positively or negatively to the structural carbon content, respectively.

Only few and low correlations were specific to the intra-genotypic level (Table 2). More significant correlations occurred at the inter-genotypic level and one was common to HL and LL: the higher the structural carbon content, the lower the seed number. Antagonisms were found between fruit fresh weight and dry matter content or structural carbon content under HL. Dry matter content was positively related to fruit cracking and leaf number and negatively to the fruit development duration under HL. Dry matter content was also positively correlated to plant height under LL. The sugar content (SUGfw) was positively related to fruit cracking and leaf number under HL and to specific leaf weight under LL. It was also negatively linked to fruit development duration and to cell expansion via cell size under LL. Moreover, a significant negative correlation was found between cell number and cell size only under HL (r=-0.58, data not shown), suggesting a competition for assimilates among cells. Finally, the structural carbon content was negatively linked to fruit development duration and to cuticular conductance under HL while positively to leaf number under LL.

QTL for plant development, fruit weight and composition traits

Linkage map of introgressions

The S. chmielewskii LA1840 population provided by Keygene had been developed using AFLP markers (http://www.keygene.com, unpublished data). In order to enhance the rate of introgression breeding, facilitate marker-assisted selection of new IL and comparisons between function maps of tomato and potato, within the framework of a large European project (http://www.eu-sol.net/), the S. chmielewskii IL population along with other four interspecific tomato mapping populations and one potato mapping population are being anchored to a common set of COSII markers. The other five mapping populations are: a potato diploid population F1840 [START_REF] Gebhardt | Comparative mapping between potato (Solanum tuberosum) and Arabidopsis thaliana reveals structurally conserved domains and ancient duplications in the potato genome[END_REF], the S. pennellii LA716 IL [START_REF] Eshed | An introgression line population of Lycopersicon pennellii in the cultivated tomato enables the identification and fine mapping of yield-associated QTL[END_REF], the S. habrochaites LA1777 IL [START_REF] Monforte | Development of a set of near isogenic and backcross recombinant inbred lines containing most of the Lycopersicon hirsutum genome in a L. esculentum genetic background : a tool for gene mapping and gene discovery[END_REF], the S. neorickii LA2133 backcross inbred lines (BIL) [START_REF] Fulton | Advanced backcross QTL analysis of a Lycopersicon esculentum x Lycopersicon parviflorum cross[END_REF], and the S. cheesmaniae LA483 recombinant inbred lines (RIL) [START_REF] Paran | Recombinant inbred lines for genetic mapping in tomato[END_REF].

The S. chmielewskii IL population has been anchored to 130 COSII markers and 3 SSR spanning a total of 1459 cM with an average interval of 11 cM between adjacent markers. One hundred and five COSII markers are common to at least four of the six mapping populations. Grandillo et al. were thus allowed to make comparisons among QTL mapping studies (unpublished data).

Fourteen genotypes (C1a, C3a, C3d, C4c, C4d, C6e, C7a, C7b, C7d, C8a, C8c, C9d, C10b, C12d) were identified as carrying a single S. chmielewskii introgression, covering around 30% of the genome, three genotypes carried two introgressions (C3c, C8e, C11b), one genotype carried three introgressions (C9a) and two genotypes carried four introgressions (C5b, C9c) (Fig. 2). Altogether, the 20 IL covered approximately 50% of the S. chmielewskii genome. The genetic length of each introgression varied from 4% to 45% of the chromosome length. Some of the lines still contained heterozygous fragments, either as single introgression (C4c) or as additional introgressions, on chromosomes 2 (C3c), 4 (C5b), and 11 (C9a, C9c).

Quantitative trait locus analysis

Tables 3, 4 and 5 present QTL detected separately under each fruit load condition, and according to the model based on the two-year experiment. Positive or negative QTL correspond to a location where the allele of S. chmielewskii increased or decreased the trait, respectively, compared to Moneyberg. Excluding QTL for cuticular conductance which were only identified in 2007, 84 QTL were detected in 2006 versus 90 in 2007. Sixty-one QTL were common to both years. Sixty eight QTL were detected under HL versus 74 under LL. Only 30% of these QTL were detected whatever the fruit load and hereafter called 'stable'. When stable, the sign of the QTL was the same under both fruit loads, except for QTL of cuticular conductance.

QTL for fruit weight and composition traits

Fourteen QTL were detected for FW, half being stable, the other half being detected only under LL (Table 3). All stable QTL were negative, with two genotypes (C11b and C3c) having a strong effect (fruits were 50% smaller than M). A multiple introgressed genotype (C5b) carried a positive QTL for FW under LL. Ten QTL were identified for DMC, all with positive effect and six were stable. Concerning sugar contents, six QTL for SUGfw were detected: two QTL were stable, and only two had negative effects. On the contrary, all the 11 QTL identified for SUGdw had low negative effects. They were twice more numerous under LL than under HL and C11b was the only genotype carrying a stable QTL whatever the fruit load.

QTL for fruit physiological and plant developmental traits

Five QTL with low effects were detected for fruit development duration and three of them were stable under both fruit loads (Table 4). Eight negative QTL for seed number were identified under HL, and half of them were not detected under LL. QTL for cell traits (number and size) were mainly detected under HL (nine QTL under HL against three under LL), and a QTL for cell number was stable (C11b). All QTL for cell number had negative effects while for cell size, two positive and one negative QTL were identified. The cuticular conductance is the trait for which the number of detected QTL was the largest, but also the trait for which a surprising behaviour was observed between the two fruit loads. Under HL, all QTL had negative effect (except C9c), while under LL all QTL had positive effect, and nine stable QTL had opposite effects under HL and LL.

This could be explained by the fact that Moneyberg had one of the highest values for conductance under HL, and the lowest under LL (data not shown). Concerning fruit cracking, C4d carried a stable positive QTL. The five other QTL were only detected under LL and had negative effects.

Only few QTL were detected for plant development traits (from one to five QTL per trait) as QTL analysis was only carried out on ten genotypes (Table 5). Among them, less than the half were stable.

Co-localizations between QTL for fruit weight and composition and QTL for fruit physiology or plant development

One objective of the present study was to look for co-localizations between QTL for fruit fresh weight and composition and QTL for physiological parameters under both fruit load conditions. In order to avoid additive or epistatic effects due to multiple introgressions, only QTL colocalizations found on genotypes carrying a single introgression were shown on Fig. 3. Eleven clusters of QTL for fruit fresh weight and composition traits were found. Five regions carrying QTL for fruit weight with negative effects of S. chmielewskii alleles carried also QTL for dry matter content with an opposite allele effect, in accordance with the inter-genotypic correlation. In the same way, three co-localizations between QTL for fruit fresh weight and sugar content (SUGfw) were identified with opposite allele effects, while two others were identified with the same negative allele effects. Finally, three co-localizations were identified between QTL for fruit fresh weight and QTL for sugar content (SUGdw), with the same negative allele effects. Several co-localizations between QTL for fruit fresh weight and QTL for cell number or for seed number (both 3 clusters) were detected with same sign effects. Co-localizations between QTL for sugar content (SUGfw) and QTL for cell size were expected according to correlations. A region located on the bottom of chromosome 9 carried such co-localization with opposite allele effects (C9d). All QTL for fruit cracking co-localized with QTL for fruit fresh weight. Three of them had the same allele effects except on chromosome 4. Some co-localizations between QTL for dry matter content and QTL for plant development did not always correspond to the direction of correlations (Table 2). For example, dry matter content was positively correlated to plant height at the inter-genotypic level and QTL analysis revealed one colocalization for these two traits with similar allele effects (C4d), and two others with opposite allele effects (C3a and C9d). Some co-localizations were also identified between QTL for sugar content (SUGfw) and QTL for specific leaf weight or leaf number (C3a and C4d).

Discussion

Relationships between fruit fresh weight and composition

The influence of carbon availability on the relationships between fruit weight and sugar content was studied via fruit thinning. In peach fruits [START_REF] Morandi | Carbohydrate availability affects growth and metabolism in peach fruit[END_REF] or in apple fruits [START_REF] Link | Significance of flower and fruit thinning on fruit quality[END_REF], a lower fruit load increased simultaneously fruit fresh weight and sugar content. In tomato, some fruit thinning experiments were carried out and showed that fruit load reduction led to increase fruit fresh and dry weights [START_REF] Gautier | Modulation of competition between fruits and leaves by flower pruning and water fogging, and consequences on tomato leaf and fuit growth[END_REF][START_REF] Bertin | Analysis of the tomato fruit growth response to temperature and plant fruit load in relation to cell division, cell expansion and DNA endoreduplication[END_REF][START_REF] Baldet | The expression of cell proliferation-related genes in early developing flowers is affected by a fruit load reduction in tomato plants[END_REF], but none of them has dealt with effects on fruit sugar content. In this tomato population, most of the genotypes (15 genotypes over 21) reacted to fruit thinning by increasing the dry matter content and, for only two of them, by also increasing the carbohydrate allocation to sugar metabolism. For the six other genotypes, sugar content relative to fresh weight was unchanged whatever the fruit load. Finally, in some cases, even if the dry matter content increased, the carbohydrate allocation to sugars was stable or decreased, leading to no change in sugar content. The only significant correlation between fruit fresh weight and composition was found at the inter-genotypic level, under high fruit load conditions (Table 2) and concerned the antagonism between fruit fresh weight and dry matter content, associated with a negative relation between fruit fresh weight and structural carbon content, but not between fruit fresh weight and sugar content. From a genetic point of view, antagonism between fruit fresh weight and soluble solids content has been observed during tomato improvement under high load conditions, and is mainly due to co-localizations of QTL with antagonistic effects [START_REF] Bernacchi | Advanced backcross QTL analysis in tomato. I. Identification of QTLs for traits of agronomic importance from Lycopersicon hirsutum[END_REF][START_REF] Chen | Mapping of QTLs for lycopene and other fruit traits in a Lycopersicon esculentum x L. pimpinellifolium cross and comparison of QTLs across tomato species[END_REF][START_REF] Saliba-Colombani | Genetic analysis of organoleptic quality in fresh market tomato. 1. Mapping QTLs for physical and chemical traits[END_REF][START_REF] Lecomte | Fine mapping of QTLs of chromosome 2 affecting the fruit architecture and composition of tomato[END_REF]). In the present population, markers used to locate the introgressed fragments were common to four tomato mapping populations. By this way, it was possible to compare QTL from the S. chmielewskii population with the previous works. Co-localizations between QTL for fruit fresh weight and QTL for dry matter content or sugar content on fresh weight basis were detected with opposite effects (Fig. 3), and corresponded to regions that were already identified in other progenies under high fruit load [START_REF] Grandillo | Identifying the loci responsible for natural variation in fruit size and shape in tomato[END_REF][START_REF] Causse | Traditional and enhanced breeding for quality traits in tomato[END_REF]. These colocalizations of QTL for fruit fresh weight with QTL for dry matter content with antagonistic effects were common to both fruit loads (C9d), specific to high load (C3a) or specific to low load (C3d, C4d and C8a) in the S. chmielewskii population. Moreover, two co-localizations were identified between QTL for fruit fresh weight and QTL for sugar content relative to fresh weight, with same negative allele effects under high load (C10b) and under low load (C12d). These regions could be involved in carbon allocation to cell structures or to sugar metabolism, because they also co-localized with QTL for sugar content relative to dry weight. In this population, we did not find any stable QTL with positive effects for sugar content which was not associated with a negative QTL for fruit fresh weight.

Processes underlying fruit weight and composition

Cuticle properties are involved in water loss by transpiration [START_REF] Schönherr | Water permeability of isolated cuticular membranes: the effect of cuticular waxes on diffusion of water[END_REF][START_REF] Becker | Water permeability of plant cuticles: permeance, diffusion and partition coefficients[END_REF][START_REF] Kerstiens | Cuticular water permeability and its physiological significance[END_REF]. Concerning fruit cracking, it usually occurs when fruit growth rate is high [START_REF] Christensen | Cracking in cherries: VI. Cracking susceptibility in relation to the growth rythm of the fruit[END_REF]. It was thus expected that under low fruit load, when fruit growth was faster, QTL for fruit fresh weight co-localized with QTL for fruit cracking. Genotype C4d showed an extreme phenotype even under HL and carried a major mutation for this trait. It probably corresponded to the mutation Cwp1 described by [START_REF] Hovav | The identification of a gene (Cwp1), silenced during Solanum evolution, which causes cuticle microfissuring and dehydration when expressed in tomato fruit[END_REF]. Three other QTL for fruit cracking were mapped and co-localized with QTL for fruit fresh weight with similar effects (Fig. 3), while no such link was found between fruit cracking and fruit fresh weight through correlation analysis (Table 2). An increase in cuticular conductance or in fruit cracking was expected to increase the dry matter content or the sugar content (SUGfw). The hypothesis of a concentration of dry matter by water loss was confirmed at the inter-genotypic level under HL condition, as fruit cracking was positively correlated to dry matter content or sugar content. To avoid interaction with fruit cracking, measurement of cuticular conductance was performed 21 days after anthesis, during fruit cellular expansion, whereas all other variables were analyzed at the red ripe stage. As cuticular conductance is known to decrease throughout fruit development [START_REF] Gibert | Modelling the effect of fruit growth on surface conductance to water vapour diffusion[END_REF], conclusions may be biased. Numerous QTL for cuticular conductance were particularly susceptible to fruit load change as S. chmielewskii alleles carried negative effects under HL and positive effects under LL. This is probably linked to the extreme behaviour of Moneyberg changing from HL to LL. No other trait behaved similarly, making the interpretation of these co-localizations difficult.

Cell number is a determinant factor of fruit sink strength, usually determined during early stages of tomato fruit development [START_REF] Bohner | Effects of fruit set sequence and defoliation on cell number, cell size and hormone levels of tomato fruits (Lycopersicon esculentum Mill.) within a truss[END_REF][START_REF] Joubes | Molecular and biochemical characterization of the involvement of cyclin-dependent kinase A during the early development of tomato fruit[END_REF]. In agreement with the literature [START_REF] Higashi | Histological analysis of fruit development between two melon (Cucumis melo L. reticulatus) genotypes setting a different size of fruit[END_REF][START_REF] Jullien | Effect of pulp cell number and assimilate availability on dry matter accumulation rate in a banana fruit [ Musa sp. AAA group 'Grande Naine' (Cavendish subgroup)[END_REF][START_REF] Bertin | Analysis of the tomato fruit growth response to temperature and plant fruit load in relation to cell division, cell expansion and DNA endoreduplication[END_REF], the fruit fresh weight was positively correlated with cell number and was present both at intra-genotypic and inter-genotypic levels whatever the fruit load (Table 2). This relationship was also confirmed by the numerous colocalizations under both fruit loads between QTL for fruit fresh weight and QTL for cell number, with similar effects (Fig. 3). The increase in fruit fresh weight and cell division in the absence of carbohydrate competition has also been described by [START_REF] Baldet | The expression of cell proliferation-related genes in early developing flowers is affected by a fruit load reduction in tomato plants[END_REF] and was linked to the regulation of fw2.2, a cell cycle-control gene. Cell expansion was not directly linked to the fruit fresh weight as no significant correlation was found between cell size and fruit fresh weight and QTL co-localizations of these two traits were either with similar allele effects (C9d) or with opposite effects (C12d). Moreover, at the inter-genotypic level and under high fruit load condition, a negative correlation between cell number and cell size suggested the existence of competition for carbohydrates among cells. This competition may reduce the cell growth potential, resulting in smaller cells [START_REF] Bertin | Analysis of the tomato fruit growth response to temperature and plant fruit load in relation to cell division, cell expansion and DNA endoreduplication[END_REF][START_REF] Tsukaya | Mechanism of leaf-shape determination[END_REF].

Positive correlations between fruit fresh weight and seed number were probably due to the effect of the latter on the sink strength of the fruit [START_REF] Nitsch | Hormonal factors in growth and development[END_REF]. QTL co-localizations of these two traits were similar to previous QTL mapping experiments carried out on a S. pimpinellifolium F 2 population (Van der [START_REF] Van Der Knaap | The making of a bell pepper-shaped tomato fruit: identification of loci controlling fruit morphology in Yellow Stuffer tomato[END_REF] or on a S. peruvianum BC 3 population [START_REF] Fulton | QTL analysis of an advanced backcross of Lycopersicon peruvianum to the cultivated tomato and comparisons with QTLs found in other wild species[END_REF]. A co-localization was common to the S. pimpinellifolium F 2 population and to the present population at the bottom of chromosome 6. The role of seeds in the sink strength could be related to hormonal signalling, hormones being strongly implicated in the control of cell division rate and sustenance [START_REF] Gillaspy | Fruits : a developmental perspective[END_REF].

Plant development traits are commonly described in ecophysiological studies, and were used here to estimate source strength. The main hypothesis was that the higher the number of leaves, or the leaf area, the higher the sugar production via photosynthesis, and thus the higher the fruit dry matter or sugar contents. Positive correlations between leaf number and dry matter or sugar content (SUGfw) comforted this hypothesis. A balance between plant vigour and distribution of assimilates to the fruit was thus found at the whole plant level, as already mentioned by [START_REF] Vaast | Fruit thinning and shade improve bean characteristics and beverage quality of coffee (Coffea arabica L.) under optimal conditions[END_REF]. The hypothesis was also corroborated at the genetic level, as co-localizations between QTL for leaf area or leaf number and QTL for dry matter content or sugar content were identified with same allele effects (C4d and C9d).

Conclusion

The present paper aimed at assessing the implication of various processes in the relationships between fruit weight and its composition. Our results suggested that these relationships could be mainly related to sink strength through cell division whose intensity was modulated by fruit load.

An antagonism between fruit fresh weight and dry matter content was only detected at the intergenotypic level, in condition of competition for assimilates. This study also revealed different behaviours of genotypes with respect to changes in fruit load and it was consequently not possible to deduce a general scenario for the whole population. Moreover, a lot of QTL had different effects according to fruit load, suggesting that carbohydrate supply can strongly interact with the genome, probably via sugar or hormonal sensing. Although co-localizations of QTL can hide either pleiotropy or genetic linkage, this work could contribute to help in choice of candidate genes as physiological hypotheses linked to quality traits were formulated at the genetic level. QTL with similar behaviours under both fruit loads could also be interesting targets for breeding programs as they would be more likely stable under various environments. 25 1 Table 1: Percentage of variation attributable to the effects of genotype (G), fruit load (FL), year (Y), genotype x fruit load interaction (G x FL), 2 genotype x year interaction (G x Y), year x fruit load interaction (Y x FL) by analysis of variance. When fruit load effect was significant, the 3 percentage of variation from high load (HL) to low load (LL) ( FL ) was calculated (equation 4 in Material and Methods). 4  FL / / 5 'ns' : not significant; *: significant at the 0.05 probability level; **: significant at the 0.001 probability level; ***: significant at the 0.0001 probability level (http://www.sgn.cornell.edu/). At the bottom of each chromosome, the total genetic length is mentioned. If the genotype carried a single homozygous introgression, its location is black-filled (for example C1a). If the genotype carried a single heterozygous introgression, its location is grey-filled (C4c). If the genotype carried multiple homozygous introgressions, its name referred to the largest introgression and the locations of all the introgressions are white-filled. The locations of the other introgressions are indicated at the bottom of the fragment. For example, the C3c carried its main introgression on chromosome 3 and a smaller one on chromosome 2. If the genotype carried multiple heterozygous introgressions, introgressions are hatched. chmielewskii alleles had negative or positive effects on the trait, respectively.

Fruit weight and composition

  of the interaction between genotype i and fruit load j ik ) ( is the effect of the interaction between genotype i and year k; jk ) ( 

Figure 1 :

 1 Figure 1: Effect of fruit load on several traits.

Figure 2 :

 2 Figure 2: Linkage map showing the locations of the 110 COS markers whose names begin with

Figure 3 :

 3 Figure 3: Genetic map of the QTL detected on genotypes carrying a single introgressed

Table 2 : Inter-genotypic (Inter) and Intra-genotypic (Intra) correlation coefficients for traits relative to weight, composition, physiology of the fruit, and plant development, under high load (HL) and low load (LL).

 2 Correlation coefficients significant at the 0.05 probability level were highlighted in grey.

	HL	LL	HL	LL	HL	LL	HL	LL	HL	LL	HL	LL	HL	LL	HL	LL	HL	LL	HL	LL	HL	LL	HL	LL

Table 4 : QTL characteristics for fruit physiological traits.

 4 Genotypes for which significant differences were detected at the 0.05 probability level, in the model taking into account the two-year experiment. Chromosomes carrying introgressions are indicated for each genotype. Effects under HL or under LL are expressed as the average percentage of difference between the genotype and Moneyberg under high load and low load, respectively, over two years. means that the QTL was not significant. 2 'na' means that QTL stability cannot be deduced as the cuticular conductance was only measured in 3

	Under each fruit load, the year when QTL was significant was indicated as '06 / 07' when the QTL
	was detected whatever the year, '06' or '07' when it was significant only in 2006 or 2007,
	respectively.

Acknowledgements

We are grateful to the greenhouse experimental team and to Yolande Carretero for taking care of the plants. We thank Jean-Claude L'Hotel and Michel Pradier for their technical support during harvests. Many thanks to Beatrice Brunel and Esther Pelpoir for managing cell counting and to Emilie Rubio for sugar analyses. Many thanks also to Rebecca Stevens for English revising.

Keygene, The Netherlands is acknowledged for providing seeds of the tomato population. This work was funded by the European EU-SOL Project PL 016214-2 and Marion Prudent was supported by a grant from INRA and Région Provence Alpes Côte d'Azur (France).

* : traits for which correlation coefficients were calculated only among 10 genotypes Table 3: QTL characteristics for fruit weight and composition traits. 1 Genotypes for which significant differences were detected at the 0.05 probability level, in the model 2 taking into account the two-year experiment. Chromosomes carrying introgressions are indicated for 3 each genotype. Effects under HL or under LL are expressed as the average percentage of difference 4 between the genotype and Moneyberg under high load and low load, respectively, over two years. 5

Under each fruit load, the year when QTL was significant was indicated as '06 / 07' when the QTL 6 was detected whatever the year, '06' or '07 '