Characterization of the USDA Cucurbita pepo, C. moschata, and C. maxima germplasm collections.

The Cucurbita genus is home to a number of economically and culturally important species. We present the analysis of genotype data generated through genotyping-by-sequencing of the USDA germplasm collections of Cucurbita pepo, C. moschata, and C. maxima. These collections include a mixture of wild, landrace, and cultivated specimens from all over the world. Roughly 1,500 - 32,000 high-quality single nucleotide polymorphisms (SNPs) were called in each of the collections, which ranged in size from 314 to 829 accessions. Genomic analyses were conducted to characterize the diversity in each of the species. Analysis revealed extensive structure corresponding to a combination of geographical origin and morphotype/market class. Genome-wide associate studies (GWAS) were conducted using both historical and contemporary data. Signals were observed for several traits, but the strongest was for the bush (Bu) gene in C. pepo. Analysis of genomic heritability, together with population structure and GWAS results, was used to demonstrate a close alignment of seed size in C. pepo, maturity in C. moschata, and plant habit in C. maxima with genetic subgroups. These data represent a large, valuable collection of sequenced Cucurbita that can be used to direct the maintenance of genetic diversity, for developing breeding resources, and to help prioritize whole-genome re-sequencing.

Genetic Resources and Vulnerabilities of Major Cucurbit Crops.

The Cucurbitaceae family provides numerous important crops including watermelons (Citrullus lanatus), melons (Cucumis melo), cucumbers (Cucumis sativus), and pumpkins and squashes (Cucurbita spp.). Centers of domestication in Africa, Asia, and the Americas were followed by distribution throughout the world and the evolution of secondary centers of diversity. Each of these crops is challenged by multiple fungal, oomycete, bacterial, and viral diseases and insects that vector disease and cause feeding damage. Cultivated varieties are constrained by market demands, the necessity for climatic adaptations, domestication bottlenecks, and in most cases, limited capacity for interspecific hybridization, creating narrow genetic bases for crop improvement. This analysis of crop vulnerabilities examines the four major cucurbit crops, their uses, challenges, and genetic resources. ex situ germplasm banks, the primary strategy to preserve genetic diversity, have been extensively utilized by cucurbit breeders, especially for resistances to biotic and abiotic stresses. Recent genomic efforts have documented genetic diversity, population structure, and genetic relationships among accessions within collections. Collection size and accessibility are impacted by historical collections, current ability to collect, and ability to store and maintain collections. The biology of cucurbits, with insect-pollinated, outcrossing plants, and large, spreading vines, pose additional challenges for regeneration and maintenance. Our ability to address ongoing and future cucurbit crop vulnerabilities will require a combination of investment, agricultural, and conservation policies, and technological advances to facilitate collection, preservation, and access to critical Cucurbitaceae diversity.

The Evolutionary History of Wild, Domesticated, and Feral Brassica oleracea (Brassicaceae).

Understanding the evolutionary history of crops, including identifying wild relatives, helps to provide insight for conservation and crop breeding efforts. Cultivated Brassica oleracea has intrigued researchers for centuries due to its wide diversity in forms, which include cabbage, broccoli, cauliflower, kale, kohlrabi, and Brussels sprouts. Yet, the evolutionary history of this species remains understudied. With such different vegetables produced from a single species, B. oleracea is a model organism for understanding the power of artificial selection. Persistent challenges in the study of B. oleracea include conflicting hypotheses regarding domestication and the identity of the closest living wild relative. Using newly generated RNA-seq data for a diversity panel of 224 accessions, which represents 14 different B. oleracea crop types and nine potential wild progenitor species, we integrate phylogenetic and population genetic techniques with ecological niche modeling, archaeological, and literary evidence to examine relationships among cultivars and wild relatives to clarify the origin of this horticulturally important species. Our analyses point to the Aegean endemic B. cretica as the closest living relative of cultivated B. oleracea, supporting an origin of cultivation in the Eastern Mediterranean region. Additionally, we identify several feral lineages, suggesting that cultivated plants of this species can revert to a wild-like state with relative ease. By expanding our understanding of the evolutionary history in B. oleracea, these results contribute to a growing body of knowledge on crop domestication that will facilitate continued breeding efforts including adaptation to changing environmental conditions.

