The USDA cucumber (Cucumis sativus L.) collection: genetic diversity, population structure, genome-wide association studies, and core collection development

Genome-wide association studies to assess genetic factors controlling cucumber resistance to CABYV and CMV in crop fields and the attractiveness for their Aphis gossypii vector

Cucumber crops face high pressure from pathogens, including various viral species. Mapping quantitative trait loci (QTL) for vegetable resistance to viruses has primarily been conducted after mechanical inoculation in controlled environments, but not in crop field conditions. Moreover, viruses that cannot be mechanically inoculated, e.g. the cucurbit aphid-borne yellows virus (CABYV), have been overlooked in resistance studies. Here, we aimed to identify QTLs reducing epidemics of two prevalent cucumber viruses: CABYV and the cucumber mosaic virus (CMV). We evaluated the resistance of 256 elite cucumber lines and landraces in crop field conditions by screening for the presence of both viruses six-times during the season. We mapped twelve QTLs reducing CABYV epidemics and seven QTLs reducing CMV epidemics by combining multiloci genome-wide association studies and local score approach analyses. We also examined the attractiveness of this cucumber panel for Aphis gossypii, a major cucumber virus vector. We identified five QTLs that reduced the attractiveness, including one co-localizing with a QTL reducing CABYV epidemics. Interestingly, some accessions deemed CMV-resistant after mechanical inoculation in controlled environments showed high infection rates in crop field conditions. Only one QTL for CMV resistance was detected in both conditions, indicating that these phenotypes are regulated by independent QTLs. Local linkage disequilibrium study findings suggested that certain QTLs reducing epidemics were introduced early into elite lines through serendipity or selection. QTLs could be pyramided with other low-effect QTLs through genomic selection to obtain cucumber cultivars with enhanced field resistance to viruses.

Clade V MLO Genes Are Negative Modulators of Cucumber Defence Response to Meloidogyne incognita

Cucumber production is seriously constrained by Meloidogyne incognita. Because no resistance resources to the pathogen have been reported, disabling susceptibility genes may represent a novel breeding strategy to introduce resistance against this nematode in cucumber. Here, we studied the clade V MLO genes for their involvement in the interaction between cucumber and M. incognita. Our results showed that Arabidopsis clade V MLO mutants were resistant to M. incognita. Cucumber has three clade V MLO genes, CsaMLO1, CsaMLO8 and CsaMLO11, with upregulated expression upon inoculation with M. incognita. Heterologous overexpression of CsaMLO1, CsaMLO8 and CsaMLO11 in Arabidopsis mutants restored susceptibility to varying degrees. Silencing and knockout of individual clade V MLO genes in cucumber reduced susceptibility to M. incognita. The cucumber CRISPR mutants produced similar fruits as the wild type (WT) did. Although the yields of two single mutants (M111 and M112) and two double mutants (M81 M111 and M81 M112) were reduced compared to WT, the yields of M81 and M82 were not decreased. In summary, clade V MLO genes function as susceptibility genes for M. incognita in cucumber. Among them, CsaMLO8 may be the most promising candidate for M. incognita resistance breeding in cucumber.

A single nucleotide substitution introducing premature stop codon within CsTFL1 explains the determinate-2 phenotype in cucumber (Cucumis sativus L.)

The determinate growth habit of plants reduces the number of internodes and shortens the main stem by terminating the shoot apical meristem through a transition to inflorescence. Understanding the genetic basis of this habit can help optimize crop yield and cultivation technology for vegetable breeding. This study aimed to identify the determinate-2 (de-2) gene responsible for the determinate growth habit in the W-sk cucumber line. Termination of the main stem in the W-sk line occurred between 14 and 23 internodes, depending on cultivation conditions. Resequencing of the W-sk genome identified a novel SNP in the cucumber TERMINAL FLOWER1 (CsTFL1) gene, explaining the de-2 phenotype. This was verified with a CAPS-T marker cosegregation with determinate growth in the F2:3 population, and this polymorphism is unique among genotyped indeterminate cucumber cultivars or breeding lines. Crossing the W-sk line with the G421 line with the determinate (de) gene confirmed the allelism of both genes. An SNP in CsTFL1 in the W-sk line introduced a premature stop codon, resulting in the putative deletion of 13 amino acids, possibly causing determinate growth habit. Overall, this study provides insights into the genetic basis of cucumber plant growth architecture and advances in cucumber breeding.

Genome-Wide Identification of the Cyclic Nucleotide-Gated Ion Channel Gene Family and Expression Profiles Under Low-Temperature Stress in Luffa cylindrica L

Cyclic nucleotide-gated ion channels (CNGCs) are cell membrane channel proteins for calcium ions. They have been reported to play important roles in survival and in the responses to environmental factors in various plants. However, little is known about the CNGC family and its functions in luffa (Luffa cylindrica L.). In this study, a bioinformatics-based method was used to identify members of the CNGC gene family in L. cylindrica. In total, 20 LcCNGCs were detected, and they were grouped into five subfamilies (I, II, Ⅲ, IV-a, and IV-b) in a phylogenetic analysis with CNGCs from Arabidopsis thaliana (20 AtCNGCs) and Momordica charantia (17 McCNGCs). The 20 LcCNGC genes were unevenly distributed on 11 of the 13 chromosomes in luffa, with none on Chromosomes 1 and 5. The members of each subfamily encoded proteins with highly conserved functional domains. An evolutionary analysis of CNGCs in luffa revealed three gene losses and a motif deletion. An examination of gene replication events during evolution indicated that two tandemly duplicated gene pairs were the primary driving force behind the evolution of the LcCNGC gene family. PlantCARE analyses of the LcCNGC promoter regions revealed various cis-regulatory elements, including those responsive to plant hormones (abscisic acid, methyl jasmonate, and salicylic acid) and abiotic stresses (light, drought, and low temperature). The presence of these cis-acting elements suggested that the encoded CNGC proteins may be involved in stress responses, as well as growth and development. Transcriptome sequencing (RNA-seq) analyses revealed tissue-specific expression patterns of LcCNGCs in various plant parts (roots, stems, leaves, flowers, and fruit) and the upregulation of some LcCNGCs under low-temperature stress. To confirm the accuracy of the RNA-seq data, 10 cold-responsive LcCNGC genes were selected for verification by quantitative real-time polymerase chain reaction (RT-qPCR) analysis. Under cold conditions, LcCNGC4 was highly upregulated (>50-fold increase in its transcript levels), and LcCNGC3, LcCNGC6, and LcCNGC13 were upregulated approximately 10-fold. Our findings provide new information about the evolution of the CNGC family in L. cylindrica and provide insights into the functions of the encoded CNGC proteins.

A genome-wide association study reveals molecular mechanism underlying powdery mildew resistance in cucumber

BACKGROUND

Powdery mildew is a disease with one of the most substantial impacts on cucumber production globally. The most efficient approach for controlling powdery mildew is the development of genetic resistance; however, few genes associated with inherent variations in cucumber powdery mildew resistance have been identified as of yet.

RESULTS

In this study, we re-sequence 299 cucumber accessions, which are divided into four geographical groups. A genome-wide association study identifies 50 sites significantly associated with natural variations in powdery mildew resistance. Linkage disequilibrium analysis further divides these 50 sites into 32 linkage disequilibrium blocks containing 41 putative genes. Virus-induced gene silencing and gene expression analysis implicate CsGy5G015960, which encodes a phosphate transporter, as the candidate gene regulating powdery mildew resistance. On the basis of the resequencing data, we generate five CsGy5G015960 haplotypes, identifying Hap.1 as the haplotype most likely associated with powdery mildew resistance. In addition, we determine that a 29-bp InDel in the 3' untranslated region of CsGy5G015960 is responsible for mRNA stability. Overexpression of CsGy5G015960Hap.1 in the susceptible line enhances powdery mildew resistance and phosphorus accumulation. Further comparative RNA-seq analysis demonstrates that CsGy5G015960Hap.1 may regulate cucumber powdery mildew resistance by maintaining a higher H2O2 level through the depletion of multiple class III peroxidases.

CONCLUSIONS

Here we identify a candidate powdery mildew-resistant gene in cucumber using GWAS. The identified gene may be a promising target for molecular breeding and genetic engineering in cucumber to enhance powdery mildew resistance.

