Differential Proteomic Analysis of Anthers between Cytoplasmic Male Sterile and Maintainer Lines in Capsicum annuum L

Cytoplasmic male sterility (CMS), widely used in the production of hybrid seeds, is a maternally inherited trait resulting in a failure to produce functional pollen. In order to identify some specific proteins associated with CMS in pepper, two-dimensional gel electrophoresis (2-DE) was applied to proteomic analysis of anthers/buds between a CMS line (designated NA3) and its maintainer (designated NB3) in Capsicum annuum L. Thirty-three spots showed more than 1.5-fold in either CMS or its maintainer. Based on mass spectrometry, 27 spots representing 23 distinct proteins in these 33 spots were identified. Proteins down-regulated in CMS anthers/buds includes ATP synthase D chain, formate dehydrogenase, alpha-mannosidas, RuBisCO large subunit-binding protein subunit beta, chloroplast manganese stabilizing protein-II, glutathione S-transferase, adenosine kinase isoform 1T-like protein, putative DNA repair protein RAD23-4, putative caffeoyl-CoA 3-O-methyltransferase, glutamine synthetase (GS), annexin Cap32, glutelin, allene oxide cyclase, etc. In CMS anthers/buds, polyphenol oxidase, ATP synthase subunit beta, and actin are up-regulated. It was predicted that male sterility in NA3 might be related to energy metabolism turbulence, excessive ethylene synthesis, and suffocation of starch synthesis. The present study lays a foundation for future investigations of gene functions associated with pollen development and cytoplasmic male sterility, and explores the molecular mechanism of CMS in pepper.


Introduction
Chili pepper (Capsicum spp.) is one of the most widely cultivated vegetables or spice crops around the world. Heterosis in pepper has been shown to increase fruit yield by up to 30% [1]. Currently, the breeding of hybrid seeds depends on hand pollination, which is costly and difficult to ensure seed purity. Selecting male sterile lines, with good combining ability, are an efficient measure to breeders in commercial hybrid seed production.
Cytoplasmic male sterility (CMS) is caused by an interaction between nuclear and mitochondrial genes. Unlike genetic male sterility, CMS belongs to maternal inheritance and is a product of mitochondrial DNA (mtDNA) evolution that gives birth to new genes for better or for worse. Compared with the male fertile mtDNA, those CMS-associated genes are usually some novel open reading frames (ORFs) or chimeric mitochondrial genes, which cause male sterility in the female lines by expressing novel RNA and/or polypeptides [2,3]. In addition, nuclear restorer (Rf) genes in the male parent reversely affect the mutant mitochondrial genes and restore pollen viability in their F 1 hybrid [4][5][6]. Therefore, this CMS/Rf system not only represents a valuable tool in large scale production of commercial hybrid seeds but is also an excellent model to investigate the interaction between nuclear and mitochondrial genomes.
Male sterile cytoplasm of pepper was first documented by Martin and Grawford (1951) and characterized by Peterson (accession USDA P.I. 164835) [7,8], which is the only known source of male sterile cytoplasm until now. In order to explore the molecular mechanism of pepper CMS, orf456 or orf507 had been identified as a candidate gene [9,10]. Orf456 contains a novel ORF at the 3'-end of the cox II gene, which is co-transcribed with the cox II gene encoding a subunit of cytochrome c oxidase in mitochondria and inducing male sterility in Arabidopsis thaliana [9]. Gulyas found that a cytosine (C) deletion at the original stop codon of orf456 leads to continue to extend its translation to a new termination codon downstream, and finally produces a total length of 507 bp RNA, named orf507 [10]. However, they did not reveal the full mechanism of the CMS phenomenon in chili pepper. Recently, Luo et al. (2013) had identified the molecular basis of male sterility in the rice CMS-WA system. They thought that a new mitochondrial gene, WA532, accumulates preferentially in the anther tapetum, thereby inhibiting COX11 function in peroxide metabolism and triggering premature tapetal programmed cell death (PCD) and consequent pollen abortion [6].
As the functional molecules of living cells, proteins play important roles in the development of cell, organ or tissue by the interaction with other molecules or modifications. To study their functions, proteomics is one of the most powerful molecular tools in describing the proteomes of different organelles or tissues [11][12][13]. In recent years, the proteomic approach has been also applied to explore anther development and pollen reproduction of plants, such as in rice [11], tomato [12], Brassica napus [13], wolfberry [14], tobacco [15], etc. However, there is little proteomic information available for answering the unresolved mysteries and clarifying the mechanism of CMS in pepper.
In our study, a comparative proteomic approach was used to identify differentially expressed proteins of anthers during reproductive development of CMS in cytoplasmic male sterile and maintainer Lines in Capsicum annuum L. The objective was to identify some specific proteins related to CMS, and further elucidate their possible biological roles and their potential effects on anther development and pollen fertility. This lays a foundation for the understanding of molecular mechanism of CMS at the proteomic level, and makes the better use of CMS/Rf system in the production of pepper hybrid seeds.