Genotyping by sequencing for SNP marker development in onion.

Onion (Allium cepa) is not highly tractable for development of molecular markers due to its large (16 gigabases per 1C) nuclear genome. Single nucleotide polymorphisms (SNPs) are useful for genetic characterization and marker-aided selection of onion because of codominance and common occurrence in elite germplasm. We completed genotyping by sequencing (GBS) to identify SNPs in onion using 46 F2 plants, parents of the F2 plants (Ailsa Craig 43 and Brigham Yellow Globe 15-23), two doubled haploid (DH) lines (DH2107 and DH2110), and plants from 94 accessions in the USDA National Plant Germplasm System (NPGS). SNPs were called using the TASSEL 3.0 Universal Network Enabled Analysis (UNEAK) bioinformatics pipeline. Sequences from the F2 and DH plants were used to construct a pseudo-reference genome against which genotypes from all accessions were scored. Quality filters were used to identify a set of 284 high quality SNPs, which were placed onto an existing genetic map for the F2 family. Accessions showed a moderate level of diversity (mean He = 0.341) and evidence of inbreeding (mean F = 0.592). GBS is promising for SNP discovery in onion, although lack of a reference genome required extensive custom scripts for bioinformatics analyses to identify high quality markers.

Population Structure and Phylogenetic Relationships in a Diverse Panel of Brassica rapa L.

The crop species Brassica rapa L. has significant economic importance around the world. However, the global distribution and complex evolutionary history of the species has made investigating its genetic population structure difficult. Crop domestication and improvement has resulted in extreme phenotypic diversity and subspecies that are used for oilseed, food for human consumption, and fodder for livestock. These subspecies include the oilseed morphotypes. oleifera (turnip rape), ssp. dichotoma (brown sarson/toria), ssp. trilocularis (yellow sarson); ssp. rapa (turnip); and Asian leafy vegetables ssp. pekinensis (Chinese cabbage), ssp. chinensis (bok choy), ssp. nipposinica (mizuna/mibuna), ssp. rapifera (rapini/broccoli rabe), ssp. narinosa (tatsoi), ssp parachinensis (choy sum), and ssp. perviridis (komatsuna). To date, studies have had insufficient sampling to determine the relationship of all morphotypes, especially oilseed morphotypes, and questions remain over the contribution of morphotype and geographic origin to population structure. We used genotyping-by-sequencing to score 18,272 single nucleotide polymorphism markers in a globally diverse panel of 333 B. rapa National Plant Germplasm System accessions that included 10 recognized subspecies. Our population genetic and phylogenetic analyses were broadly congruent and revealed five subpopulations that were largely reflective of morphotype and geography. These subpopulations were 1. European turnips/oilseed, 2. Asian turnips/oilseed, 3. yellow/brown sarson (ssp. trilocularis and ssp. dichotoma), 4. Chinese cabbage (ssp. pekinensis), and 5. bok choy, choy sum, and tatsoi (ssp. chinensis, ssp. parachinensis, ssp. narinosa). Additionally, we found evidence of polyphyly and/or paraphyly, particularly for oilseed morphotypes (ssp. oleifera and ssp. dichotoma) and turnips. The results of this study have provided improved resolution to the genetic and phylogenetic relationships of subspecies within the species B. rapa. Understanding of these relationships is key to the future genetic study and improvement of this globally important crop species.

Multiple loss-of-function 5-O-glucosyltransferase alleles revealed in Vitis vinifera, but not in other Vitis species.

