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.

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”.

The majority of disease resistance (R) genes identified to date in plants encode a nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domain containing protein. Additional domains such as coiled-coil (CC) and TOLL/interleukin-1 receptor (TIR) domains can also be present. In the recently sequenced Solanum tuberosum group phureja genome we used HMM models and manual curation to annotate 435 NBS-encoding R gene homologs and 142 NBS-derived genes that lack the NBS domain. Highly similar homologs for most previously documented Solanaceae R genes were identified. A surprising ∼41% (179) of the 435 NBS-encoding genes are pseudogenes primarily caused by premature stop codons or frameshift mutations. Alignment of 81.80% of the 577 homologs to S. tuberosum group phureja pseudomolecules revealed non-random distribution of the R-genes; 362 of 470 genes were found in high density clusters on 11 chromosomes.

Naturally occurring autopolyploid species, such as the autotetraploid potato Solanum tuberosum, face a variety of challenges during meiosis. These include proper pairing, recombination and correct segregation of multiple homologous chromosomes, which can form complex multivalent configurations at metaphase I, and in turn alter allelic segregation ratios through double reduction. Here, we present a reference map of meiotic stages in diploid and tetraploid S. tuberosum using fluorescence in situ hybridisation (FISH) to differentiate individual meiotic chromosomes 1 and 2. A diploid-like behaviour at metaphase I involving bivalent configurations was predominant in all three tetraploid varieties. The crossover frequency per bivalent was significantly reduced in the tetraploids compared with a diploid variety, which likely indicates meiotic adaptation to the autotetraploid state. Nevertheless, bivalents were accompanied by a substantial frequency of multivalents, which varied by variety and by chromosome (7–48%). We identified possible sites of synaptic partner switching, leading to multivalent formation, and found potential defects in the polymerisation and/or maintenance of the synaptonemal complex in tetraploids. These findings demonstrate the rise of S. tuberosum as a model for autotetraploid meiotic recombination research and highlight constraints on meiotic chromosome configurations and chiasma frequencies as an important feature of an evolved autotetraploid meiosis.

No abstract

Potato is a drought-sensitive crop whose global sustainable production is threatened by alterations in water availability. Whilst ancestral Solanum tuberosum Andigenum landraces retain wild drought tolerance mechanisms, their molecular bases remain poorly understood. In this study, an aeroponic growth system was established to investigate stress responses in leaf and root of two Andigenum varieties with contrasting drought tolerance. Comparative transcriptome analysis revealed widespread differences in the response of the two varieties at early and late time points of exposure to drought stress and in the recovery after rewatering. Major differences in the response of the two varieties occurred at the early time point, suggesting the speed of response is crucial. In the leaves and roots of the tolerant variety, we observed rapid upregulation of ABA-related genes, which did not occur until later in the susceptible variety and indicated not only more effective ABA synthesis and mobilization, but more effective feedback regulation to limit detrimental effects of too much ABA. Roots of both varieties showed differential expression of genes involved in cell wall reinforcement and remodeling to maintain cell wall strength, hydration and growth under drought stress, including genes involved in lignification and wall expansion, though the response was stronger in the tolerant variety. Such changes in leaf and root may help to limit water losses in the tolerant variety, while limiting the reduction in photosynthetic rate. These findings provide insights into molecular bases of drought tolerance mechanisms and pave the way for their reintroduction into modern cultivars with improved resistance to drought stress and yield stability under drought conditions.

Globodera pallida is a white potato cyst nematode present in the Andes, which causes huge losses to Peruvian farmers. An RNA-seq analysis allowed the identification of candidate genes that could mediate resistance against this pathogen. Two varieties, “María Huanca” (Solanum andigena) clone resistant (CIP 279142.12) and “Chimbina Colorada” (Solanum chaucha) (CIP 701013) clone susceptible to G. pallida, were used to identify differentially expressed genes. Total RNA from roots was extracted 72 hours post inoculation with second stage juveniles. Sequencing was done using the Illumina Hiseq 2500 platform. Reads were screened for quality issues and then mapped to the reference potato genome (clone DM1-3516 R44 v4.03). Here, we report 27717 and 27750 genes expressed in the resistant and susceptible variety respectively. The comparative analysis of expression identified 100 candidate genes. 91 genes were associated with resistance to G. pallida with Fold Change ≥ 2 (p <0.05). The remaining 9 R genes had Fold Change ≤ 1. We show differences in the expression of an NBS-LRR protein similar to Gro1-8, genes linked to late blight and TMV virus resistance.

e361 size ranges can easily be pooled and loaded in a single capillary. However, the resulting complexity of the pooling schemes and the amount of data generated require an automated workflow for acquisition and processing of fragment files. Moreover, it is well-known that fragment sizes, compared between two genotyping laboratories using capillary electrophoresis (CE), can differ considerably because different laboratories may use different CE machines and running conditions.Methods: Data processing: The software handles any pooling strategy (combination of dyes and size-compatible VNTRs) which can be parsed from the file name. Thus, each VNTR is defined by a pool, dye and expected size range, defined by the repeat length, offset and expected copy range. Using that information the software will automatically screen (in batch) all theoretical ranges for each VNTR and report on the peak(s) found (or not) in those ranges.Results: The resulting VNTR information is stored in integer-type character sets where each VNTR represents one character. For comparisons within the same laboratory or between laboratories with compatible instruments and procedures, the software allows the user to handle calculations based on an expected band size range, depending on the known offset, repeat length, expected copy number variation and a user-defined tolerance. In order to deal with a possible experimental size shift linked to the CE system used, a custom 'mapping' tool was developed, allowing observed sizes for a specific instrument to be mapped to real sizes and exact copy numbers.Conclusion: Data analysis: VNTR data can be analyzed as categorical characters (each different copy number is a different allele) or as quantitative characters. In the latter case, the larger the difference between copy numbers, the less related the organisms are considered. The Minimum Spanning Tree algorithm applied on VNTR data in BioNumerics has proven to be extremely useful for epidemiological study and population genetics.