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.

Genotype-by-environment interactions and local adaptation shape selection in the US National Chip Processing Trial.

KEY MESSAGE: We find evidence of selection for local adaptation and extensive genotype-by-environment interaction in the potato National Chip Processing Trial (NCPT). We present a novel method for dissecting the interplay between selection, local adaptation and environmental response in plant breeding schemes. Balancing local adaptation and the desire for widely adapted cultivars is challenging for plant breeders and makes genotype-by-environment interactions (GxE) an important target of selection. Selecting for GxE requires plant breeders to evaluate plants across multiple environments. One way breeders have accomplished this is to test advanced materials across many locations. Public potato breeders test advanced breeding material in the National Chip Processing Trial (NCPT), a public-private partnership where breeders from ten institutions submit advanced chip lines to be evaluated in up to ten locations across the country. These clones are genotyped and phenotyped for important agronomic traits. We used these data to interrogate the NCPT for GxE. Further, because breeders submitting clones to the NCPT select in a relatively small geographic range for the first 3 years of selection, we examined these data for evidence of incidental selection for local adaptation, and the alleles underlying it, using an environmental genome-wide association study (envGWAS). We found genomic regions associated with continuous environmental variables and discrete breeding programs, as well as regions of the genome potentially underlying GxE for yield.

Molecular characterization for food safety assessment of a genetically modified late blight resistant potato: an unusual case.

Standard food safety assessments of genetically modified crops require a thorough molecular characterization of the novel DNA as inserted into the plant that is intended for commercialization, as well as a comparison of agronomic and nutritional characteristics of the genetically modified to the non-modified counterpart. These characterization data are used to identify any unintended changes in the inserted DNA or in the modified plant that would require assessment for safety in addition to the assessment of the intended modification. An unusual case of an unintended effect discovered from the molecular characterization of a genetically modified late blight resistant potato developed for growing in Bangladesh and Indonesia is presented here. Not only was a significant portion of the plasmid vector backbone DNA inserted into the plant along with the intended insertion of an R-gene for late blight resistance, but the inserted DNA was split into two separate fragments and inserted into two separate chromosomes. One fragment carries the R-gene and the other fragment carries the NPTII selectable marker gene and the plasmid backbone DNA. The implications of this for the food safety assessment of this late blight resistant potato are considered.

Genome diversity of tuber-bearing Solanum uncovers complex evolutionary history and targets of domestication in the cultivated potato.

Cultivated potatoes (Solanum tuberosum L.), domesticated from wild Solanum species native to the Andes of southern Peru, possess a diverse gene pool representing more than 100 tuber-bearing relatives (Solanum section Petota). A diversity panel of wild species, landraces, and cultivars was sequenced to assess genetic variation within tuber-bearing Solanum and the impact of domestication on genome diversity and identify key loci selected for cultivation in North and South America. Sequence diversity of diploid and tetraploid Stuberosum exceeded any crop resequencing study to date, in part due to expanded wild introgressions following polyploidy that captured alleles outside of their geographic origin. We identified 2,622 genes as under selection, with only 14-16% shared by North American and Andean cultivars, showing that a limited gene set drove early improvement of cultivated potato, while adaptation of upland (Stuberosum group Andigena) and lowland (S. tuberosum groups Chilotanum and Tuberosum) populations targeted distinct loci. Signatures of selection were uncovered in genes controlling carbohydrate metabolism, glycoalkaloid biosynthesis, the shikimate pathway, the cell cycle, and circadian rhythm. Reduced sexual fertility that accompanied the shift to asexual reproduction in cultivars was reflected by signatures of selection in genes regulating pollen development/gametogenesis. Exploration of haplotype diversity at potato's maturity locus (StCDF1) revealed introgression of truncated alleles from wild species, particularly Smicrodontum in long-day-adapted cultivars. This study uncovers a historic role of wild Solanum species in the diversification of long-day-adapted tetraploid potatoes, showing that extant natural populations represent an essential source of untapped adaptive potential.

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.

Genome Reduction Uncovers a Large Dispensable Genome and Adaptive Role for Copy Number Variation in Asexually Propagated Solanum tuberosum.

Clonally reproducing plants have the potential to bear a significantly greater mutational load than sexually reproducing species. To investigate this possibility, we examined the breadth of genome-wide structural variation in a panel of monoploid/doubled monoploid clones generated from native populations of diploid potato (Solanum tuberosum), a highly heterozygous asexually propagated plant. As rare instances of purely homozygous clones, they provided an ideal set for determining the degree of structural variation tolerated by this species and deriving its minimal gene complement. Extensive copy number variation (CNV) was uncovered, impacting 219.8 Mb (30.2%) of the potato genome with nearly 30% of genes subject to at least partial duplication or deletion, revealing the highly heterogeneous nature of the potato genome. Dispensable genes (>7000) were associated with limited transcription and/or a recent evolutionary history, with lower deletion frequency observed in genes conserved across angiosperms. Association of CNV with plant adaptation was highlighted by enrichment in gene clusters encoding functions for environmental stress response, with gene duplication playing a part in species-specific expansions of stress-related gene families. This study revealed unique impacts of CNV in a species with asexual reproductive habits and how CNV may drive adaption through evolution of key stress pathways.

Extensive genome heterogeneity leads to preferential allele expression and copy number-dependent expression in cultivated potato.

Relative to homozygous diploids, the presence of multiple homologs or homeologs in polyploids affords greater tolerance to mutations that can impact genome evolution. In this study, we describe sequence and structural variation in the genomes of six accessions of cultivated potato (Solanum tuberosum L.), a vegetatively propagated autotetraploid and their impact on the transcriptome. Sequence diversity was high with a mean single nucleotide polymorphisms (SNP) rate of approximately 1 per 50 bases suggestive of high levels of allelic diversity. Additive gene expression was observed in leaves (3605 genes) and tubers (6156 genes) that contrasted the preferential allele expression of between 2180 and 3502 and 3367 and 5270 genes in the leaf and tuber transcriptome, respectively. Preferential allele expression was significantly associated with evolutionarily conserved genes suggesting selection of specific alleles of genes responsible for biological processes common to angiosperms during the breeding selection process. Copy number variation was rampant with between 16 098 and 18 921 genes in each cultivar exhibiting duplication or deletion. Copy number variable genes tended to be evolutionarily recent, lowly expressed, and enriched in genes that show increased expression in response to biotic and abiotic stress treatments suggestive of a role in adaptation. Gene copy number impacts on gene expression were detected with 528 genes having correlations between copy number and gene expression. Collectively, these data suggest that in addition to allelic variation of coding sequence, the heterogenous nature of the tetraploid potato genome contributes to a highly dynamic transcriptome impacted by allele preferential and copy number-dependent expression effects.

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

BACKGROUND: Genome-wide single nucleotide polymorphism (SNP) markers coupled with allele dosage information has emerged as a powerful tool for studying complex traits in cultivated autotetraploid potato (Solanum tuberosum L., 2n = 4× = 48). To date, this approach has been effectively applied to the identification of quantitative trait loci (QTLs) underlying highly heritable traits such as disease resistance, but largely unexplored for traits with complex patterns of inheritance. RESULTS: In this study, an F1 tetraploid russet mapping population (162 individuals) was evaluated for multiple quantitative traits over two years and two locations to identify QTLs associated with tuber sugar concentration, processing quality, vine maturity, and other high-value agronomic traits. We report the linkage maps for the 12 potato chromosomes and the QTL location with corresponding genetic models and candidate SNPs explaining the highest phenotypic variation for tuber quality and maturity related traits. Significant QTLs for tuber glucose concentration and tuber fry color were detected on chromosomes 4, 5, 6, 10, and 11. Collectively, these QTLs explained between 24 and 46% of the total phenotypic variation for tuber glucose and fry color, respectively. The QTL on chromosome 10 was associated with apoplastic invertases, with 'Premier Russet' contributing the favorable allele for fry processing quality. On chromosome 5, minor-effect QTLs for tuber glucose concentration and fry color co-localized with various major-effect QTLs, including vine maturity, growth habit, tuber shape, early blight (Altenaria tenuis), and Verticillium wilt (Verticillium spp.). CONCLUSIONS: Linkage analysis and QTL mapping in a russet mapping population (A05141) using SNP dosage information successfully identified favorable alleles and candidate SNPs for resistance to the accumulation of tuber reducing sugars. These novel markers have a high potential for the improvement of tuber processing quality. Moreover, the discovery of different genetic models for traits with overlapping QTLs at the maturity locus clearly suggests an independent genetic control.

Automated tetraploid genotype calling by hierarchical clustering.

KEY MESSAGE: 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.

Treatment of potato tubers with the synthetic cytokinin 1-(α-ethylbenzyl)-3-nitroguanidine results in rapid termination of endodormancy and induction of transcripts associated with cell proliferation and growth.

Perennial plants undergo repression of meristematic activity in a process called dormancy. Dormancy is a complex metabolic process with implications for plant breeding and crop yield. Endodormancy, a specific subclass of dormancy, is characteristic of internal physiological mechanisms resulting in growth suppression. In this study, we examine transcriptional changes associated with the natural cessation of endodormancy in potato tuber meristems and in endodormant tubers treated with the cytokinin analog 1-(α-ethylbenzyl)-3-niroguanidine (NG), which terminates dormancy. RNA-sequencing was used to examine transcriptome changes between endodormant and non-dormant meristems from four different harvest years. A total of 35,091 transcripts were detected with 2132 differentially expressed between endodormant and non-dormant tuber meristems. Endodormant potato tubers were treated with the synthetic cytokinin NG and transcriptome changes analyzed using RNA-seq after 1, 4, and 7 days following NG exposure. A comparison of natural cessation of dormancy and NG-treated tubers demonstrated that by 4 days after NG exposure, potato meristems exhibited transcriptional profiles similar to the non-dormant state with elevated expression of multiple histones, a variety of cyclins, and other genes associated with proliferation and cellular replication. Three homologues encoding for CYCD3 exhibited elevated expression in both non-dormant and NG-treated potato tissues. These results suggest that NG terminates dormancy and induces expression cell cycle-associated transcripts within 4 days of treatment.

T-DNA characterization of genetically modified 3-R-gene late blight-resistant potato events with a novel procedure utilizing the Samplix Xdrop® enrichment technology.

Before the commercialization of genetically modified crops, the events carrying the novel DNA must be thoroughly evaluated for agronomic, nutritional, and molecular characteristics. Over the years, polymerase chain reaction-based methods, Southern blot, and short-read sequencing techniques have been utilized for collecting molecular characterization data. Multiple genomic applications are necessary to determine the insert location, flanking sequence analysis, characterization of the inserted DNA, and determination of any interruption of native genes. These techniques are time-consuming and labor-intensive, making it difficult to characterize multiple events. Current advances in sequencing technologies are enabling whole-genomic sequencing of modified crops to obtain full molecular characterization. However, in polyploids, such as the tetraploid potato, it is a challenge to obtain whole-genomic sequencing coverage that meets the regulatory approval of the genetic modification. Here we describe an alternative to labor-intensive applications with a novel procedure using Samplix Xdrop® enrichment technology and next-generation Nanopore sequencing technology to more efficiently characterize the T-DNA insertions of four genetically modified potato events developed by the Feed the Future Global Biotech Potato Partnership: DIA_MSU_UB015, DIA_MSU_UB255, GRA_MSU_UG234, and GRA_MSU_UG265 (derived from regionally important varieties Diamant and Granola). Using the Xdrop® /Nanopore technique, we obtained a very high sequence read coverage within the T-DNA and junction regions. In three of the four events, we were able to use the data to confirm single T-DNA insertions, identify insert locations, identify flanking sequences, and characterize the inserted T-DNA. We further used the characterization data to identify native gene interruption and confirm the stability of the T-DNA across clonal cycles. These results demonstrate the functionality of using the Xdrop® /Nanopore technique for T-DNA characterization. This research will contribute to meeting regulatory safety and regulatory approval requirements for commercialization with small shareholder farmers in target countries within our partnership.

A survey of the Sli gene in wild and cultivated potato.

Inbred-hybrid breeding of diploid potatoes necessitates breeding lines that are self-compatible. One way of incorporating self-compatibility into incompatible cultivated potato (Solanum tuberosum) germplasm is to introduce the S-locus inhibitor gene (Sli), which functions as a dominant inhibitor of gametophytic self-incompatibility. To learn more about Sli diversity and function in wild species relatives of cultivated potato, we obtained Sli gene sequences that extended from the 5'UTR to the 3'UTR from 133 individuals from 22 wild species relatives of potato and eight diverse cultivated potato clones. DNA sequence alignment and phylogenetic trees based on genomic and protein sequences show that there are two highly conserved groups of Sli sequences. DNA sequences in one group contain the 533 bp insertion upstream of the start codon identified previously in self-compatible potato. The second group lacks the insertion. Three diploid and four polyploid individuals of wild species collected from geographically disjointed localities contained Sli with the 533 bp insertion. For most of the wild species clones examined, however, Sli did not have the insertion. Phylogenetic analysis indicated that Sli sequences with the insertion, in wild species and in cultivated clones, trace back to a single origin. Some diploid wild potatoes that have Sli with the insertion were self-incompatible and some wild potatoes that lack the insertion were self-compatible. Although there is evidence of positive selection for some codon positions in Sli, there is no evidence of diversifying selection at the gene level. In silico analysis of Sli protein structure did not support the hypothesis that amino acid changes from wild-type (no insertion) to insertion-type account for changes in protein function. Our study demonstrated that genetic factors besides the Sli gene must be important for conditioning a switch in the mating system from self-incompatible to self-compatible in wild potatoes.

Segregation and efficacy of the cry1Ia1 gene for control of potato tuberworm in four populations of cultivated potato.

The potato (Solanum tuberosum L.) cultivar 'SpuntaG2' contains a single copy of the Bacillus thuringiensis (Bt) cry1Ia1 gene and controls potato tuberworm (Phthorimaea operculella Zeller, Lepidoptera: Gelechiidae). Two potato cultivars and two breeding lines were crossed with SpuntaG2 creating four populations used to study cry1Ia1 segregation and efficacy. The cry1Ia1 gene segregated in each of the four populations with a 1:1 ratio. All progeny that were polymerase chain reaction positive for the cry1Ia1 gene had no surviving larvae and no leaf mining in detached leaf assays after 72 h. These results support previous evidence that SpuntaG2 carries a single copy of the cry1Ia1 gene and that transmission of the transgene from parent to progeny is not restricted and follows expected Mendelian segregation ratios. Based on detached leaf assays, the efficacy of the cry1Ia1 gene is retained through sexual transmission. If the SpuntaG2 cry1Ia1 insertion event is deregulated for commercial use, SpuntaG2 could be used for conventional breeding and the progeny carrying the SpuntaG2 event would also be available for commercial use.

The Effect of Self-Compatibility Factors on Interspecific Compatibility in Solanum Section Petota.

The relationships of interspecific compatibility and incompatibility in Solanum section Petota are complex. Inquiry into these relationships in tomato and its wild relatives has elucidated the pleiotropic and redundant function of S-RNase and HT which tandemly and independently mediate both interspecific and intraspecific pollen rejection. Our findings presented here are consistent with previous work conducted in Solanum section Lycopersicon showing that S-RNase plays a central role in interspecific pollen rejection. Statistical analyses also demonstrated that HT-B alone is not a significant factor in these pollinations; demonstrating the overlap in gene function between HT-A and HT-B, as HT-A, was present and functional in all genotypes used. We were not able to replicate the general absence of prezygotic stylar barriers observable in S. verrucosum, which has been attributed to the lack of S-RNase, indicating that other non-S-RNase factors play a significant role. We also demonstrated that Sli played no significant role in these interspecific pollinations, directly conflicting with previous research. It is possible that S. chacoense as a pollen donor is better able to bypass stylar barriers in 1EBN species such as S. pinnatisectum. Consequently, S. chacoense may be a valuable resource in accessing these 1EBN species regardless of Sli status.

Genome Report: Genome sequence of 1S1, a transformable and highly regenerable diploid potato for use as a model for gene editing and genetic engineering.

Availability of readily transformable germplasm, as well as efficient pipelines for gene discovery are notable bottlenecks in the application of genome editing in potato. To study and introduce traits such as resistance against biotic and abiotic factors, tuber quality traits and self-fertility, model germplasm that is amenable to gene editing and regeneration is needed. Cultivated potato is a heterozygous autotetraploid and its genetic redundancy and complexity makes studying gene function challenging. Genome editing is simpler at the diploid level, with fewer allelic variants to consider. A readily transformable diploid potato would be further complemented by genomic resources that could aid in high throughput functional analysis. The heterozygous Solanum tuberosum Group Phureja clone 1S1 has a high regeneration rate, self-fertility, desirable tuber traits and is amenable to Agrobacterium-mediated transformation. We leveraged its amenability to Agrobacterium-mediated transformation to create a Cas9 constitutively expressing line for use in viral vector-based gene editing. To create a contiguous genome assembly, a homozygous doubled monoploid of 1S1 (DM1S1) was sequenced using 44 Gbp of long reads generated from Oxford Nanopore Technologies (ONT), yielding a 736 Mb assembly that encoded 31,145 protein-coding genes. The final assembly for DM1S1 represents a nearly complete genic space, shown by the presence of 99.6% of the genes in the Benchmarking Universal Single Copy Orthologs (BUSCO) set. Variant analysis with Illumina reads from 1S1 was used to deduce its alternate haplotype. These genetic and genomic resources provide a toolkit for applications of genome editing in both basic and applied research of potato.

HT-B and S-RNase CRISPR-Cas9 double knockouts show enhanced self-fertility in diploid Solanum tuberosum.

The Gametophytic Self-Incompatibility (GSI) system in diploid potato (Solanum tuberosum L.) poses a substantial barrier in diploid potato breeding by hindering the generation of inbred lines. One solution is gene editing to generate self-compatible diploid potatoes which will allow for the generation of elite inbred lines with fixed favorable alleles and heterotic potential. The S-RNase and HT genes have been shown previously to contribute to GSI in the Solanaceae family and self-compatible S. tuberosum lines have been generated by knocking out S-RNase gene with CRISPR-Cas9 gene editing. This study employed CRISPR-Cas9 to knockout HT-B either individually or in concert with S-RNase in the diploid self-incompatible S. tuberosum clone DRH-195. Using mature seed formation from self-pollinated fruit as the defining characteristic of self-compatibility, HT-B-only knockouts produced little or no seed. In contrast, double knockout lines of HT-B and S-RNase displayed levels of seed production that were up to three times higher than observed in the S-RNase-only knockout, indicating a synergistic effect between HT-B and S-RNase in self-compatibility in diploid potato. This contrasts with compatible cross-pollinations, where S-RNase and HT-B did not have a significant effect on seed set. Contradictory to the traditional GSI model, self-incompatible lines displayed pollen tube growth reaching the ovary, yet ovules failed to develop into seeds indicating a potential late-acting self-incompatibility in DRH-195. Germplasm generated from this study will serve as a valuable resource for diploid potato breeding.

Cold stress induces differential gene expression of retained homeologs in Camelina sativa cv Suneson.

Camelina sativa (L.) Crantz, a member of the Brassicaceae, has potential as a biofuel feedstock which is attributable to the production of fatty acids in its seeds, its fast growth cycle, and low input requirements. While a genome assembly is available for camelina, it was generated from short sequence reads and is thus highly fragmented in nature. Using long read sequences, we generated a chromosome-scale, highly contiguous genome assembly (644,491,969 bp) for the spring biotype cultivar 'Suneson' with an N50 contig length of 12,031,512 bp and a scaffold N50 length of 32,184,682 bp. Annotation of protein-coding genes revealed 91,877 genes that encode 133,355 gene models. We identified a total of 4,467 genes that were significantly up-regulated under cold stress which were enriched in gene ontology terms associated with "response to cold" and "response to abiotic stress". Coexpression analyses revealed multiple coexpression modules that were enriched in genes differentially expressed following cold stress that had putative functions involved in stress adaptation, specifically within the plastid. With access to a highly contiguous genome assembly, comparative analyses with Arabidopsis thaliana revealed 23,625 A. thaliana genes syntenic with 45,453 Suneson genes. Of these, 24,960 Suneson genes were syntenic to 8,320 A. thaliana genes reflecting a 3 camelina homeolog to 1 Arabidopsis gene relationship and retention of all three homeologs. Some of the retained triplicated homeologs showed conserved gene expression patterns under control and cold-stressed conditions whereas other triplicated homeologs displayed diverged expression patterns revealing sub- and neo-functionalization of the homeologs at the transcription level. Access to the chromosome-scale assembly of Suneson will enable both basic and applied research efforts in the improvement of camelina as a sustainable biofuel feedstock.

