Skip to main content

Cytoplasmic genome types of European potatoes and their effects on complex agronomic traits

Abstract

Background

Various wild species germplasm has been used in European potato breeding since the first introduction of potato (Solanum tuberosum L.) to Europe. As the plant cytoplasmic genome including chloroplast and mitochondrial genomes is transmitted only through the maternal parent, cytoplasmic markers are useful tools in breeding programs to determine cytoplasmic genome types and to trace maternal ancestors. The potato cytoplasmic genome can be distinguished into six distinct types (M, P, A, W, T, and D). Male sterility was found in genotypes with S. demissum-derived D-type cytoplasm and S. stoloniferum-derived W/γ-type cytoplasm. These wild species were frequently used to incorporate superior pathogen resistance genes. As a result, the percentage of these two types is increasing unintentionally in the European germplasm pool. Other than cytoplasmic male sterility, little is known about effects of the cytoplasmic genome on complex agronomic traits in potato.

Result

The cytoplasm types of 1,217 European potato cultivars and breeding clones were determined with type specific DNA markers. Most frequent were T- (59.4 %), D- (27.4 %), and W- (12.2 %) type cytoplasm, while A- (0.7 %) and M-type cytoplasm (0.3 %) was rare and P-type cytoplasm was absent. When comparing varieties with breeding clones, the former showed a relatively higher frequency of T-type and lower frequency of D- and W-type cytoplasm. Correlation analysis of cytoplasm types and agronomic data showed that W/γ-type cytoplasm was correlated with increased tuber starch content and later plant maturity. Correlation with quantitative resistance to late blight was observed for D-type and M-type cytoplasm. Both cytoplasm types had a positive effect on resistance.

Conclusion

This study revealed and quantified the cytoplasmic diversity in the European potato germplasm pool. Knowledge of cytoplasm type is important for maintaining genetic diversity and managing the male sterility problem in breeding programs. This is the first comprehensive study to show correlations of distinct cytoplasmic genomes with complex agronomic traits in potato. Correlations particularly with tuber starch content and resistance to late blight provided new knowledge on cytoplasmic effects on these important traits, which can be exploited for genetic improvement of potato.

Background

In a plant cell, multiple copies of chloroplast and mitochondrial DNA co-exist with one copy of nuclear DNA. Thus, the expression and function of nuclear genes should be affected in various ways by chloroplast and mitochondrial genes. A typical phenomenon caused by an interaction between nuclear and mitochondrial genes is cytoplasmic male sterility frequently found in crop species [1]. Chloroplast and mitochondrial DNA compose the cytoplasmic genome and are maternally inherited in most angiosperms [2, 3]. Various effects of different cytoplasmic genomes on agronomic traits have been demonstrated in classical work in Triticum and Aegilops using cytoplasm substitution lines [4].

Potato (Solanum tuberosum L.) is a crop that possesses different types of cytoplasmic genomes within a cultivar group. Since its first introduction into Europe in the sixteenth century, many diseases have threatened potato cultivation. In order to overcome them, new genetic resources have been frequently utilized from cultivated Andean landraces and wild potato species [5, 6]. The first introduction into Europe is thought to have been S. tuberosum L. ssp. andigena (referred to as S. tuberosum L. Andigenum Group by Spooner et al. [7]. In this article we tentatively use the taxonomy of Hawkes [8]). Later Chilean forms (S. tuberosum L. ssp. tuberosum) were introduced which became the ancestors of modern cultivars improved for short stolons, early vine maturity and high tuber yield in Europe and North America [5, 6, 8, 9]. The first dramatic change in the cultivar spectrum happened with the late blight epidemics caused by Phytophthora infestans, resulting in the Irish Famine (1845–1847). Soon after the Famine, modern potato breeding started using only a few Chilean ssp. tuberosum clones including the cultivar Rough Purple Chili and an ancestor of Alte Daber [5, 6, 10]. In 1906, Salaman [11] detected resistance to late blight in S. × edinense Berth. and crossed it with ssp. tuberosum. In 1908, the Mexican hexaploid species S. demissum Lindl. was introduced, which then started the breeding for late blight resistance with S. demissum germplasm in Germany [5]. The resultant pentaploid hybrid could be backcrossed easily with S. tuberosum. More beneficial tetraploid hybrids were obtained by crossing the diploid cultivated species S. phureja Juz. et Buk. with S. demissum [5, 12]. These efforts resulted in more than 80 % of modern cultivars of Germany carrying genes from S. demissum [5]. Many aphid transmitted viruses also damage severely potato cultivation. Plants infected by Potato virus Y (PVY) and Potato virus X (PVX) are badly harmed. In Europe, extreme resistance to PVY was first detected by Stelzer in S. stoloniferum Schltdl. [13]. A dominant major gene Ry sto was found by Ross [14] and Cockerham [15]. A major gene Rx acl for resistance to PVX was found in S. acaule Bitt. [15, 16]. These resistance genes, sometimes found in complex hybrids such as those from a cross (S. acaule × S. stoloniferum) × S. tuberosum, were introduced into many parental lines in Germany. Potato cyst nematodes have raged throughout Europe probably since the mid-nineteenth century. The gene H1 for extreme resistance to the potato cyst nematode Globodera rostochiensis (Woll.) was detected in ssp. andigena accession CPC 1673 in the Commonwealth Potato Collection [1719]. Resistance to G. pallida pathotype Pa2/3 was introgressed into breeding lines from the wild potato species S. vernei Bitt. et Wittm. and S. spegazzinii Bitt.. Nowadays, most modern German cultivars have one or more nematode resistance genes from these species. After a seriously raging wart epidemic in 1910, caused by the fungus Synchytrium endobioticum, genes for resistance to S. endobioticum were found in some wild species [20, 21] and have been used for a long time for successful prevention [22, 23]. As briefly summarized above, resistance breeding in Europe has an intricate history and is based on the frequent use of wild species germplasm. Consequently, various cytoplasmic genomes are expected to be present in European potatoes. However, it is difficult to identify cytoplasmic genomes by tracing back the maternal lineage, often because no information is available about which parent was used as the female parent in a cross. Sometimes also no pedigree record is available for breeding clones.

Comparing reciprocal hybrid populations, it has long been known that the cytoplasm of S. tuberosum ssp. tuberosum is different from that of S. tuberosum ssp. andigena, the former inducing higher percentage of tuberization, higher tuber yield, higher tuber numbers, and earlier vine maturity [2428]. Several cytoplasmic genomes were distinguished among cultivated potatoes and its related wild species [29]. Polymerase chain reaction (PCR)-based markers were developed that distinguish S. tuberosum ssp. tuberosum-type chloroplast DNA from the other chloroplast DNA types [30, 31]. Lössel et al. [30] also developed PCR markers that distinguish three mitochondrial DNA types (α-, β-, and γ-types). Recently, Hosaka and Sanetomo [32] developed a simpler and more informative technique using a set of five cytoplasmic markers (four chloroplast and one mitochondrial DNA markers), which differentiate six potato cytoplasm types: M, P, A, W, T, and D. The P- and A-type cytoplasm and the T- and D-type cytoplasm are relatively distinct types within the M- and W-type cytoplasm, respectively, each of which has diverse cytoplasmic variations [33, 34]. Andean cultivated potatoes evolved from ancestral wild species with M or M-derived type cytoplasm, while all other wild species not involved in the origin of cultivated potatoes have W or W-derived type cytoplasm [32]. The A-type cytoplasm is most prevalent in S. tuberosum ssp. andigena, while the T-type cytoplasm is most prevalent in S. tuberosum ssp. tuberosum. The P-type cytoplasm was introduced from S. phureja [35], while the D-type cytoplasm was introduced from S. demissum into the common potato gene pool [36]. This cytoplasm type definition system is validated only among cultivated potatoes and their close wild relatives [37]. Many wild species have specific cytoplasmic genomes [34], all of which are assigned as W-type cytoplasm. In order to distinguish the S. stoloniferum- derived cytoplasm carried by many modern varieties, an additional mitochondrial marker ALM_4/ALM_5, developed by Lössel et al. [30] is needed, by which three mitochondrial types, α-, β-, and γ-types can be distinguished. The S. stoloniferum-derived cytoplasm is characterized as W/γ subtype [30].

Cytoplasmic male sterility was found in interspecific crosses in potato, realizing the existence of different cytoplasm, as has been known in other crop species [1, 38, 39]. The common potato cytoplasm induces various types of intrinsic sterility [40]. Cultivars carrying the PVY resistance gene Ry sto exhibit complete male sterility caused by interaction with mitochondrial DNA of S. stoloniferum [41, 42], which is characterized as W/γ subtype cytoplasm [30, 43]. Sterile pollen grains clumped together in tetrads, so it was called “tetrad sterility” [44] or “lobed sterility” [45]. The same type of sterility was also observed with S. verrucosum-derived cytoplasm [44]. Pentaploid F1 hybrids can be easily obtained when S. demissum is crossed as a female parent with S. tuberosum. The hybrids produce normal-looking pollen, however, they are non-functional as male parents. The progeny produced by continued backcrossing with the pollen of S. tuberosum can be used only as female parents, although these plants produce abundant stainable pollen [46]. Thus, the S. demissum-derived cytoplasm is also associated with functionally male sterile pollen. Once S. stoloniferum or S. demissum cytoplasm are incorporated into parental lines, they can be used only as female parents. Continued infiltration of the potato gene pool by these cytoplasm would worsen male sterility problems as warned by Provan et al. [47] and Hosaka and Sanetomo [32].

