SNPs in genes functional in starch-sugar interconversion associate with natural variation of tuber starch and sugar content of potato (Solanum tuberosum L.)

The effect of marker types and density on genomic prediction and GWAS of key performance traits in tetraploid potato

Genomic prediction and genome-wide association studies are becoming widely employed in potato key performance trait QTL identifications and to support potato breeding using genomic selection. Elite cultivars are tetraploid and highly heterozygous but also share many common ancestors and generation-spanning inbreeding events, resulting from the clonal propagation of potatoes through seed potatoes. Consequentially, many SNP markers are not in a 1:1 relationship with a single allele variant but shared over several alleles that might exert varying effects on a given trait. The impact of such redundant "diluted" predictors on the statistical models underpinning genome-wide association studies (GWAS) and genomic prediction has scarcely been evaluated despite the potential impact on model accuracy and performance. We evaluated the impact of marker location, marker type, and marker density on the genomic prediction and GWAS of five key performance traits in tetraploid potato (chipping quality, dry matter content, length/width ratio, senescence, and yield). A 762-offspring panel of a diallel cross of 18 elite cultivars was genotyped by sequencing, and markers were annotated according to a reference genome. Genomic prediction models (GBLUP) were trained on four marker subsets [non-synonymous (29,553 SNPs), synonymous (31,229), non-coding (32,388), and a combination], and robustness to marker reduction was investigated. Single-marker regression GWAS was performed for each trait and marker subset. The best cross-validated prediction correlation coefficients of 0.54, 0.75, 0.49, 0.35, and 0.28 were obtained for chipping quality, dry matter content, length/width ratio, senescence, and yield, respectively. The trait prediction abilities were similar across all marker types, with only non-synonymous variants improving yield predictive ability by 16%. Marker reduction response did not depend on marker type but rather on trait. Traits with high predictive abilities, e.g., dry matter content, reached a plateau using fewer markers than traits with intermediate-low correlations, such as yield. The predictions were unbiased across all traits, marker types, and all marker densities >100 SNPs. Our results suggest that using non-synonymous variants does not enhance the performance of genomic prediction of most traits. The major known QTLs were identified by GWAS and were reproducible across exonic and whole-genome variant sets for dry matter content, length/width ratio, and senescence. In contrast, minor QTL detection was marker type dependent.

Pho1a (plastid starch phosphorylase) is duplicated and essential for normal starch granule phenotype in tubers of Solanum tuberosum L

Reserve starch from seeds and tubers is a crucial plant product for human survival. Much research has been devoted to quantitative and qualitative aspects of starch synthesis and its relation to abiotic factors of importance in agriculture. Certain aspects of genetic factors and enzymes influencing carbon assimilation into starch granules remain elusive after many decades of research. Starch phosphorylase (Pho) can operate, depending on metabolic conditions, in a synthetic and degradative pathway. The plastidial form of the enzyme is one of the most highly expressed genes in potato tubers, and the encoded product is imported into starch-synthesizing amyloplasts. We identified that the genomic locus of a Pho1a-type starch phosphorylase is duplicated in potato. Our study further shows that the enzyme is of importance for a normal starch granule phenotype in tubers. Null mutants created by genome editing display rounded starch granules in an increased number that contained a reduced ratio of apparent amylose in the starch.

A physical map of traits of agronomic importance based on potato and tomato genome sequences

Potato, tomato, pepper, and eggplant are worldwide important crop and vegetable species of the Solanaceae family. Molecular linkage maps of these plants have been constructed and used to map qualitative and quantitative traits of agronomic importance. This research has been undertaken with the vision to identify the molecular basis of agronomic characters on the one hand, and on the other hand, to assist the selection of improved varieties in breeding programs by providing DNA-based markers that are diagnostic for specific agronomic characters. Since 2011, whole genome sequences of tomato and potato became available in public databases. They were used to combine the results of several hundred mapping and map-based cloning studies of phenotypic characters between 1988 and 2022 in physical maps of the twelve tomato and potato chromosomes. The traits evaluated were qualitative and quantitative resistance to pathogenic oomycetes, fungi, bacteria, viruses, nematodes, and insects. Furthermore, quantitative trait loci for yield and sugar content of tomato fruits and potato tubers and maturity or earliness were physically mapped. Cloned genes for pathogen resistance, a few genes underlying quantitative trait loci for yield, sugar content, and maturity, and several hundred candidate genes for these traits were included in the physical maps. The comparison between the physical chromosome maps revealed, in addition to known intrachromosomal inversions, several additional inversions and translocations between the otherwise highly collinear tomato and potato genomes. The integration of the positional information from independent mapping studies revealed the colocalization of qualitative and quantitative loci for resistance to different types of pathogens, called resistance hotspots, suggesting a similar molecular basis. Synteny between potato and tomato with respect to genomic positions of quantitative trait loci was frequently observed, indicating eventual similarity between the underlying genes.

Investigating the genetic components of tuber bruising in a breeding population of tetraploid potatoes

BACKGROUND

Tuber bruising in tetraploid potatoes (Solanum tuberosum) is a trait of economic importance, as it affects tubers' fitness for sale. Understanding the genetic components affecting tuber bruising is a key step in developing potato lines with increased resistance to bruising. As the tetraploid setting renders genetic analyses more complex, there is still much to learn about this complex phenotype. Here, we used capture sequencing data on a panel of half-sibling populations from a breeding programme to perform a genome-wide association analysis (GWAS) for tuber bruising. In addition, we collected transcriptomic data to enrich the GWAS results. However, there is currently no satisfactory method to represent both GWAS and transcriptomics analysis results in a single visualisation and to compare them with existing knowledge about the biological system under study.

RESULTS

