Genomic Selection for Late Blight and Common Scab Resistance in Tetraploid Potato (Solanum tuberosum)

Genomic prediction for potato (Solanum tuberosum) quality traits improved through image analysis

Potato (Solanum tuberosum L.) is the most widely grown vegetable in the world. Consumers and processors evaluate potatoes based on quality traits such as shape and skin color, making these traits important targets for breeders. Achieving and evaluating genetic gain is facilitated by precise and accurate trait measures. Historically, quality traits have been measured using visual rating scales, which are subject to human error and necessarily lump individuals with distinct characteristics into categories. Image analysis offers a method of generating quantitative measures of quality traits. In this study, we use TubAR, an image-analysis R package, to generate quantitative measures of shape and skin color traits for use in genomic prediction. We developed and compared different genomic models based on additive and additive plus non-additive relationship matrices for two aspects of skin color, redness, and lightness, and two aspects of shape, roundness, and length-to-width ratio for fresh market red and yellow potatoes grown in Minnesota between 2020 and 2022. Similarly, we used the much larger chipping potato population grown during the same time to develop a multi-trait selection index including roundness, specific gravity, and yield. Traits ranged in heritability with shape traits falling between 0.23 and 0.85, and color traits falling between 0.34 and 0.91. Genetic effects were primarily additive with color traits showing the strongest effect (0.47), while shape traits varied based on market class. Modeling non-additive effects did not significantly improve prediction models for quality traits. The combination of image analysis and genomic prediction presents a promising avenue for improving potato quality traits.

Potato soup: analysis of cultivated potato gene bank populations reveals high diversity and little structure

Cultivated potatoes are incredibly diverse, ranging from diploid to pentaploid and encompass four different species. They are adapted to disparate environments and conditions and carry unique alleles for resistance to pests and pathogens. Describing how diversity is partitioned within and among these populations is essential to understanding the potato genome and effectively utilizing landraces in breeding. This task is complicated by the difficulty of making comparisons across cytotypes and extensive admixture within section petota. We genotyped 730 accessions from the US Potato genebank including wild diploids and cultivated diploids and tetraploids using Genotype-by-sequencing. This data set allowed us to interrogate population structure and diversity as well as generate core subsets which will support breeders in efficiently screening genebank material for biotic and abiotic stress resistance alleles. We found that even controlling for ploidy, tetraploid material exhibited higher observed and expected heterozygosity than diploid accessions. In particular group chilotanum material was the most heterozygous and the only taxa not to exhibit any inbreeding. This may in part be because group chilotanum has a history of introgression not just from wild species, but landraces as well. All group chilotanum, exhibits introgression from group andigenum except clones from Southern South America near its origin, where the two groups are not highly differentiated. Moving north, we do not observe evidence for the same level of admixture back into group andigenum. This suggests that extensive history of admixture is a particular characteristic of chilotanum.

Genotyping-by-sequencing targets genic regions and improves resolution of genome-wide association studies in autotetraploid potato

KEY MESSAGE

De novo genotyping in potato using methylation-sensitive GBS discovers SNPs largely confined to genic or gene-associated regions and displays enhanced effectiveness in estimating LD decay rates, population structure and detecting GWAS associations over 'fixed' SNP genotyping platform. Study also reports the genetic architectures including robust sequence-tagged marker-trait associations for sixteen important potato traits potentially carrying higher transferability across a wider range of germplasm. This study deploys recent advancements in polyploid analytical approaches to perform complex trait analyses in cultivated tetraploid potato. The study employs a 'fixed' SNP Infinium array platform and a 'flexible and open' genome complexity reduction-based sequencing method (GBS, genotyping-by-sequencing) to perform genome-wide association studies (GWAS) for several key potato traits including the assessment of population structure and linkage disequilibrium (LD) in the studied population. GBS SNPs discovered here were largely confined (~ 90%) to genic or gene-associated regions of the genome demonstrating the utility of using a methylation-sensitive restriction enzyme (PstI) for library construction. As compared to Infinium array SNPs, GBS SNPs displayed enhanced effectiveness in estimating LD decay rates and discriminating population subgroups. GWAS using a combined set of 30,363 SNPs identified 189 unique QTL marker-trait associations (QTL-MTAs) covering all studied traits. The majority of the QTL-MTAs were from GBS SNPs potentially illustrating the effectiveness of marker-dense de novo genotyping platforms in overcoming ascertainment bias and providing a more accurate correction for different levels of relatedness in GWAS models. GWAS also detected QTL 'hotspots' for several traits at previously known as well as newly identified genomic locations. Due to the current study exploiting genome-wide genotyping and de novo SNP discovery simultaneously on a large tetraploid panel representing a greater diversity of the cultivated potato gene pool, the reported sequence-tagged MTAs are likely to have higher transferability across a wider range of potato germplasm and increased utility for expediting genomics-assisted breeding for the several complex traits studied.

Multienvironment genomic prediction in tetraploid potato

Multienvironment genomic prediction was applied to tetraploid potato using 147 potato varieties, tested for 2 years, in 3 locations representative of 3 distinct regions in Europe. Different prediction scenarios were investigated to help breeders predict genotypic performance in the regions from one year to the next, for genotypes that were tested this year (scenario 1), as well as new genotypes (scenario 3). In scenario 2, we predicted new genotypes for any one of the 6 trials, using all the information that is available. The choice of prediction model required assessment of the variance-covariance matrix in a mixed model that takes into account heterogeneity of genetic variances and correlations. This was done for each analyzed trait (tuber weight, tuber length, and dry matter) where examples of both limited and higher degrees of heterogeneity was observed. This explains why dry matter did not need complex multienvironment modeling to combine environments and increase prediction ability, while prediction in tuber weight, improved only when models were flexible enough to capture the heterogeneous variances and covariances between environments. We also found that the prediction abilities in a target trial condition decreased, if trials with a low genetic correlation to the target were included when training the model. Genomic prediction in tetraploid potato can work once there is clarity about the prediction scenario, a suitable training set is created, and a multienvironment prediction model is chosen based on the patterns of G×E indicated by the genetic variances and covariances.

Genomic selection in plant breeding: Key factors shaping two decades of progress