Single nucleotide polymorphism discovery in elite North American potato germplasm.

BACKGROUND: Current breeding approaches in potato rely almost entirely on phenotypic evaluations; molecular markers, with the exception of a few linked to disease resistance traits, are not widely used. Large-scale sequence datasets generated primarily through Sanger Expressed Sequence Tag projects are available from a limited number of potato cultivars and access to next generation sequencing technologies permits rapid generation of sequence data for additional cultivars. When coupled with the advent of high throughput genotyping methods, an opportunity now exists for potato breeders to incorporate considerably more genotypic data into their decision-making. RESULTS: To identify a large number of Single Nucleotide Polymorphisms (SNPs) in elite potato germplasm, we sequenced normalized cDNA prepared from three commercial potato cultivars: 'Atlantic', 'Premier Russet' and 'Snowden'. For each cultivar, we generated 2 Gb of sequence which was assembled into a representative transcriptome of ~28-29 Mb for each cultivar. Using the Maq SNP filter that filters read depth, density, and quality, 575,340 SNPs were identified within these three cultivars. In parallel, 2,358 SNPs were identified within existing Sanger sequences for three additional cultivars, 'Bintje', 'Kennebec', and 'Shepody'. Using a stringent set of filters in conjunction with the potato reference genome, we identified 69,011 high confidence SNPs from these six cultivars for use in genotyping with the Infinium platform. Ninety-six of these SNPs were used with a BeadXpress assay to assess allelic diversity in a germplasm panel of 248 lines; 82 of the SNPs proved sufficiently informative for subsequent analyses. Within diverse North American germplasm, the chip processing market class was most distinct, clearly separated from all other market classes. The round white and russet market classes both include fresh market and processing cultivars. Nevertheless, the russet and round white market classes are more distant from each other than processing are from fresh market types within these two groups. CONCLUSIONS: The genotype data generated in this study, albeit limited in number, has revealed distinct relationships among the market classes of potato. The SNPs identified in this study will enable high-throughput genotyping of germplasm and populations, which in turn will enable more efficient marker-assisted breeding efforts in potato.

Genetic diversity and population structure of advanced clones selected over forty years by a potato breeding program in the USA.

Knowledge regarding genetic diversity and population structure of breeding materials is essential for crop improvement. The Texas A&M University Potato Breeding Program has a collection of advanced clones selected and maintained in-vitro over a 40-year period. Little is known about its genetic makeup and usefulness for the current breeding program. In this study, 214 potato clones were genotyped with the Infinium Illumina 22 K V3 Potato Array. After filtering, a total of 10,106 single nucleotide polymorphic (SNP) markers were used for analysis. Heterozygosity varied by SNP, with an overall average of 0.59. Three groups of tetraploid clones primarily based on potato market classes, were detected using STRUCTURE software and confirmed by discriminant analysis of principal components. The highest coefficient of differentiation observed between the groups was 0.14. Signatures of selection were uncovered in genes controlling potato flesh and skin color, length of plant cycle and tuberization, and carbohydrate metabolism. A core set of 43 clones was obtained using Core Hunter 3 to develop a sub-collection that retains similar genetic diversity as the whole population, minimize redundancies, and facilitates long-term conservation of genetic resources. The comprehensive molecular characterization of our breeding clone bank collection contributes to understanding the genetic diversity of existing potato resources. This analysis could be applied to other breeding programs and assist in the selection of parents, fingerprinting, protection, and management of the breeding collections.