Genome-wide association study and genomic prediction of white rust resistance in USDA GRIN spinach germplasm

Genome and GWAS analyses for soybean cyst nematode resistance in USDA world-wide common bean (Phaseolus vulgaris) germplasm

Soybean cyst nematode (SCN), Heterodera glycines, has become a significant threat in common bean (Phaseolus vulgaris) production, particularly in regions like the upper Midwest USA. Host genetic resistance offers an effective and environmentally friendly approach to managing SCN. This study aimed to conduct a genome-wide association study (GWAS) and genomic prediction for resistance to SCN HG Types 7 (race 6), 2.5.7 (race 5), and 1.3.6.7 (race 14) using 0.7 million whole-genome resequencing-generated SNPs in 354 USDA worldwide common bean germplasm accessions. Among these, 26 lines exhibited resistance to all three HG types, with a female index (FI) of less than 10. Four QTL regions on chromosomes (Chr) 2, 3, 6, and 10 were associated with resistance to HG Type 7; four regions on Chrs 2, 6, 9, and 11 were associated with resistance to HG Type 2.5.7; and three regions on Chrs 2, 6, and 10 were associated with resistance to HG Type 1.3.6.7. Cross-prediction revealed high prediction ability (PA) of 75% (r-value) for resistance to each of the three HG types. However, low PA was observed for SCN resistance through across-population prediction between the two domestications, Mesoamerican and Andean common bean accessions. Yet, using a population of mixed Mesoamerican and Andean accessions as a training set showed a high PA to predict either sub-population. This study provides SNP markers for marker-assisted selection and high PA for genomic selection in common bean molecular breeding, enabling the selection of lines and plants with high SCN resistance. Moreover, the study observed high PA for resistance among the three HG types. Interestingly, the most highly associated SNP markers and QTL for SCN resistance varied between the two domestications, and SCN resistance is more associated with the Mesoamerican domestication than the Andean domestication. This result suggests that resistance to SCN in common bean may be related to domestication rather than co-evolution with SCN.

Integration of machine learning and genome-wide association study to explore the genomic prediction accuracy of agronomic trait in oats (Avena sativa L.)

Machine learning (ML) has garnered significant attention for its potential to enhance the accuracy of genomic predictions (GPs) in various economic crops with the use of complete genomic information. Genome-wide association studies (GWAS) are widely used to pinpoint trait-related causal variant loci in genomes. However, the simultaneous integration of both methods for crop genome prediction necessitates further research. In this study, we integrated ML and GWAS to assess the efficiency of GP for seven key agronomic traits in 195 oat (Avena sativa) cultivars from major oat-growing regions around the world. A total of 94 trait-associated single nucleotide polymorphisms were identified through the GWAS study. GP studies were conducted using the classical model genomic best linear unbiased prediction (GBLUP) and six ML models. GBLUP performed poorly in predicting all traits except flag leaf width, while none of the ML models consistently provided the best prediction accuracy across all traits. The prediction accuracy of the GWAS-derived markers was better than that of the use of genome-wide markers, and plant height had the highest prediction rate at 100 GWAS-derived markers, and the rest of the traits for which more markers were required. These results play an important role in advancing the use of GP in small oat breeding programs by optimizing the prediction rate of GP and reducing the number of markers, confirming that high prediction rates can be achieved with smaller datasets.

Genome-Wide Association Study and Genomic Prediction of Soybean Mosaic Virus Resistance

