Different genetic architectures underlie crop responses to the same pathogen: the {Helianthus annuus * Phoma macdonaldii} interaction case for black stem disease and premature ripening

Background Phoma macdonaldii has been reported as the causal agent of black stem disease (BS) and premature ripening (PR) on sunflower. PR is considered as the most widespread and detrimental disease on sunflower in France. While genetic variability and QTL mapping for partial resistance of sunflower to stem, collar and roots attacks have been reported on plantlets in controlled conditions, this work aims to describe the genetic variability in a subset of a sunflower lines, and for the first time to map QTL involved in PR resistance evaluated in field conditions using controlled inoculation. Results An efficient and reliable method for inoculation used in field experiments induced stem base necrosis on up to 98% of all plants. A significant genetic variability for PR resistance in the field was detected among the 20 inbred lines of the core collection tested across the two years. For QTL mapping, the PR resistance evaluation was performed on two recombinant inbred lines (RIL) populations derived from the crosses XRQxPSC8 and FUxPAZ2 in two different years. QTL analyses were based on a newly developed consensus genetic map comprising 1007 non-redundant molecular markers. In each of the two RIL populations, different QTL involved in PR partial sunflower resistance were detected. The most significant QTL were detected 49 days post infection (DPI) on LG10 (LOD 7.7) and on LG7 (LOD 12.1) in the XRQxPSC8 and FUxPAZ2 RIL population, respectively. In addition, different QTL were detected on both populations for PR resistance measured between 14 and 35 DPI. In parallel, the incidence of natural attack of P. macdonaldii resulting in BS disease was recorded, showing that in these populations, the genetic of resistance to both diseases is not governed by the same factors. Conclusion This work provides the first insights on the genetic architecture of sunflower PR resistance in the field. Moreover, the separate studies of symptoms on different organs and in time series allowed the identification of a succession of genetic components involved in the sunflower resistance to PR and BS diseases caused by Phoma macdonaldii along the development of the {plant * pathogen} interaction. Electronic supplementary material The online version of this article (10.1186/s12870-017-1116-1) contains supplementary material, which is available to authorized users.


Background
The sunflower crop is faced to several diseases caused by fungi and oomycetes in all the regions where it is cultivated. The use of resistant sunflower varieties is an efficient way to control the diseases, and resistance to these diseases remains a major target for sunflower breeding [1]. During the last two decades, the premature death [2] or premature ripening (PR) induced by Phoma macdonaldii Boerema (teleomorph: Leptosphaeria lindquistii) became the most severe and widespread sunflower disease in France, and could be partly responsible for the yield stagnation around 2.5 t.ha −1 . P. macdonaldii is also responsible for black stem disease (BS). However, the damages of PR on seed yield appeared greater than those of black stem [3]. As the chemical control of the disease remains difficult, potentially dangerous for the environment and becomes less socially acceptable, the development of even partially resistant varieties is an important breeding objective. Some genetic variability has been described for sunflower resistance to BS disease [4,5] and PR [3]. Several studies dealing with the genetic control of sunflower resistance to Phoma attacks on petiole stem base and roots of seedlings in growth chamber have demonstrated the quantitative character of this resistance [6][7][8][9][10][11]. Cytological observations of a susceptible and of a more resistant inbred line showed that the development of fungal hyphae within the stele was affected in the more resistant genotype, suggesting the involvement of a plant compound in the defense [12]. In addition, candidate genes could be identified through transcriptomic studies of the sunflower * P. macdonaldii interaction (e.g. sunflower-likelipase, MYB-related transcription factor regulating PAL2, a key enzyme involved in the phenylpropanoid pathway) [13]. However, in these studies, the pathosystem was developed on two-leaf-stage plantlets that is not fully representative of natural attacks. In field or in greenhouse, P. macdonaldii has been found to cause higher damages when contamination occurs at the star bud phenological stage (E1) than at earlier phenological stages [14]. Recently, we demonstrated the clear role of aerial Phoma infection in PR compared with soilborne inoculums [15]. In addition, we showed that artificial inoculation at the stem base with pycniospores or mycelium of P. macdonaldii could be used for screening genotypes showing a substantial level of resistance to PR. Using this phenotyping protocol in field, we investigated in this study the phenotypical variability and the genetic architecture of premature ripening and black stem resistance along plant development in Helianthus.

Sunflower genetic material
In 2009, the resistance level to PR of 42 genotypes was evaluated in a field trial in Auzeville-Tolosane: 40 lines were selected from the Helianthus annuus core collection of 48 lines [16] and two additional lines were included due to their high level of resistance in 2007 preliminary observations: (Tub-1709-1)-1-6A is derived by selfing from the USDA accession TUB-1709-1 which results from an introgression of H. tuberosus [17], and 97B7 is an INRA line derived from a cross involving H. argophyllus. According to the results, a representative subset of 21 lines was chosen to confirm their resistance level in a second trial in 2010 (Table 1). Two sets of RIL were used for QTL mapping: a) a subset of 117 F8 lines from the "INEDI" RIL population which was obtained by single seeds descent from a cross between the lines XRQ and PSC8 [18,19], b) a subset of 113 F7-F10 lines from the "FUxPAZ2" RIL population, derived by single seed descent from a cross between the lines FU and PAZ2 [8,19]. These two RIL populations were chosen because of the difference for the resistance level to premature ripening of the parental lines. Indeed, in a 2008 previous trial, the XRQ and PSC8 lines were evaluated in the field using the same protocol and the same strain of P. macdonaldii than those used in this work: at phenological stage M1.2-M1.3 (43 days after contamination), 27% of plants were prematurely ripened for XRQ against 77% for PSC8. In the core-collection trial in 2010, FU was one of the most susceptible lines to premature ripening (PR AUDPC = 13.73) and PAZ2 one of the most resistant (PR AUDPC = 1.61). The four parental lines confirmed therefore their differences and the corresponding RIL populations were chosen for these experiments in 2010 and 2011. If XRQ and PAZ2 appeared to be more resistant than PSC8 and FU, they don't represent extreme behaviors towards P. macdonaldii. The availability of recombinant lines was also a criterion for choosing this material.

Experimental design
Four main field experiments were carried out to evaluate the resistance level of sunflower lines, from the three genetic designs, against PR over 3 years (2009,2010 and 2011) at INRA, Auzeville-Tolosane, near Toulouse (Haute Garonne, South-West, France). Before sowing, N fertilization was applied on each trial (60 kg.ha −1 ). The crop was sown at the beginning of May in 2009 and at the beginning of April in 2010 and 2011. The experiments were conducted in a randomized block design, with two replications in 2009 and three replications in 2010 for core collection trials, and with two replications for XRQxPSC8 (2010) and FUxPAZ2 (2011) RIL populations trials. In these two last experiments, the parental lines of the RIL populations were included as controls. Each plot consisted of 3 rows of about 25 plants. Plant density was 7 plants.m −2 after thinning. According to rainfall, irrigation (30 mm) was performed before inoculation in 2011. As Phomopsis helianthi is regularly present in Auzeville-Tolosane, Corbel (fenpropimorph, 0.8 l ha −1 , BASF) was applied in June to control specifically this disease in all experiments. This fungicide does not affect the development of Phoma macdonaldii.

Artificial plant inoculation
A single P. macdonaldii monopycniospore strain (MPH2) was used in all experiments. It was isolated from natural severe stem base lesions observed in 2006 in a sunflower crop located in Montgaillard-Lauragais (Haute-Garonne, France) and selected for its severe aggressiveness on a commercial cultivar Heliasol RM (KWS AG) susceptible to P. macdonaldii. Mycelium conservation, inoculum production and mycelium inoculation were performed as previously described [15]. Inoculum was produced on potato dextrose agar (PDA Difco 39 g.l −1 ,150 mg of streptomycin, pH 6) and grown at 25°C for 10 days in the dark. Mycelium inoculation was carried out at bud stage (E1 to E5) on 15 uniform plants per plot. A 6 mm diameter disk of PDA with mycelium was placed at the stem base of each plant and immediately covered with a damp cotton and an aluminum foil to prevent dehydration and left for 5 days. No natural attack at the stem base was observed before artificial inoculation. Table 1 Description of the 21 lines tested for their resistance to Premature Ripening. 'Core collection ID' indicates the line code used in this study and as previously used [16], 'Type' indicates if lines restore the male sterility on PET1 cytoplasm (R) or maintain the PET1 cytoplasmic sterility (B), 'Origin' indicates the line pedigree

Disease assessment
Disease symptoms were observed weekly on the 15 plants per replication, up to physiological maturity on the four experiments. Development of necrosis at the stem base and premature ripening (PR) induced by the P. macdonaldii artificial inoculation were assessed in each experiment from 14 days post-inoculation (DPI) to 56 or 63 DPI with an assessment every week. The disease was scored using a 0-4 scale: 0 = healthy plant, 1 = less than 3/4 of the stem base circumference black, 2 = necrosis girdling the stem base, 3 = all leaves wilted but the stem green, 4 = plant completely dry (Fig. 1). A plant was defined as reaching the PR stage when it is completely dry before physiological maturity, with necrosis girdling the stem base.
As natural infection caused by P. macdonaldii occurred from adjacent wheat fields and induced BS symptoms on the four trials, a disease assessment was also performed on BS disease, using a 0-2 scale: 0 = healthy plant, 1 = plant with isolated spots on the stem, 2 = plants with coalescent spots, and according to the same timing than for PR assessment (once a week, at the same time than PR assessment). These symptoms are expected to be induced by primary inoculum, as the secondary cycle of de P. macdonaldii was never observed on the French territory, contrary to what is described by other authors [20,21].
The developmental stage of each genotype was recorded as previously described [22] at each disease assessment.
All plants affected during the cropping season by other fungal diseases (Phomopsis stem canker, Verticillium wilt, Alternaria leaf spot and blight) were excluded from the PR symptoms assessment: the PR symptoms were recorded only on the remaining stem base infected plants in order to avoid any confusion between senescence due to Alternaria sp. or due to P. macdonaldii. Alternaria leaf spot and blight was assessed in July 2010 on XRQxPSC8 RIL population (two replications) using a 0-3 scale: 0 = healthy plant, 1 = foliar symptoms in the lower part of the plant, 2 = foliar symptoms on the whole plant, 3 = blighted plant.
In order to avoid any artifact due to the experimenter effect, a same experimenter never carried out two successive disease assessments on the same plot.

Statistical analyses on disease assessment data and other phenotypic traits
For each experiment devoted to disease assessment and each disease assessment date or developmental stage, a mean disease score of each plot and each genotype was calculated and the genotype effect on BS and PR resistance was assessed by analysis of variance on mean disease scores according to a general linear model (GLM procedure, SAS software, SAS Institute Inc.).
The BS and PR Area Under Disease Progress Curves (AUDPC) are estimated as previously described [23]: where y i = mean BS disease score or percentage of premature ripened plants at the i th observation, t = time (days) after inoculation at the i th observation, and n = total number of observations. Broad-sense heritability was estimated according to the following formula: where σ 2 G is the genetic variance (MSg -MSgr)/rn, and σe 2 is the environmental variance (MSe), n is the number of plants, and r is the number of replicates.

