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

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.

Expression patterns of genes encoding plasma membrane aquaporins during fruit development in cucumber (Cucumis sativus L.)

Aquaporins are membrane channels precisely regulating water movement through cell membranes in most living organisms. Despite the advances in the physiology of fruit development, their participation during fruit development in cucumber still barely understood. In this paper, the expressions of 12 genes encoding plasma membrane intrinsic proteins (PIPs) were analyzed during cucumber fruit development in our work. Based on the homology search with known PIPs from rice, Arabidopsis and strawberry, 12 cucumber PIP genes subfamily members were identified. Cellular localization assays indicated that CsPIPs were localized in the plasma membrane. The qRT-PCR analysis of CsPIPs showed that 12 CsPIPs were differentially expressed during fruit development. These results suggest that 12 genes encoding plasma membrane intrinsic proteins (CsPIPs) play very important roles in cucumber life cycle and the data generated will be helpful in understanding their precise roles during fruit development in cucumber.

QTL mapping of cucumber fruit flesh thickness by SLAF-seq

Cucumber is an agriculturally and economically important vegetable crop worldwide. Fruit flesh thickness is an important trait for cucumber and also a central determinant of yield, yet little is known about the underlying mechanism of this trait. In this study, bulked segregant analysis (BSA) combined with specific length amplified fragment sequencing (SLAF-seq) was applied to finely map the gene that underlies fruit flesh thickness in cucumber. A 0.19-Mb-long quantitative trait locus on chromosome 2 controlling fruit flesh thickness (QTL fft2.1) was identified and further confirmed by simple sequence repeat (SSR) marker-based classical QTL mapping in 138 F₂ individuals. Gene prediction of this 0.19-Mb region identified 20 genes. Quantitative RT-PCR revealed higher expression levels of Csa2 M058670.1 (SET domain protein-lysine methyltransferase) in D8 (thick fruit flesh parent) compared with that in XUE1 (thin fruit flesh parent) during fruit development. Sequence alignment analysis of Csa2M058670.1 from thick and thin fruit flesh cucumber lines revealed a 4-bp deletion mutation in the promoter region of this candidate gene, which may result in the loss of Csa2M058670.1 activation in thin fruit flesh lines. The data presented herein suggest that Csa2M058670.1 is a possible candidate gene for controlling flesh thickness in cucumber.

QTL mapping in multiple populations and development stages reveals dynamic quantitative trait loci for fruit size in cucumbers of different market classes

QTL analysis in multi-development stages with different QTL models identified 12 consensus QTLs underlying fruit elongation and radial growth presenting a dynamic view of genetic control of cucumber fruit development. Fruit size is an important quality trait in cucumber (Cucumis sativus L.) of different market classes. However, the genetic and molecular basis of fruit size variations in cucumber is not well understood. In this study, we conducted QTL mapping of fruit size in cucumber using F2, F2-derived F3 families and recombinant inbred lines (RILs) from a cross between two inbred lines Gy14 (North American picking cucumber) and 9930 (North China fresh market cucumber). Phenotypic data of fruit length and diameter were collected at three development stages (anthesis, immature and mature fruits) in six environments over 4 years. QTL analysis was performed with three QTL models including composite interval mapping (CIM), Bayesian interval mapping (BIM), and multiple QTL mapping (MQM). Twenty-nine consistent and distinct QTLs were detected for nine traits from multiple mapping populations and QTL models. Synthesis of information from available fruit size QTLs allowed establishment of 12 consensus QTLs underlying fruit elongation and radial growth, which presented a dynamic view of genetic control of cucumber fruit development. Results from this study highlighted the benefits of QTL analysis with multiple QTL models and different mapping populations in improving the power of QTL detection. Discussion was presented in the context of domestication and diversifying selection of fruit length and diameter, marker-assisted selection of fruit size, as well as identification of candidate genes for fruit size QTLs in cucumber.

RNA-Seq-based transcriptome profiling of early nitrogen deficiency response in cucumber seedlings provides new insight into the putative nitrogen regulatory network

Nitrogen (N) is both an important macronutrient and a signal for plant growth and development. However, the early regulatory mechanism of plants in response to N starvation is not well understood, especially in cucumber, an economically important crop that normally consumes excessive N during production. In this study, the early time-course transcriptome response of cucumber leaves under N deficiency was monitored using RNA sequencing (RNA-Seq). More than 23,000 transcripts were examined in cucumber leaves, of which 364 genes were differentially expressed in response to N deficiency. Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database, gene ontology (GO) and protein-protein interaction analysis, 64 signaling-related N-deficiency-responsive genes were identified. Furthermore, the potential regulatory mechanisms of anthocyanin accumulation, Chl decline and cell wall remodeling were assessed at the transcription level. Increased ascorbic acid synthesis was identified in cucumber seedlings and fruit under N-deficient conditions, and a new corresponding regulatory hypothesis has been proposed. A data cross-comparison between model plants and cucumber was made, and some common and specific N-deficient response mechanisms were found in the present study. Our study provides novel insights into the responses of cucumber to nitrogen starvation at the global transcriptome level, which are expected to be highly useful for dissecting the N response pathways in this major vegetable and for improving N fertilization practices.

Sequencing, assembly, annotation, and gene expression: novel insights into the hormonal control of carrot root development revealed by a high-throughput transcriptome

Previous studies have indicated that hormonal control is essential for plant root growth. The root of the carrot is an edible vegetable with a high nutritional value. However, molecular mechanisms underlying hormone-mediated root growth of carrot have not been illustrated. Therefore, the present study collected carrot root samples from four developmental stages, and performed transcriptome sequencing to understand the molecular functions of plant hormones in carrot root growth. A total of 160,227 transcripts were generated from our transcriptome, which were assembled into 32,716 unigenes with an average length of 1,453 bp. A total of 4,818 unigenes were found to be differentially expressed between the four developmental stages. In total, 87 hormone-related differentially expressed genes were identified, and the roles of the hormones are extensively discussed. Our results suggest that plant hormones may regulate carrot root growth in a phase-dependent manner, and these findings will provide valuable resources for future research on carrot root development.

An acorn squash (Cucurbita pepo ssp. ovifera) fruit and seed transcriptome as a resource for the study of fruit traits in Cucurbita

Acorn squash (Cucurbita pepo) is an iconic fall vegetable in the United States, known for its unique fruit shape and also prized for its culinary properties. Little is known about the metabolism that underlies the development of fruit quality attributes such as color, sweetness, texture and nutritional qualities in acorn squash, or any other winter squash grown worldwide. To provide insight into winter squash fruit and seed development and add to the genomic resources in the Cucurbita genus, RNA sequencing was used to generate an acorn squash fruit and seed transcriptome from the cultivar Sweet REBA at critical points throughout fruit development. 141 838 600 high-quality paired-end Illumina reads were assembled into 55 949 unigenes. 85% of unigenes with predicted open reading frames had homology with previously identified genes and over 62% could be functionally annotated. Comparison with the watermelon and cucumber genomes provided confirmation that the unigenes are full-length and comprehensive, covering an average of 90% of the coding sequence of their homologs and 72% of the cucumber and watermelon exomes. Key candidate genes associated with carotenoid and carbohydrate metabolism were identified toward a resource for winter squash fruit quality trait dissection. This transcriptome represents a major advance in C. pepo genomics, providing significant new sequence information and revealing the repertoire of genes expressed throughout winter squash fruit and seed development. Future studies on the genetic basis of fruit quality and future breeding efforts will be enhanced by tools and insights developed from this resource.

Transcriptomic analysis reveals the roles of microtubule-related genes and transcription factors in fruit length regulation in cucumber (Cucumis sativus L.)

Cucumber (Cucumis sativus L.) fruit is a type of fleshy fruit that is harvested immaturely. Early fruit development directly determines the final fruit length and diameter, and consequently the fruit yield and quality. Different cucumber varieties display huge variations of fruit length, but how fruit length is determined at the molecular level remains poorly understood. To understand the genes and gene networks that regulate fruit length in cucumber, high throughout RNA-Seq data were used to compare the transcriptomes of early fruit from two near isogenic lines with different fruit lengths. 3955 genes were found to be differentially expressed, among which 2368 genes were significantly up-regulated and 1587 down-regulated in the line with long fruit. Microtubule and cell cycle related genes were dramatically activated in the long fruit, and transcription factors were implicated in the fruit length regulation in cucumber. Thus, our results built a foundation for dissecting the molecular mechanism of fruit length control in cucumber, a key agricultural trait of significant economic importance.

Molecular mapping reveals structural rearrangements and quantitative trait loci underlying traits with local adaptation in semi-wild Xishuangbanna cucumber (Cucumis sativus L. var. xishuangbannanesis Qi et Yuan)

Comparative genetic mapping revealed the origin of Xishuangbanna cucumber through diversification selection after domestication. QTL mapping provided insights into the genetic basis of traits under diversification selection during crop evolution. The Xishuangbanna cucumber, Cucumis sativus L. var. xishuangbannanesis Qi et Yuan (XIS), is a semi-wild landrace from the tropical southwest China with some unique traits that are very useful for cucumber breeding, such as tolerance to low light, large fruit size, heavy fruit weight, and orange flesh color in mature fruits. In this study, using 124 recombinant inbred lines (RILs) derived from the cross of the XIS cucumber with a cultivated cucumber inbred line, we developed a linkage map with 269 microsatellite (or simple sequence repeat) markers which covered 705.9 cM in seven linkage groups. Comparative analysis of orders of common marker loci or marker-anchored draft genome scaffolds among the wild (C. sativus var. hardwickii), semi-wild, and cultivated cucumber genetic maps revealed that the XIS cucumber shares major chromosomal rearrangements in chromosomes 4, 5, and 7 between the wild and cultivated cucumbers suggesting that the XIS cucumber originated through diversifying selection after cucumber domestication. Several XIS-specific minor structural changes were identified in chromosomes 1 and 6. QTL mapping with the 124 RILs in four environments identified 13 QTLs for domestication and diversifying selection-related traits including 2 for first female flowering time (fft1.1, fft6.1), 5 for mature fruit length (fl1.1, fl3.1, fl4.1, fl6.1, and fl7.1), 3 for fruit diameter (fd1.1, fd4.1, and fd6.1), and 3 for fruit weight (fw2.1, fw4.1, and fw6.1). Six of the 12 QTLs were consistently detected in all four environments. Among the 13 QTLs, fft1.1, fl1.1, fl3.1, fl7.1, fd4.1, and fw6.1 were major-effect QTLs for respective traits with each explaining at least 10 % of the observed phenotypic variations. Results from this study provide insights into the cytological and genetic basis of crop evolution leading to the XIS cucumber. The molecular markers associated with the QTLs should be useful in exploring the XIS cucumber genetic resources for cucumber breeding.

Transcriptome profiling reveals roles of meristem regulators and polarity genes during fruit trichome development in cucumber (Cucumis sativus L.)

Trichomes are epidermal hair-like structures that function in plant defence against biotic and abiotic stresses. Extensive studies have been performed on foliar trichomes development in Arabidopsis and tomato, but the molecular mechanism of fruit trichome formation remains elusive. Cucumber fruit is covered with trichomes (spines) that directly affect the appearance and quality of cucumber products. Here, we characterized the fruit spine development in wild-type (WT) cucumber and a spontaneous mutant, tiny branched hair (tbh). Our data showed that the cucumber trichome was multicellular and non-glandular, with malformed organelles and no endoreduplication. Fruit spine development was generally homogenous and marked by a rapid base expansion stage. Trichomes in the tbh mutant were tiny and branched, with increased density and aberrant cell shape. Transcriptome profiling indicated that meristem-related genes were highly enriched in the upregulated genes in the tbh versus the WT, as well as in WT spines after versus before base expansion, and that polarity regulators were greatly induced during spine base expansion. Quantitative reverse transcription PCR and in situ hybridization confirmed the differential expression of CUP-SHAPED COTYLEDON3 (CUC3) and SHOOT MERISTEMLESS (STM) during spine development. Therefore, cucumber trichomes are morphologically different from those of Arabidopsis and tomato, and their development may be regulated by a distinct pathway involving meristem genes and polarity regulators.

Dynamic changes in the date palm fruit proteome during development and ripening

Date palm (Phoenix dactylifera) is an economically important fruit tree in the Middle East and North Africa and is characterized by large cultivar diversity, making it a good model for studies on fruit development and other important traits. Here in gel comparative proteomics combined with tandem mass spectrometry were used to study date fruit development and ripening. Total proteins were extracted using a phenol-based protocol. A total of 189 protein spots were differentially regulated (p≤0.05). The identified proteins were classified into 14 functional categories. The categories with the most proteins were 'disease and defense' (16.5%) and 'metabolism' (15.4%). Twenty-nine proteins have not previously been identified in other fleshy fruits and 64 showed contrasting expression patterns in other fruits. Abundance of most proteins with a role in abiotic stress responses increased during ripening with the exception of heat shock proteins. Proteins with a role in anthocyanin biosynthesis, glycolysis, tricarboxylic acid cycle and cell wall degradation were upregulated particularly from the onset of ripening and during ripening. In contrast, expression of pentose phosphate- and photosynthesis-related proteins decreased during fruit maturation. Although date palm is considered a climacteric species, the analysis revealed downregulation of two enzymes involved in ethylene biosynthesis, suggesting an ethylene-independent ripening of 'Barhi' fruits. In summary, this proteomics study provides insights into physiological processes during date fruit development and ripening at the systems level and offers a reference proteome for the study of regulatory mechanisms that can inform molecular and biotechnological approaches to further improvements of horticultural traits including fruit quality and yield.

Characterization and expression profiling of cucumber kinesin genes during early fruit development: revealing the roles of kinesins in exponential cell production and enlargement in cucumber fruit

Rapid cell division and expansion in early fruit development are important phases for cucumber fruit yield and quality. Kinesin proteins are microtubule-based motors responsible for modulating cell division and enlargement. In this work, the candidate kinesin genes involved in rapid cell division and expansion during cucumber fruit development were investigated. The morphological and cellular changes during early fruit development were compared in four cucumber genotypes with varied fruit size. The correlation between the expression profiles of cucumber kinesin genes and cellular changes in fruit was investigated. Finally, the biochemical characteristics and subcellular localizations of three candidate kinesins were studied. The results clarified the morphological and cellular changes during early cucumber fruit development. This study found that CsKF2-CsKF6 were positively correlated with rapid cell production; CsKF1 and CsKF7 showed a strongly positive correlation with rapid cell expansion. The results also indicated that CsKF1 localized to the plasma membrane of fast-expanding fruit cells, that CsKF2 might play a role in fruit chloroplast division, and that CsKF3 is involved in the function or formation of phragmoplasts in fruit telophase cells. The results strongly suggest that specific fruit-enriched kinesins are specialized in their functions in rapid cell division and expansion during cucumber fruit development.

De novo assembly and characterization of the transcriptome, and development of SSR markers in wax gourd (Benicasa hispida)

BACKGROUND

Wax gourd is a widely used vegetable of Cucuribtaceae, and also has important medicinal and health values. However, the genomic resources of wax gourd were scarcity, and only a few nucleotide sequences could be obtained in public databases.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we examined transcriptome in wax gourd. More than 44 million of high quality reads were generated from five different tissues of wax gourd using Illumina paired-end sequencing technology. Approximately 4 Gbp data were generated, and de novo assembled into 65,059 unigenes, with an N50 of 1,132 bp. Based on sequence similarity search with known protein database, 36,070 (55.4%) showed significant similarity to known proteins in Nr database, and 24,969 (38.4%) had BLAST hits in Swiss-Prot database. Among the annotated unigenes, 14,994 of wax gourd unigenes were assigned to GO term annotation, and 23,977 were found to have COG classifications. In addition, a total of 18,713 unigenes were assigned to 281 KEGG pathways. Furthermore, 6,242 microsatellites (simple sequence repeats) were detected as potential molecular markers in wax gourd. Two hundred primer pairs for SSRs were designed for validation of the amplification and polymorphism. The result showed that 170 of the 200 primer pairs were successfully amplified and 49 (28.8%) of them exhibited polymorphisms.

CONCLUSION/SIGNIFICANCE

Our study enriches the genomic resources of wax gourd and provides powerful information for future studies. The availability of this ample amount of information about the transcriptome and SSRs in wax gourd could serve as valuable basis for studies on the physiology, biochemistry, molecular genetics and molecular breeding of this important vegetable crop.