Genome-Wide Association Study of Cuticle and Lipid Droplet Properties of Cucumber (Cucumis sativus L.) Fruit

The fruit surface is a critical first line of defense against environmental stress. Overlaying the fruit epidermis is the cuticle, comprising a matrix of cutin monomers and waxes that provides protection and mechanical support throughout development. The epidermal layer of the cucumber (Cucumis sativus L.) fruit also contains prominent lipid droplets, which have recently been recognized as dynamic organelles involved in lipid storage and metabolism, stress response, and the accumulation of specialized metabolites. Our objective was to genetically characterize natural variations for traits associated with the cuticle and lipid droplets in cucumber fruit. Phenotypic characterization and genome-wide association studies (GWAS) were performed using a resequenced cucumber core collection accounting for >96% of the allelic diversity present in the U.S. National Plant Germplasm System collection. The collection was grown in the field, and fruit were harvested at 16-20 days post-anthesis, an age when the cuticle thickness and the number and size of lipid droplets have stabilized. Fresh fruit tissue sections were prepared to measure cuticle thickness and lipid droplet size and number. The collection showed extensive variation for the measured traits. GWAS identified several QTLs corresponding with genes previously implicated in cuticle or lipid biosynthesis, including the transcription factor SHINE1/WIN1, as well as suggesting new candidate genes, including a potential lipid-transfer domain containing protein found in association with isolated lipid droplets.

Genome-wide association analysis reveal candidate genes and haplotypes related to root weight in cucumber (Cucumis sativus L.)

Background

The plant root system is critical for the absorption of water and nutrients, and have a direct influence on growth and yield. In cucumber, a globally consumed crop, the molecular mechanism of root development remains unclear, and this has implications for developing stress tolerant varieties. This study sought to determine the genetic patterns and related genes of cucumber root weight. A core cucumber germplasms population was used to do the GWAS analysis in three environments.

Results

Here, we investigated four root-weight related traits including root fresh weight (RFW), root dry weight (RDW), ratio of root dry weight to root fresh weight (RDFW) and the comprehensive evaluation index, D-value of root weight (DRW) deduced based on the above three traits for the core germplasm of the cucumber global repository. According to the D-value, we identified 21 and 16 accessions with light and heavy-root, respectively. We also found that the East Asian ecotype accessions had significantly heavier root than other three ecotypes. The genome-wide association study (GWAS) for these four traits reveals that 4 of 10 significant loci (gDRW3.1, gDRW3.2, gDRW4.1 and gDRW5.1) were repeatedly detected for at least two traits. Further haplotype and expression analysis for protein-coding genes positioned within these 4 loci between light and heavy-root accessions predicted five candidate genes (i.e., Csa3G132020 and Csa3G132520 both encoding F-box protein PP2-B1 for gDRW3.1, Csa3G629240 encoding a B-cell receptor-associated protein for gDRW3.2, Csa4G499330 encodes a GTP binding protein for gDRW4.1, and Csa5G286040 encodes a proteinase inhibitor for gDRW5.1).

Conclusions

We conducted a systematic analysis of the root genetic basis and characteristics of cucumber core germplasms population. We detected four novel loci, which regulate the root weight in cucumber. Our study provides valuable candidate genes and haplotypes for the improvement of root system in cucumber breeding.

Genetic characterization of cucumber genetic resources in the NARO Genebank indicates their multiple dispersal trajectories to the East

KEY MESSAGE

Genotyping-by-sequencing of 723 worldwide cucumber genetic resources revealed that cucumbers were dispersed eastward via at least three distinct routes, one to Southeast Asia and two from different directions to East Asia. The cucumber (Cucumis sativus) is an economically important vegetable crop cultivated and consumed worldwide. Despite its popularity, the manner in which cucumbers were dispersed from their origin in South Asia to the rest of the world, particularly to the east, remains a mystery due to the lack of written records. In this study, we performed genotyping-by-sequencing (GBS) on 723 worldwide cucumber accessions, mainly deposited in the Japanese National Agriculture and Food Research Organization (NARO) Genebank, to characterize their genetic diversity, relationships, and population structure. Analyses based on over 60,000 genome-wide single-nucleotide polymorphisms identified by GBS revealed clear genetic differentiation between Southeast and East Asian populations, suggesting that they reached their respective region independently, not progressively. A deeper investigation of the East Asian population identified two subpopulations with different fruit characteristics, supporting the traditional classification of East Asian cucumbers into two types thought to have been introduced by independent routes. Finally, we developed a core collection of 100 accessions representing at least 93.2% of the genetic diversity present in the entire collection. The genetic relationships and population structure, their associations with geographic distribution and phenotypic traits, and the core collection presented in this study are valuable resources for elucidating the dispersal history and promoting the efficient use and management of genetic resources for research and breeding in cucumber.

Lipid-Related Domestication Accounts for the Extreme Cold Sensitivity of Semiwild and Tropic Xishuangbanna Cucumber (Cucumis sativus L. var. xishuangbannanesis)

Xishuangbanna (XIS) cucumber (Cucumis sativus L. var. xishuangbannanesis) is a semiwild variety originating from low latitude tropic areas, and therefore shows extreme cold sensitivity and heat tolerance. Here, we mapped the quantitative trait loci (QTLs) that control the cold sensitivity and heat tolerance of XIS cucumber seedlings. Using bulked segregant analysis (BSA), we identified three QTLs (HTT1.1, HTT3.1, and HTT3.2, with a total length of 11.98 Mb) for heat tolerance and two QTLs (LTT6.1 and LTT6.2, with a total length of 8.74 Mb) for cold sensitivity. The QTL LTT6.1 was then narrowed down to a length of 641 kb by using kompetitive allele-specific PCR (KASP) markers. Based on structural variants (SVs) and single-nucleotide polymorphisms (SNPs), we found the LTT6.1 is covered by a high divergent region including a 50 kb deletion in the XIS49 genome, which affects the gene structure of lipase abhydrolase domain containing 6 (ABHD6, Csa_6G032560). Accordingly, there is a very big difference in lipid composition, but not in other osmoprotectants like free amino acids and fatty acids, between XIS49 and cultivated cucumber CL. Moreover, we calculated the composite likelihood ratio (CLR) and identified selective sweeps from 115 resequencing data, and found that lipid- and fatty-acid-related processes are major aspects in the domestication of the XIS group cucumber. LTT6.1 is a particularly special region positioned nearby lipid-related selective sweeps. These studies above suggested that the lipid-related domestication of XIS cucumbers should account for their extreme cold sensitivity.

Identification of QTL associated with resistance to Phytophthora fruit rot in cucumber (Cucumis sativus L.)

Phytophthora fruit rot (PFR) caused by the soilborne oomycete pathogen, Phytophthora capsici, can cause severe yield loss in cucumber. With no resistant variety available, genetic resources are needed to develop resistant varieties. The goal of this work was to identify quantitative trait loci (QTL) associated with resistance to PFR using multiple genomic approaches and populations. Two types of resistances have been identified: age-related resistance (ARR) and young fruit resistance. ARR occurs at 12-16 days post pollination (dpp), coinciding with the end of exponential fruit growth. A major QTL for ARR was discovered on chromosome 3 and a candidate gene identified based on comparative transcriptomic analysis. Young fruit resistance, which is observed during the state of rapid fruit growth prior to commercial harvest, is a quantitative trait for which multiple QTL were identified. The largest effect QTL, qPFR5.1, located on chromosome 5 was fine mapped to a 1-Mb region. Genome-wide association studies (GWAS) and extreme-phenotype genome-wide association study (XP-GWAS) for young fruit resistance were also performed on a cucumber core collection representing > 96% of the genetic diversity of the USDA cucumber germplasm. Several SNPs overlapped with the QTL identified from QTL-seq analysis on biparental populations. In addition, novel SNPs associated with the resistance were identified from the germplasm. The resistant alleles were found mostly in accessions from India and South Asia, the center of diversity for cucumber. The results from this work can be applied to future disease resistance studies and marker-assisted selection in breeding programs.

ShinyCore: An R/Shiny program for establishing core collection based on single nucleotide polymorphism data

BACKGROUND

Managing and investigating all available genetic resources are challenging. As an alternative, breeders and researchers use core collection-a representative subset of the entire collection. A good core is characterized by high genetic diversity and low repetitiveness. Among the several available software, GenoCore uses a coverage criterion that does not require computationally expensive distance-based metrics.