Soybean mosaic virus (SMV), a pathogen responsible for inducing leaf mosaic or necrosis symptoms, significantly compromises soybean seed yield and quality. According to the classification system in the United States, SMV is categorized into seven distinct strains (G1 to G7). In this study, we performed a genome-wide association study (GWAS) in GAPIT3 using four analytical models (MLM, MLMM, FarmCPU, and BLINK) on 218 soybean accessions. We identified 22 SNPs significantly associated with G1 resistance across chromosomes 1, 2, 3, 12, 13, 17, and 18. Notably, a major quantitative trait locus (QTL) spanning 873 kb (29.85-30.73 Mb) on chromosome 13 exhibited strong association with SMV G1 resistance, including the four key SNP markers: Gm13_29459954_ss715614803, Gm13_29751552_ss715614847, Gm13_30293949_ss715614951, and Gm13_30724301_ss715615024. Within this QTL, four candidate genes were identified: Glyma.13G194100, Glyma.13G184800, Glyma.13G184900, and Glyma.13G190800 (3Gg2). The genomic prediction (GP) accuracies ranged from 0.60 to 0.83 across three GWAS-derived SNP sets using five models, demonstrating the feasibility of GP for SMV-G1 resistance. These findings could provide a useful reference in soybean breeding targeting SMV-G1 resistance.

Improving genomic prediction of vitamin C content in spinach using GWAS-derived markers

Vitamin C (VC), also known as ascorbic acid and ascorbate, is a water-soluble antioxidant in plants that promotes skin health and immune function in humans. Spinach (Spinacia oleracea L.) is a leafy green widely consumed for its health benefits. Recent reports have shown that nutritional content, including VC, can be improved in spinach. However, due to its complex inheritance, new selection methods are needed to improve selection for cultivar development. In this study, single nucleotide polymorphism (SNP) markers identified by genome-wide association (GWAS) were used for genomic prediction (GP) to improve prediction accuracy (PA) for VC content in spinach. A set of 147,977 SNPs generated from whole genome resequencing was used for GWAS in a panel of 347 spinach genotypes by six GWAS models. Sixty-two SNP markers distributed on all six spinach chromosomes were associated with VC content. PA for the selection of VC content was estimated with fourteen random SNP sets across seven GP models. The results indicated that the PA can be > 40% after using 1,000 or more SNPs in six of the seven models tested; using GWAS-derived significant SNP markers PA increases to a high r-value up to 0.7 when using 62 associated SNP markers in Bayes ridge regression (BRR) model. Upon validation, identified accessions with high VC and high PA genomic selection model can be used in spinach breeding programs to develop high VC content cultivars.

Improvement of crop production in controlled environment agriculture through breeding

Controlled environment agriculture (CEA) represents one of the fastest-growing sectors of horticulture. Production in controlled environments ranges from highly controlled indoor environments with 100% artificial lighting (vertical farms or plant factories) to high-tech greenhouses with or without supplemental lighting, to simpler greenhouses and high tunnels. Although food production occurs in the soil inside high tunnels, most CEA operations use various hydroponic systems to meet crop irrigation and fertility needs. The expansion of CEA offers promise as a tool for increasing food production in and near urban systems as these systems do not rely on arable agricultural land. In addition, CEA offers resilience to climate instability by growing inside protective structures. Products harvested from CEA systems tend to be of high quality, both internal and external, and are sought after by consumers. Currently, CEA producers rely on cultivars bred for production in open-field agriculture. Because of high energy and other production costs in CEA, only a limited number of food crops have proven themselves to be profitable to produce. One factor contributing to this situation may be a lack of optimized cultivars. Indoor growing operations offer opportunities for breeding cultivars that are ideal for these systems. To facilitate breeding these specialized cultivars, a wide range of tools are available for plant breeders to help speed this process and increase its efficiency. This review aims to cover breeding opportunities and needs for a wide range of horticultural crops either already being produced in CEA systems or with potential for CEA production. It also reviews many of the tools available to breeders including genomics-informed breeding, marker-assisted selection, precision breeding, high-throughput phenotyping, and potential sources of germplasm suitable for CEA breeding. The availability of published genomes and trait-linked molecular markers should enable rapid progress in the breeding of CEA-specific food crops that will help drive the growth of this industry.

The whole-genome dissection of root system architecture provides new insights for the genetic improvement of alfalfa (Medicago sativa L.)

