Genetic mapping with an inbred line-derived F2 population in potato

Tuber transcriptome analysis reveals a novel WRKY transcription factor StWRKY70 potentially involved in potato pigmentation

Tuber flesh pigmentation, conferred by the presence of secondary metabolite anthocyanins, is one of many key agronomic traits for potato tubers. Although several genes of potato anthocyanin biosynthesis have been reported, transcription factors (TFs) contributing to tuber flesh pigmentation are still not fully understood. In this study, transcriptomic profiling of diploid potato accessions with or without tuber flesh pigmentation was conducted and genes of the anthocyanin biosynthesis pathway were found significantly enriched within the 1435 differentially expressed genes (DEGs). Weighted Gene Co-expression Network Analysis (WGCNA) and connectivity analysis pinpointed a subset of 173 genes closely related to the key biosynthetic gene StDFR. Of the eight transcription factors in the subset, group III WRKY StWRKY70, was chosen for showing high connectivity to StDFR and ten other anthocyanin biosynthetic genes and homology to known WRKYs of anthocyanin pathway. The transient activation assay showed StWRKY70 predominantly stimulated the expression of StDFR and StANS as well as the accumulation of anthocyanins by enhancing the function of the MYB transcription factor StAN1. Furthermore, the interaction between StWRKY70 and StAN1 was verified by Y2H and BiFC. Our analysis discovered a new transcriptional activator StWRKY70 which potentially involved in tuber flesh pigmentation, thus may lay the foundation for deciphering how the WRKY-MYB-bHLH-WD40 (WRKY-MBW) complex regulate the accumulation of anthocyanins and provide new strategies to breed for more nutritious potato varieties with enhanced tuber flesh anthocyanins.

Development of near homozygous lines for diploid hybrid TPS breeding in potatoes

Diploid inbred-based F1 hybrid True Potato Seed (DHTPS) breeding is a novel technique to transform potato breeding and cultivation across the globe. Significant efforts are being made to identify elite diploids, dihaploids and develop diploid inbred lines for heterosis exploitation in potatoes. Self-incompatibility is the first obstacle for developing inbred lines in diploid potatoes, which necessitates the introgression of a dominant S locus inhibitor gene (Sli) for switching self-incompatibility to self-compatibility. We evaluated a set of 357 diploid clones in different selfing generations for self-compatibility and degree of homozygosity using Kompetitive Allele Specific PCR (KASP) Single Nucleotide Polymorphism (SNP) markers. A subset of 10 KASP markers of the Sli candidate region on chromosome 12 showed an association with the phenotype for self-compatibility. The results revealed that the selected 10 KASP markers for the Sli gene genotype could be deployed for high throughput rapid screening of self-compatibility in diploid populations and to identify new sources of self-compatibility. The homozygosity assessed through 99 KASP markers distributed across all the chromosomes of the potato genome was 20-78 % in founder diploid clones, while different selfing generations, i.e., S0, S1, S2 and S3 observed 36.1-80.4, 56.9-82.8, 59.5-85.4 and 73.7-87.8 % average homozygosity, respectively. The diploid plants with ∼80 % homozygosity were also observed in the first selfing generation, which inferred that homozygosity assessment in the early generations itself could identify the best plants with high homozygosity to speed up the generation of diploid inbred lines.

Molecular dissection of an intronic enhancer governing cold-induced expression of the vacuolar invertase gene in potato

Potato (Solanum tuberosum) is the third most important food crop in the world. Potato tubers must be stored at cold temperatures to minimize sprouting and losses due to disease. However, cold temperatures strongly induce the expression of the potato vacuolar invertase gene (VInv) and cause reducing sugar accumulation. This process, referred to as "cold-induced sweetening," is a major postharvest problem for the potato industry. We discovered that the cold-induced expression of VInv is controlled by a 200 bp enhancer, VInvIn2En, located in its second intron. We identified several DNA motifs in VInvIn2En that bind transcription factors involved in the plant cold stress response. Mutation of these DNA motifs abolished VInvIn2En function as a transcriptional enhancer. We developed VInvIn2En deletion lines in both diploid and tetraploid potato using clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated nuclease 9 (Cas9)-mediated gene editing. VInv transcription in cold-stored tubers was significantly reduced in the deletion lines. Interestingly, the VInvIn2En sequence is highly conserved among distantly related Solanum species, including tomato (Solanum lycopersicum) and other non-tuber-bearing species. We conclude that the VInv gene and the VInvIn2En enhancer have adopted distinct roles in the cold stress response in tubers of tuber-bearing Solanum species.

A k-mer-based bulked segregant analysis approach to map seed traits in unphased heterozygous potato genomes

High-throughput sequencing-based methods for bulked segregant analysis (BSA) allow for the rapid identification of genetic markers associated with traits of interest. BSA studies have successfully identified qualitative (binary) and quantitative trait loci (QTLs) using QTL mapping. However, most require population structures that fit the models available and a reference genome. Instead, high-throughput short-read sequencing can be combined with BSA of k-mers (BSA-k-mer) to map traits that appear refractory to standard approaches. This method can be applied to any organism and is particularly useful for species with genomes diverged from the closest sequenced genome. It is also instrumental when dealing with highly heterozygous and potentially polyploid genomes without phased haplotype assemblies and for which a single haplotype can control a trait. Finally, it is flexible in terms of population structure. Here, we apply the BSA-k-mer method for the rapid identification of candidate regions related to seed spot and seed size in diploid potato. Using a mixture of F1 and F2 individuals from a cross between 2 highly heterozygous parents, candidate sequences were identified for each trait using the BSA-k-mer approach. Using parental reads, we were able to determine the parental origin of the loci. Finally, we mapped the identified k-mers to a closely related potato genome to validate the method and determine the genomic loci underlying these sequences. The location identified for the seed spot matches with previously identified loci associated with pigmentation in potato. The loci associated with seed size are novel. Both loci are relevant in future breeding toward true seeds in potato.

Linkage mapping of root shape traits in two carrot populations

This study investigated the genetic basis of carrot root shape traits using composite interval mapping in two biparental populations (n = 119 and n = 128). The roots of carrot F2:3 progenies were grown over 2 years and analyzed using a digital imaging pipeline to extract root phenotypes that compose market class. Broad-sense heritability on an entry-mean basis ranged from 0.46 to 0.80 for root traits. Reproducible quantitative trait loci (QTL) were identified on chromosomes 2 and 6 on both populations. Colocalization of QTLs for phenotypically correlated root traits was also observed and coincided with previously identified QTLs in published association and linkage mapping studies. Individual QTLs explained between 14 and 27% of total phenotypic variance across traits, while four QTLs for length-to-width ratio collectively accounted for up to 73% of variation. Predicted genes associated with the OFP-TRM (OVATE Family Proteins-TONNEAU1 Recruiting Motif) and IQD (IQ67 domain) pathway were identified within QTL support intervals. This observation raises the possibility of extending the current regulon model of fruit shape to include carrot storage roots. Nevertheless, the precise molecular mechanisms through which this pathway operates in roots characterized by secondary growth originating from cambium layers remain unknown.

QTL identified that influence tuber length-width ratio, degree of flatness, tuber size, and specific gravity in a russet-skinned, tetraploid mapping population

Potato tuber shape, size, and specific gravity are important agronomic traits in the russet market class of potatoes with an impact on quality, consistency, and product recovery of processed foods such as French fries. Therefore, identifying genetic regions associated with the three traits through quantitative trait locus/loci (QTL) analysis is a crucial process in the subsequent development of marker-assisted selection for use in potato breeding programs. QTL analysis was conducted on a tetraploid mapping population consisting of 190 individuals derived from the cross between two russet-skinned parents, Palisade Russet and the breeding clone ND028673B-2Russ. Field data collected over a 2-year period and used in the QTL analyses included tuber length-width and width-depth ratios that were obtained using a digital caliper. The width-depth ratio provided an assessment of the "flatness" of a tuber, which is of importance in potato processing. To cross-validate the accuracy and differences among tuber shape measurement methods, a trained evaluator also assessed the identical tubers based on 1-5 scale (compressed to long) visual assessment method. Furthermore, the weights of analyzed tubers and specific gravities were also collected during the phenotyping process for each mapping clone. A major tuber shape QTL was consistently observed on chromosome 10 with both the length-width ratio and visual assessments. On chromosome 4, a significant QTL for tuber shape from the visual assessment phenotypic data was also detected. Additionally, a tuber shape-related QTL on chromosome 6 was also detected from the length-width ratio data from 2020. Chromosome 2 was also identified as having a significant QTL for the width-depth ratio, which is of importance in influencing the flatness of a tuber. One significant QTL for tuber weight (i.e., tuber size) was observed on chromosome 5, and a significant QTL for specific gravity was found on chromosome 3. These significant and major QTL should be useful for developing marker-assisted selection for more efficient potato breeding.

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.

Reproducible Quantitative Trait Loci for Resistance to Soft Rot Caused by Dickeya dianthicola Derived from the Wild Potato Solanum microdontum (PI 458355) Are Located on Chromosomes 1, 3, and 5

The potato wild relative Solanum microdontum is a breeder-friendly source of genetic resistance to soft rot. Our objectives were to (i) identify loci associated with soft rot resistance in S. microdontum germplasm and (ii) develop bi-parental populations in a self-compatible S. tuberosum genetic background to recover segregating F2 progenies, construct a linkage map, and identify quantitative trait loci (QTLs). Under objective (i), tubers from 103 S. microdontum genotypes from the United States Potato Genebank were inoculated with a high virulence strain of Dickeya dianthicola, and lesion size was measured after a 24-h incubation period at 30°C. Association analysis using 3,490 polymorphic Infinium array SNP markers identified soft rot resistance loci on chromosomes 1, 2, 3, 5, 7, 8, 11, and 12. Under objective (ii), a resistant S. microdontum accession PI 458355 was crossed with a highly fertile, self-compatible, diploid S. tuberosum pollen parent (PI 654351) to generate segregating F2 populations. Composite interval mapping was conducted using a genetic linkage map with 970 GBS-based SNP markers. Reproducible QTLs were detected on chromosomes 1, 3, and 5, explaining 11, 13, and 23% of the phenotypic variation, respectively. Homozygous S. microdontum alleles at the QTL on chromosome 5 and heterozygous or homozygous S. microdontum alleles at QTLs on chromosomes 1 and 3 significantly decrease lesion size compared with the homozygous S. tuberosum parent. The germplasm created in these studies provides a resource for studying traits from S. microdontum, and we can use the advanced F2 selections for future potato improvement. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Transcriptomic analysis of Chinese yam (Dioscorea polystachya Turcz.) variants indicates brassinosteroid involvement in tuber development

Dioscorea is an important but underutilized genus of flowering plants that grows predominantly in tropical and subtropical regions. Several species, known as yam, develop large underground tubers and aerial bulbils that are used as food. The Chinese yam (D. polystachya Turcz.) is one of the few Dioscorea species that grows well in temperate regions and has been proposed as a climate-resilient crop to enhance food security in Europe. However, the fragile, club-like tubers are unsuitable for mechanical harvesting, which is facilitated by shorter and thicker storage organs. Brassinosteroids (BRs) play a key role in plant cell division, cell elongation and proliferation, as well as in the gravitropic response. We collected RNA-Seq data from the head, middle and tip of two tuber shape variants: F60 (long, thin) and F2000 (short, thick). Comparative transcriptome analysis of F60 vs. F2000 revealed 30,229 differentially expressed genes (DEGs), 1,393 of which were differentially expressed in the growing tip. Several DEGs are involved in steroid/BR biosynthesis or signaling, or may be regulated by BRs. The quantification of endogenous BRs revealed higher levels of castasterone (CS), 28-norCS, 28-homoCS and brassinolide in F2000 compared to F60 tubers. The highest BR levels were detected in the growing tip, and CS was the most abundant (439.6 ± 196.41 pmol/g in F2000 and 365.6 ± 112.78 pmol/g in F60). Exogenous 24-epi-brassinolide (epi-BL) treatment (20 nM) in an aeroponic system significantly increased the width-to-length ratio (0.045 ± 0.002) compared to the mock-treated plants (0.03 ± 0.002) after 7 weeks, indicating that exogenous epi-BL produces shorter and thicker tubers. In this study we demonstrate the role of BRs in D. polystachya tuber shape, providing insight into the role of plant hormones in yam storage organ development. We found that BRs can influence tuber shape in Chinese yam by regulating the expression of genes involved cell expansion. Our data can help to improve the efficiency of Chinese yam cultivation, which could provide an alternative food source and thus contribute to future food security in Europe.

Converting Hybrid Potato Breeding Science into Practice

Research on diploid hybrid potato has made fast advances in recent years. In this review we give an overview of the most recent and relevant research outcomes. We define different components needed for a complete hybrid program: inbred line development, hybrid evaluation, cropping systems and variety registration. For each of these components the important research results are discussed and the outcomes and issues that merit further study are identified. We connect fundamental and applied research to application in a breeding program, based on the experiences at the breeding company Solynta. In the concluding remarks, we set hybrid breeding in a societal perspective, and we identify bottlenecks that need to be overcome to allow successful adoption of hybrid potato.

GWAS in tetraploid potato: identification and validation of SNP markers associated with glycoalkaloid content

Genome-wide association studies (GWAS) are a useful tool to unravel the genetic architecture of complex traits, but the results can be difficult to interpret. Population structure, genetic heterogeneity, and rare alleles easily result in false positive or false negative associations. This paper describes the analysis of a GWAS panel combined with three bi-parental mapping populations to validate GWAS results, using phenotypic data for steroidal glycoalkaloid (SGA) accumulation and the ratio (SGR) between the two major glycoalkaloids α-solanine and α-chaconine in potato tubers. SGAs are secondary metabolites in the Solanaceae family, functional as a defence against various pests and pathogens and in high quantities toxic for humans. With GWAS, we identified five quantitative trait loci (QTL) of which Sga1.1, Sgr8.1, and Sga11.1 were validated, but not Sga3.1 and Sgr7.1. In the bi-parental populations, Sga5.1 and Sga7.1 were mapped, but these were not identified with GWAS. The QTLs Sga1.1, Sga7.1, Sgr7.1, and Sgr8.1 co-localize with genes GAME9, GAME 6/GAME 11, SGT1, and SGT2, respectively. For other genes involved in SGA synthesis, no QTLs were identified. The results of this study illustrate a number of pitfalls in GWAS of which population structure seems the most important. We also show that introgression breeding for disease resistance has introduced new haplotypes to the gene pool involved in higher SGA levels in certain pedigrees. Finally, we show that high SGA levels remain unpredictable in potato but that α-solanine/α-chaconine ratio has a predictable outcome with specific SGT1 and SGT2 haplotypes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01344-2.

Fruit Morphology and Ripening-Related QTLs in a Newly Developed Introgression Line Collection of the Elite Varieties 'Védrantais' and 'Piel de Sapo'

Melon is an economically important crop with widely diverse fruit morphology and ripening characteristics. Its diploid sequenced genome and multiple genomic tools make this species suitable to study the genetic architecture of fruit traits. With the development of this introgression line population of the elite varieties 'Piel de Sapo' and 'Védrantais', we present a powerful tool to study fruit morphology and ripening traits that can also facilitate characterization or pyramidation of QTLs in inodorous melon types. The population consists of 36 lines covering almost 98% of the melon genome, with an average of three introgressions per chromosome and segregating for multiple fruit traits: morphology, ripening and quality. High variability in fruit morphology was found within the population, with 24 QTLs affecting six different traits, confirming previously reported QTLs and two newly detected QTLs, FLQW5.1 and FWQW7.1. We detected 20 QTLs affecting fruit ripening traits, six of them reported for the first time, two affecting the timing of yellowing of the rind (EYELLQW1.1 and EYELLQW8.1) and four at the end of chromosome 8 affecting aroma, abscission and harvest date (EAROQW8.3, EALFQW8.3, ABSQW8.3 and HARQW8.3). We also confirmed the location of several QTLs, such as fruit-quality-related QTLs affecting rind and flesh appearance and flesh firmness.

QTL analysis of tuber shape in a diploid potato population

Tuber shape is one of the most important traits for potato breeding. Since poor or irregular shape increases the difficulty of handling and processing, researching the inheritance of potato tuber shape for potato breeding is highly important. To efficiently identify QTL for tuber shape, a diploid potato population (PM7) was generated by self-pollinated M6 (S. chacoense). A QTL TScha6 for tuber shape was identified by the QTL-seq approach at 50.91-59.93 Mb on chromosome 6 in the potato DM reference genome. To confirm TScha6, four SSR and twenty CAPS markers around the QTL were developed and the TScha6 was narrowed down to an interval of ~ 1.85 Mb. The CAPS marker C6-58.27_665 linked to TScha6 was then used to screen 86 potato cultivars and advanced breeding lines. The tuber length/width (LW) ratio was significantly correlated with the presence/absence of C6-58.27_665, and the correlation coefficient was r = 0.55 (p < 0.01). These results showed that C6-58.27_665 could be applied in marker-assisted selection (MAS) for tuber shape breeding in the future. Our research sets the important stage for the future cloning of the tuber shape gene and utilities of the marker in the breeding program.

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.

Molecular Approaches to Overcome Self-Incompatibility in Diploid Potatoes

There has been an increased interest in true potato seeds (TPS) as planting material because of their advantages over seed tubers. TPS produced from a tetraploid heterozygous bi-parental population produces non-uniform segregating progenies, which have had limited uniformity in yield and quality in commercial cultivation, and, thus, limited success. Inbreeding depression and self-incompatibility hamper the development of inbred lines in both tetraploid and diploid potatoes, impeding hybrid development efforts. Diploid potatoes have gametophytic self-incompatibility (SI) controlled by S-locus, harboring the male-dependent S-locus F-box (SLF/SFB) and female-dependent Stylar-RNase (S-RNase). Manipulation of these genes using biotechnological tools may lead to loss of self-incompatibility. Self-compatibility can also be achieved by the introgression of S-locus inhibitor (Sli) found in the self-compatible (SC) natural mutants of Solanum chacoense. The introgression of Sli through conventional breeding methods has gained much success. Recently, the Sli gene has been cloned from diverse SC diploid potato lines. It is expressed gametophytically and can overcome the SI in different diploid potato genotypes through conventional breeding or transgenic approaches. Interestingly, it has a 533 bp insertion in its promoter elements, a MITE transposon, making it a SC allele. Sli gene encodes an F-box protein PP2-B10, which consists of an F-box domain linked to a lectin domain. Interaction studies have revealed that the C-terminal region of Sli interacts with most of the StS-RNases, except StS-RNase 3, 9, 10, and 13, while full-length Sli cannot interact with StS-RNase 3, 9, 11, 13, and 14. Thus, Sli may play an essential role in mediating the interactions between pollen and stigma and function like SLFs to interact with and detoxify the S-RNases during pollen tube elongation to confer SC to SI lines. These advancements have opened new avenues in the diploid potato hybrid.

Breeding Diploid F1 Hybrid Potatoes for Propagation from Botanical Seed (TPS): Comparisons with Theory and Other Crops

This paper reviews the progress and the way ahead in diploid F1 hybrid potato breeding by comparisons with expectations from the theory of inbreeding and crossbreeding, and experiences from other diploid outbreeding crops. Diploid potatoes can be converted from an outbreeding species, in which self-pollination is prevented by a gametophytic self-incompatibility system, into one where self-pollination is possible, either through a dominant self-incompatibility inhibitor gene (Sli) or knockout mutations in the incompatibility locus. As a result, diploid F1 hybrid breeding can be used to produce genetically uniform potato cultivars for propagation from true potato seeds by crossing two near-homozygous inbred lines, derived from a number of generations of self-pollination despite inbreeding depression. Molecular markers can be used to detect and remove deleterious recessive mutations of large effect, including those in tight repulsion linkage. Improvements to the inbred lines can be made by introducing and stacking genes and chromosome segments of large desirable effect from wild relatives by backcrossing. Improvements in quantitative traits require a number of cycles of inbreeding and crossbreeding. Seed production can be achieved by hand pollinations. F1 hybrid planting material can be delivered to farmers as true seeds or young plants, and mini-tubers derived from true seeds.

StWRKY13 promotes anthocyanin biosynthesis in potato (Solanum tuberosum) tubers

