Population genomics identifies genetic signatures of carrot domestication and improvement and uncovers the origin of high-carotenoid orange carrots

Here an improved carrot reference genome and resequencing of 630 carrot accessions were used to investigate carrot domestication and improvement. The study demonstrated that carrot was domesticated during the Early Middle Ages in the region spanning western Asia to central Asia, and orange carrot was selected during the Renaissance period, probably in western Europe. A progressive reduction of genetic diversity accompanied this process. Genes controlling circadian clock/flowering and carotenoid accumulation were under selection during domestication and improvement. Three recessive genes, at the REC, Or and Y2 quantitative trait loci, were essential to select for the high α- and β-carotene orange phenotype. All three genes control high α- and β-carotene accumulation through molecular mechanisms that regulate the interactions between the carotenoid biosynthetic pathway, the photosynthetic system and chloroplast biogenesis. Overall, this study elucidated carrot domestication and breeding history and carotenoid genetics at a molecular level.

The Gynandropsis gynandra genome provides insights into whole-genome duplications and the evolution of C4 photosynthesis in Cleomaceae

Gynandropsis gynandra (Cleomaceae) is a cosmopolitan leafy vegetable and medicinal plant, which has also been used as a model to study C4 photosynthesis due to its evolutionary proximity to C3 Arabidopsis (Arabidopsis thaliana). Here, we present the genome sequence of G. gynandra, anchored onto 17 main pseudomolecules with a total length of 740 Mb, an N50 of 42 Mb and 30,933 well-supported gene models. The G. gynandra genome and previously released genomes of C3 relatives in the Cleomaceae and Brassicaceae make an excellent model for studying the role of genome evolution in the transition from C3 to C4 photosynthesis. Our analyses revealed that G. gynandra and its C3 relative Tarenaya hassleriana shared a whole-genome duplication event (Gg-α), then an addition of a third genome (Th-α, +1×) took place in T. hassleriana but not in G. gynandra. Analysis of syntenic copy number of C4 photosynthesis-related gene families indicates that G. gynandra generally retained more duplicated copies of these genes than C3T. hassleriana, and also that the G. gynandra C4 genes might have been under positive selection pressure. Both whole-genome and single-gene duplication were found to contribute to the expansion of the aforementioned gene families in G. gynandra. Collectively, this study enhances our understanding of the polyploidy history, gene duplication and retention, as well as their impact on the evolution of C4 photosynthesis in Cleomaceae.

Genetics of destemming in pepper: A step towards mechanical harvesting

Introduction: The majority of peppers in the US for fresh market and processing are handpicked, and harvesting can account for 20–50% of production costs. Innovation in mechanical harvesting would increase availability; lower the costs of local, healthy vegetable products; and perhaps improve food safety and expand markets. Most processed peppers require removal of pedicels (stem and calyx) from the fruit, but lack of an efficient mechanical process for this operation has hindered adoption of mechanical harvest. In this paper, we present characterization and advancements in breeding green chile peppers for mechanical harvesting. Specifically, we describe inheritance and expression of an easy-destemming trait derived from the landrace UCD-14 that facilitates machine harvest of green chiles.

Methods: A torque gauge was used for measuring bending forces similar to those of a harvester and applied to two biparental populations segregating for destemming force and rate. Genotyping by sequencing was used to generate genetic maps for quantitative trait locus (QTL) analyses.

Results: A major destemming QTL was found on chromosome 10 across populations and environments. Eight additional population and/or environment-specific QTL were also identified. Chromosome 10 QTL markers were used to help introgress the destemming trait into jalapeño-type peppers. Low destemming force lines combined with improvements in transplant production enabled mechanical harvest of destemmed fruit at a rate of 41% versus 2% with a commercial jalapeńo hybrid. Staining for the presence of lignin at the pedicel/fruit boundary indicated the presence of an abscission zone and homologs of genes known to affect organ abscission were found under several QTL, suggesting that the easy-destemming trait may be due to the presence and activation of a pedicel/fruit abscission zone.

Conclusion: Presented here are tools to measure the easy-destemming trait, its physiological basis, possible molecular pathways, and expression of the trait in various genetic backgrounds. Mechanical harvest of destemmed mature green chile fruits was achieved by combining easy-destemming with transplant management.

GRF-GIF chimeric proteins enhance in vitro regeneration and Agrobacterium-mediated transformation efficiencies of lettuce (Lactuca spp.)

KEY MESSAGE

GRF-GIF chimeric proteins from multiple source species enhance in vitro regeneration in both wild and cultivated lettuce. In addition, they enhance regeneration in multiple types of lettuce including butterheads, romaines, and crispheads. The ability of plants to regenerate in vitro has been exploited for use in tissue culture systems for plant propagation, plant transformation, and genome editing. The success of in vitro regeneration is often genotype dependent and continues to be a bottleneck for Agrobacterium-mediated transformation and its deployment for improvement of some crop species. Manipulation of transcription factors that play key roles in plant development such as BABY BOOM, WUSCHEL, and GROWTH-REGULATING FACTORs (GRFs) has improved regeneration and transformation efficiencies in several plant species. Here, we compare the efficacy of GRF-GIF gene fusions from multiple species to boost regeneration efficiency and shooting frequency in four genotypes of wild and cultivated lettuce (Lactuca spp. L.). In addition, we show that GRF-GIFs with mutated miRNA 396 binding sites increase regeneration efficiency and shooting frequency when compared to controls. We also present a co-transformation strategy for increased transformation efficiency and recovery of transgenic plants harboring a gene of interest. This strategy will enhance the recovery of transgenic plants of other lettuce genotypes and likely other crops in the Compositae family.

High-quality chromosome-scale genomes facilitate effective identification of large structural variations in hot and sweet peppers

Pepper (Capsicum annuum) is an important vegetable crop that has been subjected to intensive breeding, resulting in limited genetic diversity, especially for sweet peppers. Previous studies have reported pepper draft genome assemblies using short read sequencing, but their capture of the extent of large structural variants (SVs), such as presence-absence variants (PAVs), inversions, and copy-number variants (CNVs) in the complex pepper genome falls short. In this study, we sequenced the genomes of representative sweet and hot pepper accessions by long-read and/or linked-read methods and advanced scaffolding technologies. First, we developed a high-quality reference genome for the sweet pepper cultivar 'Dempsey' and then used the reference genome to identify SVs in 11 other pepper accessions and constructed a graph-based pan-genome for pepper. We annotated an average of 42 972 gene families in each pepper accession, defining a set of 19 662 core and 23 115 non-core gene families. The new pepper pan-genome includes informative variants, 222 159 PAVs, 12 322 CNVs, and 16 032 inversions. Pan-genome analysis revealed PAVs associated with important agricultural traits, including potyvirus resistance, fruit color, pungency, and pepper fruit orientation. Comparatively, a large number of genes are affected by PAVs, which is positively correlated with the high frequency of transposable elements (TEs), indicating TEs play a key role in shaping the genomic landscape of peppers. The datasets presented herein provide a powerful new genomic resource for genetic analysis and genome-assisted breeding for pepper improvement.

CarrotOmics: a genetics and comparative genomics database for carrot (Daucus carota)

