Qingfei Li1,2
Weili Guo1,2Bihua Chen1,2Feifei Pan1,2Helian Yang1,2Junguo Zhou1,2Guangyin Wang1,2Xinzheng Li1,2*- 1College of Horticulture and Landscape, Henan Institute of Science and Technology, Xinxiang, China
- 2Henan Province Engineering Research Center of Horticultural Plant Resource Utilization and Germplasm Enhancement, Xinxiang, China
The number and proportion of female flowers per plant can directly influence the yield and economic benefits of cucurbit crops. Ethephon is often used to induce female flowers in cucurbits. However, the mechanism through which it affects floral sex differentiation in pumpkin is unknown. We found that the application of ethephon on shoot apical meristem of pumpkin at seedling stage significantly increased the number of female flowers and expedited the appearance of the first female flower. These effects were further investigated by transcriptome and hormone analyses of plants sprayed with ethephon. A total of 647 differentially expressed genes (DEGs) were identified, among which 522 were upregulated and 125 were downregulated. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analysis indicated that these genes were mainly enriched in plant hormone signal transduction and 1-aminocyclopropane-1-carboxylate oxidase (ACO). The results suggests that ethylene is a trigger for multiple hormone signaling, with approximately 4.2% of the identified DEGs involved in ethylene synthesis and multiple hormone signaling. Moreover, ethephon significantly reduced the levels of jasmonic acid (JA), jasmonoyl-L-isoleucine (JA-ILE), and para-topolin riboside (pTR) but increased the levels of 3-indoleacetamide (IAM). Although the level of 1-aminocyclopropanecarboxylic acid was not changed, the expression of ACO genes, which code for the enzyme catalyzing the key rate-limiting step in ethylene production, was significantly upregulated after ethephon treatment. The results indicate that the ethephon affects the transcription of ethylene synthesis and signaling genes, and other hormone signaling genes, especially auxin responsive genes, and modulates the levels of auxin, jasmonic acid, and cytokinin (CK), which may together contribute to femaleness.
Introduction
Pumpkin (Cucurbita moschata Duch.) is cultivated worldwide and is popular for its high nutritional and medicinal value (Lu et al., 2019). As a monoecious plant, pumpkin is a typical material for exploration of floral sex differentiation. The sex differentiation in the flowers of pumpkin includes three stages. First, male flowers appear at the base nodes of the plant, and almost no female flowers are generated at this stage. In the second stage, female flowers alternate with male flowers; usually, one female flower appears after several male flowers. In the third stage, which occurs at the end of the blooming season, female flowers appear continuously, but these are not suitable for generating fruits. Overall, male flowers are much more abundant than female flowers in pumpkin. However, for higher yield, it is necessary to have more female flowers per plant.
Sexual differentiation of flowers is mainly influenced by sex determining genes, hormones, and environmental factors. Sex expression in cucurbitaceae crops can be affected by several phytohormones, such as ethylene, auxin, gibberellin, cytokinin, abscisic acid (ABA), brassinosteroid, and salicylic acid (SA; Byers et al., 1972; Rudich and Halevy, 1974; Trebitsh et al., 1987; Menéndez et al., 2009; Pimenta Lange and Lange, 2016; Zhang et al., 2017). Among these, ethylene is the main regulator of sex determination in cucurbitaceae (Manzano et al., 2011).
Sex differentiation in cucumber, which has been studied more intensively among cucurbitaceae crops, is mainly determined by the F (CsACS1G), M (CsACS2), and A (CsACS11) genes. These genes encode 1-aminocyclopropane-1-carboxylate synthase (ACS), a key rate-limiting enzyme in the biosynthesis of ethylene (Pan et al., 2018). The F gene promotes femaleness (Mibus and Tatlioglu, 2004; Knopf and Trebitsh, 2006), the M gene inhibits the development of stamens (Yamasaki et al., 2001; Saito et al., 2007), and the A gene is an androecious gene (plants with mutations in CsACS11 do not have female flowers; Boualem et al., 2015). 1-aminocyclopropane-1-carboxylate oxidase (ACO) is another key enzyme in the ethylene biosynthesis pathway (Adams and Yang, 1979; Houben and Van de Poel, 2019). Organ-specific overexpression of cucumber CsACO2 was reported to arrest the development of stamens in Arabidopsis. Among floral organs, stamens are the most sensitive to exogenous ethylene, and their development can be arrested by endogenous ethylene for inducing female flowers (Duan et al., 2008). Recently, it was shown that ACO is expressed in the carpel primordia and is required for the development of carpel in cucumber. Cucumber plants having a mutation in CsACO2 bear only male flowers because of impaired enzymatic activity of ACO and reduced emission of ethylene. In addition, a transcription factor, CsWIP1, which is negatively correlated with the formation of female flowers, could repress the expression of CsACO2 by binding to its promoter (Chen et al., 2016). These findings suggest that ACO is indispensable for the development of female flowers.
Besides the ethylene biosynthesis genes (ACS, ACO), many genes related to ethylene signaling have recently been reported to be involved in sex differentiation. Sex expression is a complex process. In a comparative transcriptome analysis of shoot apices from male, female, and hermaphroditic lines of cucumber, hormone synthesis and signaling, and ion homeostasis, which is important for ethylene perception and signaling, were found to be involved in sex differentiation (Pawełkowicz et al., 2019). CpETR1A and CpETR2B are ethylene receptor genes in Cucurbita pepo, which control the ethylene response; mutants in CpETR1A and CpETR2B are ethylene-insensitive and exhibit conversion from monoecy to andromonoecy (García et al., 2019). The ethylene-receptor gene, CsETR1, expressed in the pistil primordia, is involved in the arrest of stamen development by inducing DNA damage in primordial anthers of female flowers (Hao et al., 2003; Yamasaki et al., 2003; Duan et al., 2008; Wang et al., 2010). The mRNA levels of Cs-ETR2 and CS-ERS were significantly enhanced after ethrel application and were decreased upon application of an ethylene inhibitor (Yamasaki et al., 2000). We previously performed a comparative analysis of the transcriptomes of aborted and normal pistils, and showed that ethylene signal transduction genes are implicated in the development of pistils in C. moschata (Li et al., 2020).
Flower development involves interaction of multiple hormones. Besides ethylene, auxin is another important regulator of flower development, and the signaling pathways of the two are suggested to crosstalk (Muday et al., 2012). Auxin response factors (ARFs) were shown to be indispensable for the development of pistils in Japanese apricot (Song et al., 2015). CpARF is highly expressed during the early stages of flower development. Many auxin response elements (AuxREs) are present in the promoters of ethylene signaling (CpETR) and biosynthesis (CpACS, CpACO) genes (Liu et al., 2015). Ethylene stimulates the biosynthesis of auxin and its transport toward the elongation zone in the root tip; root cell expansion is, therefore, a result of coordinated actions of ethylene and auxin pathways (Vanstraelen and Benková, 2012).
Ethephon is a widely used ethylene-releasing agent in agriculture. When applied to plants, ethephon acts via release of ethylene, which can interfere in the growth process of plants (Cooke and Randall, 1968; Sun et al., 2015). In the field production of cucurbits, ethephon is used to induce more female flowers to obtain higher yields. Treatment with 50mg/L ethrel at the third leaf stage has long been known to increase the femaleness of cucumber (Iwahori et al., 1970). The number of female flowers per plant was increased to varying degrees upon treatment with ethrel in the concentration range from 50 to 250mg/L (Jin et al., 2011). As described above, in cucumber, ethylene may arrest the development of stamens and promote the generation of female flowers. In Ficus carica, ethephon treatment induced changes at the molecular level; for example, it significantly upregulated ACO and ARF and downregulated most of the ERF and PAL genes (Cui et al., 2021). Although ethephon has been used to regulate floral sex differentiation in the production of cucurbits for many years, little is known about the underlying regulatory mechanisms in pumpkin. Whether ethephon affects floral sex differentiation in pumpkin by regulating gene transcription associated with floral development or by modulating the endogenous levels of hormones and which hormones respond to ethephon is unclear. In the present study, the application of ethephon was observed to induce more female flowers in pumpkin and expedited the appearance of the first female flower. To further explore the regulatory mechanism behind the effect of ethephon in sex differentiation, we investigated the changes in the transcriptome and levels of endogenous hormones in the shoot apical meristem after external application of ethephon. Our results indicate that the early flowering and the greater number of female flowers might be the result of increased level of the ACO transcript and alteration in the expression of hormone signaling genes and endogenous hormones levels.
Materials and Methods
Plant Materials and Sample Preparation
An inbred line of pumpkin (C. moschata Duch), 009-1, was used as the plant material. The seeds were sterilized with hot water (55°C) for 10–15min, kept soaked for 4–5h, and germinated using soaked germination paper, in petri dishes at 28°C. When 90% of the seeds were germinated, they were transferred into a mixed matrix (peat:vermiculite:perlite, 3:1:0.5, v/v) in pots (diameter, 34cm) and grown under natural light.
Ethephon Application, Plant Sampling, and Floral Sex Differentiation
Shoot apical meristems were treated twice with ethephon (100mg/L; E8021, Solarbio, Beijing, China) containing 0.1% Tween-20, every 2days, when the third true leaf of seedlings unfolded. Tween-20 (0.1%) was used as a blank control. The concentration of ethephon was chosen according to Yang et al. (2015). After the first exposure to ethephon for 4h, the shoot apical meristems were immediately frozen in liquid nitrogen and stored at −80°C for RNA-sequencing (RNA-seq) and quantitative real-time PCR (qRT-PCR) validation. The results of a preliminary experiment showed no significant effect on gene expression within 4h of ethephon treatment. Hormone levels were analyzed a day after the second treatment. The shoot apical meristems from three plants were pooled as one biological replicate for both ethephon and control treatments. RNA-seq analysis and hormone quantification were performed using three biological replicates. The tissue samples were immediately frozen in liquid nitrogen and then stored at −80°C.
At the flowering stage, the number of female and male flowers within 20 nodes and the node at which the first female flower occurred were determined. The nodes of plants were marked with a red line after every observation, to prevent missing or duplicating the count of flowers. Each group had three plants, each with three biological replicates.
RNA Sequencing and Analysis
RNA was isolated from shoot apices subjected to ethephon and control treatments using the TRIzol™ reagent (Invitrogen, Carlsbad, CA, United States). The concentration of RNA samples was measured using NanoDrop DU8000 (Thermo, CA, United States), and their purity and integrity were evaluated by agarose gel electrophoresis and using the Agilent 2100 system (Agilent Technologies, CA, United States). Qualified RNA samples were used for construction of cDNA libraries as described previously (Li et al., 2020). The high-quality libraries were used for paired-end sequencing (2×150bp) on the Illumina NovaSeq 6000 System.
The raw reads were processed to filter out the adaptor sequences and low-quality reads (more than 50% bases with SQ≤20 in one read and with more than 10% N bases). The Q20, Q30, GC-content, and sequence duplication level of the clean data were calculated. The clean reads were then mapped to the C. moschata genome by using TopHat2 (Trapnell et al., 2012; Kim et al., 2013; Sun et al., 2017), allowing up to one mismatch. The DESeq R package (1.10.1) was used to identify the differentially expressed genes (DEGs) (Love et al., 2014). The fragments per kilobase of exon per million fragments mapped (FPKM) method was used to estimate the expression levels of genes. Values of p were adjusted using the Benjamini and Hochberg’s method for controlling the false discovery rate (FDR; Benjamini and Hochberg, 1995). Genes with an adjusted value of p<0.05 and a fold change ≥1.5 based on three biological replicates were considered differentially expressed.
Functional annotation of genes was based on the following databases: the NCBI non-redundant (Nr) protein sequences, Swiss–Prot, clusters of orthologous groups of proteins (KOG), protein family (Pfam), gene ontology (GO), and the Kyoto encyclopedia of genes and genomes (KEGG).
Functional Enrichment Analysis
Gene ontology enrichment analysis of DEGs was implemented using the GOSeq R package (Young et al., 2010), in which gene length bias was corrected, and a p value of DEGs ≤0.05 was considered as significantly enriched. KEGG pathway enrichment analysis of DEGs was performed using the software KOBAS (Mao et al., 2005). Pathways with their Benjamini and Hochberg adjusted values of p≤0.05 were defined as significantly enriched by DEGs.
Quantitative Real-Time RT-PCR
The first-strand cDNA was obtained using the PrimeScript™ RT Master Mix (Perfect Real Time) Reagent Kit (Takara, Dalian, China). The qRT-PCR was carried out using a Bio-Rad IQ5 instrument (Foster City, CA, United States), as follows: 95°C for 40s and 40cycles of 95°C for 5s and 61°C for 30s. ACTIN was used as an internal control. The primers used for qRT-PCR are listed in Supplementary Material. The expression levels were calculated using the 2-ΔΔCt method (Livak and Schmittgen, 2001). The expression of each gene was determined using three biological and three technical replicates. Furthermore, the correlation analysis and the Pearson correlation coefficient between the log2 (fold change) values obtained in qRT-PCR and RNA-seq were calculated using the IBM SPSS statistics 22 software.
Quantification of Hormones
The quantification of endogenous auxin, abscisic acid (ABA), jasmonic acid (JA), cytokinin (CK), gibberellic acid (GA), salicylic acid (SA), and 1-aminocyclopropane 1-carboxylic acid (ACC) was performed. Fresh shoot apical meristems (50mg) were frozen in liquid nitrogen and extracted with 1ml methanol/water/formic acid (15:4:1, V/V/V). The extracts were evaporated to dryness under nitrogen gas, reconstituted in 100μl 80% (V/V) methanol, and filtered through a 0.22μm filter for LC-MS analysis. The extracts were analyzed using a UPLC-ESI-MS/MS system (UPLC, ExionLC™ AD, MS, Applied Biosystems 6500 Triple Quadrupole). The analytical conditions were as follows: for LC: column, Waters ACQUITY UPLC HSS T3 C18 (1.8μm, 100mm×2.1mm i.d.); solvent system, water with 0.04% acetic acid (A), acetonitrile with 0.04% acetic acid (B); flow rate, 0.35ml/min; temperature, 40°C; injection volume, 2μl; for MS/MS: AB 6500+ QTRAP® LC–MS/MS System, equipped with an ESI Turbo Ion-Spray interface, operating in both positive and negative ion modes and controlled using the Analyst 1.6 software (AB Sciex). Solutions containing different concentrations (0.01, 0.05, 0.1, 0.5, 1, 5, 10, 50, 100, 200, and 500ng/ml) of 1-aminocyclopropanecarboxylic acid, indole-3-acetyl-L-aspartic acid, 3-indoleacetamide (IAM), para-topolin riboside (pTR), jasmonic acid, jasmonoyl-L-isoleucine (JA-ILE), cis(+)-12-oxophytodienoic acid, and gibberellin A15 were used to generate standard curves for each hormone.
Results
Effect of Ethephon on Floral Sex Differentiation in Pumpkin
The floral sex differentiation in pumpkin plants treated with ethephon was investigated. Ethephon treatment significantly expedited the appearance of the first female flower from node 13.29±1.89 (in the control) to node 8.29±1.50. The number of female flowers within 20 nodes in ethephon-treated plants was significantly higher (3.86±0.99) than in the control (2.00±0.58). The effect of ethephon treatment on the number of male flowers was not obvious (Figures 1A,B).
Figure 1. Effect of ethephon on floral sex differentiation of pumpkin. (A) Node of the first female flower and the number of female flower and male flower in ethephon-treated and control plants. Double asterisk indicates a significant difference between control and ethephon-treated plants at 1% level as determined by Duncan’s new multiple range test. (B) The schematic diagram of sex differentiation of ethephon-treated and control plants.
Statistical Analysis of RNA-seq Data and DEGs Induced by Ethephon
An average of 45,929,596 reads was obtained in the sequencing of the cDNA library prepared from each sample, with an average Q30 quality score≥93.30% and average Q20 quality score≥97.52%. The alignment of filtered reads with the C. moschata genome sequence revealed an average mapping percentage of 95.06% (Table 1). A total of 647 DEGs, including 522 upregulated and 125 downregulated genes, were identified by comparing the transcriptomes of the shoot apical meristem of ethephon-treated and control plants (Figure 2).
Table 1. Statistical analysis of reads mapped to the reference genome after rRNA filtering.
Figure 2. Analysis of differentially expressed genes (DEGs) between ethephon-treated and control plants. (A) Numbers of upregulated and downregulated DEGs. Red and green colors indicate up- and downregulated transcripts, respectively. (B) Volcano map of DEGs.
Functional Enrichment of DEGs
Gene ontology analysis indicated that the annotated genes were enriched in three major functional categories: biological processes, cellular components, and molecular functions. In the biological process category, most of the transcripts were enriched in metabolic processes, cellular process, and single-organism process (Figure 3A). In this category, the most significantly enriched GO terms included response to ethylene (GO:0009723), abscisic acid-activated signaling pathway (GO:0009738), and chitin catabolic process (GO:0006032; Figure 3B). In the cellular component category, DEGs were mainly enriched in cell, cell part, organelle, and membrane (Figure 3A), among which the SCF ubiquitin ligase complex (GO:0019005) was the most significantly enriched GO term (Figure 3B). In the molecular function category, catalytic activity and binding were highly enriched GO terms (Figure 3A), and the most significantly enriched GO terms were ACO and chitin binding (Figure 3B). More DEGs were enriched in the biological processes and cellular components categories, and relatively fewer DEGs were enriched in the molecular functions category.
Figure 3. Gene ontology (GO) analysis of DEGs between control and ethephon-treated shoot apical meristem of pumpkin. (A) GO categories. (B) Enrichment analysis of DEGs in biological processes, cellular components, and molecular functions category. The number of DEGs enriched in each term is marked on the left.
Kyoto encyclopedia of genes and genomes pathway analysis was performed to investigate the pathways that responded to ethephon treatment. A total of 116 DEGs were mapped to 72 KEGG pathways. The top 20 pathways in KEGG enrichment analysis are as shown in Figure 4. Among them, DEGs were significantly enriched in cysteine and methionine metabolism pathways and plant hormone signal transduction pathway (Figure 5). Notably, all the 12 DEGs enriched in the cysteine and methionine metabolism pathways were upregulated. Five of them were annotated as predicted ACO genes (Table 2). Twenty DEGs, including 19 upregulated and one downregulated genes, were enriched in plant hormone signal transduction pathway. Among them, 10 DEGs were involved in ethylene response and four were involved in auxin response and induction, according to the annotations in the Nr and Swiss-Prot databases (Table 2). Except for one auxin-induced gene, all other DEGs significantly enriched in the plant hormone signal transduction and cysteine and methionine metabolism pathways were upregulated upon ethephon treatment.
Figure 4. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis of DEGs.
Figure 5. Heatmap diagram of expression levels of DEGs in significant enrichment KEGG pathways. (A) DEGs enriched in plant hormone signal transduction pathway. (B) DEGs enriched in cysteine and methionine metabolism pathway. The red and green colors indicate high and low expression levels (log2 FPKM), respectively.
Table 2. List of DEGs enriched in plant hormone signal transduction pathway (No. 1–20) and cysteine and methionine metabolism pathway (No. 21–32).
Validation of RNA-Seq Results by qRT-PCR
The validation of the transcriptome data was done by qRT-PCR analysis. Twelve DEGs were chosen for validation. The expression levels of these genes determined by qRT-PCR were in good agreement with the RNA-seq data, with relative coefficient, R2=0.804 (Figures 6A,B), and the Pearson correlation coefficient, R=0.897 (p<0.0001). These results indicate the reliability of the RNA-seq analysis performed in this study.
Figure 6. Quantitative real-time PCR (qRT-PCR) verification. (A) Relative expression levels of 12 selected DEGs between ethephon-treated and control plants. Error bars indicate the SEs. (B) Correlation analysis of RNA-seq log2 (fold change) and qRT-PCR log2 (fold change) of 12 selected genes: AD1-0960 (CmoCh09G010960.1), ETHYLENE INSENSITIVE 3 (EIN3)-3100 (CmoCh05G013100.1), EIN3-4690 (CmoCh09G004690.1), AUX-IAA-6330 (CmoCh10G006330.1), EIN3-2320 (CmoCh12G012320.1), AUX-22D-1500 (CmoCh15G011500.1), AUX-SAUR-7170 (CmoCh19G007170.1), ACO-640 (CmoCh02G000640.1), ACO-5880 (CmoCh03G005880.1), ACO-3260 (CmoCh07G003260.1), ACO-150 (CmoCh11G020150.1), and MK-3620 (CmoCh16G003620.1).
Changes in the Endogenous Levels of Hormones Upon Ethephon Treatment
The endogenous levels of hormones, viz., auxin, ABA, JA, CK, GA, SA, and ACC, in the shoot apical meristem of pumpkin plants subjected to ethephon and control treatments were compared to assess the effects of ethephon. A fold change≥2 was considered to indicate differential hormone levels (Table 3). Ethephon treatment decreased the levels of JA, JA-ILE, and pTR and significantly increased the levels of 3-IAM (Table 3). In addition, the levels of cis(+)-12-oxophytodienoic acid, gibberellin A15, and indole-3-acetyl-L-aspartic acid were slightly decreased, while those of 1-aminocyclopropanecarboxylic acid were not effected by the treatment.
Table 3. Phytohormone levels in the shoot apical meristem of control and ethephon-treated plants.
Discussion
Ethephon is a plant growth regulator that is mainly used in field production for regulating floral sex differentiation (Papadopoulou et al., 2005; Martínez et al., 2013), promoting fruit ripening (Cui et al., 2021), and breaking dormancy (Corbineau et al., 2014). The effects of ethephon on the sex differentiation phenotype of flowers in zucchini, cucumber, and watermelon are well-known, but the underlying mechanism remains unclear. Moreover, there has been no systematic research on the effect of ethephon on the floral sex expression of pumpkin. In the present study, 100mg/l ethephon was used to treat shoot apical meristem of pumpkin at the seedling stage. The ethephon treatment significantly advanced the appearance of the first female flowers and significantly increased the number of female flowers within 20 nodes. To investigate the mechanism underlying these effects, the changes in the transcriptome and hormone levels upon ethephon treatment were analyzed.
Transcriptome analysis indicated an upregulation of seven DEGs annotated as ACO, which encode a key enzyme in ethylene biosynthesis. Study have found that CsACO2 mutants of cucumber bear only male flowers because of impairment of the enzymatic activity of ACO and reduced emission of ethylene (Chen et al., 2016). ACO can sometimes be rate limiting in ethylene biosynthesis (Houben and Van de Poel, 2019). Overexpression of ACO from Vitis vinifera in tomato was reported to increase the rate of ethylene release (Cai et al., 2013). The results of the above studies indicate that ACO is indispensable for the stable emission of ethylene and the development of female flowers. The upregulation of ACO in ethephon-treated plants mean that they may have a greater capacity to produce more ethylene.
Ethylene response, including ethylene receptor, ethylene insensitive 3, and ethylene-responsive transcription factor, was found to be the most significantly enriched GO term (Figure 3B). Ethephon treatment upregulated the expression of ethylene response-related genes, which might be due to the increased release of ethylene. ETR controls the ethylene response, and mutations in CpETR1A and CpETR2B result in ethylene-insensitivity and conversion of monoecy into andromonoecy (García et al., 2019). In cucumber, CsETR1 is localized in the pistil primordia and is involved in arresting the development of stamen in female flowers (Hao et al., 2003; Yamasaki et al., 2003; Duan et al., 2008; Wang et al., 2010). The application of ethephon increased the transcript levels of ethylene receptor, ethylene-responsive transcription factor, ETHYLENE INSENSITIVE 3 (EIN3)-binding F-box protein 1, auxin-responsive, indole-3-acetic acid-amido synthetase, and abscisic acid receptor genes. These results for pumpkin are in agreement with those of a previous study in which ACOs and IAA-amino acid hydrolase, indole-3-acetic acid-amido synthetase, auxin-responsive protein, and ABA receptor genes were found to be upregulated upon ethephon treatment (Cui et al., 2021).
Ethylene treatment induces a cascade of regulatory events. In the present study, we observed that five of the 20 DEGs included in plant hormone signal transduction pathways were annotated as EIN3-binding F-box protein and EIN3-like 1 protein, which were upregulated upon ethephon treatment. EIN3 and EIN3-LIKE1 (EIL1) are key transcription factors for ethylene signaling. EIN3 accumulates in the nucleus in the presence of ethylene. In contrast, in the absence of ethylene, EIN3 is negatively regulated and constantly degraded in plant cells (Cho and Yoo, 2015). The complex regulation of the activation of EIN3 and EIL1 in response to ethylene involves triggering of primary transcription through EIN3-binding sites in the promoters of ETHYLENE RESPONSE FACTOR1 (ERF1), EBF2, ERS1, ERS2, and ETR2, that finally activates the ethylene biosynthesis genes (ACOs and ACSes) Vandenbussche et al., 2012; Cho and Yoo, 2015; Dolgikh et al., 2019). The upregulation of genes coding EIN3-binding F-box protein and EIN3-like 1 protein may contribute to the changes of transcriptional levels of ETR and ACO genes in ethephon-treated plants. Besides the classical mechanism of ethylene-induced stabilization of EIN3/EIL1, JA can release EIN3/EIL1 from repression by accelerating the degradation of JAZ, thereby inducing ethylene responses (Zhu et al., 2011). In our study, the levels of JA and JA-ILE were significantly decreased; however, whether this is related to the upregulation of EIN3/EIL1 needs to be investigated.
Auxin and ethylene act synergistically to regulate the growth and development of plants. An increase in the levels of ethylene was reported to elevate the auxin response, monitored using auxin-inducible reporters, in the root elongation zone (Muday et al., 2012). Ethylene may positively regulate auxin synthesis. The levels of free IAA were reported to increase in the root tip upon treatment with ACC (100mM; Růžička et al., 2007). We observed that the levels of the auxin-responsive, indole-3-acetic acid-amido synthetase transcript, and IAM were increased upon ethephon treatment. These results support the notion that ethylene enhances the transport of auxin in the elongation zone, which leads to elevated IAA levels (Muday et al., 2012).
In summary, ethephon treatment enhances the synthesis of ethylene by increasing the expression of ACO genes, and promotes ethylene signaling, which may further crosstalks with an upregulation of auxin responsive genes and the increase of auxin levels. The regulatory network of ethephon in female flowers of pumpkins is fairly complex, which interferes transcriptional and hormonal levels, and in this process ethylene and EIN3 play the intermediate pivotal role. This study have confirmed the effects of ethephon on floral sex differentiation in pumpkin and have presented the mechanism through which ethephon promotes femaleness. Further genetic and biochemical analysis would help in clarifying the regulatory mechanism of ethephon on floral sex differentiation.
Data Availability Statement
The raw transcriptome data have been deposited in the National Center for Biotechnology Information (NCBI) SRA database (accession number: PRJNA736171).
Author Contributions
QL conducted the experiments and data analysis. WG, BC, FP, and HY assisted with plant material preparation. XL and QL conceived the project. QL wrote the original manuscript. JZ, GW, and XL supervised the studies and reviewed the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Key Science and Technology Program of Henan Province (212102110128), the Major Public Welfare Projects of Henan Province (201300111300), and the National Natural Science Foundation of China (U2004161 and 31401876).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.715487/full#supplementary-material
References
Adams, D. O., and Yang, S. F. (1979). Ethylene biosynthesis: identification of 1-aminocyclopropane-1-carboxylic acid as an intermediate in the conversion of methionine to ethylene. Proc. Natl. Acad. Sci. U. S. A. 76, 170–174. doi: 10.1073/pnas.76.1.170
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B Stat. Methodol. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Boualem, A., Troadec, C., Camps, C., Lemhemdi, A., Morin, H., Sari, M. A., et al. (2015). A cucurbit androecy gene reveals how unisexual flowers develop and dioecy emerges. Science 350, 688–691. doi: 10.1126/science.aac8370
Byers, R. E., Baker, L. R., Sell, H. M., Herner, R. C., and Dilley, D. R. (1972). Ethylene: a natural regulator of sex expression of Cucumis melo L. Proc. Natl. Acad. Sci. 69, 717–720. doi: 10.1073/pnas.69.3.717
Cai, T., Ye, X., Ma, Y. C., Wang, X. L., Chen, F., and Cheng, Z. M. (2013). Over-expression of 1-aminocyclopropane-1-carboxylic acid oxidase genes (ACO) from Vitis vinifera in tomato (solanum locopersicum) and its effects on ethylene release rates. J. Agric. Biotechnol. 21, 1037–1044. doi: 10.3969/j.issn.1674-7968.2013.09.004
Chen, H. M., Sun, J. J., Li, S., Cui, Q. Z., Zhang, H. M., Xin, F. J., et al. (2016). An ACC oxidase gene essential for cucumber carpel development. Mol. Plant 9, 1315–1327. doi: 10.1016/j.molp.2016.06.018
Cho, Y. H., and Yoo, S. D. (2015). Novel connections and gaps in ethylene signaling from the ER membrane to the nucleus. Front. Plant Sci. 5:733. doi: 10.3389/FPLS.2014.00733
Cooke, A. R., and Randall, D. I. (1968). 2-Haloethanephosphonic acids as ethylene releasing agents for the induction of flowering in pineapples. Nature 218, 974–975. doi: 10.1038/218974a0
Corbineau, F., Xia, Q., Bailly, C., and El-Maarouf-Bouteau, H. (2014). Ethylene, a key factor in the regulation of seed dormancy. Front. Plant Sci. 5:539. doi: 10.3389/fpls.2014.00539
Cui, Y. Y., Zhai, Y. L., Flaishman, M., Li, J. P., Chen, S. W., Zheng, C. L., et al. (2021). Ethephon induces coordinated ripening acceleration and divergent coloration responses in fig (Ficus carica L.) flowers and receptacles. Plant Mol. Biol. 105, 347–364. doi: 10.1007/s11103-020-01092-x
Dolgikh, V. A., Pukhovaya, E. M., and Zemlyanskaya, E. V. (2019). Shaping ethylene response: the role of EIN3/EIL1 transcription factors. Front. Plant Sci. 10:1030. doi: 10.3389/fpls.2019.01030
Duan, Q. H., Wang, D. H., Xu, Z. H., and Bai, S. N. (2008). Stamen development in Arabidopsis is arrested by organ-specific overexpression of a cucumber ethylene synthesis gene CsACO2. Planta 228, 537–543. doi: 10.1007/s00425-008-0756-7
García, A., Aguado, E., Martínez, C., Loska, D., Beltrán, S., Valenzuela, J. L., et al. (2019). The ethylene receptors CpETR1A and CpETR2B cooperate in the control of sex determination in Cucurbita pepo. J. Exp. Bot. 71, 154–167. doi: 10.1093/jxb/erz417
Hao, Y. J., Wang, D. H., Peng, Y. B., Bai, S. L., Xu, L. Y., Li, Y. Q., et al. (2003). DNA damage in the early primordial anther is closely correlated with stamen arrest in the female flower of cucumber (Cucumis sativus L.). Planta 217, 888–895. doi: 10.1007/s00425-003-1064-x
Houben, M., and Van de Poel, B. (2019). 1-Aminocyclopropane-1-carboxylic acid oxidase (ACO): the enzyme that makes the plant hormone ethylene. Front. Plant Sci. 10:695. doi: 10.3389/fpls.2019.00695
Iwahori, S., Lyons, J. M., and Smith, O. E. (1970). Sex expression in cucumber plants as affected by 2-Chloroethylphosphonic acid, ethylene, and growth regulators. Plant Physiol. 46, 412–415. doi: 10.1104/pp.46.3.412
Jin, X. X., Qin, Z. W., Wu, T., and Zhou, X. Y. (2011). Identification of ethylene-responsive genes in ethrel-treated shoot apices of cucumber by suppression subtractive hybridization. Plant Mol. Biol. Report. 29, 875–884. doi: 10.1007/s11105-011-0304-7
Kim, D., Pertea, G., Trapnell, C., Pimentel, H., Kelley, R., and Salzberg, S. L. (2013). Tophat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14:R36. doi: 10.1186/gb-2013-14-4-r36
Knopf, R. R., and Trebitsh, T. (2006). The female-specific CS-ACS1g gene of cucumber. A case of gene duplication and recombination between the non-sex-specific 1-aminocyclopropane-1-carboxylate synthase gene and a branched-chain amino acid transaminase gene. Plant Cell Physiol. 47, 1217–1228. doi: 10.1093/pcp/pcj092
Li, Q. F., Zhang, L., Pan, F. F., Guo, W. L., Chen, B. H., Yang, H. L., et al. (2020). Transcriptomic analysis reveals ethylene signal transduction genes involved in pistil development of pumpkin. PeerJ 8:e9677. doi: 10.7717/peerj.9677
Liu, K. D., Yuan, C. C., Li, H. L., Lin, W. H., Yang, Y. J., Shen, C. J., et al. (2015). Genome-wide identification and characterization of auxin response factor (ARF) family genes related to flower and fruit development in papaya (Carica papaya L.). BMC Genomics 16:901. doi: 10.1186/s12864-015-2182-0
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8
Lu, A. X., Yu, M. G., Fang, Z. Y., Xiao, B., Guo, L., Wang, W. M., et al. (2019). Preparation of the controlled acid hydrolysates from pumpkin polysaccharides and their antioxidant and antidiabetic evaluation. Int. J. Biol. Macromol. 121, 261–269. doi: 10.1016/j.ijbiomac.2018.09.158
Manzano, S., Martínez, C., Megías, Z., Gómez, P., Garrido, D., and Jamilena, M. (2011). The role of ethylene and brassinosteroids in the control of sex expression and flower development in Cucurbita pepo. Plant Growth Regul. 65, 213–221. doi: 10.1007/s10725-011-9589-7
Mao, X. Z., Cai, T., Olyarchuk, J., and Wei, L. P. (2005). Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 21, 3787–3793. doi: 10.1093/bioinformatics/bti430
Martínez, C., Manzano, S., Megías, Z., Garrido, D., Pico, B., and Jamilena, M. (2013). Involvement of ethylene biosynthesis and signalling in fruit set and early fruit development in zucchini squash (Cucurbita pepo L.). BMC Plant Biol. 13:139. doi: 10.1186/1471-2229-13-139
Menéndez, V., Revilla, M. A., Fal, M. A., and Fernández, H. (2009). The effect of cytokinins on growth and sexual organ development in the gametophyte of Blechnum spicant L. Plant Cell Tissue Organ Cult. 96, 245–250. doi: 10.1007/s11240-008-9481-y
Mibus, H., and Tatlioglu, T. (2004). Molecular characterization and isolation of the F/f gene for femaleness in cucumber (Cucumis sativus L.). Theor. Appl. Genet. 109, 1669–1676. doi: 10.1007/s00122-004-1793-7
Muday, G. K., Rahman, A., and Binder, B. M. (2012). Auxin and ethylene: collaborators or competitors? Trends Plant Sci. 17, 181–195. doi: 10.1016/j.tplants.2012.02.001
Pan, J., Wang, G., Wen, H., Du, H., Lian, H., He, H., et al. (2018). Differential gene expression caused by the F and M loci provides insight into ethylene-mediated female flower differentiation in cucumber. Front. Plant Sci. 9:1091. doi: 10.3389/fpls.2018.01091
Papadopoulou, E., Little, H. A., Hammar, S. A., and Grumet, R. (2005). Effect of modified endogenous ethylene production on sex expression, bisexual flower development and fruit production in melon (Cucumis melo L.). Sex. Plant Reprod. 18, 131–142. doi: 10.1007/s00497-005-0006-0
Pawełkowicz, M., Pryszcz, L., Skarzyńska, A., Wóycicki, R., Posyniak, K., Rymuszka, J., et al. (2019). Comparative transcriptome analysis reveals new molecular pathways for cucumber genes related to sex determination. Plant Reprod. 32, 193–216. doi: 10.1007/s00497-019-00362-z
Pimenta Lange, M. J., and Lange, T. (2016). Ovary-derived precursor gibberellin A9 is essential for female flower development in cucumber. Development 143, 4425–4429. doi: 10.1242/dev.135947
Rudich, J., and Halevy, A. H. (1974). Involvement of abscisic acid in the regulation of sex expression in the cucumber. Plant Cell Physiol. 15, 635–642. doi: 10.1093/oxfordjournals.pcp.a075048
Růžička, K., Ljung, K., Vanneste, S., Podhorská, R., Beeckman, T., Friml, J., et al. (2007). Ethylene regulates root growth through effects on auxin biosynthesis and transport-dependent auxin distribution. Plant Cell 19, 2197–2212. doi: 10.1105/tpc.107.052126
Saito, S., Fujii, N., Miyazawa, Y., Yamasaki, S., Matsuura, S., Mizusawa, H., et al. (2007). Correlation between development of female flower buds and expression of the CS-ACS2 gene in cucumber plants. J. Exp. Bot. 58, 2897–2907. doi: 10.1093/jxb/erm141
Song, J., Gao, Z. H., Huo, X. M., Sun, H. L., Xu, Y. S., Shi, T., et al. (2015). Genome-wide identification of the auxin response factor (ARF) gene family and expression analysis of its role associated with pistil development in Japanese apricot (Prunus mume Sieb. Et Zucc). Acta Physiol. Plant. 37:145. doi: 10.1007/s11738-015-1882-z
Sun, Z. Y., Zhang, T. J., Su, J. Q., Soon Chow, W., Liu, J. Q., Chen, L. L., et al. (2015). A novel role of ethephon in controlling the noxious weed Ipomoea cairica (Linn.) sweet. Sci. Rep. 5:11372. doi: 10.1038/srep11372
Sun, H. H., Wu, S., Zhang, G. Y., Jiao, C., Guo, S. G., Ren, Y., et al. (2017). Karyotype stability and unbiased fractionation in the paleo-allotetraploid Cucurbita genomes. Molecular Plant. 10, 1293–1306. doi: 10.1016/j.molp.2017.09.003
Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578. doi: 10.1038/nprot.2012.016
Trebitsh, T., Rudich, J., and Riov, J. (1987). Auxin, biosynthesis of ethylene and sex expression in cucumber (Cucumis sativus). Plant Growth Regul. 5, 105–113. doi: 10.1007/BF00024738
Vandenbussche, F., Vaseva, I., Vissenberg, K., and Van Der Straeten, D. (2012). Ethylene in vegetative development: a tale with a riddle. New Phytol. 194, 895–909. doi: 10.1111/j.1469-8137.2012.04100.x
Vanstraelen, M., and Benková, E. (2012). Hormonal interactions in the regulation of plant development. Annu. Rev. Cell Dev. Biol. 28, 463–487. doi: 10.1146/annurev-cellbio-101011-155741
Wang, D. H., Li, F., Duan, Q. H., Han, T., Xu, Z. H., and Bai, S. N. (2010). Ethylene perception is involved in female cucumber flower development. Plant J. 61, 862–872. doi: 10.1111/j.1365-313X.2009.04114.x
Yamasaki, S., Fujii, N., Matsuura, S., Mizusawa, H., and Takahashi, H. (2001). The M locus and ethylene-controlled sex determination in andromonoecious cucumber plants. Plant Cell Physiol. 42, 608–619. doi: 10.1093/pcp/pce076
Yamasaki, S., Fujii, N., and Takahashi, H. (2000). The ethylene-regulated expression of CS-ETR2 and CS-ERS genes in cucumber plants and their possible involvement with sex expression in flowers. Plant Cell Physiol. 41, 608–616. doi: 10.1093/pcp/41.5.608
Yamasaki, S., Fujii, N., and Takahashi, H. (2003). Characterization of ethylene effects on sex determination in cucumber plants. Sex. Plant Reprod. 16, 103–111. doi: 10.1007/s00497-003-0183-7
Yang, X. X., Qu, S. P., Yang, G. X., Wang, Y. Y., Xu, W. L., and Cui, C. S. (2015). Effect of ethephon, ethylene inhibitor on Cucurbita maxima sex expression. Northern Hortic. 18, 30–34. doi: 10.11937/bfyy.201518008
Young, M. D., Wakefield, M. J., Smyth, G. K., and Oshlack, A. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11:R14. doi: 10.1186/gb-2010-11-2-r14
Zhang, Y., Zhao, G. Y., Li, Y. S., Mo, N., Zhang, J., and Liang, Y. (2017). Transcriptomic analysis implies that GA regulates sex expression via ethylene-dependent and ethylene-independent pathways in cucumber (Cucumis sativus L.). Front. Plant Sci. 8:10. doi: 10.3389/fpls.2017.00010
Keywords: endogenous hormone, ethephon, pumpkin, sex differentiation, yield
Citation: Li Q, Guo W, Chen B, Pan F, Yang H, Zhou J, Wang G and Li X (2021) Transcriptional and Hormonal Responses in Ethephon-Induced Promotion of Femaleness in Pumpkin. Front. Plant Sci. 12:715487. doi: 10.3389/fpls.2021.715487
Edited by:
Wusirika Ramakrishna, Central University of Punjab, IndiaReviewed by:
Rachana D. Bhardwaj, Punjab Agricultural University, IndiaMarshall Keyster, University of the Western Cape, South Africa
Muhammad Azher Nawaz, University of Sargodha, Pakistan
Copyright © 2021 Li, Guo, Chen, Pan, Yang, Zhou, Wang and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xinzheng Li, lxz2283@126.com