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.

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

In planta processing and glycosylation of a nematode CLAVATA3/ENDOSPERM SURROUNDING REGION-like effector and its interaction with a host CLAVATA2-like receptor to promote parasitism.

Like other biotrophic plant pathogens, plant-parasitic nematodes secrete effector proteins into host cells to facilitate infection. Effector proteins that mimic plant CLAVATA3/ENDOSPERM SURROUNDING REGION-related (CLE) proteins have been identified in several cyst nematodes, including the potato cyst nematode (PCN); however, the mechanistic details of this cross-kingdom mimicry are poorly understood. Plant CLEs are posttranslationally modified and proteolytically processed to function as bioactive ligands critical to various aspects of plant development. Using ectopic expression coupled with nanoliquid chromatography-tandem mass spectrometry analysis, we show that the in planta mature form of proGrCLE1, a multidomain CLE effector secreted by PCN during infection, is a 12-amino acid arabinosylated glycopeptide (named GrCLE1-1Hyp4,7g) with striking structural similarity to mature plant CLE peptides. This glycopeptide is more resistant to hydrolytic degradation and binds with higher affinity to a CLAVATA2-like receptor (StCLV2) from potato (Solanum tuberosum) than its nonglycosylated forms. We further show that StCLV2 is highly up-regulated at nematode infection sites and that transgenic potatoes with reduced StCLV2 expression are less susceptible to PCN infection, indicating that interference of the CLV2-mediated signaling pathway confers nematode resistance in crop plants. These results strongly suggest that phytonematodes have evolved to utilize host cellular posttranslational modification and processing machinery for the activation of CLE effectors following secretion into plant cells and highlight the significance of arabinosylation in regulating nematode CLE effector activity. Our finding also provides evidence that multidomain CLEs are modified and processed similarly to single-domain CLEs, adding new insight into CLE maturation in plants.

Applying network and genetic analysis to the potato metabolome.

Compositional traits in potato [Solanum tuberosum L.] are economically important but genetically complex, often controlled by many loci of small effect; new methods need to be developed to accelerate analysis and improvement of such traits, like chip quality. In this study, we used network analysis to organize hundreds of metabolic features detected by mass spectrometry into groups, as a precursor to genetic analysis. 981 features were condensed into 44 modules; module eigenvalues were used for genetic mapping and correlation analysis with phenotype data collected by the Solanaceae Coordinated Agricultural Project. Half of the modules were associated with at least one SNP according to GWAS; 11 of those modules were also significantly correlated with chip color. Within those modules features associated with chipping provide potential targets for selection in addition to selection for reduced glucose. Loci associated with module eigenvalues were not evenly distributed throughout the genome but were instead clustered on chromosomes 3, 7, and 8. Comparison of GWAS on single features and modules of clustered features often identified the same SNPs. However, features with related chemistries (for example, glycoalkaloids with precursor/product relationships) were not found to be near neighbors in the network analysis and did not share common SNPs from GWAS. Instead, the features within modules were often structurally disparate, suggesting that linkage disequilibrium complicates network analyses in potato. This result is consistent with recent genomic studies of potato showing that chromosomal rearrangements that create barriers to recombination are common in cultivated germplasm.

Phased, chromosome-scale genome assemblies of tetraploid potato reveal a complex genome, transcriptome, and predicted proteome landscape underpinning genetic diversity.

Cultivated potato is a clonally propagated autotetraploid species with a highly heterogeneous genome. Phased assemblies of six cultivars including two chromosome-scale phased genome assemblies revealed extensive allelic diversity, including altered coding and transcript sequences, preferential allele expression, and structural variation that collectively result in a highly complex transcriptome and predicted proteome, which are distributed across the homologous chromosomes. Wild species contribute to the extensive allelic diversity in tetraploid cultivars, demonstrating ancestral introgressions predating modern breeding efforts. As a clonally propagated autotetraploid that undergoes limited meiosis, dysfunctional and deleterious alleles are not purged in tetraploid potato. Nearly a quarter of the loci bore mutations are predicted to have a high negative impact on protein function, complicating breeder's efforts to reduce genetic load. The StCDF1 locus controls maturity, and analysis of six tetraploid genomes revealed that 12 allelic variants of StCDF1 are correlated with maturity in a dosage-dependent manner. Knowledge of the complexity of the tetraploid potato genome with its rampant structural variation and embedded deleterious and dysfunctional alleles will be key not only to implementing precision breeding of tetraploid cultivars but also to the construction of homozygous, diploid potato germplasm containing favorable alleles to capitalize on heterosis in F1 hybrids.

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 analysis of pigmented tuber flesh in potato.