CarrotOmics (https://carrotomics.org/) is a comprehensive database for carrot (Daucus carota L.) breeding and research. CarrotOmics was developed using resources available at the MainLab Bioinformatics core (https://www.bioinfo.wsu.edu/) and is implemented using Tripal with Drupal modules. The database delivers access to download or visualize the carrot reference genome with gene predictions, gene annotations and sequence assembly. Other genomic resources include information for 11 224 genetic markers from 73 linkage maps or genotyping-by-sequencing and descriptions of 371 mapped loci. There are records for 1601 Apiales species (or subspecies) and descriptions of 9408 accessions from 11 germplasm collections representing more than 600 of these species. Additionally, 204 Apiales species have phenotypic information, totaling 28 517 observations from 10 041 biological samples. Resources on CarrotOmics are freely available, search functions are provided to find data of interest and video tutorials are available to describe the search functions and genomic tools. CarrotOmics is a timely resource for the Apiaceae research community and for carrot geneticists developing improved cultivars with novel traits addressing challenges including an expanding acreage in tropical climates, an evolving consumer interested in sustainably grown vegetables and a dynamic environment due to climate change. Data from CarrotOmics can be applied in genomic-assisted selection and genetic research to improve basic research and carrot breeding efficiency.

Database URL

https://carrotomics.org/

Chromosome-scale assembly of the Moringa oleifera Lam. genome uncovers polyploid history and evolution of secondary metabolism pathways through tandem duplication

The African Orphan Crops Consortium (AOCC) selected the highly nutritious, fast growing and drought tolerant tree crop moringa (Moringa oleifera Lam.) as one of the first of 101 plant species to have its genome sequenced and a first draft assembly was published in 2019. Given the extensive uses and culture of moringa, often referred to as the multipurpose tree, we generated a significantly improved new version of the genome based on long-read sequencing into 14 pseudochromosomes equivalent to n = 14 haploid chromosomes. We leveraged this nearly complete version of the moringa genome to investigate main drivers of gene family and genome evolution that may be at the origin of relevant biological innovations including agronomical favorable traits. Our results reveal that moringa has not undergone any additional whole-genome duplication (WGD) or polyploidy event beyond the gamma WGD shared by all core eudicots. Moringa duplicates retained following that ancient gamma events are also enriched for functions commonly considered as dosage balance sensitive. Furthermore, tandem duplications seem to have played a prominent role in the evolution of specific secondary metabolism pathways including those involved in the biosynthesis of bioactive glucosinolate, flavonoid, and alkaloid compounds as well as of defense response pathways and might, at least partially, explain the outstanding phenotypic plasticity attributed to this species. This study provides a genetic roadmap to guide future breeding programs in moringa, especially those aimed at improving secondary metabolism related traits.

Ancestral Chromosomes for Family Peronosporaceae Inferred from a Telomere-to-Telomere Genome Assembly of Peronospora effusa

Downy mildew disease of spinach, caused by the oomycete Peronospora effusa, causes major losses to spinach production. In this study, the 17 chromosomes of P. effusa were assembled telomere-to-telomere, using Pacific Biosciences high-fidelity reads. Of these, 16 chromosomes are complete and gapless; chromosome 15 contains one gap bridging the nucleolus organizer region. This is the first telomere-to-telomere genome assembly for an oomycete. Putative centromeric regions were identified on all chromosomes. This new assembly enables a reevaluation of the genomic composition of Peronospora spp.; the assembly was almost double the size and contained more repeat sequences than previously reported for any Peronospora species. Genome fragments consistently underrepresented in six previously reported assemblies of P. effusa typically encoded repeats. Some genes annotated as encoding effectors were organized into multigene clusters on several chromosomes. Putative effectors were annotated on 16 of the 17 chromosomes. The intergenic distances between annotated genes were consistent with compartmentalization of the genome into gene-dense and gene-sparse regions. Genes encoding putative effectors were enriched in gene-sparse regions. The near-gapless assembly revealed apparent horizontal gene transfer from Ascomycete fungi. Gene order was highly conserved between P. effusa and the genetically oriented assembly of the oomycete Bremia lactucae; high levels of synteny were also detected with Phytophthora sojae. Extensive synteny between phylogenetically distant species suggests that many other oomycete species may have similar chromosome organization. Therefore, this assembly provides the foundation for genomic analyses of diverse oomycetes.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Chromosome evolution and the genetic basis of agronomically important traits in greater yam

The nutrient-rich tubers of the greater yam, Dioscorea alata L., provide food and income security for millions of people around the world. Despite its global importance, however, greater yam remains an orphan crop. Here, we address this resource gap by presenting a highly contiguous chromosome-scale genome assembly of D. alata combined with a dense genetic map derived from African breeding populations. The genome sequence reveals an ancient allotetraploidization in the Dioscorea lineage, followed by extensive genome-wide reorganization. Using the genomic tools, we find quantitative trait loci for resistance to anthracnose, a damaging fungal pathogen of yam, and several tuber quality traits. Genomic analysis of breeding lines reveals both extensive inbreeding as well as regions of extensive heterozygosity that may represent interspecific introgression during domestication. These tools and insights will enable yam breeders to unlock the potential of this staple crop and take full advantage of its adaptability to varied environments.

The spinach YY genome reveals sex chromosome evolution, domestication, and introgression history of the species

Background

Spinach (Spinacia oleracea L.) is a dioecious species with an XY sex chromosome system, but its Y chromosome has not been fully characterized. Our knowledge about the history of its domestication and improvement remains limited.

Results

A high-quality YY genome of spinach is assembled into 952 Mb in six pseudo-chromosomes. By a combination of genetic mapping, Genome-Wide Association Studies, and genomic analysis, we characterize a 17.42-Mb sex determination region (SDR) on chromosome 1. The sex chromosomes of spinach evolved when an insertion containing sex determination genes occurred, followed by a large genomic inversion about 1.98 Mya. A subsequent burst of SDR-specific repeats (0.1–0.15 Mya) explains the large size of this SDR. We identify a Y-specific gene, NRT1/PTR 6.4 which resides in this insertion, as a strong candidate for the sex determination or differentiation factor. Resequencing of 112 spinach genomes reveals a severe domestication bottleneck approximately 10.87 Kya, which dates the domestication of spinach 7000 years earlier than the archeological record. We demonstrate that a strong selection signal associated with internode elongation and leaf area expansion is associated with domestication of edibility traits in spinach. We find that several strong genomic introgressions from the wild species Spinacia turkestanica and Spinacia tetrandra harbor desirable alleles of genes related to downy mildew resistance, frost resistance, leaf morphology, and flowering-time shift, which likely contribute to spinach improvement.

Conclusions

Analysis of the YY genome uncovers evolutionary forces shaping nascent sex chromosome evolution in spinach. Our findings provide novel insights about the domestication and improvement of spinach.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02633-x.

Genomic Resources to Guide Improvement of the Shea Tree

A defining component of agroforestry parklands across Sahelo-Sudanian Africa (SSA), the shea tree (Vitellaria paradoxa) is central to sustaining local livelihoods and the farming environments of rural communities. Despite its economic and cultural value, however, not to mention the ecological roles it plays as a dominant parkland species, shea remains semi-domesticated with virtually no history of systematic genetic improvement. In truth, shea's extended juvenile period makes traditional breeding approaches untenable; but the opportunity for genome-assisted breeding is immense, provided the foundational resources are available. Here we report the development and public release of such resources. Using the FALCON-Phase workflow, 162.6 Gb of long-read PacBio sequence data were assembled into a 658.7 Mbp, chromosome-scale reference genome annotated with 38,505 coding genes. Whole genome duplication (WGD) analysis based on this gene space revealed clear signatures of two ancient WGD events in shea's evolutionary past, one prior to the Astrid-Rosid divergence (116-126 Mya) and the other at the root of the order Ericales (65-90 Mya). In a first genome-wide look at the suite of fatty acid (FA) biosynthesis genes that likely govern stearin content, the primary determinant of shea butter quality, relatively high copy numbers of six key enzymes were found (KASI, KASIII, FATB, FAD2, FAD3, and FAX2), some likely originating in shea's more recent WGD event. To help translate these findings into practical tools for characterization, selection, and genome-wide association studies (GWAS), resequencing data from a shea diversity panel was used to develop a database of more than 3.5 million functionally annotated, physically anchored SNPs. Two smaller, more curated sets of suggested SNPs, one for GWAS (104,211 SNPs) and the other targeting FA biosynthesis genes (90 SNPs), are also presented. With these resources, the hope is to support national programs across the shea belt in the strategic, genome-enabled conservation and long-term improvement of the shea tree for SSA.

A chromosome-level Amaranthus cruentus genome assembly highlights gene family evolution and biosynthetic gene clusters that may underpin the nutritional value of this traditional crop

Traditional crops have historically provided accessible and affordable nutrition to millions of rural dwellers but have been neglected, with most modern agricultural systems over-reliant on a small number of internationally traded crops. Traditional crops are typically well-adapted to local agro-ecological conditions and many are nutrient-dense. They can play a vital role in local food systems through enhanced nutrition (particularly where diets are dominated by starch crops), food security and livelihoods for smallholder farmers, and a climate-resilient and biodiverse agriculture. Using short-read, long-read and phased sequencing technologies, we generated a high-quality chromosome-level genome assembly for Amaranthus cruentus, an under-researched crop with micronutrient- and protein-rich leaves and gluten-free seed, but lacking improved varieties, with respect to productivity and quality traits. The 370.9 Mb genome demonstrates a shared whole genome duplication with a related species, Amaranthus hypochondriacus. Comparative genome analysis indicates chromosomal loss and fusion events following genome duplication that are common to both species, as well as fission of chromosome 2 in A. cruentus alone, giving rise to a haploid chromosome number of 17 (versus 16 in A. hypochondriacus). Genomic features potentially underlying the nutritional value of this crop include two A. cruentus-specific genes with a likely role in phytic acid synthesis (an anti-nutrient), expansion of ion transporter gene families, and identification of biosynthetic gene clusters conserved within the amaranth lineage. The A. cruentus genome assembly will underpin much-needed research and global breeding efforts to develop improved varieties for economically viable cultivation and realization of the benefits to global nutrition security and agrobiodiversity.

An anchored chromosome-scale genome assembly of spinach improves annotation and reveals extensive gene rearrangements in euasterids

Spinach (Spinacia oleracea L.) is a member of the Caryophyllales family, a basal eudicot asterid that consists of sugar beet (Beta vulgaris L. subsp. vulgaris), quinoa (Chenopodium quinoa Willd.), and amaranth (Amaranthus hypochondriacus L.). With the introduction of baby leaf types, spinach has become a staple food in many homes. Production issues focus on yield, nitrogen-use efficiency and resistance to downy mildew (Peronospora effusa). Although genomes are available for the above species, a chromosome-level assembly exists only for quinoa, allowing for proper annotation and structural analyses to enhance crop improvement. We independently assembled and annotated genomes of the cultivar Viroflay using short-read strategy (Illumina) and long-read strategies (Pacific Biosciences) to develop a chromosome-level, genetically anchored assembly for spinach. Scaffold N50 for the Illumina assembly was 389 kb, whereas that for Pacific BioSciences was 4.43 Mb, representing 911 Mb (93% of the genome) in 221 scaffolds, 80% of which are anchored and oriented on a sequence-based genetic map, also described within this work. The two assemblies were 99.5% collinear. Independent annotation of the two assemblies with the same comprehensive transcriptome dataset show that the quality of the assembly directly affects the annotation with significantly more genes predicted (26,862 vs. 34,877) in the long-read assembly. Analysis of resistance genes confirms a bias in resistant gene motifs more typical of monocots. Evolutionary analysis indicates that Spinacia is a paleohexaploid with a whole-genome triplication followed by extensive gene rearrangements identified in this work. Diversity analysis of 75 lines indicate that variation in genes is ample for hypothesis-driven, genomic-assisted breeding enabled by this work.

Genome sequence and genetic diversity analysis of an under-domesticated orphan crop, white fonio (Digitaria exilis)

Background

Digitaria exilis, white fonio, is a minor but vital crop of West Africa that is valued for its resilience in hot, dry, and low-fertility environments and for the exceptional quality of its grain for human nutrition. Its success is hindered, however, by a low degree of plant breeding and improvement.

Findings

We sequenced the fonio genome with long-read SMRT-cell technology, yielding a ∼761 Mb assembly in 3,329 contigs (N50, 1.73 Mb; L50, 126). The assembly approaches a high level of completion, with a BUSCO score of >99%. The fonio genome was found to be a tetraploid, with most of the genome retained as homoeologous duplications that differ overall by ∼4.3%, neglecting indels. The 2 genomes within fonio were found to have begun their independent divergence ∼3.1 million years ago. The repeat content (>49%) is fairly standard for a grass genome of this size, but the ratio of Gypsy to Copia long terminal repeat retrotransposons (∼6.7) was found to be exceptionally high. Several genes related to future improvement of the crop were identified including shattering, plant height, and grain size. Analysis of fonio population genetics, primarily in Mali, indicated that the crop has extensive genetic diversity that is largely partitioned across a north-south gradient coinciding with the Sahel and Sudan grassland domains.

Conclusions

We provide a high-quality assembly, annotation, and diversity analysis for a vital African crop. The availability of this information should empower future research into further domestication and improvement of fonio.

Genetic diversity, SNP-trait associations and genomic selection accuracy in a west African collection of Kersting's groundnut [Macrotyloma geocarpum(Harms) Maréchal & Baudet]

Understanding the mechanisms governing complex traits variation is a requirement for efficient crop improvement. In this study, the molecular characterization, marker-trait associations and the possibility for genomic selection in a collection of 281 Kersting's groundnut accessions were carried out. The diversity panel was phenotyped using an Alpha lattice design with two replicates in two contrasting environments. Accessions were genotyped using genotyping by sequencing technology. Genome-wide association analyses were performed between single nucleotide polymorphism markers and yield-related traits across tested environments. SNP markers were used to calculate the observed (Ho) and expected heterozygosity (He), and the total gene diversity (Ht). Genetic differentiation among accessions across ecological regions of origin was analysed. Our results revealed 493 quality SNPs of which 113 had a minor allele frequency>0.05, a total gene diversity of 0.43 and average Ho and He values of 0.04 and 0.22, respectively. Four clusters, highly differentiated by seed coat colour (Fst = 0.79), were identified. The population structure analysis showed two subpopulations with high differentiation across ecological regions (Fst = 0.37). The GWAS revealed 10 significant marker-trait associations, of which six SNPs were consistent across environments. The genomic selection through cross-validation showed moderate to high prediction accuracies for leaflet length, seed dimension traits, 100 seed weight, days to 50% flowering and days to maturity. This demonstrates the existence of genetic variability within Kersting's groundnut and shows the potential for the improvement of the species. The findings also provide a first insight into the phenotype-to-genotype relationships in Kersting's groundnut, using SNP markers.

Transcriptional analyses of differential cultivars during resistant and susceptible interactions with Peronospora effusa, the causal agent of spinach downy mildew

Downy mildew of spinach is caused by the obligate oomycete pathogen, Peronospora effusa. The disease causes significant economic losses, especially in the organic sector of the industry where the use of synthetic fungicides is not permitted for disease control. New pathotypes of this pathogen are increasingly reported which are capable of breaking resistance. In this study, we took advantage of new spinach genome resources to conduct RNA-seq analyses of transcriptomic changes in leaf tissue of resistant and susceptible spinach cultivars Solomon and Viroflay, respectively, at an early stage of pathogen establishment (48 hours post inoculation, hpi) to a late stage of symptom expression and pathogen sporulation (168 hpi). Fold change differences in gene expression were recorded between the two cultivars to identify candidate genes for resistance. In Solomon, the hypersensitive inducible genes such as pathogenesis-related gene PR-1, glutathione-S-transferase, phospholipid hydroperoxide glutathione peroxidase and peroxidase were significantly up-regulated uniquely at 48 hpi and genes involved in zinc finger CCCH protein, glycosyltransferase, 1-aminocyclopropane-1-carboxylate oxidase homologs, receptor-like protein kinases were expressed at 48 hpi through 168 hpi. The types of genes significantly up-regulated in Solomon in response to the pathogen suggests that salicylic acid and ethylene signaling pathways mediate resistance. Furthermore, many genes involved in the flavonoid and phenylpropanoid pathways were highly expressed in Viroflay compared to Solomon at 168 hpi. As anticipated, an abundance of significantly down-regulated genes was apparent at 168 hpi, reflecting symptom development and sporulation in cultivar Viroflay, but not at 48 hpi. In the pathogen, genes encoding RxLR-type effectors were expressed during early colonization of cultivar Viroflay while crinkler-type effector genes were expressed at the late stage of the colonization. Our results provide insights on gene expression in resistant and susceptible spinach-P. effusa interactions, which can guide future studies to assess candidate genes necessary for downy mildew resistance in spinach.

In Vitro Bioaccessibility of Carotenoids and Chlorophylls in a Diverse Collection of Spinach Accessions and Commercial Cultivars

Spinach, a nutrient-dense, green-leafy vegetable, is a rich source of carotenoid and chlorophyll bioactives. While the content of bioactives is known to vary with the genotype, variation in bioaccessibility is unknown. Bioaccessibility was explored in 71 greenhouse-grown spinach genotypes in fall and spring 2018/2019. Spinach was phenotyped for its greenness, leaf texture, leaf shape, and SPAD chlorophyll content. Postharvest, spinach was washed, blanched, and homogenized prior to assessment of bioactive bioaccessibility using a novel high-throughput in vitro digestion model followed by high-performance liquid chromatography with a photodiode array detector analysis. There was a significant variation in the bioaccessible content for all bioactives (p < 0.05), except for chlorophyll b (p = 0.063) in spring-grown spinach. The correlation coefficients of bioaccessible contents between seasons reveal that lutein (r = 0.52) and β-carotene (r = 0.55) were correlated to a greater extent than chlorophyll a (r = 0.38) and chlorophyll b (r = 0.19). The results suggest that carotenoid and chlorophyll bioaccessible contents may vary based on spinach genotypes and may be stable across seasons.

Breeding Crops for Enhanced Food Safety

An increasing global population demands a continuous supply of nutritious and safe food. Edible products can be contaminated with biological (e.g., bacteria, virus, protozoa), chemical (e.g., heavy metals, mycotoxins), and physical hazards during production, storage, transport, processing, and/or meal preparation. The substantial impact of foodborne disease outbreaks on public health and the economy has led to multidisciplinary research aimed to understand the biology underlying the different contamination processes and how to mitigate food hazards. Here we review the knowledge, opportunities, and challenges of plant breeding as a tool to enhance the food safety of plant-based food products. First, we discuss the significant effect of plant genotypic and phenotypic variation in the contamination of plants by heavy metals, mycotoxin-producing fungi, and human pathogenic bacteria. In addition, we discuss the various factors (i.e., temperature, relative humidity, soil, microbiota, cultural practices, and plant developmental stage) that can influence the interaction between plant genetic diversity and contaminant. This exposes the necessity of a multidisciplinary approach to understand plant genotype × environment × microbe × management interactions. Moreover, we show that the numerous possibilities of crop/hazard combinations make the definition and identification of high-risk pairs, such as Salmonella-tomato and Escherichia coli-lettuce, imperative for breeding programs geared toward improving microbial safety of produce. Finally, we discuss research on developing effective assays and approaches for selecting desirable breeding germplasm. Overall, it is recognized that although breeding programs for some human pathogen/toxin systems are ongoing (e.g., Fusarium in wheat), it would be premature to start breeding when targets and testing systems are not well defined. Nevertheless, current research is paving the way toward this goal and this review highlights advances in the field and critical points for the success of this initiative that were discussed during the Breeding Crops for Enhanced Food Safety workshop held 5-6 June 2019 at University of California, Davis.

African Orphan Crops Consortium (AOCC): status of developing genomic resources for African orphan crops

The African Orphan Crops Consortium (AOCC) successfully initiated the ambitious genome sequencing project of 101 African orphan crops/trees with 6 genomes sequenced, 6 near completion, and 20 currently in progress. Addressing stunting, malnutrition, and hidden hunger through nutritious, economic, and resilient agri-food system is one of the major agricultural challenges of this century. As sub-Saharan Africa harbors a large portion of the severely malnourished population, the African Orphan Crops Consortium (AOCC) was established in 2011 with an aim to reduce stunting and malnutrition by providing nutritional security through improving locally adapted nutritious, but neglected, under-researched or orphan African food crops. Foods from these indigenous or naturalized crops and trees are rich in minerals, vitamins, and antioxidant, and are an integral part of the dietary portfolio and cultural, social, and economic milieu of African farmers. Through stakeholder consultations supported by the African Union, 101 African orphan and under-researched crop species were prioritized to mainstream into African agri-food systems. The AOCC, through a network of international-regional-public-private partnerships and collaborations, is generating genomic resources of three types, i.e., reference genome sequence, transcriptome sequence, and re-sequencing 100 accessions/species, using next-generation sequencing (NGS) technology. Furthermore, the University of California Davis African Plant Breeding Academy under the AOCC banner is training 150 lead African scientists to breed high yielding, nutritious, and climate-resilient (biotic and abiotic stress tolerant) crop varieties that meet African farmer and consumer needs. To date, one or more forms of sequence data have been produced for 60 crops. Reference genome sequences for six species have already been published, 6 are almost near completion, and 19 are in progress.

Association between vitamin content, plant morphology and geographical origin in a worldwide collection of the orphan crop Gynandropsis gynandra (Cleomaceae)

The variability in nutrient content and morphology in Gynandropsis gynandra is associated with the geographic origin of the accessions and provides a basis for breeding for higher levels of vitamin C, carotenoids or tocopherols in higher-yielding cultivars. We examined the variation in carotenoids, tocopherols and ascorbic acid as well as morphological traits in a worldwide germplasm of 76 accessions of the orphan leafy vegetable Gynandropsis gynandra (Cleomaceae) using greenhouse experiments and high-performance liquid chromatography analysis. The levels of carotenoids and tocopherols accumulating in the leaves varied significantly across accessions and were linked with the geographical origin and morphological variation. The main carotenoids included lutein, β-carotene, α-carotene and violaxanthin. A twofold to threefold variation was observed for these compounds. The main tocopherols detected were α-tocopherol and γ-tocopherol with a 20-fold variation. A ninefold variation in vitamin C concentration and independent of geographical origin was observed. Overall, the accessions were grouped into three clusters based on variation in nutrient content and morphology. West African accessions were short plants with small leaves and with high tocopherol contents and relatively low carotenoid contents, Asian accessions were short plants with broad leaves and with relatively low carotenoid and high tocopherol contents, while East-Southern African plants were tall with high contents of both carotenoids and chlorophylls and low tocopherol contents. Carotenoids were positively correlated with plant height as well as foliar and floral traits but negatively correlated with tocopherols. The absence of a significant correlation between vitamin C and other traits indicated that breeding for high carotenoids or tocopherols content may be coupled with improved leaf yield and vitamin C content. Our study provides baseline information on the natural variation available for traits of interest for breeding for enhanced crop yield and nutrient content in Gynandropsis gynandra.

The zinc-finger transcription factor CcLOL1 controls chloroplast development and immature pepper fruit color in Capsicum chinense and its function is conserved in tomato

Chloroplast development and chlorophyll content in the immature fruit has a major impact on the morphology and quality in pepper (Capsicum spp.) fruit. Two major quantitative trait loci (QTLs), pc1 and pc10 that affect chlorophyll content in the pepper fruit by modulation of chloroplast compartment size were previously identified in chromosomes 1 and 10, respectively. The pepper homolog of GOLDEN2-LIKE transcription factor (CaGLK2) has been found as underlying pc10, similar to its effect on tomato chloroplast development. In the present study, we identified the pepper homolog of the zinc-finger transcription factor LOL1 (LSD ONE LIKE1; CcLOL1) as the gene underlying pc1. LOL1 has been identified in Arabidopsis as a positive regulator of programmed cell death and we report here on its role in controlling fruit development in the Solanaceae in a fruit-specific manner. The light-green C. chinense parent used for QTL mapping was found to carry a null mutation in CcLOL1. Verification of the function of the gene was done by generating CRISPR/Cas9 knockout mutants of the orthologous tomato gene resulting in light-green tomato fruits, indicating functional conservation of the orthologous genes in controlling chlorophyll content in the Solanaceae. Transcriptome profiling of light and dark-green bulks differing for pc1, showed that the QTL affects multiple photosynthesis and oxidation-reduction associated genes in the immature green fruit. Allelic diversity of three known genes CcLOL1, CaGLK2, and CcAPRR2 that influence pepper immature fruit color, was found to be associated with variation in chlorophyll content primarily in C. chinense.

Diversity in the Bioaccessibility of Carotenoid and Chlorophyll Compounds in 69 Spinach Genotypes (P06-007-19)

Objectives

Spinach is one of the most nutrient dense and popular green leafy vegetables. Known as a rich source of pro-vitamin A carotenoids and chlorophylls, the variation in content of these bioactives has been associated with differences in genetics, environment and processing. However, beyond processing, factors affecting their bioavailability remain relatively unknown. Establishing the presence of diversity in bioaccessibility, a phenotypical surrogate for bioavailability, of phytochemicals from public germplasm collections and commercial cultivars would provide critical information to breeders seeking to improve the nutritional value of the material used in their programs.

Methods

69 spinach accessions from both the USDA Germplasm Resource Information Network (GRIN) and the Center for Genetic Resources (CGN), the Netherlands, and selected commercial varieties were greenhouse-grown during the Fall 2018 season at Piedmont Research Station, NC, harvested, washed, blanched (100°C, 2 min) and homogenized (30 sec). Carotenoid and chlorophyll bioaccessibility were determined from spinach homogenates formulated with 5% canola oil using a high-throughput, three-phase in vitro digestion followed by centrifugation and filtration to isolate the micellar fraction. Transfer of lutein (LUT), β-carotene (BC) and pheophytin A (Phe A) and B (Phe B) from spinach to micellar fractions was quantified by LC to determine bioaccessibility.

Results

LUT and BC were well recovered through in vitro digestion while chlorophylls were quantitatively converted to corresponding pheophytins. Relative bioaccessibility for all phytochemicals was found to be significantly different across genotypes, as determined by ANOVA analysis. LUT bioaccessibiltiy ranged from 9.4–30.6% (P = 0.01), BC (12.9–42.5%; P < 0.001), PheA (24.5–51%; P = < 0.001), and PheB (19.8–46.9% P = < 0.001). Absolute bioaccessible content ranged from 0.6–8.2 and 0.9–3.7 mg/100 g fw spinach for LUT and BC, respectively.

Conclusions

Results suggest that carotenoid and chlorophyll bioaccessibility may vary based on spinach genotype. The relationship of bioaccessibility to quality traits and potential interactions with agronomic (GxE) and processing (GxP) conditions remain to be explored.

Funding Sources

Foundation for Food and Agriculture Research.

The draft genomes of five agriculturally important African orphan crops

BACKGROUND

The expanding world population is expected to double the worldwide demand for food by 2050. Eighty-eight percent of countries currently face a serious burden of malnutrition, especially in Africa and south and southeast Asia. About 95% of the food energy needs of humans are fulfilled by just 30 species, of which wheat, maize, and rice provide the majority of calories. Therefore, to diversify and stabilize the global food supply, enhance agricultural productivity, and tackle malnutrition, greater use of neglected or underutilized local plants (so-called orphan crops, but also including a few plants of special significance to agriculture, agroforestry, and nutrition) could be a partial solution.

RESULTS

Here, we present draft genome information for five agriculturally, biologically, medicinally, and economically important underutilized plants native to Africa: Vigna subterranea, Lablab purpureus, Faidherbia albida, Sclerocarya birrea, and Moringa oleifera. Assembled genomes range in size from 217 to 654 Mb. In V. subterranea, L. purpureus, F. albida, S. birrea, and M. oleifera, we have predicted 31,707, 20,946, 28,979, 18,937, and 18,451 protein-coding genes, respectively. By further analyzing the expansion and contraction of selected gene families, we have characterized root nodule symbiosis genes, transcription factors, and starch biosynthesis-related genes in these genomes.

CONCLUSIONS

These genome data will be useful to identify and characterize agronomically important genes and understand their modes of action, enabling genomics-based, evolutionary studies, and breeding strategies to design faster, more focused, and predictable crop improvement programs.

Nitrogen fixation in a landrace of maize is supported by a mucilage-associated diazotrophic microbiota

Plants are associated with a complex microbiota that contributes to nutrient acquisition, plant growth, and plant defense. Nitrogen-fixing microbial associations are efficient and well characterized in legumes but are limited in cereals, including maize. We studied an indigenous landrace of maize grown in nitrogen-depleted soils in the Sierra Mixe region of Oaxaca, Mexico. This landrace is characterized by the extensive development of aerial roots that secrete a carbohydrate-rich mucilage. Analysis of the mucilage microbiota indicated that it was enriched in taxa for which many known species are diazotrophic, was enriched for homologs of genes encoding nitrogenase subunits, and harbored active nitrogenase activity as assessed by acetylene reduction and ¹⁵N₂ incorporation assays. Field experiments in Sierra Mixe using ¹⁵N natural abundance or ¹⁵N-enrichment assessments over 5 years indicated that atmospheric nitrogen fixation contributed 29%–82% of the nitrogen nutrition of Sierra Mixe maize.

Nitrogen is an essential nutrient for plants, and for many nonlegume crops, the requirement for nitrogen is primarily met by the use of inorganic fertilizers. These fertilizers are produced from fossil fuel by energy-intensive processes that are estimated to use 1% to 2% of the total global energy supply and produce an equivalent share of greenhouse gases. Because maize (Zea mays L.) is a significant recipient of nitrogen fertilization, a research goal for decades has been to identify or engineer mechanisms for biological fixation of atmospheric nitrogen in association with this crop. We hypothesized that isolated indigenous landraces of maize grown using traditional practices with little or no fertilizer might have evolved strategies to improve plant performance under low-nitrogen nutrient conditions. Here, we show that for one such maize landrace grown in nitrogen-depleted fields near Oaxaca, Mexico, 29%–82% of the plant nitrogen is derived from atmospheric nitrogen. High levels of nitrogen fixation are supported, at least in part, by the abundant production of a sugar-rich mucilage associated with aerial roots that provides a home to a complex nitrogen-fixing microbiome.

Topological Data Analysis as a Morphometric Method: Using Persistent Homology to Demarcate a Leaf Morphospace

Current morphometric methods that comprehensively measure shape cannot compare the disparate leaf shapes found in seed plants and are sensitive to processing artifacts. We explore the use of persistent homology, a topological method applied as a filtration across simplicial complexes (or more simply, a method to measure topological features of spaces across different spatial resolutions), to overcome these limitations. The described method isolates subsets of shape features and measures the spatial relationship of neighboring pixel densities in a shape. We apply the method to the analysis of 182,707 leaves, both published and unpublished, representing 141 plant families collected from 75 sites throughout the world. By measuring leaves from throughout the seed plants using persistent homology, a defined morphospace comparing all leaves is demarcated. Clear differences in shape between major phylogenetic groups are detected and estimates of leaf shape diversity within plant families are made. The approach predicts plant family above chance. The application of a persistent homology method, using topological features, to measure leaf shape allows for a unified morphometric framework to measure plant form, including shapes, textures, patterns, and branching architectures.

Comparative transcriptomics and genomic patterns of discordance in Capsiceae (Solanaceae)

The integration of genomics and phylogenetics allows new insight into the structure of gene tree discordance, the relationships among gene position, gene history, and rate of evolution, as well as the correspondence of gene function, positive selection, and gene ontology enrichment across lineages. We explore these issues using the tribe Capsiceae (Solanaceae), which is comprised of the genera Lycianthes and Capsicum (peppers). In combining the annotated genomes of Capsicum with newly sequenced transcriptomes of four species of Lycianthes and Capsicum, we develop phylogenies for 6747 genes, and construct a backbone species tree using both concordance and explicit phylogenetic network approaches. We quantify phylogenetic discordance among individual gene trees, measure their rates of synonymous and nonsynonymous substitution, and test whether they were positively selected along any branch of the phylogeny. We then map these genes onto the annotated Capsicum genome and test whether rates of evolution, gene history, and gene ontology vary significantly with gene position. We observed substantial discordance among gene trees. A bifurcating species tree placing Capsicum within a paraphyletic Lycianthes was supported over all phylogenetic networks. Rates of synonymous and nonsynonymous substitution varied 41-fold and 130-fold among genes, respectively, and were significantly lower in pericentromeric regions. We found that results of concordance tree analyses vary depending on the subset of genes used, and that genes within the pericentromeric regions only capture a portion of the observed discordance. We identified 787 genes that have been positively selected throughout the diversification history of Capsiceae, and discuss the importance of these genes as targets for investigation of economically important traits in the domesticated peppers.

Quantitative Trait Loci Controlling Fruit Size and Other Horticultural Traits in Bell Pepper (Capsicum annuum)

Bell pepper ( L.) is a group of fruit vegetables that has large variation in fruit shape, fruit size, and horticultural traits. Using unadapted sources of germplasm to bring in novel alleles while maintaining favorable quality and horticultural traits is challenging for breeding in pepper. A genetic map with 318 loci from genotype-by-sequencing (GBS) and single nucleotide polymorphism assays was generated from a recombinant inbred line population derived from a cultivated bell-type 'Maor' and a landrace highly resistant to , 'Criollo de Morelos-334'. Forty-nine quantitative trait loci (QTLs) were detected for fruit, leaf, and horticultural traits with the scantwo permutation and stepwiseqtl methods from R/qtl. With the availability of a pepper reference genome and GBS data, candidate genes for pepper organ size and other horticultural traits were predicted. , , and genes were candidate genes for controlling organ sizes on chromosome 1, 2, and 3, respectively. Two candidate genes controlling trichome formation in pepper are located at chromosome 10: and . The locus on chromosome 10, which encodes a member of the R2R3 MYB-domain family of proteins, has a function in anthocyanin accumulation. These QTL results and the candidate genes for each trait emphasize the genetic basis of the important traits for breeding with unadapted parents in bell pepper.

A roadmap for breeding orphan leafy vegetable species: a case study of Gynandropsis gynandra (Cleomaceae)

Despite an increasing awareness of the potential of "orphan" or unimproved crops to contribute to food security and enhanced livelihoods for farmers, coordinated research agendas to facilitate production and use of orphan crops by local communities are generally lacking. We provide an overview of the current knowledge on leafy vegetables with a focus on Gynandropsis gynandra, a highly nutritious species used in Africa and Asia, and highlight general and species-specific guidelines for participatory, genomics-assisted breeding of orphan crops. Key steps in genome-enabled orphan leafy vegetables improvement are identified and discussed in the context of Gynandropsis gynandra breeding, including: (1) germplasm collection and management; (2) product target definition and refinement; (3) characterization of the genetic control of key traits; (4) design of the 'process' for cultivar development; (5) integration of genomic data to optimize that 'process'; (6) multi-environmental participatory testing and end-user evaluation; and (7) crop value chain development. The review discusses each step in detail, with emphasis on improving leaf yield, phytonutrient content, organoleptic quality, resistance to biotic and abiotic stresses and post-harvest management.

Reference quality assembly of the 3.5-Gb genome of Capsicum annuum from a single linked-read library

Linked-Read sequencing technology has recently been employed successfully for de novo assembly of human genomes, however, the utility of this technology for complex plant genomes is unproven. We evaluated the technology for this purpose by sequencing the 3.5-gigabase (Gb) diploid pepper (Capsicum annuum) genome with a single Linked-Read library. Plant genomes, including pepper, are characterized by long, highly similar repetitive sequences. Accordingly, significant effort is used to ensure that the sequenced plant is highly homozygous and the resulting assembly is a haploid consensus. With a phased assembly approach, we targeted a heterozygous F₁ derived from a wide cross to assess the ability to derive both haplotypes and characterize a pungency gene with a large insertion/deletion. The Supernova software generated a highly ordered, more contiguous sequence assembly than all currently available C. annuum reference genomes. Over 83% of the final assembly was anchored and oriented using four publicly available de novo linkage maps. A comparison of the annotation of conserved eukaryotic genes indicated the completeness of assembly. The validity of the phased assembly is further demonstrated with the complete recovery of both 2.5-Kb insertion/deletion haplotypes of the PUN1 locus in the F₁ sample that represents pungent and nonpungent peppers, as well as nearly full recovery of the BUSCO2 gene set within each of the two haplotypes. The most contiguous pepper genome assembly to date has been generated which demonstrates that Linked-Read library technology provides a tool to de novo assemble complex highly repetitive heterozygous plant genomes. This technology can provide an opportunity to cost-effectively develop high-quality genome assemblies for other complex plants and compare structural and gene differences through accurate haplotype reconstruction.

Regions Underlying Population Structure and the Genomics of Organ Size Determination in Capsicum annuum

Fruits, as an important part of the human diet, have been under strong selection during domestication. In general, continued directed selection has led to varieties having larger fruit with greater shape variation and tremendous increases in fruit mass. Common cultivated peppers ( L.) are found in a wide range of sizes and shapes. Analysis of genetic relatedness and population structure has shown that the large-fruited, nonpungent types have reduced diversity and comprise a highly structured group. To explore this population structure, a statistical method for detecting fixation within subpopulations was applied to a set of 21 pungent and 19 nonpungent lines that represent the pepper breeding germplasm. We have identified 17 blocks within the pepper genome that are conserved among nonpungent large-fruited varieties. To determine if these regions were fixed by selection on fruit size or pungency, quantitative trait loci (QTLs) from seven studies along with capsaicin biosynthesis genes and homologs of organ size regulatory genes were mapped onto the current pepper genome assembly. Of the 17 fixed regions, 14 overlapped with fruit size or shape QTLs. There were seven putative organ size regulators and seven capsaicin biosynthetic genes within these regions. This work defines genomic regions that underly structure within the nonpungent pepper germplasm and QTLs or genes that may have been selected for during the development of large-fruited nonpungent pepper varieties.

Early fruiting in Synsepalum dulcificum (Schumach. & Thonn.) Daniell juveniles induced by water and inorganic nutrient management

Background. The miracle plant, Synsepalum dulcificum (Schumach. & Thonn.) Daniell is a native African orphan crop species that has recently received increased attention due to its promise as a sweetener and source of antioxidants in both the food and pharmaceutical industries. However, a major obstacle to the species' widespread utilization is its relatively slow growth rate and prolonged juvenile period. Method. In this study, we tested twelve treatments made up of various watering regimes and exogenous nutrient application (nitrogen, phosphorus and potassium, at varying dosages) on the relative survival, growth, and reproductive development of 15-months-old S. dulcificum juveniles. Results. While the plants survived under most tested growing conditions, nitrogen application at doses higher than 1.5 g [seedling] ^-1 was found to be highly detrimental, reducing survival to 0%. The treatment was found to affect all growth traits, and juveniles that received a combination of nitrogen, phosphorus, and potassium (each at a rate of 1.5 g [seedling] ^-1), in addition to daily watering, exhibited the most vegetative growth. The simple daily provision of adequate water was found to greatly accelerate the transition to reproductive maturity in the species (from >36 months to an average of 23 months), whereas nutrient application affected the length of the reproductive phase within a season, as well as the fruiting intensity. Conclusions. This study highlights the beneficial effect of water supply and fertilization on both vegetative and reproductive growth in S. dulcificum. Water supply appeared to be the most important factor unlocking flowering in the species, while the combination of nitrogen, phosphorus and potassium at the dose of 1.5 g (for all) consistently exhibited the highest performance for all growth and yield traits. These findings will help intensify S. dulcificum's breeding and horticultural development.

Diversity analysis of cotton (Gossypium hirsutum L.) germplasm using the CottonSNP63K Array

BACKGROUND

Cotton germplasm resources contain beneficial alleles that can be exploited to develop germplasm adapted to emerging environmental and climate conditions. Accessions and lines have traditionally been characterized based on phenotypes, but phenotypic profiles are limited by the cost, time, and space required to make visual observations and measurements. With advances in molecular genetic methods, genotypic profiles are increasingly able to identify differences among accessions due to the larger number of genetic markers that can be measured. A combination of both methods would greatly enhance our ability to characterize germplasm resources. Recent efforts have culminated in the identification of sufficient SNP markers to establish high-throughput genotyping systems, such as the CottonSNP63K array, which enables a researcher to efficiently analyze large numbers of SNP markers and obtain highly repeatable results. In the current investigation, we have utilized the SNP array for analyzing genetic diversity primarily among cotton cultivars, making comparisons to SSR-based phylogenetic analyses, and identifying loci associated with seed nutritional traits.

RESULTS

The SNP markers distinctly separated G. hirsutum from other Gossypium species and distinguished the wild from cultivated types of G. hirsutum. The markers also efficiently discerned differences among cultivars, which was the primary goal when designing the CottonSNP63K array. Population structure within the genus compared favorably with previous results obtained using SSR markers, and an association study identified loci linked to factors that affect cottonseed protein content.

CONCLUSIONS

Our results provide a large genome-wide variation data set for primarily cultivated cotton. Thousands of SNPs in representative cotton genotypes provide an opportunity to finely discriminate among cultivated cotton from around the world. The SNPs will be relevant as dense markers of genome variation for association mapping approaches aimed at correlating molecular polymorphisms with variation in phenotypic traits, as well as for molecular breeding approaches in cotton.

Use of microsatellite markers for the assessment of bambara groundnut breeding system and varietal purity before genome sequencing

Maximizing the research output from a limited investment is often the major challenge for minor and underutilized crops. However, such crops may be tolerant to biotic and abiotic stresses and are adapted to local, marginal, and low-input environments. Their development through breeding will provide an important resource for future agricultural system resilience and diversification in the context of changing climates and the need to achieve food security. The African Orphan Crops Consortium recognizes the values of genomic resources in facilitating the improvement of such crops. Prior to beginning genome sequencing there is a need for an assessment of line varietal purity and to estimate any residual heterozygosity. Here we present an example from bambara groundnut (Vigna subterranea (L.) Verdc.), an underutilized drought tolerant African legume. Two released varieties from Zimbabwe, identified as potential genotypes for whole genome sequencing (WGS), were genotyped with 20 species-specific SSR markers. The results indicate that the cultivars are actually a mix of related inbred genotypes, and the analysis allowed a strategy of single plant selection to be used to generate non-heterogeneous DNA for WGS. The markers also confirmed very low levels of heterozygosity within individual plants. The application of a pre-screen using co-dominant microsatellite markers is expected to substantially improve the genome assembly, compared to a cultivar bulking approach that could have been adopted.

DNA Sequence Evolution and Rare Homoeologous Conversion in Tetraploid Cotton

Allotetraploid cotton species are a vital source of spinnable fiber for textiles. The polyploid nature of the cotton genome raises many evolutionary questions as to the relationships between duplicated genomes. We describe the evolution of the cotton genome (SNPs and structural variants) with the greatly improved resolution of 34 deeply re-sequenced genomes. We also explore the evolution of homoeologous regions in the AT- and DT-genomes and especially the phenomenon of conversion between genomes. We did not find any compelling evidence for homoeologous conversion between genomes. These findings are very different from other recent reports of frequent conversion events between genomes. We also identified several distinct regions of the genome that have been introgressed between G. hirsutum and G. barbadense, which presumably resulted from breeding efforts targeting associated beneficial alleles. Finally, the genotypic data resulting from this study provides access to a wealth of diversity sorely needed in the narrow germplasm of cotton cultivars.

A HapMap leads to a Capsicum annuum SNP infinium array: a new tool for pepper breeding

The Capsicum genus (Pepper) is a part of the Solanacae family. It has been important in many cultures worldwide for its key nutritional components and uses as spices, medicines, ornamentals and vegetables. Worldwide population growth is associated with demand for more nutritionally valuable vegetables while contending with decreasing resources and available land. These conditions require increased efficiency in pepper breeding to deal with these imminent challenges. Through resequencing of inbred lines we have completed a valuable haplotype map (HapMap) for the pepper genome based on single-nucleotide polymorphisms (SNP). The identified SNPs were annotated and classified based on their gene annotation in the pepper draft genome sequence and phenotype of the sequenced inbred lines. A selection of one marker per gene model was utilized to create the PepperSNP16K array, which simultaneously genotyped 16 405 SNPs, of which 90.7% were found to be informative. A set of 84 inbred and hybrid lines and a mapping population of 90 interspecific F2 individuals were utilized to validate the array. Diversity analysis of the inbred lines shows a distinct separation of bell versus chile/hot pepper types and separates them into five distinct germplasm groups. The interspecific population created between Tabasco (C. frutescens chile type) and P4 (C. annuum blocky type) produced a linkage map with 5546 markers separated into 1361 bins on twelve 12 linkage groups representing 1392.3 cM. This publically available genotyping platform can be used to rapidly assess a large number of markers in a reproducible high-throughput manner for pepper. As a standardized tool for genetic analyses, the PepperSNP16K can be used worldwide to share findings and analyze QTLs for important traits leading to continued improvement of pepper for consumers. Data and information on the array are available through the Solanaceae Genomics Network.

New Insights on Eggplant/Tomato/Pepper Synteny and Identification of Eggplant and Pepper Orthologous QTL

Eggplant, pepper, and tomato are the most exploited berry-producing vegetables within the Solanaceae family. Their genomes differ in size, but each has 12 chromosomes which have undergone rearrangements causing a redistribution of loci. The genome sequences of all three species are available but differ in coverage, assembly quality and percentage of anchorage. Determining their syntenic relationship and QTL orthology will contribute to exploit genomic resources and genetic data for key agronomic traits. The syntenic analysis between tomato and pepper based on the alignment of 34,727 tomato CDS to the pepper genome sequence, identified 19,734 unique hits. The resulting synteny map confirmed the 14 inversions and 10 translocations previously documented, but also highlighted 3 new translocations and 4 major new inversions. Furthermore, each of the 12 chromosomes exhibited a number of rearrangements involving small regions of 0.5-0.7 Mbp. Due to high fragmentation of the publicly available eggplant genome sequence, physical localization of most eggplant QTL was not possible, thus, we compared the organization of the eggplant genetic map with the genome sequence of both tomato and pepper. The eggplant/tomato syntenic map confirmed all the 10 translocations but only 9 of the 14 known inversions; on the other hand, a newly detected inversion was recognized while another one was not confirmed. The eggplant/pepper syntenic map confirmed 10 translocations and 8 inversions already detected and suggested a putative new translocation. In order to perform the assessment of eggplant and pepper QTL orthology, the eggplant and pepper sequence-based markers located in their respective genetic map were aligned onto the pepper genome. GBrowse in pepper was used as reference platform for QTL positioning. A set of 151 pepper QTL were located as well as 212 eggplant QTL, including 76 major QTL (PVE ≥ 10%) affecting key agronomic traits. Most were confirmed to cluster in orthologous chromosomal regions. Our results highlight that the availability of genome sequences for an increasing number of crop species and the development of "ultra-dense" physical maps provide new and key tools for detailed syntenic and orthology studies between related plant species.

A Long-Read Transcriptome Assembly of Cotton (Gossypium hirsutum L.) and Intraspecific Single Nucleotide Polymorphism Discovery

Upland cotton (Gossypium hirsutum L.) has a narrow germplasm base, which constrains marker development and hampers intraspecific breeding. A pressing need exists for high-throughput single nucleotide polymorphism (SNP) markers that can be readily applied to germplasm in breeding and breeding-related research programs. Despite progress made in developing new sequencing technologies during the past decade, the cost of sequencing remains substantial when one is dealing with numerous samples and large genomes. Several strategies have been proposed to lower the cost of sequencing for multiple genotypes of large-genome species like cotton, such as transcriptome sequencing and reduced-representation DNA sequencing. This paper reports the development of a transcriptome assembly of the inbred line Texas Marker-1 (TM-1), a genetic standard for cotton, its usefulness as a reference for RNA sequencing (RNA-seq)-based SNP identification, and the availability of transcriptome sequences of four other cotton cultivars. An assembly of TM-1 was made using Roche 454 transcriptome reads combined with an assembly of all available public expressed sequence tag (EST) sequences of TM-1. The TM-1 assembly consists of 72,450 contigs with a total of 70 million bp. Functional predictions of the transcripts were estimated by alignment to selected protein databases. Transcriptome sequences of the five lines, including TM-1, were obtained using an Illumina Genome Analyzer-II, and the short reads were mapped to the TM-1 assembly to discover SNPs among the five lines. We identified >14,000 unfiltered allelic SNPs, of which ∼3,700 SNPs were retained for assay development after applying several rigorous filters. This paper reports availability of the reference transcriptome assembly and shows its utility in developing intraspecific SNP markers in upland cotton.

Development of a 63K SNP Array for Cotton and High-Density Mapping of Intraspecific and Interspecific Populations of Gossypium spp

High-throughput genotyping arrays provide a standardized resource for plant breeding communities that are useful for a breadth of applications including high-density genetic mapping, genome-wide association studies (GWAS), genomic selection (GS), complex trait dissection, and studying patterns of genomic diversity among cultivars and wild accessions. We have developed the CottonSNP63K, an Illumina Infinium array containing assays for 45,104 putative intraspecific single nucleotide polymorphism (SNP) markers for use within the cultivated cotton species Gossypium hirsutum L. and 17,954 putative interspecific SNP markers for use with crosses of other cotton species with G. hirsutum. The SNPs on the array were developed from 13 different discovery sets that represent a diverse range of G. hirsutum germplasm and five other species: G. barbadense L., G. tomentosum Nuttal × Seemann, G. mustelinum Miers × Watt, G. armourianum Kearny, and G. longicalyx J.B. Hutchinson and Lee. The array was validated with 1,156 samples to generate cluster positions to facilitate automated analysis of 38,822 polymorphic markers. Two high-density genetic maps containing a total of 22,829 SNPs were generated for two F2 mapping populations, one intraspecific and one interspecific, and 3,533 SNP markers were co-occurring in both maps. The produced intraspecific genetic map is the first saturated map that associates into 26 linkage groups corresponding to the number of cotton chromosomes for a cross between two G. hirsutum lines. The linkage maps were shown to have high levels of collinearity to the JGI G. raimondii Ulbrich reference genome sequence. The CottonSNP63K array, cluster file and associated marker sequences constitute a major new resource for the global cotton research community.

BAC-End Sequence-Based SNP Mining in Allotetraploid Cotton (Gossypium) Utilizing Resequencing Data, Phylogenetic Inferences, and Perspectives for Genetic Mapping

A bacterial artificial chromosome library and BAC-end sequences for cultivated cotton (Gossypium hirsutum L.) have recently been developed. This report presents genome-wide single nucleotide polymorphism (SNP) mining utilizing resequencing data with BAC-end sequences as a reference by alignment of 12 G. hirsutum L. lines, one G. barbadense L. line, and one G. longicalyx Hutch and Lee line. A total of 132,262 intraspecific SNPs have been developed for G. hirsutum, whereas 223,138 and 470,631 interspecific SNPs have been developed for G. barbadense and G. longicalyx, respectively. Using a set of interspecific SNPs, 11 randomly selected and 77 SNPs that are putatively associated with the homeologous chromosome pair 12 and 26, we mapped 77 SNPs into two linkage groups representing these chromosomes, spanning a total of 236.2 cM in an interspecific F2 population (G. barbadense 3-79 × G. hirsutum TM-1). The mapping results validated the approach for reliably producing large numbers of both intraspecific and interspecific SNPs aligned to BAC-ends. This will allow for future construction of high-density integrated physical and genetic maps for cotton and other complex polyploid genomes. The methods developed will allow for future Gossypium resequencing data to be automatically genotyped for identified SNPs along the BAC-end sequence reference for anchoring sequence assemblies and comparative studies.

Development and bin mapping of gene-associated interspecific SNPs for cotton (Gossypium hirsutum L.) introgression breeding efforts

BACKGROUND

Cotton (Gossypium spp.) is the largest producer of natural fibers for textile and is an important crop worldwide. Crop production is comprised primarily of G. hirsutum L., an allotetraploid. However, elite cultivars express very small amounts of variation due to the species monophyletic origin, domestication and further bottlenecks due to selection. Conversely, wild cotton species harbor extensive genetic diversity of prospective utility to improve many beneficial agronomic traits, fiber characteristics, and resistance to disease and drought. Introgression of traits from wild species can provide a natural way to incorporate advantageous traits through breeding to generate higher-producing cotton cultivars and more sustainable production systems. Interspecific introgression efforts by conventional methods are very time-consuming and costly, but can be expedited using marker-assisted selection.

RESULTS

Using transcriptome sequencing we have developed the first gene-associated single nucleotide polymorphism (SNP) markers for wild cotton species G. tomentosum, G. mustelinum, G. armourianum and G. longicalyx. Markers were also developed for a secondary cultivated species G. barbadense cv. 3-79. A total of 62,832 non-redundant SNP markers were developed from the five wild species which can be utilized for interspecific germplasm introgression into cultivated G. hirsutum and are directly associated with genes. Over 500 of the G. barbadense markers have been validated by whole-genome radiation hybrid mapping. Overall 1,060 SNPs from the five different species have been screened and shown to produce acceptable genotyping assays.

CONCLUSIONS

This large set of 62,832 SNPs relative to cultivated G. hirsutum will allow for the first high-density mapping of genes from five wild species that affect traits of interest, including beneficial agronomic and fiber characteristics. Upon mapping, the markers can be utilized for marker-assisted introgression of new germplasm into cultivated cotton and in subsequent breeding of agronomically adapted types, including cultivar development.

Development and bin mapping of gene-associated interspecific SNPs for cotton (Gossypium hirsutum L.) introgression breeding efforts

Background

Cotton (Gossypium spp.) is the largest producer of natural fibers for textile and is an important crop worldwide. Crop production is comprised primarily of G. hirsutum L., an allotetraploid. However, elite cultivars express very small amounts of variation due to the species monophyletic origin, domestication and further bottlenecks due to selection. Conversely, wild cotton species harbor extensive genetic diversity of prospective utility to improve many beneficial agronomic traits, fiber characteristics, and resistance to disease and drought. Introgression of traits from wild species can provide a natural way to incorporate advantageous traits through breeding to generate higher-producing cotton cultivars and more sustainable production systems. Interspecific introgression efforts by conventional methods are very time-consuming and costly, but can be expedited using marker-assisted selection.

Results

Using transcriptome sequencing we have developed the first gene-associated single nucleotide polymorphism (SNP) markers for wild cotton species G. tomentosum, G. mustelinum, G. armourianum and G. longicalyx. Markers were also developed for a secondary cultivated species G. barbadense cv. 3–79. A total of 62,832 non-redundant SNP markers were developed from the five wild species which can be utilized for interspecific germplasm introgression into cultivated G. hirsutum and are directly associated with genes. Over 500 of the G. barbadense markers have been validated by whole-genome radiation hybrid mapping. Overall 1,060 SNPs from the five different species have been screened and shown to produce acceptable genotyping assays.

Conclusions

This large set of 62,832 SNPs relative to cultivated G. hirsutum will allow for the first high-density mapping of genes from five wild species that affect traits of interest, including beneficial agronomic and fiber characteristics. Upon mapping, the markers can be utilized for marker-assisted introgression of new germplasm into cultivated cotton and in subsequent breeding of agronomically adapted types, including cultivar development.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-945) contains supplementary material, which is available to authorized users.

CaGLK2 regulates natural variation of chlorophyll content and fruit color in pepper fruit

We provide multiple evidences that CaGLK2 underlies a quantitative trait locus controlling natural variation in chlorophyll content and immature fruit color of pepper via modulating chloroplast compartment size. Pepper fruit quality is attributed to a variety of traits, affecting visual appearance, flavor, chemical composition and nutritional value. Among the quality traits, fruit color is of primary importance because the pigments that confer color are associated with nutrition, health and flavor. Although gene models have been proposed for qualitative aspects of fruit color, large natural variation in quantitative pigment content and fruit color exists in pepper. However, its genetic basis is largely unknown which hampers its utilization for plant improvement. We studied the role of GLK2, a GOLDEN2-like transcription factor that regulates chloroplast development in controlling natural variation for chlorophyll content and immature fruit color of pepper. The role of GLK2 in regulating fruit development has been studied previously in tomato using ectopic expression and the uniform ripening mutant analyses. However, pepper provides a unique opportunity to further study the function of this gene because of the wide natural variation of fruit colors in this species. Segregation, sequencing and expression analyses indicated that pepper GLK2 (CaGLK2) corresponds to the recently reported pc10 QTL that controls chloroplast development and chlorophyll content in pepper. CaGLK2 exerts its effect on chloroplast compartment size predominantly during immature fruit development. We show that the genetic background, sequence variation and expression pattern confer a complex and multi-level regulation of CaGLK2 and fruit color in Capsicum. The positive effect on fruit quality predominantly at the green stage conferred by CaGLK2 can be utilized to breed green pepper varieties with improved nutritional values and taste.

