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

Genotyping-by-sequencing and weighted gene co-expression network analysis of genes responsive against Potato virus Y in commercial potato cultivars

Potato is considered a key component of the global food system and plays a vital role in strengthening world food security. A major constraint to potato production worldwide is the Potato Virus Y (PVY), belonging to the genus Potyvirus in the family of Potyviridae. Selective breeding of potato with resistance to PVY pathogens remains the best method to limit the impact of viral infections. Understanding the genetic diversity and population structure of potato germplasm is important for breeders to improve new cultivars for the sustainable use of genetic materials in potato breeding to PVY pathogens. While, genetic diversity improvement in modern potato breeding is facing increasingly narrow genetic basis and the decline of the genetic diversity. In this research, we performed genotyping-by-sequencing (GBS)-based diversity analysis on 10 commercial potato cultivars and weighted gene co-expression network analysis (WGCNA) to identify candidate genes related to PVY-resistance. WGCNA is a system biology technique that uses the WGCNA R software package to describe the correlation patterns between genes in multiple samples. In terms of consumption, these cultivars are a high rate among Iranian people. Using population structure analysis, the 10 cultivars were clustered into three groups based on the 118343 single nucleotide polymorphisms (SNPs) generated by GBS. Read depth ranged between 5 and 18. The average data size and Q30 of the reads were 145.98 Mb and 93.63%, respectively. Based on the WGCNA and gene expression analysis, the StDUF538, StGTF3C5, and StTMEM161A genes were associated with PVY resistance in the potato genome. Further, these three hub genes were significantly involved in defense mechanism where the StTMEM161A was involved in the regulation of alkalization apoplast, the StDUF538 was activated in the chloroplast degradation program, and the StGTF3C5 regulated the proteins increase related to defense in the PVY infected cells. In addition, in the genetic improvement programs, these hub genes can be used as genetic markers for screening commercial cultivars for PVY resistance. Our survey demonstrated that the combination of GBS-based genetic diversity germplasm analysis and WGCNA can assist breeders to select cultivars resistant to PVY as well as help design proper crossing schemes in potato breeding.

Identification of QTL associated with plant vine characteristics and infection response to late blight, early blight, and Verticillium wilt in a tetraploid potato population derived from late blight-resistant Palisade Russet

Potato late blight (causal agent Phytophthora infestans) is a disease of potatoes with economic importance worldwide. Control is primarily through field monitoring and the application of fungicides. Control of late blight with fungicides and host plant resistance is difficult, with documented cases of such control measures failing with the advent of new pathotypes of P. infestans. To better understand host plant resistance and to develop more durable late blight resistance, Quantitative Trait Locus/Loci (QTL) analysis was conducted on a tetraploid mapping population derived from late blight-resistant potato cultivar Palisade Russet. Additionally, QTL analyses for other traits such as Verticillium wilt and early blight resistance, vine size and maturity were performed to identify a potential relationship between multiple traits and prepare genetic resources for molecular markers useful in breeding programs. For this, one hundred ninety progenies from intercrossing Palisade Russet with a late blight susceptible breeding clone (ND028673B-2Russ) were assessed. Two parents and progenies were evaluated over a two-year period for response to infection by the US-8 genotype of P. infestans in inoculated field screenings in Corvallis, Oregon. In Aberdeen, Idaho, the same mapping population was also evaluated for phenotypic response to early blight and Verticillium wilt, and vine size and maturity in a field over a two-year period. After conducting QTL analyses with those collected phenotype data, it was observed that chromosome 5 has a significant QTL for all five traits. Verticillium wilt and vine maturity QTL were also observed on chromosome 1, and vine size QTL was also found on chromosomes 3 and 10. An early blight QTL was also detected on chromosome 2. The QTL identified in this study have the potential for converting into breeder-friendly molecular markers for marker-assisted selection.

Potato tonoplast sugar transporter 1 controls tuber sugar accumulation during postharvest cold storage