Interest in anthocyanin-pigmented potato tuber flesh is increasing. To genetically map and characterize loci that influence this trait, diploid potato clone 10618-01, which has partially pigmented flesh, was crossed with diploid 320-02, which has white flesh. Almost all progeny exhibited purple coloration in the flesh, with some clones having only a small percentage of tissue pigmented, other clones having most tissue pigmented, and the majority of clones showing intermediate color phenotypes. The two parents and 228 progeny were genotyped with 493 AFLP, 8 CAPS, and 13 SSR markers. QTLs influencing extent of flesh pigmentation were detected on chromosomes 5, 8, and 9. The potato homolog of Petunia an1, a basic helix-loop-helix (bHLH) transcriptional regulator of anthocyanin biosynthesis, was found to co-localize with the QTL on chromosome 9. A CAPS marker based on this gene was used to evaluate a collection of 21 tetraploid potato clones with highly or fully pigmented red or purple flesh, as well as 53 cultivars with white or yellow flesh. All 21 pigmented-flesh clones shared a marker allele that was present in only 21 of the 53 white and yellow clones, suggesting that a common bHLH allele contributes toward, although it is clearly not sufficient for, highly or fully pigmented tuber flesh in cultivated potato.

The potato developer (D) locus encodes an R2R3 MYB transcription factor that regulates expression of multiple anthocyanin structural genes in tuber skin.

A dominant allele at the D locus (also known as I in diploid potato) is required for the synthesis of red and purple anthocyanin pigments in tuber skin. It has previously been reported that D maps to a region of chromosome 10 that harbors one or more homologs of Petunia an2, an R2R3 MYB transcription factor that coordinately regulates the expression of multiple anthocyanin biosynthetic genes in the floral limb. To test whether D acts similarly in tuber skin, RT-PCR was used to evaluate the expression of flavanone 3-hydroxylase (f3h), dihydroflavonol 4-reductase (dfr) and flavonoid 3',5'-hydroxylase (f3'5'h). All three genes were expressed in the periderm of red- and purple-skinned clones, while dfr and f3'5'h were not expressed, and f3h was only weakly expressed, in white-skinned clones. A potato cDNA clone with similarity to an2 was isolated from an expression library prepared from red tuber skin, and an assay developed to distinguish the two alleles of this gene in a diploid potato clone known to be heterozygous Dd. One allele was observed to cosegregate with pigmented skin in an F(1) population of 136 individuals. This allele was expressed in tuber skin of red- and purple-colored progeny, but not in white tubers, while other parental alleles were not expressed in white or colored tubers. The allele was placed under the control of a doubled 35S promoter and transformed into the light red-colored cultivar Désirée, the white-skinned cultivar Bintje, and two white diploid clones known to lack the functional allele of D. Transformants accumulated pigment in tuber skin, as well as in other tissues, including young foliage, flower petals, and tuber flesh.

Genetic Variance Partitioning and Genome-Wide Prediction with Allele Dosage Information in Autotetraploid Potato.

As one of the world's most important food crops, the potato (Solanum tuberosum L.) has spurred innovation in autotetraploid genetics, including in the use of SNP arrays to determine allele dosage at thousands of markers. By combining genotype and pedigree information with phenotype data for economically important traits, the objectives of this study were to (1) partition the genetic variance into additive vs. nonadditive components, and (2) determine the accuracy of genome-wide prediction. Between 2012 and 2017, a training population of 571 clones was evaluated for total yield, specific gravity, and chip fry color. Genomic covariance matrices for additive (G), digenic dominant (D), and additive × additive epistatic (G#G) effects were calculated using 3895 markers, and the numerator relationship matrix (A) was calculated from a 13-generation pedigree. Based on model fit and prediction accuracy, mixed model analysis with G was superior to A for yield and fry color but not specific gravity. The amount of additive genetic variance captured by markers was 20% of the total genetic variance for specific gravity, compared to 45% for yield and fry color. Within the training population, including nonadditive effects improved accuracy and/or bias for all three traits when predicting total genotypic value. When six F1 populations were used for validation, prediction accuracy ranged from 0.06 to 0.63 and was consistently lower (0.13 on average) without allele dosage information. We conclude that genome-wide prediction is feasible in potato and that it will improve selection for breeding value given the substantial amount of nonadditive genetic variance in elite germplasm.

Software for Genome-Wide Association Studies in Autopolyploids and Its Application to Potato.

Genome-wide association studies (GWAS) are widely used in diploid species to study complex traits in diversity and breeding populations, but GWAS software tailored to autopolyploids is lacking. The objectives of this research were to (i) develop an R package for autopolyploids based on the + mixed model, (ii) validate the software with simulated data, and (iii) analyze a diversity panel of tetraploid potatoes. A unique feature of the R package, called GWASpoly, is its ability to model different types of polyploid gene action, including additive, simplex dominant, and duplex dominant. Using a simulated tetraploid population, we confirmed our hypothesis that statistical power is higher when the assumed gene action in the GWAS model matches the gene action at unobserved quantitative trait loci (QTL). Thirteen traits were analyzed in the Solanaceae Coordinated Agricultural Project (SolCAP) potato diversity panel and, consistent with previous studies, significant QTL for tuber shape and eye depth co-localized on chromosome 10. For the other traits, only marginally significant QTL were detected, most likely due to insufficient statistical power: for simulated traits with a heritability () of 0.3, the median genome-wide power was only 0.01. Our results indicate that both marker density and population size were limiting factors for GWAS with the SolCAP panel.

Application of DNA markers linked to the potato H1 gene conferring resistance to pathotype Ro1 of Globodera rostochiensis.

Ninety-one potato genotypes (cultivars and breeding lines) selected as resistant or susceptible to pathotype Ro1 of Globodera rostochiensis were screened for the presence of two PCR markers, 0.14 and 0.76 kb in length. Both PCR markers were linked with the H1 gene, located at the distal end of the long arm of chromosome V, and were present in 88 to 100% of the resistant cultivars and breeding lines. The 0.76 kb PCR marker was detected in all resistant genotypes and in approximately 86% of susceptible breeding lines as well as in all susceptible cultivars. The 0.14 kb marker was detected in 88% of resistant breeding lines and in 94% of resistant cultivars. Most of the susceptible genotypes tested (91% of cultivars, but only 50% of breeding lines) did not show the presence of the 0.14 kb marker. We conclude that the 0.14 kb H1 marker is likely to be useful for the proper selection of potato genotypes resistant to the Ro1 pathotype of G. rostochiensis.

Inheritance and genetic mapping of tuber eye depth in cultivated diploid potatoes.

Tuber eye depth of the potato (Solanum tuberosum L.) is an important trait for the processing quality and appearance of potatoes. In the present study, we used a cultivated diploid potato family (12601) of 107 plants to dissect the mode of inheritance and to map the gene(s) controlling the trait. The family segregated for both eye depth (deep vs shallow) and tuber shape (round vs long) traits. The deep eye (Eyd) phenotype was found to be associated with round tubers (Ro) in most progeny clones. Further evaluation of this population with molecular markers including simple sequence repeats, amplified fragment length polymorphism, and sequence-characterized amplified regions revealed that the primary locus for eye depth is located on chromosome 10. This map location was confirmed by evaluating a second diploid family (12586). The results of this study led to the following conclusions: (1) there is a major locus controlling the eye depth trait; (2) deep eye (Eyd) is dominant to shallow (eyd); (3) the Eyd/eyd locus is located on chromosome 10; and (4) the Eyd/eyd locus is closely linked with the major locus for tuber shape (Ro/ro), at a distance of about 4 cM.

The potato P locus codes for flavonoid 3',5'-hydroxylase.

The potato P locus is required for the production of blue/purple anthocyanin pigments in any tissue of the potato plant such as tubers, flowers, or stems. We have previously reported, based on RFLP mapping in tomato, that the gene coding for the anthocyanin biosynthetic enzyme flavonoid 3',5'-hydroxylase (f3'5'h) maps to the same region of the tomato genome as P maps in potato. To further evaluate this association a Petunia f3'5'h gene was used to screen a potato cDNA library prepared from purple-colored flowers and stems. Six positively hybridizing cDNA clones were sequenced and all appeared to be derived from a single gene that shares 85% sequence identity at the amino acid level with Petunia f3'5'h. The potato gene cosegregated with purple tuber color in a diploid F1 sub-population of 37 purple and 25 red individuals and was found to be expressed in tuber skin only in the presence of the anthocyanin regulatory locus I. A potato f3'5'h cDNA clone was placed under the control of a doubled CaMV 35S promoter and introduced into the red-skinned cultivar 'Desiree'. Tuber and stem tissues that are colored red in Desiree were purple in nine of 17 independently transformed lines.