Genomic selection, the application of genomic prediction (GP) models to select candidate individuals, has significantly advanced in the past two decades, effectively accelerating genetic gains in plant breeding. This article provides a holistic overview of key factors that have influenced GP in plant breeding during this period. We delved into the pivotal roles of training population size and genetic diversity, and their relationship with the breeding population, in determining GP accuracy. Special emphasis was placed on optimizing training population size. We explored its benefits and the associated diminishing returns beyond an optimum size. This was done while considering the balance between resource allocation and maximizing prediction accuracy through current optimization algorithms. The density and distribution of single-nucleotide polymorphisms, level of linkage disequilibrium, genetic complexity, trait heritability, statistical machine-learning methods, and non-additive effects are the other vital factors. Using wheat, maize, and potato as examples, we summarize the effect of these factors on the accuracy of GP for various traits. The search for high accuracy in GP-theoretically reaching one when using the Pearson's correlation as a metric-is an active research area as yet far from optimal for various traits. We hypothesize that with ultra-high sizes of genotypic and phenotypic datasets, effective training population optimization methods and support from other omics approaches (transcriptomics, metabolomics and proteomics) coupled with deep-learning algorithms could overcome the boundaries of current limitations to achieve the highest possible prediction accuracy, making genomic selection an effective tool in plant breeding.

Gene expression analysis of potato drought-responsive genes under drought stress in potato (Solanum tuberosum L.) cultivars

The potato (Solanum tuberosum L.), an important field crop consumed extensively worldwide, is adversely affected by abiotic stress factors especially drought. Therefore, it is vital to understand the genetic mechanism under drought stress to decrease loose of yield and quality . This trial aimed to screen drought-responsive gene expressions of potato and determine the drought-tolerant potato cultivar. The trial pattern is a completely randomized block design (CRBD) with four replications under greenhouse conditions. Four cultivars (Brooke, Orwell, Vr808, Shc909) were irrigated with four different water regimes (control and three stress conditions), and the gene expression levels of 10 potato genes were investigated. The stress treatments as follows: Control = 100% field capacity; slight drought = 75% field capacity; moderate drought = 50% field capacity, and severe drought 25% field capacity. To understand the gene expression under drought stress in potato genotypes, RT-qPCR analysis was performed and results showed that the genes most associated with drought tolerance were the StRD22 gene, MYB domain transcription factor, StERD7, Sucrose Synthase (SuSy), ABC Transporter, and StDHN1. The StHSP100 gene had the lowest genetic expression in all cultivars. Among the cultivars, the Orwell exhibited the highest expression of the StRD22 gene under drought stress. Overall, the cultivar with the highest gene expression was the Vr808, closely followed by the Brooke cultivar. As a result, it was determined that potato cultivars Orwell, Vr808, and Brooke could be used as parents in breeding programs to develop drought tolerant potato cultivars.

Variation and genetic basis of mineral content in potato tubers and prospects for genomic selection

Malnutrition is a major public health concern in many parts of the world. Among other nutrients, minerals are necessary in the human diet. Potato tubers are a good source of minerals; they contribute 18% of the recommended dietary allowance of potassium; 6% of copper, phosphorus, and magnesium; and 2% of calcium and zinc. Increased public interest in improving the nutritional value of foods has prompted the evaluation of mineral content in tubers of advanced genotypes from the Texas A&M Potato Breeding Program and the investigation of the genetics underlying mineral composition in tubers. The objectives of this study were to i) assess phenotypic variation for mineral content in tubers of advanced potato genotypes, ii) identify genomic regions associated with tuber mineral content, and iii) obtain genomic-estimated breeding values. A panel of 214 advanced potato genotypes and reference varieties was phenotyped in three field environments in Texas for the content of 12 minerals in tubers and genotyped using the Infinium Illumina 22K V3 single nucleotide polymorphism (SNP) Array. There was significant variation between potato genotypes for all minerals evaluated except iron. As a market group, red-skinned potatoes had the highest amount of minerals, whereas russets had the lowest mineral content. Reds had significantly higher P, K, S, and Zn than russets and significantly higher P and Mg than chippers. Russets had significantly higher Ca, Mg, and Na than chippers. However, the chippers had significantly higher K than the russets. A genome-wide association study for mineral content using GWASpoly identified three quantitative trait loci (QTL) associated with potassium and manganese content on chromosome 5 and two QTL associated with zinc content on chromosome 7. The loci identified will contribute to a better understanding of the genetic basis of mineral content in potatoes. Genomic-estimated breeding values for mineral macro and micronutrients in tubers obtained with StageWise will guide the selection of parents and the advancement of genotypes in the breeding program to increase mineral content in potato tubers.

Identification of QTL associated with plant vine characteristics and infection response to late blight, early blight, and Verticillium wilt in a tetraploid potato population derived from late blight-resistant Palisade Russet

Potato late blight (causal agent Phytophthora infestans) is a disease of potatoes with economic importance worldwide. Control is primarily through field monitoring and the application of fungicides. Control of late blight with fungicides and host plant resistance is difficult, with documented cases of such control measures failing with the advent of new pathotypes of P. infestans. To better understand host plant resistance and to develop more durable late blight resistance, Quantitative Trait Locus/Loci (QTL) analysis was conducted on a tetraploid mapping population derived from late blight-resistant potato cultivar Palisade Russet. Additionally, QTL analyses for other traits such as Verticillium wilt and early blight resistance, vine size and maturity were performed to identify a potential relationship between multiple traits and prepare genetic resources for molecular markers useful in breeding programs. For this, one hundred ninety progenies from intercrossing Palisade Russet with a late blight susceptible breeding clone (ND028673B-2Russ) were assessed. Two parents and progenies were evaluated over a two-year period for response to infection by the US-8 genotype of P. infestans in inoculated field screenings in Corvallis, Oregon. In Aberdeen, Idaho, the same mapping population was also evaluated for phenotypic response to early blight and Verticillium wilt, and vine size and maturity in a field over a two-year period. After conducting QTL analyses with those collected phenotype data, it was observed that chromosome 5 has a significant QTL for all five traits. Verticillium wilt and vine maturity QTL were also observed on chromosome 1, and vine size QTL was also found on chromosomes 3 and 10. An early blight QTL was also detected on chromosome 2. The QTL identified in this study have the potential for converting into breeder-friendly molecular markers for marker-assisted selection.

Genome-wide association mapping and genomic prediction for late blight and potato cyst nematode resistance in potato (Solanum tuberosum L.)

Potatoes are an important source of food for millions of people worldwide. Biotic stresses, notably late blight and potato cyst nematodes (PCN) pose a major threat to potato production worldwide, and knowledge of genes controlling these traits is limited. A genome-wide association mapping study was conducted to identify the genomic regulators controlling these biotic stresses, and the genomic prediction accuracy was worked out using the GBLUP model of genomic selection (GS) in a panel of 222 diverse potato accessions. The phenotype data on resistance to late blight and two PCN species (Globodera pallida and G. rostochiensis) were recorded for three and two consecutive years, respectively. The potato panel was genotyped using genotyping by sequencing (GBS), and 1,20,622 SNP markers were identified. A total of 7 SNP associations for late blight resistance, 9 and 11 for G. pallida and G. rostochiensis, respectively, were detected by additive and simplex dominance models of GWAS. The associated SNPs were distributed across the chromosomes, but most of the associations were found on chromosomes 5, 10 and 11, which have been earlier reported as the hotspots of disease-resistance genes. The GS prediction accuracy estimates were low to moderate for resistance to G. pallida (0.04-0.14) and G. rostochiensis (0.14-0.21), while late blight resistance showed a high prediction accuracy of 0.42-0.51. This study provides information on the complex genetic nature of these biotic stress traits in potatoes and putative SNP markers for resistance breeding.

Breeding strategies for late blight resistance in potato crop: recent developments

Late blight (LB) is a serious disease that affects potato crop and is caused by Phytophthora infestans. Fungicides are commonly used to manage this disease, but this practice has led to the development of resistant strains and it also poses serious environmental and health risks. Therefore, breeding for resistance development can be the most effective strategies to control late blight. Various Solanum species have been utilized as a source of resistance genes to combat late blight disease. Several potential resistance genes and quantitative resistance loci (QRLs) have been identified and mapped through the application of molecular techniques. Furthermore, molecular markers closely linked to resistance genes or QRLs have been utilized to hasten the breeding process. However, the use of single-gene resistance can lead to the breakdown of resistance within a short period. To address this, breeding programs are now being focused on development of durable and broad-spectrum resistant cultivars by combining multiple resistant genes and QRLs using advanced molecular breeding tools such as marker-assisted selection (MAS) and cis-genic approaches. In addition to the strategies mentioned earlier, somatic hybridization has been utilized for the development and characterization of interspecific somatic hybrids. To further broaden the scope of late blight resistance breeding, approaches such as genomic selection, RNAi silencing, and various genome editing techniques can be employed. This study provides an overview of recent advances in various breeding strategies and their applications in improving the late blight resistance breeding program.

Efficient Genomic Prediction of Yield and Dry Matter in Hybrid Potato

There is an ongoing endeavor within the potato breeding sector to rapidly adapt potato from a clonal polyploid crop to a diploid hybrid potato crop. While hybrid breeding allows for the efficient generation and selection of parental lines, it also increases breeding program complexity and results in longer breeding cycles. Over the past two decades, genomic prediction has revolutionized hybrid crop breeding through shorter breeding cycles, lower phenotyping costs, and better population improvement, resulting in increased genetic gains for genetically complex traits. In order to accelerate the genetic gains in hybrid potato, the proper implementation of genomic prediction is a crucial milestone in the rapid improvement of this crop. The authors of this paper set out to test genomic prediction in hybrid potato using current genotyped material with two alternative models: one model that predicts the general combining ability effects (GCA) and another which predicts both the general and specific combining ability effects (GCA+SCA). Using a training set comprising 769 hybrids and 456 genotyped parental lines, we found that reasonable a prediction accuracy could be achieved for most phenotypes with both zero common parents (ρ=0.36-0.61) and one (ρ=0.50-0.68) common parent between the training and test sets. There was no benefit with the inclusion of non-additive genetic effects in the GCA+SCA model despite SCA variance contributing between 9% and 19% of the total genetic variance. Genotype-by-environment interactions, while present, did not appear to affect the prediction accuracy, though prediction errors did vary across the trial's targets. These results suggest that genomically estimated breeding values on parental lines are sufficient for hybrid yield prediction.

Optimal implementation of genomic selection in clone breeding programs-Exemplified in potato: I. Effect of selection strategy, implementation stage, and selection intensity on short-term genetic gain

Genomic selection (GS) is used in many animal and plant breeding programs to enhance genetic gain for complex traits. However, its optimal integration in clone breeding programs, for example potato, that up to now relied on phenotypic selection (PS) requires further research. In this study, we performed computer simulations based on an empirical genomic dataset of tetraploid potato to (i) investigate under a fixed budget how the weight of GS relative to PS, the stage of implementing GS, the correlation between an auxiliary trait and the target trait, the variance components, and the prediction accuracy affect the genetic gain of the target trait, (ii) determine the optimal allocation of resources maximizing the genetic gain of the target trait, and (iii) make recommendations to breeders how to implement GS in clone and especially potato breeding programs. In our simulation results, any selection strategy involving GS had a higher short-term genetic gain for the target trait than Standard-PS. In addition, we showed that implementing GS in consecutive selection stages can largely enhance short-term genetic gain and recommend the breeders to implement GS at single hills and A clone stages. Furthermore, we observed for selection strategies involving GS that the optimal allocation of resources maximizing the genetic gain of the target trait differed considerably from those typically used in potato breeding programs and, thus, require the adjustment of the selection and phenotyping intensities. The trends are described in our study. Therefore, our study provides new insight for breeders regarding how to optimally implement GS in a commercial potato breeding program to improve the short-term genetic gain for their target trait.

Resequencing worldwide spinach germplasm for identification of field resistance QTLs to downy mildew and assessment of genomic selection methods

Downy mildew, commercially the most important disease of spinach, is caused by the obligate oomycete Peronospora effusa. In the past two decades, new pathogen races have repeatedly overcome the resistance used in newly released cultivars, urging the need for more durable resistance. Commercial spinach cultivars are bred with major R genes to impart resistance to downy mildew pathogens and are effective against some pathogen races/isolates. This work aimed to evaluate the worldwide USDA spinach germplasm collections and commercial cultivars for resistance to downy mildew pathogen in the field condition under natural inoculum pressure and conduct genome wide association analysis (GWAS) to identify resistance-associated genomic regions (alleles). Another objective was to evaluate the prediction accuracy (PA) using several genomic prediction (GP) methods to assess the potential implementation of genomic selection (GS) to improve spinach breeding for resistance to downy mildew pathogen. More than four hundred diverse spinach genotypes comprising USDA germplasm accessions and commercial cultivars were evaluated for resistance to downy mildew pathogen between 2017-2019 in Salinas Valley, California and Yuma, Arizona. GWAS was performed using single nucleotide polymorphism (SNP) markers identified via whole genome resequencing (WGR) in GAPIT and TASSEL programs; detected 14, 12, 5, and 10 significantly associated SNP markers with the resistance from four tested environments, respectively; and the QTL alleles were detected at the previously reported region of chromosome 3 in three of the four experiments. In parallel, PA was assessed using six GP models and seven unique marker datasets for field resistance to downy mildew pathogen across four tested environments. The results suggest the suitability of GS to improve field resistance to downy mildew pathogen. The QTL, SNP markers, and PA estimates provide new information in spinach breeding to select resistant plants and breeding lines through marker-assisted selection (MAS) and GS, eventually helping to accumulate beneficial alleles for durable disease resistance.

Towards smart and sustainable development of modern berry cultivars in Europe

Fresh berries are a popular and important component of the human diet. The demand for high-quality berries and sustainable production methods is increasing globally, challenging breeders to develop modern berry cultivars that fulfill all desired characteristics. Since 1994, research projects have characterized genetic resources, developed modern tools for high-throughput screening, and published data in publicly available repositories. However, the key findings of different disciplines are rarely linked together, and only a limited range of traits and genotypes has been investigated. The Horizon2020 project BreedingValue will address these challenges by studying a broader panel of strawberry, raspberry and blueberry genotypes in detail, in order to recover the lost genetic diversity that has limited the aroma and flavor intensity of recent cultivars. We will combine metabolic analysis with sensory panel tests and surveys to identify the key components of taste, flavor and aroma in berries across Europe, leading to a high-resolution map of quality requirements for future berry cultivars. Traits linked to berry yields and the effect of environmental stress will be investigated using modern image analysis methods and modeling. We will also use genetic analysis to determine the genetic basis of complex traits for the development and optimization of modern breeding technologies, such as molecular marker arrays, genomic selection and genome-wide association studies. Finally, the results, raw data and metadata will be made publicly available on the open platform Germinate in order to meet FAIR data principles and provide the basis for sustainable research in the future.

The Past, Present, and Future of Host Plant Resistance in Cotton: An Australian Perspective

Cotton is a key global fiber crop. However, yield potential is limited by the presence of endemic and introduced pests and diseases. The introduction of host plant resistance (HPR), defined as the purposeful use of resistant crop cultivars to reduce the impact of pests and diseases, has been a key breeding target for the Commonwealth Scientific and Industrial Research Organisation (CSIRO) cotton breeding program. The program has seen success in releasing cultivars resistant to Bacterial blight, Verticillium wilt, Fusarium wilt, and Cotton bunchy top. However, emerging biotic threats such as Black root rot and secondary pests, are becoming more frequent in Australian cotton production systems. The uptake of tools and breeding methods, such as genomic selection, high throughput phenomics, gene editing, and landscape genomics, paired with the continued utilization of sources of resistance from Gossypium germplasm, will be critical for the future of cotton breeding. This review celebrates the success of HPR breeding activities in the CSIRO cotton breeding program and maps a pathway for the future in developing resistant cultivars.

Genetic parameters and selection in full-sib families of tall fescue using best linear unbiased prediction (BLUP) analysis

BACKGROUND

Better understanding of genetic structure of economic traits is crucial for identification and selection of superior genotypes in specific breeding programs. Best linear unbiased prediction (BLUP) is the most efficient method in this regard, which is poorly used in forage plant breeding. The present study aimed to assess genetic variation, estimate genetic parameters, and predict breeding values of five essential traits in full sib families (recognized by EST-SSR markers) of tall fescue using REML/BLUP procedure.

METHOD

Forty-two full-sib families of tall fescue (included of 120 individual genotypes), recognized by EST-SSR markers along with twenty-one their corresponding parental genotypes were assessed for biomass production and agro-morphological traits at three harvests (spring, summer, and autumn) in the field during 4 years (2017-2020).

RESULTS

Considerable genotypic variability was observed for all traits. Low narrow-sense heritability (h2n) for dry forage yield (DFY) at three harvest indicates that non-additive gene actions may play an important role in the inheritance of this trait. Higher h2n of yield related traits and flowering time and also significant genetic correlation of these traits with forage yield, suggests that selection based on these traits may lead to indirect genetic improvement of DFY.

CONCLUSION

Our results showed the adequacy of REML/BLUP procedure for identification and selection of preferable parental genotypes and progenies with higher breeding values for future breeding programs such as variety development in tall fescue. Parental genotypes 21 M, 1 M, and 20 L were identified as superior and stable genotypes and could also produce the best hybrid combinations when they were mostly used as maternal parent.

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-Based Genotype × Environment Prediction Enhances Potato (Solanum tuberosum L.) Improvement Using Pseudo-Diploid and Polysomic Tetraploid Modeling

Potato breeding must improve its efficiency by increasing the reliability of selection as well as identifying a promising germplasm for crossing. This study shows the prediction accuracy of genomic-estimated breeding values for several potato (Solanum tuberosum L.) breeding clones and the released cultivars that were evaluated at three locations in northern and southern Sweden for various traits. Three dosages of marker alleles [pseudo-diploid (A), additive tetrasomic polyploidy (B), and additive-non-additive tetrasomic polyploidy (C)] were considered in the genome-based prediction models, for single environments and multiple environments (accounting for the genotype-by-environment interaction or G × E), and for comparing two kernels, the conventional linear, Genomic Best Linear Unbiased Prediction (GBLUP) (GB), and the non-linear Gaussian kernel (GK), when used with the single-kernel genetic matrices of A, B, C, or when employing two-kernel genetic matrices in the model using the kernels from B and C for a single environment (models 1 and 2, respectively), and for multi-environments (models 3 and 4, respectively). Concerning the single site analyses, the trait with the highest prediction accuracy for all sites under A, B, C for model 1, model 2, and for GB and GK methods was tuber starch percentage. Another trait with relatively high prediction accuracy was the total tuber weight. Results show an increase in prediction accuracy of model 2 over model 1. Non-linear Gaussian kernel (GK) did not show any clear advantage over the linear kernel GBLUP (GB). Results from the multi-environments had prediction accuracy estimates (models 3 and 4) higher than those obtained from the single-environment analyses. Model 4 with GB was the best method in combination with the marker structure B for predicting most of the tuber traits. Most of the traits gave relatively high prediction accuracy under this combination of marker structure (A, B, C, and B-C), and methods GB and GK combined with the multi-environment with G × E model.

Germplasm, Breeding, and Genomics in Potato Improvement of Biotic and Abiotic Stresses Tolerance

Potato is one of the most important food crops in the world. Late blight, viruses, soil and tuber-borne diseases, insect-pests mainly aphids, whiteflies, and potato tuber moths are the major biotic stresses affecting potato production. Potato is an irrigated and highly fertilizer-responsive crop, and therefore, heat, drought, and nutrient stresses are the key abiotic stresses. The genus Solanum is a reservoir of genetic diversity, however, a little fraction of total diversity has been utilized in potato breeding. The conventional breeding has contributed significantly to the development of potato varieties. In recent years, a tremendous progress has been achieved in the sequencing technologies from short-reads to long-reads sequence data, genomes of Solanum species (i.e., pan-genomics), bioinformatics and multi-omics platforms such as genomics, transcriptomics, proteomics, metabolomics, ionomics, and phenomics. As such, genome editing has been extensively explored as a next-generation breeding tool. With the available high-throughput genotyping facilities and tetraploid allele calling softwares, genomic selection would be a reality in potato in the near future. This mini-review covers an update on germplasm, breeding, and genomics in potato improvement for biotic and abiotic stress tolerance.

The power of genomic estimated breeding values for selection when using a finite population size in genetic improvement of tetraploid potato

Potato breeding relies heavily on visual phenotypic scoring for clonal selection. Obtaining robust phenotypic data can be labor intensive and expensive, especially in the early cycles of a potato breeding program where the number of genotypes is very large. We have investigated the power of genomic estimated breeding values (GEBVs) for selection from a limited population size in potato breeding. We collected genotypic data from 669 tetraploid potato clones from all cycles of a potato breeding program, as well as phenotypic data for eight important breeding traits. The genotypes were partitioned into a training and a test population distinguished by cycle of selection in the breeding program. GEBVs for seven traits were predicted for individuals from the first stage of the breeding program (T1) which had not undergone any selection, or individuals selected at least once in the field (T2). An additional approach in which GEBVs were predicted within and across full-sib families from unselected material (T1) was tested for four breeding traits. GEBVs were obtained by using a Bayesian Ridge Regression model estimating single marker effects and phenotypic data from individuals at later stages of selection of the breeding program. Our results suggest that, for most traits included in this study, information from individuals from later stages of selection cannot be utilized to make selections based on GEBVs in earlier clonal generations. Predictions of GEBVs across full-sib families yielded similarly low prediction accuracies as across generations. The most promising approach for selection using GEBVs was found to be making predictions within full-sib families.

Genome-Wide Association Mapping of Late Blight Tolerance Trait in Potato (Solanum tuberosum L.)

Uncovering the genetic basis and optimizing the late blight tolerance trait in potatoes (Solanum tuberosum L.) are crucial for potato breeding. Late blight disease is one of the most significant diseases hindering potato production. The traits of late blight tolerance were evaluated for 284 potato cultivars to identify loci significantly associated with the late blight tolerance trait. Of all, 37 and 15 were the most tolerant to disease, and 107 and 30 were the most susceptible. A total of 22,489 high-quality single-nucleotide polymorphisms and indels were identified in 284 potato cultivars. All the potato cultivars were clustered into eight subgroups using population structure analysis and principal component analysis, which were consistent with the results of the phylogenetic tree analysis. The average genetic diversity for all 284 potato cultivars was 0.216, and the differentiation index of each subgroup was 0.025–0.149. Genome-wide linkage disequilibrium (LD) analysis demonstrated that the average LD was about 0.9 kb. A genome-wide association study using a mixed linear model identified 964 loci significantly associated with the late blight tolerance trait. Fourteen candidate genes for late blight tolerance traits were identified, including genes encoding late blight tolerance protein, chitinase 1, cytosolic nucleotide-binding site–leucine-rich repeat tolerance protein, protein kinase, ethylene-responsive transcription factor, and other potential plant tolerance-related proteins. This study provides novel insights into the genetic architecture of late blight tolerance traits and will be helpful for late blight tolerance in potato breeding.

Quantitative Trait Locus Mapping for Common Scab Resistance in a Tetraploid Potato Full-Sib Population

Despite the negative impact of common scab (Streptomyces spp.) on the potato industry, little is known about the genetic architecture of resistance to this bacterial disease in the crop. We evaluated a mapping population (∼150 full sibs) derived from a cross between two tetraploid potatoes ('Atlantic' × B1829-5) in three environments (MN11, PA11, ME12) under natural common scab pressure. Three measures to common scab reaction, namely percentage of scabby tubers and disease area and lesion indices, were found to be highly correlated (>0.76). Because of the large environmental effect, heritability values were zero for all three traits in MN11, but moderate to high in PA11 and ME12 (∼0.44 to 0.79). We identified a single quantitative trait locus (QTL) for lesion index in PA11, ME12, and joint analyses on linkage group 3, explaining ∼22 to 30% of the total variation. The identification of QTL haplotypes and candidate genes contributing to disease resistance can support genomics-assisted breeding approaches in the crop.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Understanding the Effectiveness of Genomic Prediction in Tetraploid Potato

Use of genomic prediction (GP) in tetraploid is becoming more common. Therefore, we think it is the right time for a comparison of GP models for tetraploid potato. GP models were compared that contrasted shrinkage with variable selection, parametric vs. non-parametric models and different ways of accounting for non-additive genetic effects. As a complement to GP, association studies were carried out in an attempt to understand the differences in prediction accuracy. We compared our GP models on a data set consisting of 147 cultivars, representing worldwide diversity, with over 39 k GBS markers and measurements on four tuber traits collected in six trials at three locations during 2 years. GP accuracies ranged from 0.32 for tuber count to 0.77 for dry matter content. For all traits, differences between GP models that utilised shrinkage penalties and those that performed variable selection were negligible. This was surprising for dry matter, as only a few additive markers explained over 50% of phenotypic variation. Accuracy for tuber count increased from 0.35 to 0.41, when dominance was included in the model. This result is supported by Genome Wide Association Study (GWAS) that found additive and dominance effects accounted for 37% of phenotypic variation, while significant additive effects alone accounted for 14%. For tuber weight, the Reproducing Kernel Hilbert Space (RKHS) model gave a larger improvement in prediction accuracy than explicitly modelling epistatic effects. This is an indication that capturing the between locus epistatic effects of tuber weight can be done more effectively using the semi-parametric RKHS model. Our results show good opportunities for GP in 4x potato.

Global multi-environment resistance QTL for foliar late blight resistance in tetraploid potato with tropical adaptation

The identification of environmentally stable and globally predictable resistance to potato late blight is challenged by the clonal and polyploid nature of the crop and the rapid evolution of the pathogen. A diversity panel of tetraploid potato germplasm bred for multiple resistance and quality traits was genotyped by genotyping by sequencing (GBS) and evaluated for late blight resistance in three countries where the International Potato Center (CIP) has established breeding work. Health-indexed, in vitro plants of 380 clones and varieties were distributed from CIP headquarters and tuber seed was produced centrally in Peru, China, and Ethiopia. Phenotypes were recorded following field exposure to local isolates of Phytophthora infestans. QTL explaining resistance in four experiments conducted across the three countries were identified in chromosome IX, and environment-specific QTL were found in chromosomes III, V, and X. Different genetic models were evaluated for prediction ability to identify best performing germplasm in each and all environments. The best prediction ability (0.868) was identified with the genomic best linear unbiased predictors (GBLUPs) when using the diploid marker data and QTL-linked markers as fixed effects. Genotypes with high levels of resistance in all environments were identified from the B3, LBHT, and B3-LTVR populations. The results show that many of the advanced clones bred in Peru for high levels of late blight resistance maintain their resistance in Ethiopia and China, suggesting that the centralized selection strategy has been largely successful.

Genomic Selection in an Outcrossing Autotetraploid Fruit Crop: Lessons From Blueberry Breeding

Blueberry (Vaccinium corymbosum and hybrids) is a specialty crop with expanding production and consumption worldwide. The blueberry breeding program at the University of Florida (UF) has greatly contributed to expanding production areas by developing low-chilling cultivars better adapted to subtropical and Mediterranean climates of the globe. The breeding program has historically focused on recurrent phenotypic selection. As an autopolyploid, outcrossing, perennial, long juvenile phase crop, blueberry breeding cycles are costly and time consuming, which results in low genetic gains per unit of time. Motivated by applying molecular markers for a more accurate selection in the early stages of breeding, we performed pioneering genomic selection studies and optimization for its implementation in the blueberry breeding program. We have also addressed some complexities of sequence-based genotyping and model parametrization for an autopolyploid crop, providing empirical contributions that can be extended to other polyploid species. We herein revisited some of our previous genomic selection studies and showed for the first time its application in an independent validation set. In this paper, our contribution is three-fold: (i) summarize previous results on the relevance of model parametrizations, such as diploid or polyploid methods, and inclusion of dominance effects; (ii) assess the importance of sequence depth of coverage and genotype dosage calling steps; (iii) demonstrate the real impact of genomic selection on leveraging breeding decisions by using an independent validation set. Altogether, we propose a strategy for using genomic selection in blueberry, with the potential to be applied to other polyploid species of a similar background.

MultiGWAS: An integrative tool for Genome Wide Association Studies in tetraploid organisms

The genome-wide association studies (GWASs) are essential to determine the genetic bases of either ecological or economic phenotypic variation across individuals within populations of the model and nonmodel organisms. For this research question, the GWAS replication testing different parameters and models to validate the results' reproducibility is common. However, straightforward methodologies that manage both replication and tetraploid data are still missing. To solve this problem, we designed the MultiGWAS, a tool that does GWAS for diploid and tetraploid organisms by executing in parallel four software packages, two designed for polyploid data (GWASpoly and SHEsis) and two designed for diploid data (GAPIT and TASSEL). MultiGWAS has several advantages. It runs either in the command line or in a graphical interface; it manages different genotype formats, including VCF. Moreover, it allows control for population structure, relatedness, and several quality control checks on genotype data. Besides, MultiGWAS can test for additive and dominant gene action models, and, through a proprietary scoring function, select the best model to report its associations. Finally, it generates several reports that facilitate identifying false associations from both the significant and the best-ranked association Single Nucleotide Polymorphisms (SNPs) among the four software packages. We tested MultiGWAS with public tetraploid potato data for tuber shape and several simulated data under both additive and dominant models. These tests demonstrated that MultiGWAS is better at detecting reliable associations than using each of the four software packages individually. Moreover, the parallel analysis of polyploid and diploid software that only offers MultiGWAS demonstrates its utility in understanding the best genetic model behind the SNP association in tetraploid organisms. Therefore, MultiGWAS probed to be an excellent alternative for wrapping GWAS replication in diploid and tetraploid organisms in a single analysis environment.

The recombination landscape and multiple QTL mapping in a Solanum tuberosum cv. 'Atlantic'-derived F1 population

There are many challenges involved with the genetic analyses of autopolyploid species, such as the tetraploid potato, Solanum tuberosum (2n = 4x = 48). The development of new analytical methods has made it valuable to re-analyze an F₁ population (n = 156) derived from a cross involving 'Atlantic', a widely grown chipping variety in the USA. A fully integrated genetic map with 4285 single nucleotide polymorphisms, spanning 1630 cM, was constructed with MAPpoly software. We observed that bivalent configurations were the most abundant ones (51.0~72.4% depending on parent and linkage group), though multivalent configurations were also observed (2.2~39.2%). Seven traits were evaluated over four years (2006-8 and 2014) and quantitative trait loci (QTL) mapping was carried out using QTLpoly software. Based on a multiple-QTL model approach, we detected 21 QTL for 15 out of 27 trait-year combination phenotypes. A hotspot on linkage group 5 was identified with co-located QTL for maturity, plant yield, specific gravity, and internal heat necrosis resistance evaluated over different years. Additional QTL for specific gravity and dry matter were detected with maturity-corrected phenotypes. Among the genes around QTL peaks, we found those on chromosome 5 that have been previously implicated in maturity (StCDF1) and tuber formation (POTH1). These analyses have the potential to provide insights into the biology and breeding of tetraploid potato and other autopolyploid species.

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.

A Bioinformatics Pipeline to Identify a Subset of SNPs for Genomics-Assisted Potato Breeding

Modern potato breeding methods following a genomic-led approach provide means for shortening breeding cycles and increasing breeding efficiency across selection cycles. Acquiring genetic data for large breeding populations remains expensive. We present a pipeline to reduce the number of single nucleotide polymorphisms (SNPs) to lower the cost of genotyping. First, we reduced the number of individuals to be genotyped with a high-throughput method according to the multi-trait variation as defined by principal component analysis of phenotypic characteristics. Next, we reduced the number of SNPs by pruning for linkage disequilibrium. By adjusting the square of the correlation coefficient between two adjacent loci, we obtained reduced subsets of SNPs. We subsequently tested these SNP subsets by two methods; (1) a genome-wide association study (GWAS) for marker identification, and (2) genomic selection (GS) to predict genomic estimated breeding values. The results indicate that both GWAS and GS can be done without loss of information after SNP reduction. The pipeline allows for creating custom SNP subsets to cover all variation found in any particular breeding population. Low-throughput genotyping will reduce the genotyping cost associated with large populations, thereby making genomic breeding methods applicable to large potato breeding populations by reducing genotyping costs.

Prediction based on estimated breeding values using genealogy for tuber yield and late blight resistance in auto-tetraploid potato (Solanum tuberosum L.)

Estimated breeding values using best linear unbiased prediction (BLUP) through pedigree relationship can enhance selection efficiency and save time as well as resources in autotetraploid potato breeding program. Here, we used historical preliminary yield evaluation trials data of 469–619 breeding lines for tuber yield and late blight resistance to estimate heritability and BLUP based breeding values modelling auto-tetraploid inheritance in mixed model analysis. The pedigree file had a depth of 3–4 generations with total 370 individuals including 111 founders. Heritability estimates varied from 0.15 for marketable tuber yield to 0.47 for late blight resistance computed using A matrix. The prediction accuracy for total tuber yield, marketable tuber yield and late blight resistance (AUDPC) was 0.53 ± 0.02, 0.44 ± 0.02 and 0.81 ± 0.01, respectively. The prediction accuracy was highest for late blight resistance and moderate for total and marketable tuber yield. The prediction bias measured as regression of observed phenotype values on predicted values for late blight resistance was almost nil in comparison to total and marketable tuber yield. Moderate to high prediction accuracies for tuber yields and late blight resistance suggest the selection of genotypes based on EBVs in Indian potato breeding programme for higher genetic gain.

Identification and functional characterisation of late blight resistance polymorphic genes in Russet Burbank potato cultivar

In plants, the biosynthesis of the phenylpropanoid, flavonoid and fatty acid pathway monomers, polymers and conjugated metabolites play a vital role in disease resistance. These are generally deposited to reinforce cell walls to contain the pathogen to the site of infection. Identification of sequence variants in genes that biosynthesise these resistance metabolites can explain the mechanisms of disease resistance. The resistant and susceptible genotypes inoculated with Phytophthora infestans were RNA sequenced to identify the single nucleotide polymorphisms (SNPs) and insertion/deletion (InDel) variations. The SNPs/InDels were annotated and classified into different categories based on their effect on gene functions. In the selected 25 biosynthetic genes overlapping 39 transcripts, a total of 52 SNPs/InDels were identified in the protein-coding (CDS) regions. These were categorised as deleterious based on prediction of their effects on protein structure and function. The SNPs/InDels data obtained in this study can be used in genome editing to enhance late blight resistance in Russet Burbank and other potato cultivars.

Making the most of all data: Combining non-genotyped and genotyped potato individuals with HBLUP

Using genomic information to predict phenotypes can improve the accuracy of estimated breeding values and can potentially increase genetic gain over conventional breeding. In this study, we investigated the prediction accuracies achieved by best linear unbiased prediction (BLUP) for nine potato phenotypic traits using three types of relationship matrices pedigree ABLUP, genomic GBLUP, and a hybrid matrix (H) combining pedigree and genomic information (HBLUP). Deep pedigree information was available for >3000 different potato breeding clones evaluated over four years. Genomic relationships were estimated from >180,000 informative SNPs generated using a genotyping-by-sequencing transcriptome (GBS-t) protocol for 168 cultivars, many of which were parents of clones. Two validation scenarios were implemented, namely "Genotyped Cultivars Validation" (a subset of genotyped lines as validation set) and "Non-genotyped 2009 Progenies Validation". Most of the traits showed moderate to high narrow sense heritabilities (range 0.22-0.72). In the Genotyped Cultivars Validation, HBLUP outperformed ABLUP on prediction accuracies for all traits except early blight, and outperformed GBLUP for most of the traits except tuber shape, tuber eye depth and boil after-cooking darkening. This is evidence that the in-depth relationship within the H matrix could potentially result in better prediction accuracy in comparison to using A or G matrix individually. The prediction accuracies of the Non-genotyped 2009 Progenies Validation were comparable between ABLUP and HBLUP, varying from 0.17-0.70 and 0.18-0.69, respectively. Better prediction accuracy and less bias in prediction using HBLUP is of practical utility to breeders as all breeding material is ranked on the same scale leading to improved selection decisions. In addition, our approach provides an economical alternative to utilize historic breeding data with current genotyped individuals in implementing genomic selection.

Investigation of Streptomyces scabies Causing Potato Scab by Various Detection Techniques, Its Pathogenicity and Determination of Host-Disease Resistance in Potato Germplasm

Streptomyces scabies is a Gram-positive bacterial pathogen that causes common scab disease to several crops, particularly in the potato. It is a soil borne pathogen, a very devastating scab pathogen and difficult to manage in the field. Streptomyces has several species that cause common scab such as S. scabiei, S. acidiscabies, S. europaeiscabiei, S. luridiscabiei, S. niveiscabiei, S. puniciscabiei, S. reticuliscabiei, S. stelliscabiei, S. turgidiscabies, S. ipomoeae. Common scab disease harmfully affects potato economic and market value due to the presence of black spots on the tuber. Owing to its genetic diversity and pathogenicity, the determination of pathogen presence in potato fields is still challenging. In this study, S. scabies genetic diversity was measured by surveying five potato-growing areas of Pakistan during the growing season 2019. A total of 50 Streptomyces isolates, including S. scabies, S. acidiscabies, S. griseoflavus were isolated and identified based on morphologic, biochemical and molecular analysis. Virulent confirmation assays confirmed ten virulent strains of Streptomyces spp. On the potato cultivars Cardinal and Santee. Among the Streptomyces species, S. scabies showed the highest scab index, followed by S. acidiscabies and S. griseoflavus by exhibiting the scab-like lesions on potato tubers. Ten potato cultivars were screened against these virulent isolates of Streptomyces. The Faisalabad white variety showed the highest scab index followed By Cardinal, Tourag, Kuroda, Santee, Lady Rosetta, Asterix, Diamant, Faisalabad red and Sadaf. Moreover, genetic diversity and pathogenicity of Streptomyces spp. on potato tubers were also likely diverse in different geographical regions and also potato cultivars. This study represents a contribution to understanding the local interaction between potatoes and Streptomyces spp. in Pakistan. It will aid in supporting a solution for the management of this pathogen around the world.

Comparative transcriptome expression analysis in susceptible and resistant potato (Solanum tuberosum) cultivars to common scab (Streptomyces scabies) revealed immune priming responses in the incompatible interaction

Common scab disease in potato has become a widespread issue in major potato production areas, leading to increasing economic losses. Varietal resistance is seen as a viable and long-term scab management strategy. However, the genes and mechanisms of varietal resistance are unknown. In the current study, a comparative RNA transcriptome sequencing and differential gene signaling and priming sensitization studies were conducted in two potato cultivars that differ by their response to common scab (Streptomyces scabies), for unraveling the genes and pathways potentially involved in resistance within this pathosystem. We report on a consistent and contrasted gene expression pattern from 1,064 annotated genes differentiating a resistant (Hindenburg) and a susceptible (Green Mountain) cultivars, and identified a set of 273 co-regulated differentially expressed genes in 34 pathways that more likely reflect the genetic differences of the cultivars and metabolic mechanisms involved in the scab pathogenesis and resistance. The data suggest that comparative transcriptomic phenotyping can be used to predict scab lesion phenotype in breeding lines using mature potato tuber. The study also showed that the resistant cultivar, Hindenburg, has developed and maintained a capacity to sense and prime itself for persistent response to scab disease over time, and suggests an immune priming reaction as a mechanism for induced-resistance in scab resistant potato cultivars. The set of genes identified, described, and discussed in the study paves the foundation for detailed characterizations towards tailoring and designing procedures for targeted gene knockout through gene editing and phenotypic evaluation.

SeqBreed: a python tool to evaluate genomic prediction in complex scenarios

Background

Genomic prediction (GP) is a method whereby DNA polymorphism information is used to predict breeding values for complex traits. Although GP can significantly enhance predictive accuracy, it can be expensive and difficult to implement. To help design optimum breeding programs and experiments, including genome-wide association studies and genomic selection experiments, we have developed SeqBreed, a generic and flexible forward simulator programmed in python3.

Results

SeqBreed accommodates sex and mitochondrion chromosomes as well as autopolyploidy. It can simulate any number of complex phenotypes that are determined by any number of causal loci. SeqBreed implements several GP methods, including genomic best linear unbiased prediction (GBLUP), single-step GBLUP, pedigree-based BLUP, and mass selection. We illustrate its functionality with Drosophila genome reference panel (DGRP) sequence data and with tetraploid potato genotype data.

