Transcriptome analyses of early cucumber fruit growth identifies distinct gene modules associated with phases of development

Parthenocarpy, a pollination-independent fruit set mechanism to ensure yield stability

Fruit development is essential for flowering plants' reproduction and a significant food source. Climate change threatens fruit yields due to its impact on pollination and fertilization processes, especially vulnerable to extreme temperatures, insufficient light, and pollinator decline. Parthenocarpy, the development of fruit without fertilization, offers a solution, ensuring yield stability in adverse conditions and enhancing fruit quality. Parthenocarpic fruits not only secure agricultural production but also exhibit improved texture, appearance, and shelf life, making them desirable for food processing and other applications. Recent research unveils the molecular mechanisms behind parthenocarpy, implicating transcription factors (TFs), noncoding RNAs, and phytohormones such as auxin, gibberellin (GA), and cytokinin (CK). Here we review recent findings, construct regulatory models, and identify areas for further research.

RNA-seq-based transcriptome profiling of early fruit development in Chieh-qua and analysis of related transcription factors

Chieh-qua (Benincasa hispida Cogn. var. Chieh-qua How.) fruit development starts post pollination. With the continuous expansion of the fruit, the soluble solid content of the fruit decreases. Because there are no reports on the early development of Chieh-qua fruit, this study compared fruit transcriptomes at 0-, 3-, and 7 day post pollination (dpp). 104,747 unigenes were assembled from clean reads and compared using six public databases for similarity searching. Compared with those of 0 dpp (C), there were differences in the expression of 12,982 and 6541 genes in the fruit tissue at 3 dpp and 7 dpp, respectively. Compared with 3 dpp (B), there were 14,314 differentially expressed genes in the fruit at 7 dpp (A). Based on the analysis of transcription factors, 213 nucleotides in the MYB superfamily were identified; among them, 94 unigenes of the MYB superfamily were differentially expressed at the three stages. In the pairwise comparison of differential expression, eight unigenes (Gene_id: TRINITY_DN32880_c1_g2, TRINITY_DN35142_c2_g2, TRINITY_DN32454_c11_g6, TRINITY_DN34105_c2_g7, TRINITY_DN32758_c3_g3, TRINITY_DN33604_c4_g10, TRINITY_DN34466_c3_g1, TRINITY_DN35924_c3_g2) were homologous to those of MYB59, MYB-GT3b, MYB18, MYB4, MYB108, MYB306, MYB340, and MYB-bHLH13. These unigenes differed significantly among the three stages. Furthermore, MYB59 and MYB18 exhibited higher expression at 7 dpp. MYB4, MYB-GT3b, MYB108, and MYB306 showed the highest expression levels in fruits at 3 dpp. In addition, MYB340 and MYB-bHLH13 showed higher expression levels during the unpollinated stage. MYB59, MYB-GT3b, MYB18, MYB4, MYB108, MYB306, MYB340, and MYB-bHLH13 may play crucial roles in Chieh-qua fruit development, defense, and blossoming. This study provides a basis for further investigation of MYB superfamily genes involved in early fruit expansion in chieh-qua.

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

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

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.

Characterization of Transcriptome Dynamics during Early Fruit Development in Olive (Olea europaea L.)

In the olive (Olea europaea L.), an economically leading oil crop worldwide, fruit size and yield are determined by the early stages of fruit development. However, few detailed analyses of this stage of fruit development are available. This study offers an extensive characterization of the various processes involved in early olive fruit growth (cell division, cell cycle regulation, and cell expansion). For this, cytological, hormonal, and transcriptional changes characterizing the phases of early fruit development were analyzed in olive fruit of the cv. 'Picual'. First, the surface area and mitotic activity (by flow cytometry) of fruit cells were investigated during early olive fruit development, from 0 to 42 days post-anthesis (DPA). The results demonstrate that the cell division phase extends up to 21 DPA, during which the maximal proportion of 4C cells in olive fruits was reached at 14 DPA, indicating that intensive cell division was activated in olive fruits at that time. Subsequently, fruit cell expansion lasted as long as 3 weeks more before endocarp lignification. Finally, the molecular mechanisms controlling the early fruit development were investigated by analyzing the transcriptome of olive flowers at anthesis (fruit set) as well as olive fruits at 14 DPA (cell division phase) and at 28 DPA (cell expansion phase). Sequential induction of the cell cycle regulating genes is associated with the upregulation of genes involved in cell wall remodeling and ion fluxes, and with a shift in plant hormone metabolism and signaling genes during early olive fruit development. This occurs together with transcriptional activity of subtilisin-like protease proteins together with transcription factors potentially involved in early fruit growth signaling. This gene expression profile, together with hormonal regulators, offers new insights for understanding the processes that regulate cell division and expansion, and ultimately fruit yield and olive size.

Telomere-to-telomere genome assembly of bitter melon (Momordica charantia L. var. abbreviata Ser.) reveals fruit development, composition and ripening genetic characteristics

Momordica charantia L. var. abbreviata Ser. (Mca), known as bitter gourd or bitter melon, is a Momordica variety with medicinal value and belongs to the Cucurbitaceae family. In view of the lack of genomic information on bitter gourd and other Momordica species and to promote Mca genomic research, we assembled a 295.6-Mb telomere-to-telomere (T2T) high-quality Mca genome with six gap-free chromosomes after Hi-C correction. This genome is anchored to 11 chromosomes, which is consistent with the karyotype information, and comprises 98 contigs (N50 of 25.4 Mb) and 95 scaffolds (N50 of 25.4 Mb). The Mca genome harbors 19 895 protein-coding genes, of which 45.59% constitute predicted repeat sequences. Synteny analysis revealed variations involved in fruit quality during the divergence of bitter gourd. In addition, assay for transposase-accessible chromatin by high-throughput sequencing and metabolic analysis showed that momordicosides and other substances are characteristic of Mca fruit pulp. A combined transcriptomic and metabolomic analysis revealed the mechanisms of pigment accumulation and cucurbitacin biosynthesis in Mca fruit peels, providing fundamental molecular information for further research on Mca fruit ripening. This report provides a new genetic resource for Momordica genomic studies and contributes additional insights into Cucurbitaceae phylogeny.

Chayote Fruit (Sechium edule var. virens levis) Development and the Effect of Growth Regulators on Seed Germination