Cold-induced sweetening (CIS), the undesirable sugar accumulation in cold-stored potato (Solanum tuberosum L.) tubers, is a severe postharvest issue in the potato processing industry. Although the process of sucrose hydrolysis by vacuolar invertase during potato CIS is well understood, there is limited knowledge about the transportation of sucrose from the cytosol to the vacuole during postharvest cold storage. Here, we report that among the three potato tonoplast sugar transporters (TSTs), StTST1 exhibits the highest expression in tubers during postharvest cold storage. Subcellular localization analysis demonstrates that StTST1 is a tonoplast-localized protein. StTST1 knockdown decreases reducing sugar accumulation in tubers during low-temperature storage. Compared to wild-type, potato chips produced from StTST1-silenced tubers displayed significantly lower acrylamide levels and lighter color after cold storage. Transcriptome analysis manifests that suppression of StTST1 promotes starch synthesis and inhibits starch degradation in cold-stored tubers. We further establish that the increased sucrose content in the StTST1-silenced tubers might cause a decrease in the ABA content, thereby inhibiting the ABA-signaling pathway. We demonstrate that the down-regulation of β-amylase StBAM1 in StTST1-silenced tubers might be directly controlled by ABA-responsive element-binding proteins (AREBs). Altogether, we have shown that StTST1 plays a critical role in sugar accumulation and starch metabolism regulation during postharvest cold storage. Thus, our findings provide a new strategy to improve the frying quality of cold-stored tubers and reduce the acrylamide content in potato chips.

Identification of Long-Distance Transport Signal Molecules Associated with Plant Maturity in Tetraploid Cultivated Potatoes (Solanum tuberosum L.)

Maturity is a key trait for breeders to identify potato cultivars suitable to grow in different latitudes. However, the molecular mechanism regulating maturity remains unclear. In this study, we performed a grafting experiment using the early-maturing cultivar Zhongshu 5 (Z5) and the late-maturing cultivar Zhongshu 18 (Z18) and found that abscisic acid (ABA) and salicylic acid (SA) positively regulate the early maturity of potato, while indole-3-acetic acid (IAA) negatively regulated early maturity. A total of 43 long-distance transport mRNAs are observed to be involved in early maturity, and 292 long-distance transport mRNAs involved in late maturity were identified using RNA sequencing. Specifically, StMADS18, StSWEET10C, and StSWEET11 are detected to be candidate genes for their association with potato early maturity. Metabolomic data analysis shows a significant increase in phenolic acid and flavonoid contents increased in the scion of the early-maturing cultivar Z5, but a significant decrease in amino acid, phenolic acid, and alkaloid contents increased in the scion of the late-maturing cultivar Z18. This work reveals a significant association between the maturity of tetraploid cultivated potato and long-distance transport signal molecules and provides useful data for assessing the molecular mechanisms underlying the maturity of potato plants and for breeding early-maturing potato cultivars.

Suppression of the tonoplast sugar transporter StTST3.2 improves quality of potato chips

Which sugar transporter regulates sugar accumulation in tubers is largely unknown. Accumulation of reducing sugar (RS) in potato (Solanum tuberosum L.) tubers negatively affects the quality of tubers undergoing the frying process. However, little is known about the genes involved in regulating RS content in tubers at harvest. Here, we have identified two tonoplast sugar transporter (TST) 3-type isoforms (StTST3.1 and StTST3.2) in potato. Quantitative real-time PCR results indicate that StTST3.1 and StTST3.2 possess distinct expression patterns in various potato tissues. StTST3.2 was found to be the expressed TST3-type isoform in tubers. Further subcellular localization analysis revealed that StTST3.2 was targeted to the tonoplast. Silencing of StTST3.2 in potato by stable transformation resulted in significantly lower RS content in tubers at harvest or after room temperature storage, suggesting StTST3.2 plays an important role in RS accumulation in tubers. Accordingly, compared with the unsilenced control, potato chips processed from StTST3.2-silenced tubers exhibited lighter color and dramatically decreased acrylamide production at harvest or after room temperature storage. In addition, we demonstrated that silencing of StTST3.2 has no significant effect on potato growth and development. Thus, suppression of StTST3.2 could be another effective approach for improving processing quality and decreasing acrylamide content in potato tubers.

Linkage and QTL mapping for tuber shape and specific gravity in a tetraploid mapping population of potato representing the russet market class

BACKGROUND

Tuber shape and specific gravity (dry matter) are important agronomic traits in potato processing and impact production costs, quality, and consistency of the final processed food products such as French fries and potato chips. In this study, linkage and QTL mapping were performed for these two traits to allow for the implementation of marker-assisted selection to facilitate breeding efforts in the russet market class. Two parents, Rio Grande Russet (female) and Premier Russet (male) and their 205 F1 progenies were initially phenotyped for tuber shape and specific gravity in field trials conducted in Idaho and North Carolina in 2010 and 2011, with specific gravity also being measured in Minnesota in 2011. Progenies and parents were previously genotyped using the Illumina SolCAP Infinium 8303 Potato SNP array, with ClusterCall and MAPpoly (R-packages) subsequently used for autotetraploid SNP calling and linkage mapping in this study. The 12 complete linkage groups and phenotypic data were then imported into QTLpoly, an R-package designed for polyploid QTL analyses.

RESULTS

Significant QTL for tuber shape were detected on chromosomes 4, 7, and 10, with heritability estimates ranging from 0.09 to 0.36. Significant tuber shape QTL on chromosomes 4 and 7 were specific to Idaho and North Carolina environments, respectively, whereas the QTL on chromosome 10 was significant regardless of growing environment. Single marker analyses identified alleles in the parents associated with QTL on chromosomes 4, 7, and 10 that contributed to significant differences in tuber shape among progenies. Significant QTL were also identified for specific gravity on chromosomes 1 and 5 with heritability ranging from 0.12 to 0.21 and were reflected across environments.

CONCLUSION

Fully automated linkage mapping and QTL analysis were conducted to identify significant QTL for tuber shape and dry matter in a tetraploid mapping population representing the russet market class. The findings are important for the development of molecular markers useful to potato breeders for marker-assisted selection for the long tuber shape and acceptable dry matter required by the potato industry within this important market class.

Genome-Wide Approach to Identify Quantitative Trait Loci for Drought Tolerance in Tetraploid Potato (Solanum tuberosum L.)

Drought represents a major abiotic stress factor negatively affecting growth, yield and tuber quality of potatoes. Quantitative trait locus (QTL) analyses were performed in cultivated potatoes for drought tolerance index DRYM (deviation of relative starch yield from the experimental median), tuber starch content, tuber starch yield, tuber fresh weight, selected transcripts and metabolites under control and drought stress conditions. Eight genomic regions of major interest for drought tolerance were identified, three representing standalone DRYM QTL. Candidate genes, e.g., from signaling pathways for ethylene, abscisic acid and brassinosteroids, and genes encoding cell wall remodeling enzymes were identified within DRYM QTL. Co-localizations of DRYM QTL and QTL for tuber starch content, tuber starch yield and tuber fresh weight with underlying genes of the carbohydrate metabolism were observed. Overlaps of DRYM QTL with metabolite QTL for ribitol or galactinol may indicate trade-offs between starch and compatible solute biosynthesis. Expression QTL confirmed the drought stress relevance of selected transcripts by overlaps with DRYM QTL. Bulked segregant analyses combined with next-generation sequencing (BSAseq) were used to identify mutations in genes under the DRYM QTL on linkage group 3. Future analyses of identified genes for drought tolerance will give a better insight into drought tolerance in potatoes.

The recombination landscape and multiple QTL mapping in a Solanum tuberosum cv. 'Atlantic'-derived F1 population

There are many challenges involved with the genetic analyses of autopolyploid species, such as the tetraploid potato, Solanum tuberosum (2n = 4x = 48). The development of new analytical methods has made it valuable to re-analyze an F₁ population (n = 156) derived from a cross involving 'Atlantic', a widely grown chipping variety in the USA. A fully integrated genetic map with 4285 single nucleotide polymorphisms, spanning 1630 cM, was constructed with MAPpoly software. We observed that bivalent configurations were the most abundant ones (51.0~72.4% depending on parent and linkage group), though multivalent configurations were also observed (2.2~39.2%). Seven traits were evaluated over four years (2006-8 and 2014) and quantitative trait loci (QTL) mapping was carried out using QTLpoly software. Based on a multiple-QTL model approach, we detected 21 QTL for 15 out of 27 trait-year combination phenotypes. A hotspot on linkage group 5 was identified with co-located QTL for maturity, plant yield, specific gravity, and internal heat necrosis resistance evaluated over different years. Additional QTL for specific gravity and dry matter were detected with maturity-corrected phenotypes. Among the genes around QTL peaks, we found those on chromosome 5 that have been previously implicated in maturity (StCDF1) and tuber formation (POTH1). These analyses have the potential to provide insights into the biology and breeding of tetraploid potato and other autopolyploid species.

A tetrasomic inheritance model and likelihood-based method for mapping quantitative trait loci in autotetraploid species

Dissecting the genetic architecture of quantitative traits in autotetraploid species is a methodologically challenging task, but a pivotally important goal for breeding globally important food crops, including potato and blueberry, and ornamental species such as rose. Mapping quantitative trait loci (QTLs) is now a routine practice in diploid species but is far less advanced in autotetraploids, largely due to a lack of analytical methods that account for the complexities of tetrasomic inheritance. We present a novel likelihood-based method for QTL mapping in outbred segregating populations of autotetraploid species. The method accounts properly for sophisticated features of gene segregation and recombination in an autotetraploid meiosis. It may model and analyse molecular marker data with or without allele dosage information, such as that from microarray or sequencing experiments. The method developed outperforms existing bivalent-based methods, which may fail to model and analyse the full spectrum of experimental data, in the statistical power of QTL detection, and accuracy of QTL location, as demonstrated by an intensive simulation study and analysis of data sets collected from a segregating population of potato (Solanum tuberosum). The study enables QTL mapping analysis to be conducted in autotetraploid species under a rigorous tetrasomic inheritance model.

Insight into the bZIP Gene Family in Solanum tuberosum: Genome and Transcriptome Analysis to Understand the Roles of Gene Diversification in Spatiotemporal Gene Expression and Function

The basic region-leucine zipper (bZIP) transcription factors (TFs) form homodimers and heterodimers via the coil–coil region. The bZIP dimerization network influences gene expression across plant development and in response to a range of environmental stresses. The recent release of the most comprehensive potato reference genome was used to identify 80 StbZIP genes and to characterize their gene structure, phylogenetic relationships, and gene expression profiles. The StbZIP genes have undergone 22 segmental and one tandem duplication events. Ka/Ks analysis suggested that most duplications experienced purifying selection. Amino acid sequence alignments and phylogenetic comparisons made with the Arabidopsis bZIP family were used to assign the StbZIP genes to functional groups based on the Arabidopsis orthologs. The patterns of introns and exons were conserved within the assigned functional groups which are supportive of the phylogeny and evidence of a common progenitor. Inspection of the leucine repeat heptads within the bZIP domains identified a pattern of attractive pairs favoring homodimerization, and repulsive pairs favoring heterodimerization. These patterns of attractive and repulsive heptads were similar within each functional group for Arabidopsis and S. tuberosum orthologs. High-throughput RNA-seq data indicated the most highly expressed and repressed genes that might play significant roles in tissue growth and development, abiotic stress response, and response to pathogens including Potato virus X. These data provide useful information for further functional analysis of the StbZIP gene family and their potential applications in crop improvement.

Multiple QTL Mapping in Autopolyploids: A Random-Effect Model Approach with Application in a Hexaploid Sweetpotato Full-Sib Population

In developing countries, the sweetpotato, Ipomoea batatas (L.) Lam. [Formula: see text], is an important autopolyploid species, both socially and economically. However, quantitative trait loci (QTL) mapping has remained limited due to its genetic complexity. Current fixed-effect models can fit only a single QTL and are generally hard to interpret. Here, we report the use of a random-effect model approach to map multiple QTL based on score statistics in a sweetpotato biparental population ('Beauregard' × 'Tanzania') with 315 full-sibs. Phenotypic data were collected for eight yield component traits in six environments in Peru, and jointly adjusted means were obtained using mixed-effect models. An integrated linkage map consisting of 30,684 markers distributed along 15 linkage groups (LGs) was used to obtain the genotype conditional probabilities of putative QTL at every centiMorgan position. Multiple interval mapping was performed using our R package QTLpoly and detected a total of 13 QTL, ranging from none to four QTL per trait, which explained up to 55% of the total variance. Some regions, such as those on LGs 3 and 15, were consistently detected among root number and yield traits, and provided a basis for candidate gene search. In addition, some QTL were found to affect commercial and noncommercial root traits distinctly. Further best linear unbiased predictions were decomposed into additive allele effects and were used to compute multiple QTL-based breeding values for selection. Together with quantitative genotyping and its appropriate usage in linkage analyses, this QTL mapping methodology will facilitate the use of genomic tools in sweetpotato breeding as well as in other autopolyploids.

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.

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.