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.

Potato (Solanum tuberosum L.) originates from the Andes and evolved short-day-dependent tuber formation as a vegetative propagation strategy. Here we describe the identification of a central regulator underlying a major-effect quantitative trait locus for plant maturity and initiation of tuber development. We show that this gene belongs to the family of DOF (DNA-binding with one finger) transcription factors and regulates tuberization and plant life cycle length, by acting as a mediator between the circadian clock and the StSP6A mobile tuberization signal. We also show that natural allelic variants evade post-translational light regulation, allowing cultivation outside the geographical centre of origin of potato. Potato is a member of the Solanaceae family and is one of the world's most important food crops. This annual plant originates from the Andean regions of South America. Potato develops tubers from underground stems called stolons. Its equatorial origin makes potato essentially short-day dependent for tuberization and potato will not make tubers in the long-day conditions of spring and summer in the northern latitudes. When introduced in temperate zones, wild material will form tubers in the course of the autumnal shortening of day-length. Thus, one of the first selected traits in potato leading to a European potato type is likely to have been long-day acclimation for tuberization. Potato breeders can exploit the naturally occurring variation in tuberization onset and life cycle length, allowing varietal breeding for different latitudes, harvest times and markets.

Assessment of genomic DNA sequence variation and genotype calling in autotetraploids implies the ability to distinguish among five possible alternative allele copy number states. This study demonstrates the accuracy of genotyping-by-sequencing (GBS) of a large collection of autotetraploid potato cultivars using next-generation sequencing. It is still costly to reach sufficient read depths on a genome wide scale, across the cultivated gene pool. Therefore, we enriched cultivar-specific DNA sequencing libraries using an in-solution hybridisation method (SureSelect). This complexity reduction allowed to confine our study to 807 target genes distributed across the genomes of 83 tetraploid cultivars and one reference (DM 1–3 511). Indexed sequencing libraries were paired-end sequenced in 7 pools of 12 samples using Illumina HiSeq2000. After filtering and processing the raw sequence data, 12.4 Gigabases of high-quality sequence data was obtained, which mapped to 2.1 Mb of the potato reference genome, with a median average read depth of 63× per cultivar. We detected 129,156 sequence variants and genotyped the allele copy number of each variant for every cultivar. In this cultivar panel a variant density of 1 SNP/24 bp in exons and 1 SNP/15 bp in introns was obtained. The average minor allele frequency (MAF) of a variant was 0.14. Potato germplasm displayed a large number of relatively rare variants and/or haplotypes, with 61% of the variants having a MAF below 0.05. A very high average nucleotide diversity (π = 0.0107) was observed. Nucleotide diversity varied among potato chromosomes. Several genes under selection were identified. Genotyping-by-sequencing results, with allele copy number estimates, were validated with a KASP genotyping assay. This validation showed that read depths of ∼60–80× can be used as a lower boundary for reliable assessment of allele copy number of sequence variants in autotetraploids. Genotypic data were associated with traits, and alleles strongly influencing maturity and flesh colour were identified.