Integrated analysis in bi-parental and natural populations reveals CsCLAVATA3 (CsCLV3) underlying carpel number variations in cucumber

Novel Allelic Gene Variations in CmCLAVATA3 (CmCLV3) Were Identified in a Genetic Population of Melon (Cucumis melo L.)

Carpel number (CN) is an important trait affecting the fruit size and shape of melon, which plays a crucial role in determining the overall appearance and market value. A unique non-synonymous single nucleotide polymorphism (SNP) in CmCLAVATA3 (CmCLV3) is responsible for the variation of CN in C. melo ssp. agrestis (hereafter agrestis), but it has been unclear in C. melo ssp. melo (hereafter melo). In this study, one major locus controlling the polymorphism of 5-CN (multi-CN) and 3-CN (normal-CN) in melo was identified using bulked segregant analysis (BSA-seq). This locus was then fine-mapped to an interval of 1.8 Mb on chromosome 12 using a segregating population containing 1451 progeny. CmCLV3 is still present in the candidate region. A new allele of CmCLV3, which contains five other nucleotide polymorphisms, including a non-synonymous SNP in coding sequence (CDS), except the SNP reported in agrestis, was identified in melo. A cis-trans test confirmed that the candidate gene, CmCLV3, contributes to the variation of CNs in melo. The qRT-PCR results indicate that there is no significant difference in the expression level of CmCLV3 in the apical stem between the multi-CN plants and the normal-CN plants. Overall, this study provides a genetic resource for melon fruit development research and molecular breeding. Additionally, it suggests that melo has undergone similar genetic selection but evolved into an independent allele.

Genome-Wide Identification and Characterization of CLAVATA3/EMBRYO SURROUNDING REGION (CLE) Gene Family in Foxtail Millet (Setaria italica L.)

The CLAVATA3/EMBRYO-SURROUNDING REGION (CLE) genes encode signaling peptides that play important roles in various developmental and physiological processes. However, the systematic identification and characterization of CLE genes in foxtail millet (Setaria italica L.) remain limited. In this study, we identified and characterized 41 SiCLE genes in the foxtail millet genome. These genes were distributed across nine chromosomes and classified into four groups, with five pairs resulting from gene duplication events. SiCLE genes within the same phylogenetic group shared similar gene structure and motif patterns, while 34 genes were found to be single-exon genes. All SiCLE peptides harbored the conserved C-terminal CLE domain, with highly conserved positions in the CLE core sequences shared among foxtail millet, Arabidopsis, rice, and maize. The SiCLE genes contained various cis-elements, including five plant hormone-responsive elements. Notably, 34 SiCLE genes possessed more than three types of phytohormone-responsive elements on their promoters. Comparative analysis revealed higher collinearity between CLE genes in maize and foxtail millet, which may be because they are both C4 plants. Tissue-specific expression patterns were observed, with genes within the same group exhibiting similar and specific expression profiles. SiCLE32 and SiCLE41, classified in Group D, displayed relatively high expression levels in all tissues except panicles. Most SiCLE genes exhibited low expression levels in young panicles, while SiCLE6, SiCLE24, SiCLE25, and SiCLE34 showed higher expression in young panicles, with SiCLE24 down-regulated during later panicle development. Greater numbers of SiCLE genes exhibited higher expression in roots, with SiCLE7, SiCLE22, and SiCLE36 showing the highest levels and SiCLE36 significantly down-regulated after abscisic acid (ABA) treatment. Following treatments with ABA, 6-benzylaminopurine (6-BA), and gibberellic acid 3 (GA3), most SiCLE genes displayed down-regulation followed by subsequent recovery, while jasmonic acid (JA) and indole-3-acetic acid (IAA) treatments led to upregulation at 30 min in leaves. Moreover, identical hormone treatments elicited different expression patterns of the same genes in leaves and stems. This comprehensive study enhances our understanding of the SiCLE gene family and provides a foundation for further investigations into the functions and evolution of SiCLE genes in foxtail millet.

Restructuring plant types for developing tailor-made crops

Plants have adapted to different environmental niches by fine-tuning the developmental factors working together to regulate traits. Variations in the developmental factors result in a wide range of quantitative variations in these traits that helped plants survive better. The major developmental pathways affecting plant architecture are also under the control of such pathways. Most notable are the CLAVATA-WUSCHEL pathway regulating shoot apical meristem fate, GID1-DELLA module influencing plant height and tillering, LAZY1-TAC1 module controlling branch/tiller angle and the TFL1-FT determining the floral fate in plants. Allelic variants of these key regulators selected during domestication shaped the crops the way we know them today. There is immense yield potential in the 'ideal plant architecture' of a crop. With the available genome-editing techniques, possibilities are not restricted to naturally occurring variations. Using a transient reprogramming system, one can screen the effect of several developmental gene expressions in novel ecosystems to identify the best targets. We can use the plant's fine-tuning mechanism for customizing crops to specific environments. The process of crop domestication can be accelerated with a proper understanding of these developmental pathways. It is time to step forward towards the next-generation molecular breeding for restructuring plant types in crops ensuring yield stability.

CsTRM5 regulates fruit shape via mediating cell division direction and cell expansion in cucumber

