Genome sequencing of adapted diploid potato clones.

Cultivated potato is a vegetatively propagated crop, and most varieties are autotetraploid with high levels of heterozygosity. Reducing the ploidy and breeding potato at the diploid level can increase efficiency for genetic improvement including greater ease of introgression of diploid wild relatives and more efficient use of genomics and markers in selection. More recently, selfing of diploids for generation of inbred lines for F1 hybrid breeding has had a lot of attention in potato. The current study provides genomics resources for nine legacy non-inbred adapted diploid potato clones developed at Agriculture and Agri-Food Canada. De novo genome sequence assembly using 10× Genomics and Illumina sequencing technologies show the genome sizes ranged from 712 to 948 Mbp. Structural variation was identified by comparison to two references, the potato DMv6.1 genome and the phased RHv3 genome, and a k-mer based analysis of sequence reads showed the genome heterozygosity range of 1 to 9.04% between clones. A genome-wide approach was taken to scan 5 Mb bins to visualize patterns of heterozygous deleterious alleles. These were found dispersed throughout the genome including regions overlapping segregation distortions. Novel variants of the StCDF1 gene conferring earliness of tuberization were found among these clones, which all produce tubers under long days. The genomes will be useful tools for genome design for potato breeding.

Identification of QTLs Controlling Resistance to Anthracnose Disease in Water Yam (Dioscorea alata).

Anthracnose disease caused by a fungus Colletotrichum gloeosporioides is the primary cause of yield loss in water yam (Dioscorea alata), the widely cultivated species of yam. Resistance to yam anthracnose disease (YAD) is a prime target in breeding initiatives to develop durable-resistant cultivars for sustainable management of the disease in water yam cultivation. This study aimed at tagging quantitative trait loci (QTL) for anthracnose disease resistance in a bi-parental mapping population of D. alata. Parent genotypes and their recombinant progenies were genotyped using the Genotyping by Sequencing (GBS) platform and phenotyped in two crop cycles for two years. A high-density genetic linkage map was built with 3184 polymorphic Single Nucleotide Polymorphism (NSP) markers well distributed across the genome, covering 1460.94 cM total length. On average, 163 SNP markers were mapped per chromosome with 0.58 genetic distances between SNPs. Four QTL regions related to yam anthracnose disease resistance were identified on three chromosomes. The proportion of phenotypic variance explained by these QTLs ranged from 29.54 to 39.40%. The QTL regions identified showed genes that code for known plant defense responses such as GDSL-like Lipase/Acylhydrolase, Protein kinase domain, and F-box protein. The results from the present study provide valuable insight into the genetic architecture of anthracnose resistance in water yam. The candidate markers identified herewith form a relevant resource to apply marker-assisted selection as an alternative to a conventional labor-intensive screening for anthracnose resistance in water yam.

Potato Response to Drought Stress: Physiological and Growth Basis.

Drought poses a major challenge to the production of potatoes worldwide. Climate change is predicted to further aggravate this challenge by intensifying potato crop exposure to increased drought severity and frequency. There is an ongoing effort to adapt our production systems of potatoes through the development of drought-tolerant cultivars that are appropriately engineered for the changing environment. The breeding of drought-tolerant cultivars can be approached through the identification of drought-related physiological and biochemical traits and their deployment in new potato cultivars. Thus, the main objective of this study was to develop a method to identify and characterize the drought-tolerant potato genotypes and the related key traits. To achieve this objective, first we studied 56 potato genotypes including 54 cultivars and 2 advanced breeding lines to assess drought tolerance in terms of tuber yield in the greenhouse experiment. Drought differentially reduced tuber yield in all genotypes. Based on their capacity to maintain percent tuber yield under drought relative to their well-watered controls, potato genotypes differed in their ability to tolerate drought. We then selected six genotypes, Bannock Russet, Nipigon, Onaway, Denali, Fundy, and Russet Norkotah, with distinct yield responses to drought to further examine the physiological and biochemical traits governing drought tolerance. The drought-induced reduction in tuber yield was only 15-20% for Bannock Russet and Nipigon, 44-47% for Onaway and Denali, and 83-91% for Fundy and Russet Norkotah. The tolerant genotypes, Bannock Russet and Nipigon, exhibited about a 2-3-fold increase in instantaneous water-use efficiency (WUE) under drought as compared with their well-watered controls. This stimulation was about 1.8-2-fold for moderately tolerant genotypes, Onaway and Denali, and only 1.5-fold for sensitive genotypes, Fundy, and Russet Norkotah. The differential stimulation of instantaneous WUE of tolerant and moderately tolerant genotypes vs. sensitive genotypes was accounted for by the differential suppression of the rates of photosynthesis, stomatal conductance, and transpiration rates across genotypes. Potato genotypes varied in their response to leaf protein content under drought. We suggest that the rates of photosynthesis, instantaneous WUE, and leaf protein content can be used as the selection criteria for the drought-tolerant potato genotypes.

Genetic parameters, prediction, and selection in a white Guinea yam early-generation breeding population using pedigree information.

Better understanding of the genetic control of traits in breeding populations is crucial for the selection of superior varieties and parents. This study aimed to assess genetic parameters and breeding values for six essential traits in a white Guinea yam (Dioscorea rotundata Poir.) breeding population. For this, pedigree-based best linear unbiased prediction (P-BLUP) was used. The results revealed significant nonadditive genetic variances and medium to high (.45-.79) broad-sense heritability estimates for the traits studied. The pattern of associations among the genetic values of the traits suggests that selection based on a multiple-trait selection index has potential for identifying superior breeding lines. Parental breeding values predicted using progeny performance identified 13 clones with high genetic potential for simultaneous improvement of the measured traits in the yam breeding program. Subsets of progeny were identified for intermating or further variety testing based on additive genetic and total genetic values. Selection of the top 5% progenies based on the multi-trait index revealed positive genetic gains for fresh tuber yield (t ha-1), tuber yield (kg plant-1), and average tuber weight (kg). However, genetic gain was negative for tuber dry matter content and Yam mosaic virus resistance in comparison with standard varieties. Our results show the relevance of P-BLUP for the selection of superior parental clones and progenies with higher breeding values for interbreeding and higher genotypic value for variety development in yam.

Genome analyses reveal the hybrid origin of the staple crop white Guinea yam (Dioscorea rotundata).

White Guinea yam (Dioscorea rotundata) is an important staple tuber crop in West Africa. However, its origin remains unclear. In this study, we resequenced 336 accessions of white Guinea yam and compared them with the sequences of wild Dioscorea species using an improved reference genome sequence of D. rotundata In contrast to a previous study suggesting that D. rotundata originated from a subgroup of Dioscorea praehensilis, our results suggest a hybrid origin of white Guinea yam from crosses between the wild rainforest species D. praehensilis and the savannah-adapted species Dioscorea abyssinica We identified a greater genomic contribution from D. abyssinica in the sex chromosome of Guinea yam and extensive introgression around the SWEETIE gene. Our findings point to a complex domestication scenario for Guinea yam and highlight the importance of wild species as gene donors for improving this crop through molecular breeding.

Comparative assessment of genetic diversity matrices and clustering methods in white Guinea yam (Dioscorea rotundata) based on morphological and molecular markers.

Understanding the diversity and genetic relationships among and within crop germplasm is invaluable for genetic improvement. This study assessed genetic diversity in a panel of 173 D. rotundata accessions using joint analysis for 23 morphological traits and 136,429 SNP markers from the whole-genome resequencing platform. Various diversity matrices and clustering methods were evaluated for a comprehensive characterization of genetic diversity in white Guinea yam from West Africa at phenotypic and molecular levels. The translation of the different diversity matrices from the phenotypic and genomic information into distinct groups varied with the hierarchal clustering methods used. Gower distance matrix based on phenotypic data and identity by state (IBS) distance matrix based on SNP data with the UPGMA clustering method found the best fit to dissect the genetic relationship in current set materials. However, the grouping pattern was inconsistent (r = - 0.05) between the morphological and molecular distance matrices due to the non-overlapping information between the two data types. Joint analysis for the phenotypic and molecular information maximized a comprehensive estimate of the actual diversity in the evaluated materials. The results from our study provide valuable insights for measuring quantitative genetic variability for breeding and genetic studies in yam and other root and tuber crops.

Paternity Assignment in White Guinea Yam (Dioscorea Rotundata) Half-Sib Progenies from Polycross Mating Design Using SNP Markers.