Conclusions

SeqBreed is a flexible and easy to use tool that can be used to optimize GP or genome-wide association studies. It incorporates some of the most popular GP methods and includes several visualization tools. Code is open and can be freely modified. Software, documentation, and examples are available at https://github.com/miguelperezenciso/SeqBreed.

Exploring Deep Learning for Complex Trait Genomic Prediction in Polyploid Outcrossing Species

Genomic prediction (GP) is the procedure whereby the genetic merits of untested candidates are predicted using genome wide marker information. Although numerous examples of GP exist in plants and animals, applications to polyploid organisms are still scarce, partly due to limited genome resources and the complexity of this system. Deep learning (DL) techniques comprise a heterogeneous collection of machine learning algorithms that have excelled at many prediction tasks. A potential advantage of DL for GP over standard linear model methods is that DL can potentially take into account all genetic interactions, including dominance and epistasis, which are expected to be of special relevance in most polyploids. In this study, we evaluated the predictive accuracy of linear and DL techniques in two important small fruits or berries: strawberry and blueberry. The two datasets contained a total of 1,358 allopolyploid strawberry (2n=8x=112) and 1,802 autopolyploid blueberry (2n=4x=48) individuals, genotyped for 9,908 and 73,045 single nucleotide polymorphism (SNP) markers, respectively, and phenotyped for five agronomic traits each. DL depends on numerous parameters that influence performance and optimizing hyperparameter values can be a critical step. Here we show that interactions between hyperparameter combinations should be expected and that the number of convolutional filters and regularization in the first layers can have an important effect on model performance. In terms of genomic prediction, we did not find an advantage of DL over linear model methods, except when the epistasis component was important. Linear Bayesian models were better than convolutional neural networks for the full additive architecture, whereas the opposite was observed under strong epistasis. However, by using a parameterization capable of taking into account these non-linear effects, Bayesian linear models can match or exceed the predictive accuracy of DL. A semiautomatic implementation of the DL pipeline is available at https://github.com/lauzingaretti/deepGP/.

Genomic Selection for Prediction of Fruit-Related Traits in Pepper (Capsicum spp.)

Pepper (Capsicum spp.) fruit-related traits are critical determinants of quality. These traits are controlled by quantitatively inherited genes for which marker-assisted selection (MAS) has proven insufficiently effective. Here, we evaluated the potential of genomic selection, in which genotype and phenotype data for a training population are used to predict phenotypes of a test population with only genotype data, for predicting fruit-related traits in pepper. We measured five fruit traits (fruit length, fruit shape, fruit width, fruit weight, and pericarp thickness) in 351 accessions from the pepper core collection, including 229 Capsicum annuum, 48 Capsicum baccatum, 48 Capsicum chinense, 25 Capsicum frutescens, and 1 Capsicum chacoense in 4 years at two different locations and genotyped these accessions using genotyping-by-sequencing. Among the whole core collection, considering its genetic distance and sexual incompatibility, we only included 302 C. annum complex (229 C. annuum, 48 C. chinense, and 25 C. frutescens) into further analysis. We used phenotypic and genotypic data to investigate genomic prediction models, marker density, and effects of population structure. Among 10 genomic prediction methods tested, Reproducing Kernel Hilbert Space (RKHS) produced the highest prediction accuracies (measured as correlation between predicted values and observed values) across the traits, with accuracies of 0.75, 0.73, 0.84, 0.83, and 0.82 for fruit length, fruit shape, fruit width, fruit weight, and pericarp thickness, respectively. Overall, prediction accuracies were positively correlated with the number of markers for fruit traits. We tested our genomic selection models in a separate population of recombinant inbred lines derived from two parental lines from the core collection. Despite the large difference in genetic diversity between the training population and the test population, we obtained moderate prediction accuracies of 0.32, 0.34, 0.50, and 0.48 for fruit length, fruit shape, fruit width, and fruit weight, respectively. This use of genomic selection for fruit-related traits demonstrates the potential use of core collections and genomic selection as tools for crop improvement.

pSBVB: A Versatile Simulation Tool To Evaluate Genomic Selection in Polyploid Species

Genomic Selection (GS) is the procedure whereby molecular information is used to predict complex phenotypes and it is standard in many animal and plant breeding schemes. However, only a small number of studies have been reported in horticultural crops, and in polyploid species in particular. In this paper, we have developed a versatile forward simulation tool, called polyploid Sequence Based Virtual Breeding (pSBVB), to evaluate GS strategies in polyploids; pSBVB is an efficient gene dropping software that can simulate any number of complex phenotypes, allowing a very flexible modeling of phenotypes suited to polyploids. As input, it takes genotype data from the founder population, which can vary from single nucleotide polymorphisms (SNP) chips up to sequence, a list of causal variants for every trait and their heritabilities, and the pedigree. Recombination rates between homeologous chromosomes can be specified, so that both allo- and autopolyploid species can be considered. The program outputs phenotype and genotype data for all individuals in the pedigree. Optionally, it can produce several genomic relationship matrices that consider exact or approximate genotype values. pSBVB can therefore be used to evaluate GS strategies in polyploid species (say varying SNP density, genetic architecture or population size, among other factors), or to optimize experimental designs for association studies. We illustrate pSBVB with SNP data from tetraploid potato and partial sequence data from octoploid strawberry, and we show that GS is a promising breeding strategy for polyploid species but that the actual advantage critically depends on the underlying genetic architecture. Source code, examples and a complete manual are freely available in GitHub https://github.com/lauzingaretti/pSBVB.

Evaluation of Genomic Prediction for Pasmo Resistance in Flax

Pasmo (Septoria linicola) is a fungal disease causing major losses in seed yield and quality and stem fibre quality in flax. Pasmo resistance (PR) is quantitative and has low heritability. To improve PR breeding efficiency, the accuracy of genomic prediction (GP) was evaluated using a diverse worldwide core collection of 370 accessions. Four marker sets, including three defined by 500, 134 and 67 previously identified quantitative trait loci (QTL) and one of 52,347 PR-correlated genome-wide single nucleotide polymorphisms, were used to build ridge regression best linear unbiased prediction (RR-BLUP) models using pasmo severity (PS) data collected from field experiments performed during five consecutive years. With five-fold random cross-validation, GP accuracy as high as 0.92 was obtained from the models using the 500 QTL when the average PS was used as the training dataset. GP accuracy increased with training population size, reaching values >0.9 with training population size greater than 185. Linear regression of the observed PS with the number of positive-effect QTL in accessions provided an alternative GP approach with an accuracy of 0.86. The results demonstrate the GP models based on marker information from all identified QTL and the 5-year PS average is highly effective for PR prediction.