When investigating population structure, we found that the STRUCTURE algorithm yielded greater insights than discriminant analysis of principal components (DAPC). Importantly, we found that markers with the highest (though non-significant) association scores were consistent with previous findings on tuber bruising. In addition, new genomic regions were found to be associated with tuber bruising. The GWAS results were backed by the transcriptomics differential expression analysis. The differential expression notably highlighted for the first time the role of two genes involved in cellular strength and mechanical force sensing in tuber resistance to bruising. We proposed a new visualisation, the HIDECAN plot, to integrate the results from the genomics and transcriptomics analyses, along with previous knowledge about genomic regions and candidate genes associated with the trait.

CONCLUSION

This study offers a unique genome-wide exploration of the genetic components of tuber bruising. The role of genetic components affecting cellular strength and resistance to physical force, as well as mechanosensing mechanisms, was highlighted for the first time in the context of tuber bruising. We showcase the usefulness of genomic data from breeding programmes in identifying genomic regions whose association with the trait of interest merit further investigation. We demonstrate how confidence in these discoveries and their biological relevance can be increased by integrating results from transcriptomics analyses. The newly proposed visualisation provides a clear framework to summarise of both genomics and transcriptomics analyses, and places them in the context of previous knowledge on the trait of interest.

Genome-Wide Association Study of Agronomic and Physiological Traits Related to Drought Tolerance in Potato

Potato (Solanum tuberosum L.) is often considered a water-sensitive crop and its production can be threatened by drought events, making water stress tolerance a trait of increasing interest. In this study, a panel of 144 tetraploid potato genotypes was evaluated for two consecutive years (2019 and 2020) to observe the variation of several physiological traits such as chlorophyll content and fluorescence, stomatal conductance, NDVI, and leaf area and circumference. In addition, agronomic parameters such as yield, tuber fresh weight, tuber number, starch content, dry matter and reducing sugars were determined. GGP V3 Potato array was used to genotype the population, obtaining a total of 18,259 high-quality SNP markers. Marker-trait association was performed using GWASpoly package in R software and Q + K linear mixed models were considered. This approach allowed us to identify eighteen SNP markers significantly associated with the studied traits in both treatments and years, which were related to genes with known functions. Markers related to chlorophyll content and number of tubers under control and stress conditions, and related to stomatal conductance, NDVI, yield and reducing sugar content under water stress, were identified. Although these markers were distributed throughout the genome, the SNPs associated with the traits under control conditions were found mainly on chromosome 11, while under stress conditions they were detected on chromosome 4. These results contribute to the knowledge of the mechanisms of potato tolerance to water stress and are useful for future marker-assisted selection programs.

Recent Advances in Molecular Improvement for Potato Tuber Traits

Potato is an important crop due to its nutritional value and high yield potential. Improving the quality and quantity of tubers remains one of the most important breeding objectives. Genetic mapping helps to identify suitable markers for use in the molecular breeding, and combined with transgenic approaches provides an efficient way for gaining desirable traits. The advanced plant breeding tools and molecular techniques, e.g., TALENS, CRISPR-Cas9, RNAi, and cisgenesis, have been successfully used to improve the yield and nutritional value of potatoes in an increasing world population scenario. The emerging methods like genome editing tools can avoid incorporating transgene to keep the food more secure. Multiple success cases have been documented in genome editing literature. Recent advances in potato breeding and transgenic approaches to improve tuber quality and quantity have been summarized in this review.

The transcription factor StTINY3 enhances cold-induced sweetening resistance by coordinating starch resynthesis and sucrose hydrolysis in potato

The accumulation of reducing sugars in cold-stored tubers, known as cold-induced sweetening (CIS), negatively affects potato processing quality. The starch to sugar interconversion pathways that are altered in cold-stored CIS tubers have been elucidated, but the mechanism that regulates them remains largely unknown. This study identified a CBF/DREB transcription factor (StTINY3) that enhances CIS resistance by both activating starch biosynthesis and repressing the hydrolysis of sucrose to reducing sugars in detached cold-stored tubers. Silencing StTINY3 in a CIS-resistant genotype decreased CIS resistance, while overexpressing StTINY3 in a CIS-sensitive genotype increased CIS resistance, and altering StTINY3 expression was associated with expression changes in starch resynthesis-related genes. We showed first that overexpressing StTINY3 inhibited sucrose hydrolysis by enhancing expression of the invertase inhibitor gene StInvInh2, and second that StTINY3 promoted starch resynthesis by up-regulating a large subunit of the ADP-glucose pyrophosphorylase gene StAGPaseL3, and the glucose-6-phosphate transporter gene StG6PT2. Using electrophoretic mobility shift assays, we revealed that StTINY3 is a nuclear-localized transcriptional activator that directly binds to the dehydration-responsive element/CRT cis-element in the promoters of StInvInh2 and StAGPaseL3. Taken together, these findings established that StTINY3 influences CIS resistance in cold-stored tubers by coordinately modulating the starch to sugar interconversion pathways and is a good target for improving potato processing quality.

Transcription factor OsSGL is a regulator of starch synthesis and grain quality in rice

Starch biosynthesis during rice endosperm development is important for grain quality, as it influences grain size and physico-chemical properties, which together determine rice eating quality. Cereal starch biosynthetic pathways have been comprehensively investigated; however, their regulation, especially by transcriptional repressors remains largely unknown. Here, we identified a DUF1645 domain-containing protein, STRESS_tolerance and GRAIN_LENGTH (OsSGL), that participates in regulating rice starch biosynthesis. Overexpression of OsSGL reduced total starch and amylose content in the endosperm compared with the wild type. Chromatin immunoprecipitation sequencing and RNA-seq analyses indicated that OsSGL targets the transcriptional activity of several starch and sucrose metabolism genes. In addition, ChIP-qPCR, yeast one-hybrid, EMSA and dual-luciferase assays demonstrated that OsSGL directly inhibits the expression of SUCROSE SYNTHASE 1 (OsSUS1) in the endosperm. Furthermore, OsSUS1 interacts with OsSGL to release its transcriptional repression ability. Unexpectedly, our results also show that knock down and mutation of OsSGL disrupts the starch biosynthetic pathway, causing lower starch and amylose content. Therefore, our findings demonstrate that accurate control of OsSGL homeostasis is essential for starch synthesis and grain quality. In addition, we revealed the molecular mechanism of OsSGL in regulating starch biosynthesis-related genes, which are required for grain quality.

