Dissection of the Plant Hormone Signal Transduction Network in Late Blight-Resistant Potato Genotype SD20 and Prediction of Key Resistance Genes.

Hormones play an important role in plant disease resistance and defense. Transcriptome data of late blight-resistant potato genotype SD20 treated by ethylene (ET), jasmonate (JA), salicylic acid (SA), and Phytophthora infestans CN152 was analyzed to assess the role of the ET/JA/SA regulatory network in plant disease resistance and defense and predict key resistant genes. The results show that there was significant crossover of differentially expressed genes among all treatments, and common and specific plant disease interaction genes for the ET, JA, and SA treatments were differentially expressed in the CN152 treatment. The resistance and defense genes of the potato genotype SD20 could be induced to regulate metabolic and hormone signaling pathways by alternative splicing in all treatments. Further analysis found that JA and ET pathways can work together synergistically. JA/ET and SA pathways antagonize each other to initiate the expression of calmodulin-domain protein kinases and calmodulin/calmodulin-like and RPM1-interacting protein 4 genes, and they activate HSP-mediated hypersensitive response and defense-related genes. Meanwhile, nine defense genes, including wound-responsive AP2-like factor, glutathione-s-transferase, serine/threonine-protein kinase BRI1, and Avr9/Cf-9 rapidly elicited protein genes, were obtained using weighted gene coexpression network analysis, which provided reliable targets for functional verification. This study provides a theoretical reference for the comprehensive application of plant hormones to improve resistance to potato late blight disease.

Genome-wide identification and gene expression analysis of the 14-3-3 gene family in potato (Solanum tuberosum L.).

BACKGROUND: 14-3-3 proteins are essential in regulating various biological processes and abiotic stress responses in plants. Although 14-3-3 proteins have been studied in model plants such as Arabidopsis thaliana and Oryza sativa, there is a lack of research on the 14-3-3 gene family in potatoes (Solanum tuberosum L.). RESULTS: A total of 18 14-3-3 genes encoding proteins containing a typical conserved PF00244 domain were identified by genome-wide analysis in potatoes. The St14-3-3 gene family members were unevenly distributed across the chromosomes, and gene structure analysis showed that gene length and intron number varied greatly among the members. Phylogenetic analysis of 14-3-3 proteins in potatoes and other plant species showed that they could be divided into two distinct groups (ε and non-ε). Members in the ε group tended to have similar exon-intron structures and conserved motif patterns. Promoter sequence analysis showed that the St14-3-3 gene promoters contained multiple hormone-, stress-, and light-responsive cis-regulatory elements. Synteny analysis suggested that segmental duplication events contributed to the expansion of the St14-3-3 gene family in potatoes. The observed syntenic relationships between some 14-3-3 genes from potato, Arabidopsis, and tomato suggest that they evolved from a common ancestor. RNA-seq data showed that St14-3-3 genes were expressed in all tissues of potatoes but that their expression patterns were different. qRT-PCR assays revealed that the expression levels of nearly all tested St14-3-3 genes were affected by drought, salt, and low-temperature stresses and that different St14-3-3 genes had different responses to these stresses. CONCLUSIONS: In summary, genome-wide identification, evolutionary, and expression analyses of the 14-3-3 gene family in potato were conducted. These results provide important information for further studies on the function and regulation of St14-3-3 gene family members in potatoes.

Genome architecture and tetrasomic inheritance of autotetraploid potato.

Potato (Solanum tuberosum) is the most consumed non-cereal food crop. Most commercial potato cultivars are autotetraploids with highly heterozygous genomes, severely hampering genetic analyses and improvement. By leveraging the state-of-the-art sequencing technologies and polyploid graph binning, we achieved a chromosome-scale, haplotype-resolved genome assembly of a cultivated potato, Cooperation-88 (C88). Intra-haplotype comparative analyses revealed extensive sequence and expression differences in this tetraploid genome. We identified haplotype-specific pericentromeres on chromosomes, suggesting a distinct evolutionary trajectory of potato homologous centromeres. Furthermore, we detected double reduction events that are unevenly distributed on haplotypes in 1021 of 1034 selfing progeny, a feature of autopolyploid inheritance. By distinguishing maternal and paternal haplotype sets in C88, we simulated the origin of heterosis in cultivated tetraploid with a survey of 3110 tetra-allelic loci with deleterious mutations, which were masked in the heterozygous condition by two parents. This study provides insights into the genomic architecture of autopolyploids and will guide their breeding.

