2.1 Plant materials
Adult trees of the pear varieties ‘Huanghua’ (H-H, S1S2), ‘Xizilv’ (XZL, S1S4), and ‘Nijisseiki’ (NI, S2S4) growing in the orchards of Nanjing Agricultural University, Jiangsu, China were used in this study. The leaves, pistils, and pollen were collected and stored at − 80°C for later analyses. S-RNase gene homozygotes were identified from the seedlings of the selfed progeny of ‘Huanghua’. The S1-homozygote H-7 (S1S1) and the S2-homozygote H-9 (S2S2) were used in this study.
2.2 Bud stage artificial self-pollination
Artificial self-pollination of ‘Huanghua’ at the bud stage (4, 6, and 8 days before full bloom) was carried out in 2008 and 2009. Before anthesis, two to three flowers per inflorescence were selected, emasculated, and enclosed in insect-proof bags to prevent contamination after artificial self-pollination. Fruit set was recorded 3 weeks later and fruits were collected when mature.
2.3 PCR amplification and sequence analysis of genomic S-RNases
Genomic DNA was extracted from young leaves of each cultivar using the CTAB (cetyl trimethyl ammonium bromide) method (Macklin, Shanghai, China). S-RNase alleles were identified by PCR analysis using the primer pair PF and PR [13], which targeted the C1 consensus conserved region and the non-consensus conserved region located downstream of the C2 region of S-RNases, respectively (Supplemental Table 1). Amplifications were conducted on a PTC-200 thermal cycler (Bio-Rad, Hercules, CA, USA) with the following cycle parameters: 3 min at 94 ℃ for pre-denaturation, followed by ten cycles of 30 s at 94 ℃, 30 s at 50 ℃, and 45 s at 70 ℃, then 25 cycles of 30 s at 94 ℃, 30 s at 48 ℃, and 45 s at 70 ℃, with a final extension for 7 min at 70 ℃. The products were separated by electrophoresis on a 1.5% (w/v) agarose gel. Band scoring was carried out using a standard DL2000 ladder (TaKaRa, Dalian, China). S1-RNase and S2-RNase alleles were discriminated on the basis of fragment length.
2.4 S-RNase gene quantitative real-time PCR assay
S-RNase gene transcript levels were analyzed by quantitative RT-PCR using the Rotor-Gene 6 fluorescence quantitative PCR kit (Beijing Heros Technology Company, Shanghai, China). The RT-PCR system was as follows: 5 µL 5×PCR buffer, 4 µL MgCl2 (25 mmol l− 1), 1 µL dNTP mixture (10 mmol l− 1), 0.5 µL specific primers (10 µmol l− 1), 1.25 µL 20× Eva green fluorescent dye, 1.5 U Taq enzyme (2.5 U·µL− 1), 1 µL template, and ddH2O to 25 µL. The sample was pre-degenerated at 95 ℃ for 5 min, 95 ℃ for 15 s, and 65 ℃ for 35 s (for fluorescence detection) over 40 cycles. The S1-RNase allele-specific primers were S1-F and S1-R. The S2-RNase allele-specific primers were S2-F and S2-R. A relative quantification method was used to calculate relative gene transcript levels [14]. The transcript level of each gene in each sample was compared with that of Actin (AF386514), which allowed for direct comparison of Actin gene-normalized RT-PCR ratios among samples. The primers Pactin-F and Pactin-R were used to amplify Actin (Supplemental Table 1).
2.5 Analysis of nuclear DNA content by flow cytometry
Nuclear protoplasm was isolated by chopping 100 mg young leaf tissue with a sharp scalpel in a glass Petri dish containing 500 µL LB01 lysis buffer with the following composition: 15 mM Tris, 2 mM EDTA-2Na, 80 mM KCl, 20 mM NaCl, 0.5 mM spermine, 15 mM β-mercaptoethanol, and 0.1% (v/v) Triton X-100, with a final pH of 7.4 [15]. To remove leaf debris, the suspension of released nuclei was filtered through a 21-µm nylon mesh. The nuclei were stained with propidium iodide (PI) at 100 mg·mL− 1 and incubated with 100 mg·mL− 1 RNase for 15–30 min prior to analysis. Unstained chicken red blood cell nuclei were used as an internal standard and were added to suspensions of isolated nuclei prior to staining. After isolation, nuclei were fixed in ice-cold 3:1 fixative (ethanol/glacial acetic acid) and stored in 70% (v/v) ethanol at − 20°C. The PI-stained nuclei suspensions were analyzed with an Epics XL flow cytometer (Beckman Coulter, Miami, FL, USA) and sorter equipped with an argon ion laser at λ = 514 nm adjusted to a 100-mW power output. Integral fluorescence and the height and width of fluorescence pulses emitted from the nuclei were collected through a 620 nm long-pass filter.
2.6 Preparation, concentration, and activity of S-RNase
Four grams of styles were prepared for isolating each S-RNase following our previously described method [6, 8]. The isolated S-RNase sample was stored in Eppendorf tubes at − 80 ℃. The S-RNase concentration and activity were determined by the methods of Bradford [16] and Brown and Ho [17], respectively.
2.7 Pollination test and segregation of S-alleles
To check cross-(in)compatibility in the three cultivars/lines, self-pollination and cross-pollination tests were conducted in ‘Huanghua’, ‘H-7’, and ‘H-9’. Field pollination experiments were carried out in 2014 and 2015. Before anthesis, 2–3 flowers per inflorescence were selected, emasculated, and enclosed in insect-proof bags to prevent contamination after pollination. Subsequent fruit set was recorded 3 weeks later, and fruit were collected when mature. Pistils were harvested and fixed in FAA solution (5% v/v formalin, 5% v/v acetic acid, 45% v/v ethanol) 48 h after self-pollination. Pistils were washed in distilled water three times to remove FAA, and then softened by autoclaving for 30 min at 100 kPa in 5% (w/v) sodium sulfite solution. The pistils were stained for 24 h with 0.1% (v/v) aniline blue (0.1 N potassium phosphate, pH 8.7). After the epidermis was removed, the pistils were examined under a microscope (BX60, Olympus, Tokyo, Japan) equipped with a BH2-RFL-T2 ultraviolet light source and an Osram HBO 100 W/2 high-pressure mercury lamp (Olympus, Tokyo, Japan) [18]. To study the segregation of the S-RNase gene, seedling genotypes were detected by PCR using the primer set PF/PR for S-RNase.
2.8 Pollen culture and SI challenge
H-7 pollen grains were pre-cultured for 2 h at 25 ℃ in a basal medium in the dark according to the method of [6]. The basal medium consisted of a 2-(N-morpholine)–ethanesulfonic acid (MES)-NaOH buffer supplemented with 10% (w/v) sucrose, 15% (v/v) polyethylene glycol 4000, 0.01% (w/v) H3BO3, 0.07% (w/v) Ca(NO3)2 4H2O, 0.02% (w/v) MgSO4·7H2O and 0.01% (w/v) KNO3, pH 6.0–6.5. After pre-culture, H-7 stylar S-RNase was added to the medium as an SI challenge, and H-9 stylar S-RNase was added to the medium as a compatible treatment. Medium without S-RNase was used as the control. The final activity of the S-RNases in the basal medium was 0.15 U.
2.9 Depolymerization assays and quantification of filamentous actin in vivo
An F-actin depolymerization assay was conducted and the amount of F-actin present in the samples was quantified as previously described [19]. After each reaction, the pollen samples were stabilized in 200 mM MES-buffered saline and fixed in 4% (w/v) paraformaldehyde. The samples were then stained individually with 5 mM rhodamine-phalloidin for the depolymerization assay or simultaneously with 5 mM rhodamine-phalloidin and 5 mM ethidium bromide (EB) for the quantification assay. A confocal microscope (LSM700, Zeiss, Jena, Germany) was used to examine the levels of actin depolymerization. Phalloidin fluorescence divided by EB fluorescence was used as an index of actin filament levels. The bound phalloidin and EB were eluted with methanol and quantified by spectrofluorometry with excitation and emission wavelengths of 492 and 514 nm, respectively, for phalloidin, and 513 and 615 nm, respectively, for EB.
2.10 Evaluation of programmed cell death in pollen tubes
We evaluated PCD in pollen tubes by terminal-deoxynucleotidyl transferase-mediated nick-end labeling (TUNEL) staining as previously described [20]. After fixing the samples in 4% (w/v) paraformaldehyde for 2 h, the pollen tubes were transferred to 70% (v/v) ethanol and incubated overnight at − 20 ℃. The Dead End Fluorometric TUNEL system (Promega, Beijing, China) was used to examine PCD in the samples. After the samples were washed with citrate buffer (pH 4.1), the pollen tubes were stained with 50 µg·mL− 1 4′,6-diamidino-2-phenylindole (DAPI) for 10 min at room temperature and then analyzed using a Fluoview FV1200 confocal microscope (Olympus). The TUNEL signal was detected at excitation and emission wavelengths of 540 and 620 nm, respectively, and the DAPI signal was detected at excitation and emission wavelengths of 360 and 420 nm, respectively. Positive TUNEL staining in the pollen tubes was indicative of PCD.
2.11 Protein extraction
One gram of frozen style from H-H, S1S1, and S2S2 was ground to a powder in liquid nitrogen. Each sample included three biological replicates. The dry powder was dissolved in 0.6 mL lysis buffer (1 M sucrose, 0.5 M Tris 8.0, 0.1 M KCl, 50 mM ascorbic acid, 1% NP 40, 1% NaDOC, 10 mM EDTA, 10 mM DTT, 1% protease inhibitor mixture and 1% phosphatase inhibitor mixture) and placed on ice for 10 min [21]. Then, 0.6 mL Tris-phenol was added to the samples, which were kept on ice for 10 min after intensive mixing. The samples were then centrifuged at 5500×g for 10 min at 4 ℃. The supernatant was transferred into a new tube, and four volumes of cold ammonium acetate/methanol were added to the tube, followed by incubation at − 20 ℃ overnight. After incubation, the pellet was collected by centrifugation at 1600×g for 10 min at 4 ℃, and then resuspended in cold methanol and kept at − 20 ℃ for 1 h. After centrifugation at 1600×g for 10 min at 4 ℃, the pellet was collected, resuspended in cold acetone, and washed twice. After washing, the mixture was centrifuged at 1600×g for 10 min at 4 ℃, and then the pellet collected and air-dried. Then, 0.5 mL lysis buffer (8 M Urea, 50 mM Tris 8.0, 1% NP40, 1% NaDOC, 10 mM EDTA, 5 mM DTT, 1% protease inhibitor mixture, and 1% phosphatase inhibitor mixture) was added. The samples were sonicated at 200 W for 5 min, and then centrifuged at 2000×g for 10 min at 4 ℃. The supernatant was collected, and the protein concentration was quantified using the Bio-Rad Protein Assay Kit (Bio-Rad, Hercules, CA, USA) following the manufacturer’s instructions.
2.12 iTRAQ labeling and strong cation exchange (SCX) fractionation
From every sample, 100 µg protein was digested with Trypsin Gold (Promega, Madison, WI, USA) at a protein: trypsin ratio of 30:1 for 16 h at 37°C. After trypsin digestion, the peptides were dried by vacuum centrifugation, and then dissolved in 0.5 M TEAB using 8-plex iTRAQ reagent [22]. Then, one unit of iTRAQ reagent was thawed and diluted with 24 µL isopropanol. The samples were marked with iTRAQ labels as follows: iT-S1S2, iT-S1S1, iT-S2S2. The peptides were labeled with isobaric tags at 25 ℃ with 2 h. Finally, the marked peptide mixtures were combined and dried under vacuum centrifugation.
SCX Chromatography was conducted using an LC-20AB HPLC pump system (Shimadzu, Japan). The iTRAQ labeled peptide mixture was reconstituted in 2 mL buffer A (25 mM NaH2PO4 in 25% ACN, pH 2.7) and then loaded onto an Ultremex SCX column (particle size, 5 µm). The peptides were eluted in a gradient formed by buffer A and buffer B (1 M KCl in 25% ACN, 25 mM NaH2PO4, pH 2.7), according to the following program: buffer A for 10 min at 1 mL/min, then 5–60% buffer B over 27 min, and 60–100% buffer B over 1 min. The column was washed with 100% buffer B for 1 min and then equilibrated with buffer A for 10 min before the next injection. Fractions were collected every 1 min. The elution of products was detected by measuring absorbance at 214 nm. Finally, the eluted peptides collected in 20 fractions were desalted on a Strata X C18 column and then vacuum-dried [23].
2.13 LC–ESI–MS/MS analysis based on Triple TOF 5600
Each dried peptide fraction was reconstituted in buffer A (0.1% FA, 5% ACN) and then centrifuged at 20,000×g for 10 min. The average concentration of final peptides was about 0.5 µg/µL. For the LC analysis, 10 µL of supernatant was injected into a LC-20AD nanoHPLC equipped with a 2-cm C18 trap column. The peptides were separated on an internally packed analytical column. All samples were loaded at a rate of 8 µL/min over 4 min, and separated by a gradient of buffer A and B (0.1% FA, 95% ACN) supplied at a flow rate of 300 µL/min. The gradient program was as follows: 2–35% B over 35 min, increased with a linear gradient to 60% over 5 min, then to 80% over 2 min, and held at 80% B for 4 min before returning to 5% buffer B within 1 min.
Data were acquired from the Triple TOF 5600 System equipped with a pulled quartz tip as transmitter and a Nanospray III source, with the following settings: 30 psi curtain gas, 150 ℃ interface heater temperature, 15 psi atomization gas, and 2.5 kV ion spray voltage. A RP (≥ 30,000 FWHM) was used for operating the MS for TOF-MS scanning [24]. For IDA, if the threshold was exceeded at 120 counts per s above 2 to 5+, survey scans were performed every 250 ms and 30 production scans were assembled. The total cycle time was fixed at 3.3 s and the Q2 transfer window at 100 Da was set to 100%. By monitoring the 40 GHz multi-channel TDC detector with a four-anode channel detector, four time bins were summarized at 11 kHz pulse frequency per scan. The 35 ± 5 eV sweep collision energy setting was used in conjunction with iTRAQ to adjust the rolling collision energy to perform collision-induced dissociation of all precursor ions. Finally, dynamic exclusion was set to 50% of the peak width (15 s), and precursors were refreshed from the exclusion list.
2.14 iTRAQ protein identification and quantification
The raw data files obtained from the Triple TOF 5600 system were converted into MGF files using 5600 MS converter and then the MGF files were searched. The triple TOF 5600 mass spectrometer used in this study had a high mass accuracy (less than 2 ppm), allowing for reliable peptide identification. The charge states of peptides were set to + 2 and + 3. For proteome analyses, the MS data were processed using Proteome Discoverer software with default parameters to produce a list of peaks. Data were collected using Analyst QS 2.0 software, and proteins were identified and quantified by Mascot 2.3.02. For iTRAQ quantification, peptides for quantification were automatically selected by an algorithm to calculate reporter peak area, error factor (EF), and p-value (using the default parameters in the Mascot software package). The resulting data sets were automatically corrected for deviations to eliminate any changes due to uneven mixing of samples with different markers. The sequences of Chinese white pear proteins (https://www.peargenome.njau.edu.cn/) were identified using the Mascot search engine (Matrix Science). The ratios between samples were obtained directly based on the protein abundance for any given protein. Proteins showing 1.5-fold or higher difference in abundance between the treated and control samples and a p-value of less than 0.05 were identified as differentially expressed proteins (DEPs). Quantification was performed at the peptide level according to the procedure described at http://www.matrixscience.com/help/quant_statistics_help.html. Each treatment was repeated once for the iTRAQ analysis. The mass spectrometry proteomics data were stored in the PRIDE partner repository with the data set identifier PXD0043543 in the ProteomeXchange [25].