Skip to main content
  • Research article
  • Open access
  • Published:

Comparative proteomic analysis reveals a dynamic pollen plasma membrane protein map and the membrane landscape of receptor-like kinases and transporters important for pollen tube growth and interaction with pistils in rice



The coordination of pollen tube (PT) growth, guidance and timely growth arrest and rupture mediated by PT-pistil interaction is crucial for the PT to transport sperm cells into ovules for double fertilization. The plasma membrane (PM) represents an important interface for cell–cell interaction, and PM proteins of PTs are pioneers for mediating PT integrity and interaction with pistils. Thus, understanding the mechanisms underlying these events is important for proteomics.


Using the efficient aqueous polymer two-phase system and alkali buffer treatment, we prepared high-purity PM from mature and germinated pollen of rice. We used iTRAQ quantitative proteomic methods and identified 1,121 PM-related proteins (PMrPs) (matched to 899 loci); 192 showed differential expression in the two pollen cell types, 119 increased and 73 decreased in abundance during germination. The PMrP and differentially expressed PMrP sets all showed a functional skew toward signal transduction, transporters, wall remodeling/metabolism and membrane trafficking. Their genomic loci had strong chromosome bias. We found 37 receptor-like kinases (RLKs) from 8 kinase subfamilies and 209 transporters involved in flux of diversified ions and metabolites. In combination with the rice pollen transcriptome data, we revealed that in general, the protein expression of these PMrPs disagreed with their mRNA expression, with inconsistent mRNA expression for 74% of differentially expressed PMrPs.


This study identified genome-wide pollen PMrPs, and provided insights into the membrane profile of receptor-like kinases and transporters important for pollen tube growth and interaction with pistils. These pollen PMrPs and their mRNAs showed discordant expression. This work provides resource and knowledge to further dissect mechanisms by which pollen or the PT controls PMrP abundance and monitors interactions and ion and metabolite exchanges with female cells in rice.


Pollen, consisting of the large vegetative cell (VC) and immobile male gametes (sperm cells) enclosed in the VC, represents an innovative phenotype of plants during evolution. Pollen is tolerant to dehydration and can transport the male gamete over a long distance with the help of wind and/or animals, important forces driving the distribution of plants in the land. Once landing on the stigma of pistils, pollen hydrates, germinates and gives rise to a tip-growing pollen tube (PT). The tube further journeys within the pistil and finally arrives in the synergids of embryo sacs, where it arrests growth and ruptures to release sperm cells for double fertilization. This process requires maintenance of PT integrity and PT–pistil interaction to guide the PT toward the embryo sac and regulate the timely growth arrest of the PT and rupture for immobile gamete release, thereby precisely guaranteeing the one PT to one ovule relationship [1].

The plasma membrane (PM) represents the semi-permeable barrier for selective flux of ion and metabolites across the membrane and an important interface for cell–cell interaction. PM proteins, especially receptor-like kinases (RLKs) and transporters/channels of PTs are the pioneers for mediating PT integrity and interaction with pistils [1]. In Arabidopsis, ANX1 and ANX2, members of the plant-specific RLKs of the Catharanthus roseus RLK1-like (CrRLK1L) subfamily [2], are localized in the PT PM and redundantly regulate PT growth and integrity [3, 4]. PTs in anx1anx2 double mutants show precocious rupture. Overexpression of ANXs caused PT growth inhibition [5]. RUPO, the rice CrRLK1L member, controls PT growth and integrity by interacting with K+ transporters, so a novel RLK signaling pathway mediates K+ homeostasis is required for PT growth and integrity [6]. Studies have also revealed the involvement of the leucine-rich repeat RLK (LRR-RLK) subfamily in PT growth and guidance. Tomato pollen-specific LRR-RLKs LePRK1 and LePRK2, expressed specifically in pollen, regulate PT growth and can bind STIGMA-SPECIFIC PROTEIN1 (STIG1), a small cysteine-rich protein from the pistil [7]. The Arabidopsis LRR-RLKs MDIS1, MIK1 and MIK2 form heteromers, and the complex functions as a receptor of LURE1, the defensin-like cysteine-rich peptide from synergids, to regulate PT guidance and perception [8]. Another Arabidopsis LRR-RLK, PRK6, was identified as a LURE1 receptor and functions in guiding PT tip growth [9]. Furthermore, several other female factors were identified to be involved in PT growth and/or guidance. Tobacco transmitting tissue-specific (TTS), a pistil transmitting tissue-specific arabinogalactan protein, may have roles in guiding PT growth [10]. Chemocyanin, a stigma-expressed small cell wall protein in lily, induces PT chemotropism [11]. ZmEA1, the small protein exclusively expressed in maize egg apparatus, helps guide PT growth in the short range [12]. Thus, additional RLKs are involved in sensing these female factors.

Ion fluxes across the PM and the ion gradient in PTs are well known and are crucial for PT growth [13]. In Arabidopsis, AtACA9, a Ca2+ATPase for Ca2+ efflux is required for PT growth and PT–synergid contact [14]. Arabidopsis cyclic nucleotide gated Ca2+ channel 7 (AtCNGC7) and 8 for Ca2+ influx redundantly regulate the initiation of PT tip growth [15]. AtCNGC18 is essential for PT directional growth in vitro [16]. Disruption of SPIK, an inward K+ channel in Arabidopsis, strongly reduced K+ influx, which resulted in impaired pollen germination and PT growth [17]. PTs lacking the cation/proton exchangers CHX21 and CHX23 grew down in the transmitting tract and failed to turn to the ovule [18]. The maize ZmES4, the synergid-expressed defensin-like cysteine-rich protein, has roles in opening the PT PM-localized K+ influx channel KZM1, which led to excessive influx of K+, thereby causing PT rupture [19]. Moreover, the rice receptor-like kinase RUPO–K+ transporter signaling pathway has been revealed in PTs [6].

Despite these promising findings, our knowledge of RLKs and transporters/channels (hereafter called transporters) that function in PT growth and interaction with pistils is limited, and a detailed understanding of these components at the omic-wide level and proteomic characteristics of pollen and PT PM is lacking. A systematic knowledge of RLKs and transporters in the PM is crucial for an in-depth understanding of the mechanisms underlying PT growth and interaction with pistils.

Here, we prepared PMs from mature pollen grains (MPGs) and germinated pollen grains (GPGs) and dissected PM proteins by using iTRAQ quantitative proteomics. We identified 1,121 PM-related proteins (PMrPs) (matched to 899 loci), with 192 differentially expressed during pollen germination, and revealed 37 RLKs and 209 transporters in the proteome. All PMrPs and differentially expressed PMrPs featured signal transduction, transporters, wall remodeling/metabolism and membrane trafficking functions. Further comparison of proteomic and transcriptomic data revealed that PMrPs are in general discordant with their mRNA levels, with inconsistent mRNA profiles for 74% of differentially expressed PMrPs. These results provide insights into the proteomic characteristics of pollen PM and the profile of RLKs and transporters in the membrane.


Pollen collection and in vitro germination

Rice cultivar Zhonghua 10 (Oryza sativa L. ssp. japonica) was planted under natural conditions in Beijing. Mature pollen grains (MPGs) were collected at anthesis stage by using a modified vacuum cleaner outfitted with nylon meshes. For germination experiments, fresh collected MPGs were transferred into liquid germination medium (40 mg/L H3BO3, 3 mM Ca (NO3)2 · 4H2O, 3 mg/L VB1, 10% PEG4000, 250 mM sucrose) immediately and cultured with gentle shaking at room temperature (~30 °C) for about 15 min. Under this condition, more than 90% of MPGs synchronously germinated to generate polar-growing PTs. Germinated pollen grains (GPGs) were collected by centrifugation at 1000 × g at 4 °C for 5 min. All collected MPGs and GPGs were used immediately or stored at −80 °C.

Plasma membrane preparation

MPGs and GPGs were homogenized in extracting buffer (250 mM sucrose, protease inhibitor cocktail, 1 mM EDTA, 1 mM DTT, 1 mM PMSF, and 50 mM MOPS/KOH, pH 7.8) by use of the high-speed bench top homogenizer FastPrep-24 (MP Biomedicals, USA). The homogenate was differentially centrifuged at 1,500 × g for 5 min, 12,000 × g for 20 min and then 31,000 × g for 15 min to remove cell debris, mitochondria and other organelle contaminants, respectively. The resulting supernatant was centrifuged at 100,000 × g for 1 h with use of BECKMAN Optima L-80XP (70Ti Rotor, Beckman Coulter, USA) to collect pellets (total microsomal vesicles [MSVs]). MSVs were resuspended in PM isolation buffer (250 mM sucrose, 1 mM DTT, 1 mM PMSF, and 5 mM potassium phosphate, pH 7.8) and used to enrich PM vesicles by use of an aqueous polymer two-phase system [20] of 6.5% (w/w) PEG3350 (Sigma), 6.5% (w/w) Dextran T-500 (Pharmacia), 250 mM sucrose, 5 mM KCl, 1 mM DTT, and 5 mM potassium phosphate, pH 7.8. After enrichment, the collected upper phase was diluted more than three-fold with dilution buffer (250 mM sucrose, 1 mM DTT, 1 mM PMSF, and 50 mM MOPS/KOH, pH 7.8) and centrifuged at 200,000 × g for 1 h to collect PM vesicles. PM vesicles were washed with the dilution buffer, then treated with 100 mM sodium carbonate (pH 11.5) to remove soluble proteins associated with the PM vesicles as described [21]. All procedures were carried out at 4 °C. Protein concentration was measured by Bradford assay with bovine serum albumin (BSA) as a standard.

