The Scientific World Journal

Volume 2014 (2014), Article ID 147061, 9 pages

http://dx.doi.org/10.1155/2014/147061

Research Article

Sequence Characterization and Spatiotemporal Expression Patterns of PbS ₂₆- RNase Gene in Chinese White Pear (Pyrus bretschneideri)

Lin Zhang,1,2 Baoguang Jia,¹ Feng Zou,¹ Xiaofeng Tan,¹ Min Liu,¹ Zhibo Song,¹ Yanling Zeng,¹ Nan Jiang,¹ and Deyi Yuan¹

¹Key Laboratory of Cultivation and Protection for Non-Wood Forest Trees, Ministry of Education, Central South University of Forestry and Technology, Changsha 410004, China

²Department of Plant Science and Landscape Architecture, University of Connecticut, Storrs, CT 06269, USA

Received 6 September 2013; Accepted 19 January 2014; Published 5 March 2014

Academic Editors: A. Ferrante and Q.-H. Zhu

Copyright © 2014 Lin Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Many flowering plants exhibit an important intraspecific reproductive barrier phenomenon, that is, self-incompatibility (SI), in which S-RNase genes play a significant role. To clarify the specific function of S-RNase genes in Chinese pears, the full length cDNA of PbS ₂₆- RNase was isolated by rapid amplification of cDNA ends (RACE) technology from Chinese white pear (Pyrus bretschneideri) cultivar “Hongpisu.” The cDNA sequence for PbS ₂₆- RNase was deposited in GenBank under accession number EU081888. At the amino acid level, the PbS ₂₆- RNase displayed the highest similarity (96.9%) with PcSa-RNase of P. communis, and only seven amino acid differences were present in the two S-RNases. Phylogenetic analysis of rosaceous S-RNases indicated that the PbS ₂₆- RNase clustered with maloideous S-RNases, forming a subfamily-specific not a species-specific group. The PbS ₂₆- RNase gene was specifically expressed in the style but not other tissues/organs. The expression level of the PbS ₂₆- RNase gene rapidly increased at bell balloon stage (BBS), and then it dropped after pollination. However, the abundance of the PbS ₂₆- RNase gene transcript in the style was greater after cross-pollination than after self-pollination. In addition, a method for rapidly detecting the PbS ₂₆- RNase gene was developed via allele-specific primers design. The present study could provide a scientific basis for fully clarifying the mechanism of pear SI at the molecular level.

1. Introduction

In flowering plants (also called angiosperms), an intraspecific reproductive mechanism, self-incompatibility (SI), has so far been widely known. This mechanism can prevent inbreeding and promote outcrossing [1]. One type of the SI, gametophytic SI (GSI), mainly exists in three plant families including Solanaceae, Scrophulariaceae, and Rosaceae, and it is controlled by a multiallelic S-locus [1–3]. The S-locus harbors a series of S-alleles that are divided into two types, that is, female S-genes and male S-genes. The expression products of the S-allele in the female organ (style) are S-RNases, a class of polymorphic proteins with ribonuclease activity.

Many fruit trees of Rosaceae, including apple, pear, and almond, belong to GSI-type plants. Among rosaceous members, Japanese pear (Pyrus pyrifolia) is the first species in which pistil S-allele specific proteins were successfully identified as S-RNases [3]. Afterwards, cDNAs encoding seven S-RNases (S₁ to S₇) of Japanese pear were cloned with their primary structure being determined [4, 5]. A hypervariable (HV) region, possibly associated with S-allele specific recognition, was identified [5]. Based on sequence characterization of Japanese pear S-RNase genes, a PCR-restriction fragment length polymorphism (PCR-RFLP) system was established for identifying S-genotypes of cultivars and searching new S-alleles [5]. Using this system, four additional alleles S₈, S₉, S₁₀, and Sk were discovered and cloned, and accordingly the system was modified [6–8].

Pears native to China include four main cultivated species, that is, Chinese sand pear (P. pyrifolia), Chinese white pear (P. bretschneideri), Ussurian pear (P. ussuriensis), and Xinjiang pear (P. sinkiangensis). Chinese pears were considered to be the major origin of Asian pears [9]. Knowledge of S-genotype is essential for pear production and breeding programs. The conventional pollination-based method for cultivars S-genotyping is time consuming and labor intensive. Moreover, the obtained results are often ambiguous due to environmental and physiological factors. PCR-based molecular techniques, such as PCR-RFLP and direct sequencing of amplified fragments, provide a rapid, reliable, and easily manipulated approach for pear cultivars S-genotyping. To date, a series of S-alleles have been found in Chinese pears and approximately 100 cultivars have been S-genotyped [10–15]. Of these S-alleles, only three full length cDNAs encoding S₁₃-, S₁₆-, and S₃₄-RNases were reported and analyzed for their possible function in SI reaction [16, 17]. More recently, several pollen specific polymorphic F-box genes were isolated and considered to be the candidates of male S-genes in pears [18–24]. To fully clarify the mechanism of pear SI, it is necessary to investigate the interaction of the S-RNase genes with the male S-genes.