Marker-assisted selection is an efficient breeding tool that connects genotypes with agronomic traits and pathogen resistances. Various diagnostic DNA markers are available now for potato breeding [35, 4850]. Recently, association genetics has been applied to identify diagnostic markers for quantitative traits that are controlled by multiple genes and environmental factors. For example, associations were discovered between DNA polymorphisms at individual candidate loci and complex traits such as tuber yield, starch and sugar content [5154]. Gebhardt et al. [55] genotyped a gene bank collection of 600 potato cultivars with five DNA markers linked to a previously mapped quantitative trait locus (QTL) for resistance to late blight and plant maturity. Significant association with quantitative resistance to late blight and plant maturity was detected with PCR markers derived from R1, a major gene for race specific resistance to late blight, or tightly linked to R1. The marker alleles associated with increased resistance and later maturity were traced to an introgression from S. demissum [55]. Pajerowska-Mukhtar et al. [56] tested 24 candidate loci for association with field resistance to late blight and plant maturity in a population of 184 breeding clones and found single nucleotide polymorphisms (SNPs) in the Allene Oxide Synthase 2 (StAOS2) gene associated with field resistance to late blight.

Using cytoplasmic markers, Lössel et al. [30] indicated that W/α and W/γ-type cytoplasm showed a higher tuber starch content than T/β-type cytoplasm. Apart from that, little is known about effects of the cytoplasmic genome on agronomic performance, mainly because an accurate method to distinguish cytoplasmic genomes was not available until recently.

In this study, we analyzed 1,383 tetraploid genotypes of six different populations to disclose the cytoplasmic diversity in European potato gene pool. These populations have been previously evaluated for agronomic traits such as late blight resistance, chip quality, tuber yield and starch content, plant maturity and susceptibility to tuber bruising in the context of searching for associations with nuclear markers. Correlations were investigated between different cytoplasm types and agronomic traits. The importance of cytoplasmic diversity and the correlation with some agronomic traits, especially with tuber starch content and resistance to late blight are discussed.

Results

Cytoplasm types of European potato collections

A total of 1,383 tetraploid cultivars and breeding clones of six populations (Table 1) were genotyped using multiplex PCR with the cytoplasmic markers T, S, SAC, D, and A. One genotype of population PIN184 showed a mixed pattern of T- and M-type cytoplasm, while another genotype of population CHIPS-ALL showed a mixed pattern of T- and D-type cytoplasm, probably due to DNA contaminations. These two genotypes were discarded for the further analysis. Genotypes with W-type cytoplasm were further examined using the mtDNA (mitochondrial DNA) marker ALM_4/ALM_5 which distinguished four different subtype cytoplasm: W/α, W/β, W/γ, and the fourth type. The fourth type detected in one cultivar had both 2.4 kb and 1.6 kb bands (= Type 3 banding pattern reported by Hosaka and Sanetomo [32]), which is designated as W/αβ-type cytoplasm in this article.

Table 1 Populations used and quantitative agronomic trait data evaluated previously

The T-type cytoplasm was the most prevalent in all six populations (Table 2, Additional file 1: Table S1). The GBC population consisting of 536 genotypes included many old varieties. 369 varieties and 6 breeding clones had T-type cytoplasm, being the highest frequency of T-type cytoplasm (70.0 %) among all populations. On the other hand, T-type cytoplasm was found in less than half of the genotypes (45.9 %) in the PIN184 population, which represented modern breeding materials. In contrast, D-type cytoplasm was found with the highest frequency in the PIN184 population (35.5 %) and the lowest frequency in the GBC population (20.3 %). Within W-type cytoplasm, the subtype W/γ was the most frequent. The frequency of W/γ-type cytoplasm in the EURO-CUL and GBC populations was 11.6 % and 6.0 %, respectively, which was lower than those of the other four populations (13.7 % − 25.0 %). The variety “Raisa” and two breeding clones “CIP 38 31 17 06” and “MPI79.452/14D” had the W/α-type cytoplasm. Eighteen genotypes had the W/β-type. Only the variety “Rita”, which was included in both EURO-CUL and GBC populations, had W/αβ-type cytoplasm. A-type cytoplasm was found in four genotypes in EURO-CUL, five genotypes in GBC and one genotype in PIN184. All four genotypes with M-type cytoplasm were from the PIN184 population.

Table 2 The number of genotypes and percentages with different cytoplasm types in each population

One hundred and nineteen varieties and one breeding line were duplicated in at least two populations. Of these duplicated genotypes, 104 had the identical cytoplasm type in all duplicates, whereas 16 (13 %) had different cytoplasm types (Additional file 2: Table S2), probably due to sampling errors. Discarding these duplicates, a total of 1,217 varieties and breeding clones were actually determined for the cytoplasm types: 723 (59.4 %) had T-type, 333 (27.4 %) D-type, 9 (0.7 %) A-type, 4 (0.3 %) M-type and 148 (12.2 %) had W-type cytoplasm (Table 2). None of the genotypes had the P-type cytoplasm. Of 1,217 genotypes, 694 were named cultivars, of which 480 (69.2 %) had T-type, 148 (21.3 %) D-type, 8 (1.2 %) A-type and 58 (8.4 %) had W-type cytoplasm (W/γ = 6.5 %). The remaining 523 genotypes were breeding clones, of which 243 (46.5 %) had T-type, 185 (35.4 %) D-type, 1 (0.2 %) A-type, 4 (0.8 %) M-type and 90 (17.2 %) had W-type cytoplasm (W/γ = 15.5 %). The cultivars “Tannenzapfen” and “Pink Fir Apple” and the three cultivars “Asparges”, “Corne de Bique” and “La Ratte” were treated as different cultivars, but are suspected to be identical, according to the German potato gene bank at Groß-Lüsewitz and the Potato Pedigree Database at Wageningen [57]. All five cultivars had the A-type cytoplasm in common. The cytoplasm types of all 694 European varieties, 26 named breeding clones and 497 breeding clones with ID number of each breeding company are listed in Additional file 1: Table S1.

Correlation of different cytoplasm types with quantitative agronomic traits

The different cytoplasm types were tested for correlation with 20 quantitative agronomic traits that have been evaluated in five of the six populations using one-way ANOVA (parametric) or Welch’s test (nonparametric) (Table 3). The cytoplasm type found in less than three genotypes in each population (Table 2) was omitted because it was a number too small to apply statistical analysis by ANOVA or Welch’s test.

Table 3 One-way ANOVA or Welch’s test of cytoplasm types with phenotypic scores for agronomic traits

For the significant traits, means of each pair were compared using Tukey’s test performed after ANOVA or Kruskal-Wallis test performed after Welch’s test.

In the BRUISE population (n = 203), T-, D- and W/γ-type cytoplasm were tested for correlation with the phenotypic traits BI (bruising index), SCB (starch corrected bruising), PM (plant maturity), TS (tuber shape), TSC (tuber starch content), and TY (tuber yield). For five of the six traits significant differences among cytoplasm types were found (Table 3). Starch corrected bruising was not significantly different among cytoplasm types. Genotypes with W/γ-type cytoplasm matured significantly later, and had more round tuber shape compared to those with T- or D-type cytoplasm. Compared to genotypes with T-type cytoplasm, genotypes with W/γ-type cytoplasm were more susceptible to black spot bruising, had higher tuber starch content and yield.

Chip quality was analyzed as CQA (chip quality after harvest) and CQS (chip quality after 3 months storage at 4 °C) in the CHIPS-ALL population (n = 227) and as RSC (tuber reducing sugar content) in the SUGAR40 population (n = 40). No significant difference was found among cytoplasm types for the three chip quality traits. Of the traits TSC, TY and TSY (tuber starch yield) analyzed in the CHIPS-ALL population, only TSC showed significant differences between cytoplasm types. Consistent with the observations in the BRUISE population, TSC was significantly higher in genotypes with W/γ-type cytoplasm compared to those with T- or D-type cytoplasm.

In the GBC population (n = 536), passport data for RLBF (foliage resistance to late blight) and RLBT (tuber resistance to late blight) were analyzed for correlation with cytoplasm types. Welch’s test revealed a significant difference for RLBF. Genotypes with T- and A-type cytoplasm showed lower resistance levels compared to those with D-, W/γ and W/β-type cytoplasm. PM did not show significant differences among cytoplasm types.