Although the role of WRKY transcription factors (TFs) in colour formation has been reported in several species, their function in potato (Solanum tuberosum L.) anthocyanin biosynthesis remains unclear. In this study, the potato WRKY gene StWRKY13 was isolated and characterised. Expression analysis revealed a significantly higher StWRKY13 expression in chromatic tubers than in yellow ones. Transient activation assays showed that StWRKY13 could enhance the role of StAN2 in promoting anthocyanin biosynthesis in tobacco (Nicotiana tabacum L.). Over-expressing the StWRKY13 gene promoted anthocyanin biosynthesis in potato tubers. Further investigations indicated that StWRKY13 could interact with the StCHS, StF3H, StDFR, and StANS gene promoters and significantly enhance their activities. Our findings showed that StWRKY13 could promote anthocyanin biosynthesis by activating StCHS, StF3H, StDFR, and StANS transcription in potato tubers, thereby supporting the theoretical basis for anthocyanins formation in coloured 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.

Rare instances of haploid inducer DNA in potato dihaploids and ploidy-dependent genome instability

In cultivated tetraploid potato (Solanum tuberosum), reduction to diploidy (dihaploidy) allows for hybridization to diploids and introgression breeding and may facilitate the production of inbreds. Pollination with haploid inducers (HIs) yields maternal dihaploids, as well as triploid and tetraploid hybrids. Dihaploids may result from parthenogenesis, entailing the development of embryos from unfertilized eggs, or genome elimination, entailing missegregation and the loss of paternal chromosomes. A sign of genome elimination is the occasional persistence of HI DNA in some dihaploids. We characterized the genomes of 919 putative dihaploids and 134 hybrids produced by pollinating tetraploid clones with three HIs: IVP35, IVP101, and PL-4. Whole-chromosome or segmental aneuploidy was observed in 76 dihaploids, with karyotypes ranging from 2n = 2x - 1 = 23 to 2n = 2x + 3 = 27. Of the additional chromosomes in 74 aneuploids, 66 were from the non-inducer parent and 8 from the inducer parent. Overall, we detected full or partial chromosomes from the HI parent in 0.87% of the dihaploids, irrespective of parental genotypes. Chromosomal breaks commonly affected the paternal genome in the dihaploid and tetraploid progeny, but not in the triploid progeny, correlating instability to sperm ploidy and to haploid induction. The residual HI DNA discovered in the progeny is consistent with genome elimination as the mechanism of haploid induction.

First-generation genome editing in potato using hairy root transformation

Genome editing and cis-gene breeding have rapidly accelerated crop improvement efforts, but their impacts are limited by the number of species capable of being genetically transformed. Many dicot species, including some vital potato relatives being used to accelerate breeding and genetics efforts, remain recalcitrant to standard Agrobacterium tumefaciens-based transformation. Hairy root transformation using Agrobacterium rhizogenes (A. rhizogenes) provides an accelerated approach to generating transgenic material but has been limited to analysis of hairy root clones. In this study, strains of A. rhizogenes were tested in the wild diploid potato relative Solanum chacoense, which is recalcitrant to infection by Agrobacterium tumefaciens. One strain of A. rhizogenes MSU440 emerged as being capable of delivering a T-DNA carrying the GUS marker and generating transgenic hairy root clones capable of GUS expression and regeneration to whole plants. CRISPR/Cas9 reagents targeting the potato PHYTOENE DESATURASE (StPDS) gene were expressed in hairy root clones and regenerated. We found that 64%-98% of transgenic hairy root clones expressing CRISPR/Cas9 reagents carried targeted mutations, while only 14%-30% of mutations were chimeric. The mutations were maintained in regenerated lines as stable mutations at rates averaging at 38% and were capable of germ-line transmission to progeny. This novel approach broadens the numbers of genotypes amenable to Agrobacterium-mediated transformation while reducing chimerism in primary events and accelerating the generation of edited materials.

Solyntus, the New Highly Contiguous Reference Genome for Potato (Solanum tuberosum)

With the rapid expansion of the application of genomics and sequencing in plant breeding, there is a constant drive for better reference genomes. In potato (Solanum tuberosum), the third largest food crop in the world, the related species S. phureja, designated “DM”, has been used as the most popular reference genome for the last 10 years. Here, we introduce the de novo sequenced genome of Solyntus as the next standard reference in potato genome studies. A true Solanum tuberosum made up of 116 contigs that is also highly homozygous, diploid, vigorous and self-compatible, Solyntus provides a more direct and contiguous reference then ever before available. It was constructed by sequencing with state-of-the-art long and short read technology and assembled with Canu. The 116 contigs were assembled into scaffolds to form each pseudochromosome, with three contigs to 17 contigs per chromosome. This assembly contains 93.7% of the single-copy gene orthologs from the Solanaceae set and has an N50 of 63.7 Mbp. The genome and related files can be found at https://www.plantbreeding.wur.nl/Solyntus/. With the release of this research line and its draft genome we anticipate many exciting developments in (diploid) potato research.

The origin and widespread occurrence of Sli-based self-compatibility in potato

Self-compatible (SC) diploid potatoes allow innovative potato breeding. Therefore, the Sli gene, originally described in S. chacoense, has received much attention. In elite S. tuberosum diploids, spontaneous berry set is occasionally observed. We aimed to map SC from S. tuberosum origin. Two full-sib mapping populations from non-inbred diploids were used. Bulks were composed based on both pollen tube growth and berry set upon selfing. After DNA sequencing of the parents and bulks, we generated k-mer tables. Set algebra and depth filtering were used to identify bulk-specific k-mers. Coupling and repulsion phase k-mers, transmitted from the SC parent, mapped in both populations to the distal end of chromosome 12. Intersection between the k-mers from both populations, in coupling phase with SC, exposed a shared haplotype of approximately 1.5 Mb. Subsequently, we screened read archives of potatoes and wild relatives for k-mers specific to this haplotype. The well-known SC clones US-W4 and RH89-039-16, but surprisingly, also S. chacoense clone M6 were positives. Hence, the S. tuberosum source of SC seems identical to Sli. Furthermore, the candidate region drastically reduced to 333 kb. Haplotype-specific KASP markers were designed and validated on a panel of diploid clones including another renown SC dihaploid G254. Interestingly, k-mers specific to the SC haplotype were common in tetraploid varieties. Pedigree information suggests that the SC haplotype was introduced into tetraploid varieties via the founder "Rough Purple Chili". We show that Sli is surprisingly widespread and indigenous to the cultivated gene pool of potato.

An SNP-Based High-Density Genetic Linkage Map for Tetraploid Potato Using Specific Length Amplified Fragment Sequencing (SLAF-Seq) Technology

Specific length amplified fragment sequencing (SLAF-seq) is a recently developed high-resolution strategy for the discovery of large-scale de novo genotyping of single nucleotide polymorphism (SNP) markers. In the present research, in order to facilitate genome-guided breeding in potato, this strategy was used to develop a large number of SNP markers and construct a high-density genetic linkage map for tetraploid potato. The genomic DNA extracted from 106 F1 individuals derived from a cross between two tetraploid potato varieties YSP-4 × MIN-021 and their parents was used for high-throughput sequencing and SLAF library construction. A total of 556.71 Gb data, which contained 2269.98 million pair-end reads, were obtained after preprocessing. According to bioinformatics analysis, a total of 838,604 SLAF labels were developed, with an average sequencing depth of 26.14-fold for parents and 15.36-fold for offspring of each SLAF, respectively. In total, 113,473 polymorphic SLAFs were obtained, from which 7638 SLAFs were successfully classified into four segregation patterns. After filtering, a total of 7329 SNP markers were detected for genetic map construction. The final integrated linkage map of tetraploid potato included 3001 SNP markers on 12 linkage groups, and covered 1415.88 cM, with an average distance of 0.47 cM between adjacent markers. To our knowledge, the integrated map described herein has the best coverage of the potato genome and the highest marker density for tetraploid potato. This work provides a foundation for further quantitative trait loci (QTL) location, map-based gene cloning of important traits and marker-assisted selection (MAS) of potato.

The genetic basis of inbreeding depression in potato

Inbreeding depression confers reduced fitness among the offspring of genetic relatives. As a clonally propagated crop, potato (Solanum tuberosum L.) suffers from severe inbreeding depression; however, the genetic basis of inbreeding depression in potato is largely unknown. To gain insight into inbreeding depression in potato, we evaluated the mutation burden in 151 diploid potatoes and obtained 344,831 predicted deleterious substitutions. The deleterious mutations in potato are enriched in the pericentromeric regions and are line specific. Using three F₂ populations, we identified 15 genomic regions with severe segregation distortions due to selection at the gametic and zygotic stages. Most of the deleterious recessive alleles affecting survival and growth vigor were located in regions with high recombination rates. One of these deleterious alleles is derived from a rare mutation that disrupts a gene required for embryo development. This study provides the basis for genome design of potato inbred lines.

A common genetic mechanism underlies morphological diversity in fruits and other plant organs

Shapes of edible plant organs vary dramatically among and within crop plants. To explain and ultimately employ this variation towards crop improvement, we determined the genetic, molecular and cellular bases of fruit shape diversity in tomato. Through positional cloning, protein interaction studies, and genome editing, we report that OVATE Family Proteins and TONNEAU1 Recruiting Motif proteins regulate cell division patterns in ovary development to alter final fruit shape. The physical interactions between the members of these two families are necessary for dynamic relocalization of the protein complexes to different cellular compartments when expressed in tobacco leaf cells. Together with data from other domesticated crops and model plant species, the protein interaction studies provide possible mechanistic insights into the regulation of morphological variation in plants and a framework that may apply to organ growth in all plant species.

Remarkable organ shape morphological diversity exists in fruits, vegetables and seeds. Here, the authors establish a link between OVATE Family Proteins and TONNEAU1 Recruiting Motif family proteins in the development pathway that governs fruit shape of tomato, melon, and cucumber as well as potato tuber shape.

Generation of self-compatible diploid potato by knockout of S-RNase

Re-domestication of potato into an inbred line-based diploid crop propagated by seed represents a promising alternative to traditional clonal propagation of tetraploid potato, but self-incompatibility has hindered the development of inbred lines. To address this problem, we created self-compatible diploid potatoes by knocking out the self-incompatibility gene S-RNase using the CRISPR-Cas9 system. This strategy opens new avenues for diploid potato breeding and will also be useful for studying other self-incompatible crops.