RESULTS

ShinyCore is a new method to select a core collection through two phases. The first phase uses the coverage criterion to quickly attain a fixed coverage, and the second phase uses a newly devised score (referred to as the rarity score) to further enhance diversity. It can attain a fixed coverage faster than a currently available algorithm devised for the coverage criterion, so it will benefit users who have big data. ShinyCore attains the minimum coverage specified by a user faster than GenoCore, and it then seeks to add entries with the rarest allele for each marker. Therefore, measures of genetic diversity and distance can be improved.

CONCLUSION

Although GenoCore is a fast algorithm, its implementation is difficult for those unfamiliar with R, ShinyCore can be easily implemented in Shiny with RStudio and an interactive web applet is available for those who are not familiar with programming languages.

Genetic diversity, population structure, and genome-wide association analysis of ginkgo cultivars

Ginkgo biloba is an economically valuable tree worldwide. The species has nearly become extinct during the Quaternary, which has likely resulted in reduction of its genetic variability. The genetic variability is now conserved in few natural populations in China and a number of cultivars that are, however, derived from a few ancient trees, helping the species survive in China through medieval times. Despite the recent interest in ginkgo, however, detailed knowledge of its genetic diversity, conserved in cultivated trees and cultivars, has remained poor. This limits efficient conservation of its diversity as well as efficient use of the existing germplasm resources. Here we performed genotyping-by-sequencing (GBS) on 102 cultivated germplasms of ginkgo collected to explore their genetic structure, kinship, and inbreeding prediction. For the first time in ginkgo, a genome-wide association analysis study (GWAS) was used to attempt gene mapping of seed traits. The results showed that most of the germplasms did not show any obvious genetic relationship. The size of the ginkgo germplasm population expanded significantly around 1500 years ago during the Sui and Tang dynasties. Classification of seed cultivars based on a phylogenetic perspective does not support the current classification criteria based on phenotype. Twenty-four candidate genes were localized after performing GWAS on the seed traits. Overall, this study reveals the genetic basis of ginkgo seed traits and provides insights into its cultivation history. These findings will facilitate the conservation and utilization of the domesticated germplasms of this living fossil plant.

Elucidation of genetic variation and population structure of melon genetic resources in the NARO Genebank, and construction of the World Melon Core Collection

Numerous genetic resources of major crops have been introduced from around the world and deposited in Japanese National Agriculture and Food Research Organization (NARO) Genebank. Understanding their genetic variation and selecting a representative subset ("core collection") are essential for optimal management and efficient use of genetic resources. In this study, we conducted genotyping-by-sequencing (GBS) to characterize the genetic relationships and population structure in 755 accessions of melon genetic resources. The GBS identified 39,324 single-nucleotide polymorphisms (SNPs) that are distributed throughout the melon genome with high density (one SNP/10.6 kb). The phylogenetic relationships and population structure inferred using this SNP dataset are highly associated with the cytoplasm type and geographical origin. Our results strongly support the recent hypothesis that cultivated melon was established in Africa and India through multiple independent domestication events. Finally, we constructed a World Melon Core Collection that covers at least 82% of the genetic diversity and has a wide range of geographical origins and fruit morphology. The genome-wide SNP dataset, phylogenetic relationships, population structure, and the core collection provided in this study should largely contribute to genetic research, breeding, and genetic resource preservation in melon.

Chloroplast Pan-Genomes and Comparative Transcriptomics Reveal Genetic Variation and Temperature Adaptation in the Cucumber

Although whole genome sequencing, genetic variation mapping, and pan-genome studies have been done on a large group of cucumber nuclear genomes, organelle genome information is largely unclear. As an important component of the organelle genome, the chloroplast genome is highly conserved, which makes it a useful tool for studying plant phylogeny, crop domestication, and species adaptation. Here, we have constructed the first cucumber chloroplast pan-genome based on 121 cucumber germplasms, and investigated the genetic variations of the cucumber chloroplast genome through comparative genomic, phylogenetic, haplotype, and population genetic structure analysis. Meanwhile, we explored the changes in expression of cucumber chloroplast genes under high- and low-temperature stimulation via transcriptome analysis. As a result, a total of 50 complete chloroplast genomes were successfully assembled from 121 cucumber resequencing data, ranging in size from 156,616-157,641 bp. The 50 cucumber chloroplast genomes have typical quadripartite structures, consisting of a large single copy (LSC, 86,339-86,883 bp), a small single copy (SSC, 18,069-18,363 bp), and two inverted repeats (IRs, 25,166-25,797 bp). Comparative genomic, haplotype, and population genetic structure results showed that there is more genetic variation in Indian ecotype cucumbers compared to other cucumber cultivars, which means that many genetic resources remain to be explored in Indian ecotype cucumbers. Phylogenetic analysis showed that the 50 cucumber germplasms could be classified into 3 types: East Asian, Eurasian + Indian, and Xishuangbanna + Indian. The transcriptomic analysis showed that matK were significantly up-regulated under high- and low-temperature stresses, further demonstrating that cucumber chloroplasts respond to temperature adversity by regulating lipid metabolism and ribosome metabolism. Further, accD has higher editing efficiency under high-temperature stress, which may contribute to the heat tolerance. These studies provide useful insight into genetic variation in the chloroplast genome, and established the foundation for exploring the mechanisms of temperature-stimulated chloroplast adaptation.

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.

Elucidating the role of key physio-biochemical traits and molecular network conferring heat stress tolerance in cucumber

Cucumber is an important vegetable crop grown worldwide and highly sensitive to prevailing temperature condition. The physiological, biochemical and molecular basis of high temperature stress tolerance is poorly understood in this model vegetable crop. In the present study, a set of genotypes with contrasting response under two different temperature stress (35/30°C and 40/35°C) were evaluated for important physiological and biochemical traits. Besides, expression of the important heat shock proteins (HSPs), aquaporins (AQPs), photosynthesis related genes was conducted in two selected contrasting genotypes at different stress conditions. It was established that tolerant genotypes were able to maintain high chlorophyll retention, stable membrane stability index, higher retention of water content, stability in net photosynthesis, high stomatal conductance and transpiration in combination with less canopy temperatures under high temperature stress conditions compared to susceptible genotypes and were considered as the key physiological traits associated with heat tolerance in cucumber. Accumulation of biochemicals like proline, protein and antioxidants like SOD, catalase and peroxidase was the underlying biochemical mechanisms for high temperature tolerance. Upregulation of photosynthesis related genes, signal transduction genes and heat responsive genes (HSPs) in tolerant genotypes indicate the molecular network associated with heat tolerance in cucumber. Among the HSPs, higher accumulation of HSP70 and HSP90 were recorded in the tolerant genotype, WBC-13 under heat stress condition indicating their critical role. Besides, Rubisco S, Rubisco L and CsTIP1b were upregulated in the tolerant genotypes under heat stress condition. Therefore, the HSPs in combination with photosynthetic and aquaporin genes were the underlying important molecular network associated with heat stress tolerance in cucumber. The findings of the present study also indicated negative feedback of G-protein alpha unit and oxygen evolving complex in relation to heat stress tolerance in cucumber. These results indicate that the thermotolerant cucumber genotypes enhanced physio-biochemical and molecular adaptation under high-temperature stress condition. This study provides foundation to design climate smart genotypes in cucumber through integration of favorable physio-biochemical traits and understanding the detailed molecular network associated with heat stress tolerance in cucumber.

CuGenDBv2: an updated database for cucurbit genomics

