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

Genetic and molecular regulation of fruit development in cucumber

Fruit development can be generally classified into a set of biologically sequential stages including fruit initiation, growth, and ripening. Cucumber, a globally important vegetable crop, displays two important features during fruit development: parthenocarpy at fruit initiation and prematurity at harvest for consumption. Therefore, fruit growth plays essential role for cucumber yield and quality formation, and has become the research hot spot in cucumber fruit development. Here, we describe recent advances in molecular mechanisms underlying fruit growth in cucumber, include key players and regulatory networks controlling fruit length variation, fruit neck elongation, and locule development. We also provide insights into future directions for scientific research and breeding strategies in cucumber.

Single nucleotide polymorphisms in SEPALLATA 2 underlie fruit length variation in cucurbits

Complete disruption of critical genes is generally accompanied by severe growth and developmental defects, which dramatically hinder its utilization in crop breeding. Identifying subtle changes, such as single-nucleotide polymorphisms (SNPs), in critical genes that specifically modulate a favorable trait is a prerequisite to fulfill breeding potential. Here, we found 2 SNPs in the E-class floral organ identity gene cucumber (Cucumis sativus) SEPALLATA2 (CsSEP2) that specifically regulate fruit length. Haplotype (HAP) 1 (8G2667A) and HAP2 (8G2667T) exist in natural populations, whereas HAP3 (8A2667T) is induced by ethyl methanesulfonate mutagenesis. Phenotypic characterization of 4 near-isogenic lines and a mutant line showed that HAP2 fruits are significantly longer than those of HAP1, and those of HAP3 are 37.8% longer than HAP2 fruit. The increasing fruit length in HAP1-3 was caused by a decreasing inhibitory effect on CRABS CLAW (CsCRC) transcription (a reported positive regulator of fruit length), resulting in enhanced cell expansion. Moreover, a 7638G/A-SNP in melon (Cucumis melo) CmSEP2 modulates fruit length in a natural melon population via the conserved SEP2-CRC module. Our findings provide a strategy for utilizing essential regulators with pleiotropic effects during crop breeding.

Fine genetic mapping and transcriptomic analysis revealed major gene modulating the clear stripe margin pattern of watermelon peel

The peel stripe margin pattern is one of the most important quality traits of watermelon. In this study, two contrasted watermelon lines [slb line (P1) with a clear peel stripe margin pattern and GWAS-38 line (P2) with a blurred peel stripe margin pattern] were crossed, and biparental F2 mapping populations were developed. Genetic segregation analysis revealed that a single recessive gene is modulating the main-effect genetic locus (Clcsm) of the clear stripe margin pattern of peel. Bulked segregant analysis-based sequencing (BSA-Seq) and fine genetic mapping exposed the delimited Clcsm locus to a 19.686-kb interval on chromosome 6, and the Cla97C06G126680 gene encoding the MYB transcription factor family was identified. The gene mutation analysis showed that two non-synonymous single-nucleotide polymorphism (nsSNP) sites [Chr6:28438793 (A-T) and Chr6:28438845 (A-C)] contribute to the clear peel stripe margin pattern, and quantitative real-time polymerase chain reaction (qRT-PCR) also showed a higher expression trend in the slb line than in the GWAS-38 line. Further, comparative transcriptomic analysis identified major differentially expressed genes (DEGs) in three developmental periods [4, 12, and 20 days after pollination (DAP)] of both parental lines. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses indicated highly enriched DEGs involved in metabolic processes and catalytic activity. A total of 44 transcription factor families and candidate genes belonging to the ARR-B transcription factor family are believed to regulate the clear stripe margin trait of watermelon peel. The gene structure, sequence polymorphism, and expression trends depicted significant differences in the peel stripe margin pattern of both parental lines. The ClMYB36 gene showed a higher expression trend for regulating the clear peel stripe margin of the slb line, and the ClAPRR5 gene depicted a higher expression for modulating the blurred peel stripe margin in the GWAS-38 line. Overall, our fine genetic mapping and transcriptomic analysis revealed candidate genes differentiating the clear and blurred peel stripe patterns of watermelon fruit.

The mutation of CsSUN, an IQD family protein, is responsible for the short and fat fruit (sff) in cucumber (Cucumis sativus L.)

The fruit shape of cucumber is an important agronomic trait, and mining regulatory genes, especially dominant ones, is vital for cucumber breeding. In this study, we identified a short and fat fruit mutant, named sff, from an EMS mutagenized population. Compared to the CCMC (WT), sff (MT) exhibited reduced fruit length and increased dimeter. Segregation analysis revealed that the sff phenotype is controlled by a semi-dominant single gene with dosage effects. Through map-based cloning, the SFF locus was narrowed down to a 52.6 kb interval with two SNPs (G651A and C1072T) in the second and third exons of CsaV3_1G039870, which encodes an IQD family protein, CsSUN. The G651A within the IQ domain of CsSUN was identified as the unique SNP among 114 cucumber accessions, and it was the primary cause of the functional alteration in CsSUN. By generating CsSUN knockout lines in cucumber, we confirmed that CsSUN was responsible for sff mutant phenotype. The CsSUN is localized to the plasma membrane. CsSUN exhibited the highest expression in the fruit with lower expression in sff compared to WT. Histological observations suggest that the sff mutant phenotype is due to increased transverse cell division and inhibited longitudinal cell division. Transcriptome analysis revealed that CsSUN significantly affected the expression of genes related to cell division, expansion, and auxin signal transduction. This study unveils CsSUN's crucial role in shaping cucumber fruit and offers novel insights for cucumber breeding.

Quantitative Trait Loci Mapping and Comparative Transcriptome Analysis of Fruit Weight (FW) in Watermelon (Citrullus lanatus L.)

The watermelon (Citrullus lanatus L.) holds substantial economic value as a globally cultivated horticultural crop. However, the genetic architecture of watermelon fruit weight (FW) remains poorly understood. In this study, we used sh14-11 with small fruit and N14 with big fruit to construct 100 recombinant inbred lines (RILs). Based on whole-genome resequencing (WGR), 218,127 single nucleotide polymorphisms (SNPs) were detected to construct a high-quality genetic map. After quantitative trait loci (QTL) mapping, a candidate interval of 31-38 Mb on chromosome 2 was identified for FW. Simultaneously, the bulked segregant analysis (BSA) in the F2 population corroborated the identification of the same interval, encompassing the homologous gene linked to the known FW-related gene fas. Additionally, RNA-seq was carried out across 11 tissues from sh14-11 and N14, revealing expression profiles that identified 1695 new genes and corrected the annotation of 2941 genes. Subsequent differential expression analysis unveiled 8969 differentially expressed genes (DEGs), with 354 of these genes exhibiting significant differences across four key developmental stages. The integration of QTL mapping and differential expression analysis facilitated the identification of 14 FW-related genes, including annotated TGA and NAC transcription factors implicated in fruit development. This combined approach offers valuable insights into the genetic basis of FW, providing crucial resources for enhancing watermelon cultivation.

Identification of CaPCR1, an OFP gene likely involved in pointed versus concave fruit tip regulation in pepper (Capsicum annuum L.) using recombinant inbred lines

KEY MESSAGE

CaPCR1 (Capana12g002165) was a candidate gene regulating fruit concave/pointed tip shape in pepper. The concave shape of the fruit tip in pepper plants is highly susceptible to drought and low temperature stresses, resulting in the appearance of a pointed tip fruit, which affects its commercial value. However, few studies on the process of fruit tip development and regulatory genes in pepper have been reported. Herein, the developmental process of the ovary before anthesis, especially changes in the shape of the ovary tip, was studied in detail. The results showed that the final fruit tip shape was consistent with the ovary tip shape before anthesis, and a concave tip shape gradually developed. F4 recombinant inbred lines (RILs) were constructed to map the genes regulating fruit tip shape through hybridization of the LRS and SBS pepper inbred lines. CaPCR1 (Capana12g002165), an OFP (OVATE Family Protein) family gene, was located in the candidate region on chr12. Three SNPs were found in the protein coding sequence of CaPCR1 between SBS and LRS, but only one SNP led to amino acid variation. Sequence variations, including base replacements, deletions and insertions, were also detected in the gene promoter region. The relative expression level of the CaPCR1 gene was significantly greater in the concave tip ovary than in the pointed tip ovary. qRT‒PCR analysis revealed that the CaPCR1 gene was expressed mainly in the gynoecium, placenta and green fruit pericarp, which was consistent with its function in ovary and fruit development. Taken together, these results suggested that CaPCR1 is a candidate gene involved in fruit tip shape determination in pepper.

