- Research article
- Open Access
Genetic mapping and identification of QTL for earliness in the globe artichoke/cultivated cardoon complex
BMC Research Notes volume 5, Article number: 252 (2012)
The Asteraceae species Cynara cardunculus (2n = 2x = 34) includes the two fully cross-compatible domesticated taxa globe artichoke (var. scolymus L.) and cultivated cardoon (var. altilis DC). As both are out-pollinators and suffer from marked inbreeding depression, linkage analysis has focussed on the use of a two way pseudo-test cross approach.
A set of 172 microsatellite (SSR) loci derived from expressed sequence tag DNA sequence were integrated into the reference C. cardunculus genetic maps, based on segregation among the F1 progeny of a cross between a globe artichoke and a cultivated cardoon. The resulting maps each detected 17 major linkage groups, corresponding to the species’ haploid chromosome number. A consensus map based on 66 co-dominant shared loci (64 SSRs and two SNPs) assembled 694 loci, with a mean inter-marker spacing of 2.5 cM. When the maps were used to elucidate the pattern of inheritance of head production earliness, a key commercial trait, seven regions were shown to harbour relevant quantitative trait loci (QTL). Together, these QTL accounted for up to 74% of the overall phenotypic variance.
The newly developed consensus as well as the parental genetic maps can accelerate the process of tagging and eventually isolating the genes underlying earliness in both the domesticated C. cardunculus forms. The largest single effect mapped to the same linkage group in each parental maps, and explained about one half of the phenotypic variance, thus representing a good candidate for marker assisted selection.
The Asteraceae (ex Compositae) species Cynara cardunculus L. comprises three taxa, namely the two domesticated form globe artichoke (var. scolymus) and cultivated cardoon (var. altilis), along with their common ancestor the wild cardoon (var. sylvestris). While the globe artichoke was selected for its large immature inflorescences, the cardoon was selected for its fleshy leaves and stalks. The three taxa remain fully cross-compatible with one another, and their F1 hybrids are fertile. The species complex has a highly heterozygous diploid genome (2n = 2x = 34), maintained by its cross-pollinating habit . The domesticated forms produce a variety of nutraceuticals and pharmaceutically active compounds like inulin, mono- and di-caffeoylquinic acids [2–6] and sesquiterpene lactones, which are responsible for its characteristic bitterness [7–9]. Globe artichoke contributes significantly to the Mediterranean agricultural economy in the form of an annual production of ~750Mt worth over US$500 M annually. It is also cultivated in the Americas, North Africa and China (http://faostat.fao.org).
Most of the Mediterranean globe artichoke germplasm is vegetatively propagated, and a number of varietal groups have been defined on the basis of the appearance of the inflorescence and harvesting time of the head (capitula) Flowering can be induced between autumn and spring in early flowering types by watering dormant underground shoots, whereas late flowering types flower only during spring and early summer. A common breeding target for both vegetatively and seed-propagated varieties is the promotion of earliness since inflorescences produced in the early part of the year command a higher price than those produced in the summer. Unlike globe artichoke, the cultivated cardoon is exclusively seed-propagated, and is generally handled as an annual crop. Of late it has been promoted as a source of lignocellulosic biomass [10–12] and the evidence suggests that it should be possible to derive types able to flower early, to produce stems with a high lignin content and to generate biomass with a good level of energy efficiency [13, 14]. Earliness is therefore an important trait in both domesticated forms.
The first generation of C. cardunculus marker-based genetic maps [15–17] have resulted in a cultivated cardoon map composed of nearly 200 loci (17 major linkage groups, LGs) spanning just over 10 M, and a globe artichoke one featuring 326 loci (20 major LGs) spanning about 15 M. The two maps have since been integrated on the basis of common loci with the inclusion of a number of genes involved in the synthesis of caffeoylquinic acids [18, 19]. More recently crosses between globe artichoke and its ancestor wild cardoon have generated highly segregating F1 populations exploitable as ornamentals  as well as for mapping studies .
The multi-allelism of many microsatellite (SSR) loci makes them particularly well suited as bridging markers to link independent maps. The design of SSR assays requires DNA sequence, which in globe artichoke exists at present largely in the form of expressed sequence tag (EST) sequence (http://compgenomics.ucdavis.edu/). Over 4,000 potential EST-SSR loci have been identified from this sequence resource, and the experimental testing of a sample of 300 loci showed that more than one half were informative between the parents of our two mapping populations . In the present report, we describe the integration into the globe artichoke and cultivated cardoon maps of a large number of these EST-SSR loci, and show that they can be used as bridging markers to merge the two maps. The resulting dense maps was then used to identify a number of quantitative trait loci (QTL) underlying early head production in C. cardunculus.
Results and discussion
Six of the 178 informative Cynara Expressed Microsatellite (CyEM) markers, identified by Scaglione et al. , were excluded from the analysis on the basis of excessive missing values. Of the remaining 172, 54 segregated in both parents (46 as 1:1:1:1, eight as 1:2:1) and 118 in just one of the parents (85 in globe artichoke ‘Romanesco C3’, 33 in cultivated cardoon ‘Altilis 41’). On the whole 228 microsatellite markers were available for map construction (Table 1). Eleven of these loci suffered from mild segregation distortion (χα=0.052 < χ2 ≤ χα=0.012) but just one (CyEM_58) from severe distortion (χ2 > χα=0.012). Since CyEM_58 was excluded from the mapping analysis, this left a total of 227 SSR loci. Co-dominant markers appear to be less affected by segregation distortion than dominant ones [23, 24], and this certainly was the case for C. cardunculus, where ~13% of AFLP and S-SAP loci , but only ~5% of SSRs and SNPs are distorted. Segregation distortion has been associated with statistical bias and/or with errors in genotyping, but they can also stem from a number of biological phenomena affecting meiosis, fertilization and embryogenesis  as well as the presence of null alleles. Null alleles at SSR loci are not uncommon, as they can arise where either one (or both) of the primers fail to anneal because of sequence mismatch or the deletion of the whole locus, and cause an higher apparent number of homozygotes because they can no longer be distinguished from the heterozygotes . In this situation, the options are either to disregard the affected loci, to score segregation in the same way as for a dominant marker , to attempt to redesign the primers [28, 29], or to adjust allele frequencies on the basis of a global estimate of the frequency of null alleles. As recently described by Lanteri et al.  the null alleles segregating in a Mendelian fashion were identified, thus limiting the segregation distortion in our populations. As noted previously [15, 17], although the inclusion of loci distorted at the 1% level and above increases the frequency of type I errors, it does help to maintain marker density throughout the map.
The updated ‘Romanesco C3’ map was built from 574 loci (359 AFLPs, 19 S-SAPs, 189 SSRs and seven SNPs), and the ‘Altilis 41’ one from 373 loci (246 AFLPs, 8 S-SAPs, 114 SSRs and five SNPs); of these, 78 (76 SSRs and two SNPs) were in common between the two parental genotypes. The CyEM SSRs were less informative in the cultivated cardoon than in the globe artichoke. Of the 228 assayed SSR loci 189 (83%) segregated in ‘Romanesco C3’, but only 114 (50%) in ‘Altilis 41’ (Table 1). The difference in level of heterozygosity between these parents has been remarked on before  and is thought to be a consequence of the sustained vegetative propagation used in globe artichoke, in contrast to the seed propagation applied to the cultivated cardoon, which led to a certain degree of purifying selection aimed at stabilizing production.
Globe artichoke map
The globe artichoke ‘Romanesco C3’ map (LOD 6.0) consisted of 473 loci falling into 20 LGs, each containing at least eight loci (Figure 1). The number of mapped SSRs has now risen from 46 to 185. The largest LG contained 73 loci, and the range in genetic length of the individual LGs was 34.5-140.9 cM. CyEM loci (139 markers) were mapped to all the major LGs, and their inclusion allowed the integration of six AFLP loci which previously had remained unlinked . Two LGs (C3_13 and C3_18) which were previously separated have now been merged, while LG C3_4 has been split into C3_4a and C3_4b as a consequence of more stringent LOD applied (Figure 1). LG C3_17 has increased in genetic length by 36 cM (86%), while that of LG C3_3 and C3_8 has increased by ~30% and ~20%, respectively. The map spanned 1543.8 cM, with a mean inter-marker distance of 3.40 cM, corresponding to a 3.8% increase in length over the earlier map , but in a ~28% decrease in the mean inter-marker distance. The proportion of intervals shorter than 7 cM is now 88% (previously 77%), and only six gaps of >15 cM remain. The SSRs appeared to be rather quite uniformly dispersed, although some clustering is present in the distal regions of C3_8, C3_2 and C3_17, and around the putative centromeric region, of C3_3, C3_15 and C3_20. These chromosomal regions are typically enriched for SSRs [30–35]. The relatively low marker saturation present in the distal regions of C3_3 and C3_14 presumably reflects a localized reduced level of polymorphism between the mapping parents.
Some segregation distortion was present at five of the CyEM loci (CyEM_19, _47, _70, _73 and _231; three at α = 0.05 and two at α = 0.01, Figure 1) which affected a cluster of loci on both C3_17 and C3_9. In both cases the distortion was due to an excess of the band detected in the female parent, thus it is likely to have a biological basis, rather than being due to either scoring error or chance . Biological mechanisms causing segregation distortions have been extensively studied in Drosophila , and are known to occur in many plant species [38–42]. On the other hand, the other 18 distorted loci were scattered across the genome, a common feature in the genetic maps of both plant and animal species .
By lowering the LOD threshold from 6.0 to 5.0, three pairs of LGs were merged: C3_10 with C3_5, C3_14 with C3_8 and C3_4a with C3_4b, resulting in the formation of 17 LGs (corresponding to the haploid chromosome number, Figure 1). It also allowed the inclusion of two unlinked pairs of loci (one into C3_2 and the other into C3_12) and the singlet AFLP locus e35/m46-156 (into C3_7). This generated an increase in the genetic length of the map of ~60 cM; one doublet still remains unlinked (Figure 1).
Both the goodness-of-fit of marker placement (mean χ2 contribution) and nearest neighbour fit (cM) were evaluated. Compared to the earlier ‘Romanesco C3’ map , the average mean χ2 contribution of markers across the LGs has been significantly reduced from 5.38 to 4.42 (t test at α = 0.005), highlighting the improvement in robustness. The variation in this parameter for each LG is illustrated in Figure 2, which confirms that LGs C3_1, C3_2, C3_5, C3_8, C3_10, C3_17 and C3_20 have all shown an improved goodness-of-fit. C3_12 remained largely composed of AFLP loci (only two CyEM loci were integrated) and thus its robustness was hardly improved. The mean nearest neighbour fit of the CyEM loci (24.3 ± 3.7 cM) was markedly lower (t test at α = 0.005) than that of the AFLP loci (51.0 ± 4.6 cM), confirming the desirability of including co-dominant markers to obtain reliable marker placement.
Cultivated cardoon map
The genetic map of the cultivated cardoon ‘Altilis 41’ parent (LOD 6.0) was constructed from 373 segregating loci (82 CyEM loci), of which 273 were ordered into 21 major LGs, whose length ranged from 27.1 to 125.2 cM, with the largest LG consisting of 29 loci. The result of integrating the CyEM loci was an increase in the number of major LGs from 17 to 21. This involved the recognition of four new LGs (Alt_18 to _21), the splitting of Alt_1 into two (Alt_1a and Alt_1b) and the merging of Alt_16 and Alt_1b (Figure 3). The updated ‘Altilis 41’ map included 107 SSR loci distributed across all but one (Alt_13) of the major LGs, with a total genetic length of 1485.7 cM and a mean inter-marker distance of 5.44 cM. This represents a marked increase in both length (+42%) and number of loci (+50%), together with a minor decrease in the mean inter-marker distance (−5%). The proportion of intervals smaller than 10 cM (about 80%) was not significantly reduced. Some of the LGs recorded large increases in their genetic length – for example, that of Alt_17 by 64.7 cM, Alt_14 by 58.9 cM and Alt_18 by 57.8 cM. The only LG which recorded a reduction in length was Alt_5. There was some clustering of CyEM loci in the distal region of Alt_2 and around the putative centromeric region of Alt_1b and Alt_15. Three CyEM loci (CyEM_73, _3 and _231; Figure 3) showed a degree of segregation distortion (two at α = 0.05 and one at α = 0.01), but none of these were linked to other distorted loci, similar to the other nine markers showing segregation distortion. The addition of the new SSR markers decreased the mean inter-marker distance on Alt_11 by ~60%, and some gaps in the previous map have been filled; but ten gaps of >15 cM remained, perhaps reflecting regions of genetic fixation which have arisen during cultivated cardoon domestication.
Lowering the LOD threshold to 5.0 led to the merging of four pairs of LG: Alt_1a with Alt_1b, Alt_11 with Alt_2, Alt_7 with Alt_10, and Alt_19 with Alt_13. At this level of stringency the number of LGs corresponded to the haploid chromosome (Figure 3). The lowered stringency also allowed the incorporation of two groups of three linked loci into LGs Alt_20, and Alt_5, and of one doublet into LG Alt_6. As a result, the overall length of the map was increased by 133.5 cM; one triplet and four doublets still remain unlinked (Figure 3).
The C. Cardunculus consensus map
The number of informative shared co-dominant markers was raised from 19 to 66 (64 SSRs, two SNPs), representing from one to 15 bridging markers per LG. As a result, 19 of the ‘Romanesco C3’ LGs were alignable with 20 of the ‘Altilis 41’ ones (Table 2). There was a one to one correspondence between 18 LG pairs, but C3_1 shared markers with both Alt_11 and Alt_2 (Table 2). C3_4b remained non-aligned, but did harbour a number of SSR loci which were informative for the second step of the analysis; this was not the case for Alt_13 (Table 2). The alignment was followed by the construction of a consensus map based on a LOD threshold of above 5.0 (Figure 4), which succeeded in capturing 694 loci, 227 (217 SSRs, ten SNPs) of which involved co-dominant markers. The map generated 17 LGs with a total genetic length of 1687.6 cM and a mean inter-marker spacing of 2.5 cM; consensus LG numbers (from LG I to LG XVII) have been assigned (Table 2, Figure 4). The length of each individual LG varied from 44.5 to 144.5 cM (mean 99.3 cM), with the largest containing 92 loci. Only three of the CyEM loci (the intercross locus CyEM_134 and the two ‘Altilis 41’ testcross loci CyEM_167 and _79) remained unlinked.
The consensus map, obtained from the domesticated C. cardunculus forms, was compared with the Sonnante et al.  map constructed from a cross between the var. scolymus cultivar ‘Mola’ and the var. sylvestris (wild cardoon) accession ‘Tolfa’, by considering 125 (117 SSRs, eight SNPs) common markers. The common markers identified each of the 17 LGs on the consensus map, with between two and 21 present on each LG (Table 2). Ten of the LGs aligned readily; LGs V and VII aligned with two ‘Mola’/‘Tolfa’ LGs, and LG VIII with two major groups and a triplet of markers. LGs X/XVII, and XIII/XVI each aligned with only a single ‘Mola’/‘Tolfa’ LG. In general, marker order and genetic separation were comparable, with some exceptions. It has been established that wild cardoon is more divergent from the two cultivated forms (globe artichoke and cultivated cardoon) than are the two cultivated forms with respect to one another [44, 45]. Somewhat surprisingly, therefore, over 100 SSR loci featured in the consensus map but apparently were either non-informative or remained as singlet loci in the ‘Mola’/‘Tolfa’ population.
EST-SSRs as functional markers
Putative functions can be deduced for markers derived from ESTs using homology searches with public protein databases. Annotation of mapped loci was performed via BlastX search as well as InterPro scan and GO categorisation made it possible to tag some biological functions.
A set of 17 CyEM markers were annotated with GO terms involved in the ‘response to stimulus’ (Table 3), five of which were derived from transcripts related to ‘response to cold stress’ and eight to ‘response to salt stress’ terms. In particular, the marker CyEM-42, developed from the contig CL4773Contig1 (1281 bp, 267 aminoacids)  and mapped on LG_12 of “Romanesco C3” map, showed high amino acid similarity (81%) with the Arabidopsis protein kinase PBS1 (NP_196820.1, unigene At.23518). To consider reliable orthology, a reciprocal tblastx analysis against the whole EST collection, currently available for C. cardunculus, was performed and no better alignment than that of contig CL4773 was detected. PBS1 was found to work as R gene against the bacterial pathogen Pseudomonas syringae, where its cleavage, operated by the pathogens’ effector AvrPphB, triggers the signalling cascade, generating the host response (HR) . Pseudomonas spp. together with other endophytic bacteria may affect globe artichoke plants both in field and during micropropagation  and the CyEM-42 may be likely considered a reliable marker for tagging a bacterial resistance trait in the species.
Our EST-SSR markers may be defined as functional markers with the potential to target polymorphisms in gene responsible for traits of interest and they can be also particularly useful for constructing comparative framework maps with other Asteraceae, giving the possibility to amplify ortholog genes and provide anchor loci.
The genetic basis of earliness
An evaluation of the variance for the three earliness-related traits established significant genotypic differences (P < 0.05) between ‘Romanesco C3’ and ‘Altilis 41’ (Table 4). Thus, eMH in the former was 162 days in “2009” and 178 days in “2008”, while in the latter the respective times were 218 and 223 days. All three traits varied continuously among the F1 progeny (the distribution for eMH is shown in Figure 5); no progeny was as early flowering as ‘Romanesco C3’, but a few were later flowering than ‘Altilis 41’, due to transgressive segregation. The mean eMH, eFOH and eSOH lay substantially above the mid-parent value, suggesting semi-dominance for lateness. The low global level of heterozygosity characteristic of the cultivated cardoon makes it possible that one or more of the earliness QTL are in the homozygous state in ‘Altilis 41’, so that the presence of dominant alleles for lateness may contribute to later flowering across the whole mapping population. The inter-trait correlations were similar in both seasons, with the strongest correlation linking eMH and eFOH (r > 0.80, P < 0.0001). The correlations between the two seasons were also strong, ranging from 0.64 (P < 0.0001) for eMH to 0.49 (p < 0.001) for eSOH (Table 5). Flowering and head harvesting time was a little earlier in “2009” than in “2008” (7–8 days on average), while performance was somewhat more variable in “2009” (Table 4), probably reflecting the difference between re-awakened and newly sown material. The broad sense heritability for eMH of 0.76 (Table 4) indicated the trait to be predominantly under genetic control, but the rather lower heritabilities shown by the traits eFOH and eSOH suggested that the environment is quite influential in their determination.
The KW test and SIM procedure identified, at first, six QTL regions stable across years in the developed consensus map (Figure 4). Those on LGs I, XI and XVII involved all three traits, those on LGs I and IX only eMH and eFOH, and the one on LG VII solely eMH. The seventh QTL cluster on LG II involved all three traits, but was only expressed in “2009” (Figure 4). On the whole, seven chromosomal regions scattered over six LGs of the consensus map were identified. When the ‘Romanesco C3’ and ‘Altilis 41’ maps were used separately for QTL validation, the percentage of phenotypic variance explained by some of the QTL differed from that predicted by the analysis based on the consensus map (data not shown), perhaps reflecting the structure and size of the segregating progeny and the existence of different allelic interactions . However, all seven QTL regions were detectable by applying the SIM method to the parental maps, and further analysed with the MQM procedure. QTL identified in each map and season are shown in Table 6 and graphically reported in Figure 6. Only three of the seven QTL regions were detectable in both parental maps, presumably these regions were heterozygous in both parental lines. The other four were only detectable in one of the two maps, suggesting that one parent was homozygous in the critical region (Figure 6). Across all three traits, a total of 25 QTL was detected, of which 19 were stable across both growing seasons, with the other six expressed only in “2009”.
With respect to eMH, two of the QTL were heterozygous in both parents, three only in ‘Romanesco C3’ and two just in ‘Altilis 41’. The largest effect stable eMH QTL in ‘Romanesco C3’ mapped to LG C3_1 in the neighbourhood of the SSR locus CELMS_40, named eMH.C3_1. This QTL was responsible for 38-48% of the phenotypic variation and was associated with an additive effect of 10–12 days. The other four QTL in ‘Romanesco C3’ mapped to LGs C3_9, _8, _12 and _2, and accounted individually for between 6-10% of the phenotypic variance; eMH.C3_2 was only detected in “2009”. The largest stable QTL detected in ‘Altilis 41’ (eMH.Alt_2), homologous to the ones detected in the same region of the ‘C3’ map, explained from 33-41% of the phenotypic variance and was associated with an additive effect of 9–11 days. A second QTL, eMH.Alt_4, was detectable only in “2009”, but its location suggested it to be identical with eMH.C3_2 (Figure 6). Further two minor QTL present on LGs Alt_2 and _15 accounted for, respectively, 8% and-11% of the variance. Globally, the QTL identified in the ‘Romanesco C3’ and ‘Altilis 41’ maps accounted for, respectively, 74% and 62% of the phenotypic variance for main head harvesting time in “2008”.
Six eFOH QTL were detected, three of which were represented in both parental maps, one on just the ‘Romanesco C3’ map, and the other two on just the ‘Altilis 41’ map. Five of the six QTL mapped to a region where a eMH QTL was also located, with overlapping LOD confidence intervals but with an overall lower phenotypic effect. The exception was eFOH.Alt_6 (Table 6, Figure 6). As for eMH, the largest effect QTL mapped to LGs C3_1 and Alt_2 in the neighbourhood of CyEM_223. Based on the “2009” data, the set of eFOH QTL accounted for 47% (‘Romanesco C3’) and 54% (‘Altilis 41’) of the variation.
Only four eSOH QTL were uncovered, due to the reduced heritability of this trait (hB2 = 0.54, Table 4). All four co-localized with eFOH QTLs, with an overall lower phenotypic effect (Table 6, Figure 6), with the largest effect QTL mapping to the cluster on C3_1 and Alt_2. Based on the “2009” data, the set of eSOH QTL accounted for 48% (‘Romanesco C3’) and 43% (‘Altilis 41’) of the variation.
We have reported here an extension of the C. cardunculus genetic map by introducing SSR loci sited within genic sequence. The integration of 139 of these loci has significantly improved the resolution and accuracy of the maps. Given that the genome size of the species is ~1.08Gbp , the mean equivalence between the physical and genetic length in this species is of the order of 1 cM to 670 kbp. Thus the mean physical separation of the mapped markers is around 2.2Mbp. On this basis, most gene sequences should lie within about 1Mbp of the nearest marker, although this value makes the non-valid assumption that recombination sites are randomly distributed along the length of the chromosomes.
Shortening the life cycle is seen as an important breeding goal in terms of both globe artichoke’s economic value, and the ease of exploiting cultivated cardoon as an energy crop [13, 14]. The newly developed consensus as well as the parental genetic maps can accelerate the process of tagging and eventually isolating the genes underlying earliness in both the domesticated C. cardunculus forms. We have shown that a cluster of large effect QTL resides on the homologous LGs C3_1 and Alt_2, and this clearly represents a reasonable candidate for marker assisted breeding. The critical genetic interval contains two SSR loci (CELMS_40 and CyEM_223), either of which is well-placed to serve as an indirect selection criterion for earliness. Before such a genotypic selection programme can be implemented, however, a validation of the presence and importance of the QTL needs to be conducted using different genetic backgrounds and in other relevant environments [50, 51]. To date, this study represents the first attempt to identify QTL in C. cardunculus, which is the necessary preliminary step for implementing marker-assisted selection for quantitative traits. Beyond tagging, mapping also prepares the ground for positional cloning, which will enable the molecular basis of trait variation to be identified.
Plant material and SSR analysis
The 178 informative CyEM markers identified by Scaglione et al.  were used to genotype a set of 94 F1 hybrid (randomly selected from 154 true hybrids as described by Portis et al.  from the cross’Romanesco C3’(globe artichoke; female parent) x ‘Altilis 41’ (cultivated cardoon; male parent). A 7 ng aliquot of genomic DNA from each mapping population individual was amplified in a 10 μl reaction containing 1x PCR buffer, 1 mM MgCl2, 0.5U Taq DNA polymerase (Qiagen Inc., Venlo, Netherlands), 40nM 5’-labelled (FAM, HEX or TAMRA) forward primer, 40 nM unlabelled reverse primer and 0.2 mM dNTP. A touchdown cycling regime was applied, consisting of an initial denaturation of 94°C/2.5 min, followed by nine cycles of 94°C/30s, 63°C/30s (decreasing by 0.7°C per cycle), 72°C/60s, and 30 further cycles of 94°C/30s, 57°C/30s, 72°C/60s. Where only weak amplification was achieved, the MgCl2 concentration was raised to 1.5 mM and the final annealing temperature was lowered to 55°C. The amplicons were separated on an ABI3730 capillary DNA sequencer (Applied Biosystem Inc., Foster City, CA, USA). Internal ROX-labelled GS500 size standards were included in each capillary. The output was analysed by GeneMapper v3.5 software (Applied Biosystems).
Linkage analysis and parental maps construction
The CyEM genotypes of the 94 mapping population individuals were combined with previous genotypic data based on 605 AFLP, 27 S-SAP and 56 other SSRs , along with ten SNP from genes underlying caffeoylquinic acid synthesis (reported by Comino et al.  and Menin et al. ). JoinMap v4.0  was used to generate two separate linkage maps (one for each parent) using the double pseudo-testcross mapping strategy . The markers fell into three classes: maternal testcross markers segregating only in ‘Romanesco C3’ (expected segregation ratio 1:1); paternal testcross markers segregating only in ‘Altilis 41’ (1:1); and intercross markers segregating within both parents (either 1:2:1 or 1:1:1:1). Differences between observed and expected segregation ratios were tested by χ2, and only markers deviating if at all only slightly from expectation (χα=0.12 < χ2 ≤ χα=0.012) were used for map construction and the estimation of genetic distances, when their presence did not alter surrounding marker order in the LG. Heavily distorted loci (χ2 > χα=0.012), along with those associated with 30 or more missing values, were excluded. LGs were established on the basis of an initial LOD threshold of 6.0. Locus order and distances between loci were established using the following parameter set: Rec = 0.40, LOD = 1.0, Jump = 5. Map distances were converted to cM using the Kosambi mapping function . Where a locus order discrepancy arose between a pair of parental LGs, the marker order of the ‘1:1:1:1’ segregating SSR and the marker order of SNP markers were taken as the ‘fixed order’. Once the framework maps had been established, additional loci were subsequently added and some LGs merged by lowering the LOD threshold to 5.0. On the resulting maps, loci suffering from slight segregation distortion have been identified with either one (χα=0.12 < χ2 ≤ χα=0.052) or two (χα=0.052 < χ2 ≤ χα=0.012) asterisks. The ‘Romanesco C3’ LGs are labelled LG_C3, and the ‘Altilis 41’ ones LG_Alt, using the numbering system suggested by Portis et al. .
Consensus map construction
A genotypic data set based on all the available markers was then used to construct a consensus map. Here, the loci belonging to two segregation classes ‘1:2:1’ (the same pair of alleles segregating in each parent), and ‘1:1:1:1’ (different alleles segregating in each parent) were used as ‘bridge markers’. The most likely locus order was established from a comparison of the ‘C3’, ‘Alt’ and consensus LGs, and where these differed substantially from one another, the most likely order was assumed to be one associated with the lowest χ2 value (estimating goodness-of-fit) and the lowest mean χ2 contribution for all loci. LGs were established on the basis of an initial LOD threshold of 5.0 and numbered according to C3 maps LGs order.
Annotation of mapped sequences was carried out using the Blast2Go software . The best twenty BlastX results were retrieved by querying the nr protein database at NCBI. Gene Ontology terms were retrieved accordingly to software capabilities and transferred to our sequences, adopting an annotation threshold of 55. Additional GO terms were obtained performing InterPro scan for conserved motif, using all available databases. ANNEX function was used to obtain further GO terms which are implicit on the base of electronically annotated ones.
Earliness evaluation and QTL analysis
The mapping population (154 F1 progeny of the cross ‘Romanesco C3’ (var. scolymus) x ‘Altilis 41’ (var. altilis)), along with six clones of each parental line, was cultivated at the University of Catania's experimental station (37°25’N; 15°30’E; 10 m a.s.l) and evaluated over the two growing seasons 2007–2008 (hereafter referred to as “2008”) and 2008–2009 (“2009”). In “2008”, seedlings at the three true leaf stages (about 40 days after germination) were transferred to the field in mid September, while in “2009”, the growing season was initiated by applying drip irrigation to field capacity in mid August. Earliness was scored either as the number of days between transplanting (“2008”) or awakening (“2009”) and harvesting of both the main (eMH trait) and first and second order heads (obtained from the ramification of the main stem: eFOH and eSOH traits). Population means, standard deviations, distribution histograms and trait correlations were calculated using R software . Analyses of variance were based on treating each growing season as an independent replicate. The broad sense heritability was given by the expression hB2 = σg2/(σg2 + σe2/y), where σg2 represented the genetic variance and σe2 the error variance. Correlations between traits were estimated using Pearson’s coefficient.
In the initial step of the QTL analysis, the consensus map was used to assign putative locations by performing a Kruskal-Wallis (KW) non-parametric test in conjunction with a simple interval mapping procedure (SIM) , applying the cross-pollination algorithm implemented within MapQTL v4.0 software . Next, the two separate parental maps were employed for a re-analysis based on the BC1 algorithm, using both SIM and multiple QTL mapping (MQM) . Markers lying within a putative QTL region and associated with the highest LOD score were used as co-factors. For the MQM, a backward elimination procedure was applied to select the appropriate co-factors (significantly associated with each trait at P < 0.02). The LOD thresholds for QTL significance were confirmed by a permutation test consisting of 1,000 replications, which implies a genome-wide significance level of 0.05 . Only those QTL associated with a LOD greater than either the genome-wide threshold or the LG threshold were considered. 1-LOD support intervals were determined for each LOD peak . The additive effect and the proportion of the overall phenotypic variance associated with each QTL and all QTL together were estimated from the MQM model. Linkage maps and QTL position were drawn using MapChart .
Basnitzki J, Zohary D: Breeding of seed planted artichoke. Plant Breeding Reviews. 1994, 12: 253-269.
Schutz K, Kammerer D, Carle R, Schieber A: Identification and quantification of caffeoylquinic acids and flavonolds from artichoke (Cynara scolymus L.) heads, juice, and pomace by HPLC-DAD-ESI/MSn. Journal of Agricultural and Food Chemistry. 2004, 52 (13): 4090-4096. 10.1021/jf049625x.
Lombardo S, Pandino G, Mauromicale G, Knödler M, Carle R, Schieber M: Influence of genotype, harvest time and plant part on polyphenolic composition of globe artichoke [Cynara cardunculus var. scolymus (L.) Fiori]. Food Chemistry. 2010, 119: 1175-1181. 10.1016/j.foodchem.2009.08.033.
Pandino G, Courts F, Lombardo S, Mauromicale G, Williamson G: Caffeoylquinic acids and flavonoids in the immature inflorescence of globe artichoke, wild cardoon, and cultivated cardoon. Journal of Agricultural and Food Chemistry. 2010, 58: 1026-1031. 10.1021/jf903311j.
Pandino G, Lombardo S, Mauromicale G, Williamson G: Phenolic acids and flavonoids in leaf and floral stem of cultivated and wild Cynara cardunculus L. genotypes. Food Chemistry. 2011, 126: 417-422. 10.1016/j.foodchem.2010.11.001.
Lattanzio V, Kroon PA, Linsalata V, Cardinali A: Globe artichoke: a functional food and source of nutraceutical ingredients. Journal of Functional Foods. 2009, 1 (2): 131-144. 10.1016/j.jff.2009.01.002.
Menin B, Comino C, Portis E, Moglia A, Cankar K, Bouwmeester H, Lanteri S, Beekwilder J: Genetic mapping characterization of the globe artichoke (+)-germacrene A synthase gene, encoding the first dedicated enzyme for biosynthesis of the bitter sesquiterpene lactone cynaropicrin. Plant Science. 2012, 190: 1-8.
Cravotto G, Nano G, Binello A, Spagliardi P, Seu G: Chemical and biological modification of cynaropicrin and grosheimin: a structure-bitterness relationship study. Journal of the Science of Food and Agriculture. 2005, 85: 1757-1764. 10.1002/jsfa.2180.
Fritsche J, Beindorff C, Dachtler M, Zhang H, Lammers J: Isolation, characterization and determination of minor artichoke (Cynara scolymus L.) leaf extract compounds. European Food Research and Technology. 2002, 215: 149-157. 10.1007/s00217-002-0507-0.
Foti S, Mauromicale G, Raccuia S, Fallico B, Fanella F, Maccarone E: Possible alternative utilization of Cynara spp. I. Biomass, grain yield and chemical composition of grain. Industrial Crops and Products. 1999, 10: 219-228. 10.1016/S0926-6690(99)00026-6.
Ierna A, Mauromicale G: Cynara cardunculus L. genotypes as a crop for energy purposes in a Mediterranean environment. Biomass and Bioenergy. 2010, 34: 754-760. 10.1016/j.biombioe.2010.01.018.
Encinar J, Gonzalez J, Rodriguez J, Tejedor A: Biodiesel fuels from vegetable oils: Transesterification of Cynara cardunculus L. oils with ethanol. Energy and Fuels. 2002, 16 (2): 443-450. 10.1021/ef010174h.
Acquadro A, Portis E, Longo A, Mauro R, Mauromicale G, Lanteri S: Potentiality of Cynara cardunculus L. as energy crop. Journal of Biotechnology. 2010, 150: 165-166.
Acquadro A, Portis E, Scaglione D, Mauro RP, Campion B, Falavigna A, Zaccardelli R, Ronga D, Perrone D, Mauromicale G, et al: Cynara cardunculus L.: From Vegetable to Energy Crop. Proceedings of the Joint Meeting AGI-SIBV-SIGA. 2011, Società Italiana di Genetica Agraria (ed),
Lanteri S, Acquadro A, Comino C, Mauro R, Mauromicale G, Portis E: A first linkage map of globe artichoke (Cynara cardunculus var. scolymus L.) based on AFLP, S-SAP, M-AFLP and microsatellite markers. Theoretical and Applied Genetics. 2006, 112 (8): 1532-1542. 10.1007/s00122-006-0256-8.
Acquadro A, Lanteri S, Scaglione D, Arens P, Vosman B, Portis E: Genetic mapping and annotation of genomic microsatellites isolated from globe artichoke. Theoretical and Applied Genetics. 2009, 118: 1573-1587. 10.1007/s00122-009-1005-6.
Portis E, Mauromicale G, Mauro R, Acquadro A, Scaglione D, Lanteri S: Construction of a reference molecular linkage map of globe artichoke (Cynara cardunculus var. scolymus). Theoretical and Applied Genetics. 2009, 120: 59-70. 10.1007/s00122-009-1159-2.
Comino C, Hehn A, Moglia A, Menin B, Bourgaud F, Lanteri S, Portis E: The isolation and mapping of a novel hydroxycinnamoyltransferase in the globe artichoke chlorogenic acid pathway. BMC Plant Biology. 2009, 9: 30-10.1186/1471-2229-9-30.
Menin B, Comino C, Moglia A, Dolzhenko Y, Portis E, Lanteri S: Identification and mapping of genes related to caffeoylquinic acid synthesis in Cynara cardunculus L. Plant Science. 2010, 179 (4): 338-347. 10.1016/j.plantsci.2010.06.010.
Lanteri S, Portis E, Acquadro A, Mauro RP, Mauromicale G: Morphology and SSR fingerprinting of newly developed Cynara cardunculus genotypes exploitable as ornamentals. Euphytica. 2012, 184: 311-321. 10.1007/s10681-011-0509-8.
Sonnante G, Gatto A, Morgese A, Montemurro F, Sarli G, Blanco E, Pignone D: Genetic map of artichoke x wild cardoon: toward a consensus map for Cynara cardunculus. Theoretical and Applied Genetics. 2011, 123 (7): 1215-1229. 10.1007/s00122-011-1661-1.
Scaglione D, Acquadro A, Portis E, Taylor C, Lanteri S, Knapp S: Ontology and diversity of transcript-associated microsatellites mined from a globe artichoke EST database. BMC Genomics. 2009, 10: 454-10.1186/1471-2164-10-454.
Lorieux M, Goffinet B, Perrier X, Deleon D, Lanaud C: Maximum-likelihood models for mapping genetic markers showing segregation distortion. 1. Backcross populations. Theoretical and Applied Genetics. 1995, 90: 73-80.
Lorieux M, Perrier X, Goffinet B, Lanaud C, Deleon D: Maximum-likelihood models for mapping genetic markers showing segregation distortion. 2. F2 populations. Theoretical and Applied Genetics. 1995, 90: 81-89.
Bradshaw H, Villar M, Watson B, Otto K, Stewart S, Stettler R: Molecular genetics of growth and development in Populus. 3. A genetic-linkage map of a hybrid poplar composed of RFLP, STS and RAPD markers. Theoretical and Applied Genetics. 1994, 89 (2–3): 167-178.
Pekkinen M, Varvio S, Kulju K, Karkkainen H, Smolander S, Vihera-Aarnio A, Koski V, Sillanpaa M: Linkage map of birch, Betula pendula Roth, based on microsatellites and amplified fragment length polymorphisms. Genome. 2005, 48 (4): 619-625. 10.1139/g05-031.
Rodzen J, May B: Inheritance of microsatellite loci in the white sturgeon (Acipenser transmontanus). Genome. 2002, 45: 1064-1076. 10.1139/g02-083.
Shaw P, Turan C, Wright J, O'Connell M, Carvalho G: Microsatellite DNA analysis of population structure in Atlantic herring (Clupea harengus), with direct comparison to allozyme and mtDNA RFLP analyses. Heredity. 1999, 83: 490-499. 10.1038/sj.hdy.6885860.
Van Oosterhout C, Hutchinson W, Wills D, Shipley P: MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Molecular Ecology Notes. 2004, 4: 535-538. 10.1111/j.1471-8286.2004.00684.x.
Arens P, Odinot P, Vanheusden A, Lindhout P, Vosman B: GATA repeats and GACA repeats are not evenly distributed throughout the tomato genome. Genome. 1995, 38 (1): 84-90. 10.1139/g95-010.
Bhattramakki D, Dong J, Chhabra A, Hart G: An integrated SSR and RFLP linkage map of Sorghum bicolor (L.) Moench. Genome. 2000, 43 (6): 988-1002.
Gill G, Wilcox P, Whittaker D, Winz R, Bickerstaff P, Echt C, Kent J, Humphreys M, Elborough K, Gardner R: A framework linkage map of perennial ryegrass based on SSR markers. Genome. 2006, 49 (4): 354-364. 10.1139/G05-120.
Jones E, Dupal M, Dumsday J, Hughes L, Forster J: An SSR-based genetic linkage map for perennial ryegrass (Lolium perenne L.). Theoretical and Applied Genetics. 2002, 105 (4): 577-584. 10.1007/s00122-002-0907-3.
McCouch S, Teytelman L, Xu Y, Lobos K, Clare K, Walton M, Fu B, Maghirang R, Li Z, Xing Y, et al: Development and mapping of 2240 new SSR markers for rice (Oryza sativa L.). DNA Research. 2002, 9 (6): 199-207. 10.1093/dnares/9.6.199.
Ramsay L, Macaulay M, Ivanissevich S, MacLean K, Cardle L, Fuller J, Edwards K, Tuvesson S, Morgante M, Massari A, et al: A simple sequence repeat-based linkage map of barley. Genetics. 2000, 156 (4): 1997-2005.
Fishman L, Kelly AJ, Morgan E, Willis JH: A genetic map in the Mimulus guttatus species complex reveals transmission ratio distortion due to heterospecific interactions. Genetics. 2001, 159 (4): 1701-1716.
Lyttle T: The Genetics and Evolutionary Biology of Meiotic Drive - Preface. American Naturalist. 1991, 137 (3): 281-282. 10.1086/285162.
Xu Y, Zhu L, Xiao J, Huang N, McCouch S: Chromosomal regions associated with segregation distortion of molecular markers in F-2, backcross, doubled haploid, and recombinant inbred populations in rice (Oryza sativa L). Molecular and General Genetics. 1997, 253: 535-545. 10.1007/s004380050355.
Faris J, Laddomada B, Gill B: Molecular mapping of segregation distortion loci in Aegilops tauschii. Genetics. 1998, 149: 319-327.
Lu H, Romero-Severson J, Bernardo R: Chromosomal regions associated with segregation distortion in maize. Theoretical and Applied Genetics. 2002, 105 (4): 622-628. 10.1007/s00122-002-0970-9.
Kleinhofs A, Kilian A, Maroof M, Biyashev R, Hayes P, Chen F, Lapitan N, Fenwick A, Blake T, Kanazin V, et al: A molecular, isozyme and morphological map of the barley (Hordeum vulgare) genome. Theoretical and Applied Genetics. 1993, 86 (6): 705-712.
Ky C, Barre P, Lorieux M, Trouslot P, Akaffou S, Louarn J, Charrier A, Hamon S, Noirot M: Interspecific genetic linkage map, segregation distortion and genetic conversion in coffee (Coffea sp.). Theoretical and Applied Genetics. 2000, 101: 669-676. 10.1007/s001220051529.
Taylor D, Ingvarsson P: Common features of segregation distortion in plants and animals. Genetica. 2003, 117: 27-35. 10.1023/A:1022308414864.
Lanteri S, Saba E, Cadinu M, Mallica G, Baghino L, Portis E: Amplified fragment length polymorphism for genetic diversity assessment in globe artichoke. Theoretical and Applied Genetics. 2004, 108 (8): 1534-1544. 10.1007/s00122-003-1576-6.
Portis E, Barchi L, Acquadro A, Macua J, Lanteri S: Genetic diversity assessment in cultivated cardoon by AFLP (amplified fragment length polymorphism) and microsatellite markers. Plant Breeding. 2005, 124 (3): 299-304. 10.1111/j.1439-0523.2005.01098.x.
Shao F, Merritt P, Bao Z, Innes R, Dixon J: A Yersinia effector and a Pseudomonas avirulence protein define a family of cysteine proteases functioning in bacterial pathogenesis. Cell. 2002, 109 (5): 575-588. 10.1016/S0092-8674(02)00766-3.
Penalver R, Duranvila N, Lopez MM: Characterization and pathogenicity of bacteria from shoot tips of the globe artichoke (Cynara-scolymus L). Annals of Applied Biology. 1994, 125 (3): 501-513. 10.1111/j.1744-7348.1994.tb04987.x.
Cabezas JA, Cervera MT, Ruiz-Garcia L, Carreno J, Martinez-Zapater JM: A genetic analysis of seed and berry weight in grapevine. Genome. 2006, 49 (12): 1572-1585. 10.1139/g06-122.
Marie D, Brown S: A cytometric exercise in plant DNA histograms, with 2 C-values for 70 species. Biology of the Cell. 1993, 78 (1–2): 41-51.
Melchinger AE, Utz HF, Schon GC: QTL analyses of complex traits with cross validation, bootstrapping and other biometric methods. Euphytica. 2004, 137 (1): 1-11.
Cloutier S, Ragupathy R, Niu ZX, Duguid S: SSR-based linkage map of flax (Linum usitatissimum L.) and mapping of QTLs underlying fatty acid composition traits. Molecular Breeding. 2011, 28 (4): 437-451. 10.1007/s11032-010-9494-1.
Van Ooijen JW: JoinMap® v.4, Software for the calculation of genetic linkage maps in experimental populations. 2006, Kyazma BV, Wageningen, Netherlands
Weeden N: Approaches to mapping in horticultural crops. Gressho. PM (ed) Plant genome analysis. 1994, CRC Press Boca Raton, , 57-68.
Kosambi D: The estimation of map distances from recombination values. In.Ann Eugen. 1944, 12: 172-175.
Conesa A, Gotz S, Garcia-Gomez J, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.
R: a language and environment for statistical computing Version 2.11.1. 2008, Vienna, Austria: R Foundation for Statistical Computing (ed), ,http://www.R-project.org,
Lander ES, Botstein D: Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121 (1): 185-199.
Van Ooijen JW, Boer MP, Jansen RC, Maliepaard C: MapQTL 4.0: software for the calculation of QTL positions on genetic maps. Plant Research International. 2002, , Wageningen, The Netherlands
Jansen RC, Stam P: High-resolution of quantitative traits into multiple loci via interval mapping. Genetics. 1994, 136 (4): 1447-1455.
Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138 (3): 963-971.
Van Ooijen JW: Accuracy of mapping quantitative trait loci in autogamous species. Theoretical and Applied Genetics. 1992, 84 (7–8): 803-811.
Voorrips R: MapChart: Software for the graphical presentation of linkage maps and QTLs. Journal of Heredity. 2002, 93 (1): 77-78. 10.1093/jhered/93.1.77.
Acquadro A, Portis E, Lanteri S: Isolation of microsatellite loci in artichoke (Cynara cardunculus L. var. scolymus). Molecular Ecology Notes. 2003, 3 (1): 37-39.
Acquadro A, Portis E, Lee D, Donini P, Lanteri S: Development and characterization of microsatellite markers in Cynara cardunculus L. Genome. 2005, 48 (2): 217-225. 10.1139/g04-111.
Acquadro A, Portis E, Albertini E, Lanteri S: M-AFLP-based protocol for microsatellite loci isolation in Cynara cardunculus L. (Asteraceae). Molecular Ecology Notes. 2005, 5 (2): 272-274. 10.1111/j.1471-8286.2005.00897.x.
Sonnante G, Carluccio A, De Paolis A, Pignone D: Identification of artichoke SSR markers: molecular variation and patterns of diversity in genetically cohesive taxa and wild allies. Genetic Resources and Crop Evolution. 2008, 55 (7): 1029-1046. 10.1007/s10722-008-9310-5.
This research was supported by grants from: (i) the National Science Foundation Plant Genome Research Program (No. 0421630), (ii) the Georgia Research Alliance, (iii) the University of Georgia Research Foundation, and (iv) by MIPAAF (Ministero delle Politiche Agricole, Alimentari e Forestali - Italy) through the CYNERGIA (“Costituzione e valutazione dell’adattabilita’ di genotipi di Cynara cardunculus per la produzione di biomassa e biodiesel in ambiente mediterraneo”) project .
The authors declare that they have no competing interests.
SL and SK planned and supervised the experimental work; GM and RM cultivated the C. cardunculus parental accessions and the F1 individuals and performed earliness trait evaluation over the two growing seasons; DS and AA performed the genotyping of the progenies; EP and DS carried out linkage analyses and map construction; EP performed QTL analysis; SL and EP supervised the drafting of the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Portis, E., Scaglione, D., Acquadro, A. et al. Genetic mapping and identification of QTL for earliness in the globe artichoke/cultivated cardoon complex. BMC Res Notes 5, 252 (2012). https://doi.org/10.1186/1756-0500-5-252
- Cynara cardunculus
- Linkage map