The Changing Face of Bacterial Soft-Rot Diseases

Bacterial soft rot is a disease complex caused by multiple genera of gram-negative and gram-positive bacteria, with Dickeya and Pectobacterium being the most widely studied soft-rot bacterial pathogens. In addition to soft rot, these bacteria also cause blackleg of potato, foot rot of rice, and bleeding canker of pear. Multiple Dickeya and Pectobacterium species cause the same symptoms on potato, complicating epidemiology and disease resistance studies. The primary pathogen species present in potato-growing regions differs over time and space, further complicating disease management. Genomics technologies are providing new management possibilities, including improved detection and biocontrol methods that may finally allow effective disease management. The recent development of inbred diploid potato lines is also having a major impact on studying soft-rot pathogens because it is now possible to study soft-rot disease in model plant species that produce starchy vegetative storage organs. Together, these new discoveries have changed how we face diseases caused by these pathogens.

Plant Resistance to Colorado Potato Beetle (Coleoptera: Chrysomelidae) in Diploid F2 Families Derived From Crosses Between Cultivated and Wild Potato

Colorado potato beetle, Leptinotarsa decemlineata Say (Coleoptera: Chrysomelidae), is a serious global pest of potato, Solanum tuberosum L. Management of L. decemlineata has relied heavily on insecticides, but repeated evolution of insecticide resistance has motivated the exploration and development of alternative strategies, such as plant resistance. The recent development of two diploid potato families derived from crosses between cultivated and wild potato species (S. chacoense and S. berthaultii) has provided a unique opportunity to reexamine plant traits for resistance breeding. In this 2-yr study, we surveyed select F2 clones for the induction of L. decemlineata mortality and a reduction in defoliation in no-choice feeding assays when challenged with adults and larvae from three sites in Wisconsin. We tested for an association with glandular trichome density and foliar levels of the glycoalkaloids chaconine and solanine. Several potato clones demonstrated resistance in specific feeding assays, but none excelled consistently across experiments. Mortality and defoliation generally differed significantly among L. decemlineata populations, which could be indicative of heritable variation in beetle responses to plant defenses or variation in the physiological status of the beetle populations tested. Contrary to expectations, higher trichome density increased mortality or decreased defoliation in only a few cases, and levels of mortality and defoliation were unrelated to foliar glycoalkaloid content, warranting further investigation of the defense mechanisms of resistant clones. In addition to identifying several potential L. decemlineata resistance sources, this study underscores the need to include multiple insect populations in surveys of plant resistance to this diverse pest.

QTL mapping in diploid potato by using selfed progenies of the cross S. tuberosum × S. chacoense

Usually, mapping studies in potato are performed with segregating populations from crosses between highly heterozygous diploid or tetraploid parents. These studies are hampered by a high level of genetic background noise due to the numerous segregating alleles, with a maximum of eight per locus. In the present study, we aimed to increase the mapping efficiency by using progenies from diploid inbred populations in which at most two alleles segregate. Selfed progenies were generated from a cross between S. tuberosum (D2; a highly heterozygous diploid) and S. chacoense (DS; a homozygous diploid clone) containing the self-incompatibility overcoming S locus inhibitor (Sli-gene). The Sli-gene enables self-pollination and the generation of selfed progenies. One F2 population was used to map several quality traits, such as tuber shape, flesh and skin color. Quantitative trait loci were identified for almost all traits under investigation. The identified loci partially coincided with known mapped loci and partially identified new loci. Nine F3 populations were used to validate the QTLs and monitor the overall increase in the homozygosity level.

Meta-analysis of chromosome-scale crossover rate variation in eukaryotes and its significance to evolutionary genomics

Understanding the distribution of crossovers along chromosomes is crucial to evolutionary genomics because the crossover rate determines how strongly a genome region is influenced by natural selection on linked sites. Nevertheless, generalities in the chromosome-scale distribution of crossovers have not been investigated formally. We fill this gap by synthesizing joint information on genetic and physical maps across 62 animal, plant and fungal species. Our quantitative analysis reveals a strong and taxonomically widespread reduction of the crossover rate in the centre of chromosomes relative to their peripheries. We demonstrate that this pattern is poorly explained by the position of the centromere, but find that the magnitude of the relative reduction in the crossover rate in chromosome centres increases with chromosome length. That is, long chromosomes often display a dramatically low crossover rate in their centre, whereas short chromosomes exhibit a relatively homogeneous crossover rate. This observation is compatible with a model in which crossover is initiated from the chromosome tips, an idea with preliminary support from mechanistic investigations of meiotic recombination. Consequently, we show that organisms achieve a higher genome-wide crossover rate by evolving smaller chromosomes. Summarizing theory and providing empirical examples, we finally highlight that taxonomically widespread and systematic heterogeneity in crossover rate along chromosomes generates predictable broad-scale trends in genetic diversity and population differentiation by modifying the impact of natural selection among regions within a genome. We conclude by emphasizing that chromosome-scale heterogeneity in crossover rate should urgently be incorporated into analytical tools in evolutionary genomics, and in the interpretation of resulting patterns.

