Open Access

Novel SNP markers in InvGE and SssI genes are associated with natural variation of sugar contents and frying color in Solanum tuberosum Group Phureja

  • Diana Duarte-Delgado1, 2,
  • Deissy Juyó1,
  • Christiane Gebhardt3,
  • Felipe Sarmiento4 and
  • Teresa Mosquera-Vásquez1Email author
BMC GeneticsBMC series – open, inclusive and trusted201718:23

DOI: 10.1186/s12863-017-0489-3

Received: 31 August 2016

Accepted: 3 March 2017

Published: 9 March 2017

Abstract

Background

Potato frying color is an agronomic trait influenced by the sugar content of tubers. The candidate gene approach was employed to elucidate the molecular basis of this trait in Solanum tuberosum Group Phureja, which is mainly diploid and represents an important genetic resource for potato breeding. The objective of this research was to identify novel genetic variants related with frying quality in loci with key functions in carbohydrate metabolism, with the purpose of discovering genetic variability useful in breeding programs. Therefore, an association analysis was implemented with 109 SNP markers identified in ten candidate genes.

Results

The analyses revealed four associations in the locus InvGE coding for an apoplastic invertase and one association in the locus SssI coding for a soluble starch synthase. The SNPs SssI-C 45711901 T and InvGE-C 2475454 T were associated with sucrose content and frying color, respectively, and were not found previously in tetraploid genotypes. The rare haplotype InvGE-A 2475187 C 2475295 A 2475344 was associated with higher fructose contents. Our study allowed a more detailed analysis of the sequence variation of exon 3 from InvGE, which was not possible in previous studies because of the high frequency of insertion-deletion polymorphisms in tetraploid potatoes.

Conclusion

The association mapping strategy using a candidate gene approach in Group Phureja allowed the identification of novel SNP markers in InvGE and SssI associated with frying color and the tuber sugar content measured by High Performance Liquid Chromatography (HPLC). These novel associations might be useful in potato breeding programs for improving quality traits and to increase crop genetic variability. The results suggest that some genes involved in the natural variation of tuber sugar content and frying color are conserved in both Phureja and tetraploid germplasm. Nevertheless, the associated variants in both types of germplasm were present in different regions of these genes. This study contributes to the understanding of the genetic architecture of tuber sugar contents and frying color at harvest in Group Phureja.

Keywords

Association mapping Potato frying quality Carbohydrate metabolism Candidate genes

Background

Potato frying color is influenced by the sugar content of the tubers, because the hydrolysis of sucrose by invertase is the main source of the reducing sugars glucose and fructose, which are precursors of the Maillard reaction. This reaction produces dark pigments and toxic products such as acrylamide during the non-enzymatic reaction between reducing sugars and free amino acids at high temperatures [1, 2]. The high content of reducing sugars in tubers causes a low quality of fried potato products. The biotic and abiotic stresses present during plant development, the temperatures during the vegetative growth, and the genotypic differences in carbohydrate metabolism are crucial factors that influence the sugar contents of tubers at harvest [3, 4]. Sucrose, glucose and fructose are important metabolic signals in transduction pathways that affect the expression of several classes of genes involved in all stages of plant development and related to stress responses. Therefore, the sugar content of storage organs such as potato tubers is a complex trait controlled by multiple genetic and environmental factors [58]. Consequently, it is of particular interest to identify variants of genes related to low tuber sugar content at harvest. Such variants might be useful diagnostic markers in breeding programs aimed at improving frying quality.

Association mapping is a strategy to study the genetic basis of complex traits that has advantages compared to linkage mapping for the discovery of diagnostic molecular markers [9]. This strategy allows the mapping of quantitative traits more precisely due to the use of natural populations or diversity panels, where there are higher numbers of ancient recombination events. Consequently, linkage disequilibrium (LD) decays faster than in a bi-parental segregating population. Therefore, the mapping resolution is increased because of the reorganization of the chromosomes in smaller regions [10, 11]. Association mapping is also advantageous because it analyses the effects of multiple alleles in a single experiment. This is important for the efficient development of diagnostic molecular markers for traits of interest in breeding populations, with the purpose of contributing to the selection of genotypes with the desired qualities [12].

There is a large body of knowledge on the metabolic pathways and enzymes involved in starch synthesis, degradation and transport in potato. Sixty nine functional genes in carbohydrate metabolism have been mapped to linkage maps of the 12 potato chromosomes [13]. Several of these candidate genes have been found to be associated with natural variation of frying color or reducing sugar content [7, 1417]. These studies have revealed the multi-loci genetic architecture of sugar content and frying color related to genes coding for functional enzymes in carbohydrate metabolism in tetraploid breeding populations of potato (2n = 4 × = 48). To assess the diagnostic power of genes with key functions in carbohydrate metabolism in diverse germplasm, it is necessary to search for genetic variants associated with frying quality in different potato populations, with the purpose of discovering genetic variability useful in potato breeding.

Solanum tuberosum Group Phureja consists mainly of diploid genotypes (2n = 2 × = 24) with short day adaptation and tuber sprouting at harvest [18]. In Andean countries, Group Phureja constitutes an important crop, particularly for small farmers [19]. Genotypes with round tubers, yellow flesh and skin represent an interesting potential for widespread use due to their desirable organoleptic properties, nutritional quality, and processing potential [2023]. This cultivated potato is also a valuable genetic resource for the introgression of agronomic traits into tetraploid potatoes from S. tuberosum subsp. tuberosum that are the most extensively grown worldwide [2427]. The doubled-monoploid genotype DM1-3516 R44 from Group Phureja represents the potato reference genome [28]. This genome sequence has contributed to the understanding of the genetic basis of various traits in tetraploid potato [7, 29]. This genomic tool is also important for studying agronomic traits of Group Phureja such as frying quality.