QTL analysis
Using the previous genotypic data and the methodology previously described [19], the two RIL populations were genotyped with a complementary AXIOM array (Affymetrix, USA) composed of 197,863 SNPs. For the "INEDI" population, we used a set of 32,666 SNPs that were polymorphic (with no segregating distortion in the whole population) and showed the highest quality ("PolyHighResolution") according to the automatic allele calling from Axiom Analysis Suite software. In the same manner, a set of 28,529 SNPs from the FUxPAZ2 RILpopulation was used for mapping. A set of SSR markers and additional SNPs from candidate genes sequencing Fig. 1 Correspondence between the scoring of Phoma macdonaldii PR disease and the symptoms observed on sunflower (a: score 0; b: score 1; c: score 2; d: score 3; e: score 4) was added to the "AXIOM" markers. A consensus map was then built with the CARTHAGENE software (http:// www.inra.fr/mia/T/CartaGene/, [24]) in using mergen command. In the process of obtaining robust consensus map, we used a) the markers mapped in [19], b) and from the AXIOM array, only the markers being genotyped in both populations. During the map building, only the markers with a LOD 2 points greater than 10 were selected.
For each RIL population, two different maps derived from this consensus map were used to perform QTL detection with MCQTL [25]. The phenotypic traits were recorded in the following independent trials: disease traits with artificial (PR) or natural (BS, Alternaria) infection for the "INEDI" RIL population (2010), diseases traits with artificial (PR) or natural (BS) infection for the FUxPAZ2 RIL population (2011).
A threshold corresponding to a Type I error rate of 1% at the genomewide level was used, as determined after 3000 replications of the resampling process for each trait and averaging the limit value for QTL detection across the different traits. The supports of the QTL were determined using the software MCQTL [25] with the QTL LOD value minus 2, which is expected to provide a 95% confidence. This procedure is producing longer interval supports that often reported in the literature (for example [9,10]).

Stem base necrosis and PR
In the four experiments, the mean percentage of plants showing Phoma necrosis at the stem base (disease score ≥ 1) at 14 DPI ranged from 97.9% to 98.7%. These high percentages showed the effectiveness of the inoculation method. Similar results were previously observed in field and greenhouse trials, both on adult plants [3,15].
In all four experiments, stem base necrosis increased from disease score 1 to disease score 2 during the first 35 DPI. All genotypes presented at least one plant with a girdling Phoma necrosis at the stem base (disease score = 2). The first premature ripened plants (disease score = 4) appeared in each trial at 42 DPI and between 17% and 33% of plants showed PR at 49 DPI according to the experiments. One week later, this percentage ranged between 35% to 51%. This evolution is consistent with previous results obtained on two commercial sunflower hybrids in field trials [3] and greenhouse trials [15], where the premature ripened plants started to appear at 43 DPI.

Black stem disease
Natural attacks of black stem disease occurred in the four trials. At the beginning of the assessment of disease resistance (14 days post inoculation at the stem base), the average percentage of plants with BS symptoms ranged from 0% (core-collection trial in 2009) to 74% (FUxPAZ2 trial in 2011) over the 3 years. The black stem disease progressed along plant development and the mean disease score reached 1.57 in core-collection trial in 2010, 1.50 in XRQxPSC8 RIL population experiment in 2010 and 1.93 in FUxPAZ2 RIL population experiment in 2011 respectively at the last observation (Table 2). A high variability between genotypes was observed in each experiment ( Table 2, Additional file 1).
Phenotypic variability in the core-collection Stem base necrosis and premature ripening in the corecollection The stem base necrosis due to Phoma macdonaldii infection showed a very high phenotypic variability among the sunflower lines (Fig. 2).
In the 20 lines present in the two core-collection trials, the percentage of PR plants (score = 4) at 63 DPI was on average 68% in 2009 and 55% in 2010, with a range from 0 to 100% each year. However, due to large developmental differences among the 20 lines, PR must be considered at the same developmental stage for each genotype: the M1.2-M1.3 stages allow a good evaluation of PR resistance, avoiding confusion of premature ripening with natural senescence. According to the genotypes, these stages are reached in 2010 between 35 and 63 DPI (Fig. 3).
In 2010, the PR AUDPC ranged from 0.23% PR plants.day to 17.53% PR plants.day (mean = 5.88% PR plants.day; Fig. 4 A correlation analysis (Kendall rank concordance test) between core-collection lines datas 2009 and 2010 on the mean genotype values of the premature ripening disease score at each phenological stage (in order to compare the genotypes in the same conditions) was performed. Of 13 of the 14 lines for which we observed at least 10 plants in 2009, the Tau b coefficient of Kendall is greater than 0.75 and is significant at the 1% threshold. Three of the four parental lines of the RIL populations developed for QTL analysis belong to this set of lines: FU, PAZ2 and XRQ (Table 3).

Black stem disease in the core-collection
Previous studies have shown the relevance of observing natural attacks for experiments related to the study of BS resistance variability of commercial hybrids [5] or to the evaluation of yied losses [26], with very consistent results between 2 years of experiments. In our experiment of 2010, at 14 DPI, the percentage of plants without any symptom on the stem under natural attack ranged from 14% for the line SF310 to 100% for SF061 (mean = 74%). The evolution of the percentage of plants with coalescent spots (Fig. 5) showed high differences in behavior between lines: SF310 and SF336 exhibited a very fast evolution of this symptom, with more than 80% of plants affected at 42 DPI. In contrast, lines Tub, 97B7 and SF334 had less than 30% of their plants with coalescent spots at 63 DPI. The SF061 line seemed to be very resistant to BS disease with 93% of healthy plants at the end of the experiment.
At each disease assessment, the analysis of variance of the mean BS disease score showed a highly significant genotype effect (p-value from 2.77*10 −9 to 0.00211), showing that the disease scoring was efficient to reveal genetic variability. This confirms the previous results obtained in field trials on six commercial hybrids [5] and on 54 other inbred lines [4].

Disease observations on RILs Stem base necrosis and premature ripening on RILs
In the two RIL populations, the PR disease score increased in the same manner from 14 dpi to 56 or 63 dpi ( Table 2). The mean of AUDPC for PR confirmed the difference between the parental lines (XRQ and PAZ2 more resistant than PSC8 and FU). In the two RIL populations, the developmental stage M1.2-M1.3 was reached between 42 and 49 DPI. At this developmental stage, the global mean disease score of XRQxPSC8 and FUxPAZ2 RIL populations was respectively 2.52 (CV = 32%) and 2.27 (CV = 29%) (Fig. 6). Within the two RIL populations, this mean disease score varied from 1.60 to 4.00 for XRQxPSC8 and 1.00 to 4.00 for FUxPAZ2. These data confirmed the behavior of the four parental lines: PSC8 and FU The mean disease score analysis of variance showed significant to highly significant differences between inbred lines in the XRQxPSC8 population (p-value from 0.0111 to 6.42*10 −13 ) and highly statistically significant differences between inbred lines in the FUxPAZ2 population (p-value from 1.54*10 −11 to P < 2*10 −16 ) according to the disease assessment date. In the two RIL populations, highly significant differences were also highlighted for AUDPC_PR (P < 0.0001).
More details are available in Additional file 1.
In the two RIL populations, broad-sense heritability ( Table 4) appeared moderate to large: respectively from 0.36 to 0.76 in the XRQxPSC8 population, from 0.73 to 0.86 in the FUxPAZ2 population. Although they were probably overestimated because the measurements were made in only one environment, these values show a certain reliability in the measured character which is a very important result considering the time required for the expression of the symptoms on adult plants.

Black stem disease on RILs
At the latest date of PR assessment, the global mean disease score of XRQxPSC8 and FUxPAZ2 RIL populations reached respectively 1.50 (CV = 26%) and 1.93 (CV = 9%). It increased 1.2 point between the first disease assessment (14 DPI) and the last notation for both RIL populations ( Table 2). The parental lines seemed to be different at the same time on their mean disease score at 14 DPI and on the evolution of the severity of the disease (AUPDC). At each assessment date, significant to highly significant differences between inbred lines have been highlighted in the two RIL populations except at 21 DPI in FUxPAZ2 population (Additional file 1).
In the FUxPAZ2 population, broad-sense heritability for black stem disease in our field trial ranged from 0.04 to 0.88 according to the disease assessment date and reached 0.88 for AUDPC_BS (Table 4). When the black stem disease is well established (from 42 DPI, mean value of BS disease score = 1.02), the estimated heritability values range between 0.75 and 0.88 and are close to that previously obtained (0.94) in F2-F3 families of the same population under controlled conditions on plantlets [8]. In the XRQxPSC8 population, h 2 ranged from 0.29 to 0.65 and reached 0.70 for AUDPC_BS. While these estimates are probably over-valued because we  In the "INEDI" population, different QTL were mapped depending on the type of disease (PR or BS). For the PR related traits, two different QTL were identified for early symptoms (Table 5, Fig. 7): a first QTL (Pm_PR_dpi14) was mapped on LG16 for early symptoms (14 DPI) accounting for 28.8% of the phenotypic variability, and another one was mapped on LG10 (Pm_PR_dpi14 & 21) accounting for 19.1 to 21% of the phenotypic variability. For the latest symptoms (42 to 63 dpi), one QTL (Pm_PR_dpi42 to 63, AUDPC_PR) is Each of these last QTL accounted between 15.6% and 19.9% of the phenotypic variability. The line XRQ was found carrying the resistance allele for PR related traits, and generally the susceptibility allele for the BS related traits.
In the trial of 2010, the "INEDI" RIL population was affected by Alternaria leaf spot and blight, with very highly statistically significant variation between RIL (Kruskal-Wallis test: P < 0.0001). A QTL related to susceptibility to Alternaria and explaining 52.8% of the phenotypic variation was found strictly colocalizing with the QTL of PR late symptoms traits (Pm_PR_dpi42, 49, 56 & 63) on LG 10. The parental line XRQ was also carrying the resistance allele to Alternaria and for the four QTL associated to PR resistance (Pm_PR_dpi42, 49, 56 & 63; AUDPC_PR) on this LG. One might suspect an artefact due to a mis-scoring of PR resistance; however, the plants showing Alternaria spots were discarded from the analysis of PR resistance, and the symptoms of the stem base necrosis due to Phoma macdonaldii and of the leaf spots due to Alternaria are clearly distinguishable. Therefore, we cannot exclude a co-location of the QTL involved in both diseases or even the involvement of a common genetic factor in the resistance to both fungi. Three QTL associated to resistance to Alternaria blight have been recently identified under natural infection in Argentina, and their effect has been confirmed in two genetic backgrounds. One of these QTL was located on LG 10 [27], consolidating our result.
In the FUxPAZ2 population (  28 & 35). These last QTL detected on LG13 during the early stage of disease development was found located in the same area as the well-known cluster of Resistance Gene Candidate [28] where race specific resistances to the Downy Mildew caused by Plasmopara halstedii have been mapped, and more precisely close to Pl21 [19]. In contrast, the QTL detected during the late stages of disease development (Pm_PR_dpi 42, 49 & 56) were found located on LG7.
Four QTL involved in Phoma resistance to stem attack on 19-day-old seedlings on FUxPAZ2 F2-F3 families have been previously detected [8]. These QTL were mapped on linkage groups 3, 8, 11 and 12 (16, 1, 14 and 15 in the denomination they used according to [29], respectively). Several QTL have been mapped for resistance to stem and root infection in a RIL population derived from the cross RHA266*PAC2, after contaminating the plantlets with respectively four [6] or three [9] different P. macdonaldii isolates. Both isolate-specific and isolate-non-specific QTL were identified. The QTL involved in the isolate-non-specific resistance to stem base inoculation was mapped on LG6 [6] and LG5 and LG15 [9] while other QTL with poor specificity were located on LG5, LG13 and LG15 [6] and LG5 [9]. These four linkage groups have been found carrying QTL of resistance on adult plants in our study. In a F2-F3population from a cross between two other sunflower lines, one isolate-non-specific QTL for resistance to petiole inoculation of plantlets was also mapped on LG5 [10]. Unfortunately, there are not enough common markers between these five maps to allow a more precise comparison. In our study, no major QTL was mapped which is consistent with the quantitative nature of the resistance. No common QTL was observed between the two RIL populations. However, there is a common trend: a) in each population, QTL involved in BS and in PR diseases are not the same, b) not the same QTL are involved in early and later stage of disease development. The existence of specific QTL associated to Sclerotinia sclerotiorum sunflower resistance of different organs at the adult developmental stage has also been demonstrated [8,30]. The QTL mapping work has been done from 1 year data. However, in further studies on the H. annuus * P. macdonaldii interaction [31], eight hybrids made in crossing the susceptible parental line FU and RILs carrying either resistant or susceptible alleles at the mapped QTL (LG7 and LG13) involved in PR were evaluated in 2013 and 2014 under natural infection for the development of phoma necrosis at the stem base. These hybrids, as well as their parental lines, appeared to be well classified in terms of resistance or susceptibility. This result represents at least indirect proof that some confidence could be given to our QTL.