Abundance Analysis of Differentially Expressed Proteins in the Normal and Sterile Cytoplasm
The objective of this study was to identify differentially expressed proteins between sterile and fertile anthers (i.e., NA3 and NB3 isogenic lines). Proteins were firstly separated by the isoelectric focusing (IEF) on a linear gradient ranging from pH 4.0 to 7.0. Based on the size of molecular masses, proteins were separated by SDS-PAGE and then stained with CBB G250 and silver nitrate ( Figure 1A,B). As shown in Figure 1C, some differential spots were found in the enlarged portion of gels. Only proteins detected in all two reps were recorded. Over 1200 non-redundant proteins were found in the anther proteome. By contrast, 33 different proteins were determined between A and B lines according to their volume ratios (more than 1.5-fold, p ≤ 0.05). Proteins with less than 1.5-fold change in volume were not considered in this study. Of these, nine proteins had much higher concentrations in sterile anthers/buds ( Figure 1A), 11 proteins had lower concentrations in sterile anthers/buds, and 13 proteins were not translated in sterile anthers/buds ( Figure 1B).

Identification and Functional Classification of Differentially Expressed Proteins
Thirty-three different proteins detected were excised from the preparative gels and digested with trypsin. Of these, 27 protein spots (81.82%) were successfully identified by MALDI-TOF/TOF MS analysis except for A5, A6, A8, A9, B8, and B18. Compared with NCBInr and pepper ESTdatabases, 23 distinct proteins identified were noted with ID, MW, pI value, MASCOT score, number of matched peptides, and sequence coverage ratio (%) (Tables 1,S1). In male sterile cytoplasm, A01 and A02 were identified as Polyphenol oxidase F, and ATP synthase subunit beta was showed in A03 and A04, In normal cytoplasm, multiple spots were also detected, such as ATP synthase CF1 alpha chain (B03 and B04), Glutelin (B14 and B15), glutamine synthetase (B09, B10, and B12), Actin (B06 and A07), and triose phosphate isomerase cytosolic isoform (B19 and B21). These results indicated that some proteins species with post-translational modifications may play a certain role in multiple organelles such as mitochondria, plastids, and cystol. and its maintainer NB3 line (B); some of the areas in NA3 and NB3 gels with differential spots have been enlarged below (C). Proteins were visualized by silver nitrate staining. The arrowed and numbered spots in the image were differentially expressed proteins. Molecular markers (kDa) are shown on the right. These results were repeated for three times.  In male sterile anthers/buds, higher expression proteins related to fertile anthers/buds include polyphenol oxidase F, ATP synthase subunit beta, and actin. While the down-regulated proteins have ATP synthase D chain, formate dehydrogenase, alpha-mannosidas, RuBisCO large subunit-binding protein subunit beta, chloroplast manganese stabilizing protein-II, glutathione S-transferase, adenosine kinase isoform 1T-like protein, putative DNA repair protein RAD23-4, and putative caffeoyl-CoA 3-O-methyl transferase. Proteins only showed in fertile anthers/buds contain ATP synthase CF1 alpha chain, triose phosphate isomerase (TPI), 23 kDa polypeptide of the oxygen evolving complex of photosystem II, glutamine synthetase (GS), annexin Cap32, glutelin, and allene oxide cyclase.
Based on their putative physiological functions compared with the gene annotations in NCBI and some relevant literatures, 27 identified proteins are classified into nine functional groups, such as respiration and energy pathway, carbohydrate metabolism, antioxidative reactions, photosynthesis, amino acid metabolism, singal transduction, cytoskeleton, storage protein, and others. The largest functional category of proteins drop into respiration and energy pathway (22.22%, two in sterile anthers and four in fertile anthers), followed by those proteins participating in carbohydrate metabolism (11.11%, three only expressed in maintainer anthers), photosynthesis (11.11%, three only expressed in maintainer anthers), antioxidative reactions (11.11%, two in CMS and one in maintainer anthers), and amino acid metabolism (11.11%, three only expressed in maintainer anthers).
More differentially expressed proteins were in fertile anthers/buds than in sterile anthers/buds. Polyphenol oxidase F (PPO, showed in A01 and A02 spots, Tables 1,S1) was a major protein showed significantly up-regulation in male sterile cytoplasm. Therefore, this enzyme was selected to further confirm changes using enzyme activities.