The phenotypic assessment of potato frying quality has been performed mostly by means of visual scales of frying color in linkage and association mapping studies [7, 15, 30]. However, it is desirable to improve the phenotypic evaluations in order to obtain precise quantitative data that contribute to the finding of reliable variants associated with quantitative traits of interest [31, 32]. This can be achieved through the use of chromatographic methods which quantify the precursors of the Maillard reaction. HPLC is capable of rapid, specific, sensitive, and precise measurements of sugar contents in potato tubers [33, 34] and should be implemented in mapping studies.

Accessions from the Colombian Central Collection of Group Phureja (CCC) that constitute the Working Collection of the Breeding Program at the National University of Colombia have been characterized for their frying quality at harvest. Plants were cultivated in plots in highlands over 2800 meters above the sea level, where low temperatures favored the accumulation of reducing sugars undesirable for frying purposes [35]. Previously, frying color was assessed in tubers from experimental fields in nine environments using a visual scale. From this multi-environmental data, a chip darkening index was estimated for each genotype using a Bayesian approach [35]. Recently, a HPLC method was applied to measure sucrose, glucose, and fructose contents after harvest in the same accessions in a single-environment trial. Reducing sugars and chip darkening index showed moderate significant correlations which supports a positive relationship between the two measurements and that both data sets are consistent [36].

The purpose of the research presented here was to find allelic variants associated with sugar contents and frying color in a Group Phureja population using a candidate gene approach. Therefore, the phenotypic assessments aforementioned were used to implement an association mapping strategy in order to elucidate the molecular basis of frying quality in Phureja tubers. The association analysis was performed with SNP markers identified in candidate genes with key function in carbohydrate metabolism, which had shown associations with potato frying color and sugar content in tetraploid potatoes [7, 1417, 37].

Methods

Plant material and phenotypes

A set of 108 diploid landrace accessions from the Working Collection of the Breeding Program at the National University of Colombia and four commercial diploid cultivars (Criolla Colombia, Criolla Latina, Criolla Galeras, and Criolla Guaneña) were used in this study. Chip color data collected in nine environments allowed the estimation of a multi-environmental chip darkening index in a previous study using a Bayesian approach [35]. The index was calculated through a scale that establishes equivalences between the International Potato Center color scale [38] and a chip darkening percentage [35]. An explanatory table of this scale with images is shown in [36]. The first degree with the lightest color represented a darkening percentage of 0, while the fifth degree with the darkest color represented a darkening percentage from 90 to 100. Tubers of the 112 genotypes were also assessed at harvest in a single environment for sucrose, glucose, fructose, reducing and total sugar contents, as well as for their glucose/fructose and sucrose/reducing sugars ratios [36]. Mean values of these eight variables (seven sugar measurements and chip darkening index) from each genotype were used for the independent association analysis of each trait.

Candidate gene sequencing and SNP calling

Genomic DNA was isolated with DNeasy Plant Mini Kit™ (Qiagen) from young leaf tissue. DNA was quantified with a NanoDrop 2000c spectrophotometer (Thermo Scientific) and adjusted to a concentration of 10 ng μL−1 for the Polymerase Chain Reaction (PCR) amplification. A group of ten candidate genes previously reported with association to the frying quality trait in tetraploid populations were selected for this study: InvGE [15, 17, 39], Pain-1 [15, 16, 39, 40], Stp23 [7, 15, 16], StpL [7, 15, 16], AGPaseS-a [7, 15], SssI [7, 16], GWD [7, 41], UGPase [42], Sps [16], and G6pdh [16]. Three primer pairs were designed per gene using the potato reference genome (version 4.03) [28, 43] in the SPUD database [44] (http://solanaceae.plantbiology.msu.edu/cgi-bin/gbrowse/potato/). In addition, seven primer pairs were tested in the Group Phureja accessions in candidate gene regions analyzed by Fischer et al. [37] (LapN and KT-InvInh) and Schreiber et al. [7] (BMY-8/2, PWD, AGPaseS, INV-8/2, and PGM-3).

The primers in 16 genes were tested in eight random genotypes (Criolla Latina, CCC 21, CCC 52, CCC 59, CCC 116, CCC 123, CCC 127, and CCC 132) and those with appropriate quality of the amplicon sequence and with at least four SNPs were selected for amplifying and sequencing in the whole population (Table 1). PCR conditions were the same as described by Schreiber et al. [7]. Annealing temperatures of primers used in the population are specified in Table 1.
Table 1

Loci (PGSC0003DMG*********) from the reference genome [28, 43] analyzed in 112 Solanum tuberosum Group Phureja accessions

Gene Acronym

Locus

Chr.

Primer Sequences 5′–3′

Ta (°C)

Amplicon

No. SNPs Scored

(GenBank Accesion No.)

size (bp)

(Scored previously)e

Stp23 (PHO1)a (D00520)

PGSC0003DMG400007782

III

*f-cagatatgttacatactctacc

59

963

4 (0)

r-tcattagtcacaactttatcgg

StpL (PHO1)a (X73684)

PGSC0003DMG400028382

V

*f-ttacattgcacaagcacaagc

57

982

14 (1) [67]

r-gtgtacatacaatactctatcc

SssI (Y10416)

PGSC0003DMG402018552

III

f-aacaataggaatttaccataacc

57

970

12 (0)

*r-atattcccaaacaaaacagagc

InvGE (INV-cw)a (AJ133765)

PGSC0003DMG400008943

IX

f-caattcttcgattcttcatagg

57

798

9 (8) [39]

*r-aattgaagcagatcatgtagg

Pain1 (INV)a (X70368)

PGSC0003DMG400013856

III

f-catacattacactatagatcc

56

926

5 (2) [39, 40]

*r-aattgaagcagatcatgtagg

   

f-caaaatgaatacatattaagaggb

56

707

4 (0)

   

*r-cttaagcagttgcttagagc

UGPase (D00667)

PGSC0003DMG401013333

XI

f-atgatgttctccacttaaaagc

56

805

11 (1) [42]

   

*r-tttcagattttcagaagagagg

   

*f-tgattaacgatactatacgtcc

56

931

12 (1) [42]

   

r-ttaaaacttccttatactatatgg

GWD (Y09533)

PGSC0003DMG400007677

V

f-ttctgttatctacttagttacg

56

992

5 (3) [7, 41]

   

*r-gttttatatcttgcttctttgg

BMY-8/2 (AF393847)

PGSC0003DMG400001855

VIII

f-gctactggacatggtgacagac

57

560

9 (9) [5]

   

*r-ttacatagaggtctgtcctgcttgag

PWD (AY747068)

PGSC0003DMG400016613

IX

*f-ggtctgatgatctatctgattgcc

57

871

16 (16) [7, 41]

   

r-gacatcttgaggagaaccaaactt

LapN (X77015)

PGSC0003DMG400007831

XII

f-gcttcctggtcttggctcd

60

985

8 (3) [37]

   

*r-gataggcatacgcagccaggtcagaaatcaa

Chr chromosome, Ta annealing temperature

*Primer used for amplicon sequencing

aThese parenthesis include the alternative acronyms used in the starch-sugar interconversion pathway scheme adapted by Schreiber et al. [7]

bFrom a region downstream of the Pain1 candidate gene

cFrom Schreiber et al. [7]

dFrom Fischer et al. [37]

eThe studies where the SNPs or the related amino acid changes were previously reported in tetraploid potatoes are indicated

Amplicons were visualized in 1.5% (w/v) agarose gels in 1× TAE (Tris-acetate-EDTA) buffer and stained with ethidium bromide (0.4 μg mL−1). PCR products were cleaned using ExoSAP-IT® (USB) and custom sequenced at the Max-Planck-Genome-Center Cologne using the dideoxy chain-termination sequencing method, an ABI PRISM Dye Terminator Cycle Sequencing Ready Reaction Kit, and an ABI PRISM 3730 automated DNA Sequencer (Applied Biosystems). SNPs were detected by visual examination of the sequence trace files for overlapping base calling peaks using Geneious Software version 8.0 (Biomatters). The SNPs were scored using Data Acquisition and Analysis Software (DAx) (Van Mierlo Software Consultancy); in the case of heterozygous individuals, overlapping base calling peaks were identified.

Association analysis

Association analysis was carried out considering that the Group Phureja accessions studied do not have a marked population structure and present a fast LD decay (r 2 value of 0.1 in a distance range of 1 to 2.5 Mbp in individual chromosomes) [27, 45]. Markers included in the analysis were those with a minor allele frequency higher than 0.01. A compressed mixed linear model (model 1) [46] and an enriched compressed mixed linear model (model 2) [47] were implemented with GAPIT R Package (version 1) [48] for the analyses. These models accounted for population structure and an additional parameter of group kinship that allows the reduction of the random effect of kinship among genotypes and improves the statistical power [46, 47]. Quantile-quantile (QQ) plots of the expected and observed F-test probabilities for the SNP markers were assessed to identify the appropriate model for controlling type I errors caused by population structure and familial relatedness [49]. Accordingly, we choose from the two models the one with better adjustment for each phenotypic variable. To improve the adjustment of the association models of glucose and fructose contents, we removed the three genotypes with extreme values which constituted a single group as described by Duarte-Delgado et al. [36]. The QQ plots from the selected models of each variable with significant associations are shown in the Additional file 1.

The probability p for marker-trait associations was expressed on a -log10 scale. SNP markers with a p-value − log10 > 2 were considered significantly associated with the variable studied. To support the reliability of the associations found with the mixed model approach, an analysis of variance (ANOVA) was performed using a linear model implemented in Genstat software version 17 (VSN International). Differences between phenotypic means of the genotypic classes were detected considering a level of significance of α = 0.05. Therefore, markers with significant p-values showing effect in the phenotypic variables were retained and presented as associated. LD was estimated between pairs of SNP markers that were found with significant associations in the same chromosome, using the r2 statistic [50]. Annotated potato gene models from the SPUD database allowed the location of associated markers in intronic or exonic regions [44]. Amino acid changes caused by SNPs in exonic regions were predicted using the Open Reading Frame (ORF) Finder tool from NCBI [51].

Results

The test of primers in 16 candidate genes allowed the sequence analysis of 12 amplicon sequences from ten genes in 108 Group Phureja accessions. This analysis led to the identification of 109 SNP markers, 44 of which were previously reported for tetraploid genotypes (Table 1). The position of these genes in the starch-sugar interconversion pathway can be observed in the scheme from Schreiber et al. [7]. Some enzyme acronyms that differ from those in the scheme are distinguished in the Table 1. The gen LapN is the exception, it is not shown in the pathway because it corresponds to a leucine aminopeptidase that was identified associated to the frying quality trait through a comparative proteomics approach [37]. From these candidate genes studied, five SNPs in two genes showed seven significant associations according to the mixed model and ANOVA analyses (Table 2). The positions of these SNPs are shown in the Additional file 2 in InvGE and SssI amplified sequences. The positions of the SNPs scored in the remaining eight genes studied are shown in the Additional file 3. Both additional files with the sequences indicate which markers correspond to the previously reported markers for tetraploid potatoes.
Table 2

SNP markers associated with sugar contents and chip darkening index in Solanum tuberosum Group Phureja

Gene-SNPa

Trait

Minor Allele

−log10(p-value)b

Percent variance

Frequency, %

explained and effectc

InvGE-C 2475454 T

Chip darkening index

2.0 (T)

2.4**

15.6 ↑

InvGE-G 2475344 A d

Fructose

9.4 (A)

2.2**

7.2 ↓

InvGE-C 2475295 A

Fructose

9.0 (C)

2.2***

7.0 ↓

InvGE-G 2475187 A

Fructose

9.0 (A)

2.7***

8.5 ↓

InvGE-G 2475187 A

Glucose

9.0 (A)

2.0*

6.6 ↓

InvGE-G 2475187 A

Total sugars

9.0 (A)

2.0***

8.8 ↓

SssI-C 45711901 T

Sucrose

35.5 (C)

2.0*

13.5 ↑

aThe SNP reference allele is followed by its position in the chromosome and the nucleotide from the allelic variant

bSNP markers with a − log10 (p-value) higher than two were considered with significant effect on the variables studied with the mixed model approach. The markers were confirmed with an ANOVA considering a level of significance of α = 0.05. p-values: ***, p ≤ 0.001; **, p ≤ 0.01; *, p ≤ 0.05

cThe percent variance explained expresses the effect of the minor allele in the trait. Arrows indicate the direction of the effect on the trait, upwards for a positive effect (low sugar content, therefore light chip color) and downwards for a negative effect (high sugar content, therefore dark chip color)

dBold SNP markers were in strong linkage disequilibrium, showed similar effects in the trait, and presented similar distributions in the population

The box plots in Fig. 1 show the effect of the genotypic classes from markers in the trait variation and supported the significance of the associations. The QQ-plots showed that the models adjusted for each phenotypic trait were effective in controlling type I error (Additional file 1). Model 1 was used for the analysis of chip darkening index and total sugars, while model 2 was used for the analysis of fructose, glucose, and sucrose contents.
Fig. 1

Box plots showing the allele effects of seven significant marker-trait associations: Chip darkening index (a), fructose (b-d), glucose (e), total sugars (f), and sucrose (g). Each plot is labeled with the SNP identification indicated in Table 2. Y-axis: Values for chip darkening index and sugar contents. X-axis: SNP genotypic classes

Four markers were found in the InvGE gene coding for an apoplastic invertase on potato chromosome IX. The average distance between adjacent markers was 88 base pairs (Additional file 2). Three of these markers were previously reported by Draffehn et al. [39] in the analysis of cDNA sequences from exon 3 in tetraploid cultivars and diploid genotypes from an experimental population. InvGE-G 2475344 A, InvGE-C 2475295 A, and InvGE-G 2475187 A corresponded to the markers snp688, snp639, and snp531 from the previous study, respectively. Two of the InvGE SNPs are synonymous and the other two are predicted to generate amino acid changes. InvGE-G 2475344 A caused a substitution of non-polar alanine by polar threonine, while InvGE-C 2475295 A causes a substitution of asparagine by lysine, both polar amino acids. The SNP InvGE-C 2475454 T was not found in the analysis of cDNA sequences of this exon [39]. It is therefore a novel variant found in Group Phureja which is associated with chip color.

The minor allele frequency of the SNPs in the InvGE gene was not found in the homozygous state in the population (Fig. 1af), which indicated that these alleles correspond to low-frequency variants [52]. The SNP InvGE-C 2475454 T presented the minor allele with the lowest frequency and with the greatest effect on the variable of chip darkening index (Table 2, Fig. 1a). The markers InvGE-G 2475344 A, InvGE-C 2475295 A, and InvGE-G 2475187 A (Fig. 1bf) were in strong LD with an average r2 of 0.94 ± 0.02 calculated from the pairwise values. The rare haplotype InvGE-A 2475187 C 2475295 A 2475344 showed association with increased fructose content. The SNP InvGE-G 2475187 A from the haplotype block showed multiple significant associations with fructose, glucose and total sugar contents. The corresponding minor frequency allele had a negative effect on the traits (Table 2, Fig. 1df). This LD region is supported by the similar effects of the SNPs on the variation of the traits and the similar distributions of the markers in the population (Table 2). Draffehn et al. [39] reported difficulties when sequencing exon 3 in a tetraploid mapping population due to the presence of insertion-deletion polymorphisms. Therefore, the report of this LD region in the exon 3 of InvGE and the association of the corresponding rare haplotype with increased sugar contents are novel in cultivated potato.

One SNP marker-trait association was found in an exon of the SssI gene coding for a soluble starch synthase on chromosome III (Additional file 2). In this case, the SNP represents a synonymous mutation where all three genotypic classes were observed (Table 2, Fig. 1g). The sequence flanking this marker has not been previously studied in tetraploids [7]. This SNP marker-trait association is therefore novel. The minor allele SssI-C 45711901 was associated with decreased sucrose content.

Discussion

The study of ten candidate genes with key functions in carbohydrate metabolism in Group Phureja allowed the identification of novel SNP markers that might be useful for further analyses of these genes in potato. We found SNP markers in two genes that were associated with sugar content and frying color in Group Phureja. These results show additional and novel genetic variability associated with frying quality in potato through the analysis of genes from Group Phureja which were associated in previous studies with sugar content and chip color variation in European tetraploid populations [7, 16, 17, 39]. The association analysis indicates the relevance of the InvGE and SssI loci in the variation of frying color and sugar content in diverse cultivated potato genetic backgrounds, despite the differences between Group Phureja and European tetraploid germplasm in their physiology and metabolism, ploidy level, allele frequencies and the different strategies for sampling and phenotyping the mapping populations [53].

The study of Group Phureja genotypes with a subtle population structure was appropriate to avoid spurious associations that are caused by highly stratified populations [9, 27]. Despite the current population size might be quite restricted (n = 112), the fast LD decay and the high heterozygosity found in the Phureja genotypes supports that the population is highly heterogeneous and provides good resolution for association mapping studies because it is representative of the diversity of this cultivated group [9, 27, 45]. Furthermore, the reliability of the reported associations is supported by the significant effects on the trait variation found with both mixed and linear model approaches.

InvGE and InvGF in chromosome IX are tandem duplicated genes both coding for apoplastic invertases that form a haplotype block associated with chip quality. They co-localize with the cold-sweetening QTL Sug9a, which was highly reproducible among environments [16, 39, 54, 55]. The association of SNP markers in the InvGE locus with chip darkening index, fructose, glucose, and total sugar content supports the reproducibility of the association in a diploid potato population. It also supports the correlations among these variables in Group Phureja as discussed previously [36]. These SNPs represent low-frequency variants (Table 2), thus adding evidence that such rare variants can have important effects in complex trait architecture [56].

The current study allowed the analysis of the sequence of exon 3 from InvGE in an association mapping population from Group Phureja. The sequence analysis discovered the novel SNP InvGE-C 2475454 T with a significant effect on the variation of the multi-environmental chip darkening index. The HPLC analysis in tubers harvested in an environment representative of Colombian highlands complemented the previous result by allowing the discovery of an InvGE haplotype that is associated with a negative effect on fructose content. Since a thorough analysis of this exonic region was not possible in a tetraploid association mapping population [39], it is not possible to establish whether this haplotype is typical for Group Phureja or if it is present also in tetraploid germplasm. The results also show contrasting effects on the trait variation of the marker InvGE-C 2475454 T and the adjacent haplotype, indicating that both regions are inherited independently because they do not form a single LD block. Additional research is required to validate the contrasting effects of the low-frequency variants present in these two gene regions given the outstanding positive effect of the rare allele InvGE-T 2475454 on light frying color. Consequently, it might be of relevance to generate and select genotypes homozygous for the T allele, in order to confirm the effect of this rare allele on light frying color.

The SNP markers InvGE-G 2475344 A and InvGE-C 2475295 A that were predicted to cause amino acid changes are expected to have an impact on protein function. Expression studies in tetraploids detected InvGE and InvGF transcripts in leaves and flowers but not in mature tubers [57]. In contrast, Maddison et al. [54] found InvGE expressed in the tissue surrounding tuber buds. Functional analyses of the corresponding potato alleles were not consistent with their statistical association with chip quality [57]. In fact, apoplastic invertases are relevant in regulating the import of glucose and fructose for tuber initiation and their expression might be time-dependent [58, 59]. Expression analyses have also shown that tissue expression of these invertases is genotype-dependent [57]. Because of the lack of conclusive evidence for a functional association of InvGE DNA variants with chip quality, further experiments are required to reveal the potential functional influence of apoplastic invertases on sugar content in mature tubers from Group Phureja. Nevertheless, experimental evidence suggests that the association of this locus is more likely to be indirect as a result of LD with a functional gene [57], which needs to be proven.

The possibility of a haplotype block in the distal region of chromosome IX is supported by the co-segregation of two microsatellite loci in landraces from Argentina. One locus is located in the gene InvGF while the other is in the gene SbeII that encodes a starch branching enzyme [60, 61]. Although SbeII is not predominantly expressed in tubers, antisense inhibition experiments have shown that it has a major effect on starch structure [62, 63]. Recently, an association study linked DNA polymorphisms in this gene to the degree of starch phosphorylation, a process related to starch degradation and sugar accumulation [41]. A linkage disequilibrium region on chromosome IX might contain InvGE, InvGF, and SbeII, the latter potentially having a strong functional association with sugar content in mature tubers. Therefore, it is relevant to explore in Group Phureja the presence of a haplotype block larger than the one reported here for InvGE, through the evaluation of markers in InvGF, SbeII, and other close genomic regions. The identification of associated regions in large haplotype blocks is favorable for breeding purposes because it might allow the design of diagnostic markers with high predictive value due to the low recombination rates in these regions [16].

The association of the SNP SssI-C 45711901 T reflects the effect of a marker with three genotypic classes showing a minor allele with a higher frequency than the minor alleles from markers in InvGE. Starch synthases catalyze the glycosyl transfer from ADP-glucose to glucan. The expression of SssI in tubers is low, suggesting a minor role of this particular enzyme in starch synthesis in storage organs [64]. As in the previous case, this SNP-trait association might be indirect, resulting from LD with causal variants, or direct through a potential novel role of SssI in the control of sugar content in tubers. The SNP in SssI is novel since sequence analyses in tetraploids have not reported this nucleotide variation [5]; therefore Group Phureja might be a new source of natural allelic variation in this gene with effect in low sucrose contents. The discovery of variants controlling sucrose content is relevant since the prevention of its accumulation is a strategy that can limit the subsequent accumulation of reducing sugars in tubers [65].

The marker-trait associations found explained a limited amount of the overall trait variability in Group Phureja. Therefore, selecting individuals which combine the superior alleles for these markers will not guarantee per se a better chip quality because there might be other regions influencing the trait, and the interaction with other loci might affect their expression [66]. Nevertheless, increasing the frequency of these alleles with positive effects in potato breeding populations can be a starting point for the marker-assisted improvement of frying quality, as well as for testing the reproducibility of the phenotypic effect of the associated markers in different genetic materials and environments [5]. There might be additional genes from the starch-sugar interconversion pathway associated to frying quality in Group Phureja that reflect the divergent metabolisms of both tetraploid and diploid genotypes. A next step would be the implementation of a genome wide approach to unravel additional genetic mechanisms and novel metabolic pathways underlying the frying quality trait of Group Phureja. These approaches would lead to the application of this knowledge on the genetic architecture of sugar contents and frying color to potato breeding programs.

Conclusions

The candidate gene approach allowed the identification of novel SNP markers in InvGE and SssI that were associated with frying quality and sugar contents in Group Phureja. These novel genetic associations need to be validated and afterwards might be useful in potato breeding programs. The results show that some genes related to sugar content and frying color variation are conserved in diploid and tetraploid cultivated potato germplasm. There is a lack of conclusive evidence for a functional association of both genes with frying quality. Therefore, further experiments are required to reveal the potential functional influence of the associated SNPs with sugar contents of mature tubers or to demonstrate that associations are indirect as a result of LD with functional genes. The implementation of a genome wide approach is necessary to analyze comprehensively the genomic regions that control starch-sugar equilibrium. The genome wide study is also required to assess the distribution of LD regions in all the chromosomes to establish the mapping resolution. These strategies will contribute to a broad understanding of the genetic architecture of sugar accumulation and frying color of Group Phureja.

Declarations

Acknowledgements

We acknowledge the Ministry of Environment and Sustainable Development of Colombia for the granted access to S. tuberosum Group Phureja genetic resources under the contract RGE0053. The authors would like to express their gratitude to Birgit Walkemeier and María Cecilia Delgado for their support in the development of this research and to Margarita Duarte for her help in the final edition of the manuscript.

Funding

This research was funded by the International Development Research Centre (IDRC) and Global Affairs Canada (GAC) through the Canadian International Food Security Research Fund that financed the project “Improving Potato Production for Increasing Food Security of Indigenous Communities in Colombia”. Additional funding was provided by the Max Planck Society and the National University of Colombia.

Availability of data and materials

The data supporting the results of this article are included within the article and its additional files. Additional datasets of interest are available from the corresponding author on a reasonable request.

Authors’ contributions

DDD carried out all the experiments, the data analyses, and wrote the manuscript. DJ and FS supported the analysis and interpretation of the data. CG supervised and participated in the design of the experiments with candidate genes. TMV supervised the entire research, participated in the design of the experiments and in the analysis and interpretation of the data. All authors participated in writing the manuscript and approved the final version.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Authors’ Affiliations

(1)
Faculty of Agricultural Sciences, Agronomy Department, National University of Colombia
(2)
Present address: INRES-Plant Breeding, University of Bonn
(3)
Department of Plant Breeding and Genetics, Max Planck Institute for Plant Breeding Research
(4)
Faculty of Sciences, Biology Department, National University of Colombia

References

  1. Halford NG, Curtis TY, Muttucumaru N, Postles J, Elmore JS, Mottram DS. The acrylamide problem: a plant and agronomic science issue. J Exp Bot. 2012;63:2841–51.View ArticlePubMedGoogle Scholar
  2. Isla MI, Vattuone MA, Sampietro AR. Hydrolysis of sucrose within isolated vacuoles from Solanum tuberosum L. tubers. Planta. 1998;205:601–5.View ArticlePubMedGoogle Scholar
  3. Kumar D, Singh BP, Kumar P. An overview of the factors affecting sugar content of potatoes. Ann Appl Biol. 2004;145:247–56.View ArticleGoogle Scholar
  4. Eck HA v. Genetics of morphological and tuber traits. In: Vreugdenhil D, Bradshaw J, Gebhardt C, Govers F, MacKerron DKL, Taylor MA, et al., editors. Potato biology and biotechnology: Advances and perspectives. Amsterdam: Elsevier Science B.V; 2007. p. 91–111.Google Scholar
  5. Roitsch T, González M-C. Function and regulation of plant invertases: sweet sensations. Trends Plant Sci. 2004;9:606–13.View ArticlePubMedGoogle Scholar
  6. Ruan Y-L. Sucrose metabolism: gateway to diverse carbon use and sugar signaling. Annu Rev Plant Biol. 2014;65:33–67.View ArticlePubMedGoogle Scholar
  7. 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 Genes Genom. Genet. 2014;4:1797–811.Google Scholar
  8. Sturm A, Tang GQ. The sucrose-cleaving enzymes of plants are crucial for development, growth and carbon partitioning. Trends Plant Sci. 1999;4:401–7.View ArticlePubMedGoogle Scholar
  9. Alvarez MF, Mosquera T, Blair MW. The use of association genetics approaches in plant breeding. In: Janick J, editor. Plant Breeding Reviews: Vol. 38. Hoboken: John Wiley & Sons, Inc; 2014. p. 17–68.Google Scholar
  10. Hamblin MT, Buckler ES, Jannink J-L. Population genetics of genomics-based crop improvement methods. Trends Genet. 2011;27:98–106.View ArticlePubMedGoogle Scholar
  11. Körber N, Bus A, Li J, Higgins J, Bancroft I, Higgins EE, et al. Seedling development traits in Brassica napus examined by gene expression analysis and association mapping. BMC Plant Biol. 2015;15:136.View ArticlePubMedPubMed CentralGoogle Scholar
  12. Gebhardt C, Li L, Pajerowska-Mukhtar K, Achenbach U, Sattarzadeh A, Bormann C, et al. Candidate gene approach to identify genes underlying quantitative traits and develop diagnostic markers in potato. Crop Sci. 2007;47:S 106–11.View ArticleGoogle Scholar
  13. Chen X, Salamini F, Gebhardt C. A potato molecular-function map for carbohydrate metabolism and transport. Theor Appl Genet. 2001;102:284–95.View ArticleGoogle Scholar
  14. Baldwin SJ, Dodds KG, Auvray B, Genet RA, Macknight RC, Jacobs JME. Association mapping of cold-induced sweetening in potato using historical phenotypic data. Ann Appl Biol. 2011;158:248–56.View ArticleGoogle Scholar
  15. Li L, Tacke E, Hofferbert H-R, 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.View ArticlePubMedPubMed CentralGoogle Scholar
  16. Li L, Paulo M-J, Strahwald J, Lübeck J, Hofferbert H-R, 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.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Li L, Strahwald J, Hofferbert H-R, Lubeck J, Tacke E, Junghans H, et al. DNA Variation at the invertase locus invGE/GF is associated with tuber quality traits in populations of potato breeding clones. Genetics. 2005;170:813–21.View ArticlePubMedPubMed CentralGoogle Scholar
  18. Huamán Z, Spooner DM. Reclassification of landrace populations of cultivated potatoes (Solanum sect. Petota). Am J Bot. 2002;89:947–65.View ArticlePubMedGoogle Scholar
  19. Mosquera T, Cuéllar D. La genética de la papa frente a la gota. Bogotá: Editorial Universidad Nacional de Colombia; 2013.Google Scholar
  20. Peña C, Restrepo-Sánchez L-P, Kushalappa A, Rodríguez-Molano L-E, Mosquera T, Narváez-Cuenca C-E. Nutritional contents of advanced breeding clones of Solanum tuberosum group Phureja. LWT - Food Sci. Technol. 2015;62:76–82.Google Scholar
  21. Rivera JE, Herrera A, Rodríguez LE. Evaluación sensorial en productos procesados de papa criolla (Solanum phureja) y su importancia para el fitomejoramiento. Fitotec Colomb. 2006;6:9–25.Google Scholar
  22. Rodríguez LE, Ñustez CE, Estrada N. Criolla Latina, Criolla Paisa and Criolla Colombia, new cultivars of “diploid potato” for the province of Antioquia (Colombia). Agron Colomb. 2009;27:289–303.Google Scholar
  23. Piñeros-Niño C, Narváez-Cuenca C-E, Kushalappa AC, Mosquera T. Hydroxycinnamic acids in cooked potato tubers from Solanum tuberosum group Phureja. Food Sci. Nutr. 2016; doi:10.1002/fsn3.403.
  24. Maine MJD, Carroll CP, Stewart HE, Solomon RM, Wastie RL. Disease resistance in Solanum phureja and diploid and tetraploid S. tuberosum × S. phureja hybrids. Potato Res. 1993;36:21–8.View ArticleGoogle Scholar
  25. Hayes RJ, Thill CA. Introgression of cold (4 C) chipping from 2× (2 endosperm balance number) potato species into 4× (4EBN) cultivated potato using sexual polyploidization. Am J Potato Res. 2002;79:421–31.View ArticleGoogle Scholar
  26. Machida-Hirano R. Diversity of potato genetic resources. Breed Sci. 2015;65:26–40.
  27. Juyó D, Sarmiento F, Álvarez M, Brochero H, Gebhardt C, Mosquera T. Genetic diversity and population structure in diploid potatoes of Group Phureja. Crop Sci. 2015;55:760–9.
  28. Potato Genome Sequencing Consortium, Xu X, Pan S, Cheng S, Zhang B, Mu D, et al. Genome sequence and analysis of the tuber crop potato. Nature. 2011;475:189–95.View ArticleGoogle Scholar
  29. Jupe F, Witek K, Verweij W, Sliwka J, Pritchard L, Etherington GJ, et al. Resistance gene enrichment sequencing (RenSeq) enables reannotation of the NB-LRR gene family from sequenced plant genomes and rapid mapping of resistance loci in segregating populations. Plant J Cell Mol Biol. 2013;76:530–44.View ArticleGoogle Scholar
  30. Werij JS, Furrer H, Van Eck HJ, Visser RGF, Bachem CWB. A limited set of starch related genes explain several interrelated traits in potato. Euphytica. 2012;186:501–16.View ArticleGoogle Scholar
  31. Ersoz ES, Yu J, Buckler ES. Applications of linkage disequilibrium and association mapping in crop Plants. In: Varshney RK, Tuberosa R, editors. Genomics-assisted crop improvement: Vol. Dordrecht: Springer Netherlands; 2007. p. 97–119.View ArticleGoogle Scholar
  32. Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, et al. Finding the missing heritability of complex diseases. Nature. 2009;461:747–53.View ArticlePubMedPubMed CentralGoogle Scholar
  33. McClements DJ. Analysis of carbohydrates. In: University of Massachussets. 2003. http://people.umass.edu/~mcclemen/581Carbohydrates.html. Accessed 9 Apr 2014
  34. Duarte-Delgado D, Narváez-Cuenca C-E, Restrepo-Sánchez L-P, Kushalappa A, Mosquera-Vásquez T. Development and validation of a liquid chromatographic method to quantify sucrose, glucose, and fructose in tubers of Solanum tuberosum group Phureja. J Chromatogr B. 2015;975:18–23.View ArticleGoogle Scholar
  35. Ñústez-López CE. Estudios fenotípicos y genéticos asociados a la calidad de fritura en Solanum phureja Juz et Buk. PhD thesis, National University of Colombia. 2011. http://www.bdigital.unal.edu.co/4318/. Accessed 26 Aug 2015
  36. Duarte-Delgado D, Ñústez-López C-E, Narváez-Cuenca C-E, Restrepo-Sánchez L-P, Melo SE, Sarmiento F, et al. Natural variation of sucrose, glucose and fructose contents in Colombian genotypes of Solanum tuberosum Group Phureja at harvest. J Sci Food Agric. 2016;96:4288–94.View ArticlePubMedPubMed CentralGoogle Scholar
  37. Fischer M, Schreiber L, Colby T, Kuckenberg M, Tacke E, Hofferbert H-R, 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.View ArticlePubMedPubMed CentralGoogle Scholar
  38. Bonierbale M, De Haan S, Forbes A. Procedures for standard evaluation trials of advanced potato clones. An international cooperator’s guide. Lima: International Potato Center; 2007.Google Scholar
  39. Draffehn AM, Meller S, Li L, Gebhardt C. Natural diversity of potato (Solanum tuberosum) invertases. BMC Plant Biol. 2010;10:271.View ArticlePubMedPubMed CentralGoogle Scholar
  40. Draffehn AM, Durek P, Nunes-Nesi A, Stich B, Fernie AR, Gebhardt C. Tapping natural variation at functional level reveals allele specific molecular characteristics of potato invertase Pain-1. Plant Cell Environ. 2012;35:2143–54.View ArticlePubMedGoogle Scholar
  41. Carpenter MA, Joyce NI, Genet RA, Cooper RD, Murray SR, Noble AD, et al. Starch phosphorylation in potato tubers is influenced by allelic variation in the genes encoding glucan water dikinase, starch branching enzymes I and II, and starch synthase III. Plant Genet Genomics. 2015;6:143.Google Scholar
  42. Sowokinos JR, Vigdorovich V, Abrahamsen M. Molecular cloning and sequence variation of UDP-glucose pyrophosphorylase cDNAs from potatoes sensitive and resistant to cold sweetening. J Plant Physiol. 2004;161:947–55.View ArticlePubMedGoogle Scholar
  43. Sharma SK, Bolser D, De Boer J, Sønderkær M, Amoros W, Carboni MF, et al. Construction of reference chromosome-scale pseudomolecules for potato: integrating the potato genome with genetic and physical maps. G3 Genes Genom. Genet. 2013;3:2031–47.Google Scholar
  44. Hirsch CD, Hamilton JP, Childs KL, Cepela J, Crisovan E, Vaillancourt B, et al. Spud DB: A resource for mining sequences, genotypes, and phenotypes to accelerate potato breeding. Plant Genome. 2014;7.
  45. Álvarez, M.F. Identification of molecular marker associated with polygenic resistance to Phytophthora infestans trough association mapping in Solanum tuberosum group Phureja. PhD thesis, National University of Colombia. 2014. http://www.bdigital.unal.edu.co/46589/. Accessed 12 Dec 2016
  46. Zhang Z, Ersoz E, Lai C-Q, Todhunter RJ, Tiwari HK, Gore MA, et al. Mixed linear model approach adapted for genome-wide association studies. Nat Genet. 2010;42:355–60.View ArticlePubMedPubMed CentralGoogle Scholar
  47. Li M, Liu X, Bradbury P, Yu J, Zhang Y-M, Todhunter RJ, et al. Enrichment of statistical power for genome-wide association studies. BMC Biol. 2014;12:73.View ArticlePubMedPubMed CentralGoogle Scholar
  48. Lipka AE, Tian F, Wang Q, Peiffer J, Li M, Bradbury PJ, et al. GAPIT: genome association and prediction integrated tool. Bioinforma Oxf Engl. 2012;28:2397–9.View ArticleGoogle Scholar
  49. Sukumaran S, Xiang W, Bean SR, Pedersen JF, Kresovich S, Tuinstra MR, et al. Association mapping for grain quality in a diverse Sorghum collection. Plant Genome J. 2012;5:126–35.
  50. Hill WG, Robertson A. Linkage disequilibrium in finite populations. Theor Appl Genet. 1968;38:226–31.View ArticlePubMedGoogle Scholar
  51. Wheeler DL, Church DM, Federhen S, Lash AE, Madden TL, Pontius JU, et al. Database resources of the National Center for Biotechnology. Nucleic Acids Res. 2003;31:28–33.View ArticlePubMedPubMed CentralGoogle Scholar
  52. Mitchell AA, Zwick ME, Chakravarti A, Cutler DJ. Discrepancies in dbSNP confirmation rates and allele frequency distributions from varying genotyping error rates and patterns. Bioinforma Oxf Engl. 2004;20:1022–32.View ArticleGoogle Scholar
  53. Schönhals EM, Ortega F, Barandalla L, Aragones A, Galarreta JIRD, Liao J-C, et al. Identification and reproducibility of diagnostic DNA markers for tuber starch and yield optimization in a novel association mapping population of potato (Solanum tuberosum L.). Theor Appl Genet. 2016;129:767–85.View ArticlePubMedPubMed CentralGoogle Scholar
  54. Maddison AL, Hedley PE, Meyer RC, Aziz N, Davidson D, Machray GC. Expression of tandem invertase genes associated with sexual and vegetative growth cycles in potato. Plant Mol Biol. 1999;41:741–51.View ArticlePubMedGoogle Scholar
  55. Menéndez CM, Ritter E, Schäfer-Pregl R, Walkemeier B, Kalde A, Salamini F, et al. Cold sweetening in diploid potato: mapping quantitative trait loci and candidate genes. Genetics. 2002;162:1423–34.PubMedPubMed CentralGoogle Scholar
  56. Panoutsopoulou K, Tachmazidou I, Zeggini E. In search of low-frequency and rare variants affecting complex traits. Hum Mol Genet. 2013;22:R16–21.View ArticlePubMedPubMed CentralGoogle Scholar
  57. Draffehn A. Structural and functional characterization of natural alleles of potato (Solanum tuberosum L.) invertases associated with tuber quality traits. PhD thesis, University of Cologne. 2010. http://kups.ub.uni-koeln.de/3121/. Accessed 11 May 2015
  58. Fotopoulos V. Plant invertases: structure, function and regulation of a diverse enzyme family. J Biol Res. 2005;4:127–37.Google Scholar
  59. Minhas JS, Rai VK, Saini HS. Carbohydrate metabolism during tuber initiation in potato: A transient surge in invertase activity marks the stolon to tuber transition. Potato Res. 2004;47:113–26.View ArticleGoogle Scholar
  60. Colman SL, Monti MC, Divito SB, Digilio A, Feingold SE. Marcadores funcionales asociados al endulzamiento inducido por frío en papas nativas de Argentina. Revista Latinoamericana de la Papa. 2009;15:61–5.Google Scholar
  61. Feingold S, Lloyd J, Norero N, Bonierbale M, Lorenzen J. Mapping and characterization of new EST-derived microsatellites for potato (Solanum tuberosum L.). Theor Appl Genet. 2005;111:456–66.View ArticlePubMedGoogle Scholar
  62. Jobling SA, Schwall GP, Westcott RJ, Sidebottom CM, Debet M, Gidley MJ, et al. A minor form of starch branching enzyme in potato (Solanum tuberosum L.) tubers has a major effect on starch structure: cloning and characterisation of multiple forms of SBE A. Plant J. 1999;18:163–71.View ArticlePubMedGoogle Scholar
  63. Larsson CT, Khoshnoodi J, Ek B, Rask L, Larsson H. Molecular cloning and characterization of starch-branching enzyme II from potato. Plant Mol Biol. 1998;37:505–11.View ArticlePubMedGoogle Scholar
  64. Kossmann J, Abel GJ, Springer F, Lloyd JR, Willmitzer L. Cloning and functional analysis of a cDNA encoding a starch synthase from potato (Solanum tuberosum L.) that is predominantly expressed in leaf tissue. Planta. 1999;208:503–11.View ArticlePubMedGoogle Scholar
  65. Dale MFB, Bradshaw JE. Progress in improving processing attributes in potato. Trends Plant Sci. 2003;8:310–2.View ArticlePubMedGoogle Scholar
  66. Li L, Paulo M-J, Van Eeuwijk F, Gebhardt C. Statistical epistasis between candidate gene alleles for complex tuber traits in an association mapping population of tetraploid potato. Theor Appl Genet. 2010;121:1303–10.View ArticlePubMedPubMed CentralGoogle Scholar
  67. Nader-Nieto C. Structural and functional characterization of natural alleles of potato (Solanum tuberosum L.) starch phosphorylases associated with tuber quality traits. PhD thesis, University of Cologne; 2011. http://kups.ub.uni-koeln.de/4791. Accessed 20 Aug 2013

Copyright

© The Author(s). 2017