Characterization of a resistance locus (Pfs-1) to the spinach downy mildew pathogen (Peronospora farinosa f. sp. spinaciae) and development of a molecular marker linked to Pfs-1

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.

Fine Mapping and Identification of a Candidate Gene of Downy Mildew Resistance, RPF2, in Spinach (Spinacia oleracea L.)

Downy mildew is a major threat to the economic value of spinach. The most effective approach to managing spinach downy mildew is breeding cultivars with resistance genes. The resistance allele RPF2 is effective against races 1-10 and 15 of Peronospora farinosa f. sp. Spinaciae (P. effusa) and is widely used as a resistance gene. However, the gene and the linked marker of RPF2 remain unclear, which limit its utilization. Herein, we located the RPF2 gene in a 0.61 Mb region using a BC₁ population derived from Sp39 (rr) and Sp62 (RR) cultivars via kompetitive allele specific PCR (KASP) markers. Within this region, only one R gene, Spo12821, was identified based on annotation information. The amino acid sequence analysis showed that there were large differences in the length of the LRR domain between the parents. Additionally, a molecular marker, RPF2-IN12821, was developed based on the sequence variation in the Spo12821, and the evaluation in the BC₁ population produced a 100% match with resistance/susceptibility. The finding of the study could be valuable for improving our understanding of the genetic basis of resistance against the downy mildew pathogen and breeding resistance lines in the future.

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.

Genomic analyses provide insights into spinach domestication and the genetic basis of agronomic traits

Spinach is a nutritious leafy vegetable belonging to the family Chenopodiaceae. Here we report a high-quality chromosome-scale reference genome assembly of spinach and genome resequencing of 305 cultivated and wild spinach accessions. Reconstruction of ancestral Chenopodiaceae karyotype indicates substantial genome rearrangements in spinach after its divergence from ancestral Chenopodiaceae, coinciding with high repeat content in the spinach genome. Population genomic analyses provide insights into spinach genetic diversity and population differentiation. Genome-wide association studies of 20 agronomical traits identify numerous significantly associated regions and candidate genes for these traits. Domestication sweeps in the spinach genome are identified, some of which are associated with important traits (e.g., leaf phenotype, bolting and flowering), demonstrating the role of artificial selection in shaping spinach phenotypic evolution. This study provides not only insights into the spinach evolution and domestication but also valuable resources for facilitating spinach breeding.

Spinach is a nutritious leafy vegetable growing worldwide. Here, the authors report a high-quality chromosome-scale reference genome assembly of spinach and genome resequencing of 305 accessions, and provide insights into spinach domestication and the genetic basis of agronomic traits.

High resolution mapping and candidate gene identification of downy mildew race 16 resistance in spinach

Background

Downy mildew, the most devastating disease of spinach (Spinacia oleracea L.), is caused by the oomycete Peronospora effusa [=P. farinosa f. sp. spinaciae]. The P. effusa shows race specificities to the resistant host and comprises 19 reported races and many novel isolates. Sixteen new P. effusa races were identified during the past three decades, and the new pathogen races are continually overcoming the genetic resistances used in commercial cultivars. A spinach breeding population derived from the cross between cultivars Whale and Lazio was inoculated with P. effusa race 16 in an environment-controlled facility; disease response was recorded and genotyped using genotyping by sequencing (GBS). The main objective of this study was to identify resistance-associated single nucleotide polymorphism (SNP) markers from the cultivar Whale against the P. effusa race 16.

Results

Association analysis conducted using GBS markers identified six significant SNPs (S3_658,306, S3_692697, S3_1050601, S3_1227787, S3_1227802, S3_1231197). The downy mildew resistance locus from cultivar Whale was mapped to a 0.57 Mb region on chromosome 3, including four disease resistance candidate genes (Spo12736, Spo12784, Spo12908, and Spo12821) within 2.69–11.28 Kb of the peak SNP.

Conclusions

Genomewide association analysis approach was used to map the P. effusa race 16 resistance loci and identify associated SNP markers and the candidate genes. The results from this study could be valuable in understanding the genetic basis of downy mildew resistance, and the SNP marker will be useful in spinach breeding to select resistant lines.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

Genome-wide simple sequence repeats (SSR) markers discovered from whole-genome sequence comparisons of multiple spinach accessions