Appropriate root system architecture (RSA) can improve alfalfa yield, yet its genetic basis remains largely unexplored. This study evaluated six RSA traits in 171 alfalfa genotypes grown under controlled greenhouse conditions. We also analyzed five yield-related traits in normal and drought stress environments and found a significant correlation (0.50) between root dry weight (RDW) and alfalfa dry weight under normal conditions (N_DW). A genome-wide association study (GWAS) was performed using 1 303 374 single-nucleotide polymorphisms (SNPs) to explore the relationships between RSA traits. Sixty significant SNPs (-log 10 (P) ≥ 5) were identified, with genes within the 50 kb upstream and downstream ranges primarily enriched in GO terms related to root development, hormone synthesis, and signaling, as well as morphological development. Further analysis identified 19 high-confidence candidate genes, including AUXIN RESPONSE FACTORs (ARFs), LATERAL ORGAN BOUNDARIES-DOMAIN (LBD), and WUSCHEL-RELATED HOMEOBOX (WOX). We verified that the forage dry weight under both normal and drought conditions exhibited significant differences among materials with different numbers of favorable haplotypes. Alfalfa containing more favorable haplotypes exhibited higher forage yields, whereas favorable haplotypes were not subjected to human selection during alfalfa breeding. Genomic prediction (GP) utilized SNPs from GWAS and machine learning for each RSA trait, achieving prediction accuracies ranging from 0.70 for secondary root position (SRP) to 0.80 for root length (RL), indicating robust predictive capability across the assessed traits. These findings provide new insights into the genetic underpinnings of root development in alfalfa, potentially informing future breeding strategies aimed at improving yield.

Unraveling genetic mysteries: A comprehensive review of GWAS and DNA insights in animal and plant pathosystems

DNA serves as the carrier of genetic information, with sequence variations playing a pivotal role in defining hereditary traits. Genome-Wide Association Studies (GWAS) facilitate the investigation of the links between genetic variations and phenotypes, significantly influencing biological research, particularly in animal and plant pathology. By identifying genetic markers associated with specific traits or diseases, GWAS enhances our understanding of host-pathogen interactions and improves disease-resistant breeding strategies. It has been vital in revealing the genetic basis of disease resistance, pinpointing key genes and DNA loci, which enrich genetic resources for breeding programs and deepen our knowledge of disease resistance mechanisms at the DNA level. Additionally, GWAS contributes to pathogen population genetics, facilitating a thorough exploration of pathogen virulence. Integrating GWAS with marker-assisted selection enhances breeding efficiency and precision in selecting for disease-resistant traits. While previous research has largely focused on host genetics, the genetic variation of pathogens is equally significant. Notably, reports integrating animal and plant pathosystems are still lacking. Given the importance of these systems, this review summarizes key advancements in this field, addresses current challenges, and proposes future directions, thereby offering a vital reference for ongoing research.

Hybrid Prediction in Horticulture Crop Breeding: Progress and Challenges

In the context of rapidly increasing population and diversified market demands, the steady improvement of yield and quality in horticultural crops has become an urgent challenge that modern breeding efforts must tackle. Heterosis, a pivotal theoretical foundation for plant breeding, facilitates the creation of superior hybrids through crossbreeding and selection among a variety of parents. However, the vast number of potential hybrids presents a significant challenge for breeders in efficiently predicting and selecting the most promising candidates. The development and refinement of effective hybrid prediction methods have long been central to research in this field. This article systematically reviews the advancements in hybrid prediction for horticultural crops, including the roles of marker-assisted breeding and genomic prediction in phenotypic forecasting. It also underscores the limitations of some predictors, like genetic distance, which do not consistently offer reliable hybrid predictions. Looking ahead, it explores the integration of phenomics with genomic prediction technologies as a means to elevate prediction accuracy within actual breeding programs.

Genomic selection for crop improvement in fruits and vegetables: a systematic scoping review