Fine Mapping and Candidate Gene Prediction of Tuber Shape Controlling Ro Locus Based on Integrating Genetic and Transcriptomic Analyses in Potato.

Tuber shape is one of the most important quality traits in potato appearance. Since poor or irregular shape results in higher costs for processing and influences the consumers' willingness to purchase, breeding for shape uniformity and shallow eye depth is highly important. Previous studies showed that the major round tuber shape controlling locus, the Ro locus, is located on chromosome 10. However, fine mapping and cloning of tuber shape genes have not been reported. In this study, the analyses of tissue sectioning and transcriptome sequencing showed that the developmental differences between round and elongated tuber shapes begin as early as the hook stage of the stolon. To fine map tuber shape genes, a high-density genetic linkage map of the Ro region on chromosome 10 based on a diploid segregating population was constructed. The total length of the genetic linkage map was 25.8 cM and the average marker interval was 1.98 cM. Combined with phenotypic data collected from 2014 to 2017, one major quantitative trait locus (QTL) for tuber shape was identified, which explained 61.7-72.9% of the tuber shape variation. Through the results of genotyping and phenotypic investigation of recombinant individuals, Ro was fine mapped in a 193.43 kb interval, which contained 18 genes. Five candidate genes were preliminarily predicted based on tissue sections and transcriptome sequencing. This study provides an important basis for cloning Ro gene(s).

Salicylic Acid Improves the Constitutive Freezing Tolerance of Potato as Revealed by Transcriptomics and Metabolomics Analyses.

Freezing severely impacts potato production. Deciphering the pathways and metabolites that regulate the freezing tolerance of potato is useful in cultivation and breeding for hardiness. In the present study, Solanum acaule was identified to be more freezing tolerant than S. tuberosum. Furthermore, the two genotypes before/after exposure to 4 °C for 7 d with additional -1 °C for 12 h were analysed by RNA-seq and metabolomics, and the results were compared with the previous -1 °C for 12 h. The results showed that S. acaule activated numerous genes that differed from those of S. tuberosum. Among the genes, five pathways, such as the hormone signalling pathway, which includes salicylic acid, were enriched. Further metabolomics analysis showed that the content of salicylic acid was improved in S. acaule in response to -1 °C for 12 h. Moreover, exogenous application of 0.1 mM salicylic acid to potato was shown to improve constitutive freezing tolerance and increase the expression of HSFC1. Following transcriptome and metabolome analyses, it was documented that the content of SA that increased in freezing-tolerant S. acaule after exposure to cold condition, associated with the SA signalling pathway, enhanced potato freezing tolerance, probably through HSFC1.

Genome-wide simple sequence repeat markers in potato: abundance, distribution, composition, and polymorphism.

Simple sequence repeats (SSRs) are important sources of genetic diversity and are widely used as markers in genetics and molecular breeding. In this study, we examined four potato genomes of DM1-3 516 R44 (DM) from Solanum phureja, RH89039-16 (RH) from Solanum tuberosum, M6 from Solanum chacoense and Solanum commersonii to determine SSR abundance and distribution and develop a larger list of polymorphic markers for a potentially wide range of uses for the potato community. A total of 1,734,619 SSRs were identified across the four genomes with an average of 433,655 SSRs per genome and 2.31kb per SSR. The most abundant repeat units for mono-, di-, tri-, and tetra-nucleotide SSRs were (A/T)n, (AT/AT)n, (AAT/ATT)n, and (ATAT/ATAT)n, respectively. The SSRs were most abundant (78.79%) in intergenic regions and least abundant (3.68%) in untranslated regions. On average, 168,069 SSRs with unique flanking sequences were identified in the four genomes. Further, we identified 16,245 polymorphic SSR markers among the four genomes. Experimental validation confirmed 99.69% of tested markers could generate target bands. The high-density potato SSR markers developed in this study will undoubtedly facilitate the application of SSR markers for genetic research and marker-pyramiding in potato breeding.

Transcriptome analysis reveals that exogenous ethylene activates immune and defense responses in a high late blight resistant potato genotype.

Ethylene (ET) is one of the many important signaling hormones that functions in regulating defense responses in plants. Gene expression profiling was conducted under exogenous ET application in the high late blight resistant potato genotype SD20 and the specific transcriptional responses to exogenous ET in SD20 were revealed. Analysis of differentially expressed genes (DEGs) generated a total of 1226 ET-specific DEGs, among which transcription factors, kinases, defense enzymes and disease resistance-related genes were significantly differentially expressed. GO enrichment and KEGG metabolic pathway analysis also revealed that numerous defense regulation-related genes and defense pathways were significantly enriched. These results were consistent with the interaction of SD20 and Phytophthora infestans in our previous study, indicating that exogenous ET stimulated the defense response and initiated a similar defense pathway compared to pathogen infection in SD20. Moreover, multiple signaling pathways including ET, salicylic acid, jasmonic acid, abscisic acid, auxin, cytokinin and gibberellin were involved in the response to exogenous ET, which indicates that many plant hormones work together to form a complex network to resist external stimuli in SD20. ET-induced gene expression profiling provides insights into the ET signaling transduction pathway and its potential mechanisms in disease defense systems in potato.

Aldehyde dehydrogenase plays crucial roles in response to lower temperature stress in Solanum tuberosum and Nicotiana benthamiana.

The aim of this study is to elucidate the role of ALDH2B7a during the response to lower temperature in Solanum tuberosum. This gene was found to have altered intragenic DNA methylation status in our previous reports. A total of 18 orthologs of StALDH2B7a were identified in the S. tuberosum genome, which were then divided into 8 aldehyde dehydrogenase (ALDH) subfamilies. The methylation statuses of four intragenic cytosine sites in intron 5 and exon 6 of genomic StALDH2B7a were altered by lower temperature stress, resulting in changes in the expression of StALDH2B7a. Silencing of NbALDH2C4, a homolog of StALDH2B7a in Nicotiana benthamiana, resulted in plants which were sensitive to lower temperature and accumulation of reactive oxygen species (ROS) and malondialdehyde (MDA). These data suggested that the expression of StALDH2B7a was upregulated by alteration of its intragenic cytosine methylation status during lower temperature stress, and additional StALDH2B7a enzymes scavenged excess aldehydes resulting from ROS in a response to cold stress in potato. Our study expands the understanding of the mechanisms involved in plant responses to lower temperature, and provides a new gene source to improve potato tolerance to cold stress in northern China, where lower temperature is one of the key limiting factors for crop production.

New Findings on the Resistance Mechanism of an Elite Diploid Wild Potato Species JAM1-4 in Response to a Super Race Strain of Phytophthora infestans.

Late blight is a devastating potato disease worldwide, caused by Phytophthora infestans. The P. infestans strain 2013-18-306 from Yunnan is a "supervirulent race" that overcomes all 11 known late blight resistance genes (R1 to R11) from Solanum demissum. In a previous study, we identified a diploid wild-type potato JAM1-4 (S. jamesii) with high resistance to 2013-18-306. dRenSeq analysis indicated the presence of novel R genes in JAM1-4. RNA-Seq was used to analyze the late blight resistance response genes and defense regulatory mechanisms of JAM1-4 against 2013-18-306. Gene ontology enrichment and KEGG pathway analysis showed that many disease-resistant pathways were significantly enriched. Analysis of differentially expressed genes (DEGs) revealed an active disease resistance mechanism of JAM1-4, and the essential role of multiple signal transduction pathways and secondary metabolic pathways comprised of SA-JA-ET in plant immunity. We also found that photosynthesis in JAM1-4 was inhibited to promote the immune response. Our study reveals the pattern of resistance-related gene expression in response to a super race strain of potato late blight and provides a theoretical basis for further exploration of potato disease resistance mechanisms, discovery of new late blight resistance genes, and disease resistance breeding.

Transcriptome analysis uncovers the gene expression profile of salt-stressed potato (Solanum tuberosum L.).

Potato (Solanum tuberosum L.) is an important staple food worldwide. However, its growth has been heavily suppressed by salt stress. The molecular mechanisms of salt tolerance in potato remain unclear. It has been shown that the tetraploid potato Longshu No. 5 is a salt-tolerant genotype. Therefore, in this study we conducted research to identify salt stress response genes in Longshu No. 5 using a NaCl treatment and time-course RNA sequencing. The total number of differentially expressed genes (DEGs) in response to salt stress was 5508. Based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, it was found that DEGs were significantly enriched in the categories of nucleic acid binding, transporter activity, ion or molecule transport, ion binding, kinase activity and oxidative phosphorylation. Particularly, the significant differential expression of encoding ion transport signaling genes suggests that this signaling pathway plays a vital role in salt stress response in potato. Finally, the DEGs in the salt response pathway were verified by Quantitative real-time PCR (qRT-PCR). These results provide valuable information on the salt tolerance of molecular mechanisms in potatoes, and establish a basis for breeding salt-tolerant cultivars.

Cysteine Protease Inhibitors Reduce Enzymatic Browning of Potato by Lowering the Accumulation of Free Amino Acids.

Enzymatic browning is a major issue affecting the quality of processed potato (Solanum tuberosum L.). To understand the molecular mechanism of browning, transcriptional analyses were performed by employing potatoes that differed in browning. Coexpression analysis indicated that 9 out of 15 upregulated genes in browning-less groups encoded for potato protease inhibitors (StPIs). In addition, gene otology analysis showed that the enriched terms were mainly involved in protease inhibitors. Overexpression of cysteine StPI 143 and StPI 146 individually reduced browning and lowered protease activities and tyrosine and total free amino acid (FAA) contents, but they could not decrease polyphenol oxidase activity. Moreover, supplementing exogenous tyrosine or total FAAs into transgenic potato mash to wild-type amounts promoted mash browning, browning with total FAAs, more than with tyrosine, resembling wild-type levels. These results implied that cysteine StPIs reduced browning via lowering the accumulation of FAAs in addition to tyrosine. Our findings have enriched the knowledge about the roles and mechanisms of protease inhibitors in regulating enzymatic browning of potato, which provide new ways for controlling potato browning.

Comparative Transcriptome Profiling Reveals Compatible and Incompatible Patterns of Potato Toward Phytophthora infestans.

Late blight, caused by Phytophthora infestans (P. infestans), is a devastating disease in potato worldwide. Our previous study revealed that the Solanum andigena genotype 03112-233 is resistant to P. infestans isolate 90128, but susceptible to the super race isolate, CN152. In this study, we confirmed by diagnostic resistance gene enrichment sequencing (dRenSeq) that the resistance of 03112-233 toward 90128 is most likely based on a distinct new R gene(s). To gain an insight into the mechanism that governs resistance or susceptibility in 03112-223, comparative transcriptomic profiling analysis based on RNAseq was initiated. Changes in transcription at two time points (24 h and 72 h) after inoculation with isolates 90128 or CN152 were analyzed. A total of 8,881 and 7,209 genes were differentially expressed in response to 90128 and CN152, respectively, and 1,083 differentially expressed genes (DEGs) were common to both time points and isolates. A substantial number of genes were differentially expressed in an isolate-specific manner with 3,837 genes showing induction or suppression following infection with 90128 and 2,165 genes induced or suppressed after colonization by CN152. Hierarchical clustering analysis suggested that isolates with different virulence profiles can induce different defense responses at different time points. Further analysis revealed that the compatible interaction caused higher induction of susceptibility genes such as SWEET compared with the incompatible interaction. The salicylic acid, jasmonic acid, and abscisic acid mediated signaling pathways were involved in the response against both isolates, while ethylene and brassinosteroids mediated defense pathways were suppressed. Our results provide a valuable resource for understanding the interactions between P. infestans and potato.

A high-throughput BAC end analysis protocol (BAC-anchor) for profiling genome assembly and physical mapping.

Traditional approaches for sequencing insertion ends of bacterial artificial chromosome (BAC) libraries are laborious and expensive, which are currently some of the bottlenecks limiting a better understanding of the genomic features of auto- or allopolyploid species. Here, we developed a highly efficient and low-cost BAC end analysis protocol, named BAC-anchor, to identify paired-end reads containing large internal gaps. Our approach mainly focused on the identification of high-throughput sequencing reads carrying restriction enzyme cutting sites and searching for large internal gaps based on the mapping locations of both ends of the reads. We sequenced and analysed eight libraries containing over 3 200 000 BAC end clones derived from the BAC library of the tetraploid potato cultivar C88 digested with two restriction enzymes, Cla I and Mlu I. About 25% of the BAC end reads carrying cutting sites generated a 60-100 kb internal gap in the potato DM reference genome, which was consistent with the mapping results of Sanger sequencing of the BAC end clones and indicated large differences between autotetraploid and haploid genotypes in potato. A total of 5341 Cla I- and 165 Mlu I-derived unique reads were distributed on different chromosomes of the DM reference genome and could be used to establish a physical map of target regions and assemble the C88 genome. The reads that matched different chromosomes are especially significant for the further assembly of complex polyploid genomes. Our study provides an example of analysing high-coverage BAC end libraries with low sequencing cost and is a resource for further genome sequencing studies.

