- Research article
- Open Access
Seedless mutant ‘Wuzi Ougan’ (Citrus suavissima Hort. ex Tanaka ‘seedless’) and the wild type were compared by iTRAQ-based quantitative proteomics and integratedly analyzed with transcriptome to improve understanding of male sterility
BMC Geneticsvolume 19, Article number: 106 (2018)
Bud mutation is a vital method of citrus. ‘Wuzi Ougan’ (mutant type, MT) as a bud variant of ‘Ougan’ (wild type, WT) was first found in 1996 and has become popular because of its male sterility and seedless character. Previous analysis of its cytological sections and transcriptome revealed that the abnormal microsporogenesis that occurs before the tetrad stage of anther development might be the result of down-regulated oxidation-reduction biological processes in MT. To reveal the mechanism behind the male sterility in MT at the post-transcriptional stage, proteome profiling and integrative analysis on previously obtained transcriptome and proteome data were performed in two strains.
The proteome profiling was performed by iTRAQ (isobaric Tags for relative and absolute quantitation) analysis and 6201 high-confidence proteins were identified, among which there were 487 differentially expressed proteins (DEPs) in one or more developmental stages of anthers between MT and WT. The main functional subcategories associated with the main category biological process into which the DEPs were classified were sporopollenin biosynthesis process and pollen exine formation. The enriched pathways were phenylpropanoid biosynthesis, flavonoid biosynthesis, and phenylalanine metabolism. Moreover, there were eight pathways linked in terms of being related to phenylpropanoid metabolism. Eighteen important genes related to phenylpropanoid metabolism were also analysized by qRT-PCR (quantitative real time PCR). An integrative analysis of the fold change at the transcript (log2 FPKM ratios) and protein (log1.2 iTRAQ ratios) levels was performed to reveal the consistency of gene expression at transcriptional and proteomic level. In general, the expression of genes and proteins tended to be positively correlated, in which the correlation coefficients were 0.3414 (all genes and all proteins) and 0.5686 (DEPs and according genes).
This study is the first to offer a comprehensive understanding of the gene regulation in ‘Wuzi Ougan’ and its wild type, especially during the microsporocyte to meiosis stage. Specifically, the involved genes include those in phenylpropanoid biosynthesis, flavonoid biosynthesis, and phenylalanine metabolism, as determined by integrative transcriptome and proteome analysis.
China has an abundance of citrus resources and many cultivars derived from bud mutation . Compared with other breeding methods, such as hybrid breeding, radiation-induced mutation breeding, and transgenic breeding, the advantages of bud mutation are the short breeding cycle and the fast breeding speed . Therefore, genetic engineering has been applied to improve the quality of citrus cultivars by manipulating the molecular mechanism behind bud mutation. Mechanisms of bud mutation research have been focused on DNA methylation , retrotransposon insertion [4, 5], and on the structural and expressional difference of referred genes .
‘Wuzi Ougan’ (Citrus suavissima ‘seedless’, mutant type, MT), a bud variant of ‘Ougan’ (Citrus suavissima, wild type, WT) (Fig. 1), has almost of the same excellent characteristics as WT except of seedlessness. Seedlessness is an important commercial feature for fresh and processed fruit in the citrus industry. Previously studies showed that male sterility was an important reason for seedlessness in ‘Wuzi Ougan’ [7, 8]. Failure staining on mature pollen by both KI-I2 and FDA suggested the pollen abortion . Further observation by transmission electron microscopy (TEM) found that empty pollen grains were revealed during pollen maturation [7, 8]. In addition, microspore mother cells were found to be abnormal at the tetrad stage by scanning electron microscopy (SEM)  and suggested that the problem led to male sterility might occur at meiosis or earlier. Follow the observations, CsRad51, which was the gene in charge of double strand break (DSB) formation and DNA damage repair, was detected overexpression by qRT-PCR during microsporocytes to meiosis stage of anthers in MT when compared with WT . Lacking of DSBs was reported to be involved in failure to synapsis [10, 11]. However, both excess of DSBs and AtRad51 expression could be induced by radiation .
Male Sterility induced by genes (nuclear and cytoplasmic) mutation is genetic and widely documented in higher plants . In general, male sterility is characterized by failure of pollen grain development or function, and is either controlled by nuclear genes alone (genic male sterility) or regulated by the complementary action of nuclear and cytoplasmic genes (genetic-cytoplasmic male sterility) . In citrus, male sterility occurred naturally in bud sports [7, 8, 14, 15] or synthesized artificially by somatic protoplast fusion [16,17,18,19,20,21,22]. Those mutations and hybrids/cybrids exhibited usually failures of stamen development [23, 24] and some were evaluated by focusing on anther development [7, 8, 25, 26]. In addition to those with genic male sterility, some citrus cultivars with genetic-cytoplasmic male sterility exhibited the incidence of aborted anther of hybrid seedlings fits the segregation criterion controlled by major genes [25, 27]. Therefore, both male sterile types, those are conditioned by genes from either nuclear or nuclear organelle (eg. mutations of mitochondrial genes), could be characterized by genetically expression on RNA and protein level. In addition to sophisticated cell fusion, cross hybridization is attractively alternative and facilitated approach on introducing male sterility in seedless breeding of citrus.
In this connection, exploring efforts for improve understanding of post-transcriptional and -translational profiles have increased in citrus. Recently, vast of metabolic pathways and biological processes, as well as included genes and proteins, were reported in male sterility cultivars or lines [22, 28,29,30,31,32,33]. Among those mentioned in reports, genes or proteins have been focusing on the phenylpropanoid metabolism pathways. Phenylpropanoids metabolism is comprehensive network and includes several pathways and key enzymes. PAL (Phenylalanine ammonia-lyase) catalyzes the first step in the biosynthesis of the phenylpropanoid skeleton. PAL is documented to express predominantly in anthers [34,35,36] and induce sterile pollen when reduce the activity . Flavonoids biosynthesis is leading pathway located in phenylpropanoids metabolism. CHS (chalcone synthase) and CHI (chalcone isomerase) are rate-limiting enzymes and involved in male sterility when their expression were down-regulated [38,39,40]. Lignin biosynthesis is the last steps in phenylpropanoids network. CCR and CAD are essential for the monolignol pathway and confirmed by triple mutant of ccc exhibiting male sterility and severe dwarf phenotype .
Integrative analysis has become an important method to deeply explore the origin of the phenomenon. It has been applied to studying the mechanisms of fruit ripening , analyzing the regulatory roles of transcription factors , and comparing different breeds . In the present study, proteome profiling and integrative analysis combined with the use of previous transcriptome data obtained by RNA-Seq were performed to identify the candidate pathways and genes related to male sterility in ‘Wuzi Ougan’. A focus was placed on a pathway that potentially plays a leading role in male sterility in ‘Wuzi Ougan’, namely, phenylpropanoid metabolism. The obtained results should improve our understanding of male sterility in citrus.
Sample preparation and data source
‘Wuzi Ougan’ (Citrus suavissima Hort. ex Tanaka ‘seedless’) and its wild type were planted at Jin Chao Gang farm, Wenzhou, Zhejiang, China. The materials used for iTRAQ analysis were collected from three stages (S1, S2, S3) according to the length of floral bud : Stage 1 (the flower bud of the sporogonium, 1.2 mm < floral bud length < 1.6 mm), Stage 2 (the anther of the early microsporocyte, 2.0 mm < floral bud length < 2.4 mm) and Stage 3 (the anther of microsporocyte to meiosis, 2.4 mm < floral bud length < 2.8 mm). The flower buds were too small to form anthers and we collected the flower bud in Stage 1. Three independent technical replicates were set in each stage. The materials of qRT-PCR were collected from four stages according to the length of floral bud: I (the early microsporocyte, 2.0 mm < floral bud length < 2.4 mm), II (the microsporocyte to meiosis, 2.4 mm < floral bud length < 2.8 mm), III (the Tetrad, 2.8 mm < floral bud length < 3.1 mm) and IV (the pollen maturation, 6.5 mm < floral bud length < 7.1 mm). Anthers at microsporocytes were collected from ‘Wuzi Ougan’ and its wild type respectively for RNA-Seq. Raw sequencing reads obtained by transcriptomic sequencing were deposited in the Sequence Read Archive (SRA) database in NCBI under the accession number of PRJNA430695 (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA430695).
Protein extraction, iTRAQ labeled and LC-MS/MS
TCA-acetone method as an important method was selected to extract protein . The method of FASP Digestion according to universal sample preparation method for proteome analysis . iTRAQ labeled peptides were fractionated by SCX chromatography using the AKTA Purifier system (GE Healthcare) [47,48,49,50]. LC-MS/MS analysis was performed on a Q Exactive mass spectrometer (Thermo Scientific) that was coupled to Easy nLC (Proxeon Biosystems, now Thermo Fisher Scientific) for 60 min. The mass spectrometer was operated in positive ion mode. MS/MS spectra were searched using MASCOT engine (Matrix Science, London, UK; version 2.2) embedded into Proteome Discoverer 1.4 and the parameters were set in Table. 1.
The protein sequences of differentially expressed proteins were in batches retrieved from UniProtKB database in FASTA format. The retrieved sequences were locally searched against SwissProt database Citrus clementina (https://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Cclementina) to search homologue sequences from which the functional annotation can be transferred to the studied sequences. The final iTRAQ ratios of proteins were then normalized by the median average protein ratio of the equal mix of different labeled samples . In this work, each sequence was retrieved and loaded into Blast2GO (Version 3.3.5) for GO mapping and annotation. The FASTA protein sequences of differentially changed proteins were blasted against the online Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://geneontology.org/) to retrieve their KOs and were subsequently mapped to pathways in KEGG.
Quantitative real-time PCR (qRT-PCR) analysis
Eighteen DEPs selected from iTRAQ results were analyzed by using qRT-PCR. Total RNA was extracted from anther of ‘Ougan’ and ‘Wuzi Ougan’ in four stages (I, II, III, IV). The primers were designed by the GenScript Real-time PCR (TaqMan) Primer Design Center (https://www.genscript.com) (Table. 2). The Actin gene GU911361 as a reference was used to calculate the relative fold-differences based on comparative cycle threshold (2-ΔΔCt) values .
The iTRAQ analysis
Proteins play a role in most cellular functions and processes, and proteomics offers a direct and integrated perspective on cellular processes and networks. Thus, proteomic analysis of the anthers of MT and WT was performed using iTRAQ to explain CMS. Three developmental stages (the sporogonium, the early microsporocyte, the microsporocyte to meiosis) were analyzed in technical replicates. More than 380,000 mass spectra were collected from each replicate in this study. By performing data filtering to eliminate low-scoring spectra, 25,108, 24,243, and 24,238 unique peptides were collected from the three replicates (Table 3). A total of 4297 proteins (70%) were identified in each replicate simultaneously. Only 1058 proteins uniquely present in one replicate among the confirmed proteins (Fig. 2a). Most proteins were mapped with at least two unique peptides in this study (Fig. 2b), which indicated the high reliability of our data. In the study, the presence of a lot of proteins with unknown function identified revealed the complex nature of the regulatory network involved in development.
Based on iTRAQ quantitative proteomics, we confirmed 6201 high-confidence proteins, 487 of which were DEPs (1.2-fold change, p-value < 0.05) in one or more stages (Additional file 1: Table S1). There were 126, 340, and 83 differentially expressed proteins in S1, S2, and S3 between MT and WT, respectively (Additional file 2: Table S2). Among the DEPs, only 3 proteins were shared by three stages, 56 proteins were shared by two stages, while only 79, 286, and 63 proteins had differential expression profiles in S1, S2, and S3, respectively (Fig. 3a). There were 90, 227, and 28 proteins up-regulated in MT in S1, S2, and S3, and 36, 113, and 55 proteins up-regulated in WT in S1, S2, and S3, respectively (Fig. 3b–d).
COG and clusters of DEPs
To obtain a deeper understanding of the functions of the DEPs, a cluster of orthologous groups of proteins (COG) analysis was performed. As shown in Fig. 4, for more than half of the DEPs, which didn’t have COG or they belonged to the category with unknown function. A total of 181 DEPs were mainly classified into the categories of carbohydrate transport and metabolism; translation; ribosomal structure, and biogenesis; transcription; post-translational modification, protein turnover, and chaperones; secondary metabolite biosynthesis, transport, and catabolism; general function prediction only and signal transduction mechanisms, accounting for 14, 22, 15, 17, 16, 57, and 12 proteins, respectively. In addition, COG functions were largely enriched in the categories of translation, ribosomal structure, and biogenesis; transcription. Therefore, we consider that male sterility in MT may be influenced by genetic information processing, which might be related to abnormal meiosis. The analysis of these DEPs should improve our understanding of the mechanism behind male sterility in MT.
Protein expression during the stages of anther development was quantitatively analyzed in two strains. A total of 487 significantly differentially expressed proteins were identified and analyzed using the criteria of 1.2-fold change and a p-value < 0.05 (Fig. 5a). To obtain a deep understanding of the major trends of DEPs, we used a K-means method with a Euclidean distance metric using the software MultiExperiment Viewer. The DEPs were divided into 12 clusters (Fig. 5b). Among these DEP clusters, cluster 1 was particularly up-regulated at stage 1 and subsequently down-regulated. Clusters 2–6 showed the similar patterns, which were significantly up-regulated and reached their peak levels at stage 2. In clusters 2 and 3, DEPs were down-regulated at stage 3. A variety of DEPs in clusters 8–10 were significantly down-regulated in stages 1 and 2. Compared with cluster 10, cluster 12 showed the completely opposite trend.
GO and KEGG analyses
Gene Ontology (GO), an international standardized protein functional classification system, is a significant tool to classify the functions of a lot of proteins. GO analysis has been widely applied to predict the functions of proteins in many organisms. The GO database is composed of three ontologies: molecular function (MF), cellular component (CC), and biological process (BP).
In this study, among the proteins within the category of biological process found to be differentially expressed in ‘metabolic process’ and ‘cellular metabolic process’, indicating that these proteins might intensively take part in assimilation and/or dissimilation involved in pollen development. In cellular components, ‘cell’ and ‘cell part’ were the predominant components of the category. Regarding the molecular function category, the main enriched subcategories were ‘catalytic activity’ and ‘binding’ in the key period of floral bud development (Additional file 3 Table S3).
KEGG, a major biological process database, contains seven categories: metabolism, genetic information processing, environmental information processing, cell process, biological system, human diseases, and drug development. To further reveal the metabolic pathways of MT involved in male sterility, an analysis of the pathways of DEPs was performed and predicted a total of 487 proteins in 175 metabolic pathways (Additional file 4 Table S4). There were 19 common pathways in three stages (Table 4). Among those, four were phenylalanine-related pathways, namely, phenylpropanoid biosynthesis (ko00940); phenylalanine metabolism (ko00360); stilbenoid, diarylheptanoid, and gingerol biosynthesis (ko00945); and flavonoid biosynthesis (ko00941). Another four pathways also related to phenylalanine pathways were associated with numerous DEPs: phenylalanine, tyrosine, and tryptophan biosynthesis (ko00400); ubiquinone and other terpenoid-quinone biosynthesis (ko00130); TCA cycle (ko00020); and oxidative phosphorylation (ko00190). These eight pathways related to this process are described in a metabolic network (Fig. 6). Phenylalanine was supplied to the network of the phenylpropanoid metabolism by pathway ko00400. Some phenylalanine enters ko00360, which enables the supply of fumarate and succinate to ko00020 and trans-cinnamate to ko00130. The rest of the phenylalanine enters ko00940, which supplies cinnamic acid and cinnamoyl-CoA to the downstream pathways ko00941, ko00945, and ko00130.
Expression analysis of unigenes according to their cognate DEPs by qRT-PCR
Eighteen unigenes, according to their cognate DEPs involved in phenylpropanoid metabolism identified by iTRAQ, were used for qRT-PCR analysis in terms of four developmental stages (I, II, III, IV) in MT and WT. Those genes include two PAL genes (phenylalanine ammonia-lyase), one 4CL gene (4-coumarate-CoA ligase), two CAD genes (cinnamyl alcohol dehydrogenase), one CYP98A3 gene (cytochrome 98A3), four POD genes (peroxidase), two COMT genes (caffeoyl-O-methyltransferase), one CCoAOMT gene (caffeoyl-CoA O-methyltransferase), one CHS gene (chalcone synthase), one CHI gene (chalcone isomerase), one MIF gene (migration inhibitory factor), and two other genes. The expression levels of these genes basically showed a similar trend to the pattern of their cognate proteins (Fig. 7).
Integrative analysis of transcriptome and proteome
An integrative analysis of transcriptome and proteome data was performed to identify the regulation of metabolism in microsporogenesis in MT and WT. The transcriptome data was obtained previously (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA430695) from anthers at microsporocyte stage. Therefore, the proteome data at S3 stage (microsporocyte to meiosis) were drew out to compare with the previous transcriptome. In this connection, a total of 3809 proteins were used to analyze synthetically, and 2809 of them hit their corresponding mRNAs in the transcriptome. Among these, there were 1585 (56.4%) protein–mRNA pairs that showed a consistent trend in both transcriptome and proteome (Additional file 5 Table S5). Total of 2809 protein-mRNA pairs were divided into four categories based on the mRNA level (log2 FPKM ratios) and protein level (log1.2 iTRAQ ratios) expression profiles. Group I consisted of 2632 proteins showing no change between transcriptome and proteome. There were 12 proteins and 152 proteins changed either in the transcriptome or the proteome in Group III and Group IV, respectively. There were 13 differentially expressed proteins that showed a consistent trend in transcription and proteome in Group II. In general, the expression tendencies of all genes and all proteins were positively correlated, in which the correlation coefficient was 0.3414 (Fig. 8a). Additionally, expression changes at transcript (log2 FPKM ratios) and protein (log1.2 iTRAQ ratios) levels were plotted for DEPs detected by proteome and their according genes, and the correlation coefficient was 0.5686 (Fig. 8b).
Correlative DEPs succeeded mapping to genes from RNA-Seq data were mainly involved in sporopollenin biosynthetic process, pollen exine formation, nuclear speck, and transcription coactivator activity (p-value< 0.01, Fig. 9). The enriched KEGG pathways included those related to metabolic pathways, such as phenylpropanoid biosynthesis (ko00940), plant-pathogen interaction (ko04626), sulfur metabolism (ko00920), ABC transporters (ko02010), flavonoid biosynthesis (ko00941), phenylalanine metabolism (ko00360), and starch and sucrose metabolism (ko00500) (p-value < 0.1, Fig. 10). Pathways related to phenylpropanoid metabolism were also classified into four categories (Table 5). Phenylpropanoid biosynthesis (ko00940) includes the most amount proteins and many DEPs in Group II and IV.
The development of flower organs is a highly coordinated and irreversible phenomenon that involves a series of physiological, biochemical, and organoleptic changes . In the present study, proteome profiling was employed to investigate the differences between the mutant type and its wild type by iTRAQ technologies, where there were 487 DEPs showing variations in abundance during pollen development in the mutant. Using previous transcriptome data obtained by RNA-Seq, an integrative analysis was conducted to reveal the crosstalk between the transcriptome and the proteome in the regulation of pollen development, especially regarding male sterility in the mutant. Pathway enrichment was analyzed to reveal the importance of phenylpropanoid metabolism in this context, which has been widely reported to be involved in male sterility in plants [37, 38].
Crosstalk between transcriptome and proteome responsible for male sterility
Integrative analysis of expressed genes and their cognate proteins improved our understanding of male sterility in the mutant. Four categories based on the mRNA level (log2 FPKM ratios) and protein level (log1.2 iTRAQ ratios) expression profiles were analyzed in detail. In general, there were 13 DEPs of which the expression was consistent with that of homologous DEGs in Group II (Fig. 8a), which suggests that those affected biological processes or pathways might be responsible for the male sterility in the mutant. The involved biological processes include pollen development (Ciclev10028834m), ATPase activity (Ciclev10010223m), pollen exine formation (Ciclev10010274m), hydrolase activity (Ciclev10016019m), 4-coumarate-CoA ligase (Ciclev10031133m), and lipid binding (Ciclev10006307m). The formation of pollen floral organ development. Among these DEPs and DEGs, 4CL as a key enzyme in phenylpropanoid biosynthesis. qRT-PCR verification also showed that 4CL is up-regulated in early microsporocytes (Fig. 7c). In rice, overexpression of the 4-coumarate-CoA ligase (4CL) related gene OsAAE3 resulted in an increase in the content of H2O2 and led to programmed cell death (PCD) of the tapetum, which contributed to the suppression of floret development and decreased fertility rate of anther . Male sterility is associated with premature or delayed PCD of the tapetum, which is the innermost layer of the anther and produces substantial nutrition for the development of microspore mother cell [53, 54]. Additionally, there were 12 genes assigned to Group III with expression changed at the transcript level and 152 genes changed at the peptide level (Group IV). Groups IV include many DEPs, such as peroxidase activity (Ciclev10001726m), lignin biosynthetic process (Ciclev10029158m) and cinnamyl-alcohol dehydrogenase activity (Ciclev10028718m). Lack of CCoAOMT (Ciclev10029158m) and CAD (Ciclev10028718m) impacts on lignification in the anther endothecium, which has been shown to be responsible for the failure of anther dehiscence and pollen release [41, 55]. The related genes were down-regulated in MT in the first two stages, as verified by qRT-PCR (Fig. 7d, m). Moreover, we observed the failure of anther dehiscence and pollen release in MT by stereoscope (Fig. 1c-d). POD can prevent the excessive accumulation of MDA and ROS enzymes by decomposing H2O2 into O2. Male sterility was reported to be due to dysfunction of the balance of POD content in floral organs . POD is the last key enzyme in lignin synthesis, and down-regulated of Ciclev10001726m might ultimately influence lignin content in anther in ‘Wuzi Ougan’ before the tetrad stage in this research  (Fig. 7g).
This integrative analysis was focused on DEPs, and the genes associated with were according to the background of all genes obtained from transcriptome. However, the transcript levels indicated poorly their according proteins between MT and WT. Therefore, there were some possibilities for the inconsistent expression. It was possible that the materials used for sequencing were sampled from the same trees at different years. The expression levels between years changed to some extent. In addition to biological samplings, protein levels are regulated by posttranscriptional, translational, and post-translational mechanisms, and feedback loops exist between the processes of mRNA translation and protein degradation .
Potential pathways related to male sterility in ‘Wuzi Ougan’
Genes and proteins revealed by the integrative analysis shown in Fig. 11 were particularly associated with phenylpropanoid biosynthesis (ko00940), and that these cooperated with other pathways as well as produced various secondary metabolites.
Phenylpropanoid biosynthesis, which plays center role in the network of phenylpropanoid metabolism, is an important pathway involved in standing upright, long-distance transport of water , leaf growth  and the development of floral organs . Weaken activities of the enzymes associated with this pathway was proposed to male sterility [29,30,31,32]. In this study, the expression of PALs was down-regulated in MT before the tetrad stage (Fig. 7a, b). 4CL as a specific enzyme in ko00940 plays an irreplaceable role in this metabolic process . As revealed by qRT-PCR analysis, the expression pattern of 4CL showed a sharp drop from microsporocyte to meiosis (Fig. 7c). Different CoA was difficult to produce due to the lack of 4CL (Fig. 11). The deficiency of this enzyme delayed phenylalanine usage for lignin synthesis and flavonoid biosynthesis, as well as resulted in male sterility . Coumaroyl-CoA is converted to p-coumaroyl-CoA, caffeoyl-CoA, and feruloyl-CoA via a series of reactions catalyzed by 4CL, CYP98A3, and CCoAOMT (Fig. 11). Accordingly, down-regulated expression of 4CL, CYP98A3 and CCoAOMT (Fig. 7c, f, m) in MT during the microsporocyte to meiosis, as well as shared situation of p-coumaroyl-CoA, caffeoyl-CoA, and feruloyl-CoA in pathways (including ko00940, ko00941, and ko00360), were presumably to interrupt pollen development due to kinds of CoA deficiency. COMT and CCoAOMT are two anther-specific genes, CCoAOMT mainly plays a role in vascular tissues of young stamens, while COMT acts on the endothecium and the epidermal layer of stamens . In fact, many genes involved in the network of phenylpropanoid metabolism efficiently expressed in anther or tapetum cells, such as PAL [34,35,36], CHI , CHS [38, 39], COMT and CCoAOMT . Those most widely reported sterile genes were CHS homologous genes (D5 , YY2 , LAP5/LAP6 ). Furthermore, the CHS mutant was reported to change colors of anthers, for example, from yellow to white, and led to dysfunctional male sterility . It was also reported that overexpression of the CHS gene leads to male sterility along with intense pigmentation on the surface of the anther . Therefore, it is considered that suitable gene expression is a genetic buffering mechanism to ensure floral organ function and appearance throughout development .
In addition, there were two pathways, oxidation-reduction and TCA cycle, frequently reported in male sterility plants and involved in the network of phenylpropanoid metabolism. Mitochondria are necessary organelles for cellular energy production because they participate in many metabolic pathways including the pentose phosphate pathway (ko00030), oxidative phosphorylation (ko00195), and the TCA cycle (ko00020) . Via the ko00360 pathway, phenylalanine is converted to fumarate, succinate, or trans-cinnamate, which can ensure the functioning of ko00020 and ko00195. In this study, numerous DEPs were enriched in these three pathways, which were down-regulated at the early stage of pollen development (Table 6).For example, PckA (Ciclev10025088m), as a key gene in ko00020, were down-regulated in MT (Fig. 7r), which might influence the inversion of oxaloacetate to phosphoenolpyruvate .
Three pathways (ko00940, ko00941, ko00360) were found to be enriched in the integrative analysis, which coordinately function in the regulatory network. Among them, the pathway of ko00940 and ko00360 produce substances and ATPs for microsporocytes (Fig. 11). Interestingly, another pathway, ko00400, involved in primary metabolism was not enriched (Table 4), and proposed an abundant of phenylalanine in MT. Therefore, the insufficiency or down-regulation of 4CL from microsporocyte stage to meiosis stage might result in disorder of ko00940 and ko00360, and subsequently led to dysfunction of the network of phenylpropanoid metabolism. Such important alterations of 4CL expression might contribute to the male sterility in MT.
It is extremely complicated to unravel the mechanism behind male sterility, and KEGG enrichment analysis has become essential to survey this issue [69, 70]. In the present study, we considered that disordered phenylpropanoid metabolism along with the pathways up- or down-streamed are the focus point of male sterility in ‘Wuzi Ougan’.
This paper has presented a comprehensive analysis of male sterility in ‘Wuzi Ougan’ and its wild type. Through integrative transcriptome and proteome analysis, DEGs and DEPs were identified to be particularly associated with phenylpropanoid biosynthesis, flavonoid biosynthesis, and phenylalanine metabolism. Genes were analyzed by qRT-PCR to present their according proteins related to the phenylpropanoid metabolism. This study provides a deeper understanding of the mechanism behind male sterility in citrus as well as the bud mutant.
Cytoplasmic male sterility
Cluster of Orthologous Groups of proteins
Differentially expressed genes
Differentially expressed proteins
Double strand break
Reads Per kb per Million reads
Isotope coded affinity tag labeling
Isobaric Tags for relative and absolute quantitation
Kyoto encyclopedia of genes and genomes
Programmed cell death
Real-time polymerase chain reaction
Scanning electron microscopy
Transmission electron microscopy
Zhang M, Deng XX. Advances in research of citrus cultivars selected by bud mutation and the mechanism of formation of mutated characteristics. J Fruit Sci. 2006;23(6):871–6.
Yi K, Yan ZY, Liu Z, Wang DM, Yang F, Zhang JE. Review on identification and utilization of apple sport selection. J Fruit Sci. 2006;23(5):745–9.
Xu JD, Xu HD, Xu Q, Deng XX. Characterization of DNA methylation variations during fruit development and ripening of sweet Orange. Plant Mol Biol Rep. 2015;33:1–11.
Yao JL, Dong YH, Morris BAM. Parthenocarpic apple fruit production conferred by transposon insertion mutations in a MADS-box transcription factor. P Natl Acad Sci USA. 2001;98(3):1306–11.
Kobayashi S, Goto-Yamamoto N, Hirochika H. Retrotransposon-induced mutations in grape skin color. Science. 2004;304(5673):982.
Dong QL, Yan ZY, Liu Z, Yao YX. Early ripening events caused by bud mutation in Beni shogun apple. Russ J Plant PhysL+. 2011;58(3):439–47.
Hu ZY, Zhang M, Wen QG, Wei J, Yi HL, Deng XX. Xu XH. Abnormal microspore development leads to pollen abortion in a seedless mutant of ‘Ougan’ mandarin (Citrus suavissima Hort. ex Tanaka). J Amer Soc Hort Sci. 2007;132(6):777–82.
Zhang C, Zhang M, Zhu Q, Liu ZH, Yan FH, Wu LH, Xu XH, Zhou XY, Chen X. Cytological observation of pollen development in 'Ougan'(Citrus suavissima Hort. ex Tanaka) and its seedless mutant. J Fruit Sci. 2014;31(2):265–9.
Zhu MM, Zhang C, Chang AL, Dang WY, Zhou CH, Yu DH, Wu YY, Zhang M. Expression of microsporocyte meiosis with special genes RAD51 and MS1 in Citrus suavissima ‘Seedless’. J Zhejiang A & F University. 2016;33(6):921–7.
Li WX, Yang XH, Lin ZG, Timofejeva L, Xiao R, Makaroff CA, Ma H. The atrad51c gene is required for normal meiotic chromosome synapsisand double-stranded break repair in Arabidopsis. Plant Physiol. 2005;138:965–76.
Kou Y, Chang Y, LI X, Wang S. The rice RAD51C gene is required for the meiosis of both femaleand male gametocytes and the DNA repair of somatic cells. J Exp Bot. 2012;63(4):5323–35.
Osakabe K, Yoshioka T, Ichikawa H, Toki S. Molecular cloning and characterization of Rad51-like genes fromArabidopsis thaliana. Plant Mol Biol. 2002;50(1):71–81.
Kaul ML. Genic male sterility. In: Male sterility in higher plants. Monographs on Theoretical and Applied Genetics, vol 10. Springer, Berlin, Heidelberg: Springer-Verlag; 1988.
Vardi A, Levin I, Carmi N. Induction of seedlessness in citrus. From classical techniques to emergingbiotechnological approaches. J Amer Soc Hortic Sci. 2008;133:117–26.
Iwamasa M. Studies on sterility in genus citrus with special reference to the seedlessness. Bull HortRes Sta Ser B. 1966;6:1–77.
Guo WW, Xiao SX, Deng XX. Somatic cybrid production via protoplast fusion for citrus improvement. Sci Hortic. 2013;163(5):20–6.
Yamamoto M, Matsumoto R, Okudai YY. Aborted anthers of Citrus result from gene-cytoplasmic male sterility. Sci Hotric. 1997;70(1):1–14.
Guo WW, Prasad D, Cheng YJ, Serrano P, Deng XX, Grosser JW. Targeted cybridization in citrus: transfer of Satsuma cytoplasm to seedy cultivars for potential seedlessness. Plant Cell Rep. 2004;20(10):752–8.
Cai XD, Fu J, Chen CL, Guo WW. Cybrid/hybrid plantsregenerated from somatic fusions between male sterile Satsuma mandarin and seedy tangelos. Sci Hortic. 2009;122(2):323–7.
Xiao SX, Biswas MK, Li MY, Deng XX, Xu Q, Guo WW. Production and molecular characterization of diploid and tetraploid somatic cybrid plants between male sterile Satsuma mandarin andseedy sweet orange cultivars. Plant Cell Tissue Organ Cult. 2014;116(1):81–8.
Wang L, Pan ZY, Guo WW. Proteomic analysis of leaves from a diploid cybrid produced by protoplast fusion between Satsuma mandarin and pummelo. Plant Cell Tissue Organ Cult. 2010;103(2):165–74.
Zheng BB, Wu XM, Ge XX, Deng XX, Grosser JW, Guo WW. Comparative transcript profiling of a male sterile cybrid pummelo and its fertile type revealed altered gene expression related to flower development. PLoS One. 2012;7(8):e43758.
Hanson MR, Bentolila S. Interactions of mitochondrial and nuclear genes that affect male gametophyte development. Plant Cell. 2004;16:154–69.
Sota F, Toriyama K. Genome barriers between nuclei and mitochondria exemplified by cytoplasmic male sterility. Plant Cell Physiol. 2008;49(10):1484–94.
Nakano M, Nesumi H, Yoshioka T, Yoshioka T. Segregation of plants with undeveloped anthers among hybrids derived from the seed parent, ‘kiyomi’ (Citrus unshiu × C sinensis). J Japan SocHort Sci. 2001;70:539–45.
Yamamoto M, Matsumoto R, Okudai N, Yamada Y. Aborted anthers of Citrus result from gene-cytoplasmic male sterility. Sci Hortic. 1997;70:9–14.
Yamamoto M, Okudai N, Matsumoto R. Study on the inheritance of aborted anthers in citrus using seedparents having aborted anthers. J Japan Soc Hort Sci. 1992;60:791–7.
Zheng BB, Fang YN, Pan ZY, Sun L, Deng XX, Grosser JW, Guo WW. iTRAQ-based quantitative proteomics analysis revealed alterations of carbohydrate metabolism pathways and mitochondrial proteins in a male sterile Cybrid Pummelo. J Proteome Res. 2014;13:2998–3015.
Li J, Yang S, Gai J. Transcriptome comparative analysis between the cytoplasmic male sterile line and fertile line in soybean ( Glycine max (L.) Merr.). Genes Genom. 2017;39(10):1117–27.
Rhee SJ, Seo M, Jang YJ, Cho S, Lee GP. Transcriptome profiling of differentially expressed genes in floral buds and flowers of male sterile and fertile lines in watermelon. BMC Genomics. 2015;16:914.
Sikora M, Ivancevic B, Jambrosic K. Identification of miRNAs with potential roles in regulation of anther development and male-sterility in 7B-1 male-sterile tomato mutant. BMC Genomics. 2015;16(1):878.
Ye J, Duan Y, Hu G, Geng X, Zhang G, Yan P, Liu Z, Zhang L, Song X. Identificationof candidate genes and BiosynthesisPathways related to Fertility Conversion by Wheat KTM3315A transcriptome profiling. Front Plant Sci. 2017;8:449.
Chen Z, Li HX, Ma XL, Dong HX, Ma D, Zeng HL. Proteome alterations of reverse photoperiod-sensitive genic male sterile rice (Oryza sativa L.) at fertility transformation stage. Genes Genom. 2014;36(6):711–26.
Rittscher M, Wiermann R. Occurrence of phenylalanine ammonia-Lyase (PAL) in isolated tapetum cells of Tulipa anthers. Protoplasm. 1983;118(3):219–24.
Beerhues L, Forkmann G, Schopker H, Stotz G, Wiermann R. Flavanone 3-hydroxylase and Dihydroflavonol oxygenase activities in anthers of Tulipa: The Significance of the Tapetum Fraction in Flavonoid Metabolism. J Plant Physiol. 1989;133(6):743–6.
Kehrel B, Wiermann R. Immunochemical localization of phenylalanine ammonia-lyase and chalcone synthase in anthers. Planta. 1985;163:183–90.
Kishitani S, Yomoda A, Konno N, Tanaka Y. Involvement of phenylalanine ammonia-lyase in the development of pollen in broccoli (Brassica oleracea L.). Sex Plant Reprod. 1993;6:244–8.
Zheng R, Yue SJ, Xu XY, Liu JY, Xu Q, Wang XL, Han L, Yu DY. Proteome analysis of the wild and YX-1Male sterile mutant anthers of wolfberry (Lycium barbarum L.). PLoS One. 2012;7(7):e41861.
Napoli CA, Fahy D, Wang HY, Taylor LP. White anther. A petunia mutant that abolishes pollen flavonol accumulation, induces male sterility, and is complemented by a chalcone synthase transgene. Plant Physiol. 1999;120(2):615–22.
Van Tunen AJ, Hartman SA, Mur LA, Mol JNM. Regulation of chalcone flavanone isomerase (CHI) gene expression in Petunia hybrida: the use of alternative promoters in corolla, anthers and pollen. Plant Mol Biol. 1989;12:539–51.
Thévenin J, Pollet B, Letarnec B, Saulnier L, Gissot L, Maia-Grondard A, Lapierre C, Jouanin L. The simultaneous repression of CCR and CAD, two enzymes of the lignin biosynthetic pathway, results in sterility and dwarfism in Arabidopsis thaliana. Mol Plant. 2011;4(1):70–82.
Wu J, Xu Z, Zhang Y, Chai L, Yi H, Deng XX. An integrative analysis of the transcriptome and proteome of the pulp of a spontaneous late-ripening sweet orange mutant and its wild type improves our understanding of fruit ripening in citrus. J Exp Bot. 2014;65(6):1651–71.
Shelton LM, Lister A, Walsh J, Jenkins RE, Wong MHL, Rowe C, Ricci E, Ressel L, Fang Y, Demougin P. Integrated transcriptomic and proteomic analyses uncover regulatory roles of Nrf2 in the kidney. Kidney Int. 2015;88(6):1261–73.
Yang H, Xu X, Ma H, Jiang J. Integrative analysis of transcriptomics and proteomics of skeletal muscles of the Chinese indigenous Shaziling pig compared with the Yorkshire breed. BMC Genet. 2016;17(1):80.
Novo JVJN. Plant proteomics: methods and protocols. By HervÃ© Thiellement, Michel Zivy, Catherine Damerval and ValÃ©rie Mechin (Eds.) . Biotechnol J. 2010;2(5):642.
Wisniewski J, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nat Methods. 2009;6(5):359–62.
Elias JE, Gygi SP. Target-decoy search strategy for increased confidence in large-scale protein identifications by mass spectrometry. Nat Methods. 2007;4(3):207–14.
Gan CS, Chong PK, Pham TK, Wright PC. 2007. Technical,experimental, and biological variations in isobaric tags for relativeand absolute quantitation (iTRAQ). J Proteome Res. 2007;6:821–7.
Unwin RD, Griffiths JR, Whetton AD. Simultaneous analysisof relative protein expression levels across multiple samples using iTRAQ isobaric tags with 2D nano LC-MS/MS. Nat Protoc. 2010;5:1574–82.
Noirel J, Evans C, Salim M, Mukherjee J, Ow SY, Pandhal J, Pham TK, Biggs CA, Wright PC. Methods in quantitative proteomics: setting iTRAQ on the right track. Curr Proteomics. 2011;8:17–30.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8.
Liu H, Guo ZH, Gu FW, Ke SW, Sun DY, Dong SY, Liu W, Huang M, Xiao WM, Yang GL, Liu YZ, Guo T, Wang H, Wang JF, Chen ZQ. 4-Coumarate-CoA ligase-like gene OsAAE3 negatively mediates the Rice blast resistance, floret development and lignin biosynthesis. Front Plant Sci. 2017;7:271–82.
Ku SJ, Yoon HJ, Suh HS, Chung YY. Male-sterility of thermosensitive genic male-sterile rice is associated with premature programmed cell death of the tapetum. Planta. 2003;217:559–65.
Yang XJ, Liang WQ, Chen MJ, Zhang DB, Zhao XX, Shi JX. Rice fatty acyl-CoA synthetase OsACOS12 is required for tapetum programmed cell death and male fertility. Planta. 2017;246:105–22.
Fellenberg C, Ohlen MV, Handrick V, Vogt T. The role of CCoAOMT1 and COMT1 in Arabidopsis anthers. Planta. 2012;236:51–61.
Li J, Dai X, Li L, Jiao Z, Huang Q. Metabolism of reactive oxygen species in cytoplasmic male sterility of Rice by marking upmost Pulvinus interval. Appl Biochem Biotechnol. 2015;175(3):1263–9.
Baucher M, Van Doorsselaere J, Van Montagu M, Inzé D. Molecular characterization of lignin biosynthesis genes in poplar. In: Fritig B, Legrand M, editors. Mechanisms of plant defense responses. Developments in plant pathology, vol 2. Dordrecht: Springer; 1993. p. 1993.
Vogel C, Marcotte EM. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat Rev Genet. 2012;13:227–32.
Li Y, Kim JI, Pysh L, Chapple C. Four isoforms of Arabidopsis thaliana 4-Coumarate. CoA ligase have overlapping yet distinct roles in phenylpropanoid metabolism. Plant Physiol. 2015;169:2409–21.
Lee SH, Park YJ, Park SU, Lee SW, Kim SC, Jung CS, Jang JK, Hur Y, Kim YB. Expression of genes related to Phenylpropanoid biosynthesis in different organs of Ixeris dentata var. albiflora. Molecules. 2017;22:901.
Kumar A, Ellis BE. The phenylalanine ammonia-Lyase gene family in raspberry. Structure, Expression, and Evolution. Plant physiol. 2001;127:230–9.
Jung JH, Kannan B, Dermawan H, Moxley GW, Altpeter F. Precision breeding for RNAi suppression of a major 4-coumarate:coenzyme a ligase gene improves cell wall saccharification from field grown sugarcane. Plant Mol Biol. 2016;92:505–17.
Beliaeva RG, Evdokimova LI. Variability of flavonol contents during floral morphogenesis in Papaver somniferum L. Russ J Dev Biol. 2004;35(1):10–5.
Zheng HH, Qu LJ, Liu MH, Zhang Y, Shen YP, Wei JM, Pan NS, Gu HY, Chen ZL. An anther-specific chalcone synthase-like gene D 5 related to rice pollen development. Chin Sci Bull. 2000;45(21):1921–6.
Hihara Y, Hara C, Uchimiya H. Isolation and characterization of two cDNA clones for mRNAs that are abundantly expressed in immature anthers of rice (Oryza sativa L.). Plant Mol Biol. 1996;30(6):1181–93.
Dobritsa AA, Lei ZT, Nishikawa SI, Urbanczy-Wochniak E, Huhman DV, Preuss D, Sumner LW. 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.
Li Y, Zhang T, Shen ZW, Xu Y, Li JY. Overexpression of maize anthocyanin regulatory gene Lc affects rice fertility. Biotechnol Lett. 2013;35:115–9.
Das SK, Chakrabartty PK. Characterization of a (pckA) mutant of the non-nodulating bacterium rhizobium pusense NRCPB10 induced by transposon Tn5 mutagenesis. Biotech. 2012;2(4):311–8.
Qiu LL, Jiang B, Fang J, Shen YK, Fang ZX, Saravana KRM, Yi KK, Shen CJ, Yan DL, Zheng BS. Analysis of transcriptome in hickory (Carya cathayensis) and uncover the dynamics in the hormonal signaling pathway during graft process. BMC Genomics. 2016;17:935.
Huang RM, Huang YJ, Sun ZC, Huang JQ, Wang ZJ. Transcriptome analysisof genes involved in lipid biosynthesis in the developing embryo of pecan (Carya illinoinensis). J Agr Food Chem. 2017;65(20):4223–36.
We are grateful to Prof. Wen-Wu Wu for helping to deal with the bioinformatics analysis, and two anonymous reviewers for providing useful comments on a previous version of the manuscript.
This research was financially supported by Zhejiang Province Public Welfare Technology Application Research Project (CN) (LGN18C160005); Team scientific and technical envoy project ((2013)215–122); The National Natural Science Foundation of China (31000897, 31301811); Science Technology Department of Zhejiang Province (No. 2016C02052–1); The earmarked fund for China Agriculture Research System (CARS-27). The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
Raw sequencing reads obtained by transcriptiomic sequencing were deposited in the Sequence Read Archive (SRA) database in NCBI under the accession number of PRJNA430695 (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA430695).
Ethics approval and consent to participate
All the plant materials used in this research were provided by State Key Laboratory of Subtropical Silviculture, Zhejiang Agriculture and Forestry University, China. The field experiments were conducted under local legislation and permissions.
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.
Table S1. All identification proteins in three stages. “N” indicates MT in three stages, and “Y” indicates WT in three stages. (XLSX 2862 kb)
Table S2. Annotated DEPs and their expression ratio (N/Y). “N” and “Y” indicate MT and WT, respectively. (XLS 181 kb)
Table S3. GO term of DEPs. (XLS 127 kb)
Table S4. KEGG term of DEPs. (XLS 98 kb)
Table S5. Integrated_Table_mRNA_protein_Anno. (XLSX 999 kb)