CmERF1 acts as a positive regulator of fruits and leaves growth in melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important horticultural and economic crop. ETHYLENE RESPONSE FACTOR1 (ERF1) plays an important role in regulating plant development, and the resistance to multiple biotic and abiotic stresses. In this study, developmental biology, molecular biology and biochemical assays were performed to explore the biological function of CmERF1 in melon. Abundant transcripts of CmERF1 were found in ovary at green-yellow bud (GYB) and rapid enlargement (ORE) stages. In CmERF1 promoter, the cis-regulatory elements for indoleacetic acid (IAA), methyl jasmonate (MeJA), salicylic acid (SA), abscisic acid (ABA), gibberellic acid (GA), light and low temperature responses were found. CmERF1 could be significantly induced by ethylene, IAA, MeJA, SA, ABA, and respond to continuous light and low temperature stresses in melon. Ectopic expression of CmERF1 increased the length of siliqua and carpopodium, and expanded the size of leaves in Arabidopsis. Knockdown of CmERF1 led to smaller ovary at anthesis, mature fruit and leaves in melon. In CmERF1-RNAi #2 plants, 75 genes were differently expressed compared with control, and the promoter regions of 28 differential expression genes (DEGs) contained the GCC-box (AGCCGCC) or DRE (A/GCCGAC) cis-acting elements of CmERF1. A homolog of cell division cycle protein 48 (CmCDC48) was proved to be the direct target of CmERF1 by the yeast one-hybrid assay and dual-luciferase (LUC) reporter (DLR) system. These results indicated that CmERF1 was able to promote the growth of fruits and leaves, and involved in multiple hormones and environmental signaling pathways in melon.

Heterotrimeric Gα-subunit regulates flower and fruit development in CLAVATA signaling pathway in cucumber

Flowers and fruits are the reproductive organs in plants and play essential roles in natural beauty and the human diet. CLAVATA (CLV) signaling has been well characterized as regulating floral organ development by modulating shoot apical meristem (SAM) size; however, the signaling molecules downstream of the CLV pathway remain largely unknown in crops. Here, we found that functional disruption of CsCLV3 peptide and its receptor CsCLV1 both resulted in flowers with extra organs and stumpy fruits in cucumber. A heterotrimeric G protein α-subunit (CsGPA1) was shown to interact with CsCLV1. Csgpa1 mutant plants derived from gene editing displayed significantly increased floral organ numbers and shorter and wider fruits, a phenotype resembling that of Csclv mutants in cucumber. Moreover, the SAM size was enlarged and the longitudinal cell size of fruit was decreased in Csgpa1 mutants. The expression of the classical stem cell regulator WUSCHEL (WUS) was elevated in the SAM, while the expression of the fruit length stimulator CRABS CLAW (CRC) was reduced in the fruit of Csgpa1 mutants. Therefore, the Gα-subunit CsGPA1 protein interacts with CsCLV1 to inhibit floral organ numbers but promote fruit elongation, via repressing CsWUS expression and activating CsCRC transcription in cucumber. Our findings identified a new player in the CLV signaling pathway during flower and fruit development in dicots, increasing the number of target genes for precise manipulation of fruit shape during crop breeding.

Identification and functional characterization of conserved cis-regulatory elements responsible for early fruit development in cucurbit crops

A number of cis-regulatory elements (CREs) conserved during evolution have been found to be responsible for phenotypic novelty and variation. Cucurbit crops such as cucumber (Cucumis sativus), watermelon (Citrullus lanatus), melon (Cucumis melo), and squash (Cucurbita maxima) develop fruits from an inferior ovary and share some similar biological processes during fruit development. Whether conserved regulatory sequences play critical roles in fruit development of cucurbit crops remains to be explored. In six well-studied cucurbit species, we identified 392,438 conserved noncoding sequences (CNSs), including 82,756 that are specific to cucurbits, by comparative genomics. Genome-wide profiling of accessible chromatin regions (ACRs) and gene expression patterns mapped 20,865 to 43,204 ACRs and their potential target genes for two fruit tissues at two key developmental stages in six cucurbits. Integrated analysis of CNSs and ACRs revealed 4,431 syntenic orthologous CNSs, including 1,687 cucurbit-specific CNSs that overlap with ACRs that are present in all six cucurbit crops and that may regulate the expression of 757 adjacent orthologous genes. CRISPR mutations targeting two CNSs present in the 1,687 cucurbit-specific sequences resulted in substantially altered fruit shape and gene expression patterns of adjacent NAC1 (NAM, ATAF1/2, and CUC2) and EXT-like (EXTENSIN-like) genes, validating the regulatory roles of these CNSs in fruit development. These results not only provide a number of target CREs for cucurbit crop improvement, but also provide insight into the roles of CREs in plant biology and during evolution.

A wild melon reference genome provides novel insights into the domestication of a key gene responsible for melon fruit acidity

KEY MESSAGE

A wild melon reference genome elucidates the genomic basis of fruit acidity domestication. Structural variants (SVs) have been reported to impose major effects on agronomic traits, representing a significant contributor to crop domestication. However, the landscape of SVs between wild and cultivated melons is elusive and how SVs have contributed to melon domestication remains largely unexplored. Here, we report a 379-Mb chromosome-scale genome of a wild progenitor melon accession "P84", with a contig N50 of 14.9 Mb. Genome comparison identifies 10,589 SVs between P84 and four cultivated melons with 6937 not characterized in previously analysis of 25 melon genome sequences. Furthermore, the population-scale genotyping of these SVs was determined in 1175 accessions, and 18 GWAS signals including fruit acidity, fruit length, fruit weight, fruit color and sex determination were detected. Based on these genotyped SVs, we identified 3317 highly diverged SVs between wild and cultivated melons, which could be the potential SVs associated with domestication-related traits. Furthermore, we identify novel SVs affecting fruit acidity and proposed the diverged evolutionary trajectories of CmPH, a key regulator of melon fruit acidity, during domestication and selection of different populations. These results will offer valuable resources for genomic studies and genetic improvement in melon.

Transcriptomics reveal useful resources for examining fruit development and variation in fruit size in Coccinia grandis

Introduction

The Cucurbitaceae family comprises many agronomically important members, that bear nutritious fruits and vegetables of great economic importance. Coccinia grandis, commonly known as Ivy gourd, belongs to this family and is widely consumed as a vegetable. Members of this family are known to display an impressive range of variation in fruit morphology. Although there have been studies on flower development in Ivy gourd, fruit development remains unexplored in this crop.

Methods

In this study, comparative transcriptomics of two Ivy gourd cultivars namely "Arka Neelachal Kunkhi" (larger fruit size) and "Arka Neelachal Sabuja" (smaller fruit size) differing in their average fruit size was performed. A de novo transcriptome assembly for Ivy gourd was developed by collecting fruits at different stages of development (5, 10, 15, and 20 days after anthesis i.e. DAA) from these two varieties. The transcriptome was analyzed to identify differentially expressed genes, transcription factors, and molecular markers.

Results

The transcriptome of Ivy gourd consisted of 155205 unigenes having an average contig size of 1472bp. Unigenes were annotated on publicly available databases to categorize them into different biological functions. Out of these, 7635 unigenes were classified into 38 transcription factor (TF) families, of which Trihelix TFs were most abundant. A total of 11,165 unigenes were found to be differentially expressed in both the varieties and the in silico expression results were validated through real-time PCR. Also, 98768 simple sequence repeats (SSRs) were identified in the transcriptome of Ivy gourd.

Discussion

This study has identified a number of genes, including transcription factors, that could play a crucial role in the determination of fruit shape and size in Ivy gourd. The presence of polymorphic SSRs indicated a possibility for marker-assisted selection for crop breeding in Ivy gourd. The information obtained can help select candidate genes that may be implicated in regulating fruit development and size in other fruit crops.

Nutritional Variation on Sequentially Harvested Shoots and Fruits of Lagenaria siceraria Landraces

Lagenaria siceraria (Molina) Standley, a member of the Cucurbitaceae family, is valued for its medicinal and nutritive properties. The nutrient status of sequentially harvested shoots and fruits at various growth stages in L. siceraria has not been documented to date. This study aimed to compare the nutritional status of L. siceraria landrace shoots and fruits harvested at different maturity stages. Micronutrients, macronutrients, and proximate composition of shoots and fruits were determined using inductively coupled plasma-optical emission spectrometry (ICP-OES). Data were subjected to multivariate statistical analysis. The nutrient attributes differed significantly (p < 0.05) within and among landraces at different growth stages. Correlation of nutritional traits was primarily based on shared absorption sites and comparable chemical composition. The first five principal components in shoots and fruits had 90.218 and 89.918% total variability, respectively. The micronutrients Ca, Mg, K, P, and N in shoots and the macronutrients Fe, Zn, Cu, and Al in fruits were the main contributors to variability. The biplot and dendrogram clustered landraces with comparable nutrient values. Shoot traits classified landraces into three major clusters, where clusters I and II grouped landraces with superior and inferior Mg, P, K, N, Fe, K/Ca+Mg, ADL, and protein levels at 42-63 DAS. Cluster III consisted of landraces with more Ca, Na, Mn, Zn, and Cu. However, the fruit nutrient status classified landraces into two major clusters. Cluster I comprised landrace KSC (as a singlet) with the highest Ca, P, N, Mn, Fe, Zn, and Cu contents at various stages of growth (7-21 DAA). Cluster II had landraces with higher K, C/N ratio, Na, moisture content, ash, protein, and fat. The nutritional status of shoots and fruits determined at different stages of growth is essential for selecting the best harvest time and landrace(s) for required daily nutrient intake.

Genetic mapping reveals a candidate gene CmoFL1 controlling fruit length in pumpkin

Fruit length (FL) is an important economical trait that affects fruit yield and appearance. Pumpkin (Cucurbita moschata Duch) contains a wealth genetic variation in fruit length. However, the natural variation underlying differences in pumpkin fruit length remains unclear. In this study, we constructed a F2 segregate population using KG1 producing long fruit and MBF producing short fruit as parents to identify the candidate gene for fruit length. By bulked segregant analysis (BSA-seq) and Kompetitive Allele-Specific PCR (KASP) approach of fine mapping, we obtained a 50.77 kb candidate region on chromosome 14 associated with the fruit length. Then, based on sequence variation, gene expression and promoter activity analyses, we identified a candidate gene (CmoFL1) encoding E3 ubiquitin ligase in this region may account for the variation of fruit length. One SNP variation in promoter of CmoFL1 changed the GT1CONSENSUS, and DUAL-LUC assay revealed that this variation significantly affected the promoter activity of CmoFL1. RNA-seq analysis indicated that CmoFL1 might associated with the cell division process and negatively regulate fruit length. Collectively, our work identifies an important allelic affecting fruit length, and provides a target gene manipulating fruit length in future pumpkin breeding.

Different control of resistance to two Colletotrichum orbiculare pathogenic races 0 and 1 in cucumber (Cucumis sativus L.)

KEY MESSAGE

Quantitative trait locus analysis identified independent novel loci in cucumbers responsible for resistance to races 0 and 1 of the anthracnose fungal pathogen Colletotrichum orbiculare. Cucumbers have been reported to be vulnerable to Colletotrichum orbiculare, causing anthracnose disease with significant yield loss under favorable conditions. The deployment of a single recessive Cssgr gene in cucumber breeding for anthracnose resistance was effective until a recent report on high-virulent strains infecting cucumbers in Japan conquering the resistance. QTL mapping was conducted to identify the resistance loci in the cucumber accession Ban Kyuri (G100) against C. orbiculare strains 104-T and CcM-1 of pathogenic races 0 and 1, respectively. A single dominant locus An5 was detected in the disease resistance hotspot on chromosome 5 for resistance to 104-T. Resistance to CcM-1 was governed by three loci with additive effects located on chromosomes 2 (An2) and 1 (An1.1 and An1.2). Molecular markers were developed based on variant calling between the corresponding QTL regions in the de novo assembly of the G100 genome and the publicly available cucumber genomes. Multiple backcrossed populations were deployed to fine-map An5 locus and narrow the region to approximately 222 kbp. Accumulation of An2 and An1.1 alleles displayed an adequate resistance to CcM-1 strain. This study provides functional molecular markers for pyramiding resistance loci that confer sufficient resistance against anthracnose in cucumbers.

The genetic basis of apple shape and size unraveled by digital phenotyping

Great diversity of shape, size, and skin color is observed among the fruits of different apple genotypes. These traits are critical for consumers and therefore interesting targets for breeding new apple varieties. However, they are difficult to phenotype and their genetic basis, especially for fruit shape and ground color, is largely unknown. We used the FruitPhenoBox to digitally phenotype 525 genotypes of the apple reference population (apple REFPOP) genotyped for 303,148 single nucleotide polymorphism (SNP) markers. From the apple images, 573 highly heritable features describing fruit shape and size as well as 17 highly heritable features for fruit skin color were extracted to explore genotype-phenotype relationships. Out of these features, seven principal components (PCs) and 16 features with the Pearson's correlation r < 0.75 (selected features) were chosen to carry out genome-wide association studies (GWAS) for fruit shape and size. Four PCs and eight selected features were used in GWAS for fruit skin color. In total, 69 SNPs scattered over all 17 apple chromosomes were significantly associated with round, conical, cylindrical, or symmetric fruit shapes and fruit size. Novel associations with major effect on round or conical fruit shapes and fruit size were identified on chromosomes 1 and 2. Additionally, 16 SNPs associated with PCs and selected features related to red overcolor as well as green and yellow ground color were found on eight chromosomes. The identified associations can be used to advance marker-assisted selection in apple fruit breeding to systematically select for desired fruit appearance.

Cytological, Phytohormone, and Transcriptome Analyses Provide Insights into Persimmon Fruit Shape Formation (Diospyros kaki Thunb.)

Fruit shape is an important external feature when consumers choose their preferred fruit varieties. Studying persimmon (Diospyros kaki Thunb.) fruit shape is beneficial to increasing its commodity value. However, research on persimmon fruit shape is still in the initial stage. In this study, the mechanism of fruit shape formation was studied by cytological observations, phytohormone assays, and transcriptome analysis using the long fruit and flat fruit produced by 'Yaoxianwuhua' hermaphroditic flowers. The results showed that stage 2-3 (June 11-June 25) was the critical period for persimmon fruit shape formation. Persimmon fruit shape is determined by cell number in the transverse direction and cell length in the longitudinal direction. High IAA, GA4, ZT, and BR levels may promote long fruit formation by promoting cell elongation in the longitudinal direction, and high GA3 and ABA levels may be more conducive to flat fruit formation by increasing the cell number in the transverse direction and inhibiting cell elongation in the longitudinal direction, respectively. Thirty-two DEGs related to phytohormone biosynthesis and signaling pathways and nine DEGs related to cell division and cell expansion may be involved in the persimmon fruit shape formation process. These results provide valuable information for regulatory mechanism research on persimmon fruit formation.

Fine mapping of TFL, a major gene regulating fruit length in snake gourd (Trichosanthes anguina L)

Fruit length is a crucial agronomic trait of snake gourd (Trichosanthes anguina L); however, genes associated with fruit length have not been characterised. In this study, F2 snake gourd populations were generated by crossing the inbred lines, S1 and S2 (fruit lengths: 110 and 20 cm, respectively). Subsequently, bulk segregant analysis, sequencing, and fine-mapping were performed on the F2 population to identify target genes. Our findings suggest that the fruit length of snake gourd is regulated by a major-effect regulatory gene. Mining of genes regulating fruit length in snake gourd to provide a basis for subsequent selection and breeding of new varieties. Genotype-phenotype association analysis was performed on the segregating F2 population comprising 6,000 plants; the results indicate that the target gene is located on Chr4 (61,846,126-61,865,087 bp, 18.9-kb interval), which only carries the annotated candidate gene, Tan0010544 (designated TFL). TFL belongs to the MADS-box family, one of the largest transcription factor families. Sequence analysis revealed a non-synonymous mutation of base C to G at position 202 in the coding sequence of TFL, resulting in the substitution of amino acid Gln to Glu at position 68 in the protein sequence. Subsequently, an InDel marker was developed to aid the marker-assisted selection of TFL. The TFL in the expression parents within the same period was analysed using quantitative real-time PCR; the TFL expression was significantly higher in short fruits than long fruits. Therefore, TFL can be a candidate gene for determining the fruit length in snake gourd. Collectively, these findings improve our understanding of the genetic components associated with fruit length in snake gourds, which could aid the development of enhanced breeding strategies for plant species.

A nonsynonymous mutation in BhLS, encoding an acyl-CoA N-acyltransferase leads to fruit and seed size variation in wax gourd (Benincasa hispida)

Wax gourd (Benincasa hispida (Thunb.) Cogn., 2n = 2x = 24) is an economically important vegetable crop cultivated widely in many tropical and subtropical regions, including China, India, and Japan. Both fruit and seeds are prized agronomic attributes in wax gourd breeding and production. However, the genetic mechanisms underlying these traits remain largely unexplored. In this study, we observed a strong correlation between fruit size and seed size variation in our mapping population, indicating genetic control by a single gene, BhLS, with large size being dominant over small. Through bulk segregant analysis sequencing and fine mapping with a large F2 population, we precisely located the BhLS gene within a 47.098-kb physical interval on Chromosome 10. Within this interval, only one gene, Bhi10M000649, was identified, showing homology to Arabidopsis HOOKLESS1. A nonsynonymous mutation (G to C) in the second exon of Bhi10M000649 was found to be significantly associated with both fruit and seed size variation in wax gourd. These findings collectively highlight the pleiotropic effect of the BhLS gene in regulating fruit and seed size in wax gourd. Our results offer molecular insights into the variation of fruit and seed size in wax gourd and establish a fundamental framework for breeding wax gourd cultivars with desired traits.

Multi-Omics Analysis Reveals the Distinct Features of Metabolism Pathways Supporting the Fruit Size and Color Variation of Giant Pumpkin

Pumpkin (Cucurbita maxima) is an important vegetable crop of the Cucurbitaceae plant family. The fruits of pumpkin are often used as directly edible food or raw material for a number of processed foods. In nature, mature pumpkin fruits differ in size, shape, and color. The Atlantic Giant (AG) cultivar has the world's largest fruits and is described as the giant pumpkin. AG is well-known for its large and bright-colored fruits with high ornamental and economic value. At present, there are insufficient studies that have focused on the formation factors of the AG cultivar. To address these knowledge gaps, we performed comparative transcriptome, proteome, and metabolome analysis of fruits from the AG cultivar and a pumpkin with relatively small fruit (Hubbard). The results indicate that up-regulation of gene-encoded expansins contributed to fruit cell expansion, and the increased presence of photoassimilates (stachyose and D-glucose) and jasmonic acid (JA) accumulation worked together in terms of the formation of large fruit in the AG cultivar. Notably, perhaps due to the rapid transport of photoassimilates, abundant stachyose that was not converted into glucose in time was detected in giant pumpkin fruits, implying that a unique mode of assimilate unloading is in existence in the AG cultivar. The potential molecular regulatory network of photoassimilate metabolism closely related to pumpkin fruit expansion was also investigated, finding that three MYB transcription factors, namely CmaCh02G015900, CmaCh01G018100, and CmaCh06G011110, may be involved in metabolic regulation. In addition, neoxanthin (a type of carotenoid) exhibited decreased accumulation that was attributed to the down-regulation of carotenoid biosynthesis genes in AG fruits, which may lead to pigmentation differences between the two pumpkin cultivars. Our current work will provide new insights into the potential formation factors of giant pumpkins for further systematic elucidation.

Deciphering the Genomic Characterization of the GGP Gene Family and Expression Verification of CmGGP1 Modulating Ascorbic Acid Biosynthesis in Melon Plants

Ascorbic acid (AsA), also known as vitamin C, is a well-known antioxidant found in living entities that plays an essential role in growth and development, as well as in defensive mechanisms. GDP-L-galactose phosphorylase (GGP) is a candidate gene regulating AsA biosynthesis at the translational and transcriptional levels in plants. In the current study, we conducted genome-wide bioinformatic analysis and pinpointed a single AsA synthesis rate-limiting enzyme gene in melon (CmGGP1). The protein prediction analysis depicted that the CmGGP1 protein does not have a signaling peptide or transmembrane structure and mainly functions in the chloroplast or nucleus. The constructed phylogenetic tree analysis in multispecies showed that the CmGGP1 protein has a highly conserved motif in cucurbit crops. The structural variation analysis of the CmGGP1 gene in different domesticated melon germplasms showed a single non-synonymous type-base mutation and indicated that this gene was selected by domestication during evolution. Wild-type (WT) and landrace (LDR) germplasms of melon depicted close relationships to each other, and improved-type (IMP) varieties showed modern domestication selection. The endogenous quantification of AsA content in both the young and old leaves of nine melon varieties exhibited the major differentiations for AsA synthesis and metabolism. The real-time quantitative polymerase chain reaction (qRT-PCR) analysis of gene co-expression showed that AsA biosynthesis in leaves was greater than AsA metabolic consumption, and four putative interactive genes (MELO3C025552.2, MELO3C007440.2, MELO3C023324.2, and MELO3C018576.2) associated with the CmGGP1 gene were revealed. Meanwhile, the CmGGP1 gene expression pattern was noticed to be up-regulated to varying degrees in different acclimated melons. We believe that the obtained results would provide useful insights for an in-depth genetic understanding of the AsA biosynthesis mechanism, aimed at the development of improving crop plants for melon.

Genome-Wide Identification of Vascular Plant One-Zinc-Finger Gene Family in Six Cucurbitaceae Species and the Role of CmoVOZ2 in Salt and Drought Stress Tolerance

As a plant-specific transcription factor, the vascular plant one-zinc-finger (VOZ) plays a crucial role in regulating various biological processes. In this study, a total of 17 VOZ genes in the Cucurbitaceae family were investigated using various bioinformatics methods. The 17 VOZ genes in Cucurbitaceae are distributed across 16 chromosomes. Based on the affinity of VOZ proteins to AtVOZ proteins, these 17 proteins were categorized into two groups: group I encompassed eight VOZ members, while group II comprised nine VOZ members. The expression profiles of CmoVOZs under various hormonal and abiotic stresses indicated that these genes were induced differentially by JA, ABA, GA, salt, and drought stress. Subsequently, CmoVOZ1 and CmoVOZ2 were found to be transcriptionally active, with the CmoVOZ2 protein being located mainly in the nucleus. Further experiments revealed that yeast cells expressing CmoVOZ2 gene showed increased tolerance to salt stress and drought stress. These results suggest that the VOZ gene family is not only important for plant growth and development but also that this mechanism may be universal across yeast and plants.

Microtubule-associated protein SlMAP70 interacts with IQ67-domain protein SlIQD21a to regulate fruit shape in tomato

Tomato (Solanum lycopersicum) fruit shape is related to microtubule organization and the activity of microtubule-associated proteins (MAPs). However, insights into the mechanism of fruit shape formation from a cell biology perspective remain limited. Analysis of the tissue expression profiles of different microtubule regulators revealed that functionally distinct classes of MAPs, including members of the plant-specific MICROTUBULE-ASSOCIATED PROTEIN 70 (MAP70) and IQ67 DOMAIN (IQD, also named SUN in tomato) families, are differentially expressed during fruit development. SlMAP70-1-3 and SlIQD21a are highly expressed during fruit initiation, which relates to the dramatic microtubule pattern rearrangements throughout this developmental stage of tomato fruits. Transgenic tomato lines overexpressing SlMAP70-1 or SlIQD21a produced elongated fruits with reduced cell circularity and microtubule anisotropy, while their loss-of-function mutants showed the opposite phenotype, harboring flatter fruits. Fruits were further elongated in plants coexpressing both SlMAP70-1 and SlIQD21a. We demonstrated that SlMAP70s and SlIQD21a physically interact and that the elongated fruit phenotype is likely due to microtubule stabilization induced by the SlMAP70-SlIQD21a interaction. Together, our results identify SlMAP70 proteins and SlIQD21a as important regulators of fruit elongation and demonstrate that manipulating microtubule function during early fruit development provides an effective approach to alter fruit shape.

Comparative physiochemical and transcriptomic analysis reveals the influences of cross-pollination on ovary and fruit development in pummelo (Citrus maxima)

'Shuijingmiyou' pummelo (SJ), one of the most popular fruits in Yunnan province of China, is of relatively low fruit shape (FS) quality. In this study, we compared the FS promoting effects of cross pollinations using pollens from seven pummelo varieties, and found that 'Guanximiyou' pummelo (GX) cross-pollination showed the best FS promoting effects on SJ fruits by shortening its fruit neck. To explore the underlying mechanism, physiochemical and transcriptomic differences between self- and cross-pollinated SJ ovaries (SJO and GXO) were investigated. Higher salicylic acid, gibberellin and indole acetic acid contents and superoxide dismutase, peroxidase and catalase activities, and lower polyphenol oxidase activity were determined in GXO compared with SJO. Enrichment analysis of the identified 578 differentially expressed genes (123 up-regulated and 455 down-regulated) in GXO showed that genes involved in solute transport, RNA biosynthesis, phytohormone action and cell wall organization were significantly enriched. The results obtained in this study will be helpful in understanding the influences of cross-pollination on pummelo ovary and fruit development, and can provide the basis for clarifying the underlying mechanism of cross-pollination improved fruit quality.

Molecular Mapping of Putative Genomic Regions Controlling Fruit and Seed Morphology of Watermelon

The genetic regulatory basis of qualitative and quantitative phenotypes of watermelon is being investigated in different types of molecular and genetic breeding studies around the world. In this study, biparental F2 mapping populations were developed over two experimental years, and the collected datasets of fruit and seed traits exhibited highly significant correlations. Whole-genome resequencing of comparative parental lines was performed and detected single nucleotide polymorphism (SNP) loci were converted into cleaved amplified polymorphic sequence (CAPS) markers. The screened polymorphic markers were genotyped in segregating populations and two genetic linkage maps were constructed, which covered a total of 2834.28 and 2721.45 centimorgan (cM) genetic lengths, respectively. A total of 22 quantitative trait loci (QTLs) for seven phenotypic traits were mapped; among them, five stable and major-effect QTLs (PC-8-1, SL-9-1, SWi-9-1, SSi-9-1, and SW-6-1) and four minor-effect QTLs (PC-2-1 and PC-2-2; PT-2-1 and PT-2-2; SL-6-1 and SSi-6-2; and SWi-6-1 and SWi-6-2) were observed with 3.77-38.98% PVE. The adjacent QTL markers showed a good fit marker-trait association, and a significant allele-specific contribution was also noticed for genetic inheritance of traits. Further, a total of four candidate genes (Cla97C09G179150, Cla97C09G179350, Cla97C09G180040, and Cla97C09G180100) were spotted in the stable colocalized QTLs of seed size linked traits (SL-9-1 and SWi-9-1) that showed non-synonymous type mutations. The gene expression trends indicated that the seed morphology had been formed in the early developmental stage and showed the genetic regulation of seed shape formation. Hence, we think that our identified QTLs and genes would provide powerful genetic insights for marker-assisted breeding aimed at improving the quality traits of watermelon.

Molecular Marker-Assisted Mapping, Candidate Gene Identification, and Breeding in Melon (Cucumis melo L.): A Review

Melon (Cucumis melo L.) is an important crop that is cultivated worldwide for its fleshy fruit. Understanding the genetic basis of a plant's qualitative and quantitative traits is essential for developing consumer-favored varieties. This review presents genetic and molecular advances related to qualitative and quantitative phenotypic traits and biochemical compounds in melons. This information guides trait incorporation and the production of novel varieties with desirable horticultural and economic characteristics and yield performance. This review summarizes the quantitative trait loci, candidate genes, and development of molecular markers related to plant architecture, branching patterns, floral attributes (sex expression and male sterility), fruit attributes (shape, rind and flesh color, yield, biochemical compounds, sugar content, and netting), and seed attributes (seed coat color and size). The findings discussed in this review will enhance demand-driven breeding to produce cultivars that benefit consumers and melon breeders.

Telomere-to-telomere genome assembly of melon (Cucumis melo L. var. inodorus) provides a high-quality reference for meta-QTL analysis of important traits

Melon is an important horticultural crop with extensive diversity in many horticultural groups. To explore its genomic diversity, it is necessary to assemble more high-quality complete genomes from different melon accessions. Meanwhile, a large number of QTLs have been mapped in several studies. Integration of the published QTLs onto a complete genome can provide more accurate information for candidate gene cloning. To address these problems, a telomere-to-telomere (T2T) genome of the elite melon landrace Kuizilikjiz (Cucumis melo L. var. inodorus) was de novo assembled and all the published QTLs were projected onto it in this study. The results showed that a high-quality Kuizilikjiz genome with the size of 379.2 Mb and N50 of 31.7 Mb was de novo assembled using the combination of short reads, PacBio high-fidelity long reads, Hi-C data, and a high-density genetic map. Each chromosome contained the centromere and telomeres at both ends. A large number of structural variations were observed between Kuizilikjiz and the other published genomes. A total of 1294 QTLs published in 67 studies were collected and projected onto the T2T genome. Several clustered, co-localized, and overlapped QTLs were determined. Furthermore, 20 stable meta-QTLs were identified, which significantly reduced the mapping intervals of the initial QTLs and greatly facilitated identification of the candidate genes. Collectively, the T2T genome assembly together with the numerous projected QTLs will not only broaden the high-quality genome resources but also provide valuable and abundant QTL information for cloning the genes controlling important traits in melon.

QTL mapping and transcriptomic analysis of fruit length in cucumber

A total of 151 recombinant inbred lines (RILs) were derived from the cross between 'Cucumis sativus L. hardwickii' (HW) and a cultivated Northern Chinese inbred line 'XinTaiMiCi' (XTMC). We used resequencing to construct the genetic map and analyze the genetic background of RIL population, and combined with the phenotypes of RIL population and the analysis of RNA-seq data, we located the major loci controlling the fruit length of cucumber and related analysis. A genetic map containing 600 bin markers was constructed via re-sequencing. Based on the phenotype data collected in two different seasons (spring 2021 and autumn 2022), the major quantitative trait loci (QTLs) controlling cucumber fruit length were located and their transcriptomic analysis carried out. The results revealed three QTLs (Fl2.1, Fl4.1, and Fl6.1) detected repeatedly in the two seasons, of which Fl4.1 was the dominant QTL. From the functional annotation of corresponding genes there, we discovered the gene Csa4G337340 encoding an auxin efflux carrier family protein. The expression of that gene was significantly lower in XTMC and the long-fruit RIL lines than in HW and the short-fruit RIL lines; hence, we speculated the gene could be negatively correlated with the fruit length of cucumber. Transcriptomic analysis showed that 259 differentially expressed genes (DEGs) were enriched in the plant hormone signal transduction pathway. In addition, among those DEGs, 509 transcription factors were detected, these distributed in several transcription factor gene families, such as bHLH, AP2/ErF -ERF, C2H2, and NAC. Therefore, we concluded that the major gene controlling the fruit length of cucumber is located in the interval of Fl4.1, whose gene Csa4G337340 may be involved in the negative regulation of fruit length. Further, genes related to plant hormone signal transduction and several transcription factors were also found involved in the regulation of cucumber fruit length. Our results provide a reference for the fine mapping of major genes and analyzing the mechanism of cucumber fruit length.

Phytochrome-interacting factor PIF3 integrates phytochrome B and UV-B signaling pathways to regulate gibberellin- and auxin-dependent growth in cucumber hypocotyls

In Arabidopsis, the photoreceptors phytochrome B (PhyB) and UV-B resistance 8 (UVR8) mediate light responses that play a major role in regulating photomorphogenic hypocotyl growth, but how they crosstalk to coordinate this process is not well understood. Here we report map-based cloning and functional characterization of an ultraviolet (UV)-B-insensitive, long-hypocotyl mutant, lh1, and a wild-type-like mutant, lh2, in cucumber (Cucumis sativus), which show defective CsPhyB and GA oxidase2 (CsGA20ox-2), a key gibberellic acid (GA) biosynthesis enzyme, respectively. The lh2 mutation was epistatic to lh1 and partly suppressed the long-hypocotyl phenotype in the lh1lh2 double mutant. We identified phytochrome interacting factor (PIF) CsPIF3 as playing a critical role in integrating the red/far-red and UV-B light responses for hypocotyl growth. We show that two modules, CsPhyB-CsPIF3-CsGA20ox-2-DELLA and CsPIF3-auxin response factor 18 (CsARF18), mediate CsPhyB-regulated hypocotyl elongation through GA and auxin pathways, respectively, in which CsPIF3 binds to the G/E-box motifs in the promoters of CsGA20ox-2 and CsARF18 to regulate their expression. We also identified a new physical interaction between CsPIF3 and CsUVR8 mediating CsPhyB-dependent, UV-B-induced hypocotyl growth inhibition. Our work suggests that hypocotyl growth in cucumber involves a complex interplay of multiple photoreceptor- and phytohormone-mediated signaling pathways that show both conservation with and divergence from those in Arabidopsis.

Molecular and genetic regulations of fleshy fruit shape and lessons from Arabidopsis and rice

Fleshy fruit shape is an important external quality trait influencing the usage of fruits and consumer preference. Thus, modification of fruit shape has become one of the major objectives for crop improvement. However, the underlying mechanisms of fruit shape regulation are poorly understood. In this review we summarize recent progress in the genetic basis of fleshy fruit shape regulation using tomato, cucumber, and peach as examples. Comparative analyses suggest that the OFP-TRM (OVATE Family Protein - TONNEAU1 Recruiting Motif) and IQD (IQ67 domain) pathways are probably conserved in regulating fruit shape by primarily modulating cell division patterns across fleshy fruit species. Interestingly, cucumber homologs of FRUITFULL (FUL1), CRABS CLAW (CRC) and 1-aminocyclopropane-1-carboxylate synthase 2 (ACS2) were found to regulate fruit elongation. We also outline the recent progress in fruit shape regulation mediated by OFP-TRM and IQD pathways in Arabidopsis and rice, and propose that the OFP-TRM pathway and IQD pathway coordinate regulate fruit shape through integration of phytohormones, including brassinosteroids, gibberellic acids, and auxin, and microtubule organization. In addition, functional redundancy and divergence of the members of each of the OFP, TRM, and IQD families are also shown. This review provides a general overview of current knowledge in fruit shape regulation and discusses the possible mechanisms that need to be addressed in future studies.

A Novel Introgression Line Library Derived from a Wild Melon Gives Insights into the Genetics of Melon Domestication, Uncovering New Genetic Variability Useful for Breeding

A collection of 30 melon introgression lines (ILs) was developed from the wild accession Ames 24297 (TRI) into 'Piel de Sapo' (PS) genetic background. Each IL carried an average of 1.4 introgressions from TRI, and the introgressions represented 91.4% of the TRI genome. Twenty-two ILs, representing 75% of the TRI genome, were evaluated in greenhouse (Algarrobo and Meliana) and field (Alcàsser) trials, mainly to study traits related to domestication syndrome such as fruit weight (FW) and flesh content (FFP), as well as other fruit quality traits as fruit shape (FS), flesh firmness (FF), soluble solid concentration (SSC), rind color and abscission layer. The IL collection showed an impressive variation in size-related traits, with FW ranging from 800 to 4100 g, reflecting the strong effect of the wild genome on these traits. Most of the ILs produced smaller fruits compared with PS; however, unexpectedly, the IL TRI05-2 produced bigger fruits, likely due to new epistatic interacions with the PS genetic background. In contrast, the genotypic effect for FS was smaller, and few QTLs with notable effects were detected. Interestingly, variability was also observed for FFP, FF and SSC, rind color and abscission layer formation. Genes in these introgressions are candidates for having been involved in melon domestication and diversification as well. These results confirm that the TRI IL collection is a very powerful tool for mapping traits of agronomic interest in melon, allowing the confirmation of previously reported QTLs and the identification of new ones to better understand the domestication process of this crop.

Integrative analysis of grapevine (Vitis vinifera L) transcriptome reveals regulatory network for Chardonnay quality formation

Anthocyanins, total phenols, soluble sugar and fruit shape plays a significant role in determining the distinct fruit quality and customer preference. However, for the majority of fruit species, little is known about the transcriptomics and underlying regulatory networks that control the generation of overall quality during fruit growth and ripening. This study incorporated the quality-related transcriptome data from 6 ecological zones across 3 fruit development and maturity phases of Chardonnay cultivars. With the help of this dataset, we were able to build a complex regulatory network that may be used to identify important structural genes and transcription factors that control the anthocyanins, total phenols, soluble sugars and fruit shape in grapes. Overall, our findings set the groundwork to improve grape quality in addition to offering novel views on quality control during grape development and ripening.

QTL Mapping for Ovary- and Fruit-Related Traits in Cucumis sativus-C. hystrix Introgression Line IL52

IL52 is a valuable introgression line obtained from interspecific hybridization between cultivated cucumber (Cucumis sativus L., 2n = 14) and the wild relative species C. hystrix Chakr. (2n = 24). IL52 exhibits high resistance to a number of diseases, including downy mildew, powdery mildew, and angular leaf spot. However, the ovary- and fruit-related traits of IL52 have not been thoroughly investigated. Here, we conducted quantitative trait loci (QTL) mapping for 11 traits related to ovary size, fruit size, and flowering time using a previously developed 155 F_7:8 RIL population derived from a cross between CCMC and IL52. In total, 27 QTL associated with the 11 traits were detected, distributed on seven chromosomes. These QTL explained 3.61% to 43.98% of the phenotypic variance. Notably, we identified a major-effect QTL (qOHN4.1) on chromosome 4 associated with the ovary hypanthium neck width and further delimited it into a 114-kb candidate region harboring 13 candidate genes. Furthermore, the QTL qOHN4.1 is co-localized with the QTL detected for ovary length, mature fruit length, and fruit neck length, all residing within the consensus QTL FS4.1, suggesting a plausible pleiotropic effect.

Comparative Anatomical and Transcriptomics Reveal the Larger Cell Size as a Major Contributor to Larger Fruit Size in Apricot

Fruit size is one of the essential quality traits and influences the economic value of apricots. To explore the underlying mechanisms of the formation of differences in fruit size in apricots, we performed a comparative analysis of anatomical and transcriptomics dynamics during fruit growth and development in two apricot cultivars with contrasting fruit sizes (large-fruit Prunus armeniaca 'Sungold' and small-fruit P. sibirica 'F43'). Our analysis identified that the difference in fruit size was mainly caused by the difference in cell size between the two apricot cultivars. Compared with 'F43', the transcriptional programs exhibited significant differences in 'Sungold', mainly in the cell expansion period. After analysis, key differentially expressed genes (DEGs) most likely to influence cell size were screened out, including genes involved in auxin signal transduction and cell wall loosening mechanisms. Furthermore, weighted gene co-expression network analysis (WGCNA) revealed that PRE6/bHLH was identified as a hub gene, which interacted with 1 TIR1, 3 AUX/IAAs, 4 SAURs, 3 EXPs, and 1 CEL. Hence, a total of 13 key candidate genes were identified as positive regulators of fruit size in apricots. The results provide new insights into the molecular basis of fruit size control and lay a foundation for future breeding and cultivation of larger fruits in apricot.

Mapping and identification of CsSF4, a gene encoding a UDP-N-acetyl glucosamine-peptide N-acetylglucosaminyltransferase required for fruit elongation in cucumber (Cucumis sativus L.)

The short fruit length phenotype in sf4 is caused by a SNP in Csa1G665390, which encodes an O-linked N-acetylglucosamine (GlcNAc) transferase in cucumber. Cucumber fruit is an excellent resource for studying fruit morphology due to its fast growth rate and naturally abundant morphological variations. The regulatory mechanisms underlying plant organ size and shape are important and fundamental biological questions. In this study, a short-fruit length mutant, sf4, was identified from an ethyl methanesulfonate (EMS) mutagenesis population derived from the North China-type cucumber inbred line WD1. Genetic analysis indicated that the short fruit length phenotype of sf4 was controlled by a recessive nuclear gene. The SF4 locus was located in a 116.7-kb genomic region between the SNP markers GCSNP75 and GCSNP82 on chromosome 1. Genomic and cDNA sequences analysis indicated that a single G to A transition at the last nucleotide of Csa1G665390 intron 21 in sf4 changed the splice site from GT-AG to GT-AA, resulting in a 42-bp deletion in exon 22. Csa1G665390 is presumed to be a candidate gene, CsSF4 that encodes an O-linked N-acetylglucosamine (GlcNAc) transferase (OGT). CsSF4 was highly expressed in the leaves and male flowers of wild-type cucumbers. Transcriptome analysis indicated that sf4 had alterations in expression of many genes involved in hormone response pathways, cell cycle regulation, DNA replication, and cell division, suggesting that cell proliferation-associated gene networks regulate fruit development in cucumber. Identification of CsSF4 will contribute to elucidating the function of OGT in cell proliferation and to understanding fruit elongation mechanisms in cucumber.

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.

Natural variation in CRABS CLAW contributes to fruit length divergence in cucumber

Fruit length is a key domestication trait that affects crop yield and appearance. Cucumber (Cucumis sativus) fruits vary from 5 to 60 cm in length. Despite the identification of several regulators and multiple quantitative trait loci (QTLs) underlying fruit length, the natural variation, and molecular mechanisms underlying differences in fruit length are poorly understood. Through map-based cloning, we identified a nonsynonymous polymorphism (G to A) in CRABS CLAW (CsCRC) as underlying the major-effect fruit size/shape QTL FS5.2 in cucumber. The short-fruit allele CsCRCA is a rare allele that has only been found in round-fruited semi-wild Xishuangbanna cucumbers. A near-isogenic line (NIL) homozygous for CsCRCA exhibited a 34∼39% reduction in fruit length. Introducing CsCRCG into this NIL rescued the short-fruit phenotype, and knockdown of CsCRCG resulted in shorter fruit and smaller cells. In natural cucumber populations, CsCRCG expression was positively correlated with fruit length. Further, CsCRCG, but not CsCRCA, targets the downstream auxin-responsive protein gene CsARP1 to regulate its expression. Knockout of CsARP1 produced shorter fruit with smaller cells. Hence, our work suggests that CsCRCG positively regulates fruit elongation through transcriptional activation of CsARP1 and thus enhances cell expansion. Using different CsCRC alleles provides a strategy to manipulate fruit length in cucumber breeding.

Translational research from models to crops: comparative genomics for plant breeding

The concept of translational research, which originated in the medical field in the 1980s, consists in improving the efficient transfer of research results obtained in a species (which can be considered as a model or pivot) to all the species for which these results are of interest for its improvement in Agriculture. In this context, comparative genomics is an important tool for translational research, effectively identifying genes controlling common functions between species. Editing and phenotyping tools must thus allow the functional validation of the gene conserved within the species for which the knowledge has been extrapolated, that is to say transferred, and the identification of the best alleles and associated genotypes for exploitation in current breeding programs.

Comparative transcriptomic analysis of early fruit development in eggplant (Solanum melongena L.) and functional characterization of SmOVATE5

Comparative transcriptome analysis of early fruits of long and round eggplants, SmOVATE5, is involved in regulating fruit development. Eggplant, a solanaceous crop that has undergone a long period of domestication, is one of the most important vegetables worldwide. The shape of its fruit is an important agronomic trait and consumers in different regions have different preferences. However, a limited understanding of the molecular mechanisms regulating fruit development and shape has hindered eggplant breeding. In this study, we performed morphological observations and transcriptome analysis of long- and round-fruited eggplant genotypes to understand the molecular regulation during the early development of different fruit shapes. Morphological studies revealed that the two varieties already exhibited distinctly different phenotypes at the initial stage of fruit development before flowering, with rapid fruit enlargement beginning on the sixth day after flowering. Comparative transcriptome analysis identified phytohormone-related genes that were significantly upregulated on the day of flowering, indicating they may be involved in regulating the initial stages of fruit development. Notably, SmARF1 showed a sustained upregulation pattern in both varieties, suggesting that it may promote eggplant fruit growth. In addition, several differentially expressed genes of the SUN, YABBY, and OVATE families are potentially involved in the regulation of fruit development or fruit shape. We demonstrated that the SmOVATE5 gene has a negative regulatory function suppressing plant growth and development. In conclusion, this study provides new insights into the molecular regulatory mechanisms of eggplant fruit development, and the genes identified may provide valuable references for different fruit shape breeding programs.

Variation in Shoot, Peduncle and Fruit Growth of Lagenaria siceraria Landraces

Lagenaria siceraria (Molina) Standley is a prominent food source as almost all its plant parts are edible. However, no studies have recorded the changes in shoots, peduncles and fruits during its growth. Hence, this study aimed to record changes in shoot traits and relate the peduncle to the fruit traits of L. siceraria landraces across different growth stages. Changes in shoots, peduncles and fruits during growth were compared within and among landraces using analysis of variance, correlation, principal component analysis, cluster analysis and heritability estimates. Almost all landraces had harvestable shoots at 42 days after sowing. Peduncles became shorter and wider as the fruits elongated. Shoots, peduncles and fruits correlated positively with each other. The informative principal components had a total variability of 84.488%, with a major contribution from shoot traits. The biplot and dendrogram clustered landraces with similar growth habits and the harvestable shoot and fruit attributes into three clusters, but KRI and NSRC formed singlets. Shoot width (60.2%) and peduncle length (55.2%) had high heritability estimates. The general low heritability estimates and genetic advances indicated the presence of non-additive gene action. This study is the first report on changes in harvested shoots and the relationship between peduncles and fruits during growth.

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.

BSR and Full-Length Transcriptome Approaches Identified Candidate Genes for High Seed Ratio in Camellia vietnamensis

(1) Background: C. vietnamensis is very suitable for growth in the low hilly areas of southern subtropical regions. Under appropriate conditions, the oil yield of C. vietnamensis can reach 1125 kg/ha (the existing varieties can reach 750 kg/ha). Moreover, the fruit of C. vietnamensis is large and the pericarp is thick (>5 cm). Therefore, a high seed ratio has become the main target economic trait for the breeding of C. vietnamensis. (2) Methods: A half-sibling population of C. vietnamensis plants with a combination of high and low seed ratios was constructed by crossing a C. vietnamensis female parent. Bulked segregant RNA analysis and full-length transcriptome sequencing were performed to determine the molecular mechanisms underlying a high seed ratio. (3) Results: Seed ratio is a complex quantitative trait with a normal distribution, which is significantly associated with four other traits of fruit (seed weight, seed number, fruit diameter, and pericarp thickness). Two candidate regions related to high seed ratio (HSR) were predicted. One spanned 140.8−148.4 Mb of chromosome 2 and was associated with 97 seed-yield-related candidate genes ranging in length from 278 to 16,628 bp. The other spanned 35.3−37.3 Mb on chromosome 15 and was associated with 38 genes ranging in length from 221 to 16,928 bp. Using the full-length transcript as a template, a total of 115 candidate transcripts were obtained, and 78 transcripts were predicted to be functionally annotated. The DEGs from two set pairs of cDNA sequencing bulks were enriched to cytochrome p450 CYP76F14 (KOG0156; GO:0055114, HSR4, HSR7), the gibberellin phytohormone pathway (GO:0016787, HSR5), the calcium signaling pathway (GO:0005509, HSR6), the polyubiquitin-PPAR signaling pathway (GO:0005515, HSR2, HSR3), and several main transcription factors (bZIP transcription factor, HSR1) in C. vietnamensis.

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.

Genome-Wide Identification of the SUN Gene Family in Melon (Cucumis melo) and Functional Characterization of Two CmSUN Genes in Regulating Fruit Shape Variation

Melon (Cucumis melo) is an important economic crop cultivated worldwide. A unique SUN gene family plays a crucial role in regulating plant growth and fruit development, but many SUN family genes and their function have not been well-characterized in melon. In the present study, we performed genome-wide identification and bioinformatics analysis and identified 24 CmSUN family genes that contain integrated and conserved IQ67 domain in the melon genome. Transcriptome data analysis and qRT-PCR results showed that most CmSUNs are specifically enriched in melon reproductive organs, such as young flowers and ovaries. Through genetic transformation in melons, we found that overexpression of CmSUN23-24 and CmSUN25-26-27c led to an increased fruit shape index, suggesting that they act as essential regulators in melon fruit shape variation. Subcellular localization revealed that the CmSUN23-24 protein is located in the cytoplasmic membrane. A direct interaction between CmSUN23-24 and a Calmodulin protein CmCaM5 was found by yeast two-hybrid assay, which indicated their participation in the calcium signal transduction pathway in regulating plant growth. These findings revealed the molecular characteristics, expression profile, and functional pattern of the CmSUN genes, and may provide the theoretical basis for the genetic improvement of melon fruit breeding.

Genome-wide association study provides genetic insights into natural variation in watermelon rind thickness and single fruit weight

Rind thickness and fruit weight are agronomic traits closely related to quality and yield, which have attracted much attention from consumers and breeders. However, the genetic mechanism of these two traits is still not well understood in natural populations. In this study, rind thickness and single fruit weight in 151 watermelon accessions were determined in 2019 and 2020, and genome-wide association analysis was performed by integrating phenotypic and genotype data. Abundant phenotypic variation was found in the test population, and the watermelon with thinner rind thickness tended to have smaller fruit weights. Five significant SNPs were closely associated with rind thickness on chromosome 2 by Genome-wide association study (GWAS), i.e., 32344170, 32321308, 32304738, 32328501, and 32311192. And there were 21 genes were annotated in the candidate interval, most notably, Cla97C02G044160 belonged to the MADS family, and part of the genes in this family played an important role in regulating organ size, further analysis of gene structure, gene expression level, and phylogenetic tree showed that Cla97C02G044160 was a candidate gene for regulating target traits. In conclusion, our study provides molecular insights into the natural variation of watermelon rind thickness and single fruit weight, meanwhile, providing data support for molecular marker-assisted breeding.

Construction of a density mutant collection in bitter gourd via new germplasms innovation and gene functional study

Although a few studies have elucidated the creation of bitter gourd mutants, the suitable concentration and duration of ethyl methanesulfonate (EMS) mutagenesis have not been determined. In this study, mutant collection was conducted to create new germplasms and widen genetic diversity. By employing the seeds of the inbred line Y52 as the mutagenic material, EMS as the mutagen, and the suitable mutagenic conditions for bitter gourd seeds (EMS concentration 0.2%, mutagenic time 10 h), we mutated 10,000 seeds and acquired 3223 independent M₁ lines. For the randomly selected 1000 M₂ lines, 199 M₂ lines with visible phenotypes were found, and 167 M₂ lines were mutants of fruit shape, size, and tubercles. Furthermore, fourteen dwarf, eleven leaf color, five leaf shape, and eight meristem defect mutants were discovered in this mutant collection. In addition, three lines of 1253, 2284, and 3269 represented recessive mutants crossed with Y52. Furthermore, the yellow leaf lines of 2284 and 3269 were not mutated at the same gene locus. This study constructed a mutant collection through innovative new germplasms and provided valuable resources for bitter gourd breeding and functional gene research.

QTL mapping reveals candidate genes for main agronomic traits in Luffa based on a high-resolution genetic map

Luffa is an important medicinal and edible vegetable crop of Cucurbitaceae. Strong heterosis effects and strikingly complementary characteristics were found between the two domesticated Luffa cultivars, Luffa acutangula and Luffa cylindrica. To explore the genetic basis underlying their important agronomic traits, we constructed the first interspecific high-density genetic linkage map using a BC₁ population of 110 lines derived from a cross between S1174 (Luffa acutangula) and P93075 (Luffa cylindrica). The map spanned a total of 2246.74 cM with an average distance of 0.48 cM between adjacent markers. Thereafter, a large-scale field-based quantitative trait loci (QTLs) mapping was conducted for 25 important agronomic traits and 40 significant genetic loci distributed across 11 chromosomes were detected. Notably, a vital QTL (qID2) located on chromosome 9 with a minimum distance of 23 kb was identified to be responsible for the internode diameter and explained 11% of the phenotypic variation. Lac09g006860 (LacCRWN3), encoding a nuclear lamina protein involved in the control of nuclear morphology, was the only gene harbored in qID2. Sequence alignment showed completely different promoter sequences between the two parental alleles of LacCRWN3 except for some nonsynonymous single nucleotide polymorphisms (SNPs) in exons, and the expression level in thick-stem P93075 was distinctively higher than that in thin-stem S1174. According to the natural variation analysis of a population of 183 inbred lines, two main haplotypes were found for LacCRWN3: the P93075-like and S1174-like, with the former haplotype lines exhibiting significantly thicker internode diameters than those of the latter haplotype lines. It showed that LacCRWN3, as the only CRWN3 gene in Cucurbitaceae, was the most likely candidate gene regulating the internode diameter of Luffa. Our findings will be beneficial for deciphering the molecular mechanism of key phenotypic traits and promoting maker-assisted breeding in Luffa.

Genome-Wide Identification of the SAUR Gene Family in Wax Gourd (Benincasa hispida) and Functional Characterization of BhSAUR60 during Fruit Development

The wax gourd (Benincasa hispida) is an important vegetable crop whose fruits contain nutrients and metabolites. Small auxin upregulated RNA (SAUR) genes constitute the largest early auxin-responsive gene family and regulate various biological processes in plants, although this gene family has not been studied in the wax gourd. Here, we performed genome-wide identification of the SAUR gene family in wax gourds and analyzed their syntenic and phylogenetic relationships, gene structures, conserved motifs, cis-acting elements, and expression patterns. A total of 68 SAUR (BhSAUR) genes were identified, which were distributed on nine chromosomes with 41 genes in two clusters. More than half of the BhSAUR genes were derived from tandem duplication events. The BhSAUR proteins were classified into seven subfamilies. BhSAUR gene promoters contained cis-acting elements involved in plant hormone and environmental signal responses. Further expression profiles showed that BhSAUR genes displayed different expression patterns. BhSAUR60 was highly expressed in fruits, and overexpression led to longer fruits in Arabidopsis. In addition, the plants with overexpression displayed longer floral organs and wavy stems. In conclusion, our results provide a systematic analysis of the wax gourd SAUR gene family and facilitate the functional study of BhSAUR60 during wax gourd fruit development.

Phenotypic Characterization and Fine Mapping of a Major-Effect Fruit Shape QTL FS5.2 in Cucumber, Cucumis sativus L., with Near-Isogenic Line-Derived Segregating Populations

Cucumber (Cucumis sativus L.) fruit size/shape (FS) is an important yield and quality trait that is quantitatively inherited. Many quantitative trait loci (QTLs) for fruit size/shape have been identified, but very few have been fine-mapped or cloned. In this study, through marker-assisted foreground and background selections, we developed near-isogenic lines (NILs) for a major-effect fruit size/shape QTL FS5.2 in cucumber. Morphological and microscopic characterization of NILs suggests that the allele of fs5.2 from the semi-wild Xishuangbanna (XIS) cucumber (C. s. var. xishuangbannesis) reduces fruit elongation but promotes radial growth resulting in shorter but wider fruit, which seems to be due to reduced cell length, but increased cellular layers. Consistent with this, the NIL carrying the homozygous XIS allele (fs5.2) had lower auxin/IAA contents in both the ovary and the developing fruit. Fine genetic mapping with NIL-derived segregating populations placed FS5.2 into a 95.5 kb region with 15 predicted genes, and a homolog of the Arabidopsis CRABS CLAW (CsCRC) appeared to be the most possible candidate for FS5.2. Transcriptome profiling of NIL fruits at anthesis identified differentially expressed genes enriched in the auxin biosynthesis and signaling pathways, as well as genes involved in cell cycle, division, and cell wall processes. We conclude that the major-effect QTL FS5.2 controls cucumber fruit size/shape through regulating auxin-mediated cell division and expansion for the lateral and longitudinal fruit growth, respectively. The gibberellic acid (GA) signaling pathway also plays a role in FS5.2-mediated fruit elongation.

An allelic variant in the ACS7 gene promotes primary root growth in watermelon

Gene mining in a C. lanatus × C. amarus population revealed one gene, ACS7, linked to primary root elongation in watermelon. Watermelon is a xerophytic crop characterized by a long primary root and robust lateral roots. Therefore, watermelon serves as an excellent model for studying root elongation and development. However, the genetic mechanism underlying the primary root elongation in watermelon remains unknown. Herein, through bulk segregant analysis we identified a genetic locus, qPRL.Chr03, controlling primary root length (PRL) using two different watermelon species (Citrullus lanatus and Citrullus amarus) that differ in their root architecture. Fine mapping revealed that xaa-Pro dipeptidase and 1-aminocyclopropane-1-carboxylate synthase 7 (ACS7) are candidate regulators of the primary root growth. Allelic variation in the delimited region among 193 watermelon accessions indicated that the long-root alleles might only exist in C. amarus. Interestingly, the discrepancy in PRL among the C. amarus accessions was clearly associated with a nonsynonymous single nucleotide polymorphism variant within the ACS7 gene. The ACS7 expression and ethylene levels in the primary root tips suggested that ethylene is a negative regulator of root elongation in watermelon, as supported by the application of 1-aminocyclopropane-1-carboxylate (ACC, the ethylene precursor) or 2-aminoethoxyvinyl glycine (AVG, an ACS inhibitor). To the best of our knowledge, these findings provide the first description of the genetic basis of root elongation in watermelon. The detected markers of the ACS7 gene will facilitate marker-assisted selection for the PRL trait to improve water and nutrient use efficacy in watermelon and beyond.