Chromosome rearrangements during domestication of cucumber as revealed by high-density genetic mapping and draft genome assembly

Chromosome-level genome assembly and population genomics reveals crucial selection for subgynoecy development in chieh-qua

Chieh-qua is an important cucurbit crop and very popular in South China and Southeast Asia. Despite its significance, its genetic basis and domestication history are unclear. In this study, we have successfully generated a chromosome-level reference genome assembly for the chieh-qua 'A36' using a hybrid assembly strategy that combines PacBio long reads and Illumina short reads. The assembled genome of chieh-qua is approximately 953.3 Mb in size and is organized into 12 chromosomes, with contig N50 of 6.9 Mb and scaffold N50 of 68.2 Mb. Notably, the chieh-qua genome is comparable in size to the wax gourd genome. Through gene prediction analysis, we have identified a total of 24 593 protein-coding genes in the A36 genome. Additionally, approximately 56.6% (539.3 Mb) of the chieh-qua genome consists of repetitive sequences. Comparative genome analysis revealed that chieh-qua and wax gourd are closely related, indicating a close evolutionary relationship between the two species. Population genomic analysis, employing 129 chieh-qua accessions and 146 wax gourd accessions, demonstrated that chieh-qua exhibits greater genetic diversity compared to wax gourd. We also employed the GWAS method to identify related QTLs associated with subgynoecy, an interested and important trait in chieh-qua. The MYB59 (BhiCQ0880026447) exhibited relatively high expression levels in the shoot apex of four subgynoecious varieties compared with monoecious varieties. Overall, this research provides insights into the domestication history of chieh-qua and offers valuable genomic resources for further molecular research.

A nanopore-based cucumber genome assembly reveals structural variations at two QTLs controlling hypocotyl elongation

The Xishuangbanna (XIS) cucumber (Cucumis sativus var. xishuangbannanesis) is a semiwild variety that has many distinct agronomic traits. Here, long reads generated by Nanopore sequencing technology helped assembling a high-quality genome (contig N50 = 8.7 Mb) of landrace XIS49. A total of 10,036 structural/sequence variations (SVs) were identified when comparing with Chinese Long (CL), and known SVs controlling spines, tubercles, and carpel number were confirmed in XIS49 genome. Two QTLs of hypocotyl elongation under low light, SH3.1 and SH6.1, were fine-mapped using introgression lines (donor parent, XIS49; recurrent parent, CL). SH3.1 encodes a red-light receptor Phytochrome B (PhyB, CsaV3_3G015190). A ∼4 kb region with large deletion and highly divergent regions (HDRs) were identified in the promoter of the PhyB gene in XIS49. Loss of function of this PhyB caused a super-long hypocotyl phenotype. SH6.1 encodes a CCCH-type zinc finger protein FRIGIDA-ESSENTIAL LIKE (FEL, CsaV3_6G050300). FEL negatively regulated hypocotyl elongation but it was transcriptionally suppressed by long terminal repeats retrotransposon insertion in CL cucumber. Mechanistically, FEL physically binds to the promoter of CONSTITUTIVE PHOTOMORPHOGENIC 1a (COP1a), regulating the expression of COP1a and the downstream hypocotyl elongation. These above results demonstrate the genetic mechanism of cucumber hypocotyl elongation under low light.

Genome-Wide Identification and Expression Analysis of the Cucumber FKBP Gene Family in Response to Abiotic and Biotic Stresses

The FKBP (FK506-binding protein) gene family is an important member of the PPlase protease family and plays a vital role during the processes of plant growth and development. However, no studies of the FKBP gene family have been reported in cucumber. In this study, 19 FKBP genes were identified in cucumber, which were located on chromosomes 1, 3, 4, 6, and 7. Phylogenetic analysis divided the cucumber FKBP genes into three subgroups. The FKBP genes in the same subgroup exhibited similar structures and conserved motifs. The cis-acting elements analysis revealed that the promoters of cucumber FKBP genes contained hormone-, stress-, and development-related cis-acting elements. Synteny analysis of the FKBP genes among cucumber, Arabidopsis, and rice showed that 12 kinds of syntenic relationships were detected between cucumber and Arabidopsis FKBP genes, and 3 kinds of syntenic relationships were observed between cucumber and rice FKBP genes. The tissue-specific expression analysis showed that some FKBP genes were expressed in all tissues, while others were only highly expressed in part of the 10 types of tissues. The expression profile analysis of cucumber FKBP genes under 13 types of stresses showed that the CsaV3_1G007080 gene was differentially expressed under abiotic stresses (high temperature, NaCl, silicon, and photoperiod) and biotic stresses (downy mildew, green mottle mosaic virus, Fusarium wilt, phytophthora capsica, angular leaf spot, and root-knot nematode), which indicated that the CsaV3_1G007080 gene plays an important role in the growth and development of cucumber. The interaction protein analysis showed that most of the proteins in the FKBP gene family interacted with each other. The results of this study will lay the foundation for further research on the molecular biological functions of the cucumber FKBP gene family.

Genomic analyses reveal dead-end hybridization between two deeply divergent kiwifruit species rather than homoploid hybrid speciation

Despite the importance of hybridization in evolution, the evolutionary consequence of homoploid hybridizations in plants remains poorly understood. Specially, homoploid hybridization events have been rarely documented due to a lack of genomic resources and methodological limitations. Actinidia zhejiangensis was suspected to have arisen from hybridization of Actinidia eriantha and Actinidia hemsleyana or Actinidia rufa. However, this species was very rare in nature and exhibited sympatric distribution with its potential parent species, which implied it might be a spontaneous hybrid of ongoing homoploid hybridization. Here, we illustrate the dead-end homoploid hybridization and genomic basis of isolating barriers between A. eriantha and A. hemsleyana through whole genome sequencing and population genomic analyses. Chromosome-scale genome assemblies of A. zhejiangensis and A. hemsleyana were generated. The chromosomes of A. zhejiangensis are confidently assigned to the two haplomes, and one of them originates from A. eriantha and the other originates from A. hemsleyana. Whole genome resequencing data reveal that A. zhejiangensis are mainly F1 hybrids of A. hemsleyana and A. eriantha and gene flow initiated about 0.98 million years ago, implying both strong genetic barriers and ongoing hybridization between these two deeply divergent kiwifruit species. Five inversions containing genes involved in pollen germination and pollen tube growth might account for the fertility breakdown of hybrids between A. hemsleyana and A. eriantha. Despite its distinct morphological traits and long recurrent hybrid origination, A. zhejiangensis does not initiate speciation. Collectively, our study provides new insights into homoploid hybridization in plants and provides genomic resources for evolutionary and functional genomic studies of kiwifruit.

Genome-Wide Analysis of Exocyst Complex Subunit Exo70 Gene Family in Cucumber