KEY MESSAGE: Wild and loss-of-function alleles of the 5 - O - glucosyltransferase gene responsible for synthesis of diglucoside anthocyanins in Vitis were characterized. The information aids marker development for tracking this gene in grape breeding. Anthocyanins in red grapes are present in two glycosylation states: monoglucoside (3-O-glucoside) and diglucoside (3, 5-di-O-glucoside). While monoglucoside anthocyanins are present in all pigmented grapes, diglucoside anthocyanins are rarely found in the cultivated grape species Vitis vinifera. Biochemically 3-O-glucoside anthocyanins can be converted into 3,5-di-O-glucoside anthocyanins by a 5-O-glucosyltransferase. In this study, we surveyed allelic variation of the 5-O-glucosyltransferase gene (5GT) in 70 V. vinifera ssp. vinifera cultivars, 52 V. vinifera ssp. sylvestris accessions, 23 Vitis hybrid grapes, and 22 accessions of seven other Vitis species. Eighteen 5GT alleles with apparent loss-of-function mutations, including seven premature stop codon mutations and six frameshift indel mutations, were discovered in V. vinifera, but not in the other Vitis species. A total of 36 5GT alleles without apparent loss-of-function mutations (W-type) were identified. These W-type alleles were predominantly present in wild Vitis species, although a few of them were also found in some V. vinifera accessions. We further evaluated some of these 5GT alleles in producing diglucoside anthocyanins by analyzing the content of diglucoside anthocyanins in a set of representative V. vinifera cultivars. Through haplotype network analysis we revealed that V. vinifera ssp. vinifera and its wild progenitor V. vinifera ssp. sylvestris shared many loss-of-function 5GT alleles and extensive divergence of the 5GT alleles was evident within V. vinifera. This work advances our understanding of the genetic diversity of 5GT and provides a molecular basis for future marker-assisted selection for improving this important wine quality trait.

Genetic structure of the four wild tomato species in the Solanum peruvianum s.l. species complex.

The most diverse wild tomato species Solanum peruvianum sensu lato (s.l.) has been reclassified into four separate species: Solanum peruvianum sensu stricto (s.s.), Solanum corneliomuelleri, Solanum huaylasense, and Solanum arcanum. However, reproductive barriers among the species are incomplete and this can lead to discrepancies regarding genetic identity of germplasm. We used genotyping by sequencing (GBS) of S. peruvianum s.l., Solanum neorickii, and Solanum chmielewskii to develop tens of thousands of mapped single nucleotide polymorphisms (SNPs) to analyze genetic relationships within and among species. The data set was condensed to 14,043 SNPs with no missing data across 46 sampled plants. Origins of accessions were mapped using geographical information systems (GIS). Isolation by distance, pairwise genetic distances, and number of clusters were estimated using population genetics approaches. Isolation by distance was strongly supported, especially between interspecific pairs. Eriopersicon (S. peruvianum s.s., S. corneliomuelleri, S. huaylasense) and Arcanum (S. arcanum, S. neorickii, S. chmielewskii) species groups were genetically distinct, except for S. huaylasense which showed 50% membership proportions in each group. Solanum peruvianum and S. corneliomuelleri were not significantly differentiated from each other. Many thousands of SNP markers were identified that could potentially be used to distinguish pairs of species, including S. peruvianum versus S. corneliomuelleri, if they are verified on larger numbers of samples. Diagnostic markers will be valuable for delimiting morphologically similar and interfertile species in germplasm management. Approximately 12% of the SNPs rejected a genome-wide test of selective neutrality based on differentiation among species of S. peruvianum s.l. These are candidates for more comprehensive studies of microevolutionary processes within this species complex.

A phylogenetic analysis of the grape genus (Vitis L.) reveals broad reticulation and concurrent diversification during neogene and quaternary climate change.

BACKGROUND: Grapes are one of the most economically important fruit crops. There are about 60 species in the genus Vitis. The phylogenetic relationships among these species are of keen interest for the conservation and use of this germplasm. We selected 309 accessions from 48 Vitis species,varieties, and outgroups, examined ~11 kb (~3.4 Mb total) of aligned nuclear DNA sequences from 27 unlinked genes in a phylogenetic context, and estimated divergence times based on fossil calibrations. RESULTS: Vitis formed a strongly supported clade. There was substantial support for species and less for the higher-level groupings (series). As estimated from extant taxa, the crown age of Vitis was 28 Ma and the divergence of subgenera (Vitis and Muscadinia) occurred at ~18 Ma. Higher clades in subgenus Vitis diverged 16 - 5 Ma with overlapping confidence intervals, and ongoing divergence formed extant species at 12 - 1.3 Ma. Several species had species-specific SNPs. NeighborNet analysis showed extensive reticulation at the core of subgenus Vitis representing the deeper nodes, with extensive reticulation radiating outward. Fitch Parsimony identified North America as the origin of the most recent common ancestor of extant Vitis species. CONCLUSIONS: Phylogenetic patterns suggested origination of the genus in North America, fragmentation of an ancestral range during the Miocene, formation of extant species in the late Miocene-Pleistocene, and differentiation of species in the context of Pliocene-Quaternary tectonic and climatic change. Nuclear SNPs effectively resolved relationships at and below the species level in grapes and rectified several misclassifications of accessions in the repositories. Our results challenge current higher-level classifications, reveal the abundance of genetic diversity in the genus that is potentially available for crop improvement, and provide a valuable resource for species delineation, germplasm conservation and use.