Fruit shape and size are important appearance and yield traits in cucumber, but the underlying genes and their regulatory mechanisms remain poorly understood. Here we identified a mutant with spherical fruits from an Ethyl Methane Sulfonate (EMS)-mutagenized library, named the qiu mutant. Compared with the cylindrical fruit shape in 32X (wild type), the fruit shape in qiu was round due to reduced fruit length and increased fruit diameter. MutMap analysis narrowed the candidate gene in the 6.47 MB range on Chr2, harboring the FS2.1 locus reported previously. A single-nucleotide polymorphism (SNP) (11359603) causing a truncated protein of CsaV3_2G013800, the homolog of tomato fruit shape gene SlTRM5, may underlie the fruit shape variation in the qiu mutant. Knockout of CsTRM5 by the CRISPR-Cas9 system confirmed that CsaV3_2G013800/CsTRM5 was the causal gene responsible for qiu. Sectioning analysis showed that the spherical fruit in qiu resulted mainly from increased and reduced cell division along the transverse and longitudinal directions, respectively. Meanwhile, the repressed cell expansion contributed to the decreased fruit length in qiu. Transcriptome profiling showed that the expression levels of cell-wall-related genes and abscisic acid (ABA) pathway genes were significantly upregulated in qiu. Hormone measurements indicated that ABA content was greatly increased in the qiu mutant. Exogenous ABA application reduced fruit elongation by inhibiting cell expansion in cucumber. Taken together, these data suggest that CsTRM5 regulates fruit shape by affecting cell division direction and cell expansion, and that ABA participates in the CsTRM5-mediated cell expansion during fruit elongation in cucumber.

Primary mapping of quantitative trait loci regulating multivariate horticultural phenotypes of watermelon (Citrullus lanatus L.)

Watermelon fruits exhibit a remarkable diversity of important horticultural phenotypes. In this study, we initiated a primary quantitative trait loci (QTL) mapping to identify the candidate regions controlling the ovary, fruit, and seed phenotypes. Whole genome sequencing (WGS) was carried out for two differentiated watermelon lines, and 350 Mb (96%) and 354 Mb (97%) of re-sequenced reads covered the reference de novo genome assembly, individually. A total of 45.53% non-synonymous single nucleotide polymorphism (nsSNPs) and 54.47% synonymous SNPs (sSNPs) were spotted, which produced 210 sets of novel SNP-based cleaved amplified polymorphism sequence (CAPS) markers by depicting 46.25% co-dominant polymorphism among parent lines and offspring. A biparental F_2:3 mapping population comprised of 100 families was used for trait phenotyping and CAPS genotyping, respectively. The constructed genetic map spanned a total of 2,398.40 centimorgans (cM) in length and averaged 11.42 cM, with 95.99% genome collinearity. A total of 33 QTLs were identified at different genetic positions across the eight chromosomes of watermelon (Chr-01, Chr-02, Chr-04, Chr-05, Chr-06, Chr-07, Chr-10, and Chr-11); among them, eight QTLs of the ovary, sixteen QTLs of the fruit, and nine QTLs of the seed related phenotypes were classified with 5.32-25.99% phenotypic variance explained (PVE). However, twenty-four QTLs were identified as major-effect and nine QTLs were mapped as minor-effect QTLs across the flanking regions of CAPS markers. Some QTLs were exhibited as tightly localized across the nearby genetic regions and explained the pleiotropic effects of multigenic nature. The flanking QTL markers also depicted significant allele specific contributions and accountable genes were predicted for respective traits. Gene Ontology (GO) functional enrichment was categorized in molecular function (MF), cellular components (CC), and biological process (BP); however, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were classified into three main classes of metabolism, genetic information processing, and brite hierarchies. The principal component analysis (PCA) of multivariate phenotypes widely demonstrated the major variability, consistent with the identified QTL regions. In short, we assumed that our identified QTL regions provide valuable genetic insights regarding the watermelon phenotypes and fine genetic mapping could be used to confirm them.

Identification of Bottle Gourd (Lagenaria siceraria) OVATE Family Genes and Functional Characterization of LsOVATE1

The OVATE gene family is a class of conserved transcription factors that play significant roles in plant growth, development, and abiotic stress, and also affect fruit shape in vegetable crops. Bottle gourd (Lagenaria siceraria), commonly known as calabash or gourd, is an annual climber belonging to the Cucurbitaceae family. Studies on bottle gourd OVATE genes are limited. In this study, we performed genome-wide identification of the OVATE gene family in bottle gourd, and identified a total of 20 OVATE family genes. The identified genes were unevenly distributed across 11 bottle gourd chromosomes. We also analyzed the gene homology, amino acid sequence conservation, and three-dimensional protein structure (via prediction) of the 20 OVATE family genes. We used RNA-seq data to perform expression analysis, which found 20 OVATE family genes to be differentially expressed based on spatial and temporal characteristics, suggesting that they have varying functions in the growth and development of bottle gourd. In situ hybridization and subcellular localization analysis showed that the expression characteristics of the LsOVATE1 gene, located on chromosome 7 homologous to OVATE, is a candidate gene for affecting the fruit shape of bottle gourd. In addition, RT-qPCR data from bottle gourd roots, stems, leaves, and flowers showed different spatial expression of the LsOVATE1 gene. The ectopic expression of LsOVATE1 in tomato generated a phenotype with a distinct fruit shape and development. Transgenic-positive plants that overexpressed LsOVATE1 had cone-shaped fruit, calyx hypertrophy, petal degeneration, and petal retention after flowering. Our results indicate that LsOVATE1 could serve important roles in bottle gourd development and fruit shape determination, and provide a basis for future research into the function of LsOVATE1.

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.

A SNP mutation in the CsCLAVATA1 leads to pleiotropic variation in plant architecture and fruit morphogenesis in cucumber (Cucumis sativus L.)

Plant architectures is predominantly determined by branching pattern, internode elongation, phyllotaxis, shoot determinacy and reproductive organs. Domestication or improvement of this critical agronomic trait played an important role in the breakthrough of crop yield. Here, we identified a mutant with fasciated plant architecture, named fas, from an ethyl methanesulfonate (EMS) induced mutant population in cucumber. The mutant exhibited abnormal phyllotaxy, flattened main stem, increased number of floral organs, and significantly shorter and thicker fruits. However, the molecular mechanism conferring this pleiotropic effect remains unknown. Using a map-based cloning strategy, we isolated the gene CsaV3_3G045960, encoding a leucine-rich repeat receptor-like kinase, a putative direct homolog of the Arabidopsis CLAVATA1 protein referred to as CsCLV1. Endogenous hormone assays showed that IAA and GA₃ levels in fas stems and ovaries were significantly reduced. Conformably, RNA-seq analysis showed that CsCLV1 regulates cucumber stem and ovary development by coordinating hormones and transcription factors. Our results contribute to the understanding of the function of CsCLV1 throughout the growth cycle, provide new evidence that the CLV signaling system is functionally conserved in Cucurbitaceae.

Identification and Characterization Roles of Phytoene Synthase (PSY) Genes in Watermelon Development

Phytoene synthase (PSY) plays an essential role in carotenoid biosynthesis. In this study, three ClPSY genes were identified through the watermelon genome, and their full-length cDNA sequences were cloned. The deduced proteins of the three ClPSY genes were ranged from 355 to 421 amino acid residues. Phylogenetic analysis suggested that the ClPSYs are highly conserved with bottle gourd compared to other cucurbit crops PSY proteins. Variation in ClPSY1 expression in watermelon with different flesh colors was observed; ClPSY1 was most highly expressed in fruit flesh and associated with the flesh color formation. ClPSY1 expression was much lower in the white-fleshed variety than the colored fruits. Gene expression analysis of ClPSY genes in root, stem, leaf, flower, ovary and flesh of watermelon plants showed that the levels of ClPSY2 transcripts found in leaves was higher than other tissues; ClPSY3 was dominantly expressed in roots. Functional complementation assays of the three ClPSY genes suggested that all of them could encode functional enzymes to synthesize the phytoene from Geranylgeranyl Pyrophosphate (GGPP). Some of the homologous genes clustered together in the phylogenetic tree and located in the synteny chromosome region seemed to have similar expression profiles among different cucurbit crops. The findings provide a foundation for watermelon flesh color breeding with regard to carotenoid synthesis and also provide an insight for the further research of watermelon flesh color formation.

Biparental genetic mapping reveals that CmCLAVATA3 (CmCLV3) is responsible for the variation in carpel number in melon (Cucumis melo L.)

Genetic analysis revealed that CmCLV3 is a candidate gene for the variation in melon carpel number. Carpel number (CN) is an important trait in melon. Three-CN melon fruit is oval, while 5-CN melon fruit has a round or flat shape. Herein, a genetic analysis of a population in which the CN locus was segregated indicated that 3-CN is controlled by a major dominant effective gene. Bulked segregant analysis and initial linkage mapping placed the CN locus in a 6.67 Mb region on chromosome 12, and it was narrowed to 882.19 kb with molecular markers and recombinant plants. Fine mapping with a large F₂ population containing 1026 individuals further narrowed the locus to an 83.98 kb region harboring five annotated genes. Gene structure alignment between the parental lines revealed MELO3C035640.2 (annotated as CLAVATA3, CmCLV3) as the best candidate gene for the CN trait. CmCLV3 was more highly expressed in 3- than 5-CN lines and specifically expressed in terminal buds rather than in young leaves, hypocotyls, and roots. The CmCLV3 coding region was cloned from eight 3- or 5-CN melon accessions, and a nonsynonymous SNP site was highly correlated with CN variation. This SNP site was also related to CN variations among 40 melon lines according to their resequencing data, causing a helix alteration in the CmCLV3 protein. Promoter region sequence alignment and activity analysis showed that, unlike in cucumber and tomato, CmCLV3 promoter variation and activity were not the main reasons for CN alteration. Overall, this study provides a genetic resource for melon fruit development research and molecular breeding tools for melon CN improvement.

The Formation of Hollow Trait in Cucumber (Cucumis sativus L.) Fruit Is Controlled by CsALMT2

The hollow trait is crucial for commercial quality of cucumber (Cucumis sativus L.) fruit, and its molecular regulatory mechanism is poorly understood due to its environmental sensitivity. In the previous research, we obtained the hollow and the non-hollow materials of ecotype cucumbers of South China, which were not easily affected by the external environment through a systematic breeding method. In this study, first, we proposed to use the percentage of the hollow area as the criterion to compare the hollow characteristics between two materials, and to analyze the formation mechanism of early hollow trait from the perspective of cytology. The results showed that the hollow trait occurred in the early stage of fruit development, and formed with the opening of carpel ventral zipped bi-cell layer, which formed rapidly from 2 to 4 days, and then slowed to a constant rate from 14 to 16 days. Meanwhile, the different genetic populations were constructed using these materials, and fine mapping was performed by bulked segregant analysis (BSA) and kompetitive allele specific PCR (KASP) method. The Csa1G630860 (CsALMT2), encoding protein ALMT2, was determined as a candidate gene for regulating the hollow trait in fruit. Furthermore, the expression profile of CsALMT2 was analyzed by qRT-PCR and fluorescence in situ hybridization. The expression of CsALMT2 had obvious tissue specificity, and it was abundantly expressed in the ovule development zone inside the fruit. In the hollow material of cucumber fruit, the expression of CsALMT2 was significantly downregulated. The subcellular localization in tobacco leaves indicated that CsALMT2 was distributed on the plasma membrane. In conclusion, in this study, for the first time, we found the regulatory gene of hollow trait in cucumber fruit, which laid the foundation for subsequent research on the molecular mechanism of hollow trait formation in cucumber fruit, and made it possible to apply this gene in cucumber breeding.

A cryptic variation in a member of the Ovate Family Proteins is underlying the melon fruit shape QTL fsqs8.1

The gene underlying the melon fruit shape QTL fsqs8.1 is a member of the Ovate Family Proteins. Variation in fruit morphology is caused by changes in gene expression likely due to a cryptic structural variation in this locus. Melon cultivars have a wide range of fruit morphologies. Quantitative trait loci (QTL) have been identified underlying such diversity. This research focuses on the fruit shape QTL fsqs8.1, previously detected in a cross between the accession PI 124112 (CALC, producing elongated fruit) and the cultivar 'Piel de Sapo' (PS, producing oval fruit). The CALC fsqs8.1 allele induced round fruit shape, being responsible for the transgressive segregation for this trait observed in that population. In fact, the introgression line CALC8-1, carrying the fsqs8.1 locus from CALC into the PS genetic background, produced perfect round fruit. Following a map-based cloning approach, we found that the gene underlying fsqs8.1 is a member of the Ovate Family Proteins (OFP), CmOFP13, likely a homologue of AtOFP1 and SlOFP20 from Arabidopsis thaliana and tomato, respectively. The induction of the round shape was due to the higher expression of the CALC allele at the early ovary development stage. The fsqs8.1 locus showed an important structural variation, being CmOFP13 surrounded by two deletions in the CALC genome. The deletions are present at very low frequency in melon germplasm. Deletions and single nucleotide polymorphisms in the fsqs8.1 locus could not be not associated with variation in fruit shape among different melon accessions, what indicates that other genetic factors should be involved to induce the CALC fsqs8.1 allele effects. Therefore, fsqs8.1 is an example of a cryptic variation that alters gene expression, likely due to structural variation, resulting in phenotypic changes in melon fruit morphology.

Genome-wide identification of CLE gene family and their potential roles in bolting and fruit bearing in cucumber (Cucumis sativus L.)

BACKGROUND

Signal peptides are essential for plant growth and development. In plants, biological processes including cell-cell communication, cellular proliferation and differentiation, cellular determination of self-incompatibility, and defensive responses, all depend heavily on peptide-signaling networks such as CLE (CLAVATA3/Embryo surrounding region-related). The CLEs are indispensable in different periods of plant growth and development, especially in maintaining the balance between proliferation and differentiation of stem cells in various meristematic tissues. The working system of CLE genes in cucumber, an important economical vegetable (Cucumis sativus L.), has not been fully studied yet. The distributional patterns of chromosome-level genome assembly in cucumber provide a fundamental basis for a genome-wide comparative analysis of CLE genes in such plants.

RESULTS

A total of 26 individual CLE genes were identified in Chinese long '9930' cucumber, the majority of which belong to unstable short alkaline and hydrophilic peptides. A comparative analysis showed a close relationship in the development of CLE genes among Arabidopsis thaliana, melon, and cucumber. Half of the exon-intron structures of all CsCLEs genes are single-exon genes, and motif 1, a typical CLE domain near the C-terminal functioning in signal pathways, is found in all cucumber CLE proteins but CsCLE9. The analysis of CREs (Cis-Regulatory Elements) in the upstream region of the 26 cucumber CLE genes indicates a possible relationship between CsCLE genes and certain functions of hormone response elements. Cucumber resulted closely related to Arabidopsis and melon, having seven and 15 orthologous CLE genes in Arabidopsis and melon, respectively. Additionally, the calculative analysis of a pair of orthologous genes in cucumber showed that as a part of the evolutionary process, CLE genes are undergoing a positive selection process which leads to functional differentiation. The specific expression of these genes was vigorous at the growth and development period and tissues. Cucumber gene CLV3 was overexpressed in Arabidopsis, more than half of the transformed plants in T1 generation showed the phenomena of obvious weakness of the development of growing point, no bolting, and a decreased ability of plant growth. Only two bolted strains showed that either the pod did not develop or the pod was short, and its development was significantly inferior to that in the wild type.

CONCLUSIONS

In this study, 26 CLE genes were identified in Chinese long '9930' cucumber genome. The CLE genes were mainly composed of alkaline hydrophilic unstable proteins. The genes of the CLE family were divided into seven classes, and shared close relationships with their homologs in Arabidopsis and melon. The specific expression of these genes was evaluated in different periods of growth and tissue development, and CLV3, which the representative gene of the family, was overexpressed in Arabidopsis, suggesting that it has a role in bolting and fruit bearing in cucumber.

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.

Transcriptome profiling reveals the occurrence mechanism of bisexual flowers in melon (Cucumis melo L.)

Most cultivated melons are andromonoecies in which male flowers arose both in main stem and lateral branches but bisexual flowers only emerged from the leaf axils of lateral branches. However, bisexual flowers emerged in leaf axils of main stem after ethephon treatment. Therefore, the mechanism regulating the occurrence of bisexual flowers were investigated by performing transcriptome analysis in two comparison sets: shoot apex of main stem (MA) versus that of lateral branches (LA), and shoot apex of main stem after ethephon treatment (Eth) versus control (Cont). KEGG results showed that genes involved in "plant hormone signal transduction", "MAPK signaling pathway" and "carbon metabolism" were significantly upregulated both in LA and Eth. Further, details of DEGs involved in ethylene signaling pathway were surveyed and six genes were co-upregulated in two comparison sets. Among these, CmERF1, downstream in ethylene signaling pathway, showed the most significantly difference and expressed higher in bisexual buds than that in male buds. Furthermore, fifteen DEGs were found to contain GCC box or CRT/DRE cis-element for CmERF1 in their putative promoter region, and these DEGs involved in several plant hormones signaling pathway, camalexin synthesis, carbon metabolism and plant pathogen interaction.

Resequencing of 297 melon accessions reveals the genomic history of improvement and loci related to fruit traits in melon

Domestication and improvement are two important stages in crop evolution. Melon (Cucumis melo L.) is an important vegetable crop with wide phenotypic diversity in many horticultural traits, especially fruit size, flesh thickness and aroma, which are likely the results of long-term extensive selection during its evolution. However, selective signals in domestication and improvement stages for these remarkable variations remain unclear. We resequenced 297 wild, landrace and improved melon accessions and obtained 2 045 412 high-quality SNPs. Population structure and genetic diversity analyses revealed independent and two-step selections in two subspecies of melon: ssp. melo and ssp. agrestis during melon breeding. We detected 233 (~18.35 Mbp) and 159 (~17.71 Mbp) novel potential selective signals during the improvement stage in ssp. agrestis and spp. melo, respectively. Two alcohol acyltransferase genes (CmAATs) unique to the melon genome compared with other cucurbit crops may have undergone stronger selection in ssp. agrestis for the characteristic aroma as compared with other cucurbits. Genome-wide association analysis identified eight fruit size and seven flesh thickness signals overlapping with selective sweeps. Compared with thin-skinned ssp. agrestis, thick-skinned ssp. melo has undergone a stronger selection for thicker flesh. In most melon accessions, CmCLV3 has pleiotropic effects on carpel number and fruit shape. Findings from this study provide novel insights into melon crop evolution, and new tools to advance melon breeding.

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.

Genome-based breeding approaches in major vegetable crops

Vegetable crops are major nutrient sources for humanity and have been well-cultivated since thousands of years of domestication. With the rapid development of next-generation sequencing and high-throughput genotyping technologies, the reference genome of more than 20 vegetables have been well-assembled and published. Resequencing approaches on large-scale germplasm resources have clarified the domestication and improvement of vegetable crops by human selection; its application on genetic mapping and quantitative trait locus analysis has led to the discovery of key genes and molecular markers linked to important traits in vegetables. Moreover, genome-based breeding has been utilized in many vegetable crops, including Solanaceae, Cucurbitaceae, Cruciferae, and other families, thereby promoting molecular breeding at a single-nucleotide level. Thus, genome-wide SNP markers have been widely used, and high-throughput genotyping techniques have become one of the most essential methods in vegetable breeding. With the popularization of gene editing technology research on vegetable crops, breeding efficiency can be rapidly increased, especially by combining the genomic and variomic information of vegetable crops. This review outlines the present genome-based breeding approaches used for major vegetable crops to provide insights into next-generation molecular breeding for the increasing global population.

Gene regulatory network controlling carpel number variation in cucumber

The WUSCHEL-CLAVATA3 pathway genes play an essential role in shoot apical meristem maintenance and floral organ development, and under intense selection during crop domestication. The carpel number is an important fruit trait that affects fruit shape, size and internal quality in cucumber, but the molecular mechanism remains elusive. Here, we found that CsCLV3 expression was negatively correlated with carpel number in cucumber cultivars. CsCLV3-RNAi led to increased number of petals and carpels, whereas overexpression of CsWUS resulted in more sepals, petals and carpels, suggesting that CsCLV3 and CsWUS function as a negative and a positive regulator for carpel number variation, respectively. Biochemical analyses indicated that CsWUS directly bound to the promoter of CsCLV3 and activated its expression. Overexpression of CsFUL1^A , a FRUITFULL-like MADS-box gene, resulted in more petals and carpels. CsFUL1^A can directly bind to the CsWUS promoter to stimulate its expression. Furthermore, we found that auxin participated in carpel number variation in cucumber through interaction of CsARF14 with CsWUS. Therefore, we have identified a gene regulatory pathway involving CsCLV3, CsWUS, CsFUL1^A and CsARF14 in determining carpel number variation in an important vegetable crop - cucumber.

Molecularly tagged genes and quantitative trait loci in cucumber with recommendations for QTL nomenclature

Cucumber, Cucumis sativus L. (2n = 2x = 14), is an important vegetable crop worldwide. It was the first specialty crop with a publicly available draft genome. Its relatively small, diploid genome, short life cycle, and self-compatible mating system offers advantages for genetic studies. In recent years, significant progress has been made in molecular mapping, and identification of genes and QTL responsible for key phenotypic traits, but a systematic review of the work is lacking. Here, we conducted an extensive literature review on mutants, genes and QTL that have been molecularly mapped or characterized in cucumber. We documented 81 simply inherited trait genes or major-effect QTL that have been cloned or fine mapped. For each gene, detailed information was compiled including chromosome locations, allelic variants and associated polymorphisms, predicted functions, and diagnostic markers that could be used for marker-assisted selection in cucumber breeding. We also documented 322 QTL for 42 quantitative traits, including 109 for disease resistances against seven pathogens. By alignment of these QTL on the latest version of cucumber draft genomes, consensus QTL across multiple studies were inferred, which provided insights into heritable correlations among different traits. Through collaborative efforts among public and private cucumber researchers, we identified 130 quantitative traits and developed a set of recommendations for QTL nomenclature in cucumber. This is the first attempt to systematically summarize, analyze and inventory cucumber mutants, cloned or mapped genes and QTL, which should be a useful resource for the cucurbit research community.

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.

Quantitative Trait Loci Mapping and Candidate Gene Analysis of Low Temperature Tolerance in Cucumber Seedlings

Cucumber (Cucumis sativus L.) is an economically important vegetable crop worldwide, but it is sensitive to low temperatures. Cucumber seedlings exposed to long-term low temperature stress (LT), i.e., below 20°C during the day, and 8°C at night, exhibit leaf yellowing, accelerated senescence, and reduced yield, therefore posing a threat to cucumber production. Studying the underlying mechanisms involved in LT tolerance in cucumber seedlings, and developing germplasm with improved LT-tolerance could provide fundamental solutions to the problem. In this study, an F₂ population was generated from two parental lines, "CG104" (LT-tolerant inbred line) and "CG37" (LT-sensitive inbred line), to identify loci that are responsible for LT tolerance in cucumber seedlings. Replicated phenotypic analysis of the F₂-derived F₃ family using a low-temperature injury index (LTII) suggested that the LT tolerance of cucumber seedlings is controlled by multiple genes. A genetic map of 990.8 cM was constructed, with an average interval between markers of 5.2 cM. One quantitative trait loci (QTL) named qLTT5.1 on chromosome 5, and two QTLs named qLTT6.1 and qLTT6.2 on chromosome 6 were detected. Among them, qLTT6.2 accounted for 26.8 and 24.1% of the phenotypic variation in two different experiments. Single-nucleotide polymorphism (SNP) variations within the region of qLTT6.2 were analyzed using two contrasting in silico bulks generated from the cucumber core germplasm. Result showed that 214 SNPs were distributed within the 42-kb interval, containing three candidate genes. Real-time quantitative reverse transcription PCR and sequence analysis suggested that two genes Csa6G445210, an auxin response factor, and Csa6G445230, an ethylene-responsive transmembrane protein, might be candidate genes responsible for LT tolerance in cucumber seedlings. This study furthers the understanding of the molecular mechanism underlying LT tolerance in cucumber seedlings, and provides new markers for molecular breeding.

Comparative transcriptome and coexpression network analysis of carpel quantitative variation in Paeonia rockii

BACKGROUND

Quantitative variation of floral organs in plants is caused by an extremely complex process of transcriptional regulation. Despite progress in model plants, the molecular mechanisms of quantitative variation remain unknown in woody flower plants. The Paeonia rockii originated in China is a precious woody plant with ornamental, medicinal and oil properties. There is a wide variation in the number of carpel in P. rockii, but the molecular mechanism of the variation has rarely been studied. Then a comparative transcriptome was performed among two cultivars of P. rockii with different development patterns of carpel in this study.

RESULTS

Through the next-generation and single-molecule long-read sequencing (NGS and SMLRS), 66,563 unigenes and 28,155 differentially expressed genes (DEGs) were identified in P. rockii. Then clustering pattern and weighted gene coexpression network analysis (WGCNA) indicated that 15 candidate genes were likely involved in the carpel quantitative variation, including floral organ development, transcriptional regulatory and enzyme-like factors. Moreover, transcription factors (TFs) from the MYB, WD, RING1 and LRR gene families suggested the important roles in the management of the upstream genes. Among them, PsMYB114-like, PsMYB12 and PsMYB61-like from the MYB gene family were probably the main characters that regulated the carpel quantitative variation. Further, a hypothetical model for the regulation pattern of carpel quantitative variation was proposed in which the candidate genes function synergistically the quantitative variation process.

CONCLUSIONS

We present the high-quality sequencing products in P. rockii. Our results summarize a valuable collective of gene expression profiles characterizing the carpel quantitative variation. The DEGs are candidate for functional analyses of genes regulating the carpel quantitative variation in tree peonies, which provide a precious resource that reveals the molecular mechanism of carpel quantitative variation in other woody flower crops.

CLAVATA signaling pathway genes modulating flowering time and flower number in chickpea

A combinatorial genomic strategy delineated functionally relevant natural allele of a CLAVATA gene and its marker (haplotype)-assisted introgression led to development of the early-flowering chickpea cultivars with high flower number and enhanced yield/productivity. Unraveling the genetic components involved in CLAVATA (CLV) signaling is crucial for modulating important shoot apical meristem (SAM) characteristics and ultimately regulating diverse SAM-regulated agromorphological traits in crop plants. A genome-wide scan identified 142 CLV1-, 28 CLV2- and 6 CLV3-like genes, and their comprehensive genomic constitution and phylogenetic relationships were deciphered in chickpea. The QTL/fine mapping and map-based cloning integrated with high-resolution association analysis identified SNP loci from CaCLV3_01 gene within a major CaqDTF1.1/CaqFN1.1 QTL associated with DTF (days to 50% flowering) and FN (flower number) traits in chickpea, which was further ascertained by quantitative expression profiling. Molecular haplotyping of CaCLV3_01 gene, expressed specifically in SAM, constituted two major haplotypes that differentiated the early-DTF and high-FN chickpea accessions from late-DTF and low-FN. Enhanced accumulation of transcripts of superior CaCLV3_01 gene haplotype and known flowering promoting genes was observed in the corresponding haplotype-introgressed early-DTF and high-FN near-isogenic lines (NILs) with narrow SAM width. The superior haplotype-introgressed NILs exhibited early-flowering, high-FN and enhanced seed yield/productivity without compromising agronomic performance. These delineated molecular signatures can regulate DTF and FN traits through SAM proliferation and differentiation and thereby will be useful for translational genomic study to develop early-flowering cultivars with enhanced yield/productivity.

Molecular basis of cucumber fruit domestication

Cucumber (Cucumis sativus L.) is an economically important vegetable crop that is cultivated worldwide. Compared to the wild ancestor bearing small, bitter and seedy fruit, domesticated cucumbers exhibit significant variation in fruit appearance, size and flavor. Understanding the molecular basis of domestication related traits can provide insights into fruit evolution and make crop breeding more efficient. Here we review recent advances in relating to the genetic basis of fruit morphological traits (femaleness, fruit spine, wart, size, color and carpel development) and organoleptic features (bitterness) during cucumber domestication.

STAYGREEN, STAY HEALTHY: a loss-of-susceptibility mutation in the STAYGREEN gene provides durable, broad-spectrum disease resistances for over 50 years of US cucumber production

The Gy14 cucumber (Cucumis sativus) is resistant to oomyceteous downy mildew (DM), bacterial angular leaf spot (ALS) and fungal anthracnose (AR) pathogens, but the underlying molecular mechanisms are unknown. Quantitative trait locus (QTL) mapping for the disease resistances in Gy14 and further map-based cloning identified a candidate gene for the resistant loci, which was validated and functionally characterized by spatial-temporal gene expression profiling, allelic diversity and phylogenetic analysis, as well as local association studies. We showed that the triple-disease resistances in Gy14 were controlled by the cucumber STAYGREEN (CsSGR) gene. A single nucleotide polymorphism (SNP) in the coding region resulted in a nonsynonymous amino acid substitution in the CsSGR protein, and thus disease resistance. Genes in the chlorophyll degradation pathway showed differential expression between resistant and susceptible lines in response to pathogen inoculation. The causal SNP was significantly associated with disease resistances in natural and breeding populations. The resistance allele has undergone selection in cucumber breeding. The durable, broad-spectrum disease resistance is caused by a loss-of-susceptibility mutation of CsSGR. Probably, this is achieved through the inhibition of reactive oxygen species over-accumulation and phytotoxic catabolite over-buildup in the chlorophyll degradation pathway. The CsSGR-mediated host resistance represents a novel function of this highly conserved gene in plants.

CsLFY is required for shoot meristem maintenance via interaction with WUSCHEL in cucumber (Cucumis sativus)

Cucumber (Cucumis sativus) is an agronomically important vegetable with indeterminant growth habit, in which leaves are produced from the shoot apical meristem (SAM), and unisexual flowers are generated from the leaf axils. LEAFY (LFY) and its homologs have been shown to play important roles in promoting flower development and branching. The LFY homolog gene CsLFY was cloned in cucumber. Molecular biology, developmental biology and biochemical tools were combined to explore the biological function of the LFY homologous gene CsLFY in cucumber. CsLFY was expressed in the SAM, floral meristem and floral organ primordia. Ectopic expression of CsLFY rescued the phenotype of the lfy-5 mutant in Arabidopsis. Knockdown of CsLFY by RNA interference (RNAi) led to defective shoot development and premature discontinuance of leaf initiation in cucumber. Transcription of CsWUS and putative CsLFY target genes including CsAP3 and CUM1 were significantly reduced in the CsLFY-RNAi lines. Further biochemical analyses indicated that CsLFY physically interacts with CsWUS in cucumber. These data suggested that CsLFY has a novel function in regulating shoot meristem maintenance through interaction with CsWUS, and promotes flower development via activation of CsAP3 and CUM1 in cucumber.

QTL mapping of domestication and diversifying selection related traits in round-fruited semi-wild Xishuangbanna cucumber (Cucumis sativus L. var. xishuangbannanesis)

QTL analysis revealed 11 QTL underlying flowering time and fruit size variation in the semi-wild Xishuangbanna cucumber, of which, FT6.2 and FS5.2 played the most important roles in determining photoperiod-dependent flowering time and round-fruit shape, respectively. Flowering time and fruit size are two important traits in domestication and diversifying selection in cucumber, but their genetic basis is not well understood. Here we reported QTL mapping results on flowering time and fruit size with F₂ and F_2:3 segregating populations derived from the cross between WI7200, a small fruited, early flowering primitive cultivated cucumber and WI7167, a round-fruited, later flowering semi-wild Xishuangbanna (XIS) cucumber. A linkage map with 267 microsatellite marker loci was developed with 138 F₂ plants. Phenotypic data of male and female flowering time, fruit length and diameter and three other traits (mature fruit weight and number, and seedling hypocotyl length) were collected in multiple environments. Three flowering time QTL, FT1.1, FT5.1 and FT6.2 were identified, in which FT6.2 played the most important role in conferring less photoperiod sensitive early flowering during domestication whereas FT1.1 seemed more influential in regulating flowering time within the cultivated cucumber. Eight consensus fruit size QTL distributed in 7 chromosomes were detected, each of which contributed to both longitudinal and radial growth in cucumber fruit development. Among them, FS5.2 on chromosome 5 exhibited the largest effect on the determination of round fruit shape that was characteristic of the WI7167 XIS cucumber. Possible roles of these flowering time and fruit size QTL in domestication of cucumber and crop evolution of the semi-wild XIS cucumber, as well as the genetic basis of round fruit shape in cucumber are discussed.

Round fruit shape in WI7239 cucumber is controlled by two interacting quantitative trait loci with one putatively encoding a tomato SUN homolog

QTL analysis revealed two interacting loci, FS1.2 and FS2.1, underlying round fruit shape in WI7239 cucumber; CsSUN , a homolog of tomato fruit shape gene SUN , was a candidate for FS1.2. Fruit size is an important quality and yield trait in cucumber, but its genetic basis remains poorly understood. Here we reported QTL mapping results on fruit size with segregating populations derived from the cross between WI7238 (long fruit) and WI7239 (round fruit) inbred cucumber lines. Phenotypic data of fruit length and diameter were collected at anthesis, immature and mature fruit stages in four environments. Ten major-effect QTL were detected for six traits; synthesis of information from these QTL supported two genes, FS1.2 and FS2.1, underlying fruit size variation in the examined populations. Under the two-gene model, deviation from expected segregation ratio in fruit length and diameter among segregating populations was observed, which could be explained mainly by the interactions between FS1.2 and FS2.1, and segregation distortion in the FS2.1 region. Genome-wide candidate gene search identified CsSUN, a homolog of the tomato fruit shape gene SUN, as the candidate for FS1.2. The round-fruited WI7239 had a 161-bp deletion in the first exon of CsSUN, and its expression in WI7239 was significantly lower than that in WI7238. A marker derived from this deletion was mapped at the peak location of FS1.2 in QTL analysis. Comparative analysis suggested the melon gene CmSUN-14, a homolog of CsSUN as a candidate of the fl2/fd2/fw2 QTL in melon. This study revealed the unique genetic architecture of round fruit shape in WI7239 cucumber. It also highlights the power of QTL analysis for traits with a simple genetic basis but their expression is complicated by other factors.

A valid name for the Xishuangbanna gourd, a cucumber with carotene-rich fruits

Herbarium specimens deposited in publicly accessible collections are the basis for all scientific names because only permanent specimens can be re-studied by independent researchers, the very essence of science. Re-investigations may be done with morphological, chemical, genomic, computer-tomographic, or other methods. Based on new herbarium material, I here provide a name for the Xishuangbanna gourd, a plant long cultivated in Yunnan because of its large non-bitter fruits, rich in β-carotene. Genome re-sequencing of numerous accessions has shown that this cucumber mutant is closer to Cucumis sativus var. sativus than is the wild bitter-fruited progenitor C. sativus var. hardwickii, and two dozen studies have further clarified the genetics of key traits, including pulp color, fruit shape, and flowering times. Morphological and molecular diagnoses of the new variety are provided and museum-quality specimens have been distributed to the World's major herbaria.

The chlorophyll-deficient golden leaf mutation in cucumber is due to a single nucleotide substitution in CsChlI for magnesium chelatase I subunit

The cucumber chlorophyll-deficient golden leaf mutation is due to a single nucleotide substitution in the CsChlI gene for magnesium chelatase I subunit which plays important roles in the chlorophyll biosynthesis pathway. The Mg-chelatase catalyzes the insertion of Mg(2+) into the protoporphyrin IX in the chlorophyll biosynthesis pathway, which is a protein complex encompassing three subunits CHLI, CHLD, and CHLH. Chlorophyll-deficient mutations in genes encoding the three subunits have played important roles in understanding the structure, function and regulation of this important enzyme. In an EMS mutagenesis population, we identified a chlorophyll-deficient mutant C528 with golden leaf color throughout its development which was viable and able to set fruits and seeds. Segregation analysis in multiple populations indicated that this leaf color mutation was recessively inherited and the green color showed complete dominance over golden color. Map-based cloning identified CsChlI as the candidate gene for this mutation which encoded the CHLI subunit of cucumber Mg-chelatase. The 1757-bp CsChlI gene had three exons and a single nucleotide change (G to A) in its third exon resulted in an amino acid substitution (G269R) and the golden leaf color in C528. This mutation occurred in the highly conserved nucleotide-binding domain of the CHLI protein in which chlorophyll-deficient mutations have been frequently identified. The mutant phenotype, CsChlI expression pattern and the mutated residue in the CHLI protein suggested the mutant allele in C528 is unique among mutations identified so far in different species. This golden leaf mutant not only has its potential in cucumber breeding, but also provides a useful tool in understanding the CHLI function and its regulation in the chlorophyll biosynthesis pathway as well as chloroplast development.