The Cucurbitaceae (cucurbit) family consists of about 1,000 species in 95 genera, including many economically important and popular fruit and vegetable crops. During the past several years, reference genomes have been generated for >20 cucurbit species, and variome and transcriptome profiling data have been rapidly accumulated for cucurbits. To efficiently mine, analyze and disseminate these large-scale datasets, we have developed an updated version of Cucurbit Genomics Database. The updated database, CuGenDBv2 (http://cucurbitgenomics.org/v2), currently hosts 34 reference genomes from 27 cucurbit species/subspecies belonging to 10 different genera. Protein-coding genes from these genomes have been comprehensively annotated by comparing their protein sequences to various public protein and domain databases. A novel 'Genotype' module has been implemented to facilitate mining and analysis of the functionally annotated variome data including SNPs and small indels from large-scale genome sequencing projects. An updated 'Expression' module has been developed to provide a comprehensive gene expression atlas for cucurbits. Furthermore, synteny blocks between any two and within each of the 34 genomes, representing a total of 595 pair-wise genome comparisons, have been identified and can be explored and visualized in the database.

Morphological and Genetic Diversity of Cucumber (Cucumis sativus L.) Fruit Development

Cucumber (Cucumis sativus L.) fruits, which are eaten at an immature stage of development, can vary extensively in morphological features such as size, shape, waxiness, spines, warts, and flesh thickness. Different types of cucumbers that vary in these morphological traits are preferred throughout the world. Numerous studies in recent years have added greatly to our understanding of cucumber fruit development and have identified a variety of genetic factors leading to extensive diversity. Candidate genes influencing floral organ establishment, cell division and cell cycle regulation, hormone biosynthesis and response, sugar transport, trichome development, and cutin, wax, and pigment biosynthesis have all been identified as factors influencing cucumber fruit morphology. The identified genes demonstrate complex interplay between structural genes, transcription factors, and hormone signaling. Identification of genetic factors controlling these traits will facilitate breeding for desired characteristics to increase productivity, improve shipping, handling, and storage traits, and enhance consumer-desired qualities. The following review examines our current understanding of developmental and genetic factors driving diversity of cucumber fruit morphology.

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

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

Unravelling cucumber resistance to several viruses via genome-wide association studies highlighted resistance hotspots and new QTLs

The mapping and introduction of sustainable resistance to viruses in crops is a major challenge in modern breeding, especially regarding vegetables. We hence assembled a panel of cucumber elite lines and landraces from different horticultural groups for testing with six virus species. We mapped 18 quantitative trait loci (QTL) with a multiloci genome wide association studies (GWAS), some of which have already been described in the literature. We detected two resistance hotspots, one on chromosome 5 for resistance to the cucumber mosaic virus (CMV), cucumber vein yellowing virus (CVYV), cucumber green mottle mosaic virus (CGMMV) and watermelon mosaic virus (WMV), colocalizing with the RDR1 gene, and another on chromosome 6 for resistance to the zucchini yellowing mosaic virus (ZYMV) and papaya ringspot virus (PRSV) close to the putative VPS4 gene location. We observed clear structuring of resistance among horticultural groups due to plant virus coevolution and modern breeding which have impacted linkage disequilibrium (LD) in resistance QTLs. The inclusion of genetic structure in GWAS models enhanced the GWAS accuracy in this study. The dissection of resistance hotspots by local LD and haplotype construction helped gain insight into the panel's resistance introduction history. ZYMV and CMV resistance were both introduced from different donors in the panel, resulting in multiple resistant haplotypes at same locus for ZYMV, and in multiple resistant QTLs for CMV.

The genomic signature of wild-to-crop introgression during the domestication of scarlet runner bean (Phaseolus coccineus L.)

The scarlet runner bean (Phaseolus coccineus) is one of the five domesticated Phaseolus species. It is cultivated in small-scale agriculture in the highlands of Mesoamerica for its dry seeds and immature pods, and unlike the other domesticated beans, P. coccineus is an open-pollinated legume. Contrasting with its close relative, the common bean, few studies focusing on its domestication history have been conducted. Demographic bottlenecks associated with domestication might reduce genetic diversity and facilitate the accumulation of deleterious mutations. Conversely, introgression from wild relatives could be a source of variation. Using Genotyping by Sequencing data (79,286 single-nucleotide variants) from 237 cultivated and wild samples, we evaluated the demographic history of traditional varieties from different regions of Mexico and looked for evidence of introgression between sympatric wild and cultivated populations. Traditional varieties have high levels of diversity, even though there is evidence of a severe initial genetic bottleneck followed by a population expansion. Introgression from wild to domesticated populations was detected, which might contribute to the recovery of the genetic variation. Introgression has occurred at different times: constantly in the center of Mexico; recently in the North West; and anciently in the South. Several factors are acting together to increase and maintain genetic diversity in P. coccineus cultivars, such as demographic expansion and introgression. Wild relatives represent a valuable genetic resource and have played a key role in scarlet runner bean evolution via introgression into traditional varieties.

Genome-Wide Association Studies Reveal Candidate Genes Related to Stem Diameter in Cucumber (Cucumis sativus L.)

The stem diameter, an important agronomic trait, affects cucumber growth and yield. However, no genes responsible for cucumber stem diameter have been identified yet. In this study, the stem diameter of 88 cucumber core germplasms were measured in spring 2020, autumn 2020 and autumn 2021, and a genome-wide association study (GWAS) was carried out based on the gene sequence and stem diameter of core germplasms. A total of eight loci (gSD1.1, gSD2.1, gSD3.1, gSD3.2, gSD4.1, gSD5.1, gSD5.2, and gSD6.1) significantly associated with cucumber stem diameter were detected. Of these, five loci (gSD1.1, gSD2.1, gSD3.1, gSD5.2, and gSD6.1) were repeatedly detected in two or more seasons and were considered as robust and reliable loci. Based on the linkage disequilibrium sequences of the associated SNP loci, 37 genes were selected. By further investigating the five loci via analyzing Arabidopsis homologous genes and gene haplotypes, five genes (CsaV3_1G028310, CsaV3_2G006960, CsaV3_3G009560, CsaV3_5G031320, and CsaV3_6G031260) showed variations in amino acid sequence between thick stem lines and thin stem lines. Expression pattern analyses of these genes also showed a significant difference between thick stem and thin stem lines. This study laid the foundation for gene cloning and molecular mechanism study of cucumber stem development.

Functional copy number variation of CsSHINE1 is associated with fruit skin netting intensity in cucumber, Cucumis sativus

Fruit skin netting in cucumber (Cucumis sativus) is associated with important fruit quality attributes. Two simply inherited genes H (Heavy netting) and Rs (Russet skin) control skin netting, but their molecular basis is unknown. Here, we reported map-based cloning and functional characterization of the candidate gene for the Rs locus that encodes CsSHINE1 (CsSHN1), an AP2 domain containing ethylene-responsive transcription factor protein. Comparative phenotypic analysis in near-isogenic lines revealed that fruit with netted skin had different epidermal structures from that with smooth skin including thicker cuticles, smaller, palisade-shaped epidermal and sub-epidermal cells with heavily suberized and lignified cell walls, higher peroxidase activities, which suggests multiple functions of CsSHN1 in regulating fruit skin netting and epidermal cell patterning. Among three representative cucumber inbred lines, three haplotypes at three polymorphic sites were identified inside CsSHN1: a functional copy in Gy14 (wild type) with light fruit skin netting, a copy number variant with two tandemly arrayed functional copies in WI7120 with heavy skin netting, and a loss-of-function copy in 9930 with smooth skin. The expression level of CsSHN1 in fruit exocarp of three lines was positively correlated with the skin netting intensity. Comparative analysis between cucumber and melon revealed conserved and divergent genetic mechanisms underlying fruit skin netting/reticulation that may reflect the different selection histories in the two crops. A discussion was made on genetic basis of fruit skin netting in the context of natural and artificial selections of fruit quality-related epidermal features during cucumber breeding.

NS encodes an auxin transporter that regulates the 'numerous spines' trait in cucumber (Cucumis sativus) fruit

Fruit spine is an important agronomic trait in cucumber and the "numerous spines (ns)" cucumber varieties are popular in Europe and West Asia. Although the classical genetic locus of ns was reported more than two decades ago, the NS gene has not been cloned yet. In this study, nine genetic loci for the different densities of fruit spines were identified by a genome-wide association study. Among the nine loci, fsdG2.1 was closely associated with the classical genetic locus ns, which harbors a candidate gene Csa2G264590. Overexpression of Csa2G264590 resulted in lower fruit spine density, and the knockout mutant generated by CRISPR/Cas9 displayed an increased spine density, demonstrating that the Csa2G264590 gene is NS. NS is specifically expressed in the fruit peel and spine. Genetic analysis showed that NS regulates fruit spine development independently of the tuberculate gene, Tu, which regulates spine development on tubercules; the cucumber glabrous mutants csgl1 and csgl3 are epistatic to ns. Furthermore, we found that auxin levels in the fruit peel and spine were significantly lower in the knockout mutant ns-cr. Moreover, RNA-sequencing showed that the plant hormone signal transduction pathway was enriched. Notably, most of the auxin responsive Aux/IAA family genes were downregulated in ns-cr. Haplotype analysis showed that the non-functional haplotype of NS exists exclusively in the Eurasian cucumber backgrounds. Taken together, the cloning of NS gene provides new insights into the regulatory network of fruit spine development.

Heat Stress Resistance Mechanisms of Two Cucumber Varieties from Different Regions

High temperatures affect the yield and quality of vegetable crops. Unlike thermosensitive plants, thermotolerant plants have excellent systems for withstanding heat stress. This study evaluated various heat resistance indexes of the thermotolerant cucumber (TT) and thermosensitive cucumber (TS) plants at the seedling stage. The similarities and differences between the regulatory genes were assessed through transcriptome analysis to understand the mechanisms for heat stress resistance in cucumber. The TT plants exhibited enhanced leaf status, photosystem, root viability, and ROS scavenging under high temperature compared to the TS plants. Additionally, transcriptome analysis showed that the genes involved in photosynthesis, the chlorophyll metabolism, and defense responses were upregulated in TT plants but downregulated in TS plants. Zeatin riboside (ZR), brassinosteroid (BR), and jasmonic acid (JA) levels were higher in TT plants than in TS. The heat stress increased gibberellic acid (GA) and indoleacetic acid (IAA) levels in both plant lines; however, the level of GA was higher in TT. Correlation and interaction analyses revealed that heat cucumber heat resistance is regulated by a few transcription factor family genes and metabolic pathways. Our study revealed different phenotypic and physiological mechanisms of the heat response by the thermotolerant and thermosensitive cucumber plants. The plants were also shown to exhibit different expression profiles and metabolic pathways. The heat resistant pathways and genes of two cucumber varieties were also identified. These results enhance our understanding of the molecular mechanisms of cucumber response to high-temperature stress.

Graph-based pan-genome reveals structural and sequence variations related to agronomic traits and domestication in cucumber

Structural variants (SVs) represent a major source of genetic diversity and are related to numerous agronomic traits and evolutionary events; however, their comprehensive identification and characterization in cucumber (Cucumis sativus L.) have been hindered by the lack of a high-quality pan-genome. Here, we report a graph-based cucumber pan-genome by analyzing twelve chromosome-scale genome assemblies. Genotyping of seven large chromosomal rearrangements based on the pan-genome provides useful information for use of wild accessions in breeding and genetic studies. A total of ~4.3 million genetic variants including 56,214 SVs are identified leveraging the chromosome-level assemblies. The pan-genome graph integrating both variant information and reference genome sequences aids the identification of SVs associated with agronomic traits, including warty fruits, flowering times and root growth, and enhances the understanding of cucumber trait evolution. The graph-based cucumber pan-genome and the identified genetic variants provide rich resources for future biological research and genomics-assisted breeding.

Increasing studies have suggested that single reference genome is insufficient to capture all variations in the genome. Here, the authors report a graph-based cucumber pan-genome by analyzing 12 chromosome-scale assemblies and reveal variations associated with agronomic traits and domestication.

Transcriptome analysis of ovary culture-induced embryogenesis in cucumber (Cucumis sativus L.)

Background.

Ovary culture is a useful technique used to generate double haploid (DH) cucumber (Cucumis sativus L.) plants. However, cucumber ovary culture have a low rate of embryo induction and plant regeneration. Moreover, the cucumber embryogenesis mechanism remains unclear. In this study, we explored the molecular basis of cucumber embryogenesis in order to establish a foundation for a more efficient ovary culture method. Using transcriptome sequencing, we also investigated the differential expression of genes during the embryogenesis process.

Methods.

Cytological and morphological observations have divided cucumber ovary culture into three stages: early embryo development (T0), embryo morphogenesis (T1, T2, T3 and T4), and shoot formation (T5). We selected six key time points for transcriptome sequencing and analysis: T0 (the ovules were cultured for 0 d), T1 (the ovules were cultured for 2 d), T2 (the embryos were cultured for 10 d), T3 (the embryos were cultured for 20 d), T4 (the embryos were cultured for 30 d), and T5 (the shoots after 60 d culture).

Results.

We used cytology and morphology to observe the characteristics of the cucumber’s developmental transformation during embryogenesis and plant regeneration. The differentially expressed genes(DEGs) at developmental transition points were analyzed using transcriptome sequencing. In the early embryogenesis stage, the cells expanded, which was the signal for gametophytes to switch to the sporophyte development pathway. RNA-seq revealed that when compared to the fresh unpollinated ovaries, there were 3,468 up-regulated genes in the embryos, including hormone signal transduction genes, hormone response genes, and stress-induced genes. The reported embryogenesis-related genes BBM, HSP90 and AGL were also actively expressed during this stage. In the embryo morphogenesis stage (from cell division to cotyledon-embryo formation), 480 genes that functioned in protein complex binding, microtubule binding, tetrapyrrole binding, tubulin binding and other microtubule activities were continuously up-regulated during the T1, T2, T3 and T4 time points. This indicated that the cytoskeleton structure was continuously being built and maintained by the action of microtubule-binding proteins and enzyme modification. In the shoot formation stage, 1,383 genes were up-regulated that were mainly enriched in phenylpropanoid biosynthesis, plant hormone signal transduction, phenylalanine metabolism, and starch and sucrose metabolism. These up-regualted genes included six transcription factors that contained a B3 domain, nine genes in the AP2/ERF family, and two genes encoding WUS homologous domain proteins.

Conclusions.

Evaluation of molecular gynogenesis events may contribute to a better understanding of the molecular mechanism of cucumber ovarian culture.

Identification and Validation of a Core Single-Nucleotide Polymorphism Marker Set for Genetic Diversity Assessment, Fingerprinting Identification, and Core Collection Development in Bottle Gourd

Germplasm collections are indispensable resources for the mining of important genes and variety improvement. To preserve and utilize germplasm collections in bottle gourd, we identified and validated a highly informative core single-nucleotide polymorphism (SNP) marker set from 1,100 SNPs. This marker set consisted of 22 uniformly distributed core SNPs with abundant polymorphisms, which were established to have strong representativeness and discriminatory power based on analyses of 206 bottle gourd germplasm collections and a multiparent advanced generation inter-cross (MAGIC) population. The core SNP markers were used to assess genetic diversity and population structure, and to fingerprint important accessions, which could provide an optimized procedure for seed authentication. Furthermore, using the core SNP marker set, we developed an accessible core population of 150 accessions that represents 100% of the genetic variation in bottle gourds. This core population will make an important contribution to the preservation and utilization of bottle gourd germplasm collections, cultivar identification, and marker-assisted breeding.

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.

Resistant Sources and Genetic Control of Resistance to ToLCNDV in Cucumber

Tomato leaf curl New Delhi virus (ToLCNDV) is a severe threat for cucurbit production worldwide. Resistance has been reported in several crops, but at present, there are no described accessions with resistance to ToLCNDV in cucumber (Cucumis sativus). C. sativus var. sativus accessions were mechanically inoculated with ToLCNDV and screened for resistance, by scoring symptom severity, tissue printing, and PCR (conventional and quantitative). Severe symptoms and high load of viral DNA were found in plants of a nuclear collection of Spanish landraces and in accessions of C. sativus from different geographical origins. Three Indian accessions (CGN23089, CGN23423, and CGN23633) were highly resistant to the mechanical inoculation, as well as all plants of their progenies obtained by selfing. To study the inheritance of the resistance to ToLCNDV, plants of the CGN23089 accession were crossed with the susceptible accession BGV011742, and F₁ hybrids were used to construct segregating populations (F₂ and backcrosses), which were mechanically inoculated and evaluated for symptom development and viral load by qPCR. The analysis of the genetic control fit with a recessive monogenic inheritance model, and after genotyping with SNPs distributed along the C. sativus genome, a QTL associated with ToLCNDV resistance was identified in chromosome 2 of cucumber.

Genetic Diversity Assessment and Cultivar Identification of Cucumber (Cucumis sativus L.) Using the Fluidigm Single Nucleotide Polymorphism Assay

Genetic diversity analysis and cultivar identification were performed using a core set of single nucleotide polymorphisms (SNPs) in cucumber (Cucumis sativus L.). For the genetic diversity study, 280 cucumber accessions collected from four continents (Asia, Europe, America, and Africa) by the National Agrobiodiversity Center of the Rural Development Administration in South Korea and 20 Korean commercial F1 hybrids were genotyped using 151 Fluidigm SNP assay sets. The heterozygosity of the SNP loci per accession ranged from 4.76 to 82.76%, with an average of 32.1%. Population genetics analysis was performed using population structure analysis and hierarchical clustering (HC), which indicated that these accessions were classified mainly into four subpopulations or clusters according to their geographical origins. The subpopulations for Asian and European accessions were clearly distinguished from each other (F_ST value = 0.47), while the subpopulations for Korean F1 hybrids and Asian accessions were closely related (F_ST = 0.34). The highest differentiation was observed between American and European accessions (F_ST = 0.41). Nei's genetic distance among the 280 accessions was 0.414 on average. In addition, 95 commercial F1 hybrids of three cultivar groups (Baekdadagi-, Gasi-, and Nakhap-types) were genotyped using 82 Fluidigm SNP assay sets for cultivar identification. These 82 SNPs differentiated all cultivars, except seven. The heterozygosity of the SNP loci per cultivar ranged from 12.20 to 69.14%, with an average of 34.2%. Principal component analysis and HC demonstrated that most cultivars were clustered based on their cultivar groups. The Baekdadagi- and Gasi-types were clearly distinguished, while the Nakhap-type was closely related to the Baekdadagi-type. Our results obtained using core Fluidigm SNP assay sets provide useful information for germplasm assessment and cultivar identification, which are essential for breeding and intellectual right protection in cucumber.

Complete chloroplast genomes shed light on phylogenetic relationships, divergence time, and biogeography of Allioideae (Amaryllidaceae)

Allioideae includes economically important bulb crops such as garlic, onion, leeks, and some ornamental plants in Amaryllidaceae. Here, we reported the complete chloroplast genome (cpDNA) sequences of 17 species of Allioideae, five of Amaryllidoideae, and one of Agapanthoideae. These cpDNA sequences represent 80 protein-coding, 30 tRNA, and four rRNA genes, and range from 151,808 to 159,998 bp in length. Loss and pseudogenization of multiple genes (i.e., rps2, infA, and rpl22) appear to have occurred multiple times during the evolution of Alloideae. Additionally, eight mutation hotspots, including rps15-ycf1, rps16-trnQ-UUG, petG-trnW-CCA, psbA upstream, rpl32-trnL-UAG, ycf1, rpl22, matK, and ndhF, were identified in the studied Allium species. Additionally, we present the first phylogenomic analysis among the four tribes of Allioideae based on 74 cpDNA coding regions of 21 species of Allioideae, five species of Amaryllidoideae, one species of Agapanthoideae, and five species representing selected members of Asparagales. Our molecular phylogenomic results strongly support the monophyly of Allioideae, which is sister to Amaryllioideae. Within Allioideae, Tulbaghieae was sister to Gilliesieae-Leucocoryneae whereas Allieae was sister to the clade of Tulbaghieae- Gilliesieae-Leucocoryneae. Molecular dating analyses revealed the crown age of Allioideae in the Eocene (40.1 mya) followed by differentiation of Allieae in the early Miocene (21.3 mya). The split of Gilliesieae from Leucocoryneae was estimated at 16.5 mya. Biogeographic reconstruction suggests an African origin for Allioideae and subsequent spread to Eurasia during the middle Eocene. Cool and arid conditions during the late Eocene led to isolation between African and Eurasian species. African Allioideae may have diverged to South American taxa in the late Oligocene. Rather than vicariance, long-distance dispersal is the most likely explanation for intercontinental distribution of African and South American Allioideae species.

Plant-ImputeDB: an integrated multiple plant reference panel database for genotype imputation

Genotype imputation is a process that estimates missing genotypes in terms of the haplotypes and genotypes in a reference panel. It can effectively increase the density of single nucleotide polymorphisms (SNPs), boost the power to identify genetic association and promote the combination of genetic studies. However, there has been a lack of high-quality reference panels for most plants, which greatly hinders the application of genotype imputation. Here, we developed Plant-ImputeDB (http://gong_lab.hzau.edu.cn/Plant_imputeDB/), a comprehensive database with reference panels of 12 plant species for online genotype imputation, SNP and block search and free download. By integrating genotype data and whole-genome resequencing data of plants from various studies and databases, the current Plant-ImputeDB provides high-quality reference panels of 12 plant species, including ∼69.9 million SNPs from 34 244 samples. It also provides an easy-to-use online tool with the option of two popular tools specifically designed for genotype imputation. In addition, Plant-ImputeDB accepts submissions of different types of genomic variations, and provides free and open access to all publicly available data in support of related research worldwide. In general, Plant-ImputeDB may serve as an important resource for plant genotype imputation and greatly facilitate the research on plant genetic research.

Quantitative trait loci for horticulturally important traits defining the Sikkim cucumber, Cucumis sativus var. sikkimensis

QTL mapping identified simply inherited genes and quantitative trait loci underlying morphologically characteristic traits of the Sikkim cucumber, which reveals their genetic basis during crop evolution. The data suggest the Sikkim cucumber as an ecotype of cultivated cucumber not worthy of formal taxonomic recognition. The Sikkim cucumber, Cucumis sativus var. sikkimensis, is featured with some morphological traits like black spine, brown fruit with fine and heavy netting, as well as large hollow in mature fruit. Despite its establishment as a botanical variety ~ 150 years ago, and its wide use as an important source of disease resistances in cucumber breeding, little is known about its taxonomic status and genetic basis of those characteristic traits. Here we reported QTL mapping with segregating populations derived from two Sikkim-type inbred lines, WI7088D and WI7120, and identification of 48 QTL underlying phenotypic variation for 18 horticulturally important traits. We found that the fruit spine and skin colors in the two populations were controlled by the previously cloned pleiotropic B (black spine) locus. The fruit netting in WI7088D and WI7120 was controlled by the well-known H (Heavy netting) and a novel Rs (Russet skin) locus, which was delimited to a 271-kb region on Chr5 and ~ 736-kb region on Chr1, respectively. A single major-effect QTL was detected for flowering time in each population (ft1.1 for WI7088D and ft6.2 for WI7120). Fifteen, six and five QTL were identified for fruit size, hollow size and flesh thickness variation in the two populations, respectively. No major structural changes were found between the Sikkim and cultivated cucumbers. Except for the rare allele at the Rs locus, there seem no private QTL/alleles identified from this study supporting the Sikkim cucumber as an ecotype of C. sativus, not worthy of formal taxonomic recognition.

Identification of Novel Loci and Candidate Genes for Cucumber Downy Mildew Resistance Using GWAS

Downy mildew (DM) is one of the most serious diseases in cucumber. Multiple quantitative trait loci (QTLs) for DM resistance have been detected in a limited number of cucumber accessions. In this study we applied genome-wide association analysis (GWAS) to detected genetic loci for DM resistance in a core germplasm (CG) of cucumber lines that represent diverse origins and ecotypes. Phenotypic data on responses to DM infection were collected in four field trials across three years, 2014, 2015, and 2016. With the resequencing data of these CG lines, GWAS for DM resistance was performed and detected 18 loci that were distributed on all the seven cucumber chromosomes. Of these 18 loci, only six (dmG1.4, dmG4.1, dmG4.3, dmG5.2, dmG7.1, and dmG7.2) were detected in two experiments, and were considered as loci with a stable effect on DM resistance. Further, 16 out of the 18 loci colocalized with the QTLs that were reported in previous studies and two loci, dmG2.1 and dmG7.1, were novel ones identified only in this study. Based on the annotation of homologous genes in Arabidopsis and pairwise LD correlation analysis, several candidate genes were identified as potential causal genes underlying the stable and novel loci, including Csa1G575030 for dmG1.4, Csa2G060360 for dmG2.1, Csa4G064680 for dmG4.1, Csa5G606470 for dmG5.2, and Csa7G004020 for dmG7.1. This study shows that the CG germplasm is a very valuable resource carrying known and novel QTLs for DM resistance. The potential of using these CG lines for future allele-mining of candidate genes was discussed in the context of breeding cucumber with resistance to DM.

Mapping a Partial Andromonoecy Locus in Citrullus lanatus Using BSA-Seq and GWAS Approaches

The sexual expression of watermelon plants is the result of the distribution and occurrence of male, female, bisexual and hermaphrodite flowers on the main and secondary stems. Plants can be monoecious (producing male and female flowers), andromonoecious (producing male and hermaphrodite flowers), or partially andromonoecious (producing male, female, bisexual, and hermaphrodite flowers) within the same plant. Sex determination of individual floral buds and the distribution of the different flower types on the plant, are both controlled by ethylene. A single missense mutation in the ethylene biosynthesis gene CitACS4, is able to promote the conversion of female into hermaphrodite flowers, and therefore of monoecy (genotype MM) into partial andromonoecy (genotype Mm) or andromonoecy (genotype mm). We phenotyped and genotyped, for the M/m locus, a panel of 207 C. lanatus accessions, including five inbreds and hybrids, and found several accessions that were repeatedly phenotyped as PA (partially andromonoecious) in several locations and different years, despite being MM. A cosegregation analysis between a SNV in CitACS4 and the PA phenotype, demonstrated that the occurrence of bisexual and hermaphrodite flowers in a PA line is not dependent on CitACS4, but conferred by an unlinked recessive gene which we called pa. Two different approaches were performed to map the pa gene in the genome of C. lanatus: bulk segregant analysis sequencing (BSA-seq) and genome wide association analysis studies (GWAS). The BSA-seq study was performed using two contrasting bulks, the monoecious M-bulk and the partially andromonoecious PA-bulk, each one generated by pooling DNA from 20 F2 plants. For GWAS, 122 accessions from USDA gene bank, already re-sequenced by genotyping by sequencing (GBS), were used. The combination of the two approaches indicates that pa maps onto a genomic region expanding across 32.24-36.44 Mb in chromosome 1 of watermelon. Fine mapping narrowed down the pa locus to a 867 Kb genomic region containing 101 genes. A number of candidate genes were selected, not only for their function in ethylene biosynthesis and signalling as well as their role in flower development and sex determination, but also by the impact of the SNPs and indels differentially detected in the two sequenced bulks.

The major-effect quantitative trait locus Fnl7.1 encodes a late embryogenesis abundant protein associated with fruit neck length in cucumber

Fruit neck length (FNL) is an important quality trait in cucumber because it directly affects its market value. However, its genetic basis remains largely unknown. We identified a candidate gene for FNL in cucumber using a next-generation sequencing-based bulked segregant analysis in F₂ populations, derived from a cross between Jin5-508 (long necked) and YN (short necked). A quantitative trait locus (QTL) on chromosome 7, Fnl7.1, was identified through a genome-wide comparison of single nucleotide polymorphisms between long and short FNL F₂ pools, and it was confirmed by traditional QTL mapping in multiple environments. Fine genetic mapping, sequences alignment and gene expression analysis revealed that CsFnl7.1 was the most likely candidate Fnl7.1 locus, which encodes a late embryogenesis abundant protein. The increased expression of CsFnl7.1 in long-necked Jin5-508 may be attributed to mutations in the promoter region upstream of the gene body. The function of CsFnl7.1 in FNL control was confirmed by its overexpression in transgenic cucumbers. CsFnl7.1 regulates fruit neck development by modulating cell expansion. Probably, this is achieved through the direct protein-protein interactions between CsFnl7.1 and a dynamin-related protein CsDRP6 and a germin-like protein CsGLP1. Geographical distribution differences of the FNL phenotype were found among the different cucumber types. The East Asian and Eurasian cucumber accessions were highly enriched with the long-necked and short-necked phenotypes, respectively. A further phylogenetic analysis revealed that the Fnl7.1 locus might have originated from India. Thus, these data support that the CsFnl7.1 has an important role in increasing cucumber FNL.

Recommendations for Choosing the Genotyping Method and Best Practices for Quality Control in Crop Genome-Wide Association Studies

High-throughput genotyping boosts genome-wide association studies (GWAS) in crop species, leading to the identification of single-nucleotide polymorphisms (SNPs) associated with economically important traits. Choosing a cost-effective genotyping method for crop GWAS requires careful examination of several aspects, namely, the purpose and the scale of the study, crop-specific genomic features, and technical and economic matters associated with each genotyping option. Once genotypic data have been obtained, quality control (QC) procedures must be applied to avoid bias and false signals in genotype–phenotype association tests. QC for human GWAS has been extensively reviewed; however, QC for crop GWAS may require different actions, depending on the GWAS population type. Here, we review most popular genotyping methods based on next-generation sequencing (NGS) and array hybridization, and report observations that should guide the investigator in the choice of the genotyping method for crop GWAS. We provide recommendations to perform QC in crop species, and deliver an overview of bioinformatics tools that can be used to accomplish all needed tasks. Overall, this work aims to provide guidelines to harmonize those procedures leading to SNP datasets ready for crop GWAS.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a ‘SyntenyViewer’ to view genome synteny between different cucurbit species and an ‘RNA-Seq’ module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

Gynoecy instability in cucumber (Cucumis sativus L.) is due to unequal crossover at the copy number variation-dependent Femaleness (F) locus

Cucumber, Cucumis sativus is an important vegetable crop, and gynoecy has played a critical role in yield increase of hybrid cucumber production. Cucumber has a unique genetic system for gynoecious sex expression, which is determined by the copy number variation (CNV)-based, dominant, and dosage-dependent femaleness (F) locus. However, this gynoecy expression system seems unstable since monecious plants could often be found in F-dependent gynoecious cucumber inbreds. We hypothesized that gynoecy instability (gynoecy loss) may be due to unequal crossing over (UCO) during meiosis among repeat units of the CNV. In this study, using high throughput genome resequencing, fiber-FISH and genomic qPCR analyses, we first confirmed and refined the structure of the F locus, which was a CNV of a 30.2-kb tandem repeat. Gynoecious plants contained three genes: CsACS1, CsACS1G, and CsMYB, of which CsACS1G is a duplication of CsACS1 but with a recombinant distal promoter that may contribute to gynoecy sex expression. In two large populations from self-pollinated gynoecious inbred lines, 'gynoecy loss' mutants were identified with similar mutation rates (~0.12%). We show that these monecious mutants have lost CsACS1G. In addition, we identified gynoecious lines in natural populations that carry two copies of CSACS1G. We proposed a model to explain gynoecy instability in F-dependent cucumbers, which is caused by UCO among CSACS1/G units during meiosis. The findings present a convincing case that the phenotypic variation of an economically important trait is associated with the dynamic changes of copy numbers at the F locus. This work also has important implications in cucumber breeding.

FLOWERING LOCUS T Improves Cucumber Adaptation to Higher Latitudes

Flowering time plays a crucial role in the geographical adaptation of most crops during domestication. Cucumber (Cucumis sativus) is a major vegetable crop worldwide. From its tropical origin on the southern Asian continent, cucumber has spread over a wide latitudinal cline, but the molecular mechanisms underlying this latitudinal adaptation and the expansion of domesticated cucumber are largely unclear. Here, we report the cloning of two flowering time loci from two distinct cucumber populations and show that two large deletions upstream from FLOWERING LOCUS T (FT) are associated with higher expression of FT and earlier flowering. We determined that the two large deletions are pervasive and occurred independently in Eurasian and East-Asian populations. Nucleotide diversity analysis further revealed that the FT locus region of the cucumber genome contains a signature for a selective sweep during domestication. Our results suggest that large genetic structural variations upstream from FT were selected for and have been important in the geographic spread of cucumber from its tropical origin to higher latitudes.

Constructing a Core Collection of the Medicinal Plant Angelica biserrata Using Genetic and Metabolic Data

Angelica biserrata is an important medicinal plant in Chinese traditional medicine. Its roots, which are known as Duhuo in Chinese, are broadly applied to treat inflammation, arthritis, and headache. With increasing market demand, the wild resources of A. biserrata have been overexploited, and conservation, assessment of genetic resources and breeding for this species is needed. Here, we sequenced the transcriptome of A. biserrata and developed simple sequence repeat (SSR) markers from it to construct a core collection based on 208 samples collected from Changyang-related regions. A total of 132 alleles were obtained for 17 SSR loci used with the polymorphic information content (PIC) ranging from 0.44 to 0.83. Abundant genetic diversity was inferred by Shannon's information index (1.51), observed (0.57) and expected heterozygosity (0.72). The clustering analysis resulted into two sample groups and analysis of molecular variance (AMOVA) showed only 6% genetic variation existed among populations. A further metabolic analysis of these samples revealed the main coumarin contents, such as osthole and columbianadin. According to the genetic and metabolic data, we adopted the least distance stepwise sampling strategy to construct seven preliminary core collections, of which the 20CC collection, which possessed 42 A. biserrata individuals accounting for 90.20% of the genetic diversity of the original germplasm, represented the best core collection. This study will contribute to the conservation and management of A. biserrata wild germplasm resources and provide a material basis for future selection and breeding of this medicinal plant.

Comparison Between Core Set Selection Methods Using Different Illumina Marker Platforms: A Case Study of Assessment of Diversity in Wheat

Collections of plant genetic resources stored in genebanks are an important source of genetic diversity for improvement in plant breeding programs and for conservation of natural variation. The establishment of reduced representative collections from a large set of genotypes is a valuable tool that provides cost-effective access to the diversity present in the whole set. Software like Core Hunter 3 is available to generate high quality core sets. In addition, general clustering approaches, e.g., k-medoids, are available to subdivide a large data set into small groups with maximum genetic diversity between groups. Illumina genotyping platforms are a very efficient tool for the assessment of genetic diversity of plant genetic resources. The accumulation of genotyping data over time using commercial genotyping platforms raises the question of how such huge amount of information can be efficiently used for creating core collections. In the present study, after developing a 15K wheat Infinium array with 12,908 SNPs and genotyping a set of 479 hexaploid winter wheat lines (Triticum aestivum), a larger data set was created by merging 411 lines previously genotyped with the 90K iSelect array. Overlaying the markers from the 15K and 90K arrays enabled the identification of a common set of 12,806 markers, suggesting that the 15K array is a valuable and cost-effective resource for plant breeding programs. Finally, we selected genetically diverse core sets out of these 890 wheat genotypes derived from five collections based on the common markers from the 15K and 90K SNP arrays. Two different approaches, k-medoids and Core Hunter 3 were compared,and k-medoids was identified as an efficient method for selecting small core sets out of a large collection of genotypes while retaining the genetic diversity of the original population.

Selection footprints reflect genomic changes associated with breeding efforts in 56 cucumber inbred lines

Cucumber selective breeding over recent decades has dramatically increased productivity and quality, but the genomic characterizations and changes associated with this breeding history remain unclear. Here, we analyzed the genome resequencing data of 56 artificially selected cucumber inbred lines that exhibit various phenotypes to detect trait-associated sequence variations that reflect breeding improvement. We found that the 56 cucumber lines could be assigned to group 1 and group 2, and the two groups formed a distinctive genetic structure due to the breeding history involving hybridization and selection. Differentially selected regions were identified between group 1 and group 2, with implications for genomic-selection breeding signatures. These regions included known quantitative trait loci or genes that were reported to be associated with agronomic traits. Our results advance knowledge of cucumber genomics, and the 56 selected inbred lines could be good germplasm resources for breeding.

Two types of mutations in the HEUKCHEEM gene functioning in cucumber spine color development can be used as signatures for cucumber domestication

The HEUKCHEEM gene plays an important role in spine color formation. A white spine occurs due to two mutations in HEUKCHEEM and is closely related to the regional distribution of these mutants. Mapping analysis revealed that the HEUKCHEEM gene is co-segregated with the B locus in the regulation of black spine color development in cucumber fruit. HEUKCHEEM induced the expression of the genes involved in the anthocyanin biosynthetic pathway, leading to the accumulation of anthocyanins in black spines. The transiently over-expressed HEUKCHEEM in cucumber and tobacco plants enhanced the expression of anthocyanin biosynthesis-related genes, leading to anthocyanin accumulation. However, two mutations-insertion of the 6994 bp mutator-like transposable element (MULE) sequence into the second intron and one single-nucleotide polymorphism (SNP) of C to T in the second exon of HEUKCHEEM-were identified in white spines, leading to no accumulation of anthocyanin biosynthesis-related gene transcripts and anthocyanins. Furthermore, association analysis using 104 cucumber accessions with different geographical origins revealed that the types of mutations in HEUKCHEEM are strongly linked to geographical origins. The MULE insertion is found extensively in cucumbers with white spines in East Asia and Australia. However, cucumbers with white spines in other areas could be significantly influenced by a single SNP mutation. Our results provide fundamental information on spine color development in cucumber fruits and spine color-based cucumber breeding programs.

A Functional Allele of CsFUL1 Regulates Fruit Length through Repressing CsSUP and Inhibiting Auxin Transport in Cucumber

Fruit length is a prominent agricultural trait during cucumber (Cucumis sativus) domestication and diversifying selection; however, the regulatory mechanisms of fruit elongation remain elusive. We identified two alleles of the FRUITFULL (FUL)-like MADS-box gene CsFUL1 with 3393C/A Single Nucleotide Polymorphism variation among 150 cucumber lines. Whereas CsFUL1A was specifically enriched in the long-fruited East Asian type cucumbers (China and Japan), the CsFUL1C allele was randomly distributed in cucumber populations, including wild and semiwild cucumbers. CsFUL1A knockdown led to further fruit elongation in cucumber, whereas elevated expression of CsFUL1A resulted in significantly shorter fruits. No effect on fruit elongation was detected when CsFUL1C expression was modulated, suggesting that CsFUL1A is a gain-of-function allele in long-fruited cucumber that acts as a repressor during diversifying selection of East Asian cucumbers. Furthermore, CsFUL1A binds to the CArG-box in the promoter region of SUPERMAN, a regulator of cell division and expansion, to repress its expression. Additionally, CsFUL1A inhibits the expression of auxin transporters PIN-FORMED1 (PIN1) and PIN7, resulting in decreases in auxin accumulation in fruits. Together, our work identifies an agriculturally important allele and suggests a strategy for manipulating fruit length in cucumber breeding that involves modulation of CsFUL1A expression.

Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops

The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a 'SyntenyViewer' to view genome synteny between different cucurbit species and an 'RNA-Seq' module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system.

QTL and Transcriptomic Analyses Implicate Cuticle Transcription Factor SHINE as a Source of Natural Variation for Epidermal Traits in Cucumber Fruit

The fruit surface is a unique tissue with multiple roles influencing fruit development, post-harvest storage and quality, and consumer acceptability. Serving as the first line of protection against herbivores, pathogens, and abiotic stress, the surface can vary markedly among species, cultivars within species, and developmental stage. In this study we explore developmental changes and natural variation of cucumber (Cucumis sativus L.) fruit surface properties using two cucumber lines which vary greatly for these traits and for which draft genomes and a single nucleotide polymorphism (SNP) array are available: Chinese fresh market type, Chinese Long '9930' (CL9930), and pickling type, 'Gy14'. Thin-section samples were prepared from the mid-region of fruit harvested at 0, 4, 8, 12, 16, 20, 24 and 30 days post pollination (dpp), stained with Sudan IV and evaluated for cuticle thickness, depth of wax intercalation between epidermal cells, epidermal cell size and shape, and number and size of lipid droplets. 'Gy14' is characterized by columnar shaped epidermal cells, a 2-3 fold thicker cuticular layer than CL9930, increased cuticular intercalations between cells and a larger number and larger sized lipid droplets. In both lines maximal deposition of cuticle and increase in epidermal size coincided with exponential fruit growth and was largely completed by approximately 16 dpp. Phenotyping and quantitative trait locus mapping (QTL) of fruit sampled from an F₇:F₈ Gy14 × CL9930 recombinant inbred line (RIL) population identified QTL regions on chromosomes 1, 4 and 5. Strong QTL for epidermal cell height, cuticle thickness, intercalation depth, and diameter of lipid droplets co-localized on chromosome 1. SSR markers on chromosome 1 were used to screen for recombinants in an extended RIL population to refine the QTL region. Further fine mapping by KASP assay combined with gene expression profiling suggested a small number of candidate genes. Tissue specificity, developmental analysis of expression, allelic diversity and gene function implicate the regulatory factor CsSHINE1/WIN1 as a source of natural variation for cucumber fruit epidermal traits.