Transcriptome Profiling Reveals Effects of Drought Stress on Gene Expression in Diploid Potato Genotype P3-198.

Potato (Solanum tuberosum L.) is one of the three most important food crops worldwide; however, it is strongly affected by drought stress. The precise molecular mechanisms of drought stress response in potato are not very well understood. The diploid potato genotype P3-198 has been verified to be highly resistant to drought stress. Here, a time-course experiment was performed to identify drought resistance response genes in P3-198 under polyethylene glycol (PEG)-induced stress using RNA-sequencing. A total of 1665 differentially expressed genes (DEGs) were specifically identified, and based on gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the transcription factor activity, protein kinase activity, and the plant hormone signal transduction process were significantly enriched. Annotation revealed that these DEGs mainly encode transcription factors, protein kinases, and proteins related to redox regulation, carbohydrate metabolism, and osmotic adjustment. In particular, genes encoding abscisic acid (ABA)-dependent signaling molecules were significantly differentially expressed, which revealed the important roles of the ABA-dependent signaling pathway in the early response of P3-198 to drought stress. Quantitative real-time PCR experimental verification confirmed the differential expression of genes in the drought resistance signaling pathway. Our results provide valuable information for understanding potato drought-resistance mechanisms, and also enrich the gene resources available for drought-resistant potato breeding.

Mapping and QTL Analysis of Early-Maturity Traits in Tetraploid Potato (Solanum tuberosum L.).

Early maturity is one of the most important agronomical traits in potato breeding. To identify the DNA segment that codes for early maturity, a tetraploid potato segregation population of "Zhongshu 19" × "Zhongshu 3" was genetically analyzed, using a combination of high throughput simplified genome sequencing (2b-RAD) and bulked segregant analysis (BSA). The DNA segment related to the early-maturity trait was identified at the 3.7~4.2 Mb locus on the short arm of chromosome 5. Eight molecular markers were developed, of which five were closely linked to the early-maturity trait loci. Additionally, 42 simple sequence repeats (SSR) markers were constructed based on the reference sequence of Solanum tuberosum group Phureja DM1-3 516 R44 (DM). Using the TetraploidMap software, the linkage map of chromosome 5 was constructed with 50 markers. The total map length was 172 centiMorgan (cM), with an average genetic distance of 3.44 cM. Correlating molecular and phenotypic data of the segregating population, the mapped Quantitative Trait Loci (QTL) on the short arm of chromosome 5 contributed to 33.55% of the early-maturity phenotype. The early-maturity QTL was located at 84 cM, flanked by the SSR5-85-1 and SCAR5-8 markers. The QTL was fine-mapped to 471 kb. Using DNA sequence annotation, 34 genes were identified in this region, 12 of them with unknown function. Among the other 22 annotated genes, E3 ubiquitin ligase gene PUB14 could be related to maturity and regulation of tuber formation. The constructed QTL map is a useful basic tool for the cloning of early-maturity related genes in tetraploid potatoes.

Gene Profiling in Late Blight Resistance in Potato Genotype SD20.

Late blight caused by the oomycete fungus Phytophthora infestans (Pi) is the most serious obstacle to potato (Solanum tuberosum) production in the world. A super race isolate, CN152, which was identified from Sichuan Province, China, could overcome nearly all known late blight resistance genes and caused serious damage in China. The potato genotype SD20 was verified to be highly resistant to CN152; however, the molecular regulation network underlying late blight resistance pathway remains unclear in SD20. Here, we performed a time-course experiment to systematically profile the late blight resistance response genes using RNA-sequencing in SD20. We identified 3354 differentially expressed genes (DEGs), which mainly encoded transcription factors and protein kinases, and also included four NBS-LRR genes. The late blight responsive genes showed time-point-specific induction/repression. Multi-signaling pathways of salicylic acid, jasmonic acid, and ethylene signaling pathways involved in resistance and defense against Pi in SD20. Gene Ontology and KEGG analyses indicated that the DEGs were significantly enriched in metabolic process, protein serine/threonine kinase activity, and biosynthesis of secondary metabolites. Forty-three DEGs were involved in immune response, of which 19 were enriched in hypersensitive response reaction, which could play an important role in broad-spectrum resistance to Pi infection. Experimental verification confirmed the induced expression of the responsive genes in the late blight resistance signaling pathway, such as WRKY, ERF, MAPK, and NBS-LRR family genes. Our results provided valuable information for understanding late blight resistance mechanism of potato.