SDS-PAGE and western blot analysis

Proteins were separated by 10% SDS-PAGE. For Western blot analysis, proteins in gels were electrotransferred onto a PVDF membrane (Pierce, USA) with 25 mM Tris, 192 mM glycine and 20% methanol and incubated with the primary rabbit antibodies for PM H+-ATPase (PMA2) from Nicotiana plumbaginifolia (1:5000 dilution) [22]; mitochondrial cytochrome oxidase subunit 2 (COX II) (Agrisera no. AS04053A, Sweden, 1:5000 dilution), vacuole ATPase (V-ATPase) (Agrisera no. AS07213, 1:5000 dilution), 40S ribosomal protein S14-1 (Agrisera no. AS09477, 1:3000 dilution), and ras-related protein1 (Sar) (Agrisera no. AS08326, 1:1000 dilution) from Arabidopsis thaliana; histone H1 (LOC_O s04g18090.1, Beijing Protein Innovation, China, 1:1000 dilution), DEAD-box ATP-dependent RNA helicase (eIF4a; LOC_Os02g05330.1, 1:1000 dilution), glyceraldehyde-3-phosphate dehydrogenase (GAPDH; LOC_Os04g40950.1, 1:1000 dilution), and flotillin like protein (Band_7; Beijing B&M Biotech Co., 1:1000 dilution) from Oryza sativa. Optical density of Western blot bands was quantified by using Image-Pro Plus v6.0 (Media Cybernetics, USA).

In-solution digestion, iTRAQ labeling and strong cation exchange fractionation

Protein digestion and iTRAQ labeling were performed according to the iTRAQ reagents chemistry reference guide (iTRAQ Reagents Multiplex kit, AB SCIEX) with a few modifications. Briefly, proteins (100 μg) from purified PM vesicles were supplemented with RapiGest SF surfactant (Waters, USA) at a final concentration of 0.2% (w/v) for denaturation and enzymatic digestion enhancing, then reduced with 10 mM TCEP, pH 8.0 at 56 °C for 1 h followed by alkylation with 50 mM iodoacetamide in the dark (room temperature, 45 min). Pretreated proteins were digested with trypsin at a ratio of 1:50 (w/w) (Roche) at 37 °C for 16 h, and resulting peptides were labeled with iTRAQ reagents. Experimental repeats were designed as follows: experiment 1 was a mixture of 115 tag-labeled MPG and 117 tag-labeled GPG samples; experiment 2 was a mixture of 117 tag-labeled MPG and 115 tag-labeled GPG samples. Experiment 1 and 2 were lyophilized and subjected to strong cation exchange (SCX) fractionation.

SCX fractionation was performed with an AKTA Purifier 10 HPLC system (GE Amersham Biosciences, USA). The lyophilized samples were resuspended in solvent A (5 mM ammonium chloride, 25% [v/v] acetonitrile, pH 3.0) and fractionated with a PolySULFOETHYL A column (2.1 × 200 mm, 5 μm, 300 Å, PolyLC, Columbia, MD, USA) at a flow rate of 200 μl/min through a linear gradient (0-60%, 90 min) of solvent B (500 mM ammonium chloride, 25% [v/v] acetonitrile, pH 3.0) followed by 60-100% solvent B for 10 min, and 100% solvent B for 15 min. Each separated sample was pooled to 16 fractions and lyophilized for LC-MS/MS analysis.

Nano LC-MS/MS analysis

Each SCX faction was reconstituted with 100 μL 0.1% formic acid in water. After the removal of salt on a Vydac C18 SPE cartridge by centrifugation, the desalted fraction was supplemented with 50 μL 60% acetonitrile and dried with a speed-vac. Each dried fraction was solved in 20 μL 0.1% formic acid, and 10 μL of the solution was injected for nanoLC-MS/MS analysis by use of AB SCIEX Triple TOF 5600 MS (Concord, Ontario, Canada) equipped with a splitless Eksigent nano Ultra 2D Plus nanoLC system and a cHiPLC-Nanoflex microchip system (Dublin, CA, USA). The cHiPLC-system used changeable microfluidic traps (200 μm × 5 mm) and analytical columns (75 μm × 150 mm) packed with ChromXP C18 (3 μm, 120 Å) for online separation analysis. Sample loading, trapping and desalting involved 100% of mobile phase A (2% acetonitrile, 0.2% formic acid, 98% water) at a flow rate of 2 μL/min for 10 min. Peptide elution was started with 5% mobile phase B (98% acetonitrile, 0.2% formic acid, 2% water), then the gradient increased linearly to 24% in 70 min at a flow rate of 300 nL/min. The total gradient length was 120 min. MS data acquisition was performed in the information dependent acquisition (IDA) mode. Triple TOF 5600 MS was operated with a resolving power of 30,000 (FWHM) for TOF MS scans. IDA survey scans were acquired in 250 ms with mass range of m/z 350–1250. As many as 30 product ion scans were collected for 100 ms with mass range of m/z 100–1500, if exceeding a threshold of 120 cps (counts/s) and with a charge state of +2 to +5. Dynamic exclusion was set for 18 sec. Collision energies were calculated on-the-fly for all precursor ions by using empirical equations based on mass and charge (Rolling CE on), and the Enhance iTRAQ function was turned on to improve the efficiency of the collision-induced dissociation.

Protein identification and quantitation

The raw data files (*.wiff) generated by Triple TOF 5600 were analyzed by using ProteinPilot 4.0 (revision 460, AB SCIEX), which involved two different algorithms; Paragon and Pro Group.

Paragon is a search engine that uses feature probabilities and sequence temperature values to identify peptides from MS/MS spectra [23]. Database and parameters used for searching were as follows: NCBI Oryza sativa nonredundant database (136,389 protein entries, August 2011); Sample Type-iTRAQ 4plex labeled; Cys Alkylation-Iodoacetamide; Digestion-Trypsin; Instrument-Triple TOF 5600; Quantitate; Bias Correction; Background Correction; Biological modifications. Precursor mass tolerance was 0.05 Da and fragment mass tolerance was 0.1 Da. As part of the Paragon analysis method, false discovery rate analysis was performed by searching the decoy database to assess the rate of inaccurately identified proteins.

Paragon search results were further processed by the Pro Group algorithm to determine the smallest justifiable set of detected proteins. Each detected protein has an unused protein score, a measurement of all the peptide evidence for a protein that is not better explained by a higher ranking protein, and this score is the true indicator of protein confidence. Unused protein scores 2.0, 1.3, 1.0, and 0.47 correspond to peptide confidence 99, 95, 90, and 66% respectively (score 1.3, 95% confidence was threshold of this work), as shown in the ProteinPilot Software Beta Help (AB SCIEX).

Protein quantitative analysis was also performed by use of ProteinPilot 4.0. The software calculates protein expression change ratios between different samples based on the relative intensities of iTRAQ-labeled peptides. Only ratios from the spectra that are distinct to each protein or protein form were used, to eliminate any masking of changes in expression due to peptides shared between proteins. For each protein expression change ratio reported, the program calculates a p-value that indicates the probability of randomly detecting a ratio different from 1. If an expression change ratio is extremely well determined, a real change can be detected even when the ratio is not very different from 1. To obtain a more accurate quantified result, we chose bias correction and background correction when searching the database by Paragon. Specifically, the criteria for determining MPGs and GPGs differentially expressed PM proteins are two experiments expressed, p-value ≤ 0.05 and GPG/MPG ≥ 1.50 or GPG/MPG ≤ 0.67.

In silico analysis

