Comparative transcriptome analysis provides insights into molecular mechanisms for parthenocarpic fruit development in eggplant (Solanum melongena L.)

The hormone regulatory mechanism underlying parthenocarpic fruit formation in tomato

Parthenocarpic fruits, known for their superior taste and reliable yields in adverse conditions, develop without the need for fertilization or pollination. Exploring the physiological and molecular mechanisms behind parthenocarpic fruit development holds both theoretical and practical significance, making it a crucial area of study. This review examines how plant hormones and MADS-box transcription factors control parthenocarpic fruit formation. It delves into various aspects of plant hormones-including auxin, gibberellic acid, cytokinins, ethylene, and abscisic acid-ranging from external application to biosynthesis, metabolism, signaling pathways, and their interplay in influencing parthenocarpic fruit development. The review also explores the involvement of MADS family gene functions in these processes. Lastly, we highlight existing knowledge gaps and propose directions for future research on parthenocarpy.

A Compendium for Novel Marker-Based Breeding Strategies in Eggplant

The worldwide production of eggplant is estimated at about 58 Mt, with China, India and Egypt being the major producing countries. Breeding efforts in the species have mainly focused on increasing productivity, abiotic and biotic tolerance/resistance, shelf-life, the content of health-promoting metabolites in the fruit rather than decreasing the content of anti-nutritional compounds in the fruit. From the literature, we collected information on mapping quantitative trait loci (QTLs) affecting eggplant's traits following a biparental or multi-parent approach as well as genome-wide association (GWA) studies. The positions of QTLs were lifted according to the eggplant reference line (v4.1) and more than 700 QTLs were identified, here organized into 180 quantitative genomic regions (QGRs). Our findings thus provide a tool to: (i) determine the best donor genotypes for specific traits; (ii) narrow down QTL regions affecting a trait by combining information from different populations; (iii) pinpoint potential candidate genes.

Genome-wide identification and characterization of parthenocarpic fruit set-related gene homologs in cucumber (Cucumis sativus L.)

Cucumber (Cucumis sativus L.), a major horticultural crop, in the family Cucurbitaceae is grown and consumed globally. Parthenocarpy is an ideal trait for many fruit and vegetables which produces seedless fruit desired by consumers. The seedlessness occurs when fruit develops without fertilization which can be either natural or induced. So far, a limited number of genes regulating parthenocarpic fruit set have been reported in several fruit or vegetable crops, most of which are involved in hormone biosynthesis or signalling. Although parthenocarpic cucumber has been widely used in commercial production for a long time; its genetic basis is not well understood. In this study, we retrieved thirty five parthenocarpy fruit-set related genes (PRGs) from bibliomic data in various plants. Thirty-five PRG homologs were identified in the cucumber genome via homology-based search. An in silico analysis was performed on phylogenetic tree, exon-intron structure, cis-regulatory elements in the promoter region, and conserved domains of their deduced proteins, which provided insights into the genetic make-up of parthenocarpy-related genes in cucumber. Simple sequence repeat (SSR) sequences were mined in these PRGs, and 31 SSR markers were designed. SSR genotyping identified three SSRs in two polymorphic genes. Quantitative real-time PCR of selected genes was conducted in five cucumber lines with varying degrees of parthenocarpic fruit set capacities, which revealed possible association of their expression with parthenocarpy. The results revealed that homologs CsWD40 and CsPIN-4 could be considered potential genes for determination of parthenocarpy as these genes showed parental polymorphism and differential gene expression in case of parthenocarpic and non-parthenocarpic parents.

Parthenocarpy-related genes induced by naphthalene acetic acid in oil palm interspecific O × G [Elaeis oleifera (Kunth) Cortés × Elaeis guineensis Jacq.] hybrids