White Guinea yam is mostly a dioecious outcrossing crop with male and female flowers produced on distinct plants. Fertile parents produce high fruit set in an open pollination polycross block, which is a cost-effective and convenient way of generating variability in yam breeding. However, the pollen parent of progeny from polycross mating is usually unknown. This study aimed to determine paternity in white Guinea yam half-sib progenies from polycross mating design. A total of 394 half-sib progenies from random open pollination involving nine female and three male parents was genotyped with 6602 SNP markers from DArTSeq platform to recover full pedigree. A higher proportion of expected heterozygosity, allelic richness, and evenness were observed in the half-sib progenies. A complete pedigree was established for all progenies from two families (TDr1685 and TDr1688) with 100% accuracy, while in the remaining families, paternity was assigned successfully only for 56 to 98% of the progenies. Our results indicated unequal paternal contribution under natural open pollination in yam, suggesting unequal pollen migrations or gene flow among the crossing parents. A total of 3.8% of progenies lacking paternal identity due to foreign pollen contamination outside the polycross block was observed. This study established the efficient determination of parental reconstruction and allelic contributions in the white Guinea yam half-sib progenies generated from open pollination polycross using SNP markers. Findings are useful for parental reconstruction, accurate dissection of the genetic effects, and selection in white Guinea yam breeding program utilizing polycross mating design.

Potato Tuber Necrosis Induced by Alfalfa Mosaic Virus Depends on Potato Cultivar Rather Than on Virus Strain.

Alfalfa mosaic virus (AMV) was identified as the causal agent of internal tuber necrosis in the potato cultivar Innovator in New Brunswick, Canada. Further pathological characterization of the isolate (designated as isolate CaM) was performed on six potato cultivars and one breeding clone. Upon mechanical inoculation, four cultivars (Innovator, Yukon Gold, Rochdale Gold-Dorée, and Shepody) showed needle-sized necrotic spots and increasing calico symptoms on new leaves, whereas the remaining cultivars only developed calico symptoms on new leaves. All tubers of CaM-infected Innovator and Shepody plants developed sporadic internal necrotic spots, as did ca. 23 and 8% tubers of CaM-infected Yukon Gold and Rochdale Gold-Dorée, respectively. Sequence analysis of the CP gene of CaM with AMV isolates from potato, all presumed belonging to the "non-necrotic" strain and retrieved from GenBank, indicated that CaM shared >97.1% sequence identity with all but four Egyptian isolates. At the complete genome level, phylogenetic analysis of all available sequences demonstrated that RNA 1 and RNA 3 can be grouped into three major clades each, whereas RNA 2 can be clustered into two clades. CaM and Ca175-1, an AMV isolate that was deemed non-necrotic in a previous study, had different phylogenetic clade patterns, indicating different RNA 1-RNA 2-RNA 3 haplotypes: IA-I-IB (CaM) versus Ca175-1 (IB-II-IA). Despite the difference in haplotype composition, CaM and Ca175-1 induced similar levels of internal necrosis in tubers of Innovator and its parent Shepody. The results suggest that the internal necrosis in AMV-infected tubers depends on potato cultivar rather than on AMV strain/haplotype, and CaM is just a "regular" isolate of AMV.

Verticillium dahliae Disease Resistance and the Regulatory Pathway for Maturity and Tuberization in Potato.

Kleb. is a pathogenic fungus causing wilting, chlorosis, and early dying in potato ( L.). Genetic mapping of resistance to was done using a diploid population of potato. The major quantitative trait locus (QTL) for resistance was found on chromosome 5. The gene, controlling earliness of maturity and tuberization, was mapped within the interval. Another QTL on chromosome 9 co-localized with the wilt resistance gene marker. Epistasis analysis indicated that the loci on chromosomes 5 and 9 had a highly significant interaction, and that functioned downstream of The alleles were sequenced and found to encode StCDF1.1 and StCDF1.3. Interaction between the resistance allele and the was demonstrated, but not for Genome-wide expression QTL (eQTL) analysis was performed and genes with eQTL at the and loci were both found to have similar functions involving the chloroplast, including photosynthesis, which declines in both maturity and wilt. Among the gene ontology (GO) terms that were specific to genes with eQTL at the , but not the locus, were those associated with fungal defense. These results suggest that controls fungal defense and reduces early dying in wilt through affecting genetic pathway controlling tuberization timing.

High Resolution DNA Melting Assays for Detection of Rx1 and Rx2 for High-Throughput Marker-Assisted Selection for Extreme Resistance to Potato virus X in Tetraploid Potato.

Assessment of the existing PCR-gel electrophoresis-based methods for detection of Rx1 and Rx2, the genes that independently control extreme resistance (ER) to Potato virus X (PVX), indicated that the 5Rx1F/5Rx1R primer pair led to reliable detection of Rx1, whereas the 106Rx2F/106Rx2R primer pair detected Rx2 despite some nonspecific reactions in potato clones/cultivars without Rx2. However, the methodology is time consuming and does not differentiate the absence of Rx1/Rx2 from a failed PCR reaction. A newly designed primer pair that targets Rx1 and Rx2 as well as rx1 and rx2 produced an amplicon for all alleles. When the primer pair was combined with 5Rx1F/5Rx1R, respective amplicons were produced, although they were not distinguishable by regular agarose gel electrophoresis. When subjected to a high-resolution DNA melting (HRM) assay, two distinct melting profiles for Rx1 and rx1, respectively, were detected. Triplex PCR-gel electrophoresis and -HRM assay for detection of Rx1, Rx2, and rx1/rx2 were also performed. The efficacy of the HRM assays were validated in potato cultivars/clones with known phenotypes, indicating its potential for high-throughput selection of potato clones/cultivars carrying Rx1 or Rx2. Duplex PCR-HRM assays of over 600 progeny from 12 crosses involving various parents correctly detected the presence or absence of Rx1 in each progeny, allowing accurate prediction of the phenotype. Progeny that tested positive for Rx1 by HRM exhibited ER to PVX whereas progeny that tested negative for Rx1 were susceptible to PVX infection. The genotype of each parent and the possible presence of Nx in two Rx1-possessing parents are also discussed.

Development and Validation of High-Resolution Melting Markers Derived from Rysto STS Markers for High-Throughput Marker-Assisted Selection of Potato Carrying Rysto.

Sequence analysis of the chromosome region harboring the sequence-tagged site (STS) markers YES3-3A and YES3-3B for Rysto, a gene responsible for extreme resistance to Potato virus Y (PVY) in potato, was performed in tetraploid potato 'Barbara' (Rrrr) and 'AC Chaleur' (rrrr) as well as their progeny selections. Three and two sequence variants were identified in Barbara resistant (R) selections and AC Chaleur susceptible (S) selections, respectively. Further analysis indicates that the variant with a 21-nucleotide (nt) deletion is likely the chromosome copy harboring the STS markers. Two primer pairs, one targeting the region containing a 20-nt deletion and the other targeting the region anchoring the YES3-3A reverse primer, were designed. As anticipated, pair one produced two visible fragments in Barbara-R bulk and one visible fragment in AC Chaleur-S bulk; pair two produced one visible fragment in all samples. When subjected to high-resolution melting (HRM) analysis, two distinct melting profiles for R and S samples were observed. Analysis of 147 progeny of Barbara × AC Chaleur revealed 72 and 75 progeny with R and S melting profiles, respectively, which was consistent with YES3-3A and YES3-3B assays and phenotyping analysis, thus demonstrating the potential of HRM profiles as novel molecular markers for Rysto. The efficacy of the newly developed HRM markers for high-throughput marker-assisted selection for Rysto-conferred resistance to PVY was validated further with three populations involving Barbara as the R parent.

Differences between the Bud End and Stem End of Potatoes in Dry Matter Content, Starch Granule Size, and Carbohydrate Metabolic Gene Expression at the Growing and Sprouting Stages.

Potatoes usually have the tuber bud end dominance in growth during tuber bulking and in tuber sprouting, likely using carbohydrates from the tuber stem end. We hypothesized that the tuber bud end and tuber stem end coordination in carbohydrate metabolism gene expression is different between the bulking dominance and sprouting dominance of the tuber bud end. After comparing the growing tubers at harvest from a green vine and the stage that sprouts just started to emerge after storage of tubers at room temperature, we found the following: (1) Dry matter content was higher in the tuber stem end than the tuber bud end at both stages. (2) The starch granule size was larger in the tuber bud end than in the tuber stem end. (3) The tuber bud end had higher gene expression for starch synthesis but a lower gene expression of sucrose transporters than the tuber stem end during tuber growing. (4) The tuber stem end at the sprouting stage showed more active gene expression in both starch degradation and resynthesis, suggesting more active export of carbohydrates, than the tuber bud end. The results indicate that the starch accumulation mechanism in the tuber bud end was different between field growing and post-harvest sprouting tubers and that tubers already increased dry matter and average starch granule sizes in the tuber bud end prior to the rapid growth of sprouts.

Implications of miR166 and miR159 induction to the basal response mechanisms of an andigena potato (Solanum tuberosum subsp. andigena) to salinity stress, predicted from network models in Arabidopsis.

MicroRNA (miRNA) mediated changes in gene expression by post-transcriptional modulation of major regulatory transcription factors is a potent mechanism for integrating growth and stress-related responses. Exotic plants including many traditional varieties of Andean potatoes (Solanum tuberosum subsp. andigena) are known for better adaptation to marginal environments. Stress physiological studies confirmed earlier reports on the salinity tolerance potentials of certain andigena cultivars. Guided by the hypothesis that certain miRNAs play important roles in growth modulation under suboptimal conditions, we identified and characterized salinity stress-responsive miRNA-target gene pairs in the andigena cultivar Sullu by parallel analysis of noncoding and coding RNA transcriptomes. Inverse relationships were established by the reverse co-expression between two salinity stress-regulated miRNAs (miR166, miR159) and their target transcriptional regulators HD-ZIP-Phabulosa/Phavulota and Myb101, respectively. Based on heterologous models in Arabidopsis, the miR166-HD-ZIP-Phabulosa/Phavulota network appears to be involved in modulating growth perhaps by mediating vegetative dormancy, with linkages to defense-related pathways. The miR159-Myb101 network may be important for the modulation of vegetative growth while also controlling stress-induced premature transition to reproductive phase. We postulate that the induction of miR166 and miR159 under salinity stress represents important network hubs for balancing gene expression required for basal growth adjustments.

Comparative metabolite profiling of Solanum tuberosum against six wild Solanum species with Colorado potato beetle resistance.

The Colorado potato beetle Leptinotarsa decemlineata (Say) (CPB) is a coleopteran herbivore that feeds on the foliage on Solanum species, in particular, potato. Six resistant wild Solanum species were identified, and two of these species had low levels of glycoalkaloids. Comparative analysis of the untargeted metabolite profiles of the foliage using UPLC-qTOF-MS was done to find metabolites shared between the wild species but not with Solanum tuberosum (L.) to identify resistance-related metabolites. It was found that only S. tuberosum produced the triose glycoalkaloids solanine and chaconine. Instead, the six wild species produced glycoalkaloids that shared in common tetrose sugar side chains. Additionally, there were non-glycoalkaloid metabolites associated with resistance including hydroxycoumarin and a phenylpropanoid, which were produced in all wild species but not in S. tuberosum.

New evidence on the relationship between Microsporidia and Fungi: a genome-wide analysis by DarkHorse software.

Microsporidia are a group of obligate intracellular eukaryotic parasites that infect a wide variety of species, including humans. Phylogenetic analysis indicates a relationship between the Microsporidia and the Fungi. However, most results are based on the analysis of relatively few genes. DarkHorse analysis involves the transformation of BLAST results into a lineage probability index (LPI) value and allows for the comparison of genes for an entire genome with those of other genomes. Thus, we can see which genes from the microsporidia score most closely based on the LPI with other eukaryotic organisms. In this analysis, we calculated the LPI for each gene from the genomes of 7 Microsporidia, Antonospora locustae, Enterocytozoon bieneusi, Encephalitozoon cuniculi, Encephalitozoon intestinalis, Nosema bombycis, Nosema ceranae, and Nematocida parisii, to analyze the genetic relationships between Microsporidia and other species. It was found that many (91%) genes were most closely correlated with genes from other microsporidial genomes and had the highest mean LPI (0.985), indicating a monophyletic origin of the Microsporidia. In a subsequent analysis, we excluded the other Microsporidia from the analysis to look for relationships before the divergence of Microsporidia, and found that 43% of the microsporidial genes scored highest with fungal genes, and a higher mean LPI was found with Fungi than with other kingdoms, suggesting that Microsporidia is closely related to Fungi at the genomic level. Microsporidial genes were functionally clustered based on the KOG (Eukaryotic COG) database, and the possible lineages for each gene family were discussed in concert with the DarkHorse results.

Phenotypic stability and genome-wide association study of late blight resistance in potato genotypes adapted to the tropical highlands.

Potato genotypes from a breeding population adapted to tropical highlands were analyzed for the stability of late blight resistance and also for marker-phenotype association. We harmonized the historical evaluation data, consisting of observations spanning 6 years from two field sites utilizing a resistance scale constructed by comparing the area under the disease progress curve (AUDPC) values of 172 genotypes with that of susceptible control 'Yungay'. In total, 70 potato genotypes had a coefficient of variability <0.5 and were considered stable across the environments tested. A principal component analysis demonstrated that the ensemble of experiments formed two distinct groups that reflect the stability of genotype resistance to late blight. Phytophthora infestans isolates present in the experimental fields belonged to the EC-1 clonal lineage and showed variation in virulence beyond the concept of the avirulence determined by the conventionally used R1-R11 differential set. A single-nucleotide polymorphism (SNP) marker on chromosome 9 was associated with late blight resistance and linked to instability. Genotypes with either AACC or AAAC combinations for this SNP were highly resistant only in some environments, while the genotypes with the AAAA combination had more moderate levels of resistance but were stable across environments.