Isolation and functional characterization of a glucose-6-phosphate/phosphate translocator (IbG6PPT1) from sweet potato (Ipomoea batatas (L.) Lam.)

Sweet potato (Ipomoea batatas (L.) Lam.) is a good source of carbohydrates, an excellent raw material for starch-based industries, and a strong candidate for biofuel production due to its high starch content. However, the molecular basis of starch biosynthesis and accumulation in sweet potato is still insufficiently understood. Glucose-6-phosphate/phosphate translocators (GPTs) mediate the import of glucose-6-phosphate (Glc6P) into plastids for starch synthesis. Here, we report the isolation of a GPT-encoding gene, IbG6PPT1, from sweet potato and the identification of two additional IbG6PPT1 gene copies in the sweet potato genome. IbG6PPT1 encodes a chloroplast membrane-localized GPT belonging to the GPT1 group and highly expressed in storage root of sweet potato. Heterologous expression of IbG6PPT1 resulted in increased starch content in the leaves, root tips, and seeds and soluble sugar in seeds of Arabidopsis thaliana, but a reduction in soluble sugar in the leaves. These findings suggested that IbG6PPT1 might play a critical role in the distribution of carbon sources in source and sink and the accumulation of carbohydrates in storage tissues and would be a good candidate gene for controlling critical starch properties in sweet potato.

Comparative Transcriptome Sequencing of Taro Corm Development With a Focus on the Starch and Sucrose Metabolism Pathway

Taro (Colocasia esculenta) is an important tuber crop and staple food. Taro corms have higher nutritional value and starch contents as compared to most of the other root/tuber crops. However, the growth and development of the taro rhizome have not been critically examined in terms of transcriptomic signatures in general or specific to carbohydrates (starch and sucrose) accumulation. In current study, we have conducted a comprehensive survey of transcripts in taro corms aged 1, 2, 3, 4, 5, and 8 months. In this context, we have employed a whole transcriptome sequencing approach for identification of mRNAs, CircRNAs, and miRNAs in corms and performed functional enrichment analysis of the screened differentially expressed RNAs. A total of 11,203 mRNAs, 245 CircRNAs, and 299 miRNAs were obtained from six developmental stages. The mRNAs included 139 DEGs associated with 24 important enzymes of starch and sucrose metabolism. The expression of genes encoding key enzymes of starch and sucrose metabolism pathway (GBSS, AGPase, UGPase, SP, SSS, βFRUCT and SuSy) demonstrated significant variations at the stage of 4 months (S4). A total of 191 CircRNAs were differentially expressed between the studied comparisons of growth stages and 99 of these were associated with those miRNA (or target genes) that were enriched in starch and sucrose metabolism pathway. We also identified 205 miRNAs including 46 miRNAs targeting DEGs enriched in starch and sucrose biosynthesis pathway. The results of current study provide valuable resources for future exploration of the molecular mechanisms involved in the starch properties of Taro.

Discovery of SNPs and InDels in papaya genotypes and its potential for marker assisted selection of fruit quality traits

Papaya is a tropical and climacteric fruit that is recognized for its nutritional benefits and medicinal applications. Its fruits ripen quickly and show a drastic fruit softening, leading to great post-harvest losses. To overcome this scenario, breeding programs of papaya must invest in exploring the available genetic variation to continue developing superior cultivars with improved fruit quality traits. The objective of this study was to perform a whole-genome genotyping (WGG) of papaya, predict the effects of the identified variants, and develop a list of ripening-related genes (RRGs) with linked variants. The Formosa elite lines of papaya Sekati and JS-12 were submitted to WGG with an Illumina Miseq platform. The effects of variants were predicted using the snpEff program. A total of 28,451 SNPs having Ts/Tv (Transition/Transversion) ratio of 2.45 and 1,982 small insertions/deletions (InDels) were identified. Most variant effects were predicted in non-coding regions, with only 2,104 and 138 effects placed in exons and splice site regions, respectively. A total of 106 RRGs were found to be associated with 460 variants, which may be converted into PCR markers to facilitate genetic mapping and diversity studies and to apply marker-assisted selection (MAS) for specific traits in papaya breeding programs.

eQTL mapping of the 12S globulin cruciferin gene PGCRURSE5 as a novel candidate associated with starch content in potato tubers

Tuber starch content (TSC) is a very important trait in potato (Solanum tuberosum L.). This study is the first to use expression quantitative trait loci (eQTL) mapping of transcript-derived markers for TSC in potato. Thirty-four differentially expressed genes were selected by comparing the RNA-seq data of contrasting bulked segregants. For the 11 candidate genes, we determined their relative expression levels across the segregating diploid potato population using RT-qPCR. We detected 36 eQTL as candidate genes distributed on all twelve potato chromosomes, and nine of them overlapped with QTL for TSC. Peaks for two eQTL, eAGPaseS-a and ePGRCRURSE5, were close to the corresponding loci of the large subunit of ADP-glucose pyrophosphorylase (AGPaseS-a) and the 12S globulin cruciferin gene (PGCRURSE5), respectively. The eQTL peaks for AGPaseS-a and PGRCRURSE5 explained 41.0 and 28.3% of the phenotypic variation at the transcript level. We showed the association of the DNA markers for AGPaseS-a and PGRCRURSE5 with QTL for TSC, and significant correlation between the expression level of PGRCRURSE5 and TSC. We did not observe a significant correlation between the expression level of AGPaseS-a and TSC. We concluded that the cruciferin gene PGRCRURSE5 is a novel candidate involved in the regulation of starch content in potato tubers.

Non-GMO potato lines, synthesizing increased amylose and resistant starch, are mainly deficient in isoamylase debranching enzyme