The chayote fruit is a nontraditional vegetable belonging to the Cucurbitaceae family. The fruit has an endocarpic recalcitrant seed that emerges postharvest, drastically shortening its shelf life. In this study, the changes during fruit and seed development before and after harvest (ah) are reported. Additionally, in order to investigate how growth regulators (GRs) affect seed germination, 2-cloroethylphosphonic acid (CPA) (200 µL L^-1), gibberellic acid (GA₃) (100 and 200 mg L^-1), auxin (2,4-D) (0.5 and 1.0 mM), and abscisic acid (ABA) (0.5 and 1.0 mM) were applied after harvest. The results showed that the chayote fruit reached horticultural maturity at 21 days after anthesis, with a sigmoid trend: phase I featured slow growth and high transpiration; in phase II, growth was accelerated and accumulation of endosperm was observed; and in phase III, both growth rate and transpiration were reduced, soluble sugars increased, and the seed showed 25% cotyledon development. At day 13 ah, CPA, GA₃, and 2,4-D (0.5 mM) increased seed germination, with values between 10 and 15 mm of the embryonary axis, and the treatments with 2,4-D (1 mM) and ABA (0.5 and 1.0 mM) retarded their growth (2-6 mm). This research allowed us to reveal the phenological phases and the shelf life of the chayote fruit, as well as the results of possible postharvest treatment with GRs; our results suggest that strategies to delay viviparism and prolong the shelf life of the fruit should be applied before 10 days ah, when the embryonic axis of the seed has not developed.

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.

Research Progress on the Leaf Morphology, Fruit Development and Plant Architecture of the Cucumber

Cucumber (Cucumis sativus L.) is an annual climbing herb that belongs to the Cucurbitaceae family and is one of the most important economic crops in the world. The breeding of cucumber varieties with excellent agronomic characteristics has gained more attention in recent years. The size and shape of the leaves or fruit and the plant architecture are important agronomic traits that influence crop management and productivity, thus determining the crop yields and consumer preferences. The growth of the plant is precisely regulated by both environmental stimuli and internal signals. Although significant progress has been made in understanding the plant morphological regulation of Arabidopsis, rice, and maize, our understanding of the control mechanisms of the growth and development of cucumber is still limited. This paper reviews the regulation of phytohormones in plant growth and expounds the latest progress in research regarding the genetic regulation pathways in leaf development, fruit size and shape, branching, and plant type in cucumber, so as to provide a theoretical basis for improving cucumber productivity and cultivation efficiency.

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.

Gene Expression in Zucchini Fruit Development

The study of fruit development in zucchini via gene expression has proven to be applicable in breeding programs. Phenotypic and transcriptomic studies of fruit set and parthenocarpy have been previously developed and some relevant genes have been reported. From these studies, three genotypes (MUCU-16, Whitaker, and Cavili) and six genes (CpAUX22, CpIAA4, CpIAMT-1, CpPIN5, CpCYCD6-1, and CpEXPLB1) were selected. The expression of these genes was analyzed in each genotype under three different treatments (pollination, auxin-treatment and non-treatment) during one week post anthesis. Also, a phenotyping analysis was conducted. The different nature of the samples and the genes selected allowed associations between different fruit traits and fruit development stages. There was a rapid response of CpAUX22 and CpIAA4 to the auxin treatment. Also, these genes and the CpIAMT-1 became more overexpressed in pollinated samples over time. The CpPIN5 gene increased its expression over time in all genotypes while CpCYCD6-1 was overexpressed in the early stages of fruit development in all samples. The CpEXPLB1 was highly up-regulated in non-treated samples, suggesting a relationship with fruit abortion. The overexpression of CpAUX22 and the non-overexpression of CpEXPLB1 in early stages may be associated with fruit growth in zucchini.

Cucurbitaceae genome evolution, gene function and molecular breeding

The Cucurbitaceae is one of the most genetically diverse plant families in the world. Many of them are important vegetables or medicinal plants and are widely distributed worldwide. The rapid development of sequencing technologies and bioinformatic algorithms has enabled the generation of genome sequences of numerous important Cucurbitaceae species. This has greatly facilitated research on gene identification, genome evolution, genetic variation and molecular breeding of cucurbit crops. So far, genome sequences of 18 different cucurbit species belonging to tribes Benincaseae, Cucurbiteae, Sicyoeae, Momordiceae and Siraitieae have been deciphered. This review summarizes the genome sequence information, evolutionary relationship, and functional genes associated with important agronomic traits (e.g., fruit quality). The progress of molecular breeding in cucurbit crops and prospects for future applications of Cucurbitaceae genome information are also discussed.

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

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

Identification and Expression Analysis of miR160 and Their Target Genes in Cucumber

miR160 plays a crucial role in various biological processes by regulating their target gene auxin response factor (ARF) in plants. However, little is known about miR160 and ARF in cucumber fruit expansion. Here, 4 Csa-MIR160 family members and 17 CsARFs were identified through a genome-wide search. Csa-miR160 showed a closer relationship with those in melon. Phylogenetic analysis revealed that CsARFs were divided into four classes and most of CsARFs presented a closer evolutionary relationship with those from tomato. Putative cis-elements analysis predicted that Csa-MIR160 and CsARFs were involved in light, phytohormone and stress response, which proved that they might take part in light, phytohormone and stress signaling pathway by the miR160-ARF module. In addition, CsARF5, CsARF11, CsARF13 and CsARF14 were predicted as the target genes of Csa-miR160. qRT-PCR revealed that Csa-miR160 and their target gene CsARFs were differentially expressed in differential cucumber tissues and developmental stages. Csa-miR160d was only expressed in the expanded cucumber fruit. CsARF5, CsARF11 and CsARF13 exhibited the lower expression in the expanded fruit than those in the ovary, while, CsARF14 showed the reverse trend. Our results suggested that Csa-miR160d might play a crucial role in cucumber fruit expansion by negatively targeting CsARF5, CsARF11 and CsARF13. This is the first genome-wide analysis of miR160 in cucumber. These findings provide useful information and resources for further studying the role of miR160 and ARF in cucumber fruit expansion.

Hexose transporter CsSWEET7a in cucumber mediates phloem unloading in companion cells for fruit development

In the fleshy fruit of cucumbers (Cucumis sativus L.), the phloem flow is unloaded via an apoplasmic pathway, which requires protein carriers to export sugars derived from stachyose and raffinose into the apoplasm. However, transporter(s) involved in this process remain unidentified. Here, we report that a hexose transporter, CsSWEET7a (Sugar Will Eventually be Exported Transporter 7a), was highly expressed in cucumber sink tissues and localized to the plasma membrane in companion cells of the phloem. Its expression level increased gradually during fruit development. Down-regulation of CsSWEET7a by RNA interference (RNAi) resulted in smaller fruit size along with reduced soluble sugar levels and reduced allocation of 14C-labelled carbon to sink tissues. CsSWEET7a overexpression lines showed an opposite phenotype. Interestingly, genes encoding alkaline α-galactosidase (AGA) and sucrose synthase (SUS) were also differentially regulated in CsSWEET7a transgenic lines. Immunohistochemical analysis demonstrated that CsAGA2 co-localized with CsSWEET7a in companion cells, indicating cooperation between AGA and CsSWEET7a in fruit phloem unloading. Our findings indicated that CsSWEET7a is involved in sugar phloem unloading in cucumber fruit by removing hexoses from companion cells to the apoplasmic space to stimulate the raffinose family of oligosaccharides (RFOs) metabolism so that additional sugars can be unloaded to promote fruit growth. This study also provides a possible avenue towards improving fruit production in cucumber.