Selection and validation of reference genes for gene expression analysis in switchgrass (Panicum virgatum) using quantitative real-time RT-PCR

Switchgrass (Panicum virgatum) has received a lot of attention as a forage and bioenergy crop during the past few years. Gene expression studies are in progress to improve new traits and develop new cultivars. Quantitative real time PCR (qRT-PCR) has emerged as an important technique to study gene expression analysis. For accurate and reliable results, normalization of data with reference genes is essential. In this work, we evaluate the stability of expression of genes to use as reference for qRT-PCR in the grass P. virgatum. Eleven candidate reference genes, including eEF-1α, UBQ6, ACT12, TUB6, eIF-4a, GAPDH, SAMDC, TUA6, CYP5, U2AF, and FTSH4, were validated for qRT-PCR normalization in different plant tissues and under different stress conditions. The expression stability of these genes was verified by the use of two distinct algorithms, geNorm and NormFinder. Differences were observed after comparison of the ranking of the candidate reference genes identified by both programs but eEF-1α, eIF-4a, CYP5 and U2AF are ranked as the most stable genes in the samples sets under study. Both programs discard the use of SAMDC and TUA6 for normalization. Validation of the reference genes proposed by geNorm and NormFinder were performed by normalization of transcript abundance of a group of target genes in different samples. Results show similar expression patterns when the best reference genes selected by both programs were used but differences were detected in the transcript abundance of the target genes. Based on the above research, we recommend the use of different statistical algorithms to identify the best reference genes for expression data normalization. The best genes selected in this study will help to improve the quality of gene expression data in a wide variety of samples in switchgrass.

Genetic structure and domestication of carrot (Daucus carota subsp. sativus) (Apiaceae)

PREMISE OF THE STUDY

Analyses of genetic structure and phylogenetic relationships illuminate the origin and domestication of modern crops. Despite being an important worldwide vegetable, the genetic structure and domestication of carrot (Daucus carota) is poorly understood. We provide the first such study using a large data set of molecular markers and accessions that are widely dispersed around the world. •

METHODS

Sequencing data from the carrot transcriptome were used to develop 4000 single nucleotide polymorphisms (SNPs). Eighty-four genotypes, including a geographically well-distributed subset of wild and cultivated carrots, were genotyped using the KASPar assay. •

KEY RESULTS

Analysis of allelic diversity of SNP data revealed no reduction of genetic diversity in cultivated vs. wild accessions. Structure and phylogenetic analysis indicated a clear separation between wild and cultivated accessions as well as between eastern and western cultivated carrot. Among the wild carrots, those from Central Asia were genetically most similar to cultivated accessions. Furthermore, we found that wild carrots from North America were most closely related to European wild accessions. •

CONCLUSIONS

Comparing the genetic diversity of wild and cultivated accessions suggested the absence of a genetic bottleneck during carrot domestication. In conjunction with historical documents, our results suggest an origin of domesticated carrot in Central Asia. Wild carrots from North America were likely introduced as weeds with European colonization. These results provide answers to long-debated questions of carrot evolution and domestication and inform germplasm curators and breeders on genetic substructure of carrot genetic resources.

An Ultra-High-Density, Transcript-Based, Genetic Map of Lettuce

We have generated an ultra-high-density genetic map for lettuce, an economically important member of the Compositae, consisting of 12,842 unigenes (13,943 markers) mapped in 3696 genetic bins distributed over nine chromosomal linkage groups. Genomic DNA was hybridized to a custom Affymetrix oligonucleotide array containing 6.4 million features representing 35,628 unigenes of Lactuca spp. Segregation of single-position polymorphisms was analyzed using 213 F7:8 recombinant inbred lines that had been generated by crossing cultivated Lactuca sativa cv. Salinas and L. serriola acc. US96UC23, the wild progenitor species of L. sativa The high level of replication of each allele in the recombinant inbred lines was exploited to identify single-position polymorphisms that were assigned to parental haplotypes. Marker information has been made available using GBrowse to facilitate access to the map. This map has been anchored to the previously published integrated map of lettuce providing candidate genes for multiple phenotypes. The high density of markers achieved in this ultradense map allowed syntenic studies between lettuce and Vitis vinifera as well as other plant species.

Characterization of Capsicum annuum genetic diversity and population structure based on parallel polymorphism discovery with a 30K unigene Pepper GeneChip

The widely cultivated pepper, Capsicum spp., important as a vegetable and spice crop world-wide, is one of the most diverse crops. To enhance breeding programs, a detailed characterization of Capsicum diversity including morphological, geographical and molecular data is required. Currently, molecular data characterizing Capsicum genetic diversity is limited. The development and application of high-throughput genome-wide markers in Capsicum will facilitate more detailed molecular characterization of germplasm collections, genetic relationships, and the generation of ultra-high density maps. We have developed the Pepper GeneChip® array from Affymetrix for polymorphism detection and expression analysis in Capsicum. Probes on the array were designed from 30,815 unigenes assembled from expressed sequence tags (ESTs). Our array design provides a maximum redundancy of 13 probes per base pair position allowing integration of multiple hybridization values per position to detect single position polymorphism (SPP). Hybridization of genomic DNA from 40 diverse C. annuum lines, used in breeding and research programs, and a representative from three additional cultivated species (C. frutescens, C. chinense and C. pubescens) detected 33,401 SPP markers within 13,323 unigenes. Among the C. annuum lines, 6,426 SPPs covering 3,818 unigenes were identified. An estimated three-fold reduction in diversity was detected in non-pungent compared with pungent lines, however, we were able to detect 251 highly informative markers across these C. annuum lines. In addition, an 8.7 cM region without polymorphism was detected around Pun1 in non-pungent C. annuum. An analysis of genetic relatedness and diversity using the software Structure revealed clustering of the germplasm which was confirmed with statistical support by principle components analysis (PCA) and phylogenetic analysis. This research demonstrates the effectiveness of parallel high-throughput discovery and application of genome-wide transcript-based markers to assess genetic and genomic features among Capsicum annuum.

Identification of QTLs for capsaicinoids, fruit quality, and plant architecture-related traits in an interspecific Capsicum RIL population

Quantitative trait loci (QTL) analyses in pepper are common for horticultural, disease resistance, and fruit quality traits; although none of the studies to date have used sequence-based markers associated with genes. In this study we measured plant architectural, phenological, and fruit quality traits in a pepper mapping population consisting of 92 recombinant inbred lines derived from a cross between Capsicum frutescens acc. 2814-6 and C. annuum var. NuMexRNAKY. Phenotypic measurements were correlated to loci in a high-density EST-based genetic map. In total, 96 QTL were identified for 38 traits, including 12 QTL associated with capsaicinoid levels. Twenty-one loci showed correlation among seemingly unrelated phenotypic categories, highlighting tight linkage or shared genetics between previously unassociated traits in pepper.

De novo assembly of the pepper transcriptome (Capsicum annuum): a benchmark for in silico discovery of SNPs, SSRs and candidate genes

Background

Molecular breeding of pepper (Capsicum spp.) can be accelerated by developing DNA markers associated with transcriptomes in breeding germplasm. Before the advent of next generation sequencing (NGS) technologies, the majority of sequencing data were generated by the Sanger sequencing method. By leveraging Sanger EST data, we have generated a wealth of genetic information for pepper including thousands of SNPs and Single Position Polymorphic (SPP) markers. To complement and enhance these resources, we applied NGS to three pepper genotypes: Maor, Early Jalapeño and Criollo de Morelos-334 (CM334) to identify SNPs and SSRs in the assembly of these three genotypes.

Results

Two pepper transcriptome assemblies were developed with different purposes. The first reference sequence, assembled by CAP3 software, comprises 31,196 contigs from >125,000 Sanger-EST sequences that were mainly derived from a Korean F₁-hybrid line, Bukang. Overlapping probes were designed for 30,815 unigenes to construct a pepper Affymetrix GeneChip® microarray for whole genome analyses. In addition, custom Python scripts were used to identify 4,236 SNPs in contigs of the assembly. A total of 2,489 simple sequence repeats (SSRs) were identified from the assembly, and primers were designed for the SSRs. Annotation of contigs using Blast2GO software resulted in information for 60% of the unigenes in the assembly. The second transcriptome assembly was constructed from more than 200 million Illumina Genome Analyzer II reads (80–120 nt) using a combination of Velvet, CLC workbench and CAP3 software packages. BWA, SAMtools and in-house Perl scripts were used to identify SNPs among three pepper genotypes. The SNPs were filtered to be at least 50 bp from any intron-exon junctions as well as flanking SNPs. More than 22,000 high-quality putative SNPs were identified. Using the MISA software, 10,398 SSR markers were also identified within the Illumina transcriptome assembly and primers were designed for the identified markers. The assembly was annotated by Blast2GO and 14,740 (12%) of annotated contigs were associated with functional proteins.

Conclusions

Before availability of pepper genome sequence, assembling transcriptomes of this economically important crop was required to generate thousands of high-quality molecular markers that could be used in breeding programs. In order to have a better understanding of the assembled sequences and to identify candidate genes underlying QTLs, we annotated the contigs of Sanger-EST and Illumina transcriptome assemblies. These and other information have been curated in a database that we have dedicated for pepper project.

A molecular linkage map of cultivated oat

A molecular linkage map of cultivated oat composed of 561 loci has been developed using 71 recombinant inbred lines from a cross between Avena byzantina cv. Kanota and A. sativa cv. Ogle. The loci are mainly restriction fragment length polymorphisms detected by oat cDNA clones from leaf, endosperm, and root tissue, as well as by barley leaf cDNA clones. The loci form 38 linkage groups ranging in size from 0.0 to 122.1 cM (mean, 39 cM) and consist of 2-51 loci each (mean, 14). Twenty-nine loci remain unlinked. The current map size is 1482 cM and the total size, on the basis of the number of unlinked loci, is estimated to be 2932.0 cM. This indicates that this map covers at least 50% of the cultivated oat genome. Comparisons with an A-genome diploid oat map and between linkage groups exhibiting homoeology to each other indicate that several major chromosomal rearrangements exist in cultivated oat. This map provides a tool for marker-assisted selection, quantitative trait loci analyses, and studies of genome organization in oat.

High-density SNP genotyping of tomato (Solanum lycopersicum L.) reveals patterns of genetic variation due to breeding

The effects of selection on genome variation were investigated and visualized in tomato using a high-density single nucleotide polymorphism (SNP) array. 7,720 SNPs were genotyped on a collection of 426 tomato accessions (410 inbreds and 16 hybrids) and over 97% of the markers were polymorphic in the entire collection. Principal component analysis (PCA) and pairwise estimates of F(st) supported that the inbred accessions represented seven sub-populations including processing, large-fruited fresh market, large-fruited vintage, cultivated cherry, landrace, wild cherry, and S. pimpinellifolium. Further divisions were found within both the contemporary processing and fresh market sub-populations. These sub-populations showed higher levels of genetic diversity relative to the vintage sub-population. The array provided a large number of polymorphic SNP markers across each sub-population, ranging from 3,159 in the vintage accessions to 6,234 in the cultivated cherry accessions. Visualization of minor allele frequency revealed regions of the genome that distinguished three representative sub-populations of cultivated tomato (processing, fresh market, and vintage), particularly on chromosomes 2, 4, 5, 6, and 11. The PCA loadings and F(st) outlier analysis between these three sub-populations identified a large number of candidate loci under positive selection on chromosomes 4, 5, and 11. The extent of linkage disequilibrium (LD) was examined within each chromosome for these sub-populations. LD decay varied between chromosomes and sub-populations, with large differences reflective of breeding history. For example, on chromosome 11, decay occurred over 0.8 cM for processing accessions and over 19.7 cM for fresh market accessions. The observed SNP variation and LD decay suggest that different patterns of genetic variation in cultivated tomato are due to introgression from wild species and selection for market specialization.