Solanum tuberosum potato lines with high amylose content were generated by crossing with the wild potato species Solanum sandemanii followed by repeated backcrossing to Solanum tuberosum lines. The trait, termed increased amylose (IAm), was recessive and present after three generations of backcrossing into S. tuberosum lines (6.25% S. sandemanii genes). The tubers of these lines were small, elongated and irregular with small and misshaped starch granules and high sugar content. Additional backcrossing resulted in less irregular tuber morphology, increased starch content (4.3%-9.5%) and increased amylose content (29%-37.9%) but indifferent sugar content. The amylose in the IAm starch granules was mainly located in peripheral spots, and large cavities were found in the granules. Starch pasting was suppressed, and the digestion-resistant starch (RS) content was increased. Comprehensive microarray polymer profiling (CoMPP) analysis revealed specific alterations of major pectic and glycoprotein cell wall components. This complex phenotype led us to search for candidate IAm genes exploiting its recessive trait. Hence, we sequenced genomic DNA of a pool of IAm lines, identified SNPs genome wide against the draft genome sequence of potato and searched for regions of decreased heterozygosity. Three regions, located on chromosomes 3, 7 and 10, respectively, displayed markedly less heterozygosity than average. The only credible starch metabolism-related gene found in these regions encoded the isoamylase-type debranching enzyme Stisa1. Decreased expression of mRNA (>500 fold) and reduced enzyme activity (virtually absent from IAm lines) supported Stisa1 as a candidate gene for IAm.

Exploring the chemotypes underlying important agronomic and consumer traits in cassava (Manihot esculenta crantz)

A broad diversity of phenotypes are available within the cassava germplasm collections. The phenotypes include improved nutritional, starch or culinary root quality as well as abiotic and biotic resistance properties. Some of these traits can be found naturally occurring in cassava landraces, whereas others are the result of targeted breeding efforts. For future breeding programmes it is important to know the underlying mechanisms of these desirable traits. Metabolomics can assist in the elucidation of these mechanisms by measuring the end products of the cellular processes conferring the traits of interest. The present study focused on the comparison of two or more variants of the same trait such as high and low culinary quality or resistance and susceptibility to thrips. Overall, eight different traits were assessed. Results showed that amino acids and umami compounds were associated with superior culinary attributes and the phenylpropanoid superpathway plays an important role in pest resistance. Furthermore, the data highlighted a low chemodiversity in African cassavas and that the source-sink relation was still active at the harvest stage.

Identification and impact of stable prognostic biochemical markers for cold-induced sweetening resistance on selection efficiency in potato (Solanum tuberosum L.) breeding programs

Biochemical markers for cold-induced sweetening (CIS) resistance were tested for their stability over years and their use in selection of parents for crossing to achieve high selection efficiency in potato breeding programs.

Two regulatory enzymes directly associated with reducing sugar (RS) accumulation during potato tubers cold storage were tested as a predictor for CIS resistance. These enzymes were studied in 33 potato clones from various breeding programs over four years. Clones with the presence of A-II isozymes of UDP-glucose pyrophosphorylase (UGPase) and low activity of vacuolar acid invertase (VAcInv) enzyme had increased resistance to cold-induced sweetening (CIS). Depending on the levels of these enzymes, clones were divided into class A, class B and class C. Clones categorized as class A had average RS of 0.73 mg per g FW after six months at 5.5°C storage. Class B and C had average RS of 1.15 and 3.80 mg per g FW respectively. The enzyme activity was closely associated with RS accumulation over long-term cold storage.

The biochemical markers were found to be stable over the years. Repeated-measure analysis showed 75% chance of maintaining class from one year to the next and a 25% chance of switching, No clone switched between class A and class C, even across all four years. Application of these biochemical markers can identify clones with CIS resistance early in the selection process. Biochemical markers were used to select parents for crossing and six families were established. Results showed that with both parents from class A, 95% of their offspring had desirable glucose levels and chip color, which dropped to 52% when one parent was from class A and other from class B. These results suggest that two regulatory enzymes, i.e., UGPase and VAcInv, can be used as stable prognostic biochemical markers for CIS resistance for precise parent selection resulting in progenies with significantly higher percentage of clones with acceptable processing quality.

The expression pattern of the Pho1a genes encoding plastidic starch phosphorylase correlates with the degradation of starch during fruit ripening in green-fruited and red-fruited tomato species

Genes encoding plastidic starch phosphorylase Pho1a were identified in 10 tomato species (Solanum section Lycopersicon). Pho1a genes showed higher variability in green-fruited than in red-fruited tomato species, but had an extremely low polymorphism level compared with other carbohydrate metabolism genes and an unusually low ratio of intron to exon single nucleotide polymorphisms (SNPs). In red-fruited species, Pho1a was expressed in all analysed tissues, including fruit at different developmental stages, with the highest level in mature green fruit, which is strong sink organ importing sucrose and accumulating starch. In green-fruited species Solanum peruvianum and Solanum arcanum, the Pho1a expression level was similar in mature green and ripe fruit, whereas in Solanum chmielewskii, it was higher in ripe fruit, and in Solanum habrochaites, the dynamics of fruit-specific Pho1a expression was similar to that in red-fruited tomatoes. During fruit development, in red-fruited Solanum lycopersicum, sucrose level was low, the monosaccharide content increased; in green-fruited S. peruvianum, the sucrose concentration increased and those of monosaccharides decreased. In both species, the starch content and Pho1a expression were downregulated. The evolutionary topology based on Pho1a sequences was consistent with the current division of tomatoes into red-fruited and green-fruited species, except for S. habrochaites.

Quantitative trait loci for starch-corrected chip color after harvest, cold storage and after reconditioning mapped in diploid potato