To ensure the nutritional needs of an expanding global population, it is crucial to optimize the growing capabilities and breeding values of fruit and vegetable crops. While genomic selection, initially implemented in animal breeding, holds tremendous potential, its utilization in fruit and vegetable crops remains underexplored. In this systematic review, we reviewed 63 articles covering genomic selection and its applications across 25 different types of fruit and vegetable crops over the last decade. The traits examined were directly related to the edible parts of the crops and carried significant economic importance. Comparative analysis with WHO/FAO data identified potential economic drivers underlying the study focus of some crops and highlighted crops with potential for further genomic selection research and application. Factors affecting genomic selection accuracy in fruit and vegetable studies are discussed and suggestions made to assist in their implementation into plant breeding schemes. Genetic gain in fruits and vegetables can be improved by utilizing genomic selection to improve selection intensity, accuracy, and integration of genetic variation. However, the reduction of breeding cycle times may not be beneficial in crops with shorter life cycles such as leafy greens as compared to fruit trees. There is an urgent need to integrate genomic selection methods into ongoing breeding programs and assess the actual genomic estimated breeding values of progeny resulting from these breeding programs against the prediction models.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01497-2.

Genomic prediction for agronomic traits in a diverse Flax (Linum usitatissimum L.) germplasm collection

Breeding programs require exhaustive phenotyping of germplasms, which is time-demanding and expensive. Genomic prediction helps breeders harness the diversity of any collection to bypass phenotyping. Here, we examined the genomic prediction's potential for seed yield and nine agronomic traits using 26,171 single nucleotide polymorphism (SNP) markers in a set of 337 flax (Linum usitatissimum L.) germplasm, phenotyped in five environments. We evaluated 14 prediction models and several factors affecting predictive ability based on cross-validation schemes. Models yielded significant variation among predictive ability values across traits for the whole marker set. The ridge regression (RR) model covering additive gene action yielded better predictive ability for most of the traits, whereas it was higher for low heritable traits by models capturing epistatic gene action. Marker subsets based on linkage disequilibrium decay distance gave significantly higher predictive abilities to the whole marker set, but for randomly selected markers, it reached a plateau above 3000 markers. Markers having significant association with traits improved predictive abilities compared to the whole marker set when marker selection was made on the whole population instead of the training set indicating a clear overfitting. The correction for population structure did not increase predictive abilities compared to the whole collection. However, stratified sampling by picking representative genotypes from each cluster improved predictive abilities. The indirect predictive ability for a trait was proportionate to its correlation with other traits. These results will help breeders to select the best models, optimum marker set, and suitable genotype set to perform an indirect selection for quantitative traits in this diverse flax germplasm collection.

A Genome-Wide Association Study Reveals the Genetic Mechanisms of Nutrient Accumulation in Spinach

Spinach is a significant source of vitamins, minerals, and antioxidants. These nutrients make it delicious and beneficial for human health. However, the genetic mechanism underlying the accumulation of nutrients in spinach remains unclear. In this study, we analyzed the content of chlorophyll a, chlorophyll b, oxalate, nitrate, crude fiber, soluble sugars, manganese, copper, and iron in 62 different spinach accessions. Additionally, 3,356,182 high-quality, single-nucleotide polymorphisms were found using resequencing and used in a genome-wide association study. A total of 2077 loci were discovered that significantly correlated with the concentrations of the nutritional elements. Data mining identified key genes in these intervals for four traits: chlorophyll, oxalate, soluble sugar, and Fe. Our study provides insights into the genetic architecture of nutrient variation and facilitates spinach breeding for good nutrition.

Genome-Wide Association Study and Genomic Prediction of Fusarium Wilt Resistance in Common Bean Core Collection