The availability of well-assembled genome sequences and reduced sequencing costs have enabled the resequencing of many additional accessions in several crops, thus facilitating the rapid discovery and development of simple sequence repeat (SSR) markers. Although the genome sequence of inbred spinach line Sp75 is available, previous efforts have resulted in a limited number of useful SSR markers. Identification of additional polymorphic SSR markers will support genetics and breeding research in spinach. This study aimed to use the available genomic resources to mine and catalog a large number of polymorphic SSR markers. A search for SSR loci on six chromosome sequences of spinach line Sp75 using GMATA identified a total of 42,155 loci with repeat motifs of two to six nucleotides in the Sp75 reference genome. Whole-genome sequences (30x) of additional 21 accessions were aligned against the chromosome sequences of the reference genome and in silico genotyped using the HipSTR program by comparing and counting repeat numbers variation across the SSR loci among the accessions. The HipSTR program generated SSR genotype data were filtered for monomorphic and high missing loci, and a final set of the 5986 polymorphic SSR loci were identified. The polymorphic SSR loci were present at a density of 12.9 SSRs/Mb and were physically mapped. Out of 36 randomly selected SSR loci for validation, two failed to amplify, while the remaining were all polymorphic in a set of 48 spinach accessions from 34 countries. Genetic diversity analysis performed using the SSRs allele score data on the 48 spinach accessions showed three main population groups. This strategy to mine and develop polymorphic SSR markers by a comparative analysis of the genome sequences of multiple accessions and computational genotyping of the candidate SSR loci eliminates the need for laborious experimental screening. Our approach increased the efficiency of discovering a large set of novel polymorphic SSR markers, as demonstrated in this report.

Insight into morphological and molecular variations across Iranian spinach landraces

Spinach is salt-tolerant leafy vegetable stemmed from central Asia. It has been well adapted to different climates of Iran. We aimed to study genetic diversity between several landraces as a prerequisite for crop improvement programs using molecular markers including microsatellites (SSRs) and morphological traits. Genetic diversity was studied among 22 spinach landraces using morphological and molecular tools. We developed 17 genic and genomic SSR markers based on the information acquired from NCBI resources. Morphological evaluation indicated high variability for economic traits including leaf color, leaf thickness, leaf wrinkle and seed type, across Iranian landraces. The molecular results exhibited that 12 out of 17 primer pairs successfully amplified genomic DNA with explicit bands. The results verified that genic markers were superior to genomic markers to detect polymorphism and genetic diversity. In this regard, PIC for genomic and genic SSRs was in the range of 0.4616-0.6621 and 0.5188-0.7394, respectively. Polymorphic genic SSRs were identified to be directly and indirectly involved in biotic/abiotic stresses. High degree of polymorphism, which was detected across the landraces by genic SSRs, could assist us to select 11 landraces for the second experiment. The second experiment was designed to evaluate the response of selected landraces to salinity stress. The results confirmed genetic variability among the landraces in terms of salinity tolerance. A highly diverse germplasm of Iranian spinach based on molecular and morphological characteristics along with the tolerance to oxidative stress provides an ample opportunity for plant breeders to select superior genotypes.

Genome Wide Association Studies in Multiple Spinach Breeding Populations Refine Downy Mildew Race 13 Resistance Genes

Downy mildew, caused by the oomycete Peronospora effusa, is the most economically important disease on spinach. Fourteen new races of P. effusa have been identified in the last three decades. The frequent emergence of new races of P. effusa continually overcome the genetic resistance to the pathogen. The objectives of this research were to more clearly map the downy mildew resistance locus RPF1 in spinach, to identify single nucleotide polymorphism (SNP) markers associated with the resistance, and to refine the candidate genes responsible for the resistance. Progeny from populations generated from crosses of cultivars resistant (due to RPF1) to race 13 of P. effusa (Swan, T-Bird, Squirrel, and Tonga) with race 13 susceptible cultivars (Whale and Polka) were inoculated and the downy mildew disease response determined. Association analysis was performed in TASSEL, GAPIT, PLINK, and GENESIS programs using SNP markers identified from genotyping by sequencing (GBS). Association analysis mapped the race 13 resistance loci (RPF1) to positions 0.39, 0.69, 0.94-0.98, and 1.2 Mb of chromosome 3. The associated SNPs were within 1-7 kb of the disease resistance genes Spo12784, Spo12719, Spo12905, and Spo12821, and 11-18 Kb from Spo12903. This study extended our understanding of the genetic basis of downy mildew resistance in spinach and provided the most promising candidate genes Spo12784 and Spo12903 near the RPF1 locus, to pursue functional validation. The SNP markers may be used to select for the resistant lines to improve genetic resistance against the downy mildew pathogen and in developing durably resistant cultivars.

Spinach Downy Mildew: Advances in Our Understanding of the Disease Cycle and Prospects for Disease Management

Downy mildew on spinach is caused by Peronospora effusa, an oomycete pathogen that poses a challenge to spinach production worldwide, especially in organic production. Following infection, P. effusa produces abundant amounts of asexual sporangia. Sporangia become windborne and initiate new infections locally or distantly, leading to widespread epidemics. Oospores produced from the union of opposite mating types have been observed within infected leaves and seeds and may remain viable for many years. Sexual reproduction increases the genetic diversity of P. effusa through sexual recombination, and thus, the movement of oospores on seed has likely fueled the rapid explosion of new pathotypes in different regions of the world over the past 20 years. This review summarizes recent advances in spinach downy mildew research, especially in light of the findings of oospores in contemporary commercial spinach seed lots as well as their germination. Knowledge of the role of the oospores and other aspects of the disease cycle can directly translate into new and effective disease management strategies.

Fine mapping and candidate gene screening of the downy mildew resistance gene RPF1 in Spinach

A SLAF-BSA approach was used to locate the RPF1 locus. The three most likely candidate genes were identified which provide a basic for cloning the resistance gene at the RPF1 locus. Spinach downy mildew is a globally devastating oomycete disease. The use of downy mildew resistance genes constitutes the most effective approach for disease management. Hence, the objective of the present study was to fine map the first-reported resistance locus RPF1. The resistance allele at this resistance locus was effective against races 1-7, 9, 11, 13, and 15 of Peronospora farinosa f. sp. spinaciae (P. effusa). The approach fine mapped RPF1 using specific-locus amplified fragment sequencing (SLAF-Seq) technology combined with bulked segregant analysis. A 1.72 Mb region localized on chromosome 3 was found to contain RPF1 based on association analysis. After screening recombinants with the SLAF markers within the region, the region was narrowed down to 0.89 Mb. Within this region, 14 R genes were identified based on the annotation information. To identify the genes involved in resistance, resequencing of two resistant inbred lines (12S2 and 12S3) and three susceptible inbred lines (12S1, 12S4, and 10S2) was performed. The three most likely candidate genes were identified via amino acid sequence analysis and conserved domain analysis between resistant and susceptible inbred lines. These included Spo12729, encoding a receptor-like protein, and Spo12784 and Spo12903, encoding a nucleotide-binding site and leucine-rich repeat domains. Additionally, based on the sequence variation in the three genes between the resistant and susceptible lines, molecular markers were developed for marker-assisted selection. The results could be valuable in cloning the RPF1 alleles and improving our understanding of the interaction between the host and pathogen.

Population structure, genetic diversity and downy mildew resistance among Ocimum species germplasm

BACKGROUND

The basil (Ocimum spp.) genus maintains a rich diversity of phenotypes and aromatic volatiles through natural and artificial outcrossing. Characterization of population structure and genetic diversity among a representative sample of this genus is severely lacking. Absence of such information has slowed breeding efforts and the development of sweet basil (Ocimum basilicum L.) with resistance to the worldwide downy mildew epidemic, caused by the obligate oomycete Peronospora belbahrii. In an effort to improve classification of relationships 20 EST-SSR markers with species-level transferability were developed and used to resolve relationships among a diverse panel of 180 Ocimum spp. accessions with varying response to downy mildew.

RESULTS

Results obtained from nested Bayesian model-based clustering, analysis of molecular variance and unweighted pair group method using arithmetic average (UPGMA) analyses were synergized to provide an updated phylogeny of the Ocimum genus. Three (major) and seven (sub) population (cluster) models were identified and well-supported (P < 0.001) by PhiPT (Φ_PT) values of 0.433 and 0.344, respectively. Allelic frequency among clusters supported previously developed hypotheses of allopolyploid genome structure. Evidence of cryptic population structure was demonstrated for the k1 O. basilicum cluster suggesting prevalence of gene flow. UPGMA analysis provided best resolution for the 36-accession, DM resistant k3 cluster with consistently strong bootstrap support. Although the k3 cluster is a rich source of DM resistance introgression of resistance into the commercially important k1 accessions is impeded by reproductive barriers as demonstrated by multiple sterile F1 hybrids. The k2 cluster located between k1 and k3, represents a source of transferrable tolerance evidenced by fertile backcross progeny. The 90-accession k1 cluster was largely susceptible to downy mildew with accession 'MRI' representing the only source of DM resistance.

CONCLUSIONS

High levels of genetic diversity support the observed phenotypic diversity among Ocimum spp. accessions. EST-SSRs provided a robust evaluation of molecular diversity and can be used for additional studies to increase resolution of genetic relationships in the Ocimum genus. Elucidation of population structure and genetic relationships among Ocimum spp. germplasm provide the foundation for improved DM resistance breeding strategies and more rapid response to future disease outbreaks.

New Races and Novel Strains of the Spinach Downy Mildew Pathogen Peronospora effusa

Downy mildew disease, caused by Peronospora effusa (=P. farinosa f. sp. spinaciae [Pfs]), is the most economically important disease of spinach. Current high-density fresh-market spinach production provides conducive conditions for disease development, and downy mildew frequently forces growers to harvest early owing to disease development, to cull symptomatic leaves prior to harvest, or to abandon the field if the disease is too severe. The use of resistant cultivars to manage downy mildew, particularly on increasing acreages of organic spinach production, applies strong selection pressure on the pathogen, and many new races of Pfs have been identified in recent years in spinach production areas worldwide. To monitor the virulence diversity in the Pfs population, downy mildew samples were collected from spinach production areas and tested for race identification based on the disease reactions of a standard set of international spinach differentials. Two new races (designated races 15 and 16) and eight novel strains were identified between 2013 and 2017. The disease reaction of Pfs 15 was similar to race 4, except race 4 could not overcome the resistance imparted by the RPF9 locus. Several resistance loci (RPF1, 2, 4, and 6) were effective in preventing disease caused by Pfs 15. The race Pfs 16 could overcome several resistance loci (RPF2, 4, 5, 9, and 10) but not others (RPF1, 3, 6, and 7). One novel strain (UA1014) could overcome the resistance of spinach resistant loci RPF1 to RPF7 but only infected the cotyledons and not the true leaves of certain cultivars. A new set of near-isogenic lines has been developed and evaluated for disease reactions to the new races and novel strains as differentials. None of the 360 U.S. Department of Agriculture spinach germplasm accessions tested were resistant to Pfs 16 or UA1014. A survey of isolates over several years highlighted the dynamic nature of the virulence diversity of the Pfs population. Identification of virulence diversity and evaluation of the genetics of resistance to Pfs will continue to allow for a more effective disease management strategy through resistance gene deployment.

A first linkage map and downy mildew resistance QTL discovery for sweet basil (Ocimum basilicum) facilitated by double digestion restriction site associated DNA sequencing (ddRADseq)

Limited understanding of sweet basil (Ocimum basilicum L.) genetics and genome structure has reduced efficiency of breeding strategies. This is evidenced by the rapid, worldwide dissemination of basil downy mildew (Peronospora belbahrii) in the absence of resistant cultivars. In an effort to improve available genetic resources, expressed sequence tag simple sequence repeat (EST-SSR) and single nucleotide polymorphism (SNP) markers were developed and used to genotype the MRI x SB22 F2 mapping population, which segregates for response to downy mildew. SNP markers were generated from genomic sequences derived from double digestion restriction site associated DNA sequencing (ddRADseq). Disomic segregation was observed in both SNP and EST-SSR markers providing evidence of an O. basilicum allotetraploid genome structure and allowing for subsequent analysis of the mapping population as a diploid intercross. A dense linkage map was constructed using 42 EST-SSR and 1,847 SNP markers spanning 3,030.9 cM. Multiple quantitative trait loci (QTL) model (MQM) analysis identified three QTL that explained 37-55% of phenotypic variance associated with downy mildew response across three environments. A single major QTL, dm11.1 explained 21-28% of phenotypic variance and demonstrated dominant gene action. Two minor QTL dm9.1 and dm14.1 explained 5-16% and 4-18% of phenotypic variance, respectively. Evidence is provided for an additive effect between the two minor QTL and the major QTL dm11.1 increasing downy mildew susceptibility. Results indicate that ddRADseq-facilitated SNP and SSR marker genotyping is an effective approach for mapping the sweet basil genome.

Draft genome of spinach and transcriptome diversity of 120 Spinacia accessions

Spinach is an important leafy vegetable enriched with multiple necessary nutrients. Here we report the draft genome sequence of spinach (Spinacia oleracea, 2n=12), which contains 25,495 protein-coding genes. The spinach genome is highly repetitive with 74.4% of its content in the form of transposable elements. No recent whole genome duplication events are observed in spinach. Genome syntenic analysis between spinach and sugar beet suggests substantial inter- and intra-chromosome rearrangements during the Caryophyllales genome evolution. Transcriptome sequencing of 120 cultivated and wild spinach accessions reveals more than 420 K variants. Our data suggests that S. turkestanica is likely the direct progenitor of cultivated spinach and spinach domestication has a weak bottleneck. We identify 93 domestication sweeps in the spinach genome, some of which are associated with important agronomic traits including bolting, flowering and leaf numbers. This study offers insights into spinach evolution and domestication and provides resources for spinach research and improvement.

Spinach is an economically important vegetable crop but previous genomic resources were of limited use for comparative and functional analyses. Here, Xu et al. present a high quality draft spinach genome and transcriptome data for multiple Spinacia accessions providing insight into Caryophyllales genome evolution.

De novo and comparative transcriptome analysis of cultivated and wild spinach

Spinach (Spinacia oleracea L.) is an economically important green leafy vegetable crop. In this study, we performed deep transcriptome sequencing for nine spinach accessions: three from cultivated S. oleracea, three from wild S. turkestanica and three from wild S. tetrandra. A total of approximately 100 million high-quality reads were generated, which were de novo assembled into 72,151 unigenes with a total length of 46.5 Mb. By comparing sequences of these unigenes against different protein databases, nearly 60% of them were annotated and 50% could be assigned with Gene Ontology terms. A total of 387 metabolic pathways were predicted from the assembled spinach unigenes. From the transcriptome sequencing data, we were able to identify a total of ~320,000 high-quality single nucleotide polymorphisms (SNPs). Phylogenetic analyses using SNPs as well as gene expression profiles indicated that S. turkestanica was more closely related to the cultivated S. oleracea than S. tetrandra. A large number of genes involved in responses to biotic and abiotic stresses were found to be differentially expressed between the cultivated and wild spinach. Finally, an interactive online database (http://www.spinachbase.org) was developed to allow the research community to efficiently retrieve, query, mine and analyze our transcriptome dataset.

Identification of New Races and Deviating Strains of the Spinach Downy Mildew Pathogen Peronospora farinosa f. sp. spinaciae

Spinach downy mildew disease, caused by the obligate pathogen Peronospora farinosa f. sp. spinaciae, is the most economically important spinach (Spinacia oleracea) disease. New races of this pathogen have been emerging at a rapid rate over the last 15 years. This is likely due to production changes, particularly in California, such as high-density plantings and year-round spinach production. As of 2004, 10 races of P. farinosa f. sp. spinaciae had been identified, and the spinach resistance locus RPF2 provided resistance to races 1 to 10. Based on disease reactions on a set of spinach differentials containing six hypothesized resistance loci (RPF1-RPF6), races 11, 12, 13, and 14 of P. farinosa f. sp. spinaciae were characterized based on samples collected in the past 5 years as part of this study. Race 11, identified in 2008, could overcome the resistance of spinach cultivars resistant to races 1 to 10. Spinach resistance loci RPF1, RPF3, and RPF6 provided resistance to race 11. Race 12 was identified in 2009 and could overcome the resistances of the RPF1 and RPF2 loci. The RPF3 locus was effective against race 12. Race 13 was identified in 2010 and could overcome the resistance imparted by the RPF2 and RPF3 loci, whereas the RPF1 locus was effective against race 13. Race 14 was similar to race 12 and caused identical disease responses on the standard differentials but could be distinguished from race 12 by its ability to cause disease on a number of newly released cultivars, including 'Pigeon', 'Cello', and 'Celesta'. Five novel strains of P. farinosa f. sp. spinaciae were also identified. For example, isolate UA4711 of the pathogen, collected from Spain in 2011, was able to overcome the resistance imparted by the RPF1 and RPF3 loci, while RPF2 and RPF4 were effective against this strain. A total of 116 spinach cultivars, including 103 commercial lines and 13 differential cultivars, were evaluated for resistance to race 10 and the newly designated races 11, 12, 13, and 14.