In the PIN184 population (n = 182), T, D, M, W/β and W/γ cytoplasm types were analyzed for correlation with resistance to late blight measured as MCR (maturity corrected resistance) and rAUDPC (relative area under disease progress curve). The results showed that cytoplasm type had a clearly significant effect on resistance to late blight (Table 3). Genotypes with D-type cytoplasm showed a significantly lower mean MCR value compared to those with T or W/γ-type cytoplasm (Table 3). The four genotypes with M-type cytoplasm had the lowest mean MCR value (−0.13). Thus, both M- and D-type cytoplasm were correlated with increased resistance to late blight. The same was true for rAUDPC, which showed higher resistance levels to late blight with M- and D-type cytoplasm compared to those with T-type cytoplasm. Mean PM and TSC also differed significantly among genotypes with different cytoplasm types. Genotypes with W/β-type cytoplasm had the latest maturity but the highest tuber starch content (21.8 %) compared to those with the other cytoplasm type. Genotypes with W/γ-type cytoplasm matured later and also had higher starch content (17.9 %) compared to those with T-type cytoplasm (16.0 %).

Thus, compared with the T-type cytoplasm, the D-type cytoplasm was correlated with increased foliage resistance to late blight in two independent populations (GBC and PIN184). The W/γ-type cytoplasm was correlated with higher tuber starch content in three populations (BRUISE, CHIPS-ALL and PIN184) and with later maturity in two populations (BRUISE and PIN184).

Correlation of D-type cytoplasm with nuclear gene markers for late blight resistance

The fact that genotypes with D-type cytoplasm showed a higher average level of late blight resistance compared to those with T-type cytoplasm in both the GBC and PIN184 populations (Table 3), could result from the joint introgression of D-type cytoplasm with R genes from S. demissum. We tested therefore whether the presence of the marker diagnostic for D-type cytoplasm was correlated with the presence of nuclear markers closely linked or identical with the late blight resistance genes R1, R3a and R3b.

Correlation coefficients were obtained for the D-type cytoplasmic marker with four markers either located in the R1 resistance gene (R11400, R11800) or tightly linked to R1 (CosA210 and GP179570) that have been scored in the GBC population [55], and with one R1 diagnostic marker (CosA210), and R3a and R3b gene specific markers scored in the PIN184 population [56]. None of nuclear markers for R genes showed significant correlation with the D cytoplasmic marker neither in the GBC nor the PIN184 population (Table 4).

Table 4 Pearson’s correlation coefficients between D-type cytoplasm and the nuclear markers linked with resistance genes R1, R3a and R3b

Furthermore, the frequencies of the StAOS2_A 691 C 692 haplotype were analyzed with different cytoplasm types by a Welch’s test in the PIN184 population (Table 5). The StAOS2_A 691 C 692 haplotype was associated with higher late blight resistance [56]. Significant difference was found at a 5 % level. Kruskal-Wallis test indicated that the genotypes with D- and M-type cytoplasm had higher haplotype frequencies compared to those with T-type cytoplasm.

Table 5 Welch’s test for cytoplasmic differences on allele frequencies of SNP haplotype StAOS2_A 691 C 692 in the PIN184 population

Combined effects of nuclear markers and D-type cytoplasm on late blight resistance and plant maturity

Since no correlation was found between R gene nuclear markers and D-type cytoplasm, their combined effects on late blight resistance and plant maturity were evaluated. Four marker classes +/+, +/−, −/+, and −/− were tested by Welch’s test for significant differences, where ‘plus’ indicates the presence and ‘minus’ the absence of the marker, irrespective of allele dosage.

In the GBC population, the four marker classes of all combinations of nuclear markers with D-type cytoplasm differed highly significantly for RLBF (P < 0.001) (Table 6). Two-way ANOVA was performed using the nuclear markers and D-type cytoplasm as two contributing factors. Only D-type cytoplasm was significant (P < 0.001), indicating that genotypes with the D-type cytoplasm were more resistant to foliage late blight. Interaction of the D-type cytoplasm with GP179570 was found at the 5 % significance level. For RLBT, the marker classes combining nuclear markers R11400, R11800 and CosA210 with the D-type cytoplasm were significantly different. However, only the nuclear markers were significant factors. The presence of R11400 and CosA210 indicated higher resistance to tuber late blight, while the presence of R11800 indicated lower resistance to tuber late blight. For PM, significant difference was not found in any nuclear marker and D-type cytoplasm combination.

Table 6 Interaction for effects of combinations of presence (+) and absence (−) of nuclear markers with presence (+) or absence (−) of D-type cytoplasm on resistance to late blight and plant maturity in the GBC population. In case differences were found among marker classes, two-way ANOVA was performed to explore the factors contributing the differences

In the PIN184 population, mean MCR and rAUDPC were both significantly different among marker classes combining D-type cytoplasm with either R3a, R3b, or CosA210 markers (Table 7). By two-way ANOVA, D-type cytoplasm was the significantly contributing factor in these combinations and indicated a higher level of resistance. Some additional effects of interactions of D-type cytoplasm with CosA210 for MCR were detected. Combination of D-type cytoplasm with CosA210 affected PM. D-type cytoplasm and its interaction with CosA210 were contributing factors; the presence of D-type cytoplasm resulted in later maturity, and in combination with CosA210 resulted in the latest maturity, whereas the presence of CosA210 without D-type cytoplasm resulted in the earliest maturity. The dosage classes of the SNP haplotype StAOS2_A 691 C 692 were grouped in two genotype classes, one lacking the haplotype StAOS2_A 691 C 692 and the other with the StAOS2 A 691 C 692 haplotype. The two StAOS2 marker classes were combined with presence or absence of the D-type cytoplasm and analyzed for effects on MCR, rAUDPC, and PM (Table 7). Significant differences among marker classes were found for all traits, although the number of genotypes lacking haplotype StAOS2_A 691 C 692 was small. By two-way ANOVA, StAOS2_SNP691/692 was found to be a significantly contributing factor. Haplotype StAOS2 A 691 C 692 was more resistant to late blight and later maturing. The contribution of D-type cytoplasm was not detected. However, the combination of D-type cytoplasm with haplotype StAOS2_A 691 C 692 showed significantly higher late blight resistance (MCR and rAUDPC) than StAOS2_A 691 C 692 alone (P < 0.001), whereas no difference was observed for plant maturity (P = 0.42).

Table 7 Interaction for effects of combinations of presence (+) and absence (−) of nuclear markers with presence (+) or absence (−) of D-type cytoplasm on resistance to late blight and plant maturity in the PIN184 population. In case significant differences were found among marker classes, two-way ANOVA was performed to explore the factors contributing the differences

Discussion

Cytoplasmic diversity in European potatoes

We found that T (59.4 %), D (27.4 %) and W/γ (10.4 %) were the major cytoplasm types in 694 varieties and 523 breeding clones of European potatoes. Lössl et al. [30] analyzed 144 German varieties and 140 di-haploid breeding clones and found plastid-mitochondrial types- T/β (corresponding to T-type cytoplasm) in 47 %, W/α (corresponding to D-type cytoplasm) in 40 %, and W/γ-type cytoplasm in 10 % of the analyzed genotypes. Thus, our result when using a much larger number of genotypes supports the finding of Lössl et al. [30] that T-, D-, and W/γ-type cytoplasm in this order, were the most prevalent cytoplasm types in European potatoes. However the frequencies of the respective cytoplasm types differed considerably between the two studies, as well as between varieties and breeding clones and between populations analyzed in our study (Table 2).

Previously, a collection of 488 Japanese potato germplasm including 84 varieties, 378 breeding clones and 26 landraces was investigated for the cytoplasm types. T, D, P, A, M and W types were found with frequencies of 72.1 %, 17.8 %, 6.4 %, 1.2 %, 0.2 % and 2.3 %, respectively [32]. The Japanese collection seems essentially similar to European potatoes in the sense that T-type cytoplasm was the most prevalent. However, the frequencies of D- and W-type cytoplasm were much lower in the Japanese collection compared to European potatoes.

The T-type cytoplasm is understandably predominant in European and Japanese potatoes because most varieties maternally descended from ‘Rough Purple Chili’ and a few other clones from ssp. tuberosum [6, 47]. The differences in the frequencies of D- and W-type cytoplasm likely result from a more extensive use of S. demissum-derived late blight resistance and S. stoloniferum-derived PVY resistance in German breeding programs [5, 30, 32]. As chemical control has become a standard practice over the last few decades, late blight resistance breeding had lower priority in Japanese potato breeding [58]. In addition, the S. chacoense-derived PVY resistance gene Ry chc has been used in Japanese breeding programs instead of S. stoloniferum-derived PVY resistance genes [35, 59]. For these reasons, the frequencies of D- and W/γ-type cytoplasm in Japanese potato germplasm are still lower compared to European potatoes. In contrast, the frequencies of D- and W-type cytoplasm were also much higher (38 % and 11 %, respectively, [60]) in the breeding program of the International Potato Center (CIP = Centro International de la Papa), because CIP aims to deliver pest- and disease-resistant varieties for developing countries where chemical control is not practical.