The common bean (Phaseolus vulgaris L.) is a globally cultivated leguminous crop. Fusarium wilt (FW), caused by Fusarium oxysporum f. sp. phaseoli (Fop), is a significant disease leading to substantial yield loss in common beans. Disease-resistant cultivars are recommended to counteract this. The objective of this investigation was to identify single nucleotide polymorphism (SNP) markers associated with FW resistance and to pinpoint potential resistant common bean accessions within a core collection, utilizing a panel of 157 accessions through the Genome-wide association study (GWAS) approach with TASSEL 5 and GAPIT 3. Phenotypes for Fop race 1 and race 4 were matched with genotypic data from 4740 SNPs of BARCBean6K_3 Infinium Bea Chips. After ranking the 157-accession panel and revealing 21 Fusarium wilt-resistant accessions, the GWAS pinpointed 16 SNPs on chromosomes Pv04, Pv05, Pv07, Pv8, and Pv09 linked to Fop race 1 resistance, 23 SNPs on chromosomes Pv03, Pv04, Pv05, Pv07, Pv09, Pv10, and Pv11 associated with Fop race 4 resistance, and 7 SNPs on chromosomes Pv04 and Pv09 correlated with both Fop race 1 and race 4 resistances. Furthermore, within a 30 kb flanking region of these associated SNPs, a total of 17 candidate genes were identified. Some of these genes were annotated as classical disease resistance protein/enzymes, including NB-ARC domain proteins, Leucine-rich repeat protein kinase family proteins, zinc finger family proteins, P-loopcontaining nucleoside triphosphate hydrolase superfamily, etc. Genomic prediction (GP) accuracy for Fop race resistances ranged from 0.26 to 0.55. This study advanced common bean genetic enhancement through marker-assisted selection (MAS) and genomic selection (GS) strategies, paving the way for improved Fop resistance.

Plant genome resequencing and population genomics: Current status and future prospects

Advances in DNA sequencing technology have sparked a genomics revolution, driving breakthroughs in plant genetics and crop breeding. Recently, the focus has shifted from cataloging genetic diversity in plants to exploring their functional significance and delivering beneficial alleles for crop improvement. This transformation has been facilitated by the increasing adoption of whole-genome resequencing. In this review, we summarize the current progress of population-based genome resequencing studies and how these studies affect crop breeding. A total of 187 land plants from 163 countries have been resequenced, comprising 54 413 accessions. As part of resequencing efforts 367 traits have been surveyed and 86 genome-wide association studies have been conducted. Economically important crops, particularly cereals, vegetables, and legumes, have dominated the resequencing efforts, leaving a gap in 49 orders, including Lycopodiales, Liliales, Acorales, Austrobaileyales, and Commelinales. The resequenced germplasm is distributed across diverse geographic locations, providing a global perspective on plant genomics. We highlight genes that have been selected during domestication, or associated with agronomic traits, and form a repository of candidate genes for future research and application. Despite the opportunities for cross-species comparative genomics, many population genomic datasets are not accessible, impeding secondary analyses. We call for a more open and collaborative approach to population genomics that promotes data sharing and encourages contribution-based credit policy. The number of plant genome resequencing studies will continue to rise with the decreasing DNA sequencing costs, coupled with advances in analysis and computational technologies. This expansion, in terms of both scale and quality, holds promise for deeper insights into plant trait genetics and breeding design.

Multiple-statistical genome-wide association analysis and genomic prediction of fruit aroma and agronomic traits in peaches

'Chinese Cling' is an important founder in peach breeding history due to the pleasant flavor. Genome-wide association studies (GWAS) combined with genomic selection are promising tools in fruit tree breeding, as there is a considerable time lapse between crossing and release of a cultivar. In this study, 242 peaches from Shanghai germplasm were genotyped with 145 456 single-nucleotide polymorphisms (SNPs). The six agronomic traits of fruit flesh color, fruit shape, fruit hairiness, flower type, pollen sterility, and soluble solids content, along with 14 key volatile odor compounds (VOCs), were recorded for multiple-statistical GWAS. Except the reported candidate genes, six novel genes were identified as associated with these traits. Thirty-nine significant SNPs were associated with eight VOCs. The putative candidate genes were confirmed for VOCs by RNA-seq, including three genes in the biosynthesis pathway found to be associated with linalool, soluble solids content, and cis-3-hexenyl acetate. Multiple-trait genomic prediction enhanced the predictive ability for γ-decalactone to 0.7415 compared with the single-trait model value of 0.1017. One PTS1-SSR marker was designed to predict the linalool content, and the favorable genotype 187/187 was confirmed, mainly existing in the 'Shanghai Shuimi' landrace. Overall, our findings will be helpful in determining peach accessions with the ideal phenotype and show the potential of multiple-trait genomic prediction to improve accuracy for highly correlated genetic traits. The diagnostic marker will be valuable for the breeder to bridge the gap between quantitative trait loci and marker-assisted selection for developing strong-aroma cultivars.

Skim resequencing finely maps the downy mildew resistance loci RPF2 and RPF3 in spinach cultivars whale and Lazio

Commercial production of spinach (Spinacia oleracea L.) is centered in California and Arizona in the US, where downy mildew caused by Peronospora effusa is the most destructive disease. Nineteen typical races of P. effusa have been reported to infect spinach, with 16 identified after 1990. The regular appearance of new pathogen races breaks the resistance gene introgressed in spinach. We attempted to map and delineate the RPF2 locus at a finer resolution, identify linked single nucleotide polymorphism (SNP) markers, and report candidate downy mildew resistance (R) genes. Progeny populations segregating for RPF2 locus derived from resistant differential cultivar Lazio were infected using race 5 of P. effusa and were used to study for genetic transmission and mapping analysis in this study. Association analysis performed with low coverage whole genome resequencing-generated SNP markers mapped the RPF2 locus between 0.47 to 1.46 Mb of chromosome 3 with peak SNP (Chr3_1, 221, 009) showing a LOD value of 61.6 in the GLM model in TASSEL, which was within 1.08 Kb from Spo12821, a gene that encodes CC-NBS-LRR plant disease resistance protein. In addition, a combined analysis of progeny panels of Lazio and Whale segregating for RPF2 and RPF3 loci delineated the resistance section in chromosome 3 between 1.18-1.23 and 1.75-1.76 Mb. This study provides valuable information on the RPF2 resistance region in the spinach cultivar Lazio compared to RPF3 loci in the cultivar Whale. The RPF2 and RPF3 specific SNP markers, plus the resistant genes reported here, could add value to breeding efforts to develop downy mildew resistant cultivars in the future.

Combining GS-assisted GWAS and transcriptome analysis to mine candidate genes for nitrogen utilization efficiency in Populus cathayana

BACKGROUND

Forest trees such as poplar, shrub willow, et al. are essential natural resources for sustainable and renewable energy production, and their wood can reduce dependence on fossil fuels and reduce environmental pollution. However, the productivity of forest trees is often limited by the availability of nitrogen (N), improving nitrogen use efficiency (NUE) is an important way to address it. Currently, NUE genetic resources are scarce in forest tree research, and more genetic resources are urgently needed.

RESULTS

Here, we performed genome-wide association studies (GWAS) using the mixed linear model (MLM) to identify genetic loci regulating growth traits in Populus cathayana at two N levels, and attempted to enhance the signal strength of single nucleotide polymorphism (SNP) detection by performing genome selection (GS) assistance GWAS. The results of the two GWAS analyses identified 55 and 40 SNPs that were respectively associated with plant height (PH) and ground diameter (GD), and 92 and 69 candidate genes, including 30 overlapping genes. The prediction accuracy of the GS model (rrBLUP) for phenotype exceeds 0.9. Transcriptome analysis of 13 genotypes under two N levels showed that genes related to carbon and N metabolism, amino acid metabolism, energy metabolism, and signal transduction were differentially expressed in the xylem of P. cathayana under N treatment. Furthermore, we observed strong regional patterns in gene expression levels of P. cathayana, with significant differences between different regions. Among them, P. cathayana in Longquan region exhibited the highest response to N. Finally, through weighted gene co-expression network analysis (WGCNA), we identified a module closely related to the N metabolic process and eight hub genes.

CONCLUSIONS

Integrating the GWAS, RNA-seq and WGCNA data, we ultimately identified four key regulatory genes (PtrNAC123, PtrNAC025, Potri.002G233100, and Potri.006G236200) involved in the wood formation process, and they may affect P. cathayana growth and wood formation by regulating nitrogen metabolism. This study will provide strong evidence for N regulation mechanisms, and reliable genetic resources for growth and NUE genetic improvement in poplar.

Crop germplasm: Current challenges, physiological-molecular perspective, and advance strategies towards development of climate-resilient crops