Relationship between morphological, developmental traits and resistance traits
In different crops, it has been found that environmental or crop management variation on one side (for example [32]), or genotypic variation on the other side (for example [33]) might result in a modification of morphological and developmental traits which affects the crop susceptibility to diseases. As such traits were recorded on the FU*PAZ2 RIL population in the frame of a supplementary experiment (2011, see Additional file 4), we compared the genetic architectures of morphological and developmental traitsincluding seed yield -, and disease (PR, BS) related traits. A QTL accounting for 26.7% of the phenotypic variation for the grain yield (in open pollination, per se evaluation) was found co-localized with the QTL of PR resistance on LG7, with the FU parental line carrying the allele of PR susceptibility and the positive allele for yield, collar diameter, Fig. 7 Mapping of QTL detected for all the traits in the "INEDI" RIL population. Premature Ripening Traits (Pm_PR_dpiN for Phoma macdonaldii PR score N days post inoculation, AUDPC_PR), Black Stem Disease Traits (Pm_BS_dpiN for P.macdonaldii BS score N days post inoculation), Alternaria for susceptibility score to Alternaria helianthi number of leaves and M3 developmental stage. In the genetic background covered by the FU*PAZ2 population, a significant part of the genetic variation for susceptibility to PR was therefore found positively associated with vigor related traits. Although a few explanations could be proposed today for this result (see Additional file 4), it would be highly premature to put forward any hypothesis without further research. However, this result might open insights into the functional relationship between disease resistance and phenotype related traits.

Conclusions
This work shows for the first time the quantitative behaviour of the sunflower resistance to Premature Ripening caused by P. macdonaldii. It also brings evidences that different genetic factors are implicated in the disease development, depending on the infection process leading to black stem or premature ripening diseases, and on the stage of disease development. Although phenotyping for PR resistance appears hard because it must be done on adult plants in field