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.

Genome Reduction in Tetraploid Potato Reveals Genetic Load, Haplotype Variation, and Loci Associated With Agronomic Traits.

The cultivated potato (Solanum tuberosum) has a complex genetic structure due to its autotetraploidy and vegetative propagation which leads to accumulation of mutations and a highly heterozygous genome. A high degree of heterozygosity has been considered to be the main driver of fitness and agronomic trait performance in potato improvement efforts, which is negatively impacted by genetic load. To understand the genetic landscape of cultivated potato, we constructed a gynogenic dihaploid (2n = 2x = 24) population from cv. Superior, prior to development of a high-density genetic map containing 12,753 single nucleotide polymorphisms (SNPs). Common quantitative trait loci (QTL) were identified for tuber traits, vigor and height on chromosomes 2, 4, 7, and 10, while specific QTL for number of inflorescences per plant, and tuber shape were present on chromosomes 4, 6, 10, and 11. Simplex rather than duplex loci were mainly associated with traits. In general, the Q allele (main effect) detected in one or two homologous chromosomes was associated with lower mean trait values suggesting the importance of dosage allelic effects, and the presence of up to two undesired alleles in the QTL region. Loss of heterozygosity has been associated with a lower rate of fitness, yet no correlation between the percent heterozygosity and increased fitness or agronomic performance was observed. Based upon linkage phase, we reconstructed the four homologous chromosome haplotypes of cv. Superior. revealing heterogeneity throughout the genome yet nearly duplicate haplotypes occurring among the homologs of particular chromosomes. These results suggest that the potentially deleterious mutations associated with genetic load in tetraploid potato could be mitigated by multiple loci which is consistent with the theory that epistasis complicates the identification of associations between markers and phenotypic performance.

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.

Mapping Loci That Control Tuber and Foliar Symptoms Caused by PVY in Autotetraploid Potato (Solanum tuberosum L.).

Potato tuber necrotic ringspot disease (PTNRD) is a tuber deformity associated with infection by the tuber necrotic strain of Potato virus Y (PVYNTN). PTNRD negatively impacts tuber quality and marketability, and poses a serious threat to seed and commercial potato production worldwide. PVYNTN symptoms differ in the cultivars Waneta and Pike: Waneta expresses severe PTNRD and foliar mosaic with vein and leaf necrosis, whereas Pike does not express PTNRD and mosaic is the only foliar symptom. To map loci that influence tuber and foliar symptoms, 236 F1 progeny of a cross between Waneta and Pike were inoculated with PVYNTN isolate NY090029 and genotyped using 12,808 potato SNPs. Foliar symptom type and severity were monitored for 10 wk, while tubers were evaluated for PTNRD expression at harvest and again after 60 d in storage. Pairwise correlation analyses indicate a strong association between PTNRD and vein necrosis (τ = 0.4195). QTL analyses revealed major-effect QTL on chromosomes 4 and 5 for mosaic, 4 for PTNRD, and 5 for foliar necrosis symptoms. Locating QTL associated with PVY-related symptoms provides a foundation for breeders to develop markers that can be used to eliminate potato clones with undesirable phenotypes, e.g., those likely to develop PTNRD or to be symptomless carriers of PVY.

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.

Comparative Analysis of Regions with Distorted Segregation in Three Diploid Populations of Potato.

Genes associated with gametic and zygotic selection could underlie segregation distortion, observed as alterations of expected Mendelian genotypic frequencies in mapping populations. We studied highly dense genetic maps based on single nucleotide polymorphisms to elucidate the genetic nature of distorted segregation in potato. Three intra- and interspecific diploid segregating populations were used. DRH and D84 are crosses between the sequenced doubled monoploid DM 1-3 516 R44 Solanum tuberosum Group Phureja and either RH89-039-16 S. tuberosum or 84SD22, a S. tuberosum × S. chacoense hybrid. MSX902 is an interspecific cross between 84SD22 and Ber83 S. berthaultii × 2 × species mosaic. At the 0.05 significance level, 21%, 57%, and 51% of the total markers mapped in DRH, D84, and MSX902 exhibited distorted segregation, respectively. Segregation distortion regions for DRH were located on chromosomes 9 and 12; for D84 on chromosomes 2, 3, 4, 6, 7, and 8; and on chromosomes 1, 2, 7, 9, and 12 for MSX902. In general, each population had unique segregation distortion regions and directions of distortion. Interspecific crosses showed greater levels of distorted segregation and lower recombination rates as determined from the male parents. The different genomic regions where the segregation distortion regions occurred in the three populations likely reflect unique genetic combinations producing distorted segregation.

Genetic Linkage Mapping of Economically Important Traits in Cultivated Tetraploid Potato (Solanum tuberosum L.).

The objective of this study was to construct a single nucleotide polymorphism (SNP)-based genetic map at the cultivated tetraploid level to locate quantitative trait loci (QTL) contributing to economically important traits in potato (Solanum tuberosum L.). The 156 F1 progeny and parents of a cross (MSL603) between "Jacqueline Lee" and "MSG227-2" were genotyped using the Infinium 8303 Potato Array. Furthermore, the progeny and parents were evaluated for foliar late blight reaction to isolates of the US-8 genotype of Phytophthora infestans (Mont.) de Bary and vine maturity. Linkage analyses and QTL mapping were performed using a novel approach that incorporates allele dosage information. The resulting genetic maps contained 1972 SNP markers with an average density of 1.36 marker per cM. QTL mapping identified the major source of late blight resistance in "Jacqueline Lee." The best SNP marker mapped ~0.54 Mb from a resistance hotspot on the long arm of chromosome 9. For vine maturity, the major-effect QTL was located on chromosome 5 with allelic effects from both parents. A candidate SNP marker for this trait mapped ~0.25 Mb from the StCDF1 gene, which is a candidate gene for the maturity trait. The identification of markers for P. infestans resistance will enable the introgression of multiple sources of resistance through marker-assisted selection. Moreover, the discovery of a QTL for late blight resistance not linked to the QTL for vine maturity provides the opportunity to use marker-assisted selection for resistance independent of the selection for vine maturity classifications.

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.

Integration of two diploid potato linkage maps with the potato genome sequence.

To facilitate genome-guided breeding in potato, we developed an 8303 Single Nucleotide Polymorphism (SNP) marker array using potato genome and transcriptome resources. To validate the Infinium 8303 Potato Array, we developed linkage maps from two diploid populations (DRH and D84) and compared these maps with the assembled potato genome sequence. Both populations used the doubled monoploid reference genotype DM1-3 516 R44 as the female parent but had different heterozygous diploid male parents (RH89-039-16 and 84SD22). Over 4,400 markers were mapped (1,960 in DRH and 2,454 in D84, 787 in common) resulting in map sizes of 965 (DRH) and 792 (D84) cM, covering 87% (DRH) and 88% (D84) of genome sequence length. Of the mapped markers, 33.5% were in candidate genes selected for the array, 4.5% were markers from existing genetic maps, and 61% were selected based on distribution across the genome. Markers with distorted segregation ratios occurred in blocks in both linkage maps, accounting for 4% (DRH) and 9% (D84) of mapped markers. Markers with distorted segregation ratios were unique to each population with blocks on chromosomes 9 and 12 in DRH and 3, 4, 6 and 8 in D84. Chromosome assignment of markers based on linkage mapping differed from sequence alignment with the Potato Genome Sequencing Consortium (PGSC) pseudomolecules for 1% of the mapped markers with some disconcordant markers attributable to paralogs. In total, 126 (DRH) and 226 (D84) mapped markers were not anchored to the pseudomolecules and provide new scaffold anchoring data to improve the potato genome assembly. The high degree of concordance between the linkage maps and the pseudomolecules demonstrates both the quality of the potato genome sequence and the functionality of the Infinium 8303 Potato Array. The broad genome coverage of the Infinium 8303 Potato Array compared to other marker sets will enable numerous downstream applications.

Evaluation of natural and engineered resistance mechanisms in potato against Colorado potato beetle in a no-choice field study.

The Colorado potato beetle, Leptinotarsa decemlineata Say, is the major insect pest of potato, Solanum tuberosum L., in eastern North America and is renowned for resistance development, currently resistant to >40 insecticides worldwide. Host plant resistance may assist in delaying in resistance development to insecticides. We evaluated natural host plant resistance mechanisms (glandular trichomes and Solanum chacoense Bitter-derived resistance) and engineered resistance mechanisms (Bacillus thuringiensis [Bt] Berliner cry3A and cry1Ia1) in a no-choice cage study. Six different potato lines representing four host plant resistance mechanisms were evaluated over 2 yr. Egg masses were placed in each cage (one egg mass per plant). Almost no feeding was observed in the Bt-cry3A lines, and only minor feeding was observed in the Bt-cry1Ia1 lines in either year. On the S. chacoense-derived line, there was significantly less defoliation than on either the susceptible line or the glandular trichome line in 2003. In 2004, there was significantly higher defoliation on the S. chacoense-derived line than on the susceptible line or glandular trichome line. The defoliation of the Solanum chacoense-derived line was largely due to larvae clipping the petioles, rather than consumption of the leaves. Defoliation on the glandular trichome line did not differ significantly from the defoliation of the susceptible line, suggesting glandular trichomes may not be effective in controlling larvae and preventing defoliation. This study suggested that Bt can provide high levels of resistance, but the natural resistance mechanisms tested here are variable for control of Colorado potato beetle larvae in no-choice situations.