The Ph-3 gene from Solanum pimpinellifolium encodes CC-NBS-LRR protein conferring resistance to Phytophthora infestans.

KEY MESSAGE: Ph-3 is the first cloned tomato gene for resistance to late blight and encodes a CC-NBS-LRR protein. Late blight, caused by Phytophthora infestans, is one of the most destructive diseases in tomato. The resistance (R) gene Ph-3, derived from Solanum pimpinellifolium L3708, provides resistance to multiple P. infestans isolates and has been widely used in tomato breeding programmes. In our previous study, Ph-3 was mapped into a region harbouring R gene analogues (RGA) at the distal part of long arm of chromosome 9. To further narrow down the Ph-3 interval, more recombinants were identified using the flanking markers G2-4 and M8-2, which defined the Ph-3 gene to a 26 kb region according to the Heinz1706 reference genome. To clone the Ph-3 gene, a bacterial artificial chromosome (BAC) library was constructed using L3708 and one BAC clone B25E21 containing the Ph-3 region was identified. The sequence of the BAC clone B25E21 showed that only one RGA was present in the target region. A subsequent complementation analysis demonstrated that this RGA, encoding a CC-NBS-LRR protein, was able to complement the susceptible phenotype in cultivar Moneymaker. Thus this RGA was considered the Ph-3 gene. The predicted Ph-3 protein shares high amino acid identity with the chromosome-9-derived potato resistance proteins against P. infestans (Rpi proteins).

Construction of reference chromosome-scale pseudomolecules for potato: integrating the potato genome with genetic and physical maps.

The genome of potato, a major global food crop, was recently sequenced. The work presented here details the integration of the potato reference genome (DM) with a new sequence-tagged site marker-based linkage map and other physical and genetic maps of potato and the closely related species tomato. Primary anchoring of the DM genome assembly was accomplished by the use of a diploid segregating population, which was genotyped with several types of molecular genetic markers to construct a new ~936 cM linkage map comprising 2469 marker loci. In silico anchoring approaches used genetic and physical maps from the diploid potato genotype RH89-039-16 (RH) and tomato. This combined approach has allowed 951 superscaffolds to be ordered into pseudomolecules corresponding to the 12 potato chromosomes. These pseudomolecules represent 674 Mb (~93%) of the 723 Mb genome assembly and 37,482 (~96%) of the 39,031 predicted genes. The superscaffold order and orientation within the pseudomolecules are closely collinear with independently constructed high density linkage maps. Comparisons between marker distribution and physical location reveal regions of greater and lesser recombination, as well as regions exhibiting significant segregation distortion. The work presented here has led to a greatly improved ordering of the potato reference genome superscaffolds into chromosomal "pseudomolecules".

Genome sequence and analysis of the tuber crop potato.

Potato (Solanum tuberosum L.) is the world's most important non-grain food crop and is central to global food security. It is clonally propagated, highly heterozygous, autotetraploid, and suffers acute inbreeding depression. Here we use a homozygous doubled-monoploid potato clone to sequence and assemble 86% of the 844-megabase genome. We predict 39,031 protein-coding genes and present evidence for at least two genome duplication events indicative of a palaeopolyploid origin. As the first genome sequence of an asterid, the potato genome reveals 2,642 genes specific to this large angiosperm clade. We also sequenced a heterozygous diploid clone and show that gene presence/absence variants and other potentially deleterious mutations occur frequently and are a likely cause of inbreeding depression. Gene family expansion, tissue-specific expression and recruitment of genes to new pathways contributed to the evolution of tuber development. The potato genome sequence provides a platform for genetic improvement of this vital crop.