Genome sequence of M6, a diploid inbred clone of the high-glycoalkaloid-producing tuber-bearing potato species Solanum chacoense, reveals residual heterozygosity

Cultivated potato (Solanum tuberosum L.) is a highly heterozygous autotetraploid that presents challenges in genome analyses and breeding. Wild potato species serve as a resource for the introgression of important agronomic traits into cultivated potato. One key species is Solanum chacoense and the diploid, inbred clone M6, which is self-compatible and has desirable tuber market quality and disease resistance traits. Sequencing and assembly of the genome of the M6 clone of S. chacoense generated an assembly of 825 767 562 bp in 8260 scaffolds with an N50 scaffold size of 713 602 bp. Pseudomolecule construction anchored 508 Mb of the genome assembly into 12 chromosomes. Genome annotation yielded 49 124 high-confidence gene models representing 37 740 genes. Comparative analyses of the M6 genome with six other Solanaceae species revealed a core set of 158 367 Solanaceae genes and 1897 genes unique to three potato species. Analysis of single nucleotide polymorphisms across the M6 genome revealed enhanced residual heterozygosity on chromosomes 4, 8 and 9 relative to the other chromosomes. Access to the M6 genome provides a resource for identification of key genes for important agronomic traits and aids in genome-enabled development of inbred diploid potatoes with the potential to accelerate potato breeding.

Exploring the potential and limitations of genotyping-by-sequencing for SNP discovery and genotyping in tetraploid potato

Genotyping-by-sequencing (GBS) potentially offers a cost-effective alternative for SNP discovery and genotyping. Here, we report the exploration of GBS in tetraploid potato. Both ApeKI and PstI/MspI enzymes were used for library preparation on eight diverse potato genotypes. ApeKI yielded more markers than PstI/MspI but provided a lower read coverage per marker, resulting in more missing data and limiting effective genotyping to the tetraploid mode. We then assessed the accuracy of these SNPs by comparison with SolCAP data (5824 data points in diploid mode and 3243 data points in tetraploid mode) and found the match rates between genotype calls was 90.4% and 81.3%, respectively. Imputation of missing data did not prove very accurate because of incomplete haplotype discovery, suggesting caution in setting the allowance for missing data. To further assess the quality of GBS-derived data, a genome-wide association analysis was performed for flower color on 318 clones (with ApeKI). A strong association signal on chromosome 2 was obtained with the most significant SNP located in the middle of the dihydroflavonol 4-reductase (DFR) gene. We conclude that an appropriate choice of enzyme for GBS library preparation makes it possible to obtain high-quality SNPs in potato and will be helpful for marker-assisted genomics.

Automated tetraploid genotype calling by hierarchical clustering

New software to make tetraploid genotype calls from SNP array data was developed, which uses hierarchical clustering and multiple F1 populations to calibrate the relationship between signal intensity and allele dosage. SNP arrays are transforming breeding and genetics research for autotetraploids. To fully utilize these arrays, the relationship between signal intensity and allele dosage must be calibrated for each marker. We developed an improved computational method to automate this process, which is provided as the R package ClusterCall. In the training phase of the algorithm, hierarchical clustering within an F1 population is used to group samples with similar intensity values, and allele dosages are assigned to clusters based on expected segregation ratios. In the prediction phase, multiple F1 populations and the prediction set are clustered together, and the genotype for each cluster is the mode of the training set samples. A concordance metric, defined as the proportion of training set samples equal to the mode, can be used to eliminate unreliable markers and compare different algorithms. Across three potato families genotyped with an 8K SNP array, ClusterCall scored 5729 markers with at least 0.95 concordance (94.6% of its total), compared to 5325 with the software fitTetra (82.5% of its total). The three families were used to predict genotypes for 5218 SNPs in the SolCAP diversity panel, compared with 3521 SNPs in a previous study in which genotypes were called manually. One of the additional markers produced a significant association for vine maturity near a well-known causal locus on chromosome 5. In conclusion, when multiple F1 populations are available, ClusterCall is an efficient method for accurate, autotetraploid genotype calling that enables the use of SNP data for research and plant breeding.

Genetic diversity and association mapping in the Colombian Central Collection of Solanum tuberosum L. Andigenum group using SNPs markers

The potato (Solanum tuberosum L.) is the fourth most important crop food in the world and Colombia has one of the most important collections of potato germplasm in the world (the Colombian Central Collection-CCC). Little is known about its potential as a source of genetic diversity for molecular breeding programs. In this study, we analyzed 809 Andigenum group accessions from the CCC using 5968 SNPs to determine: 1) the genetic diversity and population structure of the Andigenum germplasm and 2) the usefulness of this collection to map qualitative traits across the potato genome. The genetic structure analysis based on principal components, cluster analyses, and Bayesian inference revealed that the CCC can be subdivided into two main groups associated with their ploidy level: Phureja (diploid) and Andigena (tetraploid). The Andigena population was more genetically diverse but less genetically substructured than the Phureja population (three vs. five subpopulations, respectively). The association mapping analysis of qualitative morphological data using 4666 SNPs showed 23 markers significantly associated with nine morphological traits. The present study showed that the CCC is a highly diverse germplasm collection genetically and phenotypically, useful to implement association mapping in order to identify genes related to traits of interest and to assist future potato genetic breeding programs.