Parthenocarpy is the development without fertilization of seedless fruits. In the oil palm industry, the development of parthenocarpic fruits is considered an attractive option to increase palm oil production. Previous studies have shown the application of synthetic auxins in Elaeis guineensis, and interspecific O×G hybrids (Elaeis oleifera (Kunth) Cortés × E. guineensis Jacq.) induces parthenocarpy. The aim of this study was to identify the molecular mechanism through transcriptomics and biology system approach to responding to how the application of NAA induces parthenocarpic fruits in oil palm O×G hybrids. The transcriptome changes were studied in three phenological stages (PS) of the inflorescences: i) PS 603, pre-anthesis III, ii) PS 607, anthesis, and iii) PS 700, fertilized female flower. Each PS was treated with NAA, Pollen, and control (any application). The expression profile was studied at three separate times: five minutes (T0), 24 hours (T1), and 48 h post-treatment (T2). The RNA sequencing (RNA seq) approach was used with 27 oil palm O×G hybrids for a total of 81 raw samples. RNA-Seq showed around 445,920 genes. Numerous differentially expressed genes (DEGs) were involved in pollination, flowering, seed development, hormone biosynthesis, and signal transduction. The expression of the most relevant transcription factors (TF) families was variable and dependent on the stage and time post-treatment. In general, NAA treatment expressed differentially more genes than Pollen. Indeed, the gene co-expression network of Pollen was built with fewer nodes than the NAA treatment. The transcriptional profiles of Auxin-responsive protein and Gibberellin-regulated genes involved in parthenocarpy phenomena agreed with those previously reported in other species. The expression of 13 DEGs was validated by RT-qPCR analysis. This detailed knowledge about the molecular mechanisms involved in parthenocarpy could be used to facilitate the future development of genome editing techniques that enable the production of parthenocarpic O×G hybrid cultivars without growth regulator application.

Options for the generation of seedless cherry, the ultimate snacking product

This manuscript identifies cherry orthologues of genes implicated in the development of pericarpic fruit and pinpoints potential options and restrictions in the use of these targets for commercial exploitation of parthenocarpic cherry fruit. Cherry fruit contain a large stone and seed, making processing of the fruit laborious and consumption by the consumer challenging, inconvenient to eat 'on the move' and potentially dangerous for children. Availability of fruit lacking the stone and seed would be potentially transformative for the cherry industry, since such fruit would be easier to process and would increase consumer demand because of the potential reduction in costs. This review will explore the background of seedless fruit, in the context of the ambition to produce the first seedless cherry, carry out an in-depth analysis of the current literature around parthenocarpy in fruit, and discuss the available technology and potential for producing seedless cherry fruit as an 'ultimate snacking product' for the twenty-first century.

Programmed Cell Death May Be Involved in the Seedless Phenotype Formation of Oil Palm

Oil palm (Elaeis guineensis Jacq.) is a well-known vegetable oil-yielding crop. Seedlessness is one of the most prominent traits in oil palm due to its low processing costs and high oil content. Nevertheless, an extensive study on molecular mechanisms regulating seedless phenotype formation in oil palm is very limited so far. In this study, stigma, style, and ovary from seedless and seeded (Tenera and Pisifera) oil palm trees were used to investigate the possible mechanism. Results showed that non-pollination resulted in no fruits, and self- and cross-pollinations resulted in seedless fruits, while boron treatment had no effect on seedless phenotype formation, implying that seedless trees have incomplete self and outcrossing incompatibility. Furthermore, the transcriptome data analysis highlighted eight programmed cell death (PCD) genes and three groups of PCD-related genes: 4-coumarate-CoA ligase (4CL), S-RNase, and MADS-box. The majority of these genes were significantly up-regulated in the stigma and style of Seedless palm trees compared to Tenera and Pisifera. In addition, the co-expression network analysis confirmed the significant correlation among these genes. Moreover, two simple sequence repeats (SSR) markers (S41 and S44) were developed to identify the seedless phenotype. The up-regulation of 4CL and MADS-box TFs activated the expression of PCD genes; on the other hand, S-RNase resulted in pollen tube RNA degradation and triggered PCD. While the link between PCD and seedless phenotype formation in oil palm has not been extensively studied to date, these findings suggest a role of PCD in pollen tube lethality, leading to double fertilization failure and the seedless phenotype.

Auxin Metabolism Is Involved in Fruit Set and Early Fruit Development in the Parthenocarpic Tomato "R35-P"

Parthenocarpic tomato can set fruit and develop without pollination and exogenous hormone treatments under unfavorable environmental conditions, which is beneficial to tomato production from late fall to early spring in greenhouses. In this study, the endogenous hormones in the ovaries of the parthenocarpic tomato line "R35-P" (stigma removed or self-pollination) and the non-parthenocarpic tomato line "R35-N" (self-pollination) at four stages between preanthesis and postanthesis investigated, using high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS). A nearly twofold IAA (indoleacetic acid) content was found in "R35-P" rather than in "R35-N" at -2 and 0 days after anthesis (DAA). Except at -2 DAA, a lower ABA (abscisic acid) content was observed in Pe (stigma removed in "R35-P") compared to that in Ps (self-pollination in "R35-P") or CK (self-pollination in "R35-N"). After pollination, although the content of GA₁ (gibberellins acid 1) in CK increased, the levels of GAs (gibberellins acids) were notably low. At all four stages, a lower SA (salicylic acid) content was found in Ps and CK than in Pe, while the content and the change trend were similar in Ps and CK. The variation tendencies of JA (jasmonic acid) varied among Pe, Ps, and CK at the studied periods. Furthermore, KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analyses of transcriptomic data identified 175 differentially expressed genes (DEGs) related to plant hormone signal transduction, including 63 auxin-related genes, 27 abscisic acid-related genes, 22 ethylene-related genes, 16 cytokinin-related genes, 16 salicylic acid-related genes, 14 brassinosteroid-related genes, 13 jasmonic acid-related genes, and 4 gibberellin-related genes at -2 DAA and 0 DAA. Our results suggest that the fate of a fruit set or degeneration occurred before anthesis in tomato. Auxins, whose levels were independent of pollination and fertilization, play prominent roles in controlling a fruit set in "R35-P," and other hormones are integrated in a synergistic or antagonistic way.

Decoding the molecular mechanism of parthenocarpy in Musa spp. through protein-protein interaction network

Banana, one of the most important staple fruit among global consumers is highly sterile owing to natural parthenocarpy. Identification of genetic factors responsible for parthenocarpy would facilitate the conventional breeders to improve the seeded accessions. We have constructed Protein-protein interaction (PPI) network through mining differentially expressed genes and the genes used for transgenic studies with respect to parthenocarpy. Based on the topological and pathway enrichment analysis of proteins in PPI network, 12 candidate genes were shortlisted. By further validating these candidate genes in seeded and seedless accession of Musa spp. we put forward MaAGL8, MaMADS16, MaGH3.8, MaMADS29, MaRGA1, MaEXPA1, MaGID1C, MaHK2 and MaBAM1 as possible target genes in the study of natural parthenocarpy. In contrary, expression profile of MaACLB-2 and MaZEP is anticipated to highlight the difference in artificially induced and natural parthenocarpy. By exploring the PPI of validated genes from the network, we postulated a putative pathway that bring insights into the significance of cytokinin mediated CLAVATA(CLV)-WUSHEL(WUS) signaling pathway in addition to gibberellin mediated auxin signaling in parthenocarpy. Our analysis is the first attempt to identify candidate genes and to hypothesize a putative mechanism that bridges the gaps in understanding natural parthenocarpy through PPI network.

Transcriptome Analysis of Eggplant Root in Response to Root-Knot Nematode Infection

Eggplant (Solanum melongena L.), which belongs to the Solanaceae family, is an important vegetable crop. However, its production is severely threatened by root-knot nematodes (RKNs) in many countries. Solanum torvum, a wild relative of eggplant, is employed worldwide as rootstock for eggplant cultivation due to its resistance to soil-borne diseases such as RKNs. In this study, to identify the RKN defense mechanisms, the transcriptomic profiles of eggplant and Solanum torvum were compared. A total of 5360 differentially expressed genes (DEGs) were identified for the response to RKN infection. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that these DEGs are mainly involved in the processes of response to stimulus, protein phosphorylation, hormone signal transduction, and plant-pathogen interaction pathways. Many phytohormone-related genes and transcription factors (MYB, WRKY, and NAC) were differentially expressed at the four time points (ck, 7, 14, and 28 days post-infection). The abscisic acid signaling pathway might be involved in plant-nematode interactions. qRT-PCR validated the expression levels of some of the DEGs in eggplant. These findings demonstrate the nematode-induced expression profiles and provide some insights into the nematode resistance mechanism in eggplant.

Isolation of high-quality RNA from recalcitrant leaves of variegated and resurrection plants

Resurrection plant Ramonda serbica is a suitable model to investigate mechanisms of desiccation tolerance, while variegated Pelargonium zonale has been proven to serve as an excellent model for the metabolite allocation between sink tissue and source tissue within the same organ. However, the genomes of these plants are still not sequenced, limiting their application in molecular studies. To investigate the transcript abundance by next-generation sequencing, high-quality RNA input is required. Leaves of both P. zonale and R. serbica are rich in polyphenols that interfere with high-quality RNA extraction by common protocols. Moreover, low water content and high amount of sugars and other osmoprotectants in desiccated R. serbica leaves present the additional challenge in total RNA extraction. Here, we evaluated and compared several already established TRIzol- and CTAB-based protocols aiming to develop the efficient, simple and low-cost methods for the extraction of the satisfactory yield RNA of great purity and integrity, required for the construction of high-quality cDNA libraries. Our results show that the CTAB-based protocol (i.e. CTAB 1b) enabled the extraction of high-quality RNA from photosynthetically active and non-photosynthetically active leaf sectors of P. zonale, with high RIN values. On the other hand, TRIzol-based protocol provided a high RNA yield with low contamination and high RNA integrity even in desiccated leaves of R. serbica. We envisage that the proposed protocol would be suitable for the RNA extractions from other desiccated organs (e.g. seeds, grains, pollen grains).

Heracleum sosnowskyi seed development under the effect of exogenous application of GA3

Numerous studies have demonstrated the impact of exogenous gibberellin on fleshy fruit formation, but the effect on dry fruits is not yet well known. To test the role of gibberellin (GA₃) in dry fruit formation, we analysed the impact of exogenous GA₃ on the invasive plant Sosnowsky’s hogweed (H. sosnowskyi Manden.) seed development and germination. Treatment of GA₃ concentrations of 0.07 mM, 0.14 mM, 0.28 mM, 0.43 mM was applied to flowers at the early stage of development. Seeds were collected from treated satellite umbels. It was observed that GA₃treatment did not have a significant effect on the size of H. sosnowskyi seeds, but caused various changes in their shape. The data on semi-thin longitudinal sections of H. sosnowskyi mericarps and SEM micrographs of embryos showed that the embryos in GA₃ (0.43 mM) treated variants were at torpedo stage, while in control variants—mature embryos. The germination of seeds of each variant was estimated by burying them in the soil. Our studies indicated that GA₃ application reduced the germination of H. sosnowskyi seed from 98.0% (control) to 16.5% (GA₃ concentration 0.43 mM). It was assumed that exogenous application of GA₃ had influence on the development of dry Sosnowsky’s hogweed seeds and could be used to inhibit the spread of this invasive plant.

De novo assembly of Persea americana cv. 'Hass' transcriptome during fruit development

BACKGROUND

Avocado (Persea americana Mill.) is a basal angiosperm from the Lauraceae family. This species has a diploid genome with an approximated size of ~ 920 Mbp and produces a climacteric, fleshy and oily fruit. The flowering and fruit set are particularly prolonged processes, lasting between one to three months, generating important differences in physiological ages of the fruit within the same tree. So far there is no detailed genomic information regarding this species, being the cultivar 'Hass' especially important for avocado growers worldwide. With the aim to explore the fruit avocado transcriptome and to identify candidate biomarkers to monitore fruit development, we carried out an RNA-Seq approach during 4 stages of 'Hass' fruit development: 150 days after fruit set (DAFS), 240 DAFS, 300 DAFS (harvest) and 390 DAFS (late-harvest).

RESULTS

The 'Hass' de novo transcriptome contains 62,203 contigs (x̅=988 bp, N50 = 1050 bp). We found approximately an 85 and 99% of complete ultra-conserved genes in eukaryote and plantae database using BUSCO (Benchmarking Universal Single-Copy Orthologs) and CEGMA (Core Eukaryotic Gene Mapping Approach), respectively. Annotation was performed with BLASTx, resulting in a 58% of annotated contigs (90% of differentially expressed genes were annotated). Differentially expressed genes analysis (DEG; with False Discovery Rate ≤ 0.01) found 8672 genes considering all developmental stages. From this analysis, genes were clustered according to their expression pattern and 1209 genes show correlation with the four developmental stages.

CONCLUSIONS

Candidate genes are proposed as possible biomarkers for monitoring the development of the 'Hass' avocado fruit associated with lipid metabolism, ethylene signaling pathway, auxin signaling pathway, and components of the cell wall.

Transcriptome analysis of oil palm inflorescences revealed candidate genes for an auxin signaling pathway involved in parthenocarpy

Oil palm parthenocarpic fruits, which are produced without fertilization, can be targeted to increase oil content because the majority of the fruit is occupied by mesocarp, the part in which palm oil is stored. Consequently, gaining an understanding of the parthenocarpic mechanism would be instrumental for producing parthenocarpic oil palm. This study aims to determine effects of auxin treatment and analyze differentially expressed genes in oil palm pistils at the pollination/anthesis stage, using an RNA sequencing (RNA seq) approach. The auxin treatment caused 100% parthenocarpy when auxin was sprayed before stigmas opened. The parthenocarpy decreased to 55%, 8% and 5% when the auxin was sprayed 1, 2 and 3 days after the opening of stigmas, respectively. Oil palm plants used for RNA seq were plants untreated with auxin as controls and auxin-treated plants on the day before pollination and 1 day after pollination. The number of raw reads ranged from 8,425,859 to 11,811,166 reads, with an average size ranging from 99 to 137 base pairs (bp). When compared with the oil palm transcriptome, the mapped reads ranged from 8,179,948 to 11,320,799 reads, representing 95.85–98.01% of the oil palm matching. Based on five comparisons between RNA seq of treatments and controls, and confirmation using reverse transcription polymerase chain reaction and quantitative real-time RT-PCR expression, five candidate genes, including probable indole-3-acetic acid (IAA)-amido synthetase GH3.8 (EgGH3.8), IAA-amido synthetase GH3.1 (EgGH3.1), IAA induced ARG7 like (EgARG7), tryptophan amino transferase-related protein 3-like (EgTAA3) and flavin-containing monooxygenase 1 (EgFMO1), were differentially expressed between auxin-treated and untreated samples. This evidence suggests a pathway of parthenocarpic fruit development at the beginning of fruit development. However, more research is needed to identify which genes are definitely involved in parthenocarpy.

The making of virgin fruit: the molecular and genetic basis of parthenocarpy

Fruit set-the commitment of an angiosperm plant to develop fruit-is a key developmental process that normally occurs following successful fertilization. Parthenocarpy arises when fruit automatically develop in the absence of fertilization. This review uses parthenocarpic fruit development as a focal device through which to recapitulate and understand the molecular effectors that mediate and regulate fruit set. The review demonstrates that studies of parthenocarpy are providing vital insight into plant development, signaling and, potentially, high-value agricultural products.