The objective of this study was to map the quantitative trait loci (QTLs) for chip color after harvest (AH), cold storage (CS) and after reconditioning (RC) in diploid potato and compare them with QTLs for starch-corrected chip color. Chip color traits AH, CS, and RC significantly correlated with tuber starch content (TSC). To limit the effect of starch content, the chip color was corrected for TSC. The QTLs for chip color (AH, CS, and RC) and the starch-corrected chip color determined with the starch content after harvest (SCAH), after cold storage (SCCS) and after reconditioning (SCRC) were compared to assess the extent of the effect of starch and the location of genetic factors underlying this effect on chip color. We detected QTLs for the AH, CS, RC and starch-corrected traits on ten potato chromosomes, confirming the polygenic nature of the traits. The QTLs with the strongest effects were detected on chromosomes I (AH, 0 cM, 11.5% of variance explained), IV (CS, 43.9 cM, 12.7%) and I (RC, 49.7 cM, 14.1%). When starch correction was applied, the QTLs with the strongest effects were revealed on chromosomes VIII (SCAH, 39.3 cM, 10.8% of variance explained), XI (SCCS, 79.5 cM, 10.9%) and IV (SCRC, 43.9 cM, 10.8%). Applying the starch correction changed the landscape of QTLs for chip color, as some QTLs became statistically insignificant, shifted or were refined, and new QTLs were detected for SCAH. The QTLs on chromosomes I and IV were significant for all traits with and without starch correction.

Distinguishing CPT gene family members and vetting the sequence structure of a putative rubber synthesizing variant in Hevea brasiliensis

cis-Prenyltransferases (cis-PTs) constitute a large family of enzymes conserved during evolution and present in all domains of life. cis-PTs catalyze the cis-1,4-polymerization of isoprene units to generate isoprenoids with carbon skeletons varying from C10 (neryl pyrophosphate) to C > 10,000 (natural rubber). Though the previously reported CPTs in Hevea are designated based on sequence variations, their classification was done mostly by phylogenetic analysis using a mixture of partial as well as full length sequences often excluding the UTRs. In this context an attempt was made to reclassify the CPTs strictly based on their sequence similarity and distinguish the members putatively associated with rubber biosynthesis from the others. Extensive computational analysis was carried out on CPT sequences obtained from public resources and whole genome assemblies of Hevea. Based on the results from BLAST analysis, multiple sequence alignments of protein, nucleotide and untranslated regions, open reading frame analysis, gene prediction analysis and sequence length variations, we conclude that there exists mainly three CPTs namely RubCPT1, RubCPT2 and RubCPT3 putatively associated with rubber biosynthesis in Hevea brasiliensis. The rest were categorised as variants of dehydrodolichyl diphosphate synthase (DHDDS) involved in the synthesis of dolichols having short chain isoprenoids. Analysis of the sequence structure of the most highly expressed RubCPT1 in latex revealed the allele richness and diversity of this important variant prevailing in the popular rubber clones. Haplotypes consisting of SNPs with high degree of heterozygosity were also identified. Segregation and linkage disequilibrium analysis confirmed that recombination is the major contributor towards the generation of allelic diversity rather than point mutations. Alternatively, gene expression analysis indicated the possibility of association between specific haplotypes and RubCPT1 expression in Hevea clones which may have downstream impact up to the level of rubber production. The conclusions from this study may pave way for the identification and better understanding of CPTs directly involved with natural rubber biosynthesis in Hevea and the SNP data generated may aid in the development of molecular markers putatively associated with yield in rubber.

Comparative Proteomic Analysis Reveals That Antioxidant System and Soluble Sugar Metabolism Contribute to Salt Tolerance in Alfalfa ( Medicago sativa L.) Leaves

Soil salinity poses a serious threat to alfalfa ( Medicago sativa L.) productivity. To characterize the molecular mechanisms of salinity tolerance in Medicago, the comparative proteome of leaves from Medicago sativa cv. Zhongmu No.1 (ZM1, salt-tolerant) and Medicago truncatula cv. Jemalong A17 (A17, salt-sensitive) was performed using the iTRAQ approach. A total of 438 differentially expressed proteins (DEPs) were identified, among which 282 and 120 DEPs were specific to A17 and ZM1, respectively. In salt-tolerant ZM1, key DEPs were primarily enriched in antioxidant system, starch and sucrose metabolism, and secondary metabolism. ZM1 possessed a greater ability to remove reactive oxygen species and methylglyoxal under salt stress, as demonstrated by enhancement of the antioxidant system and secondary metabolism. Moreover, ZM1 orchestrated starch and sucrose metabolism to accumulate various soluble sugars (sucrose, maltose, glucose, and trehalose), which in turn facilitate osmotic homeostasis. Salt stress dramatically inhibited photosynthesis of A17 due to the downregulation of the light-harvesting complex and photosystem II related protein. Quantitative analyses of photochemical efficiency, antioxidant enzyme activities, hydrogen peroxide, malondialdehyde, and soluble sugar contents were consistent with the alterations predicted on the basis of DEP functions. These results shed light on our understanding of the mechanisms underlying the salt tolerance of alfalfa.

Linkage analysis and QTL mapping in a tetraploid russet mapping population of potato

Background

Genome-wide single nucleotide polymorphism (SNP) markers coupled with allele dosage information has emerged as a powerful tool for studying complex traits in cultivated autotetraploid potato (Solanum tuberosum L., 2n = 4× = 48). To date, this approach has been effectively applied to the identification of quantitative trait loci (QTLs) underlying highly heritable traits such as disease resistance, but largely unexplored for traits with complex patterns of inheritance.

Results

In this study, an F₁ tetraploid russet mapping population (162 individuals) was evaluated for multiple quantitative traits over two years and two locations to identify QTLs associated with tuber sugar concentration, processing quality, vine maturity, and other high-value agronomic traits. We report the linkage maps for the 12 potato chromosomes and the QTL location with corresponding genetic models and candidate SNPs explaining the highest phenotypic variation for tuber quality and maturity related traits. Significant QTLs for tuber glucose concentration and tuber fry color were detected on chromosomes 4, 5, 6, 10, and 11. Collectively, these QTLs explained between 24 and 46% of the total phenotypic variation for tuber glucose and fry color, respectively. The QTL on chromosome 10 was associated with apoplastic invertases, with ‘Premier Russet’ contributing the favorable allele for fry processing quality. On chromosome 5, minor-effect QTLs for tuber glucose concentration and fry color co-localized with various major-effect QTLs, including vine maturity, growth habit, tuber shape, early blight (Altenaria tenuis), and Verticillium wilt (Verticillium spp.).