Decreased defense gene expression in tolerance versus resistance to Verticillium dahliae in potato.

Verticillium dahliae Kleb., a soil-borne fungus that colonizes vascular tissues, induces wilting, chlorosis and early senescence in potato. Difference in senescence timing found in two diploid potato clones, 07506-01 and 12120-03, was studied and genetic variation in response to V. dahliae infection was identified as a causal factor. The clone, 07506-01, was infected with V. dahliae but did not develop symptoms, indicating tolerance to the pathogen. The other diploid clone, 12120-03 had low levels of pathogen with infection and moderate symptoms indicating partial resistance. 07506-01 was found to carry two susceptible alleles of the Ve2 gene and 12120-03 carried one Ve2 resistant and one susceptible allele. Infected leaves of the two clones were compared using gene expression profiling with the Potato Oligonucleotide Chip Initiative (POCI) microrarray. The results provide further evidence for differences in response of the two clones to infection with V. dahliae. Chlorophyll biosynthesis was higher in the tolerant 07506-01 compared to partially resistant 12120-03. On the other hand, expression of fungal defense genes, Ve resistance genes and defense phytohormone biosynthetic enzyme genes was decreased in 07506-01 compared to 12120-03 suggesting defense responses were suppressed in tolerance compared to resistance. Transcription factor gene expression differences pointed to the WRKY family as potential regulators of V. dahliae responses in potato.

Physiological and molecular adaptations to drought in Andean potato genotypes.

The drought stress tolerance of two Solanum tuberosum subsp. andigena landraces, one hybrid (adgxtbr) and Atlantic (S. tuberosum subsp. tuberosum) has been evaluated. Photosynthesis in the Andigena landraces during prolonged drought was maintained significantly longer than in the Tuberosum (Atlantic) line. Among the Andigena landraces, 'Sullu' (SUL) was more drought resistant than 'Negra Ojosa' (NOJ). Microarray analysis and metabolite data from leaf samples taken at the point of maximum stress suggested higher mitochondrial metabolic activity in SUL than in NOJ. A greater induction of chloroplast-localized antioxidant and chaperone genes in SUL compared with NOJ was evident. ABA-responsive TFs were more induced in NOJ compared with SUL, including WRKY1, mediating a response in SA signalling that may give rise to increased ROS. NOJ may be experiencing higher ROS levels than SUL. Metabolite profiles of NOJ were characterized by compounds indicative of stress, for example, proline, trehalose, and GABA, which accumulated to a higher degree than in SUL. The differences between the Andigena lines were not explained by protective roles of compatible solutes; hexoses and complex sugars were similar in both landraces. Instead, lower levels of ROS accumulation, greater mitochondrial activity and active chloroplast defences contributed to a lower stress load in SUL than in NOJ during drought.

Genes driving potato tuber initiation and growth: identification based on transcriptional changes using the POCI array.

The increasing amount of available expressed gene sequence data makes whole-transcriptome analysis of certain crop species possible. Potato currently has the second largest number of publicly available expressed sequence tag (EST) sequences among the Solanaceae. Most of these ESTs, plus other proprietary sequences, were combined and used to generate a unigene assembly. The set of 246,182 sequences produced 46,345 unigenes, which were used to design a 44K 60-mer oligo array (Potato Oligo Chip Initiative: POCI). In this study, we attempt to identify genes controlling and driving the process of tuber initiation and growth by implementing large-scale transcriptional changes using the newly developed POCI array. Major gene expression profiles could be identified exhibiting differential expression at key developmental stages. These profiles were associated with functional roles in cell division and growth. A subset of genes involved in the regulation of the cell cycle, based on their Gene Ontology classification, exhibit a clear transient upregulation at tuber onset indicating increased cell division during these stages. The POCI array allows the study of potato gene expression on a much broader level than previously possible and will greatly enhance analysis of transcriptional control mechanisms in a wide range of potato research areas. POCI sequence and annotation data are publicly available through the POCI database ( http://pgrc.ipk-gatersleben.de/poci ).