In this study we cloned and characterized the full length of cDNA of the PbS ₂₆-RNase gene in P. bretschneideri. We studied the spatiotemporal expression patterns of this gene and compared its transcript abundance between incompatible pollination and compatible pollination. In addition, we developed a method for rapidly and accurately detecting the PbS ₂₆-RNase gene.

2. Materials and Methods

2.1. Plant Materials

Plant materials used in this study include P. bretschneideri cultivars “Hongpisu” (S₁₂S₂₆) [13], “Hongxiangsu,” “Jinhua” (S₁₃S₁₈) [25], and 46 seedlings derived from cross of “Hongpisu” × “Jinhua” which are planted at Pear Germplasm Repository of Central South University of Forestry and Technology, Zhuzhou city, Hunan province. For cDNA cloning of PbS ₂₆-RNase gene, styles at bell stage (BS) were collected from “Hongpisu.” For expression analysis of the PbS ₂₆-RNase gene, styles of “Hongpisu” were collected at 96 h before anthesis (BA), 48 h BA, bell balloon stage (BBS), 12 h after pollination (AP), and 24 h AP, respectively. Moreover, stamens, roots, tender stems, and young leaves were also sampled for expression detection of the PbS ₂₆-RNase gene. The plant materials were sampled and immediately frozen in liquid nitrogen and stored at −80°C until used.

2.2. Isolation of Nucleic Acids

Genomic DNA was isolated from young leaves with a CTAB method as described in our previous study [22]. Total RNA was extracted using Micro-to-Midi Total RNA Purification System (Invitrogen) and following the manufacturer’s protocols. The quality of DNA and RNA was checked by electrophoresis in 1.0% TAE agarose gel.

2.3. Isolation of Full Length cDNA of PbS 26-RNase Gene

The full length cDNA of the PbS ₂₆-RNase was amplified using 3′ RACE System for Rapid Amplification of cDNA Ends (Invitrogen, Version E). Total RNA was reverse-transcribed using SuperScript II Reverse Transcriptase with Adapter Primer provided in the system under condition of 10 min at 70°C, 50 min at 42°C, 15 min at 70°C, and 20 min at 37°C. The second PCR amplification was conducted with primers PF2 (5′- TGCCTCGCTCTTGAACAAA-3′) and AUAP. The PF2 anneals to 5′ UTR of most of pear S-RNases and was used for pear S-RNase cloning [6, 16]. The AUAP primer was provided in the 3′ RACE system. The amplification program consisted of 35 cycles of 30 s at 94°C, 30 s at 59°C, 3 min at 72°C, and a final extension of 7 min at 72°C.

2.4. Controlled Pollinations

In order to elucidate the expression profiles after pollination of the PbS ₂₆-RNase gene, self-pollination (“Hongpisu” × “Hongpisu”) and cross-pollination (“Hongpisu” × “Jinhua”) tests were performed using 30 flowers with three replicates for each pollination combination. Anthers were collected from flowers at the BBS and incubated for 30 h to collect the pollens. For each pollination combination, flowers were emasculated, pollinated, and covered with paper bags to avoid contamination. In order to monitor pollen tubes growth in the style, samples were taken at 2, 8, 12, and 24 h after pollination and fixed in FAA solution. After fixation, the flowers were maintained in 70% ethanol and stored at 4°C until the time of pollen tubes growth was observed using a fluorescence microscope. According to the pollen tube growth, styles at 12 h and 24 h after pollination were collected for quantitative RT-PCR (qPCR) analysis. In addition, the percentages of fruit set and the number of full seeds were recorded four months after pollination. The fruit set percentage of below 30% was considered as incompatible [7].

2.5. Semiquantitative RT-PCR