Conclusions

Linkage analysis and QTL mapping in a russet mapping population (A05141) using SNP dosage information successfully identified favorable alleles and candidate SNPs for resistance to the accumulation of tuber reducing sugars. These novel markers have a high potential for the improvement of tuber processing quality. Moreover, the discovery of different genetic models for traits with overlapping QTLs at the maturity locus clearly suggests an independent genetic control.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0672-1) contains supplementary material, which is available to authorized users.

Genomic prediction of starch content and chipping quality in tetraploid potato using genotyping-by-sequencing

KEY MESSAGE

Genomic prediction models for starch content and chipping quality show promising results, suggesting that genomic selection is a feasible breeding strategy in tetraploid potato. Genomic selection uses genome-wide molecular markers to predict performance of individuals and allows selections in the absence of direct phenotyping. It is regarded as a useful tool to accelerate genetic gain in breeding programs, and is becoming increasingly viable for crops as genotyping costs continue to fall. In this study, we have generated genomic prediction models for starch content and chipping quality in tetraploid potato to facilitate varietal development. Chipping quality was evaluated as the colour of a potato chip after frying following cold induced sweetening. We used genotyping-by-sequencing to genotype 762 offspring, derived from a population generated from biparental crosses of 18 tetraploid parents. Additionally, 74 breeding clones were genotyped, representing a test panel for model validation. We generated genomic prediction models from 171,859 single-nucleotide polymorphisms to calculate genomic estimated breeding values. Cross-validated prediction correlations of 0.56 and 0.73 were obtained within the training population for starch content and chipping quality, respectively, while correlations were lower when predicting performance in the test panel, at 0.30-0.31 and 0.42-0.43, respectively. Predictions in the test panel were slightly improved when including representatives from the test panel in the training population but worsened when preceded by marker selection. Our results suggest that genomic prediction is feasible, however, the extremely high allelic diversity of tetraploid potato necessitates large training populations to efficiently capture the genetic diversity of elite potato germplasm and enable accurate prediction across the entire spectrum of elite potatoes. Nonetheless, our results demonstrate that GS is a promising breeding strategy for tetraploid potato.

Physical mapping of QTL for tuber yield, starch content and starch yield in tetraploid potato (Solanum tuberosum L.) by means of genome wide genotyping by sequencing and the 8.3 K SolCAP SNP array

Background

Tuber yield and starch content of the cultivated potato are complex traits of decisive importance for breeding improved varieties. Natural variation of tuber yield and starch content depends on the environment and on multiple, mostly unknown genetic factors. Dissection and molecular identification of the genes and their natural allelic variants controlling these complex traits will lead to the development of diagnostic DNA-based markers, by which precision and efficiency of selection can be increased (precision breeding).

Results

Three case-control populations were assembled from tetraploid potato cultivars based on maximizing the differences between high and low tuber yield (TY), starch content (TSC) and starch yield (TSY, arithmetic product of TY and TSC). The case-control populations were genotyped by restriction-site associated DNA sequencing (RADseq) and the 8.3 k SolCAP SNP genotyping array. The allele frequencies of single nucleotide polymorphisms (SNPs) were compared between cases and controls. RADseq identified, depending on data filtering criteria, between 6664 and 450 genes with one or more differential SNPs for one, two or all three traits. Differential SNPs in 275 genes were detected using the SolCAP array. A genome wide association study using the SolCAP array on an independent, unselected population identified SNPs associated with tuber starch content in 117 genes. Physical mapping of the genes containing differential or associated SNPs, and comparisons between the two genome wide genotyping methods and two different populations identified genome segments on all twelve potato chromosomes harboring one or more quantitative trait loci (QTL) for TY, TSC and TSY.

Conclusions

Several hundred genes control tuber yield and starch content in potato. They are unequally distributed on all potato chromosomes, forming clusters between 0.5–4 Mbp width. The largest fraction of these genes had unknown function, followed by genes with putative signalling and regulatory functions. The genetic control of tuber yield and starch content is interlinked. Most differential SNPs affecting both traits had antagonistic effects: The allele increasing TY decreased TSC and vice versa. Exceptions were 89 SNP alleles which had synergistic effects on TY, TSC and TSY. These and the corresponding genes are primary targets for developing diagnostic markers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3979-9) contains supplementary material, which is available to authorized users.

Quantitative trait loci from identification to exploitation for crop improvement

Advancement in the field of genetics and genomics after the discovery of Mendel's laws of inheritance has led to map the genes controlling qualitative and quantitative traits in crop plant species. Mapping of genomic regions controlling the variation of quantitatively inherited traits has become routine after the advent of different types of molecular markers. Recently, the next generation sequencing methods have accelerated the research on QTL analysis. These efforts have led to the identification of more closely linked molecular markers with gene/QTLs and also identified markers even within gene/QTL controlling the trait of interest. Efforts have also been made towards cloning gene/QTLs or identification of potential candidate genes responsible for a trait. Further new concepts like crop QTLome and QTL prioritization have accelerated precise application of QTLs for genetic improvement of complex traits. In the past years, efforts have also been made in exploitation of a number of QTL for improving grain yield or other agronomic traits in various crops through markers assisted selection leading to cultivation of these improved varieties at farmers' field. In present article, we reviewed QTLs from their identification to exploitation in plant breeding programs and also reviewed that how improved cultivars developed through introgression of QTLs have improved the yield productivity in many crops.

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

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_45711901T and InvGE-C_2475454T were associated with sucrose content and frying color, respectively, and were not found previously in tetraploid genotypes. The rare haplotype InvGE-A_2475187C_2475295A_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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0489-3) contains supplementary material, which is available to authorized users.