Cucumber (Cucumis sativus L.) is an important vegetable worldwide, but its yield is affected by a wide range of pathogens and pests. As the major subunit of the exocyst complex, the roles of Exo70 members have been shown in Arabidopsis and rice, but their function are unknown in cucumber. Here, we identified 18 CsExo70 members in cucumber, which were divided into three groups (Exo70.1-Exo70.3) and nine subgroups (Exo70A-Exo70I) based on the phylogenetic tree. Subsequently, systematical analyses were performed, including collinearity, gene structure, cis-acting elements, conserved motifs, expression patterns, and subcellular localization. Our results showed that CsExo70 genes were generally expressed in all tissues, and CsExo70C1 and CsExo70C2 were highly expressed in the stamen. Moreover, the expression levels of most CsExo70 genes were induced by Pseudomonas syringae pv. lachrymans (Psl) and Fusarium oxysporum f. sp. cucumerinum Owen (Foc), especially CsExo70E2 and CsExo70H3. In addition, these CsExo70s displayed similar location patterns with discrete and punctate signals in the cytoplasm. Together, our results indicate that CsExo70 members may be involved in plant development and resistance, and provide a reference for future in-depth studies of Exo70 genes in cucumber.

Understanding Transcription Factors and How They Affect Processes in Cucumber Sex Determination

Plant reproduction is a fundamental process on Earth from the perspective of biodiversity, biomass gain, and crop productivity. It is therefore important to understand the sex determination process, and many researchers are investigating the molecular basis of this phenomenon. However, information on the influence of transcription factors (TFs), genes that encode DNA-binding proteins, on this process is limited, although cucumber is a model plant in this regard. In the present study, based on RNA-seq data for differentially expressed genes (DEGs), we aimed to investigate the regulatory TFs that may influence the metabolic processes in the shoot apex containing the forming flower buds. Therefore, the annotation of the genome of the B10 cucumber line was supplemented with the assigned families of transcription factors. By performing ontology analyses of the DEGs, the processes they participate in were identified, and TFs were located among the results. In addition, TFs that have significantly overrepresented targets among DEGs were detected, and sex-specific interactome network maps were generated, indicating the regulatory TFs based on their effects on DEGs and furthermore, on the processes leading to the formation of different-sex flowers. Among the most overrepresented TF families in the sex comparisons were the NAC, bHLH, MYB, and bZIP families. An interaction network analysis indicated the most abundant families among DEGs' regulatory TFs were MYB, AP2/ERF, NAC, and bZIP, and those with the most significant impact on developmental processes were identified, namely the AP/ERF family, followed by DOF, MYB, MADS, and others. Thus, the networks' central nodes and key regulators were identified with respect to male, female, and hermaphrodite forms. Here, we proposed the first model of the regulatory network of TFs that influences the metabolism of sex development in cucumber. These findings may help us to understand the molecular genetics and functional mechanisms underlying sex determination processes.

Lateral Root Initiation in Cucumber (Cucumis sativus): What Does the Expression Pattern of Rapid Alkalinization Factor 34 (RALF34) Tell Us?

In Arabidopsis, the small signaling peptide (peptide hormone) RALF34 is involved in the gene regulatory network of lateral root initiation. In this study, we aimed to understand the nature of the signals induced by RALF34 in the non-model plant cucumber (Cucumis sativus), where lateral root primordia are induced in the apical meristem of the parental root. The RALF family members of cucumber were identified using phylogenetic analysis. The sequence of events involved in the initiation and development of lateral root primordia in cucumber was examined in detail. To elucidate the role of the small signaling peptide CsRALF34 and its receptor CsTHESEUS1 in the initial stages of lateral root formation in the parental root meristem in cucumber, we studied the expression patterns of both genes, as well as the localization and transport of the CsRALF34 peptide. CsRALF34 is expressed in all plant organs. CsRALF34 seems to differ from AtRALF34 in that its expression is not regulated by auxin. The expression of AtRALF34, as well as CsRALF34, is regulated in part by ethylene. CsTHESEUS1 is expressed constitutively in cucumber root tissues. Our data suggest that CsRALF34 acts in a non-cell-autonomous manner and is not involved in lateral root initiation in cucumber.

Insight into the Organization of the B10v3 Cucumber Genome by Integration of Biological and Bioinformatic Data

The availability of a well-organized and annotated reference genome is essential for genome research and the analysis of re-sequencing approaches. The B10v3 cucumber (Cucumis sativus L.) reference genome has been sequenced and assembled into 8035 contigs, a small fraction of which have been mapped to individual chromosomes. Currently, bioinformatics methods based on comparative homology have made it possible to re-order the sequenced contigs by mapping them to the reference genomes. The B10v3 genome (North-European, Borszczagowski line) was rearranged against the genomes of cucumber 9930 ('Chinese Long' line) and Gy14 (North American line). Furthermore, a better insight into the organization of the B10v3 genome was obtained by integrating the data available in the literature on the assignment of contigs to chromosomes in the B10v3 genome with the results of the bioinformatic analysis. The combination of information on the markers used in the assembly of the B10v3 genome and the results of FISH and DArT-seq experiments confirmed the reliability of the in silico assignment. Approximately 98% of the protein-coding genes within the chromosomes were assigned and a significant proportion of the repetitive fragments in the sequenced B10v3 genome were identified using the RagTag programme. In addition, BLAST analyses provided comparative information between the B10v3 genome and the 9930 and Gy14 data sets. This revealed both similarities and differences in the functional proteins found between the coding sequences region in the genomes. This study contributes to better knowledge and understanding of cucumber genome line B10v3.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Effects of Domestication on Plant-Microbiome Interactions

Through the process of domestication, selection is targeted on a limited number of plant traits that are typically associated with yield. As an unintended consequence, domesticated plants often perform poorly compared to their wild progenitors for a multitude of traits that were not under selection during domestication, including abiotic and biotic stress tolerance. Over the past decade, advances in sequencing technology have allowed for the rigorous characterization of host-associated microbial communities, termed the microbiome. It is now clear that nearly every conceivable plant interaction with the environment is mediated by interactions with the microbiome. For this reason, plant-microbiome interactions are an area of great promise for plant breeding and crop improvement. Here, we review the literature to assess the potential impact that domestication has had on plant-microbiome interactions and the current understanding of the genetic basis of microbiome variation to inform plant breeding efforts. Overall, we find limited evidence that domestication impacts the diversity of microbiomes, but domestication is often associated with shifts in the abundance and composition of microbial communities, including taxa of known functional significance. Moreover, genome-wide association studies and mutant analysis have not revealed a consistent set of core candidate genes or genetic pathways that confer variation in microbiomes across systems. However, such studies do implicate a consistent role for plant immunity, root traits, root and leaf exudates and cell wall integrity as key traits that control microbiome colonization and assembly. Therefore, selection on these key traits may pose the most immediate promise for enhancing plant-microbiome interactions through breeding.

Karyotyping of commercial cultivars of melon (Cucumis melo L.)

BACKGROUND

This study on cultivars of melon (Cucumis melo L.) marketed in Brazil was conducted to obtain information to be used in breeding programs of this species. Little is known about the karyotype variability among C. melo L. cultivars targeted at the consumer market. The objective of the present study was to verify the karyotype variability in eight commercial melon cultivars used in the Brazilian market.

METHODS AND RESULTS

Slides were stained with 2% Giemsa and assembled with Neomount to perform chromosomal morphometry. GC-rich heterochromatin was observed by CMA₃/DAPI staining. 5 S rDNA, centromeric satellite DNA (SatDNA), and telomeric sites were visualized using fluorescence in situ hybridization. All images were captured on an Olympus BX41 microscope equipped with a 5 M Olympus DP25 digital camera and DP2-BSW software. The cultivars showed symmetrical karyotypes with significant differences in total chromosome length and average chromosome size. Heterochromatic CMA₃⁺ blocks were observed in terminal regions related to satellites (secondary constrictions), as well as in centromeric and pericentromeric regions. A single chromosomal pair of 5 S rDNA sites was observed in all cultivars, but at distinct locations. Centromeric satellite sequences, tested for the first time in melon, revealed only centromeric sites. Telomeric sites were observed in all the chromosomes of the cultivars.

CONCLUSIONS

Karyotype variation was observed in cultivars of melon, which were analyzed for chromosomal morphology and localization of GC-rich heterochromatin, as well centromeric SatDNA, rDNA, and telomeric chromosomal markers. Hence, these cultivars can be used in future breeding programs.

Genome-wide identification and expression analysis of the cucumber PP2C gene family

Background

Type 2C protein phosphatase (PP2C) is a negative regulator of ABA signaling pathway, which plays important roles in stress signal transduction in plants. However, little research on the PP2C genes family of cucumber (Cucumis sativus L.), as an important economic vegetable, has been conducted.

Results

This study conducted a genome-wide investigation of the CsPP2C gene family. Through bioinformatics analysis, 56 CsPP2C genes were identified in cucumber. Based on phylogenetic analysis, the PP2C genes of cucumber and Arabidopsis were divided into 13 groups. Gene structure and conserved motif analysis showed that CsPP2C genes in the same group had similar gene structure and conserved domains. Collinearity analysis showed that segmental duplication events played a key role in the expansion of the cucumber PP2C genes family. In addition, the expression of CsPP2Cs under different abiotic treatments was analyzed by qRT-PCR. The results reveal that CsPP2C family genes showed different expression patterns under ABA, drought, salt, and cold treatment, and that CsPP2C3, 11–17, 23, 45, 54 and 55 responded significantly to the four stresses. By predicting the cis-elements in the promoter, we found that all CsPP2C members contained ABA response elements and drought response elements. Additionally, the expression patterns of CsPP2C genes were specific in different tissues.

Conclusions

The results of this study provide a reference for the genome-wide identification of the PP2C gene family in other species and provide a basis for future studies on the function of PP2C genes in cucumber.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08734-y.

Modified photoperiod response of CsFT promotes day neutrality and early flowering in cultivated cucumber

Map-based cloning and photoperiod response detection suggested that CsFT is the critical gene for cucumber photoperiod domestication. Photoperiod sensitivity is important for sensing seasonal changes and local adaptation. However, day-length sensitivity limits crop geographical adaptation and it should be modified during domestication. Cucumber was domesticated in southern Asia and is currently cultivated worldwide across a wide range of latitudes, but its photoperiod sensitivity and its change during cucumber domestication are unknown. Here, we confirmed wild cucumber (Hardwickii) was a short-day plant, and its flowering depends on short-day (SD) conditions, while the cultivated cucumber (9930) is a day-neutral plant that flowers independently of day length. A photoperiod sensitivity locus (ps-1) was identified by the 9930 × Hardwickii F₂ segregating populations, which span a ~ 970 kb region and contain 60 predicted genes. RNA-seq analysis showed that the critical photoperiod pathway gene FLOWERING LOCUS T (CsFT) within the ps-1 locus exhibits differential expression between 9930 and Hardwickii, which was confirmed by qRT-PCR detection. CsFT in Hardwickii was sensitive to day length and could be significantly induced by SD conditions, whereas CsFT was highly expressed in 9930 and was insensitive to day length. Moreover, the role of CsFT in promoting flowering was verified by overexpression of CsFT in Arabidopsis. We also identified the genetic variations existing in the promoter of CsFT among the different geographic cucumbers and suggest they have possible roles in photoperiod domestication. The results of this study suggest that a variation in photoperiod sensitivity of CsFT is associated with day neutrality and early flowering in cultivated cucumber and could contribute to cucumber cultivation in diverse regions throughout the world.

Recent status of Genotyping by Sequencing (GBS) Technology in cucumber (Cucumis sativus L.): a review

Current and advanced breeding tools are being used to improve economically important horticultural crops to meet the consumers' needs and preferences. Genotyping-by-sequencing (GBS) is an extremely useful tool in the investigation and analysis of the genetic diversity of different cultivars. Based on a broad range of genetic backgrounds like single nucleotide polymorphism (SNPs), GBS is known as a novel technique to facilitate the detection of quantitative trait loci (QTL) regions robustly linked with interested traits compared to genome-wide association study (GWAS) and QTL. GBS has gained popularity among breeders in recent years and it is also employed in cucumber breeding programs. Cucumbers (C. sativus L.) are monoecious, gynoecious and some of them are parthenocarpic species. Cucumber is one of the most economical and essential crops in the Cucurbitaceae family. For time immemorial, cucumber has been produced and consumed all over the world like other cucurbits. To a large extent, cultivated cucurbits are beneficial to human health for providing necessary minerals and fibers.Therefore, this review portrays the current status of advances made by using GBS and its combination with other tools in various studies of cucumber such as the use of GBS and single nucleotide polymorphism (SNP) markers, GBS and GWAS, also with QTL and marker-assisted selection (MAS) are applied to display and detect explicit genetic architecture complex traits in crops and chromosome rearrangements.Cucumber breeding programs have undoubtedly benefited from genotyping-by-sequencing. Using the GBS method, research discovered lots of new candidate genes that control various traits including spine color, fruit stalk-end color, and disease resistance in cucumber lines.

A critical review on cytogenetics of Cucurbitaceae with updates on Indian taxa

The cytogenetic relationships in the species of Cucurbitaceae are becoming immensely important to answer questions pertaining to genome evolution. Here, a simplified and updated data resource on cytogenetics of Cucurbitaceae is presented on the basis of foundational parameters (basic, zygotic and gametic chromosome numbers, ploidy, genome size, karyotype) and molecular cytogenetics. We have revised and collated our own findings on seven agriculturally important Indian cucurbit species in a comparative account with the globally published reports. Chromosome count (of around 19% species) shows nearly three-fold differences while genome size (of nearly 5% species) shows 5.84-fold differences across the species. There is no significant correlation between chromosome numbers and nuclear genome sizes. The possible trend of evolution is discussed here based on molecular cytogenetics data, especially the types and distribution of nucleolus organizer regions (NORs). The review supersedes the scopes of general chromosome databases and invites scopes for continuous updates. The offline resource serves as an exclusive toolkit for research and breeding communities across the globe and also opens scope for future establishment of web-database on Cucurbitaceae cytogenetics.

Comparative Analysis of Transposable Elements and the Identification of Candidate Centromeric Elements in the Prunus Subgenus Cerasus and Its Relatives

The subgenus Cerasus and its relatives include many crucial economic drupe fruits and ornamental plants. Repetitive elements make up a large part of complex genomes, and some of them play an important role in gene regulation that can affect phenotypic variation. However, the variation in their genomes remains poorly understood. This work conducted a comprehensive repetitive sequence identification across the draft genomes of eight taxa of the genus Prunus, including four of the Prunus subgenus Cerasus (Prunus pseudocerasus, P. avium, P. yedoensis and P. × yedoensis) as well as congeneric species (Prunus salicina, P. armeniaca, P. dulcis and P. persica). Annotation results showed high proportions of transposable elements in their genomes, ranging from 52.28% (P. armeniaca) to 61.86% (P. pseudocerasus). The most notable differences in the contents of long terminal repeat retrotransposons (LTR-RTs) and tandem repeats (TRs) were confirmed with de novo identification based on the structure of each genome, which significantly contributed to their genome size variation, especially in P. avium and P.salicina. Sequence comparisons showed many similar LTR-RTs closely related to their phylogenetic relationships, and a highly similar monomer unit of the TR sequence was conserved among species. Additionally, the predicted centromere-associated sequence was located in centromeric regions with FISH in the 12 taxa of Prunus. It presented significantly different signal intensities, even within the diverse interindividual phenotypes for Prunus tomentosa. This study provides insight into the LTR-RT and TR variation within Prunus and increases our knowledge about its role in genome evolution.

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.

Variation of Chromosome Composition in a Full-Sib Population Derived From 2x × 3x Interploidy Cross of Populus

Interploidy cross commonly results in complex chromosome number and structural variations. In our previous study, a progeny with segregated ploidy levels was produced by an interploidy cross between diploid female parent Populus tomentosa × Populus bolleana clone TB03 and triploid male parent Populus alba × Populus berolinensis 'Yinzhong'. However, the chromosome compositions of aneuploid genotypes in the progeny were still unclear. In the present study, a microsatellite DNA allele counting-peak ratios (MAC-PR) method was employed to analyze allelic configurations of each genotype to clarify their chromosome compositions, while 45S rDNA fluorescence in situ hybridization (FISH) analysis was used to reveal the mechanism of chromosome number variation. Based on the MAC-PR analysis of 47 polymorphic simple sequence repeat (SSR) markers distributed across all 19 chromosomes of Populus, both chromosomal number and structural variations were detected for the progeny. In the progeny, 26 hypo-triploids, 1 hyper-triploid, 16 hypo-tetraploids, 10 tetraploids, and 5 hyper-tetraploids were found. A total of 13 putative structural variation events (duplications and/or deletions) were detected in 12 genotypes, involved in chromosomes 3, 6, 7, 14, 15, 16, and 18. The 46.2% (six events) structural variation events occurred on chromosome 6, suggesting that there probably is a chromosome breakpoint near the SSR loci of chromosome 6. Based on calculation of the allelic information, the transmission of paternal heterozygosity in the hypo-triploids, hyper-triploid, hypo-tetraploids, tetraploids, and hyper-tetraploids were 0.748, 0.887, 0.830, 0.833, and 0.836, respectively, indicating that the viable pollen gains of the male parent 'Yinzhong' were able to transmit high heterozygosity to progeny. Furthermore, 45S rDNA-FISH analysis showed that specific-chromosome segregation feature during meiosis and chromosome appointment in normal and fused daughter nuclei of telophase II of 'Yinzhong,' which explained that the formation of aneuploids and tetraploids in the progeny could be attributed to imbalanced meiotic chromosomal segregation and division restitution of 'Yinzhong,' The data of chromosomal composition and structural variation of each aneuploid in the full-sib progeny of TB03 × 'Yinzhong' lays a foundation for analyzing mechanisms of trait variation relying on chromosome or gene dosages in Populus.

Genome-wide identification and expression analysis of the cucumber PYL gene family

Abscisic acid (ABA) is a very important hormone in plants. It regulates growth and development of plants and plays an important role in biotic and abiotic stresses. The Pyrabactin resistance 1-like (PYR/PYL) proteins play a central role in ABA signal transduction pathways. The working system of PYL genes in cucumber, an important economical vegetable (Cucumis sativus L.), has not been fully studied yet. Through bioinformatics, a total of 14 individual PYL genes were identified in Chinese long '9930' cucumber. Fourteen PYL genes were distributed on six chromosomes of cucumber, and their encoded proteins predicted to be distributed in cytoplasm and nucleus. Based on the phylogenetic analysis, the PYL genes of cucumber, Arabidopsis, rice, apple, Brachypodium distachyon and soybeancould be classified into three groups. Genetic structures and conserved domains analysis revealed that CsPYL genes in the same group have similar exons and conserved domains. By predicting cis-elements in the promoters, we found that all CsPYL members contained hormone and stress-related elements. Additionally, the expression patterns of CsPYL genes were specific in tissues. Finally, we further examined the expression of 14 CsPYL genes under ABA, PEG, salt stress. The qRT-PCR results showed that most PYL gene expression levels were up-regulated. Furthermore, with different treatments about 3h, the relative expression of PYL8 was up-regulated and more than 20 times higher than 0h. It indicated that this gene may play an important role in abiotic stress.

Fine mapping a quantitative trait locus underlying seedling resistance to gummy stem blight using a residual heterozygous lines-derived strategy in cucumber

Gummy stem blight (GSB), caused by Didymella bryoniae, is one of the most devastating diseases that severely reduces cucumber production. Developing resistant varieties would be an effective strategy to control GSB. Although several GSB-resistant QTLs have been reported, causal genes for GSB resistance have not yet been identified in cucumber. A novel loci gsb3.1 for seedling GSB resistance from the "PI 183967" genotype was previously identified in a 1.7-Mb interval on chromosome 3. In this study, we developed a residual heterozygous line-derived strategy from Recombinant Inbred Lines to perform fine mapping, and with this approach, the gsb3.1 locus was narrowed to a 38 kb interval. There were six predicted genes at the gsb3.1 locus, four of which differed in expression in the GSB-resistant compared to the susceptible lines after fungal inoculation. These candidate genes (Csa3G020050, Csa3G020060, Csa3G020090, and Csa3G020590) within the gsb3.1 locus could be helpful for the genetic study of GSB resistance and marker-assisted selection in cucumber. Phylogenetic analyses indicated that the resistant gsb3.1 allele may uniquely exist in the wild species present in the Indian group, and that nucleotide diversity was significantly reduced in cultivated accessions. Therefore, the gsb3.1 allele could be introgressed into existing commercial cultivars and combined with other resistance QTLs to provide broad-spectrum and robust GSB resistance in cucumber.

Identification of an Allelic Variant of the CsOr Gene Controlling Fruit Endocarp Color in Cucumber (Cucumis sativus L.) Using Genotyping-By-Sequencing (GBS) and Whole-Genome Sequencing

The cucumber is a major vegetable crop around the world. Fruit flesh color is an important quality trait in cucumber and flesh color mainly depends on the relative content of β-carotene in the fruits. The β-carotene serves as a precursor of vitamin A, which has dietary benefits for human health. Cucumbers with orange flesh contain a higher amount of β-carotene than white fruit flesh. Therefore, development of orange-fleshed cucumber varieties is gaining attention for improved nutritional benefits. In this study, we performed genotyping-by-sequencing (GBS) based on genetic mapping and whole-genome sequencing to identify the orange endocarp color gene in the cucumber breeding line, CS-B. Genetic mapping, genetic sequencing, and genetic segregation analyses showed that a single recessive gene (CsaV3_6G040750) encodes a chaperone DnaJ protein (DnaJ) protein at the Cucumis sativus(CsOr) locus was responsible for the orange endocarp phenotype in the CS-B line. The Or gene harbored point mutations T13G and T17C in the first exon of the coding region, resulting in serine to alanine at position 13 and isoleucine to threonine at position 17, respectively. CS-B line displayed increased β-carotene content in the endocarp tissue, corresponding to elevated expression of CsOr gene at fruit developmental stages. Identifying novel missense mutations in the CsOr gene could provide new insights into the role of Or mechanism of action for orange fruit flesh in cucumber and serve as a valuable resource for developing β-carotene-rich cucumbers varieties with increased nutritional benefits.

Molecular research progress and improvement approach of fruit quality traits in cucumber

Recent molecular studies revealed new opportunities to improve cucumber fruit quality. However, the fruit color and spine traits molecular basis remain vague despite the vast sources of genetic diversity. Cucumber is agriculturally, economically and nutritionally important vegetable crop. China produces three-fourths of the world's total cucumber production. Cucumber fruit quality depends on a number of traits such as the fruit color (peel and flesh color), spine (density, size and color), fruit shape, fruit size, defects, texture, firmness, taste, maturity stage and nutritional composition. Fruit color and spine traits determine critical quality attributes and have been the interest of researchers at the molecular level. Evaluating the molecular mechanisms of fruit quality traits is important to improve production and quality of cucumber varieties. Genes and qualitative trait locus (QTL) that are responsible for cucumber fruit color and fruit spine have been identified. The purpose of this paper is to reveal the molecular research progress of fruit color and spines as key quality traits of cucumber. The markers and genes identified so far could help for marker-assisted selection of the fruit color and spine trait in cucumber breeding and its associated nutritional improvement. Based on the previous studies, peel color and spine density as examples, we proposed a comprehensive approach for cucumber fruit quality traits improvement. Moreover, the markers and genes can be useful to facilitate cloning-mediated genetic breeding in cucumber. However, in the era of climate change, increased human population and high-quality demand of consumers, studies on molecular mechanisms of cucumber fruit quality traits are limited.

A near-complete genome sequence of mungbean (Vigna radiata L.) provides key insights into the modern breeding program