Genome wide analysis of kinesin gene family in Citrullus lanatus reveals an essential role in early fruit development

BACKGROUND

Kinesin (KIN) as a motor protein is a versatile nano-machine and involved in diverse essential processes in plant growth and development. However, the kinesin gene family has not been identified in watermelon, a valued and nutritious fruit, and yet their functions have not been characterized. Especially, their involvement in early fruit development, which directly determines the size, shape, yield and quality of the watermelon fruit, remains unclear.

RESULTS

In this study, we performed a whole-genome investigation and comprehensive analysis of kinesin genes in C. lanatus. In total, 48 kinesins were identified and categorized into 10 kinesin subfamilies groups based on phylogenetic analysis. Their uneven distribution on 11 chromosomes was revealed by distribution analysis. Conserved motif analysis showed that the ATP-binding motif of kinesins was conserved within all subfamilies, but not the microtubule-binding motif. 10 segmental duplication pairs genes were detected by the syntenic and phylogenetic approaches, which showed the expansion of the kinesin gene family in C. lanatus genome during evolution. Moreover, 5 ClKINs genes are specifically and abundantly expressed in early fruit developmental stages according to comprehensive expression profile analysis, implying their critical regulatory roles during early fruit development. Our data also demonstrated that the majority of kinesin genes were responsive to plant hormones, revealing their potential involvement in the signaling pathways of plant hormones.

CONCLUSIONS

Kinesin gene family in watermelon was comprehensively analyzed in this study, which establishes a foundation for further functional investigation of C. lanatus kinesin genes and provides novel insights into their biological functions. In addition, these results also provide useful information for understanding the relationship between plant hormone and kinesin genes in C. lanatus.

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

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

Developmentally regulated activation of defense allows for rapid inhibition of infection in age-related resistance to Phytophthora capsici in cucumber fruit

Background

Age-related resistance (ARR) is a developmentally regulated phenomenon conferring resistance to pathogens or pests. Although ARR has been observed in several host-pathogen systems, the underlying mechanisms are largely uncharacterized. In cucumber, rapidly growing fruit are highly susceptible to Phytophthora capsici but become resistant as they complete exponential growth. We previously demonstrated that ARR is associated with the fruit peel and identified gene expression and metabolomic changes potentially functioning as preformed defenses.

Results

Here, we compare the response to infection in fruit at resistant and susceptible ages using microscopy, quantitative bioassays, and weighted gene co-expression analyses. We observed strong transcriptional changes unique to resistant aged fruit 2–4 h post inoculation (hpi). Microscopy and bioassays confirmed this early response, with evidence of pathogen death and infection failure as early as 4 hpi and cessation of pathogen growth by 8–10 hpi. Expression analyses identified candidate genes involved in conferring the rapid response including those encoding transcription factors, hormone signaling pathways, and defenses such as reactive oxygen species metabolism and phenylpropanoid biosynthesis.

Conclusion

The early pathogen death and rapid defense response in resistant-aged fruit provide insight into potential mechanisms for ARR, implicating both pre-formed biochemical defenses and developmentally regulated capacity for pathogen recognition as key factors shaping age-related resistance.

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

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

Metabolic and transcriptomic analysis of two Cucurbita moschata germplasms throughout fruit development

BACKGROUND

Pumpkins (Cucurbita moschata; Cucurbitaceae) are valued for their fruits and seeds and are rich in nutrients. Carotenoids and sugar contents, as main feature of pumpkin pulp, are used to determine the fruit quality.

RESULTS

Two pumpkin germplasms, CMO-X and CMO-E, were analyzed regarding the essential quality traits such as dry weight, soluble solids, organic acids, carotenoids and sugar contents. For the comparison of fruit development in these two germplasms, fruit transcriptome was analyzed at 5 different developmental stages from 0 d to 40 d in a time course manner. Putative pathways for carotenoids biosynthesis and sucrose metabolism were developed in C. moschata fruit and homologs were identified for each key gene involved in the pathways. Gene expression data was found consistent with the accumulation of metabolites across developmental stages and also between two germplasms. PSY, PDS, ZEP, CRTISO and SUS, SPS, HK, FK were found highly correlated with the accumulation of carotenoids and sucrose metabolites, respectively, at different growth stages of C. moschata as shown by whole transcriptomic analysis. The results of qRT-PCR analysis further confirmed the association of these genes.

CONCLUSION

Developmental regulation of the genes associated with the metabolite accumulation can be considered as an important factor for the determination of C. moschata fruit quality. This research will facilitate the investigation of metabolic profiles in other cultivars.

Comparative Transcriptomic Analysis of Two Bottle Gourd Accessions Differing in Fruit Size

The bottle gourd (Lagenaria siceraria) is an important horticultural and medicinal crop with high nutritional value. This study aimed at examining the molecular regulation of fruit size in bottle gourd. We performed transcriptome sequencing of two bottle gourd cultivars differing in their fruit size. The average fruit length and weight of the cultivar Hang (39.48 cm/624.4 g) were higher than those of the cultivar USA (10.34 cm/152.8 g) at maturity. Transcriptome sequencing and assembly resulted in 89,347 unigenes. A total of 1250 differentially expressed genes (DEG) were found between the two cultivars, including 422 upregulated genes and 828 downregulated genes in Hang as compared to USA. Genes related to cell wall metabolism, phytohormones, cell cycle, and cell division showed significant differential expression between the two cultivars. DEGs encoding transcription factors (TF) from nine TF families were also identified. The ethylene response factor family was the most enriched among these families. Our study provides a basis for further investigations of the molecular regulation of fruit size in bottle gourd.

Transcriptome profiling reveals genes involved in spine development during CsTTG1-regulated pathway in cucumber (Cucumis sativus L.)

The cucumber (Cucumis sativus L.), a type of fleshy fruit, is covered with spines (multicellular trichomes), which have a crucial impact on the economic value of the crop. Previous studies have found that CsTTG1 plays important roles in the initiation and further differentiation of cucumber spines, but how spine formation is regulated at the molecular level by CsTTG1 remains poorly understood. In this study, we characterized a cucumber 35S:CsTTG1 transgenic T2 line, OE-2, which bears relatively large and long spines compared with the small and short spines of the wild type (WT). Phenotypic measurements and histological analyses revealed that this phenotypic change was attributed to significant increases in cell number and size. Comparison of ovary epidermis transcriptomes between OE-2 and WT by DGE (Digital Gene Expression) analysis identified 1241 differentially expressed genes, among which 712 genes were dramatically upregulated and 529 downregulated in the ovary epidermis of OE-2. XTH23 and Cyclin family genes were significantly activated in OE-2, and transcription factors (TFs) were found to participate in spine size regulation in OE-2. Further analyses confirmed that GA was implicated in the regulation of fruit spine development in cucumber. Thus, our study provides a foundation for dissecting the molecular regulatory networks of fruit spine control in cucumber.

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.