Assuring Potato Tuber Quality during Storage: A Future Perspective

Potatoes represent an important staple food crop across the planet. Yet, to maintain tuber quality and extend availability, there is a necessity to store tubers for long periods often using industrial-scale facilities. In this context, preserving potato quality is pivotal for the seed, fresh and processing sectors. The industry has always innovated and invested in improved post-harvest storage. However, the pace of technological change has and will continue to increase. For instance, more stringent legislation and changing consumer attitudes have driven renewed interest in creating alternative or complementary post-harvest treatments to traditional chemically reliant sprout suppression and disease control. Herein, the current knowledge on biochemical factors governing dormancy, the use of chlorpropham (CIPC) as well as existing and chemical alternatives, and the effects of pre- and post-harvest factors to assure potato tuber quality is reviewed. Additionally, the role of genomics as a future approach to potato quality improvement is discussed. Critically, and through a more industry targeted research, a better mechanistic understanding of how the pre-harvest environment influences tuber quality and the factors which govern dormancy transition should lead to a paradigm shift in how sustainable storage can be achieved.

Comparative Transcriptome Analysis Reveals Critical Function of Sucrose Metabolism Related-Enzymes in Starch Accumulation in the Storage Root of Sweet Potato

The starch properties of the storage root (SR) affect the quality of sweet potato (Ipomoea batatas (L.) Lam.). Although numerous studies have analyzed the accumulation and properties of starch in sweet potato SRs, the transcriptomic variation associated with starch properties in SR has not been quantified. In this study, we measured the starch and sugar contents and analyzed the transcriptome profiles of SRs harvested from sweet potatoes with high, medium, and extremely low starch contents, at five developmental stages [65, 80, 95, 110, and 125 days after transplanting (DAP)]. We found that differences in both water content and starch accumulation in the dry matter affect the starch content of SRs in different sweet potato genotypes. Based on transcriptome sequencing data, we assembled 112336 unigenes, and identified several differentially expressed genes (DEGs) involved in starch and sucrose metabolism, and revealed the transcriptional regulatory network controlling starch and sucrose metabolism in sweet potato SRs. Correlation analysis between expression patterns and starch and sugar contents suggested that the sugar-starch conversion steps catalyzed by sucrose synthase (SuSy) and UDP-glucose pyrophosphorylase (UGPase) may be essential for starch accumulation in the dry matter of SRs, and IbβFRUCT2, a vacuolar acid invertase, might also be a key regulator of starch content in the SRs. Our results provide valuable resources for future investigations aimed at deciphering the molecular mechanisms determining the starch properties of sweet potato SRs.

Natural variation of sucrose, glucose and fructose contents in Colombian genotypes of Solanum tuberosum Group Phureja at harvest

BACKGROUND

Potato frying quality is a complex trait influenced by sugar content in tubers. Good frying quality requires low content of reducing sugars to avoid the formation of dark pigments. Solanum tuberosum Group Phureja is a valuable genetic resource for breeding and for genetic studies. The sugar content after harvest was analyzed in a germplasm collection of Group Phureja to contribute to the understanding of the natural variation of this trait.

RESULTS

Sucrose, glucose and fructose genotypic mean values ranged from 6.39 to 29.48 g kg(-1) tuber dry weight (DW), from 0.46 to 28.04 g kg(-1) tuber DW and from 0.29 to 27.23 g kg(-1) tuber DW, respectively. Glucose/fructose and sucrose/reducing sugars ratios ranged from 1.01 to 6.67 mol mol(-1) and from 0.15 to 7.78 mol mol(-1) , respectively. Five clusters of genotypes were recognized, three of them with few genotypes and extreme phenotypic values.

CONCLUSION

Sugar content showed a wide variation, representing the available variability useful for potato breeding. The results provide a quantitative approach to analyze the frying quality trait and are consistent with frying color. The analyzed germplasm presents extreme phenotypes, which will contribute to the understanding of the genetic basis of this trait. © 2016 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Identification and reproducibility of diagnostic DNA markers for tuber starch and yield optimization in a novel association mapping population of potato (Solanum tuberosum L.)

SNPs in candidate genes Pain - 1, InvCD141 (invertases), SSIV (starch synthase), StCDF1 (transcription factor), LapN (leucine aminopeptidase), and cytoplasm type are associated with potato tuber yield, starch content and/or starch yield. Tuber yield (TY), starch content (TSC), and starch yield (TSY) are complex characters of high importance for the potato crop in general and for industrial starch production in particular. DNA markers associated with superior alleles of genes that control the natural variation of TY, TSC, and TSY could increase precision and speed of breeding new cultivars optimized for potato starch production. Diagnostic DNA markers are identified by association mapping in populations of tetraploid potato varieties and advanced breeding clones. A novel association mapping population of 282 genotypes including varieties, breeding clones and Andean landraces was assembled and field evaluated in Northern Spain for TY, TSC, TSY, tuber number (TN) and tuber weight (TW). The landraces had lower mean values of TY, TW, TN, and TSY. The population was genotyped for 183 microsatellite alleles, 221 single nucleotide polymorphisms (SNPs) in fourteen candidate genes and eight known diagnostic markers for TSC and TSY. Association test statistics including kinship and population structure reproduced five known marker-trait associations of candidate genes and discovered new ones, particularly for tuber yield and starch yield. The inclusion of landraces increased the number of detected marker-trait associations. Integration of the present association mapping results with previous QTL linkage mapping studies for TY, TSC, TSY, TW, TN, and tuberization revealed some hot spots of QTL for these traits in the potato genome. The genomic positions of markers linked or associated with QTL for complex tuber traits suggest high multiplicity and genome wide distribution of the underlying genes.

Mapping of quantitative trait loci for tuber starch and leaf sucrose contents in diploid potato