The expression patterns of the PbS ₂₆-RNase gene in different tissues/organs were studied by semiquantitative RT-PCR (semi-qPCR), in which pear Actin gene (GenBank accession number GU830958) was used as the reference gene. The PbS ₂₆-RNase gene was amplified with primer pair S26PF1 (5′-GAATACAACCTTGAATTCTACTAA-3′) and S26PR1 (5′-CACAGCTGCCATGTTTGTT-3′), which could produce a fragment of 131 bp. The pear Actin gene was amplified with primer pair ActinF (5′-TGGTGTCATGGTTGGTATGG-3′) and ActinR (5′- CAGGAGCAACACGAAGTTCA-3′), which could produce a fragment of 173 bp. The reaction volumes were 20 μL containing 1×PCR buffer, 0.5 mM dNTPs, 1.5 mM MgCl₂, 0.5 mM of each primer, 1U Taq polymerase, and the 2.0 μL synthesized cDNA first strand as template. The amplification program consisted of 35 cycles of 30 s at 94°C, 30 s at 58°C, 3 min at 72°C, and a final extension of 7 min at 72°C.

2.6. Quantitative RT-PCR

Quantitative RT-PCR (qPCR) was carried out with three replicates in each reaction using the BIO-RAD CFX system (Bio-Rad). PbS ₂₆-RNase gene-specific primers, reference gene, and its corresponding primers were the same as described in semi-qPCR. PCR was performed in a 20 μL volume containing 2 μL diluted cDNA, 250 nM each primer, and 1 × SYBR Premix Ex Taq II (TaKaRa) using the same conditions with semiquantitative RT-PCR. The results from PbS ₂₆-RNase gene-specific amplification were analyzed using the comparative Cq method which uses an arithmetic formula, 2-ΔΔCq, to obtain results for relative quantification [26].

2.7. Development of the Rapid Method for the PbS ₂₆ -RNase Gene Detection

To specifically detect the S₂₆-allele, a primer pair S26PF2 (5′-GCACAGGAAATGACCCATC-3′) and S26PR2 (5′-GGTTCGATCGAGTACGTTG-3′) was designed via multiple genomic DNA sequence alignment that would produce a fragment of 254 bp by genomic PCR. The S26PF2 was located in intron region, and the S26PR2 was located between conserved 4 region (C4) and conserved 5 region (C5). The amplification program consisted of 35 cycles of 30 s at 94°C, 30 s at 53°C, 3 min at 72°C, and a final extension of 7 min at 72°C. The amplification products were detected by 1.5% agarose gel electrophoresis.

2.8. Cloning, Sequencing, and Sequence Analysis of PCR Products

The purified PCR products were TA-ligated into pMD18-T Vector (Takara, Japan) and transformed into Escherichia coli (E. coli) strain DH5α. The positive clones are screened by the white/blue colony and identified by plasmid PCR with BcaBESTTM sequencing primer M13-47 (5′-CGCCAGGGTTTTCCCAGTCACGAC-3′) and BcaBESTTM sequencing primer RV-M (5′-GAGCGGATAACAATTTCACACAGG-3′). PCR reactions were conducted with the following cycling conditions: one initial step of 4 min at 94°C, 35 cycles of 30 s at 94°C, 30 s at 62°C, 1.5 min at 72°C, and one final extension step of 10 min at 72°C. Target clones were sequenced in both directions by Shanghai Biosune Biotechnology Ltd.