Mungbean (Vigna radiata L.), a fast-growing legume species, is an important source of carbohydrates and proteins in developing countries of Asia. Here, we constructed a near-complete genome sequence of mungbean with a scaffold N50 value of 5.2 Mb and only a 0.4% gap, with a total scaffold size of 475 Mb. We identified several misassembled pseudomolecules (Chr03, Chr04, Chr05, and Chr08) in the previous draft assembly; Chr03, Chr04, and Chr08 were assembled into one chromosome, and Chr05 was broken into two chromosomes in the improved reference genome assembly, thus providing more accurate linkage information to breeders. Additionally, using an ultra-high-resolution linkage map constructed based on resequencing data, we identified several quantitative trait loci (QTLs) and the underlying candidate genes affecting synchronous pod maturity (SPM). Mungbean homologs of two soybean ([Glycine max (L.) Merr.] flowering genes, E3 (phytochrome A) and J (early flowering 3), were identified as candidate genes for the QTLs, and the candidate genes for plant height, node number, and SPM showed critical nucleotide substitutions between the reference cultivar and other genotypes (landraces and wild accessions). Based on the analysis of genetic diversity among 276 accessions collected from 23 countries, we identified 36 selective sweep regions and observed that the overall genetic diversity of cultivars decreased to 30% of that in wild accessions postdomestication. The near-complete genome sequence of mungbean represents an important resource for genome-assisted improvement in the mungbean breeding program.

Development of SNP Markers for White Immature Fruit Skin Color in Cucumber (Cucumis sativus L.) Using QTL-seq and Marker Analyses

Despite various efforts in identifying the genes governing the white immature fruit skin color in cucumber, the genetic basis of the white immature fruit skin color is not well known. In the present study, genetic analysis showed that a recessive gene confers the white immature fruit skin-color phenotype over the light-green color of a Korean slicer cucumber. High-throughput QTL-seq combined with bulked segregation analysis of two pools with the extreme phenotypes (white and light-green fruit skin color) in an F₂ population identified two significant genomic regions harboring QTLs for white fruit skin color within the genomic region between 34.1 and 41.67 Mb on chromosome 3, and the genomic region between 12.2 and 12.7 Mb on chromosome 5. Further, nonsynonymous SNPs were identified with a significance of p < 0.05 within the QTL regions, resulting in eight homozygous variants within the QTL region on chromosome 3. SNP marker analysis uncovered the novel missense mutations in Chr3CG52930 and Chr3CG53640 genes and showed consistent results with the phenotype of light-green and white fruit skin-colored F₂ plants. These two genes were located 0.5 Mb apart on chromosome 3, which are considered strong candidate genes. Altogether, this study laid a solid foundation for understanding the genetic basis and marker-assisted breeding of immature fruit skin color in cucumber.

Genomic Prediction and the Practical Breeding of 12 Quantitative-Inherited Traits in Cucumber (Cucumis sativus L.)

Genomic prediction is an effective way for predicting complex traits, and it is becoming more essential in horticultural crop breeding. In this study, we applied genomic prediction in the breeding of cucumber plants. Eighty-one cucumber inbred lines were genotyped and 16,662 markers were identified to represent the genetic background of cucumber. Two populations, namely, diallel cross population and North Carolina II population, having 268 combinations in total were constructed from 81 inbred lines. Twelve cucumber commercial traits of these two populations in autumn 2018, spring 2019, and spring 2020 were collected for model training. General combining ability (GCA) models under five-fold cross-validation and cross-population validation were applied to model validation. Finally, the GCA performance of 81 inbred lines was estimated. Our results showed that the predictive ability for 12 traits ranged from 0.38 to 0.95 under the cross-validation strategy and ranged from -0.38 to 0.88 under the cross-population strategy. Besides, GCA models containing non-additive effects had significantly better performance than the pure additive GCA model for most of the investigated traits. Furthermore, there were a relatively higher proportion of additive-by-additive genetic variance components estimated by the full GCA model, especially for lower heritability traits, but the proportion of dominant genetic variance components was relatively small and stable. Our findings concluded that a genomic prediction protocol based on the GCA model theoretical framework can be applied to cucumber breeding, and it can also provide a reference for the single-cross breeding system of other crops.

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.

Reconstruction of ancestral karyotype illuminates chromosome evolution in the genus Cucumis

Karyotype dynamics driven by complex chromosome rearrangements constitute a fundamental issue in evolutionary genetics. The evolutionary events underlying karyotype diversity within plant genera, however, have rarely been reconstructed from a computed ancestral progenitor. Here, we developed a method to rapidly and accurately represent extant karyotypes with the genus, Cucumis, using highly customizable comparative oligo-painting (COP) allowing visualization of fine-scale genome structures of eight Cucumis species from both African-origin and Asian-origin clades. Based on COP data, an evolutionary framework containing a genus-level ancestral karyotype was reconstructed, allowing elucidation of the evolutionary events that account for the origin of these diverse genomes within Cucumis. Our results characterize the cryptic rearrangement hotspots on ancestral chromosomes, and demonstrate that the ancestral Cucumis karyotype (n = 12) evolved to extant Cucumis genomes by hybridizations and frequent lineage- and species-specific genome reshuffling. Relative to the African species, the Asian species, including melon (Cucumis melo, n = 12), Cucumis hystrix (n = 12) and cucumber (Cucumis sativus, n = 7), had highly shuffled genomes caused by large-scale inversions, centromere repositioning and chromothripsis-like rearrangement. The deduced reconstructed ancestral karyotype for the genus allowed us to propose evolutionary trajectories and specific events underlying the origin of these Cucumis species. Our findings highlight that the partitioned evolutionary plasticity of Cucumis karyotype is primarily located in the centromere-proximal regions marked by rearrangement hotspots, which can potentially serve as a reservoir for chromosome evolution due to their fragility.

Genetic regulation of shoot architecture in cucumber

Cucumber (Cucumis sativus L.) is an important vegetable crop species with great economic value. Shoot architecture determines the visual appearance of plants and has a strong impact on crop management and yield. Unlike most model plant species, cucumber undergoes vegetative growth and reproductive growth simultaneously, in which leaves are produced from the shoot apical meristem and flowers are generated from leaf axils, during the majority of its life, a feature representative of the Cucurbitaceae family. Despite substantial advances achieved in understanding the regulation of plant form in Arabidopsis thaliana, rice, and maize, our understanding of the mechanisms controlling shoot architecture in Cucurbitaceae crop species is still limited. In this review, we focus on recent progress on elucidating the genetic regulatory pathways underlying the determinant/indeterminant growth habit, leaf shape, branch outgrowth, tendril identity, and vine length determination in cucumber. We also discuss the potential of applying biotechnology tools and resources for the generation of ideal plant types with desired architectural features to improve cucumber productivity and cultivation efficiency.

Domestication of Crop Metabolomes: Desired and Unintended Consequences

The majority of the crops and vegetables of today were domesticated from their wild progenitors within the past 12 000 years. Considerable research effort has been expended on characterizing the genes undergoing positive and negative selection during the processes of crop domestication and improvement. Many studies have also documented how the contents of a handful of metabolites have been altered during human selection, but we are only beginning to unravel the true extent of the metabolic consequences of breeding. We highlight how crop metabolomes have been wittingly or unwittingly shaped by the processes of domestication, and highlight how we can identify new targets for metabolite engineering for the purpose of de novo domestication of crop wild relatives.

Polymorphism in the Chloroplast ATP Synthase Beta-Subunit Is Associated with a Maternally Inherited Enhanced Cold Recovery in Cucumber

Cucumber (Cucumis sativus L.) is a warm-season crop that is sensitive to chilling temperatures and a maternally inherited cold tolerance exists in the heirloom cultivar 'Chipper' (CH). Because the organelles of cucumber show differential transmission (maternal for chloroplast and paternal for mitochondrion), this cold tolerance is hypothesized to be chloroplast-associated. The goal of this research was to characterize the cold tolerant phenotype from CH and determine its genetic basis. Doubled haploid (DH) lines were produced from CH and cold susceptible cucumbers, reciprocal hybrids with identical nuclear genotypes were produced, and plants were subjected to cold treatments under lights at 4 °C for 5.5 h. Hybrid plants with CH as the maternal parent had significantly higher fresh and dry weights 14 days after cold treatment compared to the reciprocal hybrid, revealing an enhanced cold recovery phenotype maternally conferred by CH. Results from analyses of the nuclear transcriptome and reactive oxygen species (ROS) between reciprocal hybrids were consistent with the cold recovery phenotype. Sequencing of the chloroplast genome and transcriptome of the DH parents and reciprocal hybrids, respectively, revealed one maternally transmitted non-synonymous single nucleotide polymorphism (SNP) in the chloroplast F₁F_O-ATP synthase (CF₁F_O-ATPase) beta-subunit gene (atpB) of CH which confers an amino acid change from threonine to arginine. Protein modeling revealed that this change is located at the interface of the alpha- and beta-subunits in the CF₁F_O-ATPase complex. Polymorphisms in the CF₁F_O-ATPase complex have been associated with stress tolerances in other plants, and selection for or creation of polymorphic beta-subunit proteins by chloroplast transformation or gene editing could condition improved recovery from cold stress in plants.

A SNP of HD-ZIP I transcription factor leads to distortion of trichome morphology in cucumber (Cucumis sativus L.)

BACKGROUND

Trichomes are excellent model systems for the analysis of cell differentiation and play essential roles in plant protection. From cucumber inbred line 'WD1', we identified an EMS-induced trichome abnormally developing mutant, nps, which exhibited smaller, denser and no pyramid-shaped head trichomes.

RESULTS

Using F₂ and BC₁ populations constructed from a cross between nps and '9930', the genetic analysis showed that the nps trait is controlled by a single recessive nuclear gene. We identified CsNps by map-based cloning with 576 individuals of the F₂ population generated from the cross of nps and inbred line '9930'. The CsNps was located at a 13.4-kb genomic region on chromosome 3, which region contains three predicted genes. Sequence analysis showed that only one single nucleotide mutation (C → T) between 9930 and nps was found in the second exon of Csa3G748220, a plant-specific class I HD-Zip gene. The result of allelism test also indicated that nps is a novel allelic mutant of Mict (Micro-trichome). Thus, nps was renamed mict-L130F. By comparing the transcriptome of mict-L130F vs WD1 and 06-2 (mict) vs 06-1 (wildtype, near-isogenic line of 06-2), several potential target genes that may be related to trichome development were identified.

CONCLUSIONS

Our results demonstrate that Mict-L130F is involved in the morphogenesis of trichomes. Map-based cloning of the Mict-L130F gene could promote the study of trichome development in cucumber.

Agarose-resolvable InDel markers based on whole genome re-sequencing in cucumber

Insertion and Deletion (InDel) are common features in genomes and are associated with genetic variation. The whole-genome re-sequencing data from two parents (X1 and X2) of the elite cucumber (Cucumis sativus) hybrid variety Lvmei No.1 was used for genome-wide InDel polymorphisms analysis. Obtained sequence reads were mapped to the genome reference sequence of Chinese fresh market type inbred line ‘9930’ and gaps conforming to InDel were pinpointed. Further, the level of cross-parents polymorphism among five pairs of cucumber breeding parents and their corresponding hybrid varieties were used for evaluating hybrid seeds purity test efficiency of InDel markers. A panel of 48 cucumber breeding lines was utilized for PCR amplification versatility and phylogenetic analysis of these markers. In total, 10,470 candidate InDel markers were identified for X1 and X2. Among these, 385 markers with more than 30 nucleotide difference were arbitrary chosen. These markers were selected for experimental resolvability through electrophoresis on an Agarose gel. Two hundred and eleven (211) accounting for 54.81% of markers could be validated as single and clear polymorphic pattern while 174 (45.19%) showed unclear or monomorphic genetic bands between X1 and X2. Cross-parents polymorphism evaluation recorded 68 (32.23%) of these markers, which were designated as cross-parents transferable (CPT) InDel markers. Interestingly, the marker InDel114 presented experimental transferability between cucumber and melon. A panel of 48 cucumber breeding lines including parents of Lvmei No. 1 subjected to PCR amplification versatility using CPT InDel markers successfully clustered them into fruit and common cucumber varieties based on phylogenetic analysis. It is worth noting that 16 of these markers were predominately associated to enzymatic activities in cucumber. These agarose-based InDel markers could constitute a valuable resource for hybrid seeds purity testing, germplasm classification and marker-assisted breeding in cucumber.

QTL-seq analysis of powdery mildew resistance in a Korean cucumber inbred line

QTL mapping and RT-PCR analyses identified the CsGy5G015660 as a strong powdery mildew resistance candidate gene and natural variation of CsGy5G015660 allele was observed using 115 core germplasm. Powdery mildew (PM) is among the most serious fungal diseases encountered in the cultivation of cucurbits. The development of PM-resistant inbred lines is thus of considerable significance for cucumber breeding programs. In this study, we applied bulked segregant analysis combined with QTL-seq to identify PM resistance loci using F₂ population derived from a cross between two Korean cucumber inbred lines, PM-R (resistant) and PM-S (susceptible). Genome-wide SNP profiling using bulks of the two extreme phenotypes identified two QTLs on chromosomes 5 and 6, designated pm5.2 and pm6.1, respectively. The two PM resistance loci were validated using molecular marker-based classical QTL analysis: pm5.2 (30% R² at LOD 11) and pm6.1 (11% R² at LOD 3.2). Furthermore, reverse transcriptase-PCR analyses, using genes found to be polymorphic between PM-R and PM-S, were conducted to identify the candidate gene(s) responsible for PM resistance. We found that transcripts of the gene CsGy5G015660, encoding a putative leucine-rich repeat receptor-like serine/threonine-protein kinase (RPK2), showed specific accumulation in PM-R prior to the appearance of disease symptoms, and was accordingly considered a strong candidate gene for PM resistance. In addition, cleaved amplified polymorphic sequence markers from CsGy5G015660 were developed and used to screen 35 inbred lines. Natural variation in the CsGy5G015660 allele was also observed based on analysis of a core collection of 115 cucumber accessions. Our results provide new genetic insights for gaining a better understanding of the genetic basis of PM resistance in cucumber, and pave the way for further utilization in cucumber PM resistance breeding programs.

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.

The Formation of Fruit Quality in Cucumis sativus L

Cucumber is one of the most widely grown vegetables in China and an indispensable fresh fruit in the diet. With the development of society, the demand of people for cucumber quality is higher and higher. Therefore, cultivating high-quality cucumber varieties is one of the main goals of cucumber breeding. With the rapid development of biotechnology such as molecular marker, cucumber quality control network is becoming clear. In this review, we describe the formation mechanism of cucumber fruit quality from three aspects: (1) the commercial quality of cucumber fruit, (2) nutritional quality formation, and (3) flavor quality of cucumber fruit. In addition, the determinants of cucumber fruit quality were summarized from two aspects of genetic regulation and cultivation methods in order to provide ideas for cucumber researchers and cultivators to improve fruit quality.

The genome and transcriptome analysis of snake gourd provide insights into its evolution and fruit development and ripening

Snake gourd (Trichosanthes anguina L.), which belongs to the Cucurbitaceae family, is a popular ornamental and food crop species with medicinal value and is grown in many parts of the world. Although progress has been made in its genetic improvement, the organization, composition, and evolution of the snake gourd genome remain largely unknown. Here, we report a high-quality genome assembly for snake gourd, comprising 202 contigs, with a total size of 919.8 Mb and an N50 size of 20.1 Mb. These findings indicate that snake gourd has one of the largest genomes of Cucurbitaceae species sequenced to date. The snake gourd genome assembly harbors 22,874 protein-coding genes and 80.0% of the genome consists of repetitive sequences. Phylogenetic analysis reveals that snake gourd is closely related to sponge gourd but diverged from their common ancestor ~33-47 million years ago. The genome sequence reported here serves as a valuable resource for snake gourd genetic research and comparative genomic studies in Cucurbitaceae and other plant species. In addition, fruit transcriptome analysis reveals the candidate genes related to quality traits during snake gourd fruit development and provides a basis for future research on snake gourd fruit development and ripening at the transcript level.

Comparative genomics of muskmelon reveals a potential role for retrotransposons in the modification of gene expression

Melon exhibits substantial natural variation especially in fruit ripening physiology, including both climacteric (ethylene-producing) and non-climacteric types. However, genomic mechanisms underlying such variation are not yet fully understood. Here, we report an Oxford Nanopore-based high-grade genome reference in the semi-climacteric cultivar Harukei-3 (378 Mb + 33,829 protein-coding genes), with an update of tissue-wide RNA-seq atlas in the Melonet-DB database. Comparison between Harukei-3 and DHL92, the first published melon genome, enabled identification of 24,758 one-to-one orthologue gene pairs, whereas others were candidates of copy number variation or presence/absence polymorphisms (PAPs). Further comparison based on 10 melon genome assemblies identified genome-wide PAPs of 415 retrotransposon Gag-like sequences. Of these, 160 showed fruit ripening-inducible expression, with 59.4% of the neighboring genes showing similar expression patterns (r > 0.8). Our results suggest that retrotransposons contributed to the modification of gene expression during diversification of melon genomes, and may affect fruit ripening-inducible gene expression.

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

QTL for horticulturally important traits associated with pleiotropic andromonoecy and carpel number loci, and a paracentric inversion in cucumber

The legendary cucumber inbred line WI2757 possesses a rare combination of resistances against nine pathogens, which is an important germplasm for cucumber breeding. However, WI2757 flowers late and does not perform well under field conditions. The genetic basis for horticulturally important traits other than disease resistances in WI2757 is largely unknown. In this study, we conducted QTL mapping using F₂ and recombinant inbred line (RIL) populations from the WI2757 × True Lemon cross that were segregating for multiple traits. Phenotypic data were collected in replicated field trials across multiple years for seven traits including fruit carpel number (CN) and sex expression. A high-density SNP-based genetic map was developed with genotyping by sequencing of the RIL population, which revealed a region on chromosome 1 with strong recombination suppression. The reduced recombination in this region was due to a ~ 10-Mbp paracentric inversion in WI2757 that was confirmed with additional segregation and cytological (FISH) analyses. Thirty-six QTL were detected for flowering time, fruit length (FL), fruit diameter (FD), fruit shape (LD), fruit number (FN), CN, and powdery mildew resistance. Five moderate- or major-effect QTL for FL, FD, LD, and FN inside the inversion are likely the pleiotropic effects of the andromonoecy (m), or the cn locus. The major-effect flowering time QTL ft1.1 was also mapped inside the inversion, which seems to be different from the previously assigned delayed flowering in WI2757. Implications of these findings on the use of WI2757 in cucumber breeding are discussed.

Cucumber Fruit Size and Shape Variations Explored from the Aspects of Morphology, Histology, and Endogenous Hormones

Fruit size and shape are important qualities and yield traits in cucumber (Cucumis sativus L.), but the factors that influence fruit size and shape remain to be explored. In this study, we investigated the dynamic changes of fruit size and shape from the aspects of morphology, cellular levels and endogenous hormones for nine typical cucumber inbred lines. The results show that fruit length had a strong positive correlation to the cell number in the longitudinal section of fruit throughout the four stages of 0, 6, 12, and 30 DAA (days after anthesis). However, the significant negative correlations were found between fruit length and the fruit cell size at 12 and 30 DAA. Furthermore, fruit diameter was positively correlated to the cell number in the cross section at all the investigated fruit growth stages. The indole-3-acetic acid (IAA) content showed significant positive correlations to the fruit length at all fruit growth stages of −6, −3, 0, 3, 6, 9 and 12 DAA, but IAA content and fruit diameter showed significant negative correlations for all the stages except for at −6 DAA. The trans-zeatin riboside (tZR), zeatin (ZT), gibberellic acid (GA₃) and jasmonic acid (JA) content had a positive or negative correlation with fruit length or diameter only at certain stages. Neither fruit length nor diameter had significant correlations to abscisic acid (ABA) content. These results indicate that variations in fruit size and shape of different cucumber inbred lines mainly result from the differences in fruit cell number and endogenous IAA content. The present work is the first to propose cucumber fruit size and shape changes from the combined aspects of morphology, cellular levels, and endogenous hormones.

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.

Quantitative Trait Loci for Seed Size Variation in Cucurbits - A Review

Cucurbits (Cucurbitaceae family) include many economically important fruit vegetable crops such as watermelon, pumpkin/squash, cucumber, and melon. Seed size (SS) is an important trait in cucurbits breeding, which is controlled by quantitative trait loci (QTL). Recent advances have deciphered several signaling pathways underlying seed size variation in model plants such as Arabidopsis and rice, but little is known on the genetic basis of SS variation in cucurbits. Here we conducted literature review on seed size QTL identified in watermelon, pumpkin/squash, cucumber and melon, and inferred 14, 9 and 13 consensus SS QTL based on their physical positions in respective draft genomes. Among them, four from watermelon (ClSS2.2, ClSS6.1, ClSS6.2, and ClSS8.2), two from cucumber (CsSS4.1 and CsSS5.1), and one from melon (CmSS11.1) were major-effect, stable QTL for seed size and weight. Whole genome sequence alignment revealed that these major-effect QTL were located in syntenic regions across different genomes suggesting possible structural and functional conservation of some important genes for seed size control in cucurbit crops. Annotation of genes in the four watermelon consensus SS QTL regions identified genes that are known to play important roles in seed size control including members of the zinc finger protein and the E3 ubiquitin-protein ligase families. The present work highlights the utility of comparative analysis in understanding the genetic basis of seed size variation, which may help future mapping and cloning of seed size QTL in cucurbits.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

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.

Genetic architecture of fruit size and shape variation in cucurbits: a comparative perspective

The Cucurbitaceae family hosts many economically important fruit vegetables (cucurbits) such as cucumber, melon, watermelon, pumpkin/squash, and various gourds. The cucurbits are probably best known for the diverse fruit sizes and shapes, but little is known about their genetic basis and molecular regulation. Here, we reviewed the literature on fruit size (FS), shape (FSI), and fruit weight (FW) QTL identified in cucumber, melon, and watermelon, from which 150 consensus QTL for these traits were inferred. Genome-wide survey of the three cucurbit genomes identified 253 homologs of eight classes of fruit or grain size/weight-related genes cloned in Arabidopsis, tomato, and rice that encode proteins containing the characteristic CNR (cell number regulator), CSR (cell size regulator), CYP78A (cytochrome P450), SUN, OVATE, TRM (TONNEAU1 Recruiting Motif), YABBY, and WOX domains. Alignment of the consensus QTL with candidate gene homologs revealed widespread structure and function conservation of fruit size/shape gene homologs in cucurbits, which was exemplified with the fruit size/shape candidate genes CsSUN25-26-27a and CsTRM5 in cucumber, CmOFP1a in melon, and ClSUN25-26-27a in watermelon. In cucurbits, the andromonoecy (for 1-aminocyclopropane-1-carboxylate synthase) and the carpel number (for CLAVATA3) loci are known to have pleiotropic effects on fruit shape, which may complicate identification of fruit size/shape candidate genes in these regions. The present work illustrates the power of comparative analysis in understanding the genetic architecture of fruit size/shape variation, which may facilitate QTL mapping and cloning for fruit size-related traits in cucurbits. The limitations and perspectives of this approach are also discussed.

Analysis of Protein Synthesis in Cucumber Leaves after Inoculation with Corynespora cassiicola: A Proteomic Approach

Cucumber target leaf spot (TLS) disease caused by Corynespora cassiicola has become one of the most important fungal foliar diseases of cultivated cucumbers. However, the defense mechanisms of cucumber plants (Cucumis sativus) against C. cassiicola are still poorly understood. Here, proteins from resistant cucumber plants were analyzed using iTRAQ (isobaric tags for relative and absolute quantification) method. A total of 286 differentially expressed proteins were identified (p < 0.05, ratio > 1.2 or < 0.83) 6 and 24 h after pathogen inoculation in the resistant cultivar Jinyou 38 (the data are available via ProteomeXchange; identifier, PXD012903). Some of the early responses to C. cassiicola infection were revealed, and four factors related to the resistance of cucumber plants to TLS were discovered. First, the proteomic approach revealed modulation of signaling pathways in resistant cucumber plants in response to C. cassiicola infection. Second, the plant immune system recognizes the pathogen and initiates expression of immune response proteins, including those related to plant defense, stress response, signal transduction, cell metabolism, and redox regulation. Third, C. cassiicola activates common stress response pathways; in particular, mildew resistance locus O (MLO) proteins were found to play a crucial role in the TLS prevention. Fourth, rapid activation of the carbohydrate and secondary metabolic pathways, modification and reinforcement of cell walls, and adjustment of the apoplastic environment to the highly stressful conditions were crucial in the cucumber resistance to TLS. Overall, our data contribute to the understanding of interactions between plants and their pathogens and provide new insight into molecular processes involved in the resistance of cucumber plants to disease.

Oligo-painting and GISH reveal meiotic chromosome biases and increased meiotic stability in synthetic allotetraploid Cucumis ×hytivus with dysploid parental karyotypes

BACKGROUND

Meiosis of newly formed allopolyploids frequently encounter perturbations induced by the merging of divergent and hybridizable genomes. However, to date, the meiotic properties of allopolyploids with dysploid parental karyotypes have not been studied in detail. The allotetraploid Cucumis ×hytivus (HHCC, 2n = 38) was obtained from interspecific hybridization between C. sativus (CC, 2n = 14) and C. hystrix (HH, 2n = 24) followed by chromosome doubling. The results of this study thus offer an excellent opportunity to explore the meiotic properties of allopolyploids with dysploid parental karyotypes.

RESULTS

In this report, we describe the meiotic properties of five chromosomes (C5, C7, H1, H9 and H10) and two genomes in interspecific hybrids and C. ×hytivus (the 4th and 14th inbred family) through oligo-painting and genomic in situ hybridization (GISH). We show that 1) only two translocations carrying C5-oligo signals were detected on the chromosomes C2 and C4 of one 14th individual by the karyotyping of eight 4th and 36 14th plants based on C5- and C7-oligo painting, and possible cytological evidence was observed in meiosis of the 4th generation; 2) individual chromosome have biases for homoeologous pairing and univalent formation in F₁ hybrids and allotetraploids; 3) extensive H-chromosome autosyndetic pairings (e.g., H-H, 25.5% PMCs) were observed in interspecific F₁ hybrid, whereas no C-chromosome autosyndetic pairings were observed (e.g. C-C); 4) the meiotic properties of two subgenomes have significant biases in allotetraploids: H-subgenome exhibits higher univalent and chromosome lagging frequencies than C-subgenome; and 5) increased meiotic stability in the S₁₄ generation compared with the S₄ generation, including synchronous meiosis behavior, reduced incidents of univalent and chromosome lagging.

CONCLUSIONS

These results suggest that the meiotic behavior of two subgenomes has dramatic biases in response to interspecific hybridization and allopolyploidization, and the meiotic behavior harmony of subgenomes is a key subject of meiosis evolution in C. ×hytivus. This study helps to elucidate the meiotic properties and evolution of nascent allopolyploids with the dysploid parental karyotypes.