Transcriptome profiling reveals key genes related to astringency during cucumber fruit development

The goal of this study was to provide quantitative data on the catechin contents and underlying molecular regulatory mechanisms in cucumber during fruit development. The dynamic changes in the total catechin contents and RNA-seq-based transcriptome profiling of the flesh and peel of the cucumber cultivar 'YanBai', which is strongly astringent, were examined at three key developmental stages 3, 6 and 9 days post-pollination. The total catechin content decreased as cucumber fruit developed and was significantly lower in the flesh than in the peel. In total, 5092 and 4004 genes were found to be differently expressed in the peel and flesh, respectively. Based on a functional annotation, eight structural genes encode enzymes involved in the catechin biosynthesis pathway. Three genes encoding 4-coumarate-CoA ligases, two genes encoding chalcone isomerases, two genes encoding dihydroflavonol-4-reductase and one gene each encoding a phenylalanine ammonia-lyase, flavanone 3-hydroxylase and cinnamate 4-hydroxylase were identified as affecting the catechin content of cucumber. The transcriptome data also revealed the significance of transcription factors, including WD40-repeat proteins, MYB and bHLH, in regulating catechin biosynthesis. These findings help increase our understanding of the molecular mechanisms controlling catechin biosynthesis and astringency development in cucumber fruit.

Identification of miRNAs and Their Target Genes Involved in Cucumber Fruit Expansion Using Small RNA and Degradome Sequencing

Fruit expansion is an essential and very complex biological process. Regulatory roles of microRNAs (miRNAs) and miRNA-mRNA modules in the cucumber fruit expansion are not yet to be investigated. In this work, 1253 known and 1269 novel miRNAs were identified from nine cucumber fruit small RNA (sRNA) libraries through high-throughput sequencing. A total of 105 highly differentially expressed miRNAs were recognized in the fruit on five days post anthesis with pollination (EXP_5d) sRNA library. Further, expression patterns of 11 differentially expressed miRNAs were validated by quantitative real-time PCR (qRT-PCR). The expression patterns were similar to sRNAs sequencing data. Transcripts of 1155 sequences were predicted as target genes of differentially expressed miRNAs by degradome sequencing. Gene Ontology (GO) enrichment showed that these target genes were involved in 24 biological processes, 15 cell components and nine molecular functions. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis demonstrated that these target genes were significantly enriched in 19 pathways and the enriched KEGG pathways were associated with environmental adaptation, signal transduction and translation. Based on the functional prediction of miRNAs and target genes, our findings suggest that miRNAs have a potential regulatory role in cucumber fruit expansion by targeting their target genes, which provide important data for understanding the miRNA-mediated regulatory networks controlling fruit expansion in cucumber. Specific miRNAs could be selected for further functional research and molecular breeding in cucumber.

RNA-seq analysis and fluorescence imaging of melon powdery mildew disease reveal an orchestrated reprogramming of host physiology

The cucurbit powdery mildew elicited by Podosphaera xanthii is one of the most important limiting factors in cucurbit production. Our knowledge of the genetic and molecular bases underlying the physiological processes governing this disease is very limited. We used RNA-sequencing to identify differentially expressed genes in leaves of Cucumis melo upon inoculation with P. xanthii, using RNA samples obtained at different time points during the early stages of infection and their corresponding uninfected controls. In parallel, melon plants were phenotypically characterized using imaging techniques. We found a high number of differentially expressed genes (DEGs) in infected plants, which allowed for the identification of many plant processes that were dysregulated by the infection. Among those, genes involved in photosynthesis and related processes were found to be upregulated, whereas genes involved in secondary metabolism pathways, such as phenylpropanoid biosynthesis, were downregulated. These changes in gene expression could be functionally validated by chlorophyll fluorescence imaging and blue-green fluorescence imaging analyses, which corroborated the alterations in photosynthetic activity and the suppression of phenolic compound biosynthesis. The powdery mildew disease in melon is a consequence of a complex and multifaceted process that involves the dysregulation of many plant pathways such as primary and secondary metabolism.

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

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

A CsTu-TS1 regulatory module promotes fruit tubercule formation in cucumber

The fruit epidermal features such as the size of tubercules are important fruit quality traits for cucumber production. But the mechanisms underlying tubercule formation remain elusive. Here, tubercule size locus CsTS1 was identified by map-based cloning and was found to encode an oleosin protein. Allelic variation was identified in the promoter region of CsTS1, resulting in low expression of CsTS1 in all 22 different small-warty or nonwarty cucumber lines. High CsTS1 expression levels were closely correlated with increased fruit tubercule size among 44 different cucumber lines. Transgenic complementation and RNAi-mediated gene silencing of CsTS1 in transgenic cucumber plants demonstrated that CsTS1 positively regulates the development of tubercules. CsTS1 is highly expressed in the peel at fruit tubercule forming and enlargement stage. Auxin content and expression of three auxin signalling pathway genes were altered in the 35S:CsTS1 and CsTS1-RNAi fruit tubercules, a result that was supported by comparing the cell size of the control and transgenic fruit tubercules. CsTu, a C₂ H₂ zinc finger domain transcription factor that regulates tubercule initiation, binds directly to the CsTS1 promoter and promotes its expression. Taken together, our results reveal a novel mechanism in which the CsTu-TS1 complex promotes fruit tubercule formation in cucumber.

Cucumber ovaries inhibited by dominant fruit express a dynamic developmental program, distinct from either senescence-determined or fruit-setting ovaries

Cucurbits represent an attractive model to explore the dynamics of fruit set, whose regulation is not fully understood, despite its importance for yield determination. A fertilized ovary must integrate signals from distant plant parts and 'decide' whether to set fruit, or remain inhibited and later senesce. Here, we set out to characterize first-fruit inhibition (FFI), that is, the inhibitory effect of the first fruit on subsequent development of younger ovaries during pollination-induced and parthenocarpic fruit set. After the first fertilized ovaries set fruit, younger fertilized ovaries remained in a temporary state of inhibition. Such ovaries preserved their size and green color, and if the older fruit were removed within a 1-week reversibility window, they set fruit. The FFI effect was documented in both fertilized and parthenocarpic ovaries. We compared the gene expression profiles of pollinated ovaries (committed to set fruit) with respect to those affected by FFI, and to non-pollinated ovaries (undergoing senescence). The three fates of the ovaries were characterized by wide changes in gene expression, with several specific transcripts being up- or down-regulated in response to pollination, and to the presence of inhibitory fruit. Metabolic profiling was undertaken and integrated with the transcriptomic data in order to characterize early physiological changes that occur in post-anthesis ovaries in parthenocarpic and non-parthenocarpic genotypes. The combined results are discussed with respect to current models of fruit set and specifically with regard to FFI. Moreover, these metabolome and transcriptome data provide a valuable resource for studying ovary development and fruit set.