Evidence of cryptic introgression in tomato (Solanum lycopersicum L.) based on wild tomato species alleles.

BACKGROUND: Many highly beneficial traits (e.g. disease or abiotic stress resistance) have been transferred into crops through crosses with their wild relatives. The 13 recognized species of tomato (Solanum section Lycopersicon) are closely related to each other and wild species genes have been extensively used for improvement of the crop, Solanum lycopersicum L. In addition, the lack of geographical barriers has permitted natural hybridization between S. lycopersicum and its closest wild relative Solanum pimpinellifolium in Ecuador, Peru and northern Chile. In order to better understand patterns of S. lycopersicum diversity, we sequenced 47 markers ranging in length from 130 to 1200 bp (total of 24 kb) in genotypes of S. lycopersicum and wild tomato species S. pimpinellifolium, Solanum arcanum, Solanum peruvianum, Solanum pennellii and Solanum habrochaites. Between six and twelve genotypes were comparatively analyzed per marker. Several of the markers had previously been hypothesized as carrying wild species alleles within S. lycopersicum, i.e., cryptic introgressions. RESULTS: Each marker was mapped with high confidence (e<1 x 10-30) to a single genomic location using BLASTN against tomato whole genome shotgun chromosomes (SL2.40) database. Neighbor-joining trees showed high mean bootstrap support (86.8 ± 2.34%) for distinguishing red-fruited from green-fruited taxa for 38 of the markers. Hybridization and parsimony splits networks, genomic map positions of markers relative to documented introgressions, and historical origins of accessions were used to interpret evolutionary patterns at nine markers with putatively introgressed alleles. CONCLUSION: Of the 47 genetic markers surveyed in this study, four were involved in linkage drag on chromosome 9 during introgression breeding, while alleles at five markers apparently originated from natural hybridization with S. pimpinellifolium and were associated with primitive genotypes of S. lycopersicum. The positive identification of introgressed genes within crop species such as S. lycopersicum will help inform conservation and utilization of crop germplasm diversity, for example, facilitating the purging of undesirable linkage drag or the exploitation of novel, favorable alleles.

EST, COSII, and arbitrary gene markers give similar estimates of nucleotide diversity in cultivated tomato (Solanum lycopersicum L.).

Because cultivated tomato (Solanum lycopersicum L.) is low in genetic diversity, public, verified single nucleotide polymorphism (SNP) markers within the species are in demand. To promote marker development we resequenced approximately 23 kb in a diverse set of 31 tomato lines including TA496. Three classes of markers were sampled: (1) 26 expressed-sequence tag (EST), all of which were predicted to be polymorphic based on TA496, (2) 14 conserved ortholog set II (COSII) or unigene, and (3) ten published sequences, composed of nine fruit quality genes and one anonymous RFLP marker. The latter two types contained mostly noncoding DNA. In total, 154 SNPs and 34 indels were observed. The distributions of nucleotide diversity estimates among marker types were not significantly different from each other. Ascertainment bias of SNPs was evaluated for the EST markers. Despite the fact that the EST markers were developed using SNP prediction within a sample consisting of only one TA496 allele and one additional allele, the majority of polymorphisms in the 26 EST markers were represented among the other 30 tomato lines. Fifteen EST markers with published SNPs were more closely examined for bias. Mean SNP diversity observations were not significantly different between the original discovery sample of two lines (53 SNPs) and the 31 line diversity panel (56 SNPs). Furthermore, TA496 shared its haplotype with at least one other line at 11 of the 15 markers. These data demonstrate that public EST databases and noncoding regions are a valuable source of unbiased SNP markers in tomato.