SWATH-MS based quantitive proteomics reveal regulatory metabolism and networks of androdioecy breeding system in Osmanthus fragrans
BMC Plant Biology volume 21, Article number: 468 (2021)
The fragrant flower plant Osmanthus fragrans has an extremely rare androdioecious breeding system displaying the occurrence of males and hermaphrodites in a single population, which occupies a crucial intermediate stage in the evolutionary transition between hermaphroditism and dioecy. However, the molecular mechanism of androdioecy plant is very limited and still largely unknown.
Here, we used SWATH-MS-based quantitative approach to study the proteome changes between male and hermaphroditic O. fragrans pistils. A total of 428 proteins of diverse functions were determined to show significant abundance changes including 210 up-regulated and 218 down-regulated proteins in male compared to hermaphroditic pistils. Functional categorization revealed that the differentially expressed proteins (DEPs) primarily distributed in the carbohydrate metabolism, secondary metabolism as well as signaling cascades. Further experimental analysis showed the substantial carbohydrates accumulation associated with promoted net photosynthetic rate and water use efficiency were observed in purplish red pedicel of hermaphroditic flower compared with green pedicel of male flower, implicating glucose metabolism serves as nutritional modulator for the differentiation of male and hermaphroditic flower. Meanwhile, the entire upregulation of secondary metabolism including flavonoids, isoprenoids and lignins seem to protect and maintain the male function in male flowers, well explaining important feature of androdioecy that aborted pistil of a male flower still has a male function. Furthermore, nine selected DEPs were validated via gene expression analysis, suggesting an extra layer of post-transcriptional regulation occurs during O. fragrans floral development.
Taken together, our findings represent the first SWATH-MS-based proteomic report in androdioecy plant O. fragrans, which reveal carbohydrate metabolism, secondary metabolism and post-transcriptional regulation contributing to the androdioecy breeding system and ultimately extend our understanding on genetic basis as well as the industrialization development of O. fragrans.
Androdioecy is an exceedingly rare mating system in which males and hermaphrodites co-occur in a single population with normal sexual function . Androdioecious plants are extremely rare in nature, less than 50 plants have been reported worldwide such as Datisca glomerata , Mercurialis annua , Pseudoxandra spiritus-sancti  and Tapiscia sinensis [5, 6], drawing widespread interests by scientists. Surprisingly, such an unusual breeding system prefers to be distributed in the Oleaceae species including Phillyrea latifolia , Phillyrea angustifolia , Fraxinus lanuginose , Chionanthus retusus , Osmanthus serrulatus , O. delavayi  etc., accompanied by higher frequencies of males compared to other androdioecious species . Osmanthus fragrans Lour. belongs to Osmanthus of Oleaceae, with important dietary value and economic value, widely used in food processing industry through their fragrant flowers. O. fragrans has two sex traits, the male plant has an abortive pistil and a normally developed stamen, the other is the hermaphroditic plant with normal development of both stamen and pistil, which is a proper model to investigate androdioecy . Nevertheless, the molecular mechanism of reproductive development and sex differentiation in androdioecy plants is still largely unknown and requires more detailed investigations.
Since sex differentiation has recently become a research hotspot, proteomics has been demonstrated as a powerful tool to examine sex-related differences in plants and even the stress adaption changes from dioecious plants . For example, comparative proteomic analysis of the mitochondrial proteome between sterile rice and fertile rice revealed a clade of regulatory proteins with decreased abundance in mitochondrial complex leading to cytoplasmic male sterility in sterile rice . A gel-based proteomics study on male and female leaves of jojoba (Simmondsia chinensis) established 45 molecular protein markers for early sex differentiation . Similarity, proteomics analysis of male and female pumpkin nectar identified a total of 12 specific proteins unique to female and male nectar . The total proteome profiling of flower buds of different genders was analyzed using label-free proteomics methods, providing a basis for identifying the key proteins in the development of male and female flowers of Coccinia grandis . Interestingly, proteome changes of male and female Populus cathayana under drought stress were investigated by two-dimensional electrophoresis (2-DE) method, the upregulation of photosynthesis, energy metabolism and stress responses in male poplars explaining the observation that male poplars are more resistant to drought stress compared to female poplars . Such phenomenon also occurred in other abiotic stresses such as salt, chilling and heavy metals [21,22,23,24]. Moreover, large-scale proteomic analysis has been widely used in many other plant research as a promising tool for protein identification and gene function characterization [25, 26].
The innovative SWATH-MS (Sequential Windowed Acquisition of All Theoretical Mass Spectra) approach employs a high specificity data-independent acquisition (DIA) method coupled with a novel targeted data extraction methodology, displays the advantages of higher reproducibility, quantitative consistency and accuracy compared with gel-based and label-based or defined DDA proteomic approaches [27,28,29], particularly for detecting and quantifying low abundant proteins, which has always been the major challenges in the proteomics investigation [30, 31]. In this study, a total of 2298 proteins were quantitatively identified by SWATH-MS in O. fragrans flowers, 428 of which were found to be differentially expressed between male and hermaphroditic pistils. Alternations were observed in different metabolic processes such as glucose metabolism, flavonoids and isoprenoids metabolism, related signaling pathway, which contributes to the normal developmental growth of anthers and pollen tube. Further experimental analysis revealed that hermaphroditic O. fragrans could accumulate more carbohydrates than males during the full flowering stage. Our work represents the first SWATH-MS investigation in woody plant O. fragrans, hoping to provide a comprehensive proteome reference for the unique androdioecy breeding system from the perspective of proteomics.
Phenotype characterization for appropriate sampling and following SWATH-MS analysis
To perform a comprehensive proteomes analyses on androdioecy breeding system in Osmanthus fragrans, we used innovative SWATH-MS-based quantitative proteomic strategy to compare the proteomes of male and hermaphroditic flowers as outlined in Fig. 1, which employs a high specificity data-independent acquisition (DIA) method coupled with a novel targeted data extraction methology. The acquired DIA file recording all the complete chromatographic elution traces of peak groups was introduced into a non target analysis by the MS2-based quantification methods, which provides continuous and sufficient information in protein quantification. Hence, the quantitative analyses of peptides were supported by given extracted ion chromatograms from both MS1 and MS2 level, enabling comprehensive proteome profiling and detecting low-abundance proteins in this study.
Both male and hermaphroditic flowers of O. fragrans were sampled at the full bloom stage, and other floral organs were removed manually, leaving only the pistils. The hermaphroditic flower pistil has three complete organs: papillary stigma, shorter style, and enlarged ovary. Nevertheless, the pistil of male flowers was aborted, two carpels were separated without these three parts, accounting for pistil abortion. Simultaneously, the pistils of male flowers were significantly shorter than that of hermaphroditic flowers, 0.75 mm and 1.21 mm, respectively (Fig. 2a, b). Therefore, the pistils of male and hermaphroditic flowers were used in this proteome study by the SWATH-MS method. The correlation coefficient of the samples within the group (M and H) was generally higher than 0.95 whereas the correlation between the groups was relatively lower at around 0.6 (Fig. 2c), strongly reflected the high reliability of the appropriate sampling and suitable for subsequent analysis. In this study, a total of 2298 proteins were identified in male and hermaphroditic flower pistils (Table S1), among which 428 proteins showing significant changes in protein abundance including 210 increased and 218 decreased in the male compared to hermaphroditic pistils as shown in the volcano map (Fig. 2d).
Comprehensive inventory of proteome changes between male and hermaphroditic flowers
The functional classifications of these DEPs were carried out through Gene Ontology (GO) analysis including biological processes (BP), cellular components (CC) and molecular functions (MF). The top terms corresponding to the up-regulated DEPs of BP/CC/MF were response to cadmium ion, thylakoid, and coenzyme binding whereas the top terms corresponding to the down-regulated DEPs of BP/CC/MF were response to metal ion, plant-type cell wall and hydrolase activity, hydrolyzing O-glycosyl compounds (Fig. 3). In BP terms for up-regulated and down-regulated proteins, “response to cadmium ion” and “response to metal ion” constitute the largest number of proteins in both cases, as revealed by previous studies that O. fragrans has resistance to metal ions such as cadmium, lead and copper [32, 33]. However, there were few studies on the molecular response and resistance mechanisms of O. fragrans to heavy metals still requiring further investigations. It is also worth noting that massive down-regulated proteins are assigned with the “hydrolase activity” and “hydrolyzing O-glycosyl compounds” MF term in the male flower, implicating their potential roles for normal pistils development.
All the DEPs were mapped and enriched in top 13 pathways such as carbon metabolism, amino sugar metabolism and starch metabolism by KEGG pathway classification (Fig. 4a), which provide a quick view of the most significant pathway in the female sterility. Meanwhile, MapMAN BIN analysis was also conducted on all the DEPs (Table S2), revealing several metabolic processes including protein metabolism, RNA processing, transport, and secondary metabolism were greatly influenced during the pistil abortion of male O. fragrans flowers (Fig. 4b). In particularly, 12 of the 13 DEPs involving secondary metabolism were up-regulated while most of signaling related proteins were down-regulated proteins in the male pistils. Such completely opposite protein expression patterns in these two categories between male and hermaphroditic pistils, which may be the key response attribute to the pistil abortion of male flowers. In the following sections, those activated or repressed enzymes will be mapped to different physiological and biochemical pathways, and the relevance of diverse biological processes and molecular androdioecy characteristics of O. fragrans is discussed.
Protein-protein interaction (PPI) network analysis of DEPs
The study of protein-protein interaction (PPI) is beneficial to explore the core proteins and better understand the protein regulatory network. We used 428 DEPs for protein-protein interaction network analysis, and the two main protein interaction modules were screened by MCODE (Fig. 5). In the first module, there were 17 nodes and 117 edges, of which the key protein was ofr.gene22654, an enrichment factor protein (Fig. 5a). Its knockout mutant can significantly reduce seed setting rate in Arabidopsis, only about 40% of the wild type . In the second module, there were a total of 8 nodes and 23 edges. Among all 8 proteins, the protein ofr.gene26044 is the only down-regulated protein that plays a role in the synthesis of Jasmonic acid (Fig. 5b). As a signaling molecule, jasmonic acid is deeply involved in the growth and development of plants. Many studies have shown that the decrease of JA content caused by JA synthesis related gene mutation can lead to female infertility in tomato (Solanum lycopersicum) [35,36,37].
Correlation analysis between differentially expressed gene and protein in male and hermaphroditic pistils
To understand whether the expression of differentially expressed proteins are associated with transcriptional changes, we selected 9 genes involving a variety of processes for expression levels analysis by qRT-PCR. As shown in Fig. 6, four of the selected genes showed up- or down- regulation in accordance with the changes in abundance of the corresponding proteins identified from SWATH-MS experiment (Table S1), which are distributed in carbohydrate metabolism, secondary metabolism, signaling pathway as well as cell wall synthesis. However, two genes including ofr.gene28906 and ofr.gene18814 were expressed at the opposite level of protein expression (Fig. 6). Such difference between transcription levels and protein expression levels may be resulting from post-transcriptional regulation such as alternative splicing (AS), mRNA stability, and mRNA translation . Additionally, ofr.gene43330, ofr.gene8745 and ofr.gene41183 were all up-regulated in the proteomic data, but their gene expression was slightly up-regulated or no significant difference in the male compared to hermaphroditic pistil (Fig. 6), implicating a complex and delicate regulation stage from the transcriptome to the proteome.
Carbohydrate metabolism and TCA cycle in pistils
Fertilization of higher plants depend on the pollen tube elongation of the style toward the ovary, which was also considered to be the fastest growing plant cell requiring a huge consumption of nutrient and energy [39, 40]. Since the DEPs primarily focused on the carbohydrates and amino acids metabolism pathways, such gap probably restricted the normal developmental growth of pollen tube in male flowers, resulting in the disruption of pollen germination at stigma and complete fertilization. Consistently, down-regulation of major carbon and sugar metabolism responsible for the nutrient and energy supply were observed in male plants with incomplete pistil (Table S2).
In the starch and sucrose metabolism, UDP-glucose eventually converted to D-glucose, 4 up-regulated proteins and 7 down-regulated proteins involving this pathway were identified in male flower compared with hermaphroditic flower pistils (Fig. 7a). The up-regulated proteins promoted the conversion of UDP-glucose into other forms of carbohydrate and the down-regulated proteins inhibited the synthesis of D-glucose. Therefore, the UDP-glucose and glucose in male flower were both reduced in comparison to hermaphroditic flowers. As UDP-glucose functions as metabolic precursors of cellulose for the inner layer of pollen tube cell wall formation , the reduced contents of glucose may be the main cause for the abortive pistil in male O. fragrans flowers without complete pistil structure such as stigma, style and ovary.
The TCA cycle is a crucial pathway connection for sugar, lipid, and protein metabolism. KEGG pathway analysis and MapMAN BIN classification simultaneously revealed significant changes in the TCA cycle with 6 identical up-regulated proteins in male pistils (Fig. 4b, Fig. 7a), which well mapped to each reaction step involving the TCA cycle within the generation of different energy related products. In fact, the promoted TCA cycle accelerated the consumption of glucose, further restricting the development of pollen tube cell walls, thus affecting pistil growth in male flower through the deficiency of nutrient and energy supply.
Moreover, the full bloom stage is the main stage of flowering and pollination. In this process, the pedicel of hermaphroditic flower gradually turns purplish red whereas the male remains green (Fig. 7b), strongly implicating the accumulation of carbohydrate-related compounds in hermaphroditic flower. To further detect the carbohydrate metabolism distribution, LI-6400XT portable photosynthetic measurement system was performed to determine photosynthesis of male and hermaphroditic O. fragrans during the full flowering stage (Table S4). The results showed that the net photosynthetic rate and water use efficiency in hermaphroditic flower were higher than that of males (Fig. 7c), indicating that hermaphroditic O. fragrans could accumulate more carbohydrates with less water consumption to fulfill pollen germination, pollen tube growth and following fruit development in accordance with the results of the above proteome data.
Up-regulation of secondary metabolism for promoting male function in male flowers
By MapMAN BIN and KEGG categorizations, a large proportion of proteins involving the secondary metabolism was shown to be activated in male flowers, particularly for the biosynthesis of secondary metabolites such as flavonoids, isoprenoids and phenylpropanoids, etc. (Fig. 8). Generally, the precursors of flavonoids and isoprenoids were both p-coumaroyl CoA in the phenylpropanoids pathway, which was synthesized by chalcone synthase (CHS) and hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase (HCT), respectively . The ofr.gene25988 and ofr.gene25148 are the key enzymes for the synthesis of flavonoids and their homologs namely as less adhesive pollen 5 (LAP5) and LAP6, respectively, were essential for the formation of the pollen exine in Arabidopsis . Meanwhile, the acetoacetyl-CoA thiolase (AACT) ofr.gene40770 catalyzes the conversion of acetyl-CoA to acetoacetyl-CoA for isoprenoids metabolic pathway. Mutation of its homolog in Arabidopsis could lose the pollen coat and tapetal cells leading to male sterility . Accordingly, the flavonoids and other substances are important components of the anther or pollen tapetum to ensure the fertility of males, playing an important role in the attachment and recognition of pollination . The absence of flavonoids severely influence pollen germination and pollen tube growth without any damage to the size and structure of pollen resulting in male sterility . DEPs participating in the flavonoids and isoprenoids biosynthesis pathway were up-regulated by 9.3-fold on average in the pistil of male flowers for ensuring the fertility of pollen (Fig. 8), feasible to explain important features of androdioecy that even the pistil of the male flower has been aborted, whose pollen still remains vigorous and can pollinate with the fertile pistil of other hermaphroditic flowers. And previous studies believe that the fitness of male function in the male plant was at least twice as high as that of the hermaphroditic plant in the androdioecy breeding system in order to maintain its existence in nature [47, 48]. Therefore, those important secondary metabolites including flavonoids, isoprenoids and lignin are likely to be protective a role for maintaining male function in male flowers.
Down-regulation of signaling related proteins in pistils of male flowers
In contrast to secondary metabolism, DEPs in the signaling categorization showed reduced protein abundances changes, most of which were highly correlated with calcium signaling, receptor-like kinase (RLK) and leucine rich repeat receptor-like kinase (LRR-RLK) (Table 1; Table S2). It was well known that calcium is an important second messenger involving various signal transduction pathways, particular for stress response and reproductive growth in plants . Free calcium is essential for pollen germination, pollen tube growth, and pollen tube guidance . The dynamic changes of calcium content were investigated during the pistil development of olive (Olea europaea L.), revealing calcium gradually accumulated on the stigma and further expanded toward the style and ovary after the anthers dehiscence, implicating the growth direction of the pollen tube was determined by the concentration of calcium . Meanwhile, distribution changes of calcium regulate programmed cell death to affect the differentiation of male and female flowers in litchi (Litchi chinensis Sonn.) such as the degeneration of style cells by the absence of calcium precipitation . Therefore, down regulation of calcium signaling (ofr.gene26367, ofr.gene11138 and ofr.gene18811) probably attributed to the degeneration of pistils in male flowers. On the other hand, the CrRLK gene family has been widely demonstrated to be associated with the tip growth and cell wall integrity of pollen tubes during fertilization [53, 54]. For instance, the Arabidopsis knockout mutant feronia (FRE) belonging to CrRLK family failed to release sperm cells and impair pollen tube reception . The Buddha’s Paper Seal 1 (BUPS1) and BUPS2 were also required for the pollen tube development and disruption of them caused the swollen of pollen tube unable to complete fertilization . Recently, 26 CrRLK genes in the pear genome were identified, among which PbrCrRLK1L26 regulate pollen tube elongation and PbrCrRLK1L3 could control pollen tube rupture . Accordingly, we identified two CrRLK (ofr.gene38451 and ofr.gene2855) may modulate pollen tubes development responsible for the aborted pistils of male flowers in O. fragrans. Taken together, entire down-regulation of both calcium signaling and RLK related proteins probably severely restrict the pistils morphogenesis including pollen tube elongation, sperm cells guidance after rupture as well as carpel fusion, leading to female sterility in male flowers.
Other candidate proteins contribute to androdioecy breeding system in O. fragrans
Seeking target proteins for the differentiation of male and female flowers in O. fragrans is necessary to illuminate the molecular mechanism of androdioecy breeding system. Here, we collected 25 candidate proteins with dramatic abundance changes (Ratio > 30 or Ratio < 0.03) between M and H as shown in Table 2. Representatively, UDP-sugar pyrophosphorylase (ofr.gene8745, USP) was the largest up to 190 folds of up-regulation in the male flowers, which has been reported to be essential for the recycling of xylose and arabinose. Mutation of USP exhibited the incomplete cell wall of anther and pollen leading to sterility in Arabidopsis . Furthermore, ofr.gene58609 and ofr.gene10710 participates in the biosynthesis of lignin, which is critical for cell wall remodeling and modification in pollen development [59, 60]. Therefore, these proteins showing very high abundance in the male pistil are probably guaranteed to maintain male fertility, and improve the competitiveness and compatibility of male pollen. On the other hand, the lowest expression of candidate protein (ofr.gene35198 and ofr.gene10041) in the pistils of male flower compared to the hermaphroditic flower encodes pectin lyase-like and pectin acetylesterase activity, emerging evidences have demonstrated that pectin is the main component of the apical cell wall, which contributes to the mechanical strength of the cell wall and plays an important role in the growth of pollen tubes which may be associated with the pistil abortion of O. fragrans [61, 62].
In the androdioecy breeding system of O. fragrans, the pistils of male flowers were aborted without fully developed pistils but the pollen remains vigorous and can be pollinated whereas the pistil of the hermaphroditic flower developed normally. The present work provides a comprehensive understanding on the molecular androdioecy characteristics in O. fragrans by the advanced label-free SWATH-MS quantitative proteomics platform, which revealed significant proteome changes in carbohydrate metabolism, secondary metabolism as well as calcium signaling between M and H pistil. Accordingly, the reduced glucose metabolism cannot support normal pollen tube development leading to the aborted pistils in male flowers. On the contrary, sufficient carbohydrates accumulation in purplish red pedicel of hermaphroditic flower associated with experimental evidence of promoted net photosynthetic rate and water use efficiency in comparison to male flower (Fig. 7) further suggest glucose serves as nutritional modulator for the differentiation of male and hermaphroditic flower. Moreover, the entire upregulation of secondary metabolism including flavonoids, isoprenoids and lignin seem to protect and maintain the male function in male flowers, well explaining important feature of androdioecy that aborted pistil of a male flower still has a male function. In addition, down-regulation of calcium signaling and RLK related proteins also inhibit the pistil morphogenesis resulting in female sterility in male plants. Taken together, our work represents the first SWATH-MS-based proteomic study in androdioecy plant O. fragrans, which would provide new clues for further studies on the sex differentiation of the androdioecious O. fragrans and will extend our understanding on androdioecy breeding system.
Sample collection and preparation
Male and hermaphroditic flowers of O. fragrans were collected at their full flowering stage from the campus of Nanjing Forestry University. The development process of male and hermaphroditic O. fragrans is almost identical, and the main difference lies in the structure of pistil . Thus, other tissues such as petals, pedicels, anthers were removed manually, only pistils of male and hermaphroditic flowers were left as experimental materials to eliminate the influence of other floral tissues as much as possible. The pistil length was measured with electronic vernier caliper, and each was repeated 6 times. The pistil samples of each sex for proteome detection had three repetitions, approximately 1 g per replicate. The treated pistils were fully ground in liquid nitrogen, added with 2.5% SDS/100 mm Tris-HCl lysis buffer. After 15 min of ultrasonic treatment on ice water, the supernatants were centrifuged at 16000 g for 20 min. Acetone was added to the supernatants to precipitate the protein. After cleaning with acetone and drying in the air, 8 M Urea / 100 mm Tris-HCl solutions were added to the protein precipitation to fully dissolve the proteins. After centrifuged at 12000 g for 15 min, the supernatants were added with dithiothreitol (DTT) to the final concentration of 10 mM, and incubated at 37 °C for 1 h. Iodoacetamide (IAM) was added until the final concentration of 40 mM, and the alkylation was carried out at room temperature in a dark state to seal the sulfhydryl group. Then, 100 mM Tris-HCl solutions were added, determining protein concentration by the Bradford method, and diluted urea concentration to less than 2 M. Trypsin (50:1 protein to trypsin) was added and oscillated at 37 °C overnight. The PH value of the solution was adjusted at about 6.0, then centrifuged at 12000 g for 15 min, and desalted with C18 columns. The desalted peptide solution was dried by a centrifugal concentrator and frozen at - 20 °C for mass spectrometry detection.
Peptides samples were detected by the Triple TOF 5600 (Sciex) LC/MS system. The prepared samples were bound to the trap column and then separated by the analytical column (45 min gradient, 60 min total time). Two mobile phases were established to analyze the gradient: Buffer A-0.1% (V/V) formic acid, 5% DMSO in H2O, Buffer B-0.1% (V/V) formic acid, and 5% DMSO in acetonitrile. For SWATH scanning, one MS1 scan (ion accumulation time 250 ms, scanning range 350–1500 m/z) and 100 MS2 scans with variable windows (ion accumulation time 33 ms, scanning range 100–1800 m/z) were included in each cycle. The mass spectrum files obtained by SWATH scanning were processed by DIA-Umpire to obtain the secondary mass spectrum file that can be used for database search. TPP software was used for database retrieval, and the retrieved results were used as a spectral library, and OpenSWATH algorithm was used for SWATH targeting extraction, and a false discovery rate (FDR) of <1% was set as selection criteria . The protein quantitative intensity information obtained by SWATH analysis was applied to log2 conversion, data filling, and data normalization using the imputation algorithm in Perseus software for difference comparison and T-test analysis. Proteins of male and hermaphroditic pistils with a ratio of above 5 or below 0.2 (P < 0.05) were considered as differentially expressed proteins (DEPs) in this study.
Proteomic data analysis
The correlation of protein quantification was analyzed by corrplot R package to evaluate the reliability of SWATH quantitative proteome data. The DEPs identified were used for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis by clusterProfiler R package , and pathways of KEGG enrichment analysis results is drawn with reference to KEGG mapper . After the DEPs were compared with the homologous proteins in Arabidopsis thaliana, the MapMAN BIN system was used for functional classification (http://ppdb.tc.cornell.edu/dbsearch/searchacc.aspx). Heatmaps of DEPs between M and H were drawn by pheatmap R package . Protein-protein interaction (PPI) network analysis of DEPs was carried out using String (https://string-db.org/), and the clustering function modules were screened by the Molecular Complex Detection (MCODE) plugin and visualized by Cytoscape [67, 68].
Measurement of photosynthesis and quantitative real-time PCR analysis
The LI-6400XT portable photosynthesis measurement system was used to determine the photosynthesis of male and hermaphroditic O. fragrans at the full flowering stage according to the instructions. Four plants of each sex were selected, and 8 healthy and mature leaves of each plant were selected for measurement, and the net photosynthetic rate (Pn), water use efficiency (WUE) and other data were recorded.
Total RNA of male and hermaphroditic pistils were extracted with RNAprep Pure Plant Plus Kit (TIANGEN Biotech, Beijing) according to the manufacturer’s instructions. Then, 5 μg RNA was reversed transcribed by Evo M-MLV RT Premix for qPCR (Accurate Biotechnology, Hunan) for cDNA synthesis. SYBR Green Premix Pro Taq HS qPCR Kit (Accurate Biotechnology, Hunan) was used for quantitative real-time PCR (qRT-PCR) experiments with OfACT as a reference . qRT-PCR was performed using ABI StepOnePlus Systems (Applied Biosystems, USA), the reaction steps were as follows: 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s and 60 °C for 30 s. All primers used in this experiment were listed in Table S3.
Availability of data and materials
The data sets are included within the article and its additional files. The plant materials are available from the corresponding author on reasonable request.
Sequential Windowed Acquisition of All Theoretical Mass Spectra
Differentially expressed proteins
Kyoto Encyclopedia of Genes and Genomes
- TCA cycle:
tricarboxylic acid cycle
Hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase
Less adhesive pollen
Charlesworth D. Androdioecy and evolution of dioecy. Biol J Linn Soc. 1984;22:333–48.
Aaron L, Loren HR, Thomas SE. Functional androdioecy in the flowering plant Datisca glomerata. Nature. 1990;343:641–2.
John P. Widespread functional androdioecy in Mercurialis annua L. (Euphorbiaceae). Biol J Linn Soc. 1997;61:95–116.
Lopes JC, Chatrou LW, Mello-Silva R, Rudall PJ, Sajo MG. Phylogenomics and evolution of floral traits in the Neotropical tribe Malmeeae (Annonaceae). Mol Phylogenet Evol. 2018;118:379–91.
Yang K, Zhou XJ, Wang YY, Feng HL, Ren XL, Liu HD, et al. Carbohydrate metabolism and gene regulation during anther development in an androdioecious tree, Tapiscia sinensis. Ann Bot. 2017;120:967–77.
Xin GL, Liu JQ, Liu J, Ren XL, Du XM, Liu WZ. Anatomy and RNA-Seq reveal important gene pathways regulating sex differentiation in a functionally Androdioecious tree, Tapiscia sinensis. BMC Plant Biol. 2019;19:554.
Aronne G, Wilcock CC. Reproductive characteristics and breeding system of shrubs of the Mediterranean region. Funct Ecol. 1994;8:69–76.
Jacques L, Bertrand D. Is Phillyrea angustifolia L. (Oleaceae) an androdioecious species? Biol J Linn Soc. 1992;108:375–87.
Ishida K, Hiura T. Pollen fertility and flowering phenology in an Androdioecious tree, Fraxinus lanuginosa (Oleaceae), in Hokkaido, Japan. Int J Plant Sci. 1998;159:941–7.
Song JH, Oak MK, Hong SP. Morphological traits in an androdioecious species, Chionanthus retusus (Oleaceae). Flora. 2016;223:129–37.
Yang GD, Qian HR, Chen L, Wang XR. Reproduction system of Osmanthus serrulatus, an endemic plant to China. Sci Silva Sin. 2018;54:17–29 (in Chinese).
Duan YF, Li WH, Zheng SY, Sylveste SP, Li YF, Cai FY, et al. Functional androdioecy in the ornamental shrub Osmanthus delavayi (Oleaceae). PLoS One. 2019;14:e0221898.
Saumitou-Laprade P, Vernet P, Vassiliadis C, Hoareau Y, de Magny G, Dommée B, et al. A self-incompatibility system explains high male frequencies in an androdioecious plant. Science. 2010;327:1648–50.
Xu YC, Zhou LH, Hu SQ, Hao RM, Huang CJ, Zhao HB. The differentiation and development of pistils of hermaphrodites and pistillodes of males in androdioecious Osmanthus fragrans L. and implications for the evolution to androdioecy. Plant Syst Evol. 2014;300:843–9.
Yang L, Gong FP, Xiong EH, Wang W. Proteomics: a promising tool for research on sex-related differences in dioecious plants. Front Plant Sci. 2015;6:954.
Liu G, Tian H, Huang YQ, Hu J, Ji YX, Li SQ, et al. Alterations of mitochondrialprotein assembly and jasmonic acid biosynthesis pathway in honglian (HL)-type cytoplasmic male sterility rice. J Biol Chem. 2012;287:40051–60.
Al-Obaidi JR, Rahmad N, Hanafi NM, Halabi MF, Al-Soqeer AA. Comparative proteomic analysis of male and female plants in Jojoba (Simmondsia chinensis) leaves revealed changes in proteins involved in photosynthesis, metabolism, energy, and biotic and abiotic stresses. Acta Physiol Plant. 2017;39:179.
Chatt EC, Aderkas PV, Carter CJ, Smith D, Elliott M, Nikolau BJ. Sex-dependent variation of pumpkin (Cucurbita maxima cv. Big Max) nectar and nectaries as determined by proteomics and metabolomics. Front Plant Sci. 2018;9:860.
Devani RS, Chirmade T, Sinha S, Bendahmane A, Dholakia BB, Banerjee AK, et al. Flower bud proteome reveals modulation of sex-biased proteins potentially associated with sex expression and modification in dioecious Coccinia grandis. BMC Plant Biol. 2019;19:330.
Zhang S, Chen FG, Peng SM, Ma WJ, Korpelainen H, Li CY. Comparative physiological, ultrastructural and proteomic analyses reveal sexual differences in the responses of Populus cathayana under drought stress. Proteomics. 2010;10:2661–77.
Chen F, Zhang S, Jiang H, Korpelainen H, Li CY. Comparative proteomics analysis of salt response reveals sex-related photosynthetic inhibition by salinity in Populus cathayana cuttings. J Proteome Res. 2011;10:3944–58.
Chen F, Zhang S, Zhu GP, Korpelainen H, Li CY. Populus cathayana males are less affected than females by excess manganese: Comparative proteomic and physiological analyses. Proteomics. 2013;13:2424–37.
Zhang S, Feng LH, Jiang H, Ma WJ, Korpelainen H, Li CY. Biochemical and proteomic analyses reveal that Populus cathayana males and females have different metabolic activities under chilling stress. J Proteome Res. 2012;11:5815–26.
Zhang S, Zhang YX, Cao YC, Lei YB, Jiang H. Quantitative proteomic analysis reveals Populus cathayana females are more sensitive and respond more sophisticatedly to iron deficiency than males. J Proteome Res. 2016;15:840–50.
Henriet C, Balliau T, Aimé D, Le Signor C, Kreplak J, Zivy M, et al. Proteomics of developing pea seeds reveals a complex antioxidant network underlying the response to sulfur deficiency and water stress. J Exp Bot. 2021;72:2611–26.
Wang G, Wang G, Wang J, Du Y, Yao D, Shuai B, et al. Comprehensive proteomic analysis of developing protein bodies in maize (Zea mays) endosperm provides novel insights into its biogenesis. J Exp Bot. 2016;67:6323–35.
Bourassa S, Fournier F, Nehmé B, Kelly I, Tremblay A, Lemelin V, et al. Evaluation of iTRAQ and SWATH-MS for the quantification of proteins associated with insulin resistance in human duodenal biopsy samples. PLoS One. 2015;10:e0125934.
Ludwig C, Gillet L, Rosenberger G, Amon S, Collins BC, Aebersold R. Data-independent acquisition-based SWATH-MS for quantitative proteomics: a tutorial. Mol Syst Biol. 2018;14:e8126.
Zhu XC, Chen YP, Subramanian R. Comparison of information-dependent acquisition, SWATH, and MSAll techniques in metabolite identification study employing ultrahigh-performance liquid chromatography-quadrupole time-of-flight mass spectrometry. Anal Chem. 2014;86:1202–9.
Chen MX, Zhang YJ, Fernie AR, Liu YG, Zhu FY. SWATH-MS-based proteomics: strategies and applications in plants. Trends Biotechnol. 2020;39:433–7.
Zhu FY, Song YC, Zhang KL, Chen X, Chen MX. Quantifying plant dynamic proteomes by SWATH-based mass spectrometry. Trends Plant Sci. 2020;25:1171–2.
Wu FZ, Yang WQ, Zhang J, Zhou LQ. Effects of cadmium stress on growth and nutrient accumulation, distribution and utilization in Osmanthus fragrans var. thunbergii. Chin J Plant Ecol. 2010;34:1220–6 (in Chinese).
Qian YZ, Zhang N, Chen H, Du JL, Meng LY. Comparative study of uptake, translocation and accumulation characters in 4 varieties of Osmanthus fragrans (Thunb.) Lour. under single stress of Cd, Pb or Cu heavy mental treatment stress. J Anhui Univ (Natural Science Edition). 2018;42:102–8 (in Chinese).
Albrecht V, Ingenfeld A, Apel K. Characterization of the snowy cotyledon 1 mutant of Arabidopsis thaliana: the impact of chloroplast elongation factor G on chloroplast development and plant vitality. Plant Mol Biol. 2006;60:507–18.
Li L, Zhao YF, Mccaig BC, Wingerd BA, Wang JH, Whalon ME, et al. The tomato homolog of CORONATINE-INSENSITIVE1 is required for the maternal control of seed maturation, jasmonate-signaled defense responses, and glandular trichome development. Plant Cell. 2004;16:126–43.
Goetz S, Hellwege A, Stenzel I, Kutter C, Hauptmann V, Forner S, et al. Role of cis-12-Oxo-phytodienoic acid in tomato embryo development. Plant Physiol. 2012;158:1715–27.
Schubert R, Dobritzsch S, Gruber C, Hause G, Athmer B, Schreiber T, et al. Tomato MYB21 acts in ovules to mediate jasmonate-regulated fertility. Plant Cell. 2019;31:1043–62.
Floris M, Mahgoub H, Lanet E, Robaglia C, Menand B. Post-transcriptional regulation of gene expression in plants during abiotic stress. Int J Mol Sci. 2009;10:3168–85.
Gass N, Glagotskaia T, Mellema S, Stuurman J, Barone M, Mandel T, et al. Pyruvate decarboxylase provides growing pollen tubes with a competitive advantage in petunia. Plant Cell. 2005;17:2355–68.
Selinski J, Scheibe R. Pollen tube growth: where does the energy come from? Plant Signal Behav. 2014;9:e977200.
Persia D, Cai G, Casino CD, Faleri C, Willemse MTM, Cresti M. Sucrose synthase is associated with the cell wall of tobacco pollen tubes. Plant Physiol. 2008;147:1603–18.
Besseau S, Hoffmann L, Geoffroy P, Lapierre C, Pollet B, Legrand M. Flavonoid accumulation in Arabidopsis repressed in lignin synthesis affects auxin transport and plant growth. Plant Cell. 2007;19:148–62.
Dobritsa AA, Lei ZT, Nishikawa SI, Urbanczyk-Wochniak E, Huhman DV, Preuss D, et al. LAP5 and LAP6 encode anther-specific proteins with similarity to chalcone synthase essential for pollen exine development in Arabidopsis. Plant Physiol. 2010;153:937–55.
Jin H, Song ZH, Nikolau BJ. Reverse genetic characterization of two paralogous acetoacetyl CoA thiolase genes in Arabidopsis reveals their importance in plant growth and development. Plant J. 2012;70:1015–32.
Ning L, Lin ZW, Gu JW, Gan L, Li YH, Wang H, et al. The initial deficiency of protein processing and flavonoids biosynthesis were the main mechanisms for the male sterility induced by SX-1 in Brassica napus. BMC Genomics. 2018;19:806.
Wang LX, Lam PY, Lui ACW, Zhu FY, Chen MX, Liu HJ, et al. Flavonoids are indispensable for complete male fertility in rice. J Exp Bot. 2020;71:4715–28.
Lloyd DG. The maintenance of gynodioecy and androdioecy in angiosperms. Genetica. 1975;45:325–39.
Wolf DE, Takebayashi N. Pollen limitation and the evolution of androdioecy from dioecy. Am Nat. 2004;163:122–37.
Furuyama T, Dzelzkalns VA. A novel calcium-binding protein is expressed in Brassica pistils and anthers late in flower development. Plant Mol Biol. 1999;39:729–37.
Malhó R, Monteiro QLD, Rato C, Camacho L, Dinis A. Signalling pathways in pollen germination and tube growth. Protoplasma. 2006;228:21–30.
Krzysztof Z, Rejón JD, Suárez C, Castro AJ, Alché JD, García MIR. Whole-Organ analysis of calcium behaviour in the developing pistil of olive (Olea europaea L.) as a tool for the determination of key events in sexual plant reproduction. BMC Plant Biol. 2011;11:150.
Wang XP, Su LX, Su JW. Distribution changes of calcium and programmed cell death in the pistil of litchi (Litchi chinensis Sonn.) flower during its development. J Plant Physiol Mol Biol. 2006;32:607–16.
Boisson-Dernier A, Franck CM, Lituiev DS, Grossniklaus U. Receptor-like cytoplasmic kinase MARIS functions downstream of CrRLK1L-dependent signaling during tip growth. Proc Natl Acad Sci U S A. 2015;112:12211–6.
Nissen KS, Willats WGT, Malinovsky FG. Understanding CrRLK1L function: cell walls and growth control. Trends Plant Sci. 2016;21:516–27.
Huck N, Moore JM, Federer M, Grossniklaus U. The Arabidopsis mutant feronia disrupts the female gametophytic control of pollen tube reception. Development. 2013;130:2149–59.
Zhu L, Chu LC, Liang Y, Zhang XQ, Chen LQ, Ye D. The Arabidopsis CrRLK1L protein kinases BUPS1 and BUPS2 are required for normal growth of pollen tubes in the pistil. Plant J. 2018;95:474–86.
Kou XB, Qi KJ, Qiao X, Yin H, Liu X, Zhang SL, et al. Evolution, expression analysis, and functional verification of Catharanthus roseus RLK1-like kinase (CrRLK1L) family proteins in pear (Pyrus bretchneideri). Genomics. 2017;109:290–301.
Geserick C, Tenhaken R. UDP-sugar pyrophosphorylase is essential for arabinose and xylose recycling, and is required during vegetative and reproductive growth in Arabidopsis. Plant J. 2013;74:239–47.
Fernández-Pérez F, Pomar F, Pedreño MA, Novo-Uzal E. The suppression of AtPrx52 affects fibers but not xylem lignification in Arabidopsis by altering the proportion of syringyl units. Physiol Plantarum. 2015;154:395–406.
Joshi CP, Chiang VL. Conserved sequence motifs in plant S-adenosyl-L-methionine-dependent methyltransferases. Plant Mol Biol. 1998;37:663–74.
Bosch M, Cheung AY, Hepler PK. Pectin methylesterase, a regulator of pollen tube growth. Plant Physiol. 2005;138:1334–46.
Gou JY, Miller LM, Hou GC, Yu XH, Chen XY, Liu CJ. Acetylesterase-mediated deacetylation of pectin impairs cell elongation, pollen germination, and plant reproduction. Plant Cell. 2012;24:50–65.
Röst HL, Rosenberger G, Navarro P, Gillet L, Miladinović SM, Schubert OT, et al. OpenSWATH enables automated, targeted analysis of data-independent acquisition MS data. Nat Biotechnol. 2014;32:219–23.
Yu GC, Wang LG, Han YY, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16:284–7.
Kanehisa M, Sato Y. KEGG Mapper for inferring cellular functions from protein sequences. Protein Sci. 2020;29:28–35.
Imran QM, Hussain A, Lee SU, Mun BG, Falak N, Loake GJ, et al. Transcriptome profile of NO-induced Arabidopsis transcription factor genes suggests their putative regulatory role in multiple biological processes. Sci Rep. 2018;8:771.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4:2–2.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Zhang C, Fu JX, Wang YG, Bao ZY, Zhao HB. Identification of suitable reference genes for gene expression normalization in the quantitative real-time PCR analysis of sweet osmanthus (Osmanthus fragrans Lour.). PLoS ONE. 2015;10:e0136355.
We thank Dongdong Sun and his company for their help in the treatment of experimental materials, as well as the editor and reviewers for their professional and meticulous work.
This work was supported by the Postgraduate Research & Practice Innovation Program of Jiangsu Province (grant number KYCX19_1073), Natural Science Foundation of Jiangsu Province (grant numbers SBK2020042924, BK20200786), Major project of natural science research in colleges of Jiangsu Province (grant number 20KJA220001), China Postdoctoral Science Foundation (grant numbers 2016 M590462, 2019 M651839), Nanjing Forestry University project funding (grant number GXL2018005), Innovation Fund for Young Scholars of Nanjing Forestry University (grant number CX2019029), National Natural Science Foundation of China (grant numbers 31300558, 32071782) and the Priority Academic Program Development (PAPD) of Jiangsu High Education Institutions, Jiangsu Province, China. Founding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
All methods for the use of plants comply with relevant institutional, national, and international guidelines and legislation.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Proteins identified by SWATH-MS in this study.
Functional classification of DEPs with the MapMAN BIN system.
Primers used for qRT-PCR experiment in this study.
Photosynthesis measurements of male and hermaphroditic O. fragrans at the full flowering stage.
About this article
Cite this article
Duan, YF., Zhang, C., Zhang, M. et al. SWATH-MS based quantitive proteomics reveal regulatory metabolism and networks of androdioecy breeding system in Osmanthus fragrans. BMC Plant Biol 21, 468 (2021). https://doi.org/10.1186/s12870-021-03243-8