Evaluation of PPO and EFE Activities
The activity of PPO and ethylene forming enzyme (EFE) enzymes in CMS and maintainer line were analyzed at different stages of pollen development. As showed in Figure 2, the activities of PPO and EFE are always higher in CMS line NA3 than in maintainer line NB3. Especially, in the stage II and stage III, the activity of PPO was nearly two times higher in CMS compared with that in maintainer line ( Figure 2). In stage III, the activity of EFE was nearly three times higher in CMS than that in maintainer line.

Discussion
The use of CMS/Rf system has been widely applied to F 1 hybrid seed production of chili pepper. However, the molecular mechanism of CMS remains poorly understood. As a powerful molecular tool, proteomics has been used to study anther development and pollen reproduction in other plants [11][12][13][14][15]. In our study, the proteomic approach was performed to identify differentially expression proteins in CMS and its maintainer line to understand the molecular mechanism of pepper male sterility. Thirty-three proteins displayed differential expression levels between male sterile and normal anthers/buds. These down-regulated proteins identified using MS in CMS anthers/buds included RuBisCO large subunit-binding protein subunit beta, ATP synthase CF1 alpha chain, alpha-mannosidas, formate dehydrogenase, glutamine synthetase, caffeoyl-CoA 3-O-methyltransferase, triose phosphate isomerase, and glutathione S-transferase. The other proteins up-regulated in male sterile anthers/buds include polyphenol oxidase, ATP synthase subunit beta, and actin. The potential roles of some differential expression proteins in anther and pollen development are discussed as following.

Proteins Involved in Respiration and Energy Pathway
Six identified proteins belong to these groups of respiration and energy pathways. In particular, ATP synthase CF1 alpha chain (B03), ATP synthase CF1 alpha subunit (B04), formate dehydrogenase (B07), and ATP synthase D chain (B24) were down-regulated in the CMS line. In CMS lines of rice, Brassica napus, and wolfberry, two enzymes are also involved in respiration and energy metabolism [11,13,14].
Anther and pollen development is a complex process that depends on a series of well-coordinated metabolic and structural changes. During flower development, male meiotic cells and developing microspores require higher energy, particularly much higher mitochondrial biosynthetic products than vegetative development and growth. In maize, the number of mitochondria per cell increases 20-and 40-fold in meiocytes and tapetal cells during pollen development [16]. Microspore death was observed, when antisense of mitochondrial pyruvate dehydrogenase El down-regulated the alternative oxidase by suppressing pyruvate dehydrogenase in tobacco tapetum [17]. Many CMS-associated genes include portions of ATP synthase subunits or closely linked to normal ATP synthase subunit genes. The ATP synthase, a key enzyme for the synthesis of ATP for cellular biosynthesis, comprises three parts: F0, F1, and FA [18,19]. It has been reported that alterations of mitochondrial-encoded subunits of the F 0 F 1 -ATP synthase induce CMS in plants. For instance, in sunflower, sterile plants expressing mitochondrial ORF522 showed a specific decreased ATP synthase activity [20]. The chimeric mitochondrial ORF522 shares sequence similarity with ORFB, a plant-type ATP8, which might result in competition between two proteins leading to decreased activity of the F 0 F 1 -ATP synthase complex. In CMS-HongLian rice, sterility is associated with the expression of atp6-OrfH79, which might disturb the formation of the F 0 F 1 -ATPase complex, resulting in decreased activity of ATPase and pollen abortion [21]. The Ψatp6-2 in the pepper CMS line HW203A was also down-regulated, but it was up-regulated in the maintainer line. The corresponding F 0 F 1 -ATPase activity in the CMS line was gradually decreased along with the development of the anther, while maintainer line, F 0 F 1 -ATPase activity sharply decreased after the stage of sporogenous cell, but gradually increased following the tetrad stage [22]. One hypothesis on the mechanism of CMS is that the sterile lines are unable to meet the increased energy demand in the form of ATP during anther development, that leading to pollen abortion [20,23].