Sequence data were analyzed with DNAMAN (version 6.0.40). Prediction of isoelectric point (pI), molecular weight (Mw), and Prosite was performed by using online ExPASy proteomics tools (http://expasy.org/tools/pi_tool.html). The BLAST algorithm was used to search the NCBI GenBank (http://www.ncbi.nlm.nih.gov/) databases for homologous sequences. Motif analysis of the target gene was conducted by Motif Scan (http://myhits.isb-sib.ch/cgi-bin/motif_scan) and sequence comparisons. Multiple amino acid sequences alignment was performed using the DNAMAN software. Based on this alignment, a phylogenetic tree was generated. The following S-RNase sequences, together with the PbS ₂₆-RNase sequence, were used for the phylogenetic analysis: P. bretschneideri: PbS₁₃-RNase (DQ414812), PbS₁₆-RNase (DQ991388), and PbS₃₄-RNase (DQ414813); P. communis: PcSa-RNase (AB236430), PcSb-RNase (AB236429), PcSd-RNase (AB236427), PcSh-RNase (AB236431), PcSk-RNase (AB236432), PcSq-RNase (AB236424), and PcSr-RNase (AB236426); P. pyrifolia: PpS₁-RNase (AB002139), PpS₂-RNase (D49527), PpS₃-RNase (AB002140), PpS₄-RNase (D49528), PpS₅-RNase (D88282), PpS₆-RNase (AB002142), PpS₇-RNase (AB002143), PpS₈-RNase (AB104908), and PpS₉-RNase (AB104909); M. domestica: MdS₄-RNase (AF327223), MdS₁₀-RNase (AF327221), MdS₂₄-RNase (AF016920), MdS₂₆-RNase (AF016918), and MdS₃₁-RNase (DQ135990); P. armeniaca: ParSc-RNase (DQ422947); P. avium: PavS₁₂-RNase (AY259115), PavS₁₃-RNase (DQ385842), PavS₂₃-RNase (AY259114), PavS₂₄-RNase (AY259112), and PavS₂₅-RNase (AY259113); P. dulcis: PdSa-RNase (AB026836), PdSb-RNase (AB011469), PdSc-RNase (AB011470), and PdSd-RNase (AB011471); P. salicina: PsSa-RNase (AB026981), PsSb-RNase (AB026982), PsSc-RNase (AB084102), and PsSd-RNase (AB084103).

3. Results and Discussion

3.1. Sequence and Characterization of Full Length cDNA of the PbS 26 -RNase

Previously, the PbS ₂₆-RNase allele was identified in P. bretschneideri cultivar “Hongpisu” (S₁₂S₂₆) [13]. Its partial genomic sequence was isolated, but the functional completed cDNA has not been gained. The most common way to isolate full length cDNA includes the processes of 5′ RACE and 3′ RACE. The process of 5′ RACE is a more complicated and challenging procedure than 3′ RACE, with many more steps. The PF2 primer anneals to a conserved region in the 5′ untranslated region (UTR) of Japanese pear S-RNases [6]; accordingly the produced target gene contains initiation codon ATG in theory. Following RACE-PCR with primer PF2 and AUAP in this study, an approximately 900 bp product was yielded. Since the PF2 does not match the S₁₂-RNase [27], the 900 bp product was considered to correspond to the PbS ₂₆-RNase gene. The target clone was sequenced and confirmed as a full length cDNA encoding S₂₆-RNase by detailed sequence analysis.

The cDNA sequence for PbS ₂₆-RNase was 906 bp long with a complete open reading frame (ORF) of 684 bp encoding 228 amino acids (GenBank accession number EU081888). The putative initiation codon ATG was located at poisons 29–31. The termination codon TAA was present at positions 713–715, followed by 3′ UTR. The deduced amino acid sequences for the PbS ₂₆-RNase had similar structural characteristics of pear S-RNases, for example, a signal peptide, five conserved regions, and a HV region. The molecular weight (Mw) and the isoelectric point (pI) of the PbS ₂₆-RNase were predicted to be 26.2 kDa and 9.14, respectively, which agreed with the basic properties of pear S-RNases [7]. The PbS ₂₆-RNase showed the conserved motifs with two histidine residues, His-60 and His-116, which are essential for T2/S type RNase activity [28]. It also harbored eight cysteine residues (Cys-42, Cys-49, Cys-75, Cys-119, Cys-183, Cys-198,Cys-209, and Cys-221) that were mostly conserved in S-RNases. These eight cysteine residues could form four disulfide bridges, which plays a significant role in the formation or stabilization of their tertiary structure [7]. The PbS ₂₆-RNase also presented six potential N-glycosylation sites, of which only one, that is, Asn-, has glycans which were considered to be important for the folding and the stabilization of the core structure (Figure 1) [7].

Figure 1: Nucleotide and deduced amino acid sequences of the full length PbS ₂₆-RNase cDNA. The signal peptide, C1 to C5 domains, and HV region are shaded. Conserved cysteine residues, histidine residues essential for the RNase activity, and potential N-glycosylation sites are marked with symbols ▲, ●, and ◆, respectively, under the amino acid sequence. Asterisk indicates termination codon.

3.2. Amino Acid Sequences Comparison and Phylogenetic Analysis

The newly cloned PbS ₂₆-RNase gene allows us to better understand the evolutionary relationship of rosaceous S-RNases. The multiple amino acid sequences alignment was made using clustal X program with DNAMAN. A total of 39 complete amino acid sequences of S-RNases of eight species, including sweet cherry (P. avium), apricot (P. armeniaca), almond (P. duslics), plum (P. salicina), Japanese pear (P. pyrifolia), Chinese white pear (P. bretschneideri), European pear (P. communis), and apple (M. domestica), were used in the alignment. Among the P. bretschneideri species, the PbS ₂₆-RNase showed the lowest similarity (64.0%) with the PbS₁₃-RNase and the highest (90.4%) with the PbS₁₆-RNase. However, the PbS ₂₆-RNase showed extremely high similarities to P. communis PcSa-RNase (96.9%), which were remarkably higher than to P. bretschneideri S-RNases. Amino acid comparison between the PbS ₂₆-RNase and the PcSa-RNase showed that only seven amino acids differences were presented in the two S-RNases (Figure 2). The existence of such two S-alleles, in one case, agreed with the hypothesis that new S-alleles are generated by an accumulation of point mutations of earlier-formed alleles, which may or may not lead to differences in amino acid sequence [29]. In another case, the two genes also may generate from the same ancestor after separate mutations.

Figure 2: Alignment of amino acid sequences of the PbS ₂₆-RNase and the PcSa-RNase. Amino acid differences are indicated by asterisk.

Based on the multiple sequence alignment, a phylogenetic tree was generated (Figure 3). As indicated by the tree, the 39 S-RNases were divided into two subfamily-specific groups, that is, amygdaloideous S-RNases and maloideousS-RNases. The amygdaloideous group includes S-RNases of pears (P. pyrifolia, P. bretschneideri, and P. communis) and apple (M. domestica). The maloideous group includes those of almond (P. duslics), apricot (P. armeniaca), sweet cherry (P. avium), and plum (P. salicina). In the maloideous group, pear and apple S-RNases were closely related that they did not form species-specific subgroup. In the case of amygdaloideous group, the almond, apricot, plum, and sweet cherry S-RNases were grouped with one another and did not form species-specific subgroup either. The phylogenetic analysis of the study agreed with the proposal of Ushijima et al. [30]; namely, the divergence of S-RNase alleles predated speciation of the Rosaceae family, but it occurred shortly after the divergence into subfamilies.

Figure 3: Phylogenetic analysis of the members of S-RNase genes in rosaceous plants. The PbS ₂₆-RNase is boxed.

3.3. The PbS 26-RNase Gene Was Specifically Expressed in the Style

To investigate tissues/organs expression of PbS ₂₆-RNase gene, semi-qPCR was performed on total RNA of styles prior to pollination, stamens, roots, stems, and leaves of “Hongpisu” pear. Following RT-PCR with PbS ₂₆-RNase gene-specific primers, one fragment of expected sized (130 bp) was produced from styles, but not from the other tissues/organs (Figure 4). The results suggested that the PbS ₂₆-RNase gene was specific in pears. The expression patterns coincided with those of other S-RNase genes in Rosaceae plants such as apple [31], Japanese pear [7], sweet cherry [32], and Japanese apricot [33].

Figure 4: Semi-qPCR for the PbS ₂₆-RNase gene expression detection in different tissues/organs. 1–5: roots, stems, leaves, stamens, and styles.

3.4. Comparison of Expression Patterns of the PbS 26-RNase Gene between Self-Pollination and Cross-Pollination

Although many data show that the S-RNase genes are specific in the style, little is known about the expression profiles of the pear S-alleles during the style development, including before and after pollination. In this study, qPCR indicated that the PbS ₂₆-RNase gene showed an increased expression pattern in the developing style before pollination in “Hongpisu.” A slight transcript level of the Pb₂₆-RNase allele was detected at 96 h BA, and then it rapidly increased at 48 h BA and at BBS (Figure 7).

To further investigate the expression patterns after compatible pollination and incompatible pollination and compare their difference, controlled pollinations, that is, “Hongpisu” × “Hongpisu” and “Hongpisu” × “Jinhua,” were conducted. In self-pollination, pollen tubes were inhibited and twisted 12 h after pollination (Figure 5, a4). The self-pollination led to the low percentages of fruit set (5.8%) and seed set (on the average 3.7 seeds per fruit) (Figure 6), showing incompatible behaviors. In cross-pollination, in contrast, pollen tubes could grow to the cut end of the style, that is, the stylar base 24 h after pollination (Figure 5, b4). The cross-pollination resulted in remarkably higher percentages of fruit set (84.6%) and seed set (on the average 8.7 seeds per fruit) (Figure 6), showing compatible behaviors.

Figure 5: Fluorescence micrographs of squashed styles of “Hongpisu” harvested after different pollination treatments. (a1) Germination of pollen grains on stigma 2 h after self-pollination. (a2) Growth of pollen tubes in style transmitting tissue 6 h after self-pollination. (a3) Growth inhibition of pollen tubes in style transmitting tissue 8 h after self-pollination. (a4) Twisted pollen tubes 12 h after self-pollination. (b1) Germination of pollen grains on stigma 1 h after cross-pollination. (b2) Growth of pollen tubes in style transmitting tissue 8 h after cross-pollination. (b3) Growth of pollen tubes in style transmitting tissue 12 h after cross-pollination. (b4) Pollen tubes arrive at the base of the style 24 h after cross-pollination.

Figure 6: Percentages of fruit set and seed set of “Hongpisu” after controlled pollination.

(a)

(b)

Figure 7: QPCR analysis for the PbS ₂₆-RNase gene expression detection during styles development in “Hongpisu.”  (a) Relative expression of the PbS ₂₆-RNase gene in the style before and after self-pollination; (b) relative expression of the PbS ₂₆-RNase gene in the style before and after cross-pollination. BA: before anthesis; BBS: bell balloon stage; AP: after pollination.

After pollination, the expression level of PbS ₂₆-RNase allele gradually decreased in the style of both different pollinations (Figure 7). However, the expression level of the PbS ₂₆-RNase gene was higher after compatible pollination than after incompatible pollination, with the abundance of transcript being around 9-fold and 7-fold higher at 12 h AP and 24 AP, respectively, in cross-pollination than in self-pollination. This is the first report of the expression patterns of pear S-RNase genes in the developing style after pollination. The difference of S-RNase gene expression level between self-pollination and cross-pollination may be associated with the behaviors of different pollen S-genes. The pollen S-gene is identified as an F-box protein (SFB/SLF) which appears to act as a component of SCF complex. The SCF complex is associated with 26S proteasome degradation pathway in which non-self-S-RNases are targeted for degradation, and then the self-S-RNases are retained to specifically inhibit the growth of self-pollen tubes [34]. More recently, Xu et al. reported that two pollen specific SLF-interacting Skp1-like (SSK) proteins, PbSSK1 and PbSSK2, in P. bretschneideri could connect PbSLFs to PbCUL1 to form a putative canonical (SSK/CUL1/SLF) complex [24]. Based on the studies of Qiao et al. [34] and Xu et al. [24], we presume that pollen specific SFB₂₆/SLF₂₆ may be involved in the formation of 26S proteasome pathway and specifically degrade the non-self-S-RNases after pollination, resulting in the fact that the PbS ₂₆-RNase is reserved. Therefore, the abundance of the PbS ₂₆-RNase gene transcript in the style could be greater after cross-pollination than after self-pollination.

3.5. Development of the Rapid Method for Identifying the PbS 26-RNase Gene

In Japanese pear, a PCR-RFLP system was established for S-genotype assignments based on cloning of S ₁- to S ₉-RNase gene [5, 7]. Herein, we reported a different method for rapidly identifying S-RNase genes in Chinese pears. For the PbS ₂₆-RNase gene, a novel primer pair S26PF2 and S26PR2 was designed based on DNA sequence comparisons. Genomic DNA of “Hongxiangsu” was amplified with S26PF2 and S26PR2, and one specific fragment with expected size (254 bp) was generated. Sequencing of this fragment revealed that it corresponded to the PbS ₂₆-RNase allele. Based on the PbS ₂₆-RNase allele, the other one allele was easily identified as S₂₁ by further PCR amplification and DNA sequencing. The S26PF2 and S26PR2 primer combination was also applied in 46 progenies of “Hongpisu” (S₁₂S₂₆) × “Jinhua” (S₁₃S₁₈), and the 254 bp sized fragment was produced in 13 progenies (Figure 8), suggesting that the 13 individuals harbor the PbS ₂₆-RNase allele. These results showed that the newly designed PbS₂₆-specific primers could be used for rapid identification of the PbS ₂₆-RNase allele and will be useful in pear breeding programs.

Figure 8: Amplification of S-RNase alleles in 46 progenies of “Hongpisu” × “Jinhua” with S26PF2 and S26PR2 primers. M: 100 bp DNA ladder; lanes 1 to 46 correspond to the progeny numbers.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This study was supported by the Scientific and Technological Achievements Promotion Project of China State Forestry Administration, the National Natural Science Foundation of China (31272124), and the Senior Talent Introduction Foundation of Central South University of Forestry & Technology (104-0129).

References

B. A. McClure, V. Haring, P. R. Ebert et al., “Style self-incompatibility gene products of Nicotiana alata are ribonucleases,” Nature, vol. 342, no. 6252, pp. 955–957, 1989.

Y. Xue, R. Carpenter, H. G. Dickinson, and E. S. Coen, “Origin of allelic diversity in antirrhinum S locus RNases,” Plant Cell, vol. 8, no. 5, pp. 805–814, 1996.

T. Ishimizu, Y. Sato, T. Saito et al., “Identification and partial amino acid sequences of seven S-RNases associated with self-incompatibility of Japanese pear, Pyrus pyrifolice Nakai,” Journal of Biochemistry, vol. 120, no. 2, pp. 326–334, 1996.

H. Sassa, H. Hirano, T. Nishio, and T. Koba, “Style-specific self-compatible mutation caused by deletion of the S-RNase gene in Japanese pear (Pyrus serotina),” Plant Journal, vol. 12, no. 1, pp. 223–227, 1997.

T. Ishimizu, K. Inoue, M. Shimonaka, T. Saito, O. Terai, and S. Norioka, “PCR-based method for identifying the S-genotypes of Japanese pear cultivars,” Theoretical and Applied Genetics, vol. 98, no. 6-7, pp. 961–967, 1999.

C. Castillo, T. Takasaki, T. Saito, S. Norioka, and T. Nakanishi, “Cloning of the S₈-RNase (S₈-allele) of Japanese pear (Pyrus pyrifolia Nakai),” Plant Biotechnology, vol. 19, no. 1, pp. 1–6, 2002.

T. Takasaki, K. Okada, C. Castillo et al., “Sequence of the S₉-RNase cDNA and PCR-RFLP system for discriminating S₁-to S₉-allele in Japanese pear,” Euphytica, vol. 135, no. 2, pp. 157–167, 2004.

H. T. Kim, H. Y. Kakui, T. Koba, Y. Hirata, and H. Sassa, “Cloning of a new S-RNase and development of a PCR-RFLP system for the determination of the S-genotypes of Japanese pear,” Breeding Science, vol. 57, no. 2, pp. 159–164, 2007.

L. Bao, K. Chen, D. Zhang, Y. Cao, T. Yamamoto, and Y. Teng, “Genetic diversity and similarity of pear (Pyrus L.) cultivars native to East Asia revealed by SSR (simple sequence repeat) markers,” Genetic Resources and Crop Evolution, vol. 54, no. 5, pp. 959–971, 2007.

Y. M. Gong, C. W. Zong, H. N. Cao, et al., “Determination of S-genotypes of eleven pear cultivars growing in Yanbian region,” Journal of Fruit Science, vol. 28, no. 4, pp. 575–579, 2011.

W. Heng, S. L. Zhang, Y. Y. Zhang, et al., “Identification of S-genotypes of twelve pear cultivars by analysis of DNA sequence,” Acta Hortic Sinica, vol. 34, no. 4, pp. 853–858, 2007.

L. X. Wang, X. F. Tan, X. G. Li, et al., “Identification of S-genotypes of six pear cultivars and sequence analyses of S₄₅-RNase,” Journal of Fruit Science, vol. 27, no. 5, pp. 689–693, 2010.

T. N. Wuyun, X. F. Tan, X. G. Li, et al., “Isolation and identification of S₁₂-RNase gene of Pyrus bretschneideri,” Scientia Silvae Sinicae, vol. 42, no. 4, pp. 117–121, 2006.

X. F. Tan, L. Zhang, T. N. Wuyun et al., “Molecular identification of two new self-incompatible alleles (S-alleles) in Chinese pear (Pyrus bretschneideri),” Journal of Plant Physiology and Molecular Biology, vol. 33, no. 1, pp. 61–70, 2007.

L. Zhang, X. F. Tan, J. Zhou, et al., “Determination of S-genotypes of seven cultivars and identification of a novel S-RNase allele in Pyrus pyrifolia,” Scientia Silvae Sinicae, vol. 44, no. 10, pp. 42–48, 2008.

X. F. Tan, L. Zhang, D. Y. Yuan, et al., “Cloning and sequence analysis of full length cDNA encoding S₁₆-RNase from “Xuehua” pear (Pyrus bretschneideri Rehd.),” Journal of Zhejiang University, vol. 34, no. 2, pp. 149–157, 2008.

Y. L. Zeng, Identification of S-Genotype of Eli and the S-Alleles CDNA Cloning, Central South University of Forestry and Technology, Hunan, China, 2007.

A. Zisovich, R. A. Stern, and M. Goldway, “Identification of seven haplotype-specific SFBs in European pear (Pyrus communis) and their use as molecular markers,” Scientia Horticulturae, vol. 121, no. 1, pp. 49–53, 2009.

H. Kakui, T. Tsuzuki, T. Koba, and H. Sassa, “Polymorphism of SFBB-γ and its use for S genotyping in Japanese pear (Pyrus pyrifolia),” Plant Cell Reports, vol. 26, no. 9, pp. 1619–1625, 2007.

H. Sassa, H. Kakui, M. Miyamoto et al., “S locus F-box brothers: multiple and pollen-specific F-box genes with S haplotype-specific polymorphisms in apple and Japanese pear,” Genetics, vol. 175, no. 4, pp. 1869–1881, 2007.

H. Sassa, H. Kakui, and M. Minamikawa, “Pollen-expressed F-box gene family and mechanism of S-RNase-based gametophytic self-incompatibility (GSI) in Rosaceae,” Sexual Plant Reproduction, vol. 23, no. 1, pp. 39–43, 2010.

L. Zhang, X. F. Tan, J. Hu, et al., “Molecular cloning and characterization of PbSFBB13-gamma from Chinese white pear (Pyrus bretschneideri),” Scientia Silvae Sinicae, vol. 45, no. 11, pp. 36–43, 2009.

T. N. Wuyun, M. R. Bao, and T. Z. Li, “Identification and sequence analysis of new SFBB-gamma of Chinese pear,” Journal of Central South University of Technology, vol. 30, no. 8, pp. 5–10, 2010.

C. Xu, M. Li, J. Wu, et al., “Identification of a canonical SCF(SLF) complex involved in S-RNase-based self-incompatibility of Pyrus (Rosaceae),” Plant Molecular Biology, vol. 81, no. 3, pp. 245–257, 2013.

X. F. Tan, T. N. Wuyun, T. Nakanishi, et al., “Isolation and identification of S-RNase gene of the Pyrus bretschneideri (Chinese pear),” Journal of the Central South University of Forestry & Technology, vol. 25, no. 1, pp. 1–3, 2005.

F. J. Shao and S. F. Lu, “Genome-wide identification, molecular cloning, expression profiling and posttranscriptional regulation analysis of the Argonaute gene family in Salvia miltiorrhiza, an emerging model medicinal plant,” BMC Genomics, vol. 14, p. 512, 2013.

L. Zhang, Identification of S-Genotypes of Cultivars and Cloning of Self-Incompatibility-Associated Genes in Pear, Central South University of Forestry and Technology, Hunan, China, 2008.

Y. Kawata, F. Sakiyama, F. Hayashi, and Y. Kyogoku, “Identification of two essential histidine residues of ribonuclease T2 from Aspergillus oryzae,” European Journal of Biochemistry, vol. 187, no. 1, pp. 255–262, 1990.

C. Després, M. SabaElLeil, S. R. Rivard, D. Morse, and M. Cappadocia, “Molecular cloning of two Solanum chacoense S-alleles and a hypothesis concerning their evolution,” Sexual Plant Reproduction, vol. 7, no. 3, pp. 169–176, 1994.

K. Ushijima, H. Sassa, R. Tao et al., “Cloning and characterization of cDNAs encoding S-RNases from almond (Prunus dulcis): primary structural features and sequence diversity of the S-RNases in Rosaceae,” Molecular and General Genetics, vol. 260, no. 2-3, pp. 261–268, 1998.

W. Broothaerts, G. A. Janssens, P. Proost, and W. F. Broekaert, “cDNA cloning and molecular analysis of two self-incompatibility alleles from apple,” Plant Molecular Biology, vol. 27, no. 3, pp. 499–511, 1995.

R. Tao, H. Yamane, A. Sugiura, H. Murayama, H. Sassa, and H. Mori, “Molecular typing of S-alleles through identification, characterization and cDNA cloning for S-RNases in sweet cherry,” Journal of the American Society for Horticultural Science, vol. 124, no. 3, pp. 224–233, 1999.

H. Yaegaki, T. Shimada, T. Moriguchi, H. Hayama, T. Haji, and M. Yamaguchi, “Molecular characterization of S-RNase genes and S-genotypes in the Japanese apricot (Prunus mume Sieb.et Zucc.),” Sexual Plant Reproduction, vol. 13, no. 5, pp. 251–257, 2000.

H. Qiao, H. Wang, L. Zhao et al., “The F-box protein AhSLF-S2 physically interacts with S-RNases that may be inhibited by the ubiquitin/26S proteasome pathway of protein degradation during compatible pollination in antirrhinum,” Plant Cell, vol. 16, no. 3, pp. 582–595, 2004.