Transcriptome profiling of pumpkin (Cucurbita moschata Duch.) leaves infected with powdery mildew

Cucurbit powdery mildew (PM) is one of the most severe fungal diseases, but the molecular mechanisms underlying PM resistance remain largely unknown, especially in pumpkin (Cucurbita moschata Duch.). The goal of this study was to identify gene expression differences in PM-treated plants (harvested at 24 h and 48 h after inoculation) and untreated (control) plants of inbred line “112–2” using RNA sequencing (RNA-Seq). The inbred line “112–2” has been purified over 8 consecutive generations of self-pollination and shows high resistance to PM. More than 7600 transcripts were examined in pumpkin leaves, and 3129 and 3080 differentially expressed genes (DEGs) were identified in inbred line “112–2” at 24 and 48 hours post inoculation (hpi), respectively. Based on the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway database and GO (Gene Ontology) database, a complex regulatory network for PM resistance that may involve hormone signal transduction pathways, transcription factors and defense responses was revealed at the transcription level. In addition, the expression profiles of 16 selected genes were analyzed using quantitative RT-PCR. Among these genes, the transcript levels of 6 DEGs, including bHLH87 (Basic Helix-loop-helix transcription factor), ERF014 (Ethylene response factor), WRKY21 (WRKY domain), HSF (heat stress transcription factor A), MLO3 (Mildew Locus O), and SGT1 (Suppressor of G-Two Allele of Skp1), in PM-resistant “112–2” were found to be significantly up- or down-regulated both before 9 hpi and at 24 hpi or 48 hpi; this behavior differed from that observed in the PM-susceptible material (cultivar “Jiujiangjiaoding”). The transcriptome data provide novel insights into the response of Cucurbita moschata to PM stress and are expected to be highly useful for dissecting PM defense mechanisms in this major vegetable and for improving pumpkin breeding with enhanced resistance to PM.

Identification of Conserved and Diverse Metabolic Shift of the Stylar, Intermediate and Peduncular Segments of Cucumber Fruit during Development

Cucumber (Cucumis sativus L.) is one of the most important vegetables and contains a high content of nutritionally beneficial metabolites. However, little is known about the metabolic variations among different parts of cucumber fruit and their kinetics during growth. In this study, the dynamic metabolic profiles in the stylar end, the intermediate segment and the peduncular end of cucumber fruit during the development were investigated by employing a non-targeted metabolomics approach, where 238 metabolites were identified. Statistical analyses revealed that both development time and tissue type influenced metabolic changes, while development time seemed to exert more effects than tissue type on the cucumber fruit metabolome. The levels of the most of the detected metabolites decreased gradually, while those of some amino acids, carbohydrates and flavonoids increased across development. The metabolomes of the stylar end and the intermediate segment were similar, although all three parts of the cucumber fruit were separated from each other in orthogonal partial least squares projection to latent structures-discriminant analysis (OPLS-DA) plots. Metabolites association analysis revealed that sn-1 and sn-2 lysophospholipids are synthesized via independent pathways in cucumber fruit. In sum, this study demonstrated both conserved and diverse metabolic kinetics of three parts of cucumber fruit, which will facilitate further study of the regulation of cucumber fruit development as well as their potential applications in nutritious quality improvement of cucumber fruit.

Silicon bioavailability in exocarp of Cucumis sativus Linn

Scanning electron microscopy (SEM) and electron probe micro-analyzer (EPMA) techniques have been used to detect the silicon bioavailability in the exocarp of warty cucumber surface. Warts appear at the time of anthesis and are remnant part of spines/trichomes which on further fruit maturation abscised from the exocarp. Results of EPMA and phytolith analysis clearly revealed that the surface of exocarp (fruit) of Cucumis sativus Linn. containing warts has greater quantity of silicon as compared to the other part of the fruit. Besides silicon, some other elements were also found, on the fruit exocarp and its surrounding area. The other elements are magnesium (Mg), aluminum (Al), phosphorus (P), sulfur (S), potassium (K), calcium (Ca), iron (Fe), nickel (Ni), copper (Cu), and sodium. The percentage of silica is highest followed by Ni, Ca, Al, P, Mg, Fe, S, Cu, K, and Cl. Thus, this study clearly demonstrates that Cucumis sativus Linn. fruits which are used as salads and appetizers on daily basis are loaded with silicon and other useful elements and possess numerous health benefits.

Proteomic insight into fruit set of cucumber (Cucumis sativus L.) suggests the cues of hormone-independent parthenocarpy

BACKGROUND

Parthenocarpy is an excellent agronomic trait that enables crops to set fruit in the absence of pollination and fertilization, and therefore to produce seedless fruit. Although parthenocarpy is widely recognized as a hormone-dependent process, hormone-insensitive parthenocarpy can also be observed in cucumber; however, its mechanism is poorly understood. To improve the global understanding of parthenocarpy and address the hormone-insensitive parthenocarpy shown in cucumber, we conducted a physiological and proteomic analysis of differently developed fruits.

RESULTS

Physiological analysis indicated that the natural hormone-insensitive parthenocarpy of 'EC1' has broad hormone-inhibitor resistance, and the endogenous hormones in the natural parthenocarpy (NP) fruits were stable and relatively lower than those of the non-parthenocarpic cultivar '8419 s-1.' Based on the iTRAQ technique, 683 fruit developmental proteins were identified from NP, cytokinin-induced parthenocarpic (CP), pollinated and unpollinated fruits. Gene Ontology (GO) analysis showed that proteins detected from both set and aborted fruits were involved in similar biological processes, such as cell growth, the cell cycle, cell death and communication. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that 'protein synthesis' was the major biological process that differed between fruit set and fruit abortion. Clustering analysis revealed that different protein expression patterns were involved in CP and NP fruits. Forty-one parthenocarpy-specialized DEPs (differentially expressed proteins) were screened and divided into two distinctive groups: NP-specialized proteins and CP-specialized proteins. Furthermore, qRT-PCR and western blot analysis indicated that NP-specialized proteins showed hormone- or hormone-inhibitor insensitive expression patterns in both ovaries and seedlings.

CONCLUSIONS

In this study, the global molecular regulation of fruit development in cucumber was revealed at the protein level. Physiological and proteomic comparisons indicated the presence of hormone-independent parthenocarpy and suppression of fruit abortion in cucumber. The proteomic analysis suggested that hormone-independent parthenocarpy is regulated by hormone-insensitive proteins such as the NP-specialized proteins. Moreover, the regulation of fruit abortion suppression may be closely related to protein synthesis pathways.

A high-density linkage map and QTL mapping of fruit-related traits in pumpkin (Cucurbita moschata Duch.)

Pumpkin (Cucurbita moschata) is an economically worldwide crop. Few quantitative trait loci (QTLs) were reported previously due to the lack of genomic and genetic resources. In this study, a high-density linkage map of C. moschata was structured by double-digest restriction site-associated DNA sequencing, using 200 F2 individuals of CMO-1 × CMO-97. By filtering 74,899 SNPs, a total of 3,470 high quality SNP markers were assigned to the map spanning a total genetic distance of 3087.03 cM on 20 linkage groups (LGs) with an average genetic distance of 0.89 cM. Based on this map, both pericarp color and strip were fined mapped to a novel single locus on LG8 in the same region of 0.31 cM with phenotypic variance explained (PVE) of 93.6% and 90.2%, respectively. QTL analysis was also performed on carotenoids, sugars, tuberculate fruit, fruit diameter, thickness and chamber width with a total of 12 traits. 29 QTLs distributed in 9 LGs were detected with PVE from 9.6% to 28.6%. It was the first high-density linkage SNP map for C. moschata which was proved to be a valuable tool for gene or QTL mapping. This information will serve as significant basis for map-based gene cloning, draft genome assembling and molecular breeding.

Variation in cucumber (Cucumis sativus L.) fruit size and shape results from multiple components acting pre-anthesis and post-pollination

Morphological, QTL, and gene expression analyses indicate variation in cucumber fruit size and shape results from orientation, timing, and extent of cell division and expansion, and suggest candidate gene factors. Variation in cucumber (Cucumis sativus L.) fruit size and shape is highly quantitative, implicating interplay of multiple components. Recent studies have identified numerous fruit size and shape quantitative trait loci (QTL); however, underlying factors remain to be determined. We examined ovary and fruit development of two sequenced cucumber genotypes with extreme differences in fruit size and shape, Chinese Long '9930' (CL9930), and pickling type 'Gy14'. Differences were observed in several independent factors that can influence size and shape: ovule number, rate and period of cell division in longitudinal and cross section in ovaries and fruit, timing and rate of fruit expansion in length and diameter, and cell shape. Level and timing of expression of select fruit growth stage marker genes and candidate fruit size gene homologs associated with cucumber fruit size and shape QTL were examined from 5-day pre-anthesis to 20-day post-pollination. Our results indicate that variation in fruit size and shape results from differences in cell number and shape in longitudinal and cross section, driven in turn by differences in orientation, timing, and duration of cell division and expansion, both pre- and post-anthesis, and suggest candidate genes contributing to determination of cucumber fruit size and shape.

Transcriptomic analysis of short-fruit 1 (sf1) reveals new insights into the variation of fruit-related traits in Cucumis sativus

Fruit size is an important quality trait in different market classes of Cucumis sativus L., an economically important vegetable cultivated worldwide, but the genetic and molecular mechanisms that control fruit size are largely unknown. In this study, we isolated a natural cucumber mutant, short fruit 1 (sf1), caused by a single recessive Mendelian factor, from the North China-type inbred line CNS2. In addition to significantly decreased fruit length, other fruit-related phenotypic variations were also observed in sf1 compared to the wild-type (WT) phenotype, indicating that sf1 might have pleiotropic effects. Microscopic imaging showed that fruit cell size in sf1 was much larger than that in WT, suggesting that the short fruit phenotype in sf1 is caused by decreased cell number. Fine mapping revealed that sf1 was localized to a 174.3 kb region on chromosome 6. Similarly, SNP association analysis of bulked segregant RNA-Seq data showed increased SNP frequency in the same region of chromosome 6. In addition, transcriptomic analysis revealed that sf1 might control fruit length through the fine-tuning of cytokinin and auxin signalling, gibberellin biosynthesis and signal transduction in cucumber fruits. Overall, our results provide important information for further study of fruit length and other fruit-related features in cucumber.

Transcriptomic and metabolomic analyses of cucumber fruit peels reveal a developmental increase in terpenoid glycosides associated with age-related resistance to Phytophthora capsici

The oomycete, Phytophthora capsici, infects cucumber (Cucumis sativus L.) fruit. An age-related resistance (ARR) to this pathogen was previously observed in fruit of cultivar 'Vlaspik' and shown to be associated with the peel. Young fruits are highly susceptible, but develop resistance at ~10-12 days post pollination (dpp). Peels from resistant (16 dpp) versus susceptible (8 dpp) age fruit are enriched with genes associated with defense, and methanolic extracts from resistant age peels inhibit pathogen growth. Here we compared developing fruits from 'Vlaspik' with those of 'Gy14', a line that does not exhibit ARR. Transcriptomic analysis of peels of the two lines at 8 and 16 dpp identified 80 genes that were developmentally upregulated in resistant 'Vlaspik' 16 dpp versus 8 dpp, but not in susceptible 'Gy14' at 16 dpp. A large number of these genes are annotated to be associated with defense and/or specialized metabolism, including four putative resistance (R) genes, and numerous genes involved in flavonoid and terpenoid synthesis and decoration. Untargeted metabolomic analysis was performed on extracts from 8 and 16 dpp 'Vlaspik' and 'Gy14' fruit peels using Ultra-Performance Liquid Chromatography and Quadrupole Time-of-Flight Mass Spectrometry. Multivariate analysis of the metabolomes identified 113 ions uniquely abundant in resistant 'Vlaspik' 16 dpp peel extracts. The most abundant compounds in this group had relative mass defects consistent with terpenoid glycosides. Two of the three most abundant ions were annotated as glycosylated nor-terpenoid esters. Together, these analyses reveal potential mechanisms by which ARR to P. capsici may be conferred.

Synthesis of Pyrimethanil-Loaded Mesoporous Silica Nanoparticles and Its Distribution and Dissipation in Cucumber Plants

Mesoporous silica nanoparticles are used as pesticide carries in plants, which has been considered as a novel method to reduce the indiscriminate use of conventional pesticides. In the present work, mesoporous silica nanoparticles with particle diameters of 200-300 nm were synthesized in order to obtain pyrimethanil-loaded nanoparticles. The microstructure of the nanoparticles was observed by scanning electron microscopy. The loading content of pyrimethanil-loaded nanoparticles was investigated. After treatment on cucumber leaves, the concentrations of pyrimethanil were determined in different parts of cucumber over a period of 48 days using high performance liquid chromatography tandem mass spectrometry. It was shown that the pyrimethanil-loaded mesoporous silica nanoparticles might be more conducive to acropetal, rather than basipetal, uptake, and the dosage had almost no effect on the distribution and dissipation rate in cucumber plants. The application of the pesticide-loaded nanoparticles in leaves had a low risk of pyrimethanil accumulating in the edible part of the plant.

De novo comparative transcriptome analysis of genes involved in fruit morphology of pumpkin cultivars with extreme size difference and development of EST-SSR markers

The genetic basis of fruit size and shape was investigated for the first time in Cucurbita species and genetic loci associated with fruit morphology have been identified. Although extensive genomic resources are available at present for tomato (Solanum lycopersicum), cucumber (Cucumis sativus), melon (Cucumis melo) and watermelon (Citrullus lanatus), genomic databases for Cucurbita species are limited. Recently, our group reported the generation of pumpkin (Cucurbita pepo) transcriptome databases from two contrasting cultivars with extreme fruit sizes. In the current study we used these databases to perform comparative transcriptome analysis in order to identify genes with potential roles in fruit morphology and fruit size. Differential Gene Expression (DGE) analysis between cv. 'Munchkin' (small-fruit) and cv. 'Big Moose' (large-fruit) revealed a variety of candidate genes associated with fruit morphology with significant differences in gene expression between the two cultivars. In addition, we have set the framework for generating EST-SSR markers, which discriminate different C. pepo cultivars and show transferability to related Cucurbitaceae species. The results of the present study will contribute to both further understanding the molecular mechanisms regulating fruit morphology and furthermore identifying the factors that determine fruit size. Moreover, they may lead to the development of molecular marker tools for selecting genotypes with desired morphological traits.

Genome-wide characterization and expression profiling of the NAC genes under abiotic stresses in Cucumis sativus

The NAC (standing for no apical meristem [NAM], Arabidopsis transcription activation factor [ATAF] and cup-shaped cotyledon [CUC]) proteins pertain to one of the plant-specific transcription factor families that play important roles in plant development, abiotic stress resistance and signalling transduction. In the present study, the genomic features of the NAC genes in cucumber were analysed in depth using in silico tools. To reveal a tissue-specific, abiotic stress and hormone-responsive expression profile of CsNAC genes, RT-qPCR was performed under different treatments. Phylogenetic analyses and genome-wide annotation indicated that 82 high-confidence CsNAC genes were clustered into 13 sub-groups with uneven distribution in the cucumber genome. Furthermore, the CsNAC genes exhibited different tissue-specific expression patterns in 10 tissues under normal growth conditions, while 13 (16%) and 28 (34%) genes displayed preferential expression in roots and flowers, respectively. Moreover, CsNAC genes were more sensitive to salinity than other stresses; however, their responses were relatively rapid and transient to nutrition deprivation. Several CsNAC genes, including CsNAC35, which is an orthologue of the known stress-responsive Arabidopsis RD26, were identified as highly responsive to abiotic stresses and hormones. Overall, our findings revealed the genomic landscape and expression profiling of the CsNAC genes in response to multiple stresses and hormones, offering clues for further function analyses and molecular breeding.

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.

Effect of nitrogen deficiency on ascorbic acid biosynthesis and recycling pathway in cucumber seedlings

L-Ascorbic acid (AsA, ascorbate) is one of the most abundant natural antioxidants, and it is an important factor in the nutritional quality of cucumber. In this work, key enzymes involved in the ascorbic acid biosynthesis and recycling pathway in cucumber seedlings under nitrogen deficiency were investigated at the levels of transcription and enzyme activity. The activities of myo-inositol oxygenase (MIOX) and transcript levels of MIOXs increased dramatically, while the activities of ascorbate oxidase (AO) and glutathione reductase (GR) and transcript levels of AOs and GR2 decreased significantly in N-limited leaves, as did the ascorbate concentration, in nitrogen-deficient cucumber seedlings. The activities of other enzymes and transcript levels of other genes involved in the ascorbate recycling pathway and ascorbate synthesis pathways decreased or remained unchanged under nitrogen deficiency. These results indicate that nitrogen deficiency induced genes involved in the ascorbate-glutathione recycling and myo-inositol pathway in cucumber leaves. Thus, the AO, GR and MIOX involved in the pathways might play roles in AsA accumulation.

Rapid identification of fruit length loci in cucumber (Cucumis sativus L.) using next-generation sequencing (NGS)-based QTL analysis

The cucumber (Cucumis sativus L.) exhibits extensive variations in fruit size and shape. Fruit length is an important agronomic and domesticated trait controlled by quantitative trait loci (QTLs). Nonetheless, the underlying molecular and genetic mechanisms that determine cucumber fruit length remain unclear. QTL-seq is an efficient strategy for QTL identification that takes advantage of bulked-segregant analysis (BSA) and next-generation sequencing (NGS). In the present study, we conducted QTL mapping and QTL-seq of cucumber fruit length. QTL mapping identified 8 QTLs for immature and mature fruit length. A major-effect QTL fl3.2, which explained a maximum of 38.87% of the phenotypic variation, was detected. A genome-wide comparison of SNP profiles between two DNA bulks identified 6 QTLs for ovary length. QTLs ovl3.1 and ovl3.2 both had major effects on ovary length with a △ (SNP-index) of 0.80 (P < 0.01) and 0.74 (P < 0.01), respectively. Quantitative RT-PCR of fruit size-related homologous genes localized in the consensus QTL FL3.2 was conducted. Four candidate genes exhibited increased expression levels in long fruit genotypes. Our results demonstrated the power of the QTL-seq method in rapid QTL detection and provided reliable QTL regions for fine mapping of fruit length-related loci and for identifying candidate genes.

The WD-Repeat Protein CsTTG1 Regulates Fruit Wart Formation through Interaction with the Homeodomain-Leucine Zipper I Protein Mict

The cucumber (Cucumis sativus) fruit is covered with bloom trichomes and warts (composed of spines and tubercules), which have an important impact on the commercial value of the crop. However, little is known about the regulatory mechanism underlying their formation. Here, we reported that the cucumber WD-repeat homolog CsTTG1, which is localized in the nucleus and cytomembrane, plays an important role in the formation of cucumber fruit bloom trichomes and warts. Functional characterization of CsTTG1 revealed that it is mainly expressed in the epidermis of cucumber ovary and that its overexpression in cucumber alters the density of fruit bloom trichomes and spines, thereby promoting the warty fruit trait. Conversely, silencing CsTTG1 expression inhibits the initiation of fruit spines. Molecular and genetic analyses showed that CsTTG1 acts in parallel to Mict/CsGL1, a key trichome formation factor, to regulate the initiation of fruit trichomes, including fruit bloom trichomes and spines, and that the further differentiation of fruit spines and formation of tubercules regulated by CsTTG1 is dependent on Mict Using yeast two-hybrid assay and bimolecular fluorescence complementation assay, we determined that CsTTG1 directly interacts with Mict. Collectively, our results indicate that CsTTG1 is an important component of the molecular network that regulates fruit bloom trichome and wart formation in cucumber.

Next generation sequencing and omics in cucumber (Cucumis sativus L.) breeding directed research

In the post-genomic era the availability of genomic tools and resources is leading us to novel generation methods in plant breeding, as they facilitate the study of the genotype and its relationship with the phenotype, in particular for complex traits. In this study we have mainly concentrated on the Cucumis sativus and (but much less) Cucurbitaceae family several important vegetable crops. There are many reports on research conducted in Cucurbitaceae plant breeding programs on the ripening process, phloem transport, disease resistance, cold tolerance and fruit quality traits. This paper presents the role played by new omic technologies in the creation of knowledge on the mechanisms of the formation of the breeding features. The analysis of NGS (NGS-next generation sequencing) data allows the discovery of new genes and regulatory sequences, their positions, and makes available large collections of molecular markers. Genome-wide expression studies provide breeders with an understanding of the molecular basis of complex traits. Firstly a high density map should be created for the reference genome, then each re-sequencing data could be mapped and new markers brought out into breeding populations. The paper also presents methods that could be used in the future for the creation of variability and genomic modification of the species in question. It has been shown also the state and usefulness in breeding the chloroplastomic and mitochondriomic study.

Phloem transcriptome signatures underpin the physiological differentiation of the pedicel, stalk and fruit of cucumber (Cucumis sativus L.)

Cucumber is one of the most important vegetables grown worldwide due to its important economic and nutritional value. The cucumber fruit consists morphologically of the undesirable stalk and the tasty fruit; however, physiological differentiation of these two parts and the underlying molecular basis remain largely unknown. Here we characterized the physiological differences among the pedicel, stalk and fruit, and compared the respective phloem transcriptomes using laser capture microdissection coupled with RNA sequencing (RNA-Seq). We found that the pedicel was characterized by minor cell expansion and a high concentration of stachyose, the stalk showed rapid cell expansion and high raffinose accumulation, and the fruit featured transition from cell division to cell expansion and high levels of monosaccharides. Analyses of transcriptome data indicated that cell wall- and calcium ion binding-related genes contributed to the cell expansion in the pedicel and stalk, whereas genes implicated in cell cycle and hormone actions regulated the transition from cell division to cell expansion in the fruit. Differential sugar distribution in these three phloem-connected tissues resulted from tissue-specific sugar metabolism and transport. Enrichment of transcription factors in the stalk and fruit may facilitate nutrient accumulation in these sink organs. As such, phloem-located gene expression partially orchestrated physiological differentiation of the pedicel, stalk and fruit in cucumber. In addition, we identified 432 cucumber-unique genes and five phloem markers guiding future functional studies.

A high-density genetic map for anchoring genome sequences and identifying QTLs associated with dwarf vine in pumpkin (Cucurbita maxima Duch.)

Background

Pumpkin (Cucurbita maxima Duch.) is an economically important crop belonging to the Cucurbitaceae family. However, very few genomic and genetic resources are available for this species. As part of our ongoing efforts to sequence the pumpkin genome, high-density genetic map is essential for anchoring and orienting the assembled scaffolds. In addition, a saturated genetic map can facilitate quantitative trait locus (QTL) mapping.

Results

A set of 186 F₂ plants derived from the cross of pumpkin inbred lines Rimu and SQ026 were genotyped using the genotyping-by-sequencing approach. Using the SNPs we identified, a high-density genetic map containing 458 bin-markers was constructed, spanning a total genetic distance of 2,566.8 cM across the 20 linkage groups of C. maxima with a mean marker density of 5.60 cM. Using this map we were able to anchor 58 assembled scaffolds that covered about 194.5 Mb (71.7 %) of the 271.4 Mb assembled pumpkin genome, of which 44 (183.0 Mb; 67.4 %) were oriented. Furthermore, the high-density genetic map was used to identify genomic regions highly associated with an important agronomic trait, dwarf vine. Three QTLs on linkage groups (LGs) 1, 3 and 4, respectively, were recovered. One QTL, qCmB2, which was located in an interval of 0.42 Mb on LG 3, explained 21.4 % phenotypic variations. Within qCmB2, one gene, Cma_004516, encoding the gibberellin (GA) 20-oxidase in the GA biosynthesis pathway, had a 1249-bp deletion in its promoter in bush type lines, and its expression level was significantly increased during the vine growth and higher in vine type lines than bush type lines, supporting Cma_004516 as a possible candidate gene controlling vine growth in pumpkin.

Conclusions

A high-density pumpkin genetic map was constructed, which was used to successfully anchor and orient the assembled genome scaffolds, and to identify QTLs highly associated with pumpkin vine length. The map provided a valuable resource for gene cloning and marker assisted breeding in pumpkin and other related species. The identified vine length QTLs would help to dissect the underlying molecular basis regulating pumpkin vine growth.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2312-8) contains supplementary material, which is available to authorized users.

Exocarp Properties and Transcriptomic Analysis of Cucumber (Cucumis sativus) Fruit Expressing Age-Related Resistance to Phytophthora capsici

Very young cucumber (Cucumis sativus) fruit are highly susceptible to infection by the oomycete pathogen, Phytophthora capsici. As the fruit complete exponential growth, at approximately 10-12 days post pollination (dpp), they transition to resistance. The development of age-related resistance (ARR) is increasingly recognized as an important defense against pathogens, however, underlying mechanisms are largely unknown. Peel sections from cucumber fruit harvested at 8 dpp (susceptible) and 16 dpp (resistant) showed equivalent responses to inoculation as did whole fruit, indicating that the fruit surface plays an important role in defense against P. capsici. Exocarp from 16 dpp fruit had thicker cuticles, and methanolic extracts of peel tissue inhibited growth of P. capsici in vitro, suggesting physical or chemical components to the ARR. Transcripts specifically expressed in the peel vs. pericarp showed functional differentiation. Transcripts predominantly expressed in the peel were consistent with fruit surface associated functions including photosynthesis, cuticle production, response to the environment, and defense. Peel-specific transcripts that exhibited increased expression in 16 dpp fruit relative to 8 dpp fruit, were highly enriched (P<0.0001) for response to stress, signal transduction, and extracellular and transport functions. Specific transcripts included genes associated with potential physical barriers (i.e., cuticle), chemical defenses (flavonoid biosynthesis), oxidative stress, penetration defense, and molecular pattern (MAMP)-triggered or effector-triggered (R-gene mediated) pathways. The developmentally regulated changes in gene expression between peels from susceptible- and resistant- age fruits suggest programming for increased defense as the organ reaches full size.