Germplasm is a long-term resource management mission and investment for civilization. An estimated ∼7.4 million accessions are held in 1750 plant germplasm centres around the world; yet, only 2% of these assets have been utilized as plant genetic resources (PGRs). According to recent studies, the current food yield trajectory will be insufficient to feed the world's population in 2050. Additionally, possible negative effects in terms of crop failure because of climate change are already being experienced across the world. Therefore, it is necessary to reconciliation of research advancement and innovation of practices for further exploration of the potential of crop germplasm especially for the complex traits associated with yield such as water- and nitrogen use efficiency. In this review, we tried to address current challenges, research gaps, physiological and molecular aspects of two broad spectrum complex traits such as water- and nitrogen-use efficiency, and advanced integrated strategies that could provide a platform for combined stress management for climate-smart crop development. Additionally, recent development in technologies that are directly related to germplasm characterization was highlighted for further molecular utilization towards the development of elite varieties.

Mapping and selection of downy mildew resistance in spinach cv. whale by low coverage whole genome sequencing

Spinach (Spinacia oleracea) is a popular leafy vegetable crop and commercial production is centered in California and Arizona in the US. The oomycete Peronospora effusa causes the most important disease in spinach, downy mildew. A total of nineteen races of P. effusa are known, with more than 15 documented in the last three decades, and the regular emergence of new races is continually overcoming the genetic resistance to the pathogen. This study aimed to finely map the downy mildew resistance locus RPF3 in spinach, identify single nucleotide polymorphism (SNP) markers associated with the resistance, refine the candidate genes responsible for the resistance, and evaluate the prediction performance using multiple machine learning genomic prediction (GP) methods. Segregating progeny population developed from a cross of resistant cultivar Whale and susceptible cultivar Viroflay to race 5 of P. effusa was inoculated under greenhouse conditions to determine downy mildew disease response across the panel. The progeny panel and the parents were resequenced at low coverage (1x) to identify genome wide SNP markers. Association analysis was performed using disease response phenotype data and SNP markers in TASSEL, GAPIT, and GENESIS programs and mapped the race 5 resistance loci (RPF3) to 1.25 and 2.73 Mb of Monoe-Viroflay chromosome 3 with the associated SNP in the 1.25 Mb region was 0.9 Kb from the NBS-LRR gene SOV3g001250. The RPF3 locus in the 1.22-1.23 Mb region of Sp75 chromosome 3 is 2.41-3.65 Kb from the gene Spo12821 annotated as NBS-LRR disease resistance protein. This study extended our understanding of the genetic basis of downy mildew resistance in spinach cultivar Whale and mapped the RPF3 resistance loci close to the NBS-LRR gene providing a target to pursue functional validation. Three SNP markers efficiently selected resistance based on multiple genomic selection (GS) models. The results from this study have added new genomic resources, generated an informed basis of the RPF3 locus resistant to spinach downy mildew pathogen, and developed markers and prediction methods to select resistant lines.

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.

Genetic dissection of nitrogen induced changes in the shoot and root biomass of spinach

Efficient partitioning of above and below-ground biomass in response to nitrogen (N) is critical to the productivity of plants under sub-optimal conditions. It is particularly essential in vegetable crops like spinach with shallow root systems, a short growth cycle, and poor nitrogen use efficiency. In this study, we conducted a genome-wide association study (GWAS) to explore N-induced changes using spinach accessions with diverse genetic backgrounds. We evaluated phenotypic variations as percent changes in the shoot and root biomass in response to N using 201 spinach accessions grown in randomized complete blocks design in a soilless media under a controlled environment. A GWAS was performed for the percent changes in the shoot and root biomass in response to N in the 201 spinach accessions using 60,940 whole-genome resequencing generated SNPs. Three SNP markers, chr4_28292655, chr6_1531056, and chr6_37966006 on chromosomes 4 and 6, were significantly associated with %change in root weight, and two SNP markers, chr2_18480277 and chr4_47598760 on chromosomes 2 and 4, were significantly associated with % change shoot weight. The outcome of this study established a foundation for genetic studies needed to improve the partitioning of total biomass and provided a resource to identify molecular markers to enhance N uptake via marker-assisted selection or genomic selection in spinach breeding programs.