Most QTL for leaf sucrose content map to positions that are similar to positions of QTL for tuber starch content in diploid potato. In the present study, using a diploid potato mapping population and Diversity Array Technology (DArT) markers, we identified twelve quantitative trait loci (QTL) for tuber starch content on seven potato chromosomes: I, II, III, VIII, X, XI, and XII. The most important QTL spanned a wide region of chromosome I (42.0–104.6 cM) with peaks at 63 and 84 cM which explained 17.6 and 19.2% of the phenotypic variation, respectively. ADP-glucose pyrophosphorylase (AGPase) is the key enzyme for starch biosynthesis. The gene encoding the large subunit of this enzyme, AGPaseS-a, was localized to chromosome I at 102.3 cM and accounted for 15.2% of the variance in tuber starch content. A more than 100-fold higher expression of this gene was observed in RT-qPCR assay in plants with the marker allele AGPaseS-a1334. This study is the first to report QTL for sucrose content in potato leaves. QTL for sucrose content in leaves were located on eight potato chromosomes: I, II, III, V, VIII, IX, X and XII. In 5-week-old plants, only one QTL for leaf sucrose content was detected after 8 h of darkness; four QTL were detected after 8 h of illumination. In 11-week-old plants, 6 and 3 QTL were identified after dark and light phases, respectively. Of fourteen QTL for leaf sucrose content, eleven mapped to positions that were similar to QTL for tuber starch content. These results provide genetic information for further research examining the relationships between metabolic carbon molecule sources and sinks in potato plants.

Development and Identification of SSR Markers Associated with Starch Properties and β-Carotene Content in the Storage Root of Sweet Potato (Ipomoea batatas L.)

Sweet potato (Ipomoea batatas L.) is a nutritious food crop and, based on the high starch content of its storage root, a potential bioethanol feedstock. Enhancing the nutritional value and starch quantity of storage roots are important goals of sweet potato breeding programs aimed at developing improved varieties for direct consumption, processing, and industrial uses. However, developing improved lines of sweet potato is challenging due to the genetic complexity of this plant and the lack of genome information. Short sequence repeat (SSR) markers are powerful molecular tools for tracking important loci in crops and for molecular-based breeding strategies; however, few SSR markers and marker-trait associations have hitherto been identified in sweet potato. In this study, we identified 1824 SSRs by using a de novo assembly of publicly available ESTs and mRNAs in sweet potato, and designed 1476 primer pairs based on SSR-containing sequences. We mapped 214 pairs of primers in a natural population comprised of 239 germplasms, and identified 1278 alleles with an average of 5.972 alleles per locus and a major allele frequency of 0.7702. Population structure analysis revealed two subpopulations in this panel of germplasms, and phenotypic characterization demonstrated that this panel is suitable for association mapping of starch-related traits. We identified 32, 16, and 17 SSR markers associated with starch content, β-carotene content, and starch composition in the storage root, respectively, using association analysis and further evaluation of a subset of sweet potato genotypes with various characteristics. The SSR markers identified here can be used to select varieties with desired traits and to investigate the genetic mechanism underlying starch and carotenoid formation in the starchy roots of sweet potato.

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

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.

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

Starch phosphorylation is an important aspect of plant metabolism due to its role in starch degradation. Moreover, the degree of phosphorylation of starch determines its physicochemical properties and is therefore relevant for industrial uses of starch. Currently, starch is chemically phosphorylated to increase viscosity and paste stability. Potato cultivars with elevated starch phosphorylation would make this process unnecessary, thereby bestowing economic and environmental benefits. Starch phosphorylation is a complex trait which has been previously shown by antisense gene repression to be influenced by a number of genes including those involved in starch synthesis and degradation. We have used an association mapping approach to discover genetic markers associated with the degree of starch phosphorylation. A diverse collection of 193 potato lines was grown in replicated field trials, and the levels of starch phosphorylation at the C6 and C3 positions of the glucosyl residues were determined by mass spectrometry of hydrolyzed starch from tubers. In addition, the potato lines were genotyped by amplicon sequencing and microsatellite analysis, focusing on candidate genes known to be involved in starch synthesis. As potato is an autotetraploid, genotyping included determination of allele dosage. Significant associations (p < 0.001) were found with SNPs in the glucan water dikinase (GWD), starch branching enzyme I (SBEI) and the starch synthase III (SSIII) genes, and with a SSR allele in the SBEII gene. SNPs in the GWD gene were associated with C6 phosphorylation, whereas polymorphisms in the SBEI and SBEII genes were associated with both C6 and C3 phosphorylation and the SNP in the SSIII gene was associated with C3 phosphorylation. These allelic variants have potential as genetic markers for starch phosphorylation in potato.

Novel candidate genes AuxRP and Hsp90 influence the chip color of potato tubers

Potato (Solanum tuberosum L.) tubers exhibit significant variation in reducing sugar content directly after harvest, cold storage and reconditioning. Here, we performed QTL analysis for chip color, which is strongly influenced by reducing sugar content, in a diploid potato mapping population. Two QTL on chromosomes I and VI were detected for chip color after harvest and reconditioning. Only one region on chromosome VI was linked with cold-induced sweetening. Using the RT-PCR technique, we showed differential expression of the auxin-regulated protein (AuxRP) gene. The AuxRP transcript was presented in light chip color parental clone DG 97-952 and the RNA progeny of the bulk sample consisting of light chip color phenotypes after cold storage. This amplicon was absent in dark chip parental clone DG 08-26/39 and the RNA bulk sample of dark chip progeny. Genetic variation of AuxRP explained up to 16.6 and 15.2 % of the phenotypic variance after harvest and 3 months of storage at 4 °C, respectively. Using an alternative approach, the RDA-cDNA method was used to recognize 25 gene sequences, of which 11 could be assigned to potato chromosome VI. One of these genes, Heat-shock protein 90 (Hsp90), demonstrated higher mRNA and protein expression in RT-qPCR and western blotting assays in the dark chip color progeny bulk sample compared with the light chip color progeny bulk sample. Our study, for the first time, suggests that the AuxRP and Hsp90 genes are novel candidate genes capable of influencing the chip color of potato tubers.