It is known that clones with the D- or W/γ-type cytoplasm are functionally male sterile [30, 41, 42, 4446], so that these clones were used only as female parents, resulting in the infiltration of the common potato gene pool by these cytoplasm types [32]. The comparison of European, Japanese and Latin American gene pools demonstrates that our gene pools are being infiltrated by male sterility accompanying the D- or W/γ-type cytoplasm. The increasing frequencies of D- and W/γ- type cytoplasm enlarge the problem in designing successful mating combinations because the choice of male parents will be strictly limited, as warned previously by Provan et al. [47] and Hosaka and Sanetomo [32]. Another S. stoloniferum-derived PVY resistance gene Ry-f sto could be used instead of Ry sto because it is delivered through male fertile clones [61, 62]. Some genotypes with D-type cytoplasm also have been empirically known to be male fertile. Once the target gene(s) is transferred to genotypes with a cytoplasmic genome other than D- or W/γ-type cytoplasm, breeders could be liberated from tedious crossing activities associated with male sterility. Alternatively, a fertility-restoring gene, such as the Rt gene, which partially circumvents male sterility caused by the nuclear and T-type cytoplasm interactions [63], can be searched among genotypes with the D- or W/γ-type cytoplasm.

A prominent difference between European and Japanese potatoes was the presence of P-type cytoplasm (6.4 %) in the Japanese collection, while none was found in European potatoes. In Japan, the first diploid variety “Inca-no-mezame” was released in 2001, which has P-type cytoplasm derived from S. phureja [64]. It has an excellent taste, although it produces smaller tubers and lower yield than tetraploid cultivars. A chromosome-doubled, tetraploid clone with excellent taste was formed from “Inca-no-mezame” and combined with cyst nematode and PVY resistance genes, resulting in a breeding line “Saikai 35” [35]. Since good male fertility was recognized in genotypes with P-type cytoplasm, Saikai 35 has been extensively used to create new series of breeding clones for double cropping in Southern Japan [35]. This is why the P-type cytoplasm is increasingly present in Japanese breeding programs.

Cytoplasmic origin of European potatoes

The T-type cytoplasm was most likely derived from Chilean S. tubersoum ssp. tuberoum via a series of selfed lines “Rough Purple Chili”, “Garnet Chili”, and “Early Rose”, and a few other clones [6, 10]. All D-type cytoplasm was derived from S. demissum, whilst W-type cytoplasm was designated to a diverse group of wild potato species [32]. Within W-type cytoplasm, α-, β-, and γ-type mitochondrial DNA can be distinguished using one PCR marker (ALM_4/ALM_5, [30]). Although the W/γ subtype was assigned to genotypes with S. stoloniferum-derived cytoplasm [30], the W/γ subtype was not only found in S. stoloniferum but also in S. chacoense Bitt., S. pampasense Hawkes, S. pinnatisectum Dun., and S. vernei Bitt. et Wittm. [32]. Fifty-two of 720 publicly available genotypes (Additional file 1: Table S1) had W/γ-type cytoplasm. According to the literature [5, 30, 43] at least 19 varieties and one MPI breeding clone had either the S. stoloniferum-derived Ry sto gene or S. stoloniferum cytoplasm. Moreover, at least 21 of 25 breeding clones with W/γ-type cytoplasm in the PIN184 population were male sterile (H.-R. Hofferbert and E. Tacke, personal communication). Thus, most genotypes with W/γ-type cytoplasm were likely derived from S. stoloniferum. S. stoloniferum is a Mexican tetraploid species and is highly polymorphic [65], in which at least W/α-, W/γ-, and D-type cytoplasm were detected [30, 66, 67]. Nevertheless, varieties having the Ry sto gene have exclusively W/γ-type cytoplasm and their maternal lineages could be traced back to three S. stoloniferum accessions from which parental lines MPI 13128, MPI 46.152/1 and 43/60/96/1 were derived [43]. It is unknown whether the W/γ-type cytoplasm interacts in some way with the nuclear gene Ry sto for expression of resistance to PVY. Genotypes with another S. stoloniferum-derived PVY resistance gene Ry-f sto are male fertile [61] and have D-type cytoplasm (unpublished data). Thus, it is a further question whether tetrad-sterility is caused by the W/γ-type cytoplasm of S. stoloniferum or by a few S. stoloniferum accessions introduced to the aforementioned parental lines of a source of Ry sto .

In the past, German potato breeding used S. vernei as a source for resistance to Globodera rostochiensis [5, 68, 69], implying that some genotypes may have W/γ-type cytoplasm derived from S. vernei. According to pedigree analysis, sixteen varieties might retain S. vernei germplasm. However, twelve of those varieties had T-type cytoplasm (cvs. Amethyst, Amigo, Aula, Benol, Compagnon, Culpa, Danva, Darwina, Hydra, Krostar, Puntila and Valetta), while the remaining four had D-type cytoplasm (cvs. Astarte, Mara, Nordlicht and Roxy) (Additional file 1: Table S1). This suggests that S. vernei was used as a pollen parent in the initial cross or its hybrid progeny of a later generation.

W/αβ-type cytoplasm was found in one variety “Rita”. This type has been found in S. chacoense Bitt. [32] and S. spegazzinii (Hosaka K, personal communication). S. spegazzinii was used in the past as a source for resistance to pathotypes Ro1 and Ro5 of G. rostochiensis [5]. Thus, “Rita” might have S. spegazzinii-derived cytoplasm.

Nine genotypes had the A-type cytoplasm, which was possibly derived from S. tuberosum ssp. andigena. Four genotypes had the M-type cytoplasm (Table 2). Considering German breeding history, S. acaule might be the most probable source for this cytoplasm, although pedigrees for these breeding clones are not available. As W/α- and W/β-type cytoplasm was found in various wild species [32], we could not clarify the cytoplasmic origin for the few genotypes with W/α- or W/β-type cytoplasm.

Correlation of different cytoplasm types with agronomic traits

To the best of our knowledge, this study is the first that uses large populations to demonstrate a correlation of distinct cytoplasmic genomes with complex agronomic traits in potato. We found significant effects of cytoplasm type on various agronomic traits such as resistance to late blight and tuber bruising, plant maturity, tuber shape, starch content and yield (Table 3). No cytoplasmic difference was found for processing quality traits such as chip color and reducing sugar content (Table 3). In classical studies using reciprocal hybrids, it has been shown that the cytoplasmic genome of S. tuberosum ssp. tuberosum (mostly T-type cytoplasm) is characterized by higher percentage of tuberization, higher tuber yield, higher tuber numbers, and earlier vine maturity compared with that of S. tuberosum ssp. andigena (mostly A-type cytoplasm) [2428]. Unfortunately, in the present study, genotypes with A-type cytoplasm were too rare for statistically significant comparisons between T- and A-type cytoplasm. Nevertheless, earlier plant maturity of genotypes with T-type cytoplasm compared to other cytoplasm types was observed in three independent populations (BRUISE, GBC, and PIN184), although statistical support was obtained only in the BRUISE and PIN184 populations (Table 3). A small effect on tuber yield was only observed in the BRUISE population.

Correlation of W/γ-type cytoplasm with tuber starch content

W/γ-type cytoplasm was strongly correlated with increased tuber starch content in all three populations evaluated for TSC (Table 3). Lössl et al. [30] ranked tuber starch content according to W/γ (19 varieties) ≥ D (46 varieties) > T (79 varieties), which is in good agreement with our findings. Starch biosynthesis and degradation takes place in chloroplasts and amyloplasts and is controlled by a large number of nuclear genes. Eighteen quantitative traits loci (QTLs) for tuber starch-content were identified on all 12 potato linkage groups [70]. Efficiency of photosynthesis and carbon flux from source to sink tissues are important for the accumulation of starch in the tuber amyloplasts [71]. It is therefore conceivable that not only nuclear but also cytoplasmic factors play a functional role in starch accumulation and degradation and thereby cytoplasm type influences tuber starch content [30]. Higher tuber starch content is correlated with increased susceptibility to tuber bruising, later plant maturity and more circular tuber shape [52]. Thus, the effects of W/γ-type cytoplasm on the traits BI, PM and TS can be explained by correlation with TSC (Table 3).

Correlation of cytoplasm type with late blight resistance

Data on late blight resistance were available for the GBC and PIN184 populations [55, 56]. The two populations were rather different. The GBC population consisted of historical varieties and breeding materials that were accompanied by resistance scores obtained in various environments. The PIN184 population represented contemporary germplasm of commercial breeding programs in Germany, and all genotypes were field evaluated under similar conditions. Nevertheless reproducible effects of cytoplasm type on foliage resistance to late blight were observed in both populations. Genotypes with D-type cytoplasm had on average higher levels of foliage resistance to late blight than genotypes with T-type cytoplasm. Genotypes with M-, W/β-, or W/γ-type cytoplasm also tended to have higher foliage resistance compared to those with T-type cytoplasm (Table 3). The four genotypes with M-type cytoplasm had the highest resistance level.

For the GBC population, Gebhardt et al. [55] reported that presence of the allele-specific markers CosA210, and R11400 was associated with increased resistance of both foliage and tubers to late blight, whereas presence of R11800 was associated with increased susceptibility. CosA210 and R11400, but not R11800, were also slightly associated with later maturity. The effect of D-type cytoplasm on late blight resistance was independent from the presence of the R1 resistance gene on Chromosome V because none of the nuclear R1 markers was correlated with D-type cytoplasm (Table 4). The combination of D-type cytoplasm with the nuclear R1 markers showed that the effects of D-type cytoplasm was dominant (Table 6). This is probably because S. demissum-derived R genes other than R1 were also incorporated together with D-type cytoplasm in the GBC population. D-type cytoplasm did not affect tuber resistance to late blight (Table 6). Consequently, D-type cytoplasm was not a contributing factor to tuber resistance in the marker combinations.

In the PIN184 population, neither R1 nor R3 markers showed association with late blight resistance [56]. Accordingly, in combinations of D-type cytoplasm with these markers only the D-type cytoplasm contributed to the observed significant effects on late blight resistance, confirming the independence of the effect of D-type cytoplasm from the R genes. The late blight resistance genes R3a and R3b are tightly linked with StAOS2 (allene oxide synthase 2) on chromosome XI [56]. StAOS2 was identified as a major locus for quantitative resistance to late blight in the PIN184 population and might also have a functional role in quantitative resistance to black leg/tuber soft rot caused by bacterium Erwinia carotovora ssp. atroseptica [56, 72, 73]. The frequencies of the SNP haplotype StAOS2_A 691 C 692 which is associated with increased resistance were significantly higher in genotypes with M- and D-type cytoplasm compared to those with T-type cytoplasm (Table 5). In combination with D-type cytoplasm, the SNP haplotype StAOS2_A 691 C 692 dominated the effect of the D-type cytoplasm. However, genotypes having both D-type cytoplasm and haplotype StAOS2_A 691 C 692 had the highest level of resistance to late blight compared to the other three marker classes (Table 7), suggesting that the combination of both markers might improve field resistance to late blight more than each marker alone. The StAOS2 gene encodes a key enzyme in the biosynthesis of the defense signaling molecule jasmonic acid [74]. Plant AOS enzymes are localized in chloroplasts [7578]. Cloned StAOS2 alleles included chloroplast targeting signal peptides and were localized in potato chloroplasts [73, 79]. In general, processes in the chloroplast seem to play a role in quantitative resistance to late blight, because many nuclear genes operating in the chloroplast show constitutive higher expression levels in genotypes with quantitative resistance when compared to more susceptible genotypes [80]. Cytoplasm type might interfere with these processes either positively or negatively, thus providing as explanation for the effects of cytoplasm type on resistance.

Conclusions

Until recently, potato cytoplasmic genomes could be distinguished only by phenotypic comparisons between reciprocal cross hybrids. Molecular markers that could be used practically distinguished only between S. tuberosum ssp. tuberosum-type (T-type) and the other types [7, 9, 81, 82]. The new PCR-based technique distinguishes several potato cytoplasm types, and is simple, inexpensive and rapid, thereby facilitating large-scale cytoplasmic surveys in a short period of time [32]. Applying this technique to a large collection of European potatoes, we show that the use of M-, A- and P-type cytoplasm of Andean primitive cultivated potatoes was very limited in European potato breeding. Instead, the male-sterility-inducing D- and W/γ-type cytoplasm have been extensively used. The genetic basis of the cultivated potato could be broadened by introgression of cultivars with M- A- and P-type cytoplasm.

Thanks to the fact that phenotypic data for several important agronomic traits were available for the majority of the collection genotyped in this study for cytoplasm type, it was possible to assess for the first time whether cytoplasm type affects agronomic performance. D- and W/γ-type cytoplasm showed positive correlations with tuber starch content and foliage resistance to late blight. We identified four breeding clones with high resistance to late blight, which had the M-type cytoplasm. The M-type cytoplasm is a major cytoplasm among ancestral wild species of cultivated potatoes [34], but is rarely found in modern varieties [32]. The M-type cytoplasmic genome might have a direct function in resistance or is just co-inherited with other, unknown nuclear resistance factors. Identifying the cytoplasm types and verifying their relationships with phenotypic differences will contribute to the understanding of cytoplasmic genome diversity and function. Cytoplasm type is a novel DNA-based marker that can be used in combination with nuclear markers, for more efficient germplasm enhancement in potato.

Methods

Plant materials and phenotypic data

A total of 1,383 tetraploid genotypes consisting of 856 varieties and 527 breeding clones were determined for the cytoplasm type (Tables 1, 2 and Additional file 1: Table S1). These consisted of six populations, five of which have been used for association mapping: 205 genotypes selected for variation of susceptibility to tuber bruising [52] here referred to as the ‘BRUISE’ population, 228 genotypes selected for variation of chip quality [50, 51] referred to as the ‘CHIPS-ALL’ population, a gene bank collection of 536 genotypes with passport data for late blight resistance and plant maturity [55] referred to as the ‘GBC’ population, 40 genotypes selected for high and low processing quality [53] referred to as the ‘SUGAR40’ population, and 184 genotypes selected for variation of late blight resistance [56] referred to as the ‘PIN184’ population. The sixth population of 190 genotypes was named ‘EURO-CUL’ (European cultivars), included 140 varieties grown between 1987 and 1994 in the field at Scharnhorst (outstation of the MPI for Plant Breeding Research), 121 of which were used by Görg et al. [83] for variety identification, 10 varieties collected in the summer 2013 at the demonstration garden of the Max Planck Institute for Plant Breeding Research at Cologne, 29 varieties and three breeding clones received between 2001 and 2013 from the IPK Genebank, Germany, and the remaining 8 varieties obtained from various sources.

Evaluation data for several quantitative agronomic traits were available for these populations except EURO-CUL (Table 1). These phenotypic data were used in this study to evaluate effects of different cytoplasmic DNA types. The BRUISE population has been phenotyped for susceptibility to black spot bruising upon mechanical impact, measured as bruising index (BI) and starch corrected bruising (SCB), plant maturity (PM), tuber shape (TS), tuber starch content (TSC) and tuber yield (TY) [52]. The CHIPS-ALL population has been evaluated for chip quality after harvest (CQA) and after 3 months storage at 4 °C (CQS), TSC, tuber yield (TY) and starch yield (TSY = TSC × TY) [51]. Scores between 1 and 9 for plant maturity (PM), foliage resistance to late blight (RLBF) and tuber resistance to late blight (RLBT) were available for the GBC population [55]. This population has been genotyped with four nuclear DNA markers selected based on linkage to QTL for resistance to late blight and plant maturity for association analysis [55]. Markers CosA, and GP179, and the R1 gene for resistance to late blight [8486] are located within 1 cM on potato chromosome V in a “hot spot” for pathogen resistance [87]. Phenotypic traits available for the PIN184 population were the relative area under disease progress curve (rAUDPC) for field infestation of Phytophthora infestans, plant maturity (PM) and maturity corrected resistance (MCR) [56]. In addition, the PIN184 population has been evaluated for tuber starch content (TSC) by measuring the specific gravity. Adjusted entry means for TSC across environments were estimated as described [56]. The PIN184 population has been genotyped among others with the nuclear marker CosA, two allele specific markers derived from the R3a and R3b late blight resistance genes and SNPs at position 691 and 692 of the StAOS2 gene, which were associated with field resistance to late blight [56]. Tuber reducing sugar content (RSC) before and after cold storage has been determined in the SUGAR40 panel [53]. Note that although plant maturity (PM) has always been scored by a 1–9 scale, in the BRUISE population PM was scored from very early (score 1) to very late (score 9), while in the GBC and PIN184 populations PM was scored from very late (score 1) to very early (score 9). Where phenotypic data are based on trials in several environments (BRUISE, CHIPS-ALL, PIN184), the adjusted means were used in the present study.

Determination of cytoplasm types

Genomic DNA was extracted from 10 to 30 mg of freeze-dried leaf tissue from each EURO-CUL genotype using MagAttract 96 DNA Plant kit (Qiagen, Hilden, Germany) and the supplier’s protocol. Cytoplasmic markers T, S, SAC, A and D for cytoplasm type determination and marker detection procedures were performed as described in Hosaka and Sanetomo [32] with the following modifications. The PCR reaction was performed in 5 μl total volume, consisting of 1 μl template DNA (approximately 20 ng/μl), 0.2 mM dNTPs, 2 mM MgCl2, 1 U of Taq DNA polymerase (PEQ LAB), and 0.5 μl each of 10× forward primer mix and 10× reverse primer mix (6 μM for T, S, and SAC markers and 10 μM for A and D markers). PCR thermal conditions were as follows: 3 min at 94 °C, followed by 35 cycles of 30 s at 94 °C, 30 s at 60 °C, and 1 min and 45 s at 72 °C, and terminated with one cycle of 7 min at 72 °C. After the PCR reaction, the samples were mixed with 5 μl digestion mix including 1 μl 10× FastDigest buffer (Fermentas) and 5 U FastDigest BamHI (Fermentas). Restriction digestions were performed at 37 °C in the thermal cycler for at least 10 min. After BamHI digestion, samples with 10 μl volume were mixed with 3 μl Orange dye (Sigma) and loaded on a 3 % agarose gel in 1× TBE buffer (89 mM Tris–borate, 89 mM boric acid, and 2 mM EDTA). mtDNA types α, β, and γ were distinguished using the ALM_4/ALM_5 primers as described in [30]. 0.5 μl of 10× primer mix in the PCR reaction were replaced by 0.5 μl each of 10 μM ALM_4 and ALM_5 primers. According to Lössl et al. [30], presence of a single 2.4 kb band or a single 1.6 kb band corresponds to α-type or β-type mtDNA, respectively, while band absence indicates γ-type mtDNA.

Statistic analysis

Effects of different cytoplasm types on parametric agronomic traits (TY, TSC, TSY, SCB, rAUDPC and MCR) were analyzed by a one-way analysis of variance (one-way ANOVA) test. If significant differences were found at the 5 % level, cytoplasm types were compared using Tukey’s test. Nonparametric traits showing the difference of the group variances by Bartlett test (PM, TS, CQA, CQS, PM, RLBE, RLBT, RSC and BI) were analyzed using Welch’s test. If significant differences were found at 5 % level, each pair of cytoplasm types was compared by a Kruskal-Wallis test. Note that not all genotypes determined for cytoplasm types had phenotypic data, so that the number of genotypes (n) in the populations differed between traits (Table 3).

Nuclear markers scored in the GBC and PIN184 populations were analyzed for correlation with the presence or absence of the marker for D-type cytoplasm (Table 4). For the SNP marker StAOS2_SNP691/692, a combination of adenine at position 691 and cytosine at position 692 was regarded as the StAOS2_A 691 C 692 haplotype. Allele frequency was calculated by totalizing allele dosage of StAOS2_A 691 C 692 over all genotypes in a group and then dividing by the total number of chromosomes (=number of genotypes in the group × 4). The allele frequencies of different cytoplasm types were compared by a Welch’s analysis. Each pair was compared using Kruskal-Wallis test (Table 5). Combined effects of each nuclear marker and the D-marker on late blight resistance and plant maturity were analyzed. Four marker classes +/+, +/−, −/+, and −/− were tested by Welch’s test for significant differences, where ‘plus’ indicates the presence and ‘minus’ the absence of the marker, irrespective of allele dosage. If significant differences were found among marker classes, two-way ANOVA with the nuclear marker and the D-type cytoplasm as different factors was performed. All statistical analyses were conducted using JMP software (SAS Institute, Inc.).

Availability of supporting data

The data sets supporting the results of this article are included within the article and its additional files, Additional file 1 and Additional file 2.

Abbreviations

PVY:

Potato virus Y

PVX:

Potato virus X

PCR:

Polymerase chain reaction

QTL:

Quantitative trait locus

SNPs:

Single nucleotide polymorphisms

StAOS2 :

Allene Oxide Synthase 2

mtDNA:

Mitochondrial DNA

BI:

Bruising index

SCB:

Starch corrected bruising

PM:

Plant maturity

TS:

Tuber shape

TSC:

Tuber starch content

TY:

Tuber yield

CQA:

Chip quality after harvest

CQS:

Chip quality after 3 months storage at 4 °C

RSC:

Tuber reducing sugar content

RLBF:

Foliage resistance to late blight

RLBT:

Tuber resistance to late blight

MCR:

Maturity corrected resistance

rAUDPC:

Relative area under disease progress curve

CIP:

Centro International de la Papa (International Potato Center)

References

  1. Chase CD. Cytoplasmic male sterility: a window to the world of plant mitochondrial-nuclear interactions. Trends Genet. 2007;23:81–90.

    Article  CAS  PubMed  Google Scholar 

  2. Palmer JD, Jansen RK, Michaels HJ, Chase MW, Manhart JR. Chloroplast DNA variation and plant phylogeny. Ann Missouri Bot Gard. 1988;75:1180–206.

    Article  Google Scholar 

  3. Birky Jr CW. Uniparental inheritance of mitochondrial and chloroplast genes: mechanisms and evolution. Proc Natl Acad Sci U S A. 1995;92:11331–38.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  4. Tsunewaki K. Genetic Diversity of the Cytoplasm in Triticum and Aegilops. Tokyo: Japan Soc Promot Sci; 1980.

    Google Scholar 

  5. Ross H. Potato Breeding - Problems and Perspectives. Paul Parey: Berlin and Hamburg; 1986.

    Google Scholar 

  6. Plaisted RL, Hoopes RW. The past record and future prospects for the use of exotic potato germplasm. Am Potato J. 1989;66:603–27.

    Article  Google Scholar 

  7. Spooner DM, Núñez J, Trujillo G, Herrera MDR, Guzmán F, Ghislain M. Extensive simple sequence repeat genotyping of potato landraces supports a major reevaluation of their gene pool structure and classification. Proc Natl Acad Sci U S A. 2007;104:19398–403.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  8. Hawkes JG. The Potato: Evolution, Biodiversity and Genetic Resources. London: Belhaven Press; 1990.

    Google Scholar 

  9. Ames M, Spooner DM. DNA from herbarium specimens settles a controversy about origins of the European potato. Am J Bot. 2008;95:252–7.

    Article  CAS  PubMed  Google Scholar 

  10. Mendosa HA, Haynes FL. Genetic relationships among potato cultivars grown in the United States. HortSci. 1974;9:328–30.

    Google Scholar 

  11. Salaman RN. Breeding for immunity to late blight and other diseases in the potato. Edinburgh: Proc 7th Int Genet Congr 1939; 1941. p. 253–4.

    Google Scholar 

  12. Glendinning DR. Potato introductions and breeding up to the early 20th century. New Phytol. 1983;94:479–505.

    Article  Google Scholar 

  13. Cockerham G. Potato breeding for virus resistance. Ann Appl Biol. 1943;30:105–8.

    Article  Google Scholar 

  14. Ross H. Inheritance of extreme resistance to virus Y in Solanum stoloniferum and its hybrids with Solanum tuberosum. In: Proc Third Conf Potato Virus Diseases. 1958. p. 204–11.

    Google Scholar 

  15. Cockerham G. Genetical studies on resistance to potato viruses X and Y. Heredity. 1970;25:309–48.

    Article  Google Scholar 

  16. Ross H. Die Vererbung der “Immunität” gegen das X-Virus in tetraploidem Solanum acaule. Carylogia. 1954;6(Suppl):1128–32.

    Google Scholar 

  17. Ellenby C. Resistance to the potato-root eelworm. Nature. 1948;162:704.

    Article  Google Scholar 

  18. Ellenby C. Tuber forming species and varieties of the genus Solanum tested for resistance to the potato root eelworm Heterodera rostochiensis Wollenweber. Euphytica. 1954;3:195–202.

    Article  Google Scholar 

  19. Huijsman CA. Breeding for resistance to the potato root eelworm. II. Data on the inheritance of resistance in andigenum-tuberosum crosses obtained in 1954. Euphytica. 1955;4:133–40.

    Article  Google Scholar 

  20. Frandsen NO. Resistenzzüchtung gegen pilzliche und bakterielle Krankheiten der Kartoffel. In: Kappert H, Rudorf H, editors. Handbuch der Pflanzenzüchtung. Berlin: Paul Parey; 1958. p. 71–97.

    Google Scholar 

  21. van Soest LJM, Seidewitz L. Evaluation data on tuber-bearing Solanum species. In Inst f Pflanzenbau u Pflanzenzüchtg der FAL Braunschweig u Stichting voor Plantenveredeling. Wageningen; 1981.

  22. Schilberszky K. A new parasite causing potato wart disease (in German). Ber Dtsch Bot Ges. 1896;14:36–7.

    Google Scholar 

  23. Black W. Studies on the inheritance of resistance to wart disease (Synchytrium endobioticum (Schilb.) Perc.) in potatoes. J Genet. 1935;30:127–46.

    Article  Google Scholar 

  24. Sanford JC, Hanneman Jr RE. Reciprocal differences in the photoperiod reaction of hybrid populations in Solanum tuberosum. Am Potato J. 1979;56:531–40.

    Article  Google Scholar 

  25. Hoopes RW, Plaisted RL, Cubillos AG. Yield and fertility of reciprocal-cross Tuberosum-Andigena hybrids. Am Potato J. 1980;57:275–84.

    Article  Google Scholar 

  26. Sanford JC, Hanneman Jr RE. Large yield differences between reciprocal families of Solanum tuberosum. Euphytica. 1982;31:1–12.

    Article  Google Scholar 

  27. Hilali A, Lauer FI, Veilleux RE. Reciprocal differences between hybrids of Solanum tuberosum Groups Tuberosum (haploid) and Phureja. Euphytica. 1987;36:631–9.

    Article  Google Scholar 

  28. Maris B. Analysis of an incomplete diallel cross among three ssp. tuberosum varieties and seven long-day adapted ssp. andigena clones of the potato (Solanum tuberosum L.). Euphytica. 1989;41:163–82.

    Article  Google Scholar 

  29. Hosaka K. Who is the mother of the potato? Restriction endonuclease analysis of chloroplast DNA of cultivated potatoes. Theor Appl Genet. 1986;72:606–18.

    Article  CAS  PubMed  Google Scholar 

  30. Lössl A, Götz M, Braun A, Wenzel G. Molecular markers for cytoplasm in potato: male sterility and contribution of different plastid-mitochondrial configurations to starch production. Euphytica. 2000;116:221–30.

    Article  Google Scholar 

  31. Hosaka K. Distribution of the 241 bp deletion of chloroplast DNA in wild potato species. Am J Potato Res. 2002;79:119–23.

    Article  CAS  Google Scholar 

  32. Hosaka K, Sanetomo R. Development of a rapid identification method for potato cytoplasm and its use for evaluating Japanese collections. Theor Appl Genet. 2012;125:1237–51.

    Article  CAS  PubMed  Google Scholar 

  33. Sukhotu T, Kamijima O, Hosaka K. Nuclear and chloroplast DNA differentiation in Andean potatoes. Genome. 2004;47:46–56.

    Article  CAS  PubMed  Google Scholar 

  34. Hosaka K, Sanetomo R. Comparative differentiation in mitochondrial and chloroplast DNA among cultivated potatoes and closely related wild species. Genes Genet Syst. 2009;84:371–8.

    Article  CAS  PubMed  Google Scholar 

  35. Mori K, Mukojima N, Nakao T, Tamiya S, Sakamoto Y, Sohbaru N, et al. Germplasm release: Saikai 35, a male and female fertile breeding line carrying Solanum phureja-derived cytoplasm and potato cyst nematode resistance (H1) and Potato virus Y resistance (Ry chc ) genes. Am J Potato Res. 2012;89:63–72.

    Article  Google Scholar 

  36. Sanetomo R, Hosaka K. A maternally inherited DNA marker, descended from Solanum demissum (2n = 6x = 72) to S. tuberosum (2n = 4x = 48). Breed Sci. 2011;61:426–34.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  37. Hosaka K, Sanetomo R. Application of a PCR-based cytoplasm genotyping method for phylogenetic analysis in potato. Am J Potato Res. 2014;91:246–53.

    Article  Google Scholar 

  38. Hanson MR, Bentolila S. Interactions of mitochondrial and nuclear genes that affect male gametophyte development. Plant Cell. 2004;16:S154–69.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  39. Sofi PA, Rather AG, Wani SA. Genetic and molecular basis of cytoplasmic male sterility in maize. Commun Biometry Crop Sci. 2007;2:49–60.

    Google Scholar 

  40. Grun P, Ochoa C, Capage D. Evolution of cytoplasmic factors in tetraploid cultivated potatoes (Solanaceae). Am J Bot. 1977;64:412–20.

    Article  Google Scholar 

  41. Brown CR. Tetrad sterility: a cytoplasmic-genic male sterility attractive to bumblebees. Interlaken: Proc 9th Trien conf Eur Assn Potato Res; 1984. p. 101–2.

    Google Scholar 

  42. Ortiz R, Iwanaga M, Peloquin SJ. Male sterility and 2n pollen in 4x progenies derived from 4x × 2x and 4x × 4x crosses in potatoes. Potato Res. 1993;36:227–36.

    Article  Google Scholar 

  43. Song YS, Schwarzfischer A. Development of STS markers for selection of extreme resistance (Ry sto ) to PVY and maternal pedigree analysis of extremely resistant cultivars. Am J Potato Res. 2008;85:159–70.

    Article  CAS  Google Scholar 

  44. Abdalla MMF, Hermsen JGT. The plasmon-genic basis of pollen lobedness and tetrad sterility in Solanum verrucosum hybrids and duplicate linkage groups. Genetica. 1971;42:261–70.

    Article  Google Scholar 

  45. Grun P, Aubertin M, Radlow A. Multiple differentiation of plasmons of diploid species of Solanum. Genetics. 1962;47:1321–33.

    CAS  PubMed Central  PubMed  Google Scholar 

  46. Dionne LA. Cytoplasmic sterility in derivatives of Solanum demissum. Am Potato J. 1961;38:117–20.

    Article  Google Scholar 

  47. Provan J, Powell W, Dewar H, Bryan G, Machray GC, Waugh R. An extreme cytoplasmic bottleneck in the modern European cultivated potato (Solanum tuberosum) is not reflected in decreased levels of nuclear diversity. Proc R Soc Lond B. 1999;266:633–9.

    Article  Google Scholar 

  48. Ortega F, Lopez-Vizcon C. Application of molecular marker-assisted selection (MAS) for disease resistance in a practical potato breeding programme. Potato Res. 2012;55:1–13.

    Article  CAS  Google Scholar 

  49. Gebhardt C. Bridging the gap between genome analysis and precision breeding in potato. Trends Genet. 2013;29:248–56.

    Article  CAS  PubMed  Google Scholar 

  50. Li L, Tacke E, Hofferbert HR, Lübeck J, Strahwald J, Draffehn AM, et al. Validation of candidate gene markers for marker-assisted selection of potato cultivars with improved tuber quality. Theor Appl Genet. 2013;126:1039–52.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  51. Li L, Paulo MJ, Strahwald J, Lübeck J, Hofferbert HR, Tacke E, et al. Natural DNA variation at candidate loci is associated with potato chip color, tuber starch content, yield and starch yield. Theor Appl Genet. 2008;116:1167–81.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  52. Urbany C, Stich B, Schmidt L, Simon L, Berding H, Junghans H, et al. Association genetics in Solanum tuberosum provides new insights into potato tuber bruising and enzymatic tissue discoloration. BMC Genomics. 2011;12:7.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  53. Fischer M, Schreiber L, Colby T, Kuckenberg M, Tacke E, Hofferbert HR, et al. Novel candidate genes influencing natural variation in potato tuber cold sweetening identified by comparative proteomics and association mapping. BMC Plant Biol. 2013;13:113.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  54. Schreiber L, Nader-Nieto AC, Schönhals EM, Walkemeier B, Gebhardt C. SNPs in genes functional in starch-sugar interconversion associate with natural variation of tuber starch and sugar content of potato (Solanum tuberosum L.). G3. 2014;4:1797–811.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  55. Gebhardt C, Ballvora A, Walkemeier B, Oberhagemann P, Schüler K. Assessing genetic potential in germplasm collections of crop plants by marker-trait association: a case study for potatoes with quantitative variation of resistance to late blight and maturity type. Mol Breed. 2004;13:93–102.

    Article  CAS  Google Scholar 

  56. Pajerowska-Mukhtar KM, Stich B, Achenbach U, Ballvora A, Lübeck J, Strahwald J, et al. Single nucleotide polymorphisms in the Allene Oxide Synthase 2 gene are associated with field resistance to late blight in populations of tetraploid potato cultivars. Genetics. 2009;181:1115–27.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. van Berloo R, Hutten RCB, van Eck HJ, Visser RGF. An online potato pedigree database resource. Potato Res. 2007;50:45–57.

    Article  Google Scholar 

  58. Akino S, Takemoto D, Hosaka K. Phytophthora infestans: a review of past and current studies on potato late blight. J Gen Plant Pathol. 2013;80:24–37.

    Article  Google Scholar 

  59. Hosaka K, Hosaka Y, Mori M, Maida T, Matsunaga H. Detection of a simplex RAPD marker linked to resistance to potato virus Y in a tetraploid potato. Am J Potato Res. 2001;78:191–6.

    Article  CAS  Google Scholar 

  60. Mihovilovich E, Sanetomo R, Hosaka K, Ordoñez B, Aponte M, Bonierbale M. Cytoplasmic diversity in potato breeding: case study from the International Potato Center. Mol Breed. 2015;35:137.

    Article  Google Scholar 

  61. Flis B, Hennig J, Strzelczyk-Żyta D, Gebhardt C, Marczewski W. The Ry-f sto gene from Solanum stoloniferum for extreme resistance to Potato virus Y maps to potato chromosome XII and is diagnosed by PCR marker GP122718 in PVY resistant potato cultivars. Mol Breed. 2005;15:95–101.

    Article  CAS  Google Scholar 

  62. Song YS, Hepting L, Schweizer G, Hartl L, Wenzel G, Schwarzfischer A. Mapping of extreme resistance to PVY (Ry sto ) on chromosome XII using anther-culture-derived primary dihaploid potato lines. Theor Appl Genet. 2005;111:879–87.

    Article  CAS  PubMed  Google Scholar 

  63. Iwanaga M, Ortiz R, Cipar MS, Peloquin SJ. A restorer gene for genetic-cytoplasmic male sterility in cultivated potatoes. Am Potato J. 1991;68:19–28.

    Article  Google Scholar 

  64. Mori M, Ohara-Takada A, Umemura Y, Maida T, Kimura T, Takada N, et al. Breeding of diploid potato variety “Inca no mezame” with orange in the tuber flesh color (in Japanese). Breed Res. 2009;11:53–8.

    Google Scholar 

  65. Spooner DM, van den Berg RG, Rodríguez A, Bamberg J, Hijmans RJ, Cabrera SJL. Wild potatoes (Solanum section Petota; Solanaceae) of North and Central America. In: Anderson C, editor. Systematic Botany Monographs vol 68. USA: The American Society of Plant Taxonomists; 2004. p. 1–209.

    Google Scholar 

  66. Lössl A, Adler N, Horn R, Frei U, Wenzel G. Chondriome-type characterization of potato: mt α, β, γ, δ, ε and novel plastid-mitochondrial configurations in somatic hybrids. Theor Appl Genet. 1999;99:1–10.

    Article  Google Scholar 

  67. Sanetomo R, Hosaka K. A recombination-derived mitochondrial genome retained stoichiometrically only among Solanum verrucosum Schltdl. and Mexican polyploid wild potato species. Genet Resour Crop Evol. 2013;60:2391–404.

    Article  Google Scholar 

  68. Scurrah MM, Plaisted RL, Harrison ML. Resistance to the potato nematode Heterodera rostochiensis Woll. in clones derived from Solarium vernei. Am Potato J. 1973;50:9–18.

    Article  Google Scholar 

  69. Jacobs JME, van Eck HJ, Horsman K, Arens PFP, Verkerk-Bakker B, Jacobsen E, et al. Mapping of resistance to the potato cyst nematode Globodera rostochiensis from the wild potato species Solanum vernei. Mol Breed. 1996;2:51–60.

    Article  CAS  Google Scholar 

  70. Schäfer-Pregl R, Ritter E, Concilio L, Hesselbach J, Lovatti L, Walkemeier B, et al. Analysis of quantitative trait loci (QTLs) and quantitative trait alleles (QTAs) for potato tuber yield and starch content. Theor Appl Genet. 1998;97:834–46.

    Article  Google Scholar 

  71. Frommer WB, Sonnewald U. Molecular analysis of carbon partitioning in solanaceous species. J Expl Bot. 1995;46:587–607.

    Article  CAS  Google Scholar 

  72. Zimnoch-Guzowska E, Marczewski W, Lebecka R, Flis B, Schäfer-Pregl R, Salamini F, et al. QTL analysis of new sources of resistance to Erwinia carotovora ssp. atroseptica in potato done by AFLP, RFLP, and resistance-gene-like markers. Crop Sci. 2000;40:1156–67.

    Article  CAS  Google Scholar 

  73. Pajerowska-Mukhtar KM, Mukhtar MS, Guex N, Halim VA, Rosahl S, Somssich IE, et al. Natural variation of potato Allene Oxide Synthase 2 causes differential levels of jasmonates and pathogen resistance in Arabidopsis. Planta. 2008;228:293–306.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  74. Liechti R, Farmer EE. The jasmonic pathway. Science. 2002;296:1649–50.

    Article  CAS  PubMed  Google Scholar 

  75. Vick BA, Zimmerman DC. Pathways of fatty acid hydroperoxide metabolism in spinach leaf chloroplasts. Plant Physiol. 1987;85:1073–8.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  76. Laudert D, Pfannschmidt U, Lottspeich F, Holländer-Czytko H, Weiler EW. Cloning, molecular and functional characterization of Arabidopsis thaliana allene oxide synthase (CYP 74), the first enzyme of the octadecanoid pathway to jasmonates. Plant Mol Biol. 1996;31:323–35.

    Article  CAS  PubMed  Google Scholar 

  77. Maucher H, Hause B, Feussner I, Ziegler J, Wasternack C. Allene oxide synthases of barley (Hordeum vulgare cv. Salome): tissue specific regulation in seedling development. Plant J. 2000;21:199–213.

    Article  CAS  PubMed  Google Scholar 

  78. Froehlich JE, Itoh A, Howe GA. Tomato allene oxide synthase and fatty acid hydroperoxide lyase, two cytochrome P450s involved in oxylipin metabolism, are targeted to different membranes of chloroplast envelope. Plant Physiol. 2001;125:306–17.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  79. Farmaki T, Sanmartin M, Jimenez P, Paneque M, Sanz C, Vancanneyt G, et al. Differential distribution of the lipoxygenase pathway enzymes within potato chloroplasts. J Exp Bot. 2007;58:555–68.

    Article  CAS  PubMed  Google Scholar 

  80. Draffehn AM, Li L, Krezdorn N, Ding J, Lübeck J, Strahwald J, et al. Comparative transcript profiling by SuperSAGE identifies novel candidate genes for controlling potato quantitative resistance to late blight not compromised by late maturity. Front Plant Sci. 2013. doi:10.3389/fpls.2013.00423.

  81. Gavrilenko TA, Antonova OY, Kostina LI. Study of genetic diversity in potato cultivars using PCR analysis of organelle DNA. Russ J Genet. 2007;43:1550–5.

    Article  CAS  Google Scholar 

  82. Chimote VP, Chakrabarti SK, Pattanayak D, Pandey SK, Naik PS. Molecular analysis of cytoplasm type in Indian potato varieties. Euphytica. 2008;162:69–80.

    Article  CAS  Google Scholar 

  83. Görg R, Schachtschabel U, Ritter E, Salamini F, Gebhardt C. Discrimination among 136 tetraploid potato varieties by fingerprints using highly polymorphic DNA markers. Crop Sci. 1992;32:815–9.

    Article  Google Scholar 

  84. Leonards-Schippers C, Gieffers W, Salamini F, Gebhardt C. The R1 gene conferring race-specific resistance to Phytophthora infestans in potato is located on potato chromosome V. Mol Gen Genet. 1992;233:278–83.

    Article  CAS  PubMed  Google Scholar 

  85. Meksem K, Leister D, Peleman J, Zabeau M, Salamini F, Gebhardt C. A high-resolution map of the vicinity of the R1 locus on chromosome V of potato based on RFLP and AFLP markers. Mol Gen Genet. 1995;249:74–81.

    Article  CAS  PubMed  Google Scholar 

  86. Ballvora A, Ercolano MR, Weiß J, Meksem K, Bormann C, Oberhagemann P, et al. The R1 gene for potato resistance to late blight (Phytophthora infestans) belongs to the leucine zipper/NBS/LRR class of plant resistance genes. Plant J. 2002;30:361–71.

    Article  CAS  PubMed  Google Scholar 

  87. Gebhardt C, Valkonen JPT. Organization of genes controlling disease resistance in the potato genome. Ann Rev Phytopathol. 2001;39:79–102.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

This work was carried out when R. S. stayed as a Guest Scientist for 10 months at the Max-Planck Institute for Plant Breeding Research, Cologne, Germany, with C. G. We thank Kazuyoshi Hosaka for kindly reading and checking earlier versions of the manuscript, and Hans-Reinhard Hofferbert and Eckhard Tacke for providing pollen fertility data. This study was financially supported by Calbee Inc., Hokkaido Potato Growers Association, Kewpie Corp., KENKO Mayonnaise Co., Ltd., and Japan Snack Cereal Foods Association, and by the Max-Planck Society.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Rena Sanetomo.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RS carried out the cytoplasm type determination and correlation analysis with agronomic traits, and drafted the manuscript. CG provided populations, participated in study design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Authors’ information

RS is an Assistant Professor affiliated to Potato Germplasm Enhancement Laboratory at Obihiro University of Agriculture and Veterinary Medicine. CG is the Research Group Leader of the Potato Genome Analysis Group, Department of Plant Breeding and Genetics, Max-Planck Institute for Plant Breeding Research.

Additional files

Additional file 1: Table S1.

The cytoplasm types of all European varieties and breeding clones. Sheet 1 shows the cytoplasm types of European varieties. Sheet 2 shows the cytoplasm types of European breeding clones.

Additional file 2: Table S2.

A list of duplicated varieties with different cytoplasm types. Sixteen varieties duplicated in at least two populations had different cytoplasm types.

Rights and permissions

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sanetomo, R., Gebhardt, C. Cytoplasmic genome types of European potatoes and their effects on complex agronomic traits. BMC Plant Biol 15, 162 (2015). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-015-0545-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-015-0545-y

Keywords