Protein molecular weight (MW) and isoelectric point (pI) were calculated by using the ProParam tool of Expasy ( The matched loci IDs were obtained from the Rice Genome Annotation Project ( Cellular component, biological process, and molecular function for proteins were annotated by using gene ontology (GO, or WoLF PSORT ( Protein transmembrane domains (TMDs) were predicted by using HMMTOP 2.0 ( Modifications related to membrane localization of a protein including glycosylphosphatidyl inositol (GPI) attachment, prenylation, myristoylation, and palmitoylation were predicted by using the big-PI Predictor (, PrePS-Prenylation Prediction Suite (, N-Myristoyltransferase (NMT,, and CSS-Palm 3.0 (, respectively. Protein annotations were comprehensively evaluated by using a combination of NCBI (, RGAP 7 (, and ARAMEMNON 7.0 (

Sequence and phylogenetic analysis

Protein sequences of Arabidopsis were obtained from The Arabidopsis Information Resource and those for rice and other species were obtained from NCBI. Multiple sequence alignments involved use of BioEdit with the Clustal W method. The protein relevance and phylogenetic tree analysis involved use of MEGA4.0 software.


Preparation and purity detection of plasma membranes from pollen

To identify PM proteins, we first prepared PM vesicles from rice MPGs (Fig. 1a) and GPGs (Fig. 1b) by using the aqueous two-phase partition system followed by a high pH carbonate buffer wash (Fig. 1c). Furthermore, we evaluated the purity of the purified PMs by Western blot analysis with antibodies for PM-specific P-type H+-ATPase PMA2, mitochondrial COX II, vacuole-specific V-ATPase, ribosome protein S14-1, and nuclear protein histone H1. PMA2 was detected as two different isoforms in rice pollen: one showed increased abundance from the entire cell (EC) lysate to carbonate buffer-washed PM (CPM) and one was almost undetectable in EC lysates but was most abundant in CPM (Fig. 2a). All other marker proteins were abundant in EC lysates and/or microsome vesicles and almost undetectable in CPM. These results indicated high purity of the purified PMs with high pH carbonate buffer treatment.

Fig. 1
figure 1

Morphology of rice mature pollen grains (MPGs) and germinated pollen grains (GPGs) and the overview of workflow. a Highly dehydrated MPGs. b in vitro germinated GPGs. c Overview of the experimental scheme. Bar = 50 μm

Fig. 2
figure 2

Plasma membrane (PM) purity verification and PM-related proteins evaluation. a Western blot examination of plasma membrane enrichment. Proteins from entire cell (EC) lysates, microsomal vesicles (MSV), PM (plasma membrane vesicles) and carbonate-washed PM (CPM) were separated by 10% SDS-PAGE, transferred to PVDF membranes and detected with antibodies for PMA2, a PM marker; COXII, a mitochondrial marker; V-ATPase, a vacuole marker; S14-1, a ribosome marker; or H1, a nucleus marker. For detection of PMA2, COXII and V-ATPase, 5 μg protein was loaded per lane; for detection of S14-1 and H1, 10 μg protein was loaded per lane. b Summary of proteins with transmembrane domain (TMD) predicted with use of HMMTOP 2.0. c Venn diagram depicting the distribution of proteins with different lipid modifications. GPI, glycosylphosphatidylinositol anchor; Pre, prenylation site; Myr, myristoylation site; Pal, palmitoylation site. d Proteins predicted to have a TMD or post-translational modification (PTM) or both. e Protein subcellular locations annotated by Gene Ontology or WoLF PSORT showed that 1,121 of the 1,797 proteins with a TMD and/or PTM had PM location information

Protein identification and PM-related protein evaluation

To identify PM-related proteins and determine differences in PM proteomes between MPGs and GPGs, we digested PM proteins with trypsin by using RapiGest SF and obtained well-digested peptides (Additional file 1) for iTRAQ labeling. iTRAQ-labeled peptides were fractioned into 16 fractions by SCX chromatography to reduce sample complexity and increase the identification efficiency of low-abundant PM proteins. The UV-Time curves showed high reproducibility between biological repeated experiments (Additional file 2).

Proteins in these fractions were analyzed by reverse-phase high-performance liquid chromatography coupled with MS/MS. Under the criterion of false discovery rate (FDR) < 1%, unused ≥ 1.3 and two or more unique peptides matched, we identified 1,474 proteins with FDR 0.07% in experiment 1, and 1,284 proteins with FDR 0.08% in experiment 2 (Additional file 3). In total, 1,979 proteins were identified (matched to 1,631 loci), of which 779 were shared in both experiments, with 695 only in experiment 1 and 505 in experiment 2 (Additional file 4).

We analyzed PM location information for the identified 1,979 protein according to annotation of transmembrane domain (TMD), posttranslational modification (PTM), and subcellular location. The analysis of TMDs showed that 1,137/1,979 proteins (57.5%) had at least one TMD (1–4 TMDs for 922 proteins, 5–10 for 135, 10–20 for 80) (Fig. 2b). PTMs including GPI-anchor, prenylation, myristoylation, and palmitoylation are important in mediating PM localization of proteins and in regulating stability and function of proteins [24]. Among the 1,979 proteins, 20 (1.0%) had potential GPI-anchor motifs, 28 (1.4%) potential prenylationsites, 42 (2.1%) potential myristoylation sites and 1,610 (81.4%) potential palmitoylation sites (Fig. 2c and Additional file 4). Studies have revealed that palmitoylation plays a key role in protein sorting [25]. Ultimately, our analysis revealed 1,618/1,979 proteins (82%) with membrane-anchoring motifs (Fig. 2c). In total, 1,797/1,979 proteins (91%) were predicted to have TMDs or one or more membrane-anchoring motifs, or both, which represented potential PM proteins in pollen (Fig. 2d). Furthermore, we obtained protein subcellular location information annotated by gene ontology or WoRF PSORT; among the 1,797 proteins, 1,121 showed PM localization, with the remaining 676 having information for localization in cytoplasm, plastid, mitochondria, nuclei, ribosome, peroxisome or vacuole (Fig. 2e and Additional file 4). Therefore, we considered the 676 proteins as possible contaminants, although they or some also possibly localized in the PM. We finally revealed 1,121 proteins (matched to 899 loci) showing a strong relationship with the PM (Fig. 2e and Additional file 5) and used them for the following analysis.

Functional categories of PM-related proteins

We collected function information for the 1,121 PM-related proteins (PMrPs) from NCBI, RGAP, ARAMEMNON, and GO databases and analyzed by functional categories to gain insight into the biological process occurring in or around pollen PM. These proteins could be organized into 9 categories (Additional file 5). Overall, 67% of these proteins were in 4 groups: signal transduction (25%), transporters (19%), membrane trafficking (11%) and wall remodeling and metabolism (12%); 24% were associated with 4 other groups (cytoskeleton dynamic, protein destination, stress response, and other process); and the function of the remaining 9% was unknown (Fig. 3a). The preferred distribution of these PMrPs in signal transduction, transporters, membrane trafficking and wall dynamics and metabolism is consistent with important roles of these processes in PT growth and cell–cell interaction during fertilization.

Fig. 3
figure 3

Functional categories and network of PM-related proteins (PMrPs) and differentially expressed PMrPs. a Proportion of 1,121 PMrPs and 192 differentially expressed PMrPs in each functional category. b Number of abundance-increased and -decreased PMrPs in germinated pollen in each functional category. c PMrPs and differentially expressed PMrPs were preferentially encoded by chromosomes (chr) 1, 2 and 3. The proteins biased in chr 3 were significantly skewed toward signal transduction as well as wall remodeling and metabolism, and membrane trafficking. d a network of PMrPs. Proteins in red are increased and in blue are decreased in levels in GPGs. Detailed information for the proteins’ abbreviation were listed in Additional file 5A

Quantitative difference between MPG and GPG PMrPs

Our iTRAQ analysis showed high protein quantitative efficiency. Overall, 1,381/1,474 proteins (94%) identified in experiment 1 were quantified and 1,147/1,284 (89%) in experiment 2 were quantified (Additional file 3). Experiments 1 and 2 shared 728 proteins with quantified information. The Pearson correlation coefficient for the two independent iTRAQ experiments was 0.877 (Additional file 6), which indicates well-quantified reproducibility (Additional file 7).

Of the 1,121 PMrPs, 446 (matched to 442 loci) were reproducibly identified in the two independent experiments. Using the cut-off of fold change in expression (GPG/MPG) ≥ 1.5 or ≤0.67 and p-value ≤0.05, we revealed 192 PMrPs with significantly changed expression (matched to 192 loci) between MPGs and GPGs, with 119 abundance-increased and 73 abundance-decreased in GPGs (Additional file 8). Among these changed proteins, proteins involved in signal transduction were overrepresented (35%), with a high proportion of proteins related to transporters (17%), wall remodeling and metabolism (17%) and membrane trafficking (13%); the remaining were implicated in cytoskeleton dynamics (8%), protein destination (7%), stress response (2%) and other processes (2%) (Fig. 3a). Most of the proteins implicated in signal transduction (38/67), wall remodeling and metabolism (25/32), cytoskeleton dynamic (12/15), and protein destination (12/14) showed increased expression in GPGs; the number of abundance-increased and -decreased proteins in transporters (16 vs 17) and membrane trafficking (12 vs 13) seemed similar (Fig. 3b).

The rice genome is estimated to encode 50,000 to 60,000 genes. Annotations from RGAP 7.0 ( showed 12% of these genes distributed in chromosome (chr) 1; 10% each in chr 2, 3, and 4; 8% each in chr 5, 6, 7 and 8; 6% each in chr 9 and 10, and 7% each in chr 11 and 12. The genomic loci of pollen PMrPs and the differentially expressed PMrPs showed a significant chromosome bias. They were enriched on chr 1, 2 and 3, but not on the other chromosomes (Fig. 3c). These enriched proteins on the 3 chromosomes were significantly represented by signal transduction proteins (Fig. 3c, Additional file 9).

To validate the expression patterns of proteins detected by the iTRAQ proteomic approach, we used Western blot analysis to examine the expression of 4 proteins that were increased (eIF4a and GAPDH) and decreased (Sar) in levels in GPGs or had no change in level (Band_7) between MPGs and GPGs in iTRAQ data (Fig. 4a). Signal intensity values of Western blot bands were used for quantity analysis (Additional file 10). The expression patterns for all 4 proteins were consistent with the detection by iTRAQ analysis, with a correlation coefficient of 0.9983 (Fig. 4b), thus indicating the reliability of the iTRAQ proteomic results.

Fig. 4
figure 4

Western blot evaluation of the iTRAQ quantitative information. a Western blot analysis of the expression patterns of eIF4a (the eukaryotic initiation factor-4a, gi|115444197), GAPDH (glyceraldehyde-3-phosphate dehydrogenase, gi|115459078), Band_7 (flotillin like protein, gi|48716660) and Sar (ras-related protein, gi|115436368) in MPGs and GPGs. b Proteins examined by Western blot analysis and iTRAQ show similar tendency in expression pattern

RLKs in pollen plasma membrane

Our analysis revealed 277 PM-related components of signal transduction. These proteins are implicated in diverse signaling pathways, such as Ca2+, phospholipid, auxin, abscisic acid (ABA), gibberellic acid (GA) and phosphorylation cascades (Fig. 3d, Additional file 5). RLKs are important regulators of diverse cellular and developmental processes, such as pollen-stigma recognition [26], cell elongation [27], PT guidance [28] and rupture [4]. With a blast-based search of these predicted rice RLK sequences [29], we found 916 unique RLKs annotated in RGAP 7.0; 99 RLKs were in the reported MPG/GPG transcriptome and could be assigned to 11 subfamilies and 1 unassigned group (Additional file 11). Furthermore, 37 RLKs in our pollen PMrP dataset were assigned to the subfamilies CrRLK1L (n = 3), extensin (n = 1), leucine-rich repeat (LRR-RLK, n = 15), proline-rich extensin like receptor kinase (PERK, n = 2), receptor-like cytoplasmic kinases (RLCK, n = 13), S-domain (SD-RLK, n = 1) and unknown receptor kinase (URK, n = 1) with 1 not assigned (Table 1). These RLKs had diverse functions in different cellular processes (Additional file 11). Among the 37 pollen PM-localized RLKs, one CrRLK1L and 2 RLCKs were increased, and 2 LRR-RLKs and one RLCK were decreased in abundance in GPGs (Table 1).

Table 1 Thirty seven receptor-like kinases (RLKs) identified in plasma membrane of mature pollen grains/germinated pollen grains (MPGs/GPGs) in rice

Transporters in pollen plasma membrane

To understand the mechanisms underlying ion and metabolite flux across the rice PM, we systemically identified transporters in the pollen PM proteome. The transporter classification (TC) system was used to build a transporter classification database (TCDB, with about 10,000 representative and putative non-redundant transporters. By blast searching this database, we identified 209 transporters (matched to 161 loci) in the rice pollen PM proteome; these transporters involved 33 families (Table 2, Additional file 12). These transporters are involved in exchanges and flux across the PM of diverse inorganic ion and metabolites such as Ca2+ (n = 12), H+ (n = 38), K+ (n = 22), Cl (n = 6), Mg2+ (n = 5), sugar (n = 24), phospholipids (n = 3), amino acid/oligopeptide (n = 16), phosphate (n = 7) and sulfate (n = 3) (Table 2). In total, 34 transporters showed changed abundance during pollen germination. Abundance-increased transporters were the H+ transporters gi|218184289, gi|194033213, gi|194033219 and gi|218199814), one K+ transporter (gi|125533127), 3 sugar transporters (gi|222636644, gi|115478530 and gi|108706417), one ABC transporter (gi|218188091), 2 phospholipid transporters (gi|40253457 and gi|53793271), 2 oligopeptide transporters (gi|215697740 and gi|90265689) and 3 other transporters (gi|338817657, gi|10140720 and gi|38567827), and abundance-decreased transporters were another 7 H+ transporters (gi|125597623, gi|115469362, gi|297597907, gi|115451943, gi|115444549, gi|115437984 and gi|115465801), 3 K+ transporters (gi|115462953, gi|15128390 and gi|297722665), one Mg2+ transporter (gi|115454637), 2 sugar transporters (gi|222622219 and gi|115434360), 3 ABC transporters (gi|115485837, gi|115477865 and gi|218198932) and 2 other transporters (gi|90399194 and gi|75253347) (Additional file 12).

Table 2 Two hundred nine transporters identified in plasma membrane of MPGs/GPGs in rice

Comparison between rice pollen PM proteome and transcriptome

To evaluate the possible relation of PMrPs and their transcripts, we retrieved previously reported data for transcripts expressed in MPGs and GPGs [30]. The analysis involved 5,939 transcripts detected in MPGs and 5,945 in GPGs, for a total of 7,161 unique transcripts expressed in MPGs or/and GPGs.

Locus number comparison showed 525/899 pollen PMrPs (58.4%) with corresponding transcripts; 317 of these were pollen-preferential (cutoff at Ratio ≥ 2.0, Ratio = MAX (median (MPGs), median (GPGs))/MAX (callus cells1-3, roots1-3, leaves1-3)) (Additional file 13) and were mainly involved in signal transduction, wall remodeling and metabolism and transporters (Fig. 5a). Unexpectedly, 374/899 of the pollen PMrPs (41.6%) had no corresponding transcripts in the dataset of 7,161 transcripts. Their transcripts may be short-lived or extremely low-abundant or the protein was synthesized at early stages of pollen development and deposited for late requirement [31]. We found that 60% of transcripts for encoding pollen PMrPs were pollen-preferential, with only 44% of total MPG/GPG-expressed transcripts being pollen-preferential (Fig. 5b and c), which suggests that PMrP-encoding genes have high organ- or cell-specific expression.

Fig. 5
figure 5

Comparison of pollen PM protein and transcript levels. a 525 of the PMrPs have corresponding transcripts, and 317 of these transcripts (60%) are pollen-preferential and are mainly implicated in signal transduction. b 44% of the pollen mRNA are pollen-preferential. c 60% of the PMrPs transcripts are pollen-preferential. d No significant correlation between PMrP expression and their transcript expression. e Four expression patterns of PMrP–transcript pairs and their distribution in each functional category

Next, we analyzed the correlation of expression profiles between pollen PMrPs and their transcripts. In total, 500 of the 525 pollen PMrP–transcript pairs had quantitative information for both protein and RNA levels (Additional file 14). The expression profiles of PMrPs and their transcripts were not significantly correlated (correlation coefficient −0.1631) (Fig. 5d). When this analysis was restricted to the 192 differentially expressed PMrPs, 130 of which had corresponding transcripts, the 130 protein–transcript pairs were assigned to 4 patterns: (A) both mRNA and protein levels increased in GPGs (25 pairs), (B) both mRNA and protein levels decreased in GPGs (25 pairs), (C) mRNA levels decreased and protein levels increased in GPGs (53 pairs), and (D) mRNA levels increased and protein levels decreased in GPGs (27 pairs) (Fig. 5e, Additional file 14). In total, 74% of these differentially expressed PMrP showed discordant expression with their mRNAs. Thus, the protein expression pattern of PMrPs was not closely related to their mRNA levels omic-wide, which indicates the importance of PM proteomic studies in understanding pollen function.


We prepared high-purity pollen PM by using the aqueous polymer two-phase system and alkali buffer treatment and identified PMrPs and differentially expressed PMrPs between mature and germinated pollen by using LC-MS and iTRAQ-based quantitative proteomics approaches. This study revealed dynamic characteristics of the pollen PM proteome and a large set of RLKs and transporters in the proteome. As well, the expression pattern of PMrPs was in general inconsistent with that of corresponding mRNAs.

Pollen PM-related protein-encoding genes show biased distribution in rice chromosomes

Studies have shown that 26% of genes of the human genome [32] and 25% of genes of the Arabidopsis genome [33] encode membrane proteins. Accordingly, the rice genome, with an estimated 50,000-60,000 genes [34, 35], encodes about 13,000 membrane proteins. However, PM protein-encoding genes possibly represent a small percentage of a genome. Proteomic clues suggest that only 3% genes in a genome encode PM proteins [36], for about 750 PM proteins in Arabidopsis [37] and more than 2,500 PM proteins in humans [38]. Therefore, rice may contain about 1,500 PM proteins. Our study revealed 1,121 PMrPs (matched 899 loci) in rice pollen, which should represent nearly 60% of the predicted PM proteins.

Data analysis showed that annotated genes appear to distribute in rice chromosomes at a similar frequency, with the exception of chr 1, which has a larger number of genes than other chromosomes. However, the distribution frequency of PMrPs and differentially expressed PMrPs was significantly higher than that of genome genes in chr 1, 2 and 3, with the highest biased distribution for differentially expressed PMrPs in chr 3 (Fig. 3c). As well, these enriched PMrPs and differentially expressed PMrPs in these chromosomes showed significantly functional skew toward signal transduction (Additional file 9). Thus, the biased distribution may reflect a mechanism to coordinate these diverse components of a pathway in action, although further studies are needed to address the phenomena and the biological roles of the biased distribution of pollen PMrPs.

Pollen synthesizes a set of PM proteins specialized for pollen function

Pollen grains of most plants are metabolically quiescent when released from anthers and are tolerant of desiccation for spreading over long distances by pollinators. On loading onto the stigma, pollen germinates to generate a fast-growing PT, which journeys within the pistil via turgor-driving growth at the tip. Thus, PT growth and successful fertilization requires active membrane and wall material transportation, wall remodeling, and signal and materials communication across the PM [39].

The pollen proteome has highly represented wall remodeling and metabolism-related proteins, which differs from vegetative cell proteomes [40]. Consistent with previous findings, pollen PM-related wall remodeling and metabolism proteins were preferentially represented in both pollen PMrP and differentially expressed PMrP sets in our study. Such proteins are involved in wall synthesis, loosening and extension. For example, we revealed abundance-increased enzymes for pectin biosynthesis, such as methylesterase (PME) (gi|115452515, gi|115452623, gi|125601598, gi|115486641, and gi|108864650), and abundance-decreased enzymes for pectin degradation, such as glycosyl hydrolases (gi|115443693), PME inhibiter (gi|297721723), and polygalacturonase (gi|115437052 and gi|125555670) (Additional file 8). This finding is compatible with the importance of pectin in PT growth.

Along with the highly represented wall remodeling and metabolism PMrPs, those implicated in signal transduction, transporters and membrane trafficking were overrepresented in both pollen PMrP and differentially expressed PMrP sets. These membrane trafficking-associated PMrPs mainly involve exocytosis and endocytosis pathways, vesicle targeting and vesicle-mediated transport, which sustain the material basis for cell-wall and PM extension and thus support fast PT growth [41]. We revealed a large set of PMrPs for diverse signaling pathways such as auxin, ABA, calcium, phospholipid, and GTPase signaling; mitogen-activated protein kinase phosphorylation cascades; and receptor-like kinase pathways, which suggest that pollen function requires multiple signaling pathways and coordination of these pathways. We also revealed multiple transporters involved in flux of diverse ions and metabolites across the PM, which are consistent with the importance of the ion gradient and dynamics in PT growth and integrity maintenance (for details, see below). Together, these results indicate that these highly represented terms are compatible with the requirement of unique PT growth mechanism and function, thus providing insights into the molecular network for PT growth and function.

Previous proteome and transcriptome results showed mature and germinated pollen had similar protein and mRNA profiles [30, 40]. Thus, germination and early PT growth depend mainly on presynthesized mRNAs and proteins. Yet the impact of mRNA levels on protein expression is unclear. Here, we evaluated the association between variations in mRNA and protein expression levels. The protein expression patterns of PMrPs and their mRNA levels were not significantly correlated, with uncorrelated mRNA expression profiles for 74% of these differentially expressed PMrPs. This finding differs from the empirical conclusion that the correlation between protein and mRNA levels is generally modest [42]. Studies of mammals have revealed that protein levels are heritable molecular phenotypes and evolve under greater evolutionary constraints than mRNA levels [42, 43]. This highly discordant expression between PMrP protein and mRNA levels suggests that pollen has mechanisms to build a PM proteome for fast tip-growth and function, although further studies are needed to understand the mechanisms.

Pollen has multiple types of PM-localized RLKs

Fast-growing PTs within pistils need to constantly monitor the surface tension of the cell wall and PM and communicate with female stimuli to coordinate PT growth, integrity and timely growth arrest and rupture as they enter the ovule, but the molecular mechanisms remain largely unknown. RLKs relay extracellular and cell surface signals for initiating intracellular signaling cascades, thereby representing the central components of the signaling pathways mediating cell–cell interaction [1]. Transcriptome studies revealed 43 and 99 RLK genes expressed in Arabidopsis MPGs and rice GPGs/MPGs [4, 30], respectively, about 7 and 11% of the predicted RLK genes in the respective genomes (Additional file 11). Here, we revealed the protein levels of 37 RLKs from 8 subfamilies in rice pollen, with 6 levels changed with pollen germination (Table 1).

CrRLK1Ls have putative carbohydrate-binding domains, the malectin domain, in their extracellular regions and are considered candidates for sensing changes in the cell surface or external stimuli [44]. In Arabidopsis, THE1, HERK1 and HERK2 expressed in vegetative tissues are required for cell extension [45]. Pollen-expressed ANX1 and ANX2 are redundant for regulating PT integrity [3, 4], whereas synergid-expressed FER regulates PT reception and timely rupture [46]. We revealed three CrRLK1Ls (gi|56783691, gi|115462979 and gi|25553554) in rice pollen, the latter showing increased expression in germinated pollen (Table 1). Sequence alignment showed that they had a variable malectin domain but highly conserved kinase domains with THE1, HERK1, HERK2, ANX1, ANX2 and FER (Additional file 15). gi|56783691 and gi|115462979 are phylogenetically close to ANX1 and ANX2, but gi|25553554 is phylogenetically close to THE1 (Fig. 6). gi|25553554, termed RUPO, regulates PT integrity by interacting with K+ transporters [6]. These results suggest that multiple CrRLK1L signaling pathways orchestrate PT growth and cross-talk between PTs and female cells.

Fig. 6
figure 6

Phylogenetic relationship of the identified CrRLK1Ls in rice pollen and functionally known Arabidopsis CrRLK1Ls. The tree was created by using MEGA 4.0

The LRR-RLK subfamily contains the most members of RLKs and has 384 members in rice (312 by RGAP7.0) [29]. We revealed 15 LRR-RLKs in rice pollen PM (Table 1). To obtain functional features, we analyzed their phylogenetic relations with functionally known pollen-specific LRR-RLKs (PRKs) [47] required for pollen maturation (PiPRK1 from Petunia inflate) [48], pollen–pistil interaction (LePRK1 and LePRK2 from Lycopersicon esculentum) [28], PT elongation (AtPRK2 from Arabidopsis thaliana) [49] and PT targeting (AtPRK6, AtMDIS1, AtMIK1 and AtMIK2 from Arabidopsis thaliana) [8, 9]. Four of the 15 identified LRR-RLKs (gi|49388978, gi|52075918, gi|115477354 and gi|222640883) are in a clade with ZmPRK1 from Zea mays; in evolution, these LRR-RLKs appeared to differentiate from the dicot members AtPRP1-2 and 4–5, LePRK1-2 and PiPRK1 before the monocot–dicot split. gi|115486303 and gi|4680345 are in a clade with LePRK3, AtPRK3 and AtPRK6 and gi|125547150, gi|115480655, gi|125531685, gi|125561357 and gi|218201938 are in a clade with AtMIK1 and 2, which indicates their differentiation after the monocot–dicot split. The remaining 4 are in an independent clade without functionally known dicot homologs (Fig. 7). In fact, recent study revealed the pollen-expressed Arabidopsis LRR-RLKs MDIS1 and MIK acting as a pollen PM-localized, heteromer receptor for sensing egg-expressed attractant LURE1 to control PT reception and guidance [8]. AtPRK6 was found to be another receptor of the species-specific AtLURE1 that can trigger the Rop signaling cascade to regulate PT growth [9]. Together, these findings suggest that multiple male receptors can respond to the same female attractor within the same species.

Fig. 7
figure 7

Phylogenetic relationship of identified LRR-RLKs in rice pollen and functionally known members of this subfamily. This analysis involves 15 LRR-RLKs identified in rice pollen PM and 14 functionally known members from Arabidopsis thaliana (At), Petunia inflate (Pi), Lycopersicon esculentum (Le), and Zea mays (Zm). The tree was made by MEGA 4.0

RLCK is another large subfamily of RLKs that lack the transmembrane domain [29]. Two pollen-preferential RLCKs in Arabidopsis – LIP1 and LIP2 – are anchored to the PM of the PT tip via palmitoylation, and they can perceive the female signal – LURE1 – attraction [50]. SD-RLKs have the extracellular domains homologous to the secreted S-locus glycoproteins (SLGs) that display genotype-specific sequence polymorphisms. S receptor kinases were proposed to be responsible for the allelic specific self-incompatibility response along with SLGs [51]. We identified 13 RLCKs and one SD-RLK, with gi|222635113 and gi|222613313 increased and gi|115436274 decreased in levels in GPGs (Table 1), as candidates for pollen–pistil recognition. In addition, we revealed the RLKs Extensin, PERK and URK subfamilies, although their functions in PT–pistil interaction remain unknown. Together, our results indicate that multiple RLK signaling pathways are required for orchestrating PT growth and cross-talk between PTs and female cells and provide important candidates for further understanding the mechanisms.

Coordination of diverse ion and metabolites across the pollen PM

The gradient and homeostasis of Ca2+, H+, K+ and Cl are critical for pollen germination and PT growth [52], but the molecular mechanisms underlying the gradient and homeostasis are not fully understood. We revealed 209 transporters in rice pollen PM involved in flux across the PM of Ca2+, K+, H+, Cl, water, sugar, amino acids, oligopeptides and phospholipids (Table 2). Tip-fused Ca2+ gradient is the central player orchestrating the cross-talk of diverse signaling pathways and the ion is also important for regulating the differential plasticity of the PT apical and subapical wall [53]. Lines of evidence have shown that in Arabidopsis, transporters for Ca2+ flux were involved in PT growth, PT–synergid contact [14] and PT guidance [16]. Here, we showed 12 Ca2+ transporters in rice pollen PM from the VIC and P-ATPase superfamily (Table 2), which indicates their importance in controlling Ca2+ gradient and homeostasis in pollen.

The increased pH from the apical to the subapical region functions as an important node in regulating diverse activities in PT, such as the wall plasticity, vesicle transport and actin dynamics [54]. The pH gradient is mainly mediated by PM H+-ATPases (PMAs) that promote H+ exflux in the subapical region [55]. In rice pollen PM, we revealed 12 PMAs from the P-ATPase superfamily (Table 2), which highlights their central roles in sustaining PT alkalinity. Which transporters are responsible for H+ influx in the PT apex are unknown, but the nonspecific cation channels are generally considered responsible for apical H+ influx [52]. The rice pollen PM showed a novel cation-transporting ATPase (OsP5, gi|218196773) of the P-type ATPase superfamily, which gives some hints for PT H+ influx. Besides the PMAs, 13 VHAs and 3 H+-PPases were revealed, with the VHA subunit H (gi|218199814) significantly increased and subunit A (gi|115444549) and C (gi|115437984 and gi|115465801) decreased in levels in GPGs; they may also affect PT growth by sustaining cytosolic pH homeostasis [54].

K+ is an essential mineral for pollen germination and PT growth [56]. Studies have demonstrated that K+ channels were required for PT growth [17] and directional growth to the ovule in Arabidopsis [18] and integrity in maize [19]. We found 22 cation or specific K+ transporters (Table 2), which suggests that multiple transporters are involved in maintaining K+ homeostasis and PT turgor pressure. In addition, studies showed that PTs have a tip-fused negative Cl gradient correlated with vesicle movement and close relationship with water homeostasis [57]. Besides the functionally known anion channel SLAH3, which is activated by calcium-dependent protein kinases via phosphorylation for Cl exflux [58], we identified 3 anion channels and 3 Cl channels in rice pollen PM (Table 2), thus supplying the molecular information for insights into the Cl gradient.

Moreover, we revealed 5 Mg2+ transporters and 2 mechanosensitive channels MscS-like in pollen PM (Table 2). MscS-like 8 plays key roles in helping Arabidopsis pollen survive hypo-osmotic shock during rehydration [59]. Sucrose transporters (SUTs) and monosaccharide transporters (MSTs) are suggested to participate in pollen development, germination and PT growth but show high functional redundancy [60, 61]. The rice genome has 5 SUT genes (OsSUT1-5); only OsSUT1 and OsSUT3 transcripts were detectable in pollen [62]. We identified both OsSUT1 (gi|108706417) and 2 OsSUT3 isoforms (gi|222641552 and gi|115481924) in pollen PM (Additional file 12), which agrees with their functions in pollen germination and starch accumulation [62]. The rice genome has 65 predicted MST genes [63], but so far only 7 have been characterized [64]. We identified 15 MSTs from STP, pGlcT, AZT, INT, PLT, and ERD subfamilies in the PM, with 2 (gi|222636644, gi|115478530) increased and 2 (gi|222622219 and gi|115434360) decreased in levels in GPGs (Additional files 8 and 12 ). Several ATP-binding cassette (ABC) transporters were found to function in exine lipid-soluble precursor secretion [65]. ABC transporters consist of 2 nucleotide-binding domains (NBDs) to bind and hydrolyse ATP and 2 transmembrane domains involved in substrate recognition. According to the NBD evolutionary relationship, the transporters can be classified into 8 major subfamilies (A-G and I) [66]. ABCG transporters were well known to be expressed in tapetum, but little was known about whether they were expressed in pollen [67]. We first identified 6 ABCG transporters and also found 23 other ABC transporters in the PM (Additional file 12). We have little knowledge about the substrates and molecular functions of the 29 ABC transporters.

Finally, we identified other transporters in the rice PM (Additional file 12): (1) the pollen-specific boron efflux transporter OsBOR4, implicated in pollen germination and PT elongation [68]; (2) 3 P-type ATPases belonging to the aminophospholipid flipping (ALA) subfamily [69], which may be involved in establishing membrane lipid asymmetry [70]; (3) 9 amino acid transporters and 7 oligopeptide transporters (The types of transporters display different substrate specificities and may help pollen to sustain the rehydration process by supplying compatible solutes [71]); (4) an aquaporin that may adjust the water condition in pollen; and (5) 7 phosphate, 3 sulfate, 8 micronutrient and some other transporters for various substrates, which may provide the essential nutrients to sustain pollen activity under unclear mechanisms. Together, the large set of diversified transporters identified in pollen PM suggests that fine-tuned multiple ion and metabolite homeostasis and coordination of homeostasis are crucial for PT growth and function. Thus, we provide important molecular information for understanding the mechanisms underlying the homeostasis and coordination (Fig. 8).

Fig. 8
figure 8

Diversified RLKs and transporters are fine-tuned and coordinated in regulating pollen-pistil interaction and pollen tubes tip growth


We revealed 1,121 PMrPs and 192 differentially expressed PMrPs during rice pollen germination. Pollen presynthesizes a set of PM proteins for PT growth and interaction with pistils. The pollen PM protein levels are in general discordant with their mRNA levels, and their genes showed biased distribution in rice chromosomes. The pollen PM have a large set of RLKs mediating cell–cell interaction and transporters for flux of diversified ions and metabolites (Fig. 8). This study provides a resource to further dissect the molecular mechanisms by which pollen or the PT controls PMrP abundance and monitors interactions and ion and metabolite exchanges with female cells in rice.



Catharanthus roseus RLK1-like


Germinated pollen grains


Isobaric tags for relative and absolute quantitation


Leucine-rich repeat


Mature pollen grains


Plasma membrane


PM-related proteins


Pollen tube


Receptor-like kinases


Vegetative cell


  1. Higashiyama T, Takeuchi H. The mechanism and key molecules involved in pollen tube guidance. Annu Rev Plant Biol. 2015;66:393–413.

    Article  CAS  PubMed  Google Scholar 

  2. Nissen KS, Willats WG, Malinovsky FG. Understanding CrRLK1L function: cell walls and growth control. Trends Plant Sci. 2016;21:516–27.

    Article  CAS  PubMed  Google Scholar 

  3. Boisson-Dernier A, Roy S, Kritsas K, Grobei MA, Jaciubek M, Schroeder JI, Grossniklaus U. Disruption of the pollen-expressed FERONIA homologs ANXUR1 and ANXUR2 triggers pollen tube discharge. Development. 2009;136:3279–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Miyazaki S, Murata T, Sakurai-Ozato N, Kubo M, Demura T, Fukuda H, Hasebe M. ANXUR1 and 2, sister genes to FERONIA/SIRENE, are male factors for coordinated fertilization. Curr Biol. 2009;19:1327–31.

    Article  CAS  PubMed  Google Scholar 

  5. Boisson-Dernier A, Lituiev DS, Nestorova A, Franck CM, Thirugnanarajah S, Grossniklaus U. ANXUR receptor-like kinases coordinate cell wall integrity with growth at the pollen tube tip via NADPH oxidases. PLoS Biol. 2013;11:e1001719.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Liu L, Zheng C, Kuang B, Wei L, Yan L, Wang T. Receptor-Like Kinase RUPO interacts with potassium transporters to regulate pollen tube growth and integrity in rice. PLoS Genet. 2016;12:e1006085.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Tang W, Kelley D, Ezcurra I, Cotter R, McCormick S. LeSTIG1, an extracellular binding partner for the pollen receptor kinases LePRK1 and LePRK2, promotes pollen tube growth in vitro. Plant J. 2004;39:343–53.

    Article  CAS  PubMed  Google Scholar 

  8. Wang T, Liang L, Xue Y, Jia PF, Chen W, Zhang MX, Wang YC, Li HJ, Yang WC. A receptor heteromer mediates the male perception of female attractants in plants. Nature. 2016;531(7593):241–4.

    Article  CAS  PubMed  Google Scholar 

  9. Takeuchi H, Higashiyama T. Tip-localized receptors control pollen tube growth and LURE sensing in Arabidopsis. Nature. 2016;531:245–8.

    Article  CAS  PubMed  Google Scholar 

  10. Cheung AY, Wang H, Wu HM. A floral transmitting tissue-specific glycoprotein attracts pollen tubes and stimulates their growth. Cell. 1995;82:383–93.

    Article  CAS  PubMed  Google Scholar 

  11. Kim S, Mollet JC, Dong J, Zhang K, Park SY, Lord EM. Chemocyanin, a small basic protein from the lily stigma, induces pollen tube chemotropism. Proc Natl Acad Sci U S A. 2003;100:16125–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Marton ML, Cordts S, Broadhvest J, Dresselhaus T. Micropylar pollen tube guidance by egg apparatus 1 of maize. Science. 2005;307:573–6.

    Article  CAS  PubMed  Google Scholar 

  13. Zonia L. Spatial and temporal integration of signalling networks regulating pollen tube growth. J Exp Bot. 2010;61:1939–57.

    Article  CAS  PubMed  Google Scholar 

  14. Schiott M, Romanowsky SM, Baekgaard L, Jakobsen MK, Palmgren MG, Harper JF. A plant plasma membrane Ca2+ pump is required for normal pollen tube growth and fertilization. Proc Natl Acad Sci U S A. 2004;101:9502–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Tunc-Ozdemir M, Rato C, Brown E, Rogers S, Mooneyham A, Frietsch S, Myers CT, Poulsen LR, Malho R, Harper JF. Cyclic nucleotide gated channels 7 and 8 are essential for male reproductive fertility. PLoS One. 2013;8:e55277.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Frietsch S, Wang YF, Sladek C, Poulsen LR, Romanowsky SM, Schroeder JI, Harper JF. A cyclic nucleotide-gated channel is essential for polarized tip growth of pollen. Proc Natl Acad Sci U S A. 2007;104:14531–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Mouline K, Very AA, Gaymard F, Boucherez J, Pilot G, Devic M, Bouchez D, Thibaud JB, Sentenac H. Pollen tube development and competitive ability are impaired by disruption of a shaker K(+) channel in Arabidopsis. Genes Dev. 2002;16:339–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Lu Y, Chanroj S, Zulkifli L, Johnson MA, Uozumi N, Cheung A, Sze H. Pollen tubes lacking a pair of K+ transporters fail to target ovules in Arabidopsis. Plant Cell. 2011;23:81–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Amien S, Kliwer I, Marton ML, Debener T, Geiger D, Becker D, Dresselhaus T. Defensin-like ZmES4 mediates pollen tube burst in maize via opening of the potassium channel KZM1. PLoS Biol. 2010;8:e1000388.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Schindler J, Nothwang HG. Aqueous polymer two-phase systems: effective tools for plasma membrane proteomics. Proteomics. 2006;6:5409–17.

    Article  CAS  PubMed  Google Scholar 

  21. Fujiki Y, Hubbard AL, Fowler S, Lazarow PB. Isolation of intracellular membranes by means of sodium carbonate treatment: application to endoplasmic reticulum. J Cell Biol. 1982;93:97–102.

    Article  CAS  PubMed  Google Scholar 

  22. Morsomme P, Dambly S, Maudoux O, Boutry M. Single point mutations distributed in 10 soluble and membrane regions of the Nicotiana plumbaginifolia plasma membrane PMA2 H + −ATPase activate the enzyme and modify the structure of the C-terminal region. J Biol Chem. 1998;273:34837–42.

    Article  CAS  PubMed  Google Scholar 

  23. Shilov IV, Seymour SL, Patel AA, Loboda A, Tang WH, Keating SP, Hunter CL, Nuwaysir LM, Schaeffer DA. The paragon algorithm, a next generation search engine that uses sequence temperature values and feature probabilities to identify peptides from tandem mass spectra. Mol Cell Proteomics. 2007;6:1638–55.

    Article  CAS  PubMed  Google Scholar 

  24. Nadolski MJ, Linder ME. Protein lipidation. FEBS J. 2007;274:5202–10.

    Article  CAS  PubMed  Google Scholar 

  25. Greaves J, Chamberlain LH. Palmitoylation-dependent protein sorting. J Cell Biol. 2007;176:249–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Stein JC, Howlett B, Boyes DC, Nasrallah ME, Nasrallah JB. Molecular cloning of a putative receptor protein kinase gene encoded at the self-incompatibility locus of Brassica oleracea. Proc Natl Acad Sci U S A. 1991;88:8816–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Guo H, Li L, Ye H, Yu X, Algreen A, Yin Y. Three related receptor-like kinases are required for optimal cell elongation in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2009;106:7648–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Muschietti J, Eyal Y, McCormick S. Pollen tube localization implies a role in pollen-pistil interactions for the tomato receptor-like protein kinases LePRK1 and LePRK2. Plant Cell. 1998;10:319–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Shiu SH, Karlowski WM, Pan RS, Tzeng YH, Mayer KFX, Li WH. Comparative analysis of the receptor-like kinase family in Arabidopsis and rice. Plant Cell. 2004;16:1220–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Wei LQ, Xu WY, Deng ZY, Su Z, Xue Y, Wang T. Genome-scale analysis and comparison of gene expression profiles in developing and germinated pollen in Oryza sativa. BMC Genomics. 2010;11:338.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Holmes-Davis R, Tanaka CK, Vensel WH, Hurkman WJ, McCormick S. Proteome mapping of mature pollen of Arabidopsis thaliana. Proteomics. 2005;5:4864–84.

    Article  CAS  PubMed  Google Scholar 

  32. Fagerberg L, Jonasson K, von Heijne G, Uhlen M, Berglund L. Prediction of the human membrane proteome. Proteomics. 2010;10:1141–9.

    Article  CAS  PubMed  Google Scholar 

  33. Schwacke R, Schneider A, van der Graaff E, Fischer K, Catoni E, Desimone M, Frommer WB, Flugge UI, Kunze R. ARAMEMNON, a novel database for Arabidopsis integral membrane proteins. Plant Physiol. 2003;131:16–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Goff SA, Ricke D, Lan TH, Presting G, Wang R, Dunn M, Glazebrook J, Sessions A, Oeller P, Varma H, et al. A draft sequence of the rice genome (Oryza sativa L. ssp. japonica). Science. 2002;296:92–100.

    Article  CAS  PubMed  Google Scholar 

  35. Yu J, Hu S, Wang J, Wong GK, Li S, Liu B, Deng Y, Dai L, Zhou Y, Zhang X, et al. A draft sequence of the rice genome (Oryza sativa L. ssp. indica). Science. 2002;296:79–92.

    Article  CAS  PubMed  Google Scholar 

  36. Masson F, Rossignol M. Basic plasticity of protein expression in tobacco leaf plasma-membrane. Plant J. 1995;8:77–85.

    Article  CAS  Google Scholar 

  37. Marmagne A, Rouet MA, Ferro M, Rolland N, Alcon C, Joyard J, Garin J, Barbier-Brygoo H, Ephritikhine G. Identification of new intrinsic proteins in Arabidopsis plasma membrane proteome. Mol Cell Proteomics. 2004;3:675–91.

    Article  CAS  PubMed  Google Scholar 

  38. Natera SH, Ford KL, Cassin AM, Patterson JH, Newbigin EJ, Bacic A. Analysis of the Oryza sativa plasma membrane proteome using combined protein and peptide fractionation approaches in conjunction with mass spectrometry. J Proteome Res. 2008;7:1159–87.

    Article  CAS  PubMed  Google Scholar 

  39. Cole RA, Fowler JE. Polarized growth: maintaining focus on the tip. Curr Opin Plant Biol. 2006;9:579–88.

    Article  CAS  PubMed  Google Scholar 

  40. Dai S, Chen T, Chong K, Xue Y, Liu S, Wang T. Proteomics identification of differentially expressed proteins associated with pollen germination and tube growth reveals characteristics of germinated Oryza sativa pollen. Mol Cell Proteomics. 2007;6:207–30.

    Article  CAS  PubMed  Google Scholar 

  41. McKenna ST, Kunkel JG, Bosch M, Rounds CM, Vidali L, Winship LJ, Hepler PK. Exocytosis precedes and predicts the increase in growth in oscillating pollen tubes. Plant Cell. 2009;21:3026–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Wu L, Candille SI, Choi Y, Xie D, Jiang L, Li-Pook-Than J, Tang H, Snyder M. Variation and genetic control of protein abundance in humans. Nature. 2013;499:79–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Battle A, Khan Z, Wang SH, Mitrano A, Ford MJ, Pritchard JK, Gilad Y. Genomic variation. Impact of regulatory variation from RNA to protein. Science. 2015;347:664–7.

    Article  CAS  PubMed  Google Scholar 

  44. Lindner H, Muller LM, Boisson-Dernier A, Grossniklaus U. CrRLK1L receptor-like kinases: not just another brick in the wall. Curr Opin Plant Biol. 2012;15:659–69.

    Article  CAS  PubMed  Google Scholar 

  45. Guo H, Ye H, Li L, Yin Y. A family of receptor-like kinases are regulated by BES1 and involved in plant growth in Arabidopsis thaliana. Plant Signal Behav. 2009;4:784–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Escobar-Restrepo JM, Huck N, Kessler S, Gagliardini V, Gheyselinck J, Yang WC, Grossniklaus U. The FERONIA receptor-like kinase mediates male–female interactions during pollen tube reception. Science. 2007;317:656–60.

    Article  CAS  PubMed  Google Scholar 

  47. Kim HU, Cotter R, Johnson S, Senda M, Dodds P, Kulikauska R, Tang W, Ezcura I, Herzmark P, McCormick S. New pollen-specific receptor kinases identified in tomato, maize and Arabidopsis: the tomato kinases show overlapping but distinct localization patterns on pollen tubes. Plant Mol Biol. 2002;50:1–16.

    Article  CAS  PubMed  Google Scholar 

  48. Lee HS, Karunanandaa B, McCubbin A, Gilroy S, Kao T. PRK1, a receptor-like kinase of Petunia inflata, is essential for postmeiotic development of pollen. Plant J. 1996;9:613–24.

    Article  CAS  Google Scholar 

  49. Chang F, Gu Y, Ma H, Yang Z. AtPRK2 promotes ROP1 activation via RopGEFs in the control of polarized pollen tube growth. Mol Plant. 2013;6:1187–201.

    Article  CAS  PubMed  Google Scholar 

  50. Liu J, Zhong S, Guo X, Hao L, Wei X, Huang Q, Hou Y, Shi J, Wang C, Gu H, Qu LJ. Membrane-bound RLCKs LIP1 and LIP2 are essential male factors controlling male–female attraction in Arabidopsis. Curr Biol. 2013;23:993–8.

    Article  CAS  PubMed  Google Scholar 

  51. Stein JC, Nasrallah JB. A plant receptor-like gene, the S-locus receptor kinase of Brassica oleracea L., encodes a functional serine/threonine kinase. Plant Physiol. 1993;101:1103–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Holdaway-Clarke TL, Hepler PK. Control of pollen tube growth: role of ion gradients and fluxes. New Phytol. 2003;159:539–63.

    Article  CAS  Google Scholar 

  53. Ge LL, Tian HQ, Russell SD. Calcium function and distribution during fertilization in angiosperms. Am J Bot. 2007;94:1046–60.

    Article  CAS  PubMed  Google Scholar 

  54. Lang V, Pertl-Obermeyer H, Safiarian MJ, Obermeyer G. Pump up the volume - a central role for the plasma membrane H(+) pump in pollen germination and tube growth. Protoplasma. 2014;251:477–88.

    Article  CAS  PubMed  Google Scholar 

  55. Certal AC, Almeida RB, Carvalho LM, Wong E, Moreno N, Michard E, Carneiro J, Rodriguez-Leon J, Wu HM, Cheung AY, Feijo JA. Exclusion of a proton ATPase from the apical membrane is associated with cell polarity and tip growth in Nicotiana tabacum pollen tubes. Plant Cell. 2008;20:614–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Leigh RA, Jones RGW. A hypothesis relating critical potassium concentrations for growth to the distribution and functions of this Ion in the plant-cell. New Phytol. 1984;97:1–13.

    Article  CAS  Google Scholar 

  57. Zonia L, Cordeiro S, Tupy J, Feijo JA. Oscillatory chloride efflux at the pollen tube apex has a role in growth and cell volume regulation and is targeted by inositol 3,4,5,6-tetrakisphosphate. Plant Cell. 2002;14:2233–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Gutermuth T, Lassig R, Portes MT, Maierhofer T, Romeis T, Borst JW, Hedrich R, Feijo JA, Konrad KR. Pollen tube growth regulation by free anions depends on the interaction between the anion channel SLAH3 and calcium-dependent protein kinases CPK2 and CPK20. Plant Cell. 2013;25:4525–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Hamilton ES, Jensen GS, Maksaev G, Katims A, Sherp AM, Haswell ES. Mechanosensitive channel MSL8 regulates osmotic forces during pollen hydration and germination. Science. 2015;350:438–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Scholz-Starke J, Buttner M, Sauer N. AtSTP6, a new pollen-specific H + −monosaccharide symporter from Arabidopsis. Plant Physiol. 2003;131:70–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Rottmann T, Zierer W, Subert C, Sauer N, Stadler R. STP10 encodes a high-affinity monosaccharide transporter and is induced under low-glucose conditions in pollen tubes of Arabidopsis. J Exp Bot. 2016;67:2387–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Hirose T, Zhang Z, Miyao A, Hirochika H, Ohsugi R, Terao T. Disruption of a gene for rice sucrose transporter, OsSUT1, impairs pollen function but pollen maturation is unaffected. J Exp Bot. 2010;61:3639–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Johnson DA, Thomas MA. The monosaccharide transporter gene family in Arabidopsis and rice: a history of duplications, adaptive evolution, and functional divergence. Mol Biol Evol. 2007;24:2412–23.

    Article  CAS  PubMed  Google Scholar 

  64. Wang Y, Xiao Y, Zhang Y, Chai C, Wei G, Wei X, Xu H, Wang M, Ouwerkerk PB, Zhu Z. Molecular cloning, functional characterization and expression analysis of a novel monosaccharide transporter gene OsMST6 from rice (Oryza sativa L.). Planta. 2008;228:525–35.

    Article  CAS  PubMed  Google Scholar 

  65. Ariizumi T, Toriyama K. Genetic regulation of sporopollenin synthesis and pollen exine development. Annu Rev Plant Biol. 2011;62:437–60.

    Article  CAS  PubMed  Google Scholar 

  66. Verrier PJ, Bird D, Burla B, Dassa E, Forestier C, Geisler M, Klein M, Kolukisaoglu U, Lee Y, Martinoia E, et al. Plant ABC proteins--a unified nomenclature and updated inventory. Trends Plant Sci. 2008;13:151–9.

    Article  CAS  PubMed  Google Scholar 

  67. Wu L, Guan Y, Wu Z, Yang K, Lv J, Converse R, Huang Y, Mao J, Zhao Y, Wang Z, et al. OsABCG15 encodes a membrane protein that plays an important role in anther cuticle and pollen exine formation in rice. Plant Cell Rep. 2014;33:1881–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Tanaka N, Uraguchi S, Saito A, Kajikawa M, Kasai K, Sato Y, Nagamura Y, Fujiwara T. Roles of pollen-specific boron efflux transporter, OsBOR4, in the rice fertilization process. Plant Cell Physiol. 2013;54:2011–9.

    Article  CAS  PubMed  Google Scholar 

  69. Baxter I, Tchieu J, Sussman MR, Boutry M, Palmgren MG, Gribskov M, Harper JF, Axelsen KB. Genomic comparison of P-type ATPase ion pumps in Arabidopsis and rice. Plant Physiol. 2003;132:618–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Gomes E, Jakobsen MK, Axelsen KB, Geisler M, Palmgren MG. Chilling tolerance in Arabidopsis involves ALA1, a member of a new family of putative aminophospholipid translocases. Plant Cell. 2000;12:2441–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Schwacke R, Grallath S, Breitkreuz KE, Stransky E, Stransky H, Frommer WB, Rentsch D. LeProT1, a transporter for proline, glycine betaine, and gamma-amino butyric acid in tomato pollen. Plant Cell. 1999;11:377–92.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references


We thank Dr. Marc Boutry (University of Louvain, Louvain-la-Neuve, Belgium) for providing the antibody against plasma membrane H+-ATPase from N. plumbaginifolia. We also acknowledge Dr. Wenhai Jin and Lihai Guo (AB SCIEX Asia Pacific Application Support Centre) for technical assistance in MS analyses.


This work was supported by the Ministry of Science and Technology (grant nos. 2013CB945101 and 2012CB910504).

Availability of data and materials

The relevant supporting data can be found within the article and its additional files. The mass spectrometry data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD004994.

Authors’ contributions

Conceived and designed the experiments: TW. Performed the experiments: NY. Analyzed the data: NY TW. Contributed reagents/materials/analysis tools:TW. Wrote the paper: TW NY. Both authors read and approved the final manuscript.

Competing interests

The authors delare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Tai Wang.

Additional files

Additional file 1:

RapiGest-assisted PM protein digestion. Proteins of carbonate-washed PM proteins treated without (A, B) or with trypsin coupled with 0.2% RapiGest (C, D) were separated by 10% SDS-PAGE and stained with Coomassie brilliant blue. MPG, mature pollen grain; GPG, germinated pollen grain. 20 μg loaded per lane. (PDF 1401 kb)

Additional file 2:

Separation of iTRAQ-labeled peptides through SCX chromatography. The peptides from experiments 1 and 2 were fractionized by strong cation exchange chromatography (SCX) independently. After eluted by ammonium chloride, both of them were recombined into 16 fractions according to the 214 nm UV absorbance signal. (PDF 2055 kb)

Additional file 3:

Summary of spectra, peptides, and proteins identified in the 2 independent iTRAQ experiments. (PDF 11 kb)

Additional file 4:

All proteins identified in the experiments and their quantitative information. (XLSX 590 kb)

Additional file 5:

Plasma membrane-related proteins. (XLSX 240 kb)

Additional file 6:

Pearson analysis shows high quantitative reproducibility between the two independent iTRAQ experiments. ln [Exp 1], ln-transformed protein expression change ratio from experiment 1 (GPG/MPG = 117/115); ln [Exp 2], ln-transformed protein expression change ratio from experiment 2 (GPG/MPG = 115/117). (PDF 866 kb)

Additional file 7:

Mass spectrometry spectra show the repeatability of iTRAQ quantitative information. (PDF 634 kb)

Additional file 8:

Differentially expressed plasma membrane-related proteins. (XLSX 56 kb)

Additional file 9:

PMrPs and differentially expressed PMrPs show biased distribution in chromosomes 1, 2 and 3. A, PMrPs of each functional group distributed in rice 12 chromosomes. B, differentially expressed PMrPs of each functional group distributed in rice 12 chromosomes. (PDF 3061 kb)

Additional file 10:

Western blot validation of the iTRAQ quantitative information. (PDF 94 kb)

Additional file 11:

Summary of receptor-like kinases (RLKs) in Arabidopsis thaliana and Oryza sativa. (XLSX 26 kb)

Additional file 12:

Classification of 209 transporters identified from rice MPG and GPG plasma membrane. (XLSX 39 kb)

Additional file 13:

Qualitative comparison of rice pollen transcriptome and proteome. (XLSX 204 kb)

Additional file 14:

Quantitative comparison of pollen transcriptome and proteome. (XLSX 76 kb)

Additional file 15:

Multiple sequence alignment of CrRLK1Ls. Sequence alignments of identified rice pollen CrRLK1Ls gi|56783691, gi|115462979 and gi|25553554, and functionally known HERK1 (AT3G46290), HERK2 (AT1G30570), THE1 (AT5G54380), ANX1 (AT3G04690), ANX2 (AT5G28680) and FER (AT3G51550) from Arabidopsis. These CrRLK1Ls have a variable malectin domain (underlined in purple) and highly conserved kinase domain (in red) along with one transmembrane domain (in orange) predicted by SMART. Identical and similar amino acids are shaded in blue and yellow. (PDF 18045 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, N., Wang, T. Comparative proteomic analysis reveals a dynamic pollen plasma membrane protein map and the membrane landscape of receptor-like kinases and transporters important for pollen tube growth and interaction with pistils in rice. BMC Plant Biol 17, 2 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: