- Research article
- Open access
- Published:
Construction of a potato consensus map and QTL meta-analysis offer new insights into the genetic architecture of late blight resistance and plant maturity traits
BMC Plant Biology volume 11, Article number: 16 (2011)
Abstract
Background
Integrating QTL results from independent experiments performed on related species helps to survey the genetic diversity of loci/alleles underlying complex traits, and to highlight potential targets for breeding or QTL cloning. Potato (Solanum tuberosum L.) late blight resistance has been thoroughly studied, generating mapping data for many Rpi-genes (R-genes to Phytophthora infestans) and QTLs (quantitative trait loci). Moreover, late blight resistance was often associated with plant maturity. To get insight into the genomic organization of late blight resistance loci as compared to maturity QTLs, a QTL meta-analysis was performed for both traits.
Results
Nineteen QTL publications for late blight resistance were considered, seven of them reported maturity QTLs. Twenty-one QTL maps and eight reference maps were compiled to construct a 2,141-marker consensus map on which QTLs were projected and clustered into meta-QTLs. The whole-genome QTL meta-analysis reduced by six-fold late blight resistance QTLs (by clustering 144 QTLs into 24 meta-QTLs), by ca. five-fold maturity QTLs (by clustering 42 QTLs into eight meta-QTLs), and by ca. two-fold QTL confidence interval mean. Late blight resistance meta-QTLs were observed on every chromosome and maturity meta-QTLs on only six chromosomes.
Conclusions
Meta-analysis helped to refine the genomic regions of interest frequently described, and provided the closest flanking markers. Meta-QTLs of late blight resistance and maturity juxtaposed along chromosomes IV, V and VIII, and overlapped on chromosomes VI and XI. The distribution of late blight resistance meta-QTLs is significantly independent from those of Rpi-genes, resistance gene analogs and defence-related loci. The anchorage of meta-QTLs to the potato genome sequence, recently publicly released, will especially improve the candidate gene selection to determine the genes underlying meta-QTLs. All mapping data are available from the Sol Genomics Network (SGN) database.
Background
The number of publications reporting the mapping of QTLs (quantitative trait locus) in plants has exponentially increased since the Eighties, reaching a total of about 34,300 papers in 2010 (source: Google Scholar with key words "QTL" and "plant"). For a few species only, this huge amount of QTL data has been recorded in databases that enable quick comparison of QTL mapping results from independent experiments (e.g. Gramene for maize and rice). But for most species, QTL data accumulates in bibliography until the coming out of hot-spot genomic regions that become targets for introgression into breeding material or for cloning. To get a comprehensive understanding of the genetic control of a polygenic trait and to optimize its use in breeding, it is needed to get a complete view of the genetic architecture of the trait with the distribution of the involved loci along the genome. This synthesis can be greatly facilitated by achieving a QTL meta-analysis.
The general principle of a meta-analysis is to pool the results of several studies that address the same issue to improve the estimate of targeted parameters. Meta-analysis was first used in social and medical sciences, like epidemiology. More recently, it was applied in plant genetics to combine on a single map the genetic marker data and the QTL characteristics (location, confidence interval, effect and trait used for QTL detection) from independent QTL mapping experiments to finally estimate the optimal set of distinct consensus QTLs, called meta-QTLs. The positions of those meta-QTLs are estimated with a higher accuracy as compared to the individual QTLs in the original experiments [1]. To date, QTL meta-analyses have been achieved for traits related to plant development and plant response to environment (nutrients, abiotic and biotic stresses) in maize, wheat, rice, rapeseed, cotton, soybean, cocoa and apricot [2–18].
Statistical methods have been proposed for the meta-analysis of QTLs from several experiments. The method proposed by Goffinet and Gerber (2000) was implemented in the Biomercator software [1, 19]. It compiles the QTLs that have been projected on an existing reference map and uses the transformed Akaike classification criterion to determine the best model between one QTL, two QTLs, three QTLs etc. until the maximum number of QTLs mapped in the same region. This method was first used by Chardon et al. (2004) and by most authors until recently [2, 3, 6, 8–10, 15, 16]. Then Veyrieras et al. (2007) have extended the statistical method and implemented the new algorithms in the MetaQTL software [20]. MetaQTL notably uses a weighted least squares strategy to build the consensus map from the maps of individual studies and offers a new clustering approach based on a Gaussian mixture model to define the optimal number of QTL clusters or meta-QTLs on each chromosome that best explain the observed distribution of the individual projected QTLs. The Gaussian mixture model has shown to be flexible and robust to the non-independence of the experiments [4]. Moreover, simulations demonstrated that the number of meta-QTLs selected by the Akaike criterion is lower than the expected number with random distributions of QTLs and that it has a very low probability to happen by chance [4]. The MetaQTL software has successfully been used in wheat, maize, rice and apricot [4, 5, 12, 13, 17].
Potato (Solanum tuberosum L.) late blight resistance is typically a trait for which meta-analysis can be applied. From 1994 to 2009, 19 studies have been published on QTL mapping in different crosses and with different related species, generating a significant amount of QTL data. All these publications reflect the interest of the potato scientific community towards polygenic partial resistance to late blight. Late blight, caused by the oomycete Phytophthora infestans, is one of the most serious diseases in potato, which is the third most important food crop in the world after rice and wheat. Almost all Rpi-genes (R-genes to P. infestans) deployed in the potato fields have been rapidly overcome, while polygenic resistance appears to be a fairly efficient and durable alternative. However, it has been observed that this kind of resistance in potato is often associated with plant maturity, as most resistant plants are also the ones that mature the latest. This is a handicap for breeders and growers who aim to get early maturing plants to shorten the time of tuber production.
Attempts to get a synthetic view of the loci controlling polygenic late blight resistance in potato with comparison of their positions with maturity QTLs have already been published [21, 22]. However, because of a lack of common markers, the comparison of QTLs was achieved at a half-chromosome scale, which made the compilation imprecise. Consequently, to enhance the comparison of QTL positions coming from different mapping studies and also to refine the localization of hot-spot genomic regions, the mapping of common markers between maps is crucial.
Reference dense maps constructed with transferable markers are privileged sources of common markers. A UHD potato map containing 10,000 AFLP markers has been designed to become a reference map [23, 24]. However, the anchorage of AFLP markers is restricted to closely-related species. In addition, as the comparison is based on the comigration of the marker bands on the gel, AFLP gels are required, which does not make the comparison easy to achieve [25]. Other reference maps containing SSR and RFLP markers have been developed in potato (SSR maps [26–28]; RFLP map [29]). These markers are well defined by specific primers or a probe sequence, which makes them easily transferable from one cross to another, even between distantly related species; they are thus handy tools for map comparison.
A functional map for pathogen resistance, enriched with RGA (resistance gene analog) and DRL (defence-related locus) sequences, SNPs and InDels tightly linked or located within NBS-LRR-like genes, has been developed on the basis of two potato populations (BC9162 and F1840 [30–33]; PoMaMo database [34]). This functional map also contains CAPS, SSR and RFLP literature-derived markers, which enables the comparison with other QTL maps. However, it remains difficult to infer precisely functional locus information to QTL mapping results as QTLs often have large confidence intervals.
QTL meta-analysis thus appears here to be an adequate tool i) to narrow-down the confidence intervals of hot-spot loci where congruent late blight resistance QTLs of multiple origins map, and ii) to investigate colocalization of these loci with Rpi-genes, RGAs, DRLs and maturity QTLs as well. In this paper, we present a three-step meta-analysis process achieved with the MetaQTL software. First, we built a consensus potato map by compiling 21 QTL maps and eight reference maps. This consensus map includes common markers and specific markers tagging Rpi-genes, as well as RGA and DRL markers. Second, individual QTLs for late blight resistance and maturity were projected onto the consensus map. Third, for each trait, QTLs were clustered into meta-QTLs on the basis of the distribution of their projected positions on the consensus potato map.
Results
Bibliographic review of QTL mapping studies
The initial map set comprised a total of 37 maps divided into i) 29 QTL maps from 19 publications related to QTL detection of late blight resistance and maturity type, and ii) eight independent reference maps (without any QTL) (Table 1). Reference maps were included because they provided numerous pivotal markers, which improved connections between maps. Because of a lack of shared markers, the initial 29 QTL map set was refined to a core subset of 21 "connected" QTL maps coming from 14 publications that were included in the meta-analysis (Table 1).
The 21 "connected" QTL maps were representative of the diversity of assessments for late blight resistance and maturity, the QTL detection methods and the sources of resistance (Table 2). Resistance tests were based on disease spread on foliage in the field (FF) or in the greenhouse (FG), sporulation or necrosis spots on in vitro detached leaflets or leaf discs (LT), necrosis progression on stems (ST) and disease damage on tuber slices (TS) or whole tubers (T% or WT) in controlled conditions. Maturity type was evaluated by the number of days before flowering or senescence (MT), plant height (PH) and plant vigour (PV). QTLs were detected with different statistical detection methods according to the number of available markers, the size of the progeny and the frequency distribution profile of the raw or transformed data (non-parametric statistical tests or ANOVA, Interval Mapping, Composite Interval Mapping or Multiple QTL Mapping with permutation tests). Most of the P. infestans isolates used for late blight resistance assessments were of A1 mating type and virulent towards the 11 S. demissum Rpi-genes. However, it was difficult to say whether some of the isolates used in the different studies were the same or not. As wild tuber-bearing relatives of potato have proven to be high-potential sources of resistance, most mapping populations derived from a cross between a dihaploid S. tuberosum clone (the susceptible parent) and a clone derived from a diploid relative (the resistant parent). Two mapping populations even derived from crosses between two potato relatives (without S. tuberosum, Table 2). The parental pedigrees were sometimes quite complex. Nevertheless, the marker order in all maps was well conserved and aligned with the S. tuberosum map [35, 36]. If all known species of the parent pedigrees are taken into account, a total of 13 potato-related species were involved in the meta-analysis.
Consensus potato map
Common markers between the 21 "connected" QTL maps and eight reference maps (Table 3) made it possible the construction of a consensus map for the 12 potato chromosomes. The number of maps used to construct each consensus chromosome varied between 20 and 25 (Figure 1). The consensus potato map had a total length of 1,260 cM (Haldane) and contained a total of 2,141 markers (SSR, SSCP, CAPS, RFLP, AFLP, SNP, InDels and STS markers). Among them, 514 markers were shared by at least two different maps. There were between 28 and 58 common markers per chromosome, corresponding to 16% up to 29% of the total number of markers per chromosome. The name, map position and occurrence of each marker are given in Additional file 1 and on the SGN database [37].
QTL dataset for meta-analysis
On the basis of the 19 publications of QTL studies, a total of 211 late blight resistance QTLs and 64 maturity QTLs were collected (Table 1). However, some QTL intervals did not include the minimum of two anchor markers, which were required for their projection onto the consensus map. Thus, the QTL dataset for meta-analysis was reduced down to 144 late blight resistance QTLs and 42 maturity QTLs, coming from 14 publications. The excluded QTLs, which harboured a single common marker with the consensus map, were referred to "anchored QTLs" and indicated at this marker position in Additional file 1 but their orientation and projected confidence interval could not be determined.
As far as the QTLs included in the meta-analysis are considered, late blight and maturity QTLs spread on the 12 potato chromosomes. The number of QTLs per chromosome ranged between six and 21 for late blight resistance, and between one and eight for maturity.
For late blight resistance, R² values were available for 106 QTLs out of the 144 input QTLs and ranged between 4% (chromosome I, foliage test [38]; chromosomes V, IX, XI, XII, foliage test [39]) and 63% (chromosome X, tuber test [40]). 75% of the late blight QTLs had relatively small effects, ranging between 4% and 15%; 7% of the QTLs had large effects, ranging between 30% and 63%. Confidence intervals ranged between 3 cM (chromosome III, leaf disc test [41]) and 66 cM (chromosome VI, foliage test [42]), with a mean of 24 cM.
For maturity, R² values were available for 20 QTLs out of the 42 input QTLs and ranged between 4% (chromosomes IX and XII [39]) and 71% (chromosome V [42]). 75% of the maturity QTLs had R² values ranging between 4% and 15%; 10% of the QTLs explained more than 30% of the phenotypic variation (60% and 71% on chromosome V [42]). Confidence intervals ranged between 4 cM (chromosome XI [42]) and 61 cM (chromosome VI [42]), with a mean of 20 cM.
Meta-analysis
We determined the number of meta-QTLs per chromosome by using the modified Akaike Information Criterion (AICc) and by taking into account the consistency with the different criteria provided by the MetaQTL software (Additional file 2). Our analysis yielded a total of 32 meta-QTLs. Each meta-QTL corresponded to clusters of individual QTLs coming from different experiments. Meta-QTLs were composed of a maximum of 18 individual QTLs for late blight resistance (chromosome V) and eight individual QTLs for maturity (chromosome V). The QTL meta-analysis on the whole potato genome reduced by six-fold the initial number of late blight QTLs by passing from 144 QTLs to 24 meta-QTLs and by ca. five-fold the number of maturity QTLs by passing from 42 QTLs to eight meta-QTLs. Figure 2 presents an example of the meta-analysis steps for chromosome IV, from QTL projection on the consensus map to QTL clustering into meta-QTLs.
A graphical overview of the late blight and maturity meta-QTLs is presented on Figure 3. Late blight meta-QTLs spread on the 12 chromosomes, with one to three meta-QTLs per chromosome. Maturity meta-QTLs spread on only six chromosomes, with one or two meta-QTLs per chromosome. Other maturity QTLs were reported in literature on the other six chromosomes, but they were single in their region and no meta-QTL could be computed. Single QTLs for late blight resistance and for maturity that were excluded from the clustering step are shown in Additional file 1, with the other excluded QTLs which were anchored by a single marker to the consensus map.
The confidence intervals of late blight meta-QTLs ranged between 0.27 cM (chromosome VII) to 49.81 cM (chromosome I), with a mean of 10.25 cM (SD±10.79). The confidence intervals of maturity meta-QTLs ranged between 0.88 cM (chromosome V) to 39.28 cM (chromosome VI), with a mean of 10.67 cM (SD±12.54). With respect to the length reduction of the mean confidence interval from the published QTLs to the meta-QTLs, confidence intervals were reduced by 2.3-fold for late blight resistance and by 1.9-fold for maturity (Additional file 3).
Maturity meta-QTLs overlapped late blight meta-QTLs on chromosomes VI and XI, while there was no overlap on chromosomes IV, V, VII and VIII. However, by running meta-analysis on late blight resistance QTLs and maturity QTLs altogether under a single "super-trait", we observed that for all 12 chromosomes, maturity QTLs were always clustered together with late blight resistance QTLs in a "super meta-QTL" (data not shown). On the other way round, we observed at least one "super meta-QTL" free of maturity QTLs for 11 chromosomes; for chromosome XI only, both "super meta-QTLs" included at least one maturity QTL.
The three most consistent late blight meta-QTLs were located on chromosomes IV, V and X (MQTL_1_Late_blight of chromosome IV, MQTL_1_Late_blight of chromosome V and MQTL_2_Late_blight of chromosome X; Additional file 3). These meta-QTLs were composed of the highest number of QTLs (10 to 18 QTLs) with the largest effects (R² up to 63%, tuber test [40]). In addition, they corresponded to individual QTLs identified in different potato-related species or in plant material with complex pedigree. This result suggests that these regions could correspond to conserved resistance QTLs retrieved from several tuber-bearing Solanum species.
Candidate gene analysis
Literature reported the map positions of several Rpi-genes determining late blight resistance (reviewed in [43, 44]). However, only a few flanking markers were supplied (Rpi-genes were linked to a single marker or included in a large marker interval), which hampered the accurate location of Rpi-genes on the consensus map (Additional file 1). Due to their rough positions, it was thus not possible to say definitely whether they were included or not in the late blight meta-QTLs. Out of the 33 Rpi-genes positioned on our consensus potato map, 10 were included in the confidence intervals of late blight meta-QTLs (Table 4). One example of overlap was on chromosome IV, where the TG370-TG339 marker interval (~12 cM) containing a large NBS-LRR Rpi-gene cluster (R2-like genes) largely overlapped the meta-QTL MQTL_1_Late_blight [45]. On chromosome VI, the CT119 marker tagging the Rpi-blb2 R-gene was included in MQTL_1_Late_blight. On chromosome X, the TG422 and TG403 markers flanking the Rpi-ber2 gene were included in MQTL_2_Late_blight. However, on chromosome XI, the lack of anchor markers hindered the accurate location of the 10 Rpi-genes (Rpi-Stirling, R5 to R11, R3a and R3b). According to the flanking markers (STM5130-STM5109 for Rpi-Stirling, TG105-GP250 for R3a, TG26 for R3b and R5 to R11), only Rpi-Stirling might be included in MQTL_2_Late_blight.
Conversely, a few Rpi-genes clearly did not belong to any late blight meta-QTLs. This was the case for Rpi1 on chromosome VII and for the Rpi-vnt1, Rpi-phu1 and Rpi-mcq1/moc1 loci of chromosome IX. In three additional cases, the distinction between Rpi-genes and late blight meta-QTLs was doubtful. On chromosome V, R1 gene (BA213c14 and BA87d17 BACs) was located less than 2 cM far below the lower bound of MQTL_1_Late_blight. On chromosome VIII, the RB cluster (Rpi-blb1, Rpi-pta1, Rpi-plt1, Rpi-sto1, tagged by RB marker) was located 1 cM far up to the upper bound of MQTL_2_Late_blight [46]. On chromosome X, the Rber/Rpi-ber1 locus was located between both meta-QTLs of this chromosome, in a 3-cM interval (Additional file 1).
In total, 80 RGAs and 72 DRLs were reported on our consensus map, mainly from the PoMaMo functional map [32, 34]. Fourteen RGAs and 26 DRLs belonged to late blight meta-QTLs that covered about 20% of the consensus map (Table 4). Comparatively, 24 RGAs and nine DRLs belonged to maturity meta-QTLs that covered about 7% of the consensus map (Additional file 1). Independency Chi-2 tests indicate that the number of RGAs and Rpi-genes are under expectation in late blight meta-QTLs (p-value=0.035 under the hypothesis of independency) and over expectation in maturity meta-QTLs (p-value<0.0001). The heterogeneous distribution of RGAs and Rpi-genes corroborate the fact that they are often clustered or alleles. Conversely, the distribution of DRLs was independent on the distribution of both late blight meta-QTLs (p-value=0.323) and maturity meta-QTLs (p-value=0.909).
Discussion
A dense consensus reference potato map for map comparisons
Twenty-nine published potato maps were merged together into a single consensus map. From the information available in the publications of the genetic maps, at least three maps come from the cultivated potato species (S. tuberosum) and 23 maps from crosses between S. tuberosum and potato wild relatives (S. microdontum, S. phureja, S. sparsipilum, etc., Tables 2 and 3). Sixteen maps are already consensus maps of both parents, with generally one being a S. tuberosum clone and the other one a resistant wild potato species. This ability to compile genetic map information of S. tuberosum and its wild relatives indicates a high level of conservation of the marker order, and thus, of genomic sequences all over the genome. This stresses the very close genetic relationships of those genetic backgrounds and gives evidence of the validity to compile their deriving published QTL data produced with their maps. The genetic relationships between the cultivated potato and its wild relatives have been described in details by Spooner et al. (2008) [47].
Composed of 2,141 markers, the consensus map constructed in our study constitutes a new valuable dense reference map of potato. Marker positions are available on the SGN database, enabling map comparisons [37, 48]. This map can be used either as a source of markers for regions of interest or as an anchor reference map. However, for regions with a high density of markers, precise marker order has to be taken with care as the marker positions were calculated according to the positions of common markers from different maps and not based on recombination fractions. Thus, precise marker positions have to be checked by mapping in a large population.
One feature of this consensus map is that markers are not regularly spread along the chromosomes and tend to concentrate in the medium regions. For example, on chromosome VIII, 18 markers are spread on the top 71 cM while 167 markers are spread on the next 29 cM. This phenomenon is indeed observed on the published maps where dense regions are often assimilated to centromeric regions characterized by a small number of recombinations and consequently compressed maps. We also assume that genomic regions known to be involved in late blight resistance mostly gather in medium regions, which had been enriched with markers. Another explanation would be that distal markers generally originate from a single published map and their positions could be due to genotyping errors or skewed segregations.
A clear picture of the structural organization of late blight resistance loci on the potato genome
The synthetic potato map with meta-QTLs offers a refined overview of the structural organization of the loci of polygenic resistance to late blight in terms of number of QTLs and lengths of confidence intervals. By reducing the number of resistance loci by a factor of six (from 144 QTLs to 24 meta-QTLs), meta-analysis highlights the well-known resistance gene clusters on chromosomes IV and V, and also points out loci which had not appeared as notable in individual experiments like the loci on chromosomes I, III, VII, and XII. Twenty-four meta-QTLs summarized about 96% of the individual QTLs included in the analysis, which illustrates the power of QTL meta-analysis to combine QTLs from various studies.
One may question the validity of such QTL meta-analysis compiling information of as different species as the cultivated potato (S. tuberosum) and its wild relatives (S. stenotomum, S. berthaultii, S. bulbocastanum, etc.). However, as we explained earlier, this meta-analysis could only be performed thanks to the presence of common molecular markers mapped in a conserved order across maps of different related species. This structural tight genetic relationships of the different backgrounds sets up the hypothesis that the same genes are present in the same order in the genome across species, and that the genetic variation, if any, would take place at the allele level. The fact that the Potato Genome Sequencing Consortium is currently exploiting the high genomic similarity between S. tuberosum and S. phureja to reduce the complexity in assembly supports this hypothesis [49]. A high level of sequence conservation was also observed at the nucleotide level of the coding sequence among six Solanaceae genera (potato, tomato, pepper, petunia, tobacco and Nicotiana benthamiana) [50]. Because the number of sequence matches among different Solanaceae EST libraries was inversely correlated with the phylogenetic distance, we assume that the tuber-bearing species are also very similar at the level of expressed genes. These hypotheses have been already proposed in other genera where meta-analysis was conducted across relatives (e.g. for cotton fiber [15] and for rice blast [2]). In potato, resistance QTLs from one population frequently mapped, as far as resolution allows, in close proximity to those described in other populations. At the gene level, high sequence homology of Rpi genes were described between potato relatives. Functional homologues of the R2 resistance gene to P. infestans located on potato chromosome IV were cloned by an allele mining approach in three related species, and recognized the same effector protein of P. infestans [45]. Rpi genes and their general functions are overall well conserved across potato related species, variation being in fact limited to differences at the base pair and allele function levels [51]. At last, in a recent published study, meta-analysis was performed on populations involving different species related to bread wheat to narrow-down the interval of a QTL controlling the nitrogen use efficiency [52]. The functional underlying gene has been identified and showed to be conserved at orthologous positions in wheat species and in much further related cereals species such as rice, sorghum and maize. These different elements demonstrate that the analysis of the genetic factors controlling a trait across genomes of different related species and even of different genera of a plant family can be very powerful to perform a map-based dissection of a conserved gene that controls the same trait in several species.
In our study, the locus confidence intervals have been reduced by 2.3-fold in average. Locus accuracy has especially increased for the loci for which colocalizing QTLs are numerous, like on chromosome V where 18 colocalizing QTLs with an averaged confidence interval of 23 cM were clustered into a single meta-QTL of only 5 cM. In this way, meta-analysis refines the genomic regions of interest frequently described. This enables the determination of a set of markers for selection and a reasonable list of candidate genes when the genetic map is anchored to the annotated genome sequence. The closest flanking markers of the locus are also provided for subsequent fine-mapping in a real large population for map-based gene cloning.
However, confidence intervals of meta-QTLs have to be taken with caution, as they are the result of two successive statistical transformations (projection onto a consensus map and clustering), based themselves on the individual QTL confidence intervals. This stresses the importance of the accuracy of the initial mapping data. For our analysis, we took into account the individual QTL confidence intervals as described in the original publication when available (interval length of a certain LOD decrease). Otherwise, the confidence interval estimate was calculated with the empirical formula of Darvasi and Soller (1997) whose accuracy depends on the population size and the QTL effect [53]. To reduce the risk of giving artificially too much weight to a locus, we made a quite strict selection of the individual QTLs to be included in the meta-analysis by removing all possible redundancy (e.g. several repetitions or related traits in the same experiment). In our study, five meta-QTLs displayed confidence intervals lower than 1 cM (on chromosomes III, V, VII, VIII and XI), while the confidence interval of individual QTLs varied greatly. In these regions, the consensus map appeared condensed in comparison with the original maps; therefore, the projected confidence intervals of individual QTLs were very tight. Löffler et al. (2009) have also found in wheat a very tight meta-QTL of 0.1 cM that encompassed six QTLs only [12]. These over-reduced confidence intervals underline the necessity to validate the marker-trait association, either by association mapping or by transcriptomics when possible as performed by Norton et al. (2008) [14].
The meta-analysis implemented by the Meta-QTL software assumes that QTL experiments are independent from each other. Therefore, it could not take into account common features to several studies, such as the relatedness between mapping populations or between Phytophthora isolates that would have increased the power of the analysis. Nevertheless, by projecting QTLs on the same consensus map, meta-analysis still makes it possible to rapidly compare QTL mapping results of linked studies and highlights QTLs conferring isolate-specific resistance (same population and assessment but different isolates [54, 55]) or with tissue-specificity effect (same population and isolate but assessments on foliage and tuber [35, 40]). Another limit of the meta-analysis implemented in MetaQTL is that it does not provide the direction of the allelic effects, meaning that QTLs composing a meta-QTL can have opposite direction alleles. Consequently, we have to come back to the individual QTL data to be able to select the origin of the target favourable allele.
Polygenic late blight resistance and maturity relationships
Late blight meta-QTLs overlapped maturity meta-QTLs on chromosomes VI and XI. In addition, if we consider individual maturity QTLs excluded from the meta-analysis but anchored with a single common marker, other overlaps were presumed with late blight meta-QTLs on chromosome I, II, V, and XII, and reciprocally individual late blight QTLs overlapped maturity meta-QTLs on chromosomes IV, V, VI, VII, and XI (Additional file 1). In these cases, either pleiotropic genes might control both traits, or the resolution is not accurate enough to distinguish two closely linked genes.
The most famous association between QTLs for late blight resistance and for maturity is located on the upper part of chromosome V [39, 40, 42, 55–57]. Here, meta-analysis results show that the maturity meta-QTL MQTL_1_Maturity consisting in eight maturity QTLs reported for this chromosome was very close but distinct from the most consistent late blight meta-QTL MQTL_1_Late_blight consisting in 18 individual QTLs. This result goes rather in favour of the hypothesis that each trait would be controlled by independent genes but very closely linked. Nevertheless, this result has to be taken with care as two individual maturity QTLs that were not included in the meta-analysis were anchored to markers of the confidence interval of MQTL_1_Late_blight. Also, another individual late blight QTL was anchored to a marker of the MQTL_1_Maturity's confidence interval (Additional file 1).
Cases of clearly distinct maturity and late blight meta-QTLs were found on chromosomes IV, VII and VIII; the hypothesis of a pleiotropic gene would then be excluded for these regions. In addition, anchored individual maturity QTLs did not coincide with late blight meta-QTLs on chromosomes III and X, corroborating the independence of late blight and maturity loci in these regions as suggested in previous studies [39, 40, 55, 58]. The clearest cases of physical independency between maturity and late blight resistance QTLs could be preferential targets for introgression into elite cultivars for late blight resistance breeding.
Candidate genes for late blight resistance QTLs
Most frequent hypotheses about resistance QTLs were either defeated R-genes with residual effects or defence genes. The confidence intervals of the late blight meta-QTLs of chromosomes IV, VI, X and XI include R-genes of the NBS-LRR class (R2 cluster, Rpi-blb2, Rpi-ber2 and Rpi-Stirling respectively). The late blight resistance locus of chromosome IV was particularly well documented with the Stirling cultivar and S. microdontum cases of study [55, 59, 60]. The detection of the same locus as an R-gene or a QTL can be accounted for several factors such as the allelic form of the gene (for instance, defeated complete resistance alleles could be detected as QTLs), the composition of the pathogen isolate, the way of scoring the disease or the genetic background. Bhaskar et al. (2008) demonstrated that the resistance level conferred by the RB gene was dependent on the quantity of proteins encoded by the essential cell cycle regulator SGT1 gene. This showed the importance of the genetic background in the efficiency of R-gene-triggered disease resistance [61]. A parallel can be made with the fact that quantitative resistance is often controlled by a few large-effect QTLs in association with several minor-effect QTLs which can interact with the major QTLs to modulate the expression of the given trait (reviewed in [62]).
By analysing colocalizations between the 24 late blight meta-QTLs and the 33 Rpi-genes, 80 RGAs and 72 DRLs, we observed that 25% of late blight meta-QTLs included RGAs (33% included Rpi-genes or RGAs), and 50% included DRLs. It also appeared that the frequency of RGAs was not significantly greater inside late blight meta-QTL confidence intervals than outside, and that DRLs were neither significantly associated with late blight meta-QTLs, nor with maturity meta-QTLs. Even if our analysis was biased by the limited number of candidate genes and QTLs, our results do not favour one particular hypothesis for molecular basis of resistance QTLs rather than another, corroborating Ballini et al.'s conclusions [2]. However, the meta-analysis presents the advantage to reduce QTL confidence intervals, which contributes to increase the resolution in selecting relevant candidate genes. As an example, the StAOS2 gene encoding the potato allene oxide synthase 2 was located within the MQTL_2_Late_blight's confidence interval. Pajerowska-Mukhtar et al. (2008) showed that the natural variation of this gene was associated with a late blight resistance QTL identified by Oberhagemann et al. (1999) [63, 64]. Such congruency between meta-analysis and fine mapping results was also reported for the Vgt1 QTL in maize [20, 65].
Conclusions
In our study, we produced the first consensus map and performed the first meta-analysis dealing with both development trait and resistance to a biotic stress in potato. Through this study, we demonstrated that, as soon as a large amount of QTL data is collected from different studies and connected by common genetic markers, meta-analysis becomes a powerful tool to highlight chromosomal regions to focus further researches on and to use in breeding. To narrow-down the target loci confidence intervals, it is thus worth systematically integrating all new QTLs into meta-analysis on a regular basis. The anchorage of the new annotated potato genome sequence to meta-QTLs will especially provide interesting targets for candidate gene approach and for marker-assisted breeding [66]. Meta-analysis could also be useful for comparative QTL mapping across widely related crops of the same family, as achieved between rice and maize [67]. This opens a new type of analysis that would integrate gene evolution and functional conservation.
To improve meta-analysis, it would be necessary to integrate the relationships between parental clones across experiments, along with their pedigree, to be able to determine the donors of resistant alleles. In addition, adding information on P. infestans isolates used for resistance assessments would enlighten on the resistance spectrum mediated by meta-QTLs, which is one of the predictors of the durability of resistance to pathogens. We assume that broad-spectrum meta-QTLs probably target essential functions of the pathogen and that meta-QTLs supported by QTLs detected from several genitors or related species probably provide a selective advantage. Consequently, we presume that meta-QTLs with a broad-spectrum and retrieved from different genitors correspond to constrained genes, and could therefore be preferential targets to increase the durability of the resistance.
Methods
Consensus potato map
The construction of the consensus map was performed chromosome by chromosome. To be able to align the chromosome maps in the right orientation, a chromosome of a study should contain a minimum of two common markers with the corresponding chromosome of another study. QTL maps that did not share common markers were discarded from the construction of the consensus map. For several maps, few chromosomes were also missing, which lead to a variation of the number of input maps depending on the chromosome (Figure 1).
In case of inversions of two markers between maps, only the marker that was present in the lowest number of maps was manually removed to ensure that the most frequent common markers would be systematically retained. We repeated this process until no more inversion was observed between maps.
The ConsMap command of the MetaQTL software version 1.0 was used to create the consensus marker map [20]. The implemented method is based on a Weighted Least Square (WLS) strategy, which made it possible to compile all the input maps into a consensus map in a single step. It takes into account the distances between adjacent markers from all individual maps rescaled in Haldane unit. The size and type of the mapping population are used to estimate the map accuracy and are integrated into the compilation. Marker names and positions were provided in the input map files, along with a file specifying the synonymous names of the same markers that were mapped in different maps (Additional file 4 and on the SGN database [37]).
QTL meta-analysis
QTL meta-analysis was performed with the MetaQTL software version 1.0 [20]. MetaQTL requires a minimal set of descriptors characterizing each observed QTL: the QTL position, its confidence interval and/or its individual R² value (at least one of them is mandatory), the trait related to the QTL and the size of the QTL mapping population used for the QTL detection. The statistical method implemented in the MetaQTL software hypothesizes that the input mapping studies are independent from each other. QTL mapping studies, which were repeated in time and space, often detected redundant QTLs at the same position for the same trait. In that case, we kept only the QTL with the highest effect (R²) to avoid the attribution of a too strong weight to that QTL in the meta-analysis.
QTLProj command enabled the homothetic projection of the positions and the confidence intervals of the individual QTLs onto the consensus map. It is based on a scaling rule between the positions of the flanking markers of the QTLs on their original maps and the positions of these markers on the consensus map. The MetaQTL software first took into account the confidence interval reported in the study if available, otherwise an estimation of the confidence interval was calculated using the empirical formula proposed by Darvasi and Soller (1997) [53]: CI=530/NxR², where N is the population size and R² the QTL effect as reported in the individual study. This formula generally gives larger confidence interval than the usual interval length of LOD-1 decrease. We used trait ontology to classify and group original trait names according to their relatedness.
QTLClust command performed the clustering of the projected QTLs referring to the same trait on a given chromosome into all possible numbers of hypothetic clusters or "models", i.e., from the model consisting in only one hypothetic cluster to the model consisting in as many clusters as the total number of individual QTLs reported for the chromosome. For a given model, a Gaussian mixture approach was applied to jointly perform a quantitative clustering of the projected QTLs and estimate meta-QTL positions and confidence intervals by maximizing the likelihood of the initial QTL positions. The clustering could only be performed in the genomic regions where at least two QTLs overlapped. If QTLs were single in a genomic region (referred as outlayer QTLs in Additional file 1), they were excluded from the clustering step.
QTLModel command determined the best clustering model based on information-based criteria that were computed for each possible model: AIC (Akaike information criterion), AICc, AIC3, BIC (Bayesian information criterion) and AWE (average weight of evidence) [20]. The best model was the one which criteria values were the lowest relatively to the criteria values of the other possible models. It corresponds to the optimal number of clusters that best explain the observed QTL distribution along the consensus chromosome map. As a result, each meta-QTL position and confidence interval correspond to the consensus position of all the individual QTLs attributed to this meta-QTL, weighted by their individual accuracies and probability of being attributed to the meta-QTL.
A two-round QTL meta-analysis was adopted. First, a meta-analysis was performed by declaring late blight resistance distinct from maturity in the trait ontology; separated meta-QTLs were thus obtained for each trait in the same analysis. In the second round of QTL meta-analysis, both traits were merged into one single "super trait". The purpose of this second round was to investigate whether maturity QTLs tended to cluster with late blight resistance QTLs or not. For convenience, meta-QTLs of late blight resistance were called "late blight meta-QTLs", and meta-QTLs of maturity, vigour and plant height were called "maturity meta-QTLs".
References
Goffinet B, Gerber S: Quantitative trait loci: a meta-analysis. Genetics. 2000, 155 (1): 463-473.
Ballini E, Morel J-Bt, Droc Gt, Price A, Courtois B, Notteghem JL, Tharreau D: A genome-wide meta-analysis of rice blast resistance genes and quantitative trait loci provides new insights into partial and complete resistance. Molecular Plant-Microbe Interactions. 2008, 21 (7): 859-868. 10.1094/MPMI-21-7-0859.
Chardon F, Virlon B, Moreau L, Falque M, Joets J, Decousset L, Murigneux A, Charcosset A: Genetic architecture of flowering time in maize as inferred from quantitative trait loci meta-analysis and synteny conservation with the rice genome. Genetics. 2004, 168 (4): 2169-2185. 10.1534/genetics.104.032375.
Coque M, Martin A, Veyrieras J, Hirel B, Gallais A: Genetic variation for N-remobilization and postsilking N-uptake in a set of maize recombinant inbred lines. 3. QTL detection and coincidences. TAG Theoretical and Applied Genetics. 2008, 117 (5): 729-747. 10.1007/s00122-008-0815-2.
Courtois B, Ahmadi N, Khowaja F, Price AH, Rami JF, Frouin J, Hamelin C, Ruiz M: Rice root genetic architecture: meta-analysis from a drought QTL database. Rice. 2009, 2 (2/3): 115-128. 10.1007/s12284-009-9028-9.
Griffiths S, Simmonds J, Leverington M, Wang Y, Fish L, Sayers L, Alibert L, Orford S, Wingen L, Herry L, Faure S, Laurie D, Bilham L, Snape J: Meta-QTL analysis of the genetic control of ear emergence in elite European winter wheat germplasm. TAG Theoretical and Applied Genetics. 2009, 119 (3): 383-395. 10.1007/s00122-009-1046-x.
Guo B, Sleper DA, Lu P, Shannon JG, Nguyen HT, Arelli PR: QTLs associated with resistance to soybean cyst nematode in soybean: meta-analysis of QTL locations. Crop Science. 2006, 46 (2): 595-602. 10.2135/cropsci2005.04-0036-2.
Hanocq E, Laperche A, Jaminon O, Lainé A, Le Gouis J: Most significant genome regions involved in the control of earliness traits in bread wheat, as revealed by QTL meta-analysis. TAG Theoretical and Applied Genetics. 2007, 114 (3): 569-584. 10.1007/s00122-006-0459-z.
Khowaja F, Norton G, Courtois B, Price A: Improved resolution in the position of drought-related QTLs in a single mapping population of rice by meta-analysis. BMC Genomics. 2009, 10 (1): 276-10.1186/1471-2164-10-276.
Lanaud C, Fouet O, Clement D, Boccara M, Risterucci AM, Surujdeo-Maharaj S, Legavre T, Argout X: A meta-QTL analysis of disease resistance traits of Theobroma cacao L. Molecular Breeding. 2009, 24 (4): 361-374. 10.1007/s11032-009-9297-4.
Liu S, Hall MD, Griffey CA, McKendry AL: Meta-analysis of QTL associated with fusarium head blight resistance in wheat. Crop Science. 2009, 49 (6): 1955-1968. 10.2135/cropsci2009.03.0115.
Löffler M, Schon CC, Miedaner T: Revealing the genetic architecture of FHB resistance in hexaploid wheat (Triticum aestivum L.) by QTL meta-analysis. Molecular Breeding. 2009, 23 (3): 473-488. 10.1007/s11032-008-9250-y.
Marandel G, Salava J, Abbott A, Candresse T, Decroocq V: Quantitative trait loci meta-analysis of Plum pox virus resistance in apricot (Prunus armeniaca L.): new insights on the organization and the identification of genomic resistance factors. Molecular Plant Pathology. 2009, 10 (3): 347-360. 10.1111/j.1364-3703.2009.00535.x.
Norton GJ, Aitkenhead MJ, Khowaja FS, Whalley WR, Price AH: A bioinformatic and transcriptomic approach to identifying positional candidate genes without fine mapping: an example using rice root-growth QTLs. Genomics. 2008, 92 (5): 344-352. 10.1016/j.ygeno.2008.07.002.
Rong J, Feltus FA, Waghmare VN, Pierce GJ, Chee PW, Draye X, Saranga Y, Wright RJ, Wilkins TA, May OL, Smith CW, Gannaway JR, Wendel JF, Paterson AH: Meta-analysis of polyploid cotton QTL shows unequal contributions of subgenomes to a complex network of genes and gene clusters implicated in lint fiber development. Genetics. 2007, 176 (4): 2577-2588. 10.1534/genetics.107.074518.
Shi J, Li R, Qiu D, Jiang C, Long Y, Morgan C, Bancroft I, Zhao J, Meng J: Unraveling the Complex Trait of Crop Yield With Quantitative Trait Loci Mapping in Brassica napus. Genetics. 2009, 182 (3): 851-861. 10.1534/genetics.109.101642.
Truntzler M, Barrière Y, Sawkins M, Lespinasse D, Betran J, Charcosset A, Moreau L: Meta-analysis of QTL involved in silage quality of maize and comparison with the position of candidate genes. Theor Appl Genet. 2010, 121 (8): 1465-1482. 10.1007/s00122-010-1402-x.
Wisser RJ, Balint-Kurti PJ, Nelson RJ: The genetic architecture of disease resistance in maize: a synthesis of published studies. Phytopathology. 2006, 96 (2): 120-129. 10.1094/PHYTO-96-0120.
Arcade A, Labourdette A, Falque M, Mangin B, Chardon F, Charcosset A, Joets J: BioMercator: integrating genetic maps and QTL towards discovery of candidate genes. Bioinformatics. 2004, 20 (14): 2324-2326. 10.1093/bioinformatics/bth230.
Veyrieras JB, Goffinet B, Charcosset A: MetaQTL: a package of new computational methods for the meta-analysis of QTL mapping experiments. BMC Bioinformatics. 2007, 8 (1): 49-10.1186/1471-2105-8-49.
Simko I: Comparative analysis of quantitative trait loci for foliage resistance to Phytophthora infestans in tuber-bearing Solanum species. American Journal of Potato Research. 2002, 79 (2): 125-132. 10.1007/BF02881521.
Åšliwka J: Genetic factors encoding resistance to late blight caused by Phytophthora infestans (Mont.) de Bary on the potato genetic map. Cell Mol Biol Lett. 2004, 9 (4B): 855-867.
Isidore E, van Os H, Andrzejewski S, Bakker J, Barrena I, Bryan GJ, Caromel B, van Eck H, Ghareeb B, de Jong W, van Koert P, Lefebvre V, Milbourne D, Ritter E, van der Voort JR, Rousselle-Bourgeois F, van Vliet J, Waugh R: Toward a marker-dense meiotic map of the potato genome: lessons from linkage group I. Genetics. 2003, 165 (4): 2107-2116.
Van Os H, Andrzejewski S, Bakker E, Barrena I, Bryan GJ, Caromel B, Ghareeb B, Isidore E, de Jong W, van Koert P, Lefebvre V, Milbourne D, Ritter E, van der Voort JNAMR, Rousselle-Bourgeois F, van Vliet J, Waugh R, Visser RGF, Bakker J, van Eck HJ: Construction of a 10,000-marker ultradense genetic recombination map of potato: Providing a framework for accelerated gene isolation and a genomewide physical map. Genetics. 2006, 173 (2): 1075-1087. 10.1534/genetics.106.055871.
Rouppe van der Voort J, Wolters P, Folkertsma R, Hutten R, van Zandvoort P, Vinke H, Kanyuka K, Bendahmane A, Jacobsen E, Janssen R, Bakker J: Mapping of the cyst nematode resistance locus Gpa2 in potato using a strategy based on co-migrating AFLP markers. TAG Theoretical and Applied Genetics. 1997, 95: 874-880. 10.1007/s001220050638.
Feingold S, Lloyd J, Norero N, Bonierbale M, Lorenzen J: Mapping and characterization of new EST-derived microsatellites for potato (Solanum tuberosum L.). TAG Theoretical and Applied Genetics. 2005, 111 (3): 456-466. 10.1007/s00122-005-2028-2.
Ghislain M, Núñez J, del Rosario Herrera M, Pignataro J, Guzman F, Bonierbale M, Spooner D: Robust and highly informative microsatellite-based genetic identity kit for potato. Molecular Breeding. 2009, 23 (3): 377-388. 10.1007/s11032-008-9240-0.
Milbourne D, Meyer RC, Collins AJ, Ramsay LD, Gebhardt C, Waugh R: Isolation, characterisation and mapping of simple sequence repeat loci in potato. Mol Gen Genet. 1998, 259 (3): 233-245. 10.1007/s004380050809.
Tanksley SD, Ganal MW, Prince JP, de Vicente MC, Bonierbale MW, Broun P, Fulton TM, Giovannoni JJ, Grandillo S, Martin GB, Messeguer R, Miller JC, Miller L, Paterson AH, Pineda O, Roeder MS, Wing RA, Wu W, Young ND: High density molecular linkage maps of the tomato and potato genomes. Genetics. 1992, 132: 1141-1160.
Gebhardt C, Ritter E, Barone A, Debener T, Walkemeier B, Schachtschabel U, Kaufmann H, Thompson RD, Bonierbale MW, Ganal MW, Tanksley SD, Salamini F: RFLP maps of potato and their alignment with the homoeologous tomato genome. Theor Appl Genet. 1991, 83 (1): 49-57. 10.1007/BF00229225.
Leister D, Ballvora A, Salamini F, Gebhardt C: A PCR-based approach for isolating pathogen resistance genes from potato with potential for wide application in plants. Nature Genetics. 1996, 14 (4): 421-429. 10.1038/ng1296-421.
Rickert AM, Kim JH, Meyer S, Nagel A, Ballvora A, Oefner PJ, Gebhardt C: First-generation SNP/InDel markers tagging loci for pathogen resistance in the potato genome. Plant Biotechnol J. 2003, 1 (6): 399-410. 10.1046/j.1467-7652.2003.00036.x.
Trognitz F, Manosalva P, Gysin R, Nino-Liu D, Simon R, Herrera MD, Trognitz B, Ghislain M, Nelson R: Plant defense genes associated with quantitative resistance to potato late blight in Solanum phureja x dihaploid S. tuberosum hybrids. Molecular Plant-Microbe Interactions. 2002, 15 (6): 587-597. 10.1094/MPMI.2002.15.6.587.
The PoMaMo Database. [http://www.gabipd.org/projects/Pomamo/]
Costanzo S, Simko I, Christ BJ, Haynes KG: QTL analysis of late blight resistance in a diploid potato family of Solanum phureja x S. stenotomum. TAG Theoretical and Applied Genetics. 2005, 111 (3): 609-617. 10.1007/s00122-005-2053-1.
Villamon FG, Spooner DM, Orrillo M, Mihovilovich E, Perez W, Bonierbale M: Late blight resistance linkages in a novel cross of the wild potato species Solanum paucissectum (series Piurana). Theoretical and Applied Genetics. 2005, 111 (6): 1201-1214. 10.1007/s00122-005-0053-9.
The Sol Genomics Network (SGN) database. [http://solgenomics.net/]
Danan S, Chauvin JE, Caromel B, Moal JD, Pellé R, Lefebvre V: Major-effect QTLs for stem and foliage resistance to late blight in the wild potato relatives Solanum sparsipilum and S. spegazzinii are mapped to chromosome X. TAG Theoretical and Applied Genetics. 2009, 119 (4): 705-719. 10.1007/s00122-009-1081-7.
Bormann CA, Rickert AM, Castillo Ruiz RA, Paal J, Lubeck J, Strahwald J, Buhr K, Gebhardt C: Tagging quantitative trait loci for maturity-corrected late blight resistance in tetraploid potato with PCR-based candidate gene markers. Molecular Plant-Microbe Interactions. 2004, 17 (10): 1126-1138. 10.1094/MPMI.2004.17.10.1126.
Simko I, Costanzo S, Ramanjulu V, Christ BJ, Haynes KG: Mapping polygenes for tuber resistance to late blight in a diploid Solanum phureja x S. stenotomum hybrid population. Plant Breeding. 2006, 125: 385-389. 10.1111/j.1439-0523.2006.01232.x.
Leonards-Schippers C, Gieffers W, Schaferpregl R, Ritter E, Knapp SJ, Salamini F, Gebhardt C: Quantitative resistance to Phytophthora infestans in potato - a case-study for QTL mapping in an allogamous plant-species. Genetics. 1994, 137 (1): 67-77.
Collins A, Milbourne D, Ramsay L, Meyer R, Chatot-Balandras C, Oberhagemann P, De Jong W, Gebhardt C, Bonnel E, Waugh R: QTL for field resistance to late blight in potato are strongly correlated with maturity and vigour. Molecular Breeding. 1999, 5 (5): 387-398. 10.1023/A:1009601427062.
Danan S: Diversité structurale des locus de résistance à Phytophthora infestans chez la pomme de terre et synténie chez les Solanacées. Thèse de Montpellier SupAgro. 2009, 184-
Park TH, Vleeshouwers VGAA, Jacobsen E, Vossen Evd, Visser RGF: Molecular breeding for resistance to Phytophthora infestans (Mont.) de Bary in potato (Solanum tuberosum L.): a perspective of cisgenesis. Plant Breeding. 2009, 128 (2): 109-117. 10.1111/j.1439-0523.2008.01619.x.
Lokossou AA, Park T, Arkel Gv, Arens M, Ruyter-Spira C, Morales J, Whisson SC, Birch PRJ, Visser RGF, Jacobsen E, Vossen EAGvd: Exploiting knowledge of R/Avr genes to rapidly clone a new LZ-NBS-LRR family of late blight resistance genes from potato linkage group IV. Molecular Plant-Microbe Interactions. 2009, 22 (6): 630-641. 10.1094/MPMI-22-6-0630.
Millett B, Bradeen J: Development of allele-specific PCR and RT-PCR assays for clustered resistance genes using a potato late blight resistance transgene as a model. TAG Theoretical and Applied Genetics. 2007, 114 (3): 501-513. 10.1007/s00122-006-0449-1.
Spooner DM, Rodriguez F, Polgar Z, Ballard HE, Jansky SH: Genomic Origins of Potato Polyploids: GBSSI Gene Sequencing Data. Crop Sci. 2008, 48 (Supplement_1): S27-S36.
Mueller LA, Mills AA, Skwarecki B, Buels RM, Menda N, Tanksley SD: The SGN comparative map viewer. Bioinformatics. 2008, 24 (3): 422-423. 10.1093/bioinformatics/btm597.
The Potato Genome Sequencing Consortium (PGSC) website. [http://www.potatogenome.net/]
Rensink W, Lee Y, Liu J, Iobst S, Ouyang S, Buell CR: Comparative analyses of six solanaceous transcriptomes reveal a high degree of sequence conservation and species-specific transcripts. BMC Genomics. 2005, 6 (1): 124-10.1186/1471-2164-6-124.
Lokossou AA, Rietman H, Wang M, Krenek P, van der Schoot H, Henken B, Hoekstra R, Vleeshouwers VGAA, van der Vossen EAG, Visser RGF, Jacobsen E, Vosman B: Diversity, Distribution, and Evolution of Solanum bulbocastanum Late Blight Resistance Genes. Molecular Plant-Microbe Interactions. 2010, 23 (9): 1206-1216. 10.1094/MPMI-23-9-1206.
Quraishi UM, Abrouk M, Murat F, Pont C, Foucrier S, Desmaizieres G, Confolent C, Rivière N, Charmet G, Paux E, Murigneux A, Guerreiro L, Lafarge S, Le Gouis J, Feuillet C, Salse J: Cross-Genome Map Based Dissection of a Nitrogen Use Efficiency Ortho-metaQTL in Bread Wheat Unravels Concerted Cereal Genome Evolution. Plant J. 2011, Accepted manuscript online: 13 DEC 2010
Darvasi A, Soller M: A simple method to calculate resolving power and confidence interval of QTL map location. Behavior Genetics. 1997, 27 (2): 125-132. 10.1023/A:1025685324830.
Ewing EE, Simko I, Smart CD, Bonierbale MW, Mizubuti ESG, May GD, Fry WE: Genetic mapping from field tests of qualitative and quantitative resistance to Phytophthora infestans in a population derived from Solanum tuberosum and Solanum berthaultii. Molecular Breeding. 2000, 6 (1): 25-36. 10.1023/A:1009648408198.
Bradshaw JE, Bryan GJ, Hackett CA, McLean K, Pande B, Stewart HE, Waugh R: Dissection and analysis of quantitative disease resistance in tetraploid potato. Euphytica. 2004, 137 (1): 13-18. 10.1023/B:EUPH.0000040499.19190.fb.
Visker M, Keizer LCP, Van Eck HJ, Jacobsen ELT, Colon LT, Struik PC: Can the QTL for late blight resistance on potato chromosome 5 be attributed to foliage maturity type? . TAG Theoretical and Applied Genetics. 2003, 106: 317-325.
Sliwka J, Jakuczun H, Lebecka R, Marczewski W, Gebhardt C, Zimnoch-Guzowska E: Tagging QTLs for late blight resistance and plant maturity from diploid wild relatives in a cultivated potato (Solanum tuberosum) background. TAG Theoretical and Applied Genetics. 2007, 115 (1): 101-112. 10.1007/s00122-007-0546-9.
Visker M, Heilersig H, Kodde LP, Van de Weg WE, Voorrips RE, Struik PC, Colon LT: Genetic linkage of QTLs for late blight resistance and foliage maturity type in six related potato progenies. Euphytica. 2005, 143 (1-2): 189-199. 10.1007/s10681-005-3444-8.
Hein I, Squires J, Birch P, Bryan G: Screening wild potato accessions for resistance to the virulent allele of the Phytophthora infestans effector Avr3a. 13th International Congress on Molecular Plant-Microbe Interactions. Sorrento, Italy;2007, 259-2007
Tan MYA, Hutten RCB, Celis C, Park TH, Niks RE, Visser RGF, van Eck HJ: The Rpi-mcd1 locus from Solanum microdontum involved in resistance to Phytophthora infestans, causing a delay in infection, maps on potato chromosome 4 in a cluster of NBS-LRR genes. Molecular Plant-Microbe Interactions. 2008, 21 (7): 909-918. 10.1094/MPMI-21-7-0909.
Bhaskar P, Raasch J, Kramer L, Neumann P, Wielgus S, Austin-Phillips S, Jiang J: Sgt1, but not Rar1, is essential for the RB-mediated broad-spectrum resistance to potato late blight. BMC Plant Biology. 2008, 8 (1): 8-10.1186/1471-2229-8-8.
Lefebvre V, Chèvre AM: Tools for marking plant disease and pest resistance genes - a review. Agronomie. 1995, 15 (1): 3-19. 10.1051/agro:19950101.
Pajerowska-Mukhtar KM, Mukhtar MS, Guex N, Halim VA, Rosahl S, Somssich IE, Gebhardt C: Natural variation of potato allene oxide synthase 2 causes differential levels of jasmonates and pathogen resistance in Arabidopsis. Planta. 2008, 228 (2): 293-306. 10.1007/s00425-008-0737-x.
Oberhagemann P, Chatot-Balandras C, Schäfer-Pregl R, Wegener D, Palomino C, Salamini F, Bonnel E, Gebhardt C: A genetic analysis of quantitative resistance to late blight in potato: towards marker-assisted selection. Molecular Breeding. 1999, 5 (5): 399-415. 10.1023/A:1009623212180.
Salvi S, Tuberosa R, Chiapparino E, Maccaferri M, Veillet S, van Beuningen L, Isaac P, Edwards K, Phillips RL: Toward positional cloning of Vgt1, a QTL controlling the transition from the vegetative to the reproductive phase in maize. Plant Molecular Biology. 2002, 48 (5): 601-613. 10.1023/A:1014838024509.
Visser R, Bachem C, de Boer J, Bryan G, Chakrabati S, Feingold S, Gromadka R, van Ham R, Huang S, Jacobs J, Kuznetsov B, de Melo P, Milbourne D, Orjeda G, Sagredo B, Tang X: Sequencing the potato genome: Outline and first results to come from the elucidation of the sequence of the world's third most important food crop. American Journal of Potato Research. 2009, 86 (6): 417-429. 10.1007/s12230-009-9097-8.
Coudert Y, Périn C, Courtois B, Khong NG, Gantet P: Genetic control of root development in rice, the model cereal. Trends in Plant Science. 2010, 15 (4): 219-226. 10.1016/j.tplants.2010.01.008.
Bradshaw JE, Hackett CA, Lowe R, McLean K, Stewart HE, Tierney I, Vilaro MDR, Bryan GJ: Detection of a quantitative trait locus for both foliage and tuber resistance to late blight [Phytophthora infestans (Mont.) de Bary] on chromosome 4 of a dihaploid potato clone (Solanum tuberosum subsp tuberosum). Theoretical and Applied Genetics. 2006, 113 (5): 943-951. 10.1007/s00122-006-0353-8.
Ghislain M, Trognitz B, Herrera MD, Solis J, Casallo G, Vasquez C, Hurtado O, Castillo R, Portal L, Orrillo M: Genetic loci associated with field resistance to late blight in offspring of Solanum phureja and S. tuberosum grown under short-day conditions. Theoretical and Applied Genetics. 2001, 103 (2-3): 433-442. 10.1007/s00122-001-0545-1.
Meyer RC, Milbourne D, Hackett CA, Bradshaw JE, McNichol JW, Waugh R: Linkage analysis in tetraploid potato and association of markers with quantitative resistance to late blight (Phytophthora infestans). Mol Gen Genet. 1998, 259 (2): 150-160. 10.1007/s004380050800.
Naess SK, Bradeen JM, Wielgus SM, Haberlach GT, McGrath JM, Helgeson JP: Resistance to late blight in Solanum bulbocastanum is mapped to chromosome 8. Theoretical and Applied Genetics. 2000, 101 (5-6): 697-704. 10.1007/s001220051533.
Sandbrink JM, Colon LT, Wolters P, Stiekema WJ: Two related genotypes of Solanum microdontum carry different segregating alleles for field resistance to Phytophthora infestans. Molecular Breeding. 2000, 6 (2): 215-225. 10.1023/A:1009697318518.
Sorensen KK, Madsen MH, Kirk HG, Madsen DK, Torp AM: Linkage and quantitative trait locus mapping of foliage late blight resistance in the wild species Solanum vernei. Plant Breeding. 2006, 125 (3): 268-276. 10.1111/j.1439-0523.2006.01219.x.
Bonierbale MW, Plaisted RL, Tanksley SD: RFLP maps based on a common set of clones reveal modes of chromosomal evolution in potato and tomato. Genetics. 1988, 120 (4): 1095-1103.
Yamanaka S, Ikeda S, Imai A, Luan Y, Watanabe JA, Watanabe KN: Construction of integrated genetic map between various existing DNA markers and newly developed P450-related PBA markers in diploid potato (Solanum tuberosum). Breeding Science. 2005, 55 (2): 223-230. 10.1270/jsbbs.55.223.
Acknowledgements
This work was financially supported by the BIOEXPLOIT Integrated Project FOOD CT2005-513959 that resides under the 6th framework programme of the European Union. In addition, it was approved by the PEIFL (European fruits and vegetables innovative cluster), which is a competitiveness cluster at French national level. Sarah Danan received a PhD fellowship funded by the European BIOEXPLOIT project. Authors warmly thank Lukas Mueller for his useful critical review of this manuscript.
Author information
Authors and Affiliations
Corresponding author
Additional information
Authors' contributions
VL conceived of the initial idea and coordinated the study. All authors participated in the conception of the study. JBV originally coded the package MetaQTL and gave a significant support during the meta-analysis process. SD carried out the collection of data and applied the QTL meta-analysis. VL and SD contributed to the interpretation of the results and drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
12870_2010_801_MOESM1_ESM.XLS
Additional file 1: Description of the consensus map with meta-QTLs positions, additional individual QTLs, Rpi-genes, RGAs and DRLs. The chromosome Arabic numbers, the locus names, and the positions of the loci on the consensus map (in cM, Haldane unit) are listed in the ascending order of marker positions, from chromosome I to chromosome XII. The occurrence number of the markers across the 21 compiled maps is indicated (meta.occurrence). The "type" information indicates whether the locus is a marker (M) or a meta-QTL (Q). If the marker sequence has been previously described as being a resistance gene analog (RGA) or a defence-related locus (DRL) (information mainly retrieved from the PoMaMo database), the type "M_RGA" or "M_DRL" is mentioned; in this case, a short description about the sequence is added. Rpi-gene (R-gene to Phytophthora infestans) most probable locations are indicated by shaded areas between their two closest flanking markers, as described in literature. Unless markers were designed from known Rpi-gene sequences (as for RB gene), Rpi-genes have most of the time an approximate location on the consensus map. Overlapping locations of Rpi-genes in the same region are indicated in two different columns of the table. Rpi-genes with different names are considered as being different genes, even if they are alleles [44].The columns "qtl.ci.from" and "qtl.ci.to" indicate the range of the meta-QTL confidence interval on the consensus map. In the table, the confidence intervals of late blight meta-QTLs are framed in blue while those of maturity meta-QTLs are framed in red, and overlapping regions in purple. The column "meta-QTL trait" specifies the meta-QTL trait: late blight resistance (Late_blight) or maturity (Maturity).The columns "anchored QTL, trait" and "outlayer QTL, trait" show additional QTLs that have not been compiled in the final meta-analysis. The QTL nomenclature is described in the Figure 2 legend. For "anchored QTLs", the final number is the chromosome number alone; for "outlayer QTLs", the final number is the chromosome number juxtaposed to the QTL mapping order on the chromosome and a letter was sometimes added to distinguish colocalizing QTLs that were detected with different traits. The trait with which the QTL was detected is indicated after the comma by the trait code, as defined in the Table 2 legend."Anchored QTLs" are QTLs from non-compiled maps and they could only be anchored to the consensus map if the confidence interval comprised a common marker; they are mentioned at the position of the common marker. No confidence interval information is available."Outlayer QTLs" are QTLs from compiled maps and, as they were alone in their region on the consensus map, they eventually could not be clustered with other QTLs in a meta-QTL. Confidence intervals of these "outlayer QTLs" are represented by dotted areas. (XLS 308 KB)
12870_2010_801_MOESM2_ESM.XLS
Additional file 2: Criterion values of the meta-QTL models. The values of the AIC, AICc, AIC3, BIC and AWE criteria provided by the MetaQTL software are listed for the twelve chromosomes (first column), for both maturity and late blight resistance traits (second column), and for all possible meta-QTL models (third column, K is the number of hypothetic meta-QTLs of a model). "Delta" is the rescaled value of the criterion, it is the difference between the value of the given model and the value of the best model. The weight is the "weight of evidence" of the model, the higher it is, the more confidence we can have in the corresponding model. (XLS 134 KB)
12870_2010_801_MOESM3_ESM.XLS
Additional file 3: Meta-QTL details: number per chromosome, confidence intervals and composition in published QTLs. For each chromosome, late blight and maturity meta-QTLs are described. Meta-QTL numbers, names and confidence intervals (CI, in cM, Haldane unit) are specified. The QTL composition of each meta-QTL is detailed by giving the individual QTL names, the species origins of the mapping parents (when pedigrees are known), R² values (if specified in the original publication) and confidence intervals as described in the original publication (in their original unit, Haldane or Kosambi) with their means for each meta-QTL. The confidence intervals of the projected QTLs with their means are also indicated. The QTL nomenclature is described in the Figure 2 legend. For Villamon et al. and Leonards-Schippers et al.' studies, the original name of the QTL was added [36, 41]. When QTLs are equally shared by two meta-QTLs, they are shown in italic ('LB' for late blight resistance). The mean of the original confidence intervals of the QTLs and the number of the QTLs composing the meta-QTL are given for each meta-QTL. NA: not available. (XLS 72 KB)
12870_2010_801_MOESM4_ESM.XLS
Additional file 4: Synonymous names of markers described in the compiled published maps. Marker names that were the most often encountered for a marker and which had a familiar nomenclature were considered as standard names. Their corresponding synonymous marker names that were found in the compiled maps were collected. Only standard names were used in the final consensus map. (XLS 49 KB)
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
Danan, S., Veyrieras, JB. & Lefebvre, V. Construction of a potato consensus map and QTL meta-analysis offer new insights into the genetic architecture of late blight resistance and plant maturity traits. BMC Plant Biol 11, 16 (2011). https://doi.org/10.1186/1471-2229-11-16
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1471-2229-11-16