Proteins Involved in Carbohydrate Metabolic Pathways
Carbohydrate metabolism is one of the most basic metabolic pathways in biological metabolism. Its main physiological function is to provide required energy and carbon sources. Two proteins expressed only in NB3 anthers that were involved in carbohydrate metabolism were triose phosphate isomerase (TPI, B19, and B21) and alpha-mannosidase (B05). TPI plays an important role in the glycolytic pathway, which is the glycolytic enzyme. It can catalyze the reversible inter conversion of glyceraldehyde 3-phosphate and dihydroxyacetone phosphate [24]. Moreover, glycolysis belongs to the common method of aerobic and anaerobic respiration, which is critical for energy supply.
Carbohydrate plays a key role in anther and pollen development. It not only provides nutrition for anther development, but also affects anther and pollen development as a signal substance. Deng reported that the content of ATP and NADH decreased significantly in the maintainer line at the stage of pollen abortion, meanwhile, the expression level of TPI reduced obviously in anther [25]. When restorer genes (Rfs) were crossed to sterile lines, TPI activity and the expression level were increased and maintained at the normal levels. Therefore, we conclude that decreased abundance of these enzymes in CMS anthers down regulated sugar and starch concentrations, two molecules are necessary for biosynthesis and energy balance and induce male sterility.

Proteins Involved in Photosynthesis
Three spots (B01, B13, and B23) showed higher concentrations in fertile anthers/buds. They were identified as chloroplastic RuBisCO large subunit-binding protein subunit beta, chloroplast manganese stabilizing protein-II, and 23 kDa polypeptide of the oxygen evolving complex of PSII. Photosynthesis is dependent on the coordinated activity of the chloroplast, mitochondrial, and nuclear components [26]. Photosynthesis provides some substrates for mitochondrial respiration but it also depends on several compounds synthesized by mitochondria. In the dark, mitochondria are the main source of ATP for cellular processes, including those in the chloroplasts. Moreover, the mitochondrial ATP maintains the proton gradient across the thylakoid membrane, thus protecting the chloroplast from photo inhibition upon reillumination [27]. In the light, mitochondria provide the chloroplast with carbon skeletons for NH 4 + assimilation [28], while ATP supports various biosynthetic reactions, including the repair and recovery of photosystem II (PSII). The tobacco (Nicotianata bacum) CMSII mutant lacks the major mitochondrial NADH dehydrogenase (Complex I) and exhibits a decrease in the rate of photosynthesis, notably during dark-light transitions or when carbon fixation and photo respiration are simultaneously active [29,30]. It is possible that, due to a lack of these photosynthetic proteins, in the CMS line of pepper, a decrease in the rate of photosynthesis might lead to pollen abortion.

Proteins Involved in Antioxidative Reactions
PPO is a nuclear encoded gene whose enzyme attaches to the thylakoids of chloroplasts [31]. It can be released during ripening or senescence and associated with browning in fruits and vegetables like bananas, avocados, cocoa, and tea. In our study, two spots (A01 and A02) were identified as polyphenol oxidase (PPO), which showed much higher concentrations in sterile anthers/buds. As shown in Figure 3A, anthers in male sterile become gradually brown and senescence. Therefore, we would expect to find an abundance of this enzyme in sterile anthers. These results are consistent with premature cell death and point to the A-gene's inhibition of the cox II enzyme, which is known to cause premature tapetal programed cell death (PCD). The up-regulation of PPO in sterile anthers is a symptom of CMS. This result supports the inhibition of the cox II enzyme in our A-B molecular model.

Proteins Involved in Amino Acid Metabolism
One of the most notable differences in NA3 and NB3 anther/buds gels was the location of glutamine synthetase protein (GS, B09, B10, and B12). As these enzymes encoded by nuclear genes, they are not candidates for causes of cytoplasmic male sterility. In addition, they are targeted to chloroplast and cystosol, so the multiple spots were detected in the gel [32].
Other anther proteins with higher spot volumes in the maintainer line might be involved in signal transduction such as adenosine kinase isoform 1T-like protein (B11) and annexin Cap32 (B16).
Annexins interact in a calcium-dependent manner with membrane phospholipids, and have been proposed to be involved in a variety of cellular processes. In cabbage (Brassica oleracea L. var. capitata L.), BoAnnexin2 gene plays an important role in pollen germination [34]. Allene oxide cyclase (AOC) catalyzes the stereospecific cyclization of unstable alleneoxideto cis-(+) 12-oxo-phytodienoicacid (OPDA), which is a crucial step in the biosynthsis of jasmonic acid (JA). The AOC protein accumulated in ovules and in parenchymatic cells of vascular bundles of flower stalks [35]. JA and cis-(+)-OPDA are not only important signaling molecules in the coordination of plant response, but also play an important role in the regulation of developmental processes [36]. In Arabidopsis thaliana, JA is important for the release of pollen and elongation of filaments. JA-deficient or JA-insensitive plants are male sterile. The similar results had been reported in the moss Physcomitrella patens [37].

Plant Materials and Anther Collection
The pepper CMS line NA3 and its maintainer line NB3 with the same nuclear background were used in this study. NA3 with stable male sterile was generated by backcrossing the male sterile line of 8907A [38] as the genetic background of sterile cytoplasm, and the common chili pepper inbred line, North3, as the male parent. The backcrossing has been over 15 generations. Therefore, any differences in anther proteins between these two lines are attributed to their cytoplasms. Plants were field grown without nutrient and moisture stress at the Zengchen Experimental Station, South China Agricultural University, Guangzhou, China (23°8N, 113°17E). Flower buds were cut from NA3 and NB3 at five temporal stages that targeted phase 1 of anther development ( Figure 3F). Bulk samples were collected from 30 individuals. The whole flower buds of stage I and II, and anthers of other stages, were excised and immediately frozen in liquid nitrogen. Following, the samples are stored at −80 °C until protein extraction.

Protein Extraction and 2-DE Electrophoresis
Protein extractions were performed using a trichlor-oacetic acid (TCA)-acetone protocol with some modifications [13]. The samples were ground to a fine powder in liquid nitrogen using a pestle and mortar. The samples were transferred to a centrifuge tube and immediately suspended in 10 volumes of TCA/acetone solution (10% w/v TCA and 0.07% w/v 2-mercaptoethanol in acetone) for over-night protein precipitation at −20 °C to remove secondary metabolites. Samples were centrifuged at 4 °C for 15 min at 13,000× g, and precipitated proteins were washed twice with three volumes of cold (−20 °C) 100% acetone, incubated at −20 °C for 2 h, and then centrifuged. The pellet was air-dried and resuspended in solubilization buffer comprising 7 M urea, 2 M thiourea, 4% CHAPS, and 30 mM Tris, adjusted to pH 8.5. The solution was incubated at 4 °C for 1 h with continuous shaking, and then centrifuged at 4 °C for 15 min at 13,000 × g. The supernatant was collected and stored at −80 °C for further analysis. Protein concentrations were determined by Bradford assay (Bio-Rad, Alfred Nobel Drive Hercules, CA, USA) with BSA as the standard.
Isoelectric focusing (IEF) was carried out using an IPGphor III electrophoresis system (GE Healthcare, Pittsburgh, PA, USA) and 24 cm immobiline dry strips, pH 4-7 (GE Healthcare, Pittsburgh, PA, USA). Nine hundred micrograms of protein samples was loaded during the rehydration step (14 h). IEF was then performed by ramping to 300 V for 30 min, at 700 V for 30 min, and 1500 V for 1.5 h successively, ramping to 9000 V over 1 h, and holding at 9000 V, until a total of 52 kV was reached. Prior to the second dimension separation, the IPG strips were equilibrated in equilibration buffer (6 M urea, 30% (w/v) glycerol, 2% (w/v) SDS, 50 mM Tris-HCl, pH 8.0), first with 1% DTT and then with 2.5% iodoacetamide each for 15 min. The electrophoresed and equilibrated strips were then transferred to 12.5% vertical SDS-PAGE gels for the second dimension electrophoresis using an EttanDALT six Large Vertical System (GE Healthcare, Pittsburgh, PA, USA). SDS-PAGE was run at 2 W/gel for 45 min and then 17 W/gel until the bromphenol blue dye front reached the gel end. After electrophoresis, gels were stained with Coomassie Brilliant Blue (CBB) G-250 and silver nitrate.

Gel Imaging and Analysis
For image acquisition, the stained gels were scanned using ImageScanner III (GE Healthcare, Pittsburgh, PA, USA) at a resolution of 300 dpi and 16-bit grayscale pixel depth. The comparison between NA3 and NB3 gels was made with Image Master 2D Platinum software, Version 5.0 (GE Healthcare, Pittsburgh, PA, USA) as described in the user manual. Background subtraction and normalization were fully automatic. Minimal manual editing was performed to correct mismatched and unmatched spots between gels. The average vol.% values were calculated from three technical replicates to represent the final vol.% values of each biological replicate. NA3 and NB3 spots with 1.5-fold change or more and p < 0.05 were considered to be differentially expressed.

Trypsindigestion, Mass Spectrometry, and Protein Identification
The selected spots were excised from 2-D gels, and pellets were washed with 100 μL destaining solution (50 mM NH 4 HCO 3 in 50% ACN) in micro centrifuge tubes. This step was repeated for two to three times until the pellets were colorless. The gel pieces were washed by Milli-Q water and lyophilized, then rehydrated in digestion buffer containing 20 μg/mL of sequencing grade modified trypsin (Promega, Madison, WI, USA) in 25 mM NH 4 HCO 3 at 37 °C overnight. After brief centrifugation, the peptides were collected from supernatant, and left gel pieces were further sonicated for 10 min in 5 μL of 0.1% TFA (trifluoroacetic acid) in 50% ACN to collect the remaining peptides. The peptides from one protein spot were combined.
Mass spectrometry analysis was carried out using Autoflex speed™ MALDI-TOF-TOF (Bruker Daltonics, Bremen, Germany) tandem mass spectrometer. Parameters were set as follows: laser intensity 4500, mass range 700-3, 200 kDa, acceleration voltage 20 kV, repeat rate 200 Hz. Peptide mass finger printings (PMFs) obtained from MAIDI-TOF/MS were used to search the NCBInr protein database (http://www.ncbi.nlm.nih.gov/) and SwissProt protein databases using the MASCOT program (http://www.matrixscience.com). The following search parameters were used: monoisotopic peptide mass; 800-4000 Da; one missed cleavage per peptide; enzyme, trypsin; taxonomy, green plants; precursorion mass tolerance, 50 ppm; MS/MS fragmention mass tolerance, 0.5 Da; variable modifications, carbamidom ethylation for cysteine, and oxidation for methionine were allowed. Known contaminant ions corresponding to trypsin and keratins were excluded from the peak lists before database searching. Top six hits for each protein search were reported. Only proteins with Mascot proteins cores (based on both MS and MS/MS spectra) of 80 or more (p < 0.05), and a minimum of two matched peptides were considered to be positively identified.

PPO and EFE Activity Assay
Polyphenol oxidase (PPO) activity was determined according to the method of Leja et al. [39]. One milliliter reaction mixture contained 20 µL enzyme extract and 10 mmol/L phosphate buffer (pH 7.0). Each sample was aerated for 2 min in a small test tube followed by the addition of catechol as the substrate at a final concentration of 20 mmol/L. PPO activity was presented as the change in one unit of absorbance at 420 nm per minute per gram fresh weight of sample.
Ethylene forming enzyme (EFE) activity was determined according to the method of Wei and Kenji [40]. Zero point three grams of the buds or anthers was sealed in a vial with 2.5 mL reaction buffer (pH 7.2), 1 mmol/L ACC, 100 mmol/L MOPS. After 20 min of continuous shaking at 35 °C , 1 mL gas sample of the vial atmosphere was withdrawn by syringe and its ethylene content was determined by gas chromatography (GC-17A, SHIMADZU, Nakagyo-ku, Kyoto, Japan).

Conclusions
In summary, a comparative proteomic approach was used to identify differentially expressed proteins in developing anthers of CMS and its maintainer of Capsicum annuum L. Two-dimensional gel electrophoresis (2-DE) was used to identify 27 proteins that showed more than 1.5-fold volume ratio difference in the two lines. The results showed the breakdown of pollen development in the CMS line was associated with differential expression of certain proteins. We speculated that male sterility in NA3 might be related to energy metabolism turbulence, excessive ethylene synthesis or suffocation of starch synthesis. Therefore, this study lays the foundation for future investigation of gene function related to pollen development and cytoplasmic male sterility in pepper.