Skip Navigation


JXB Advance Access originally published online on July 9, 2008
Journal of Experimental Botany 2008 59(11):3169-3185; doi:10.1093/jxb/ern172
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Material
Right arrowOA All Versions of this Article:
59/11/3169    most recent
ern172v1
Right arrow E-letters: Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when E-letters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Tsukamoto, T.
Right arrow Articles by Iezzoni, A. F.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Tsukamoto, T.
Right arrow Articles by Iezzoni, A. F.
Agricola
Right arrow Articles by Tsukamoto, T.
Right arrow Articles by Iezzoni, A. F.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© 2008 The Author(s).
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/2.0/uk/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited. This paper is available online free of all access charges (see
http://jxb.oxfordjournals.org/open_access.html for further details)


RESEARCH PAPER

Genetic and molecular characterization of three novel S-haplotypes in sour cherry (Prunus cerasus L.)

Tatsuya Tsukamoto1 *, Daniel Potter2, Ryutaro Tao3, Cristina P. Vieira4, Jorge Vieira4 and Amy F. Iezzoni1,{dagger}

1Department of Horticulture, Michigan State University, East Lansing, MI 48824, USA
2Department of Plant Sciences, University of California, Davis, CA 95616–8780, USA
3Laboratory of Pomology, Graduate School of Agriculture, Kyoto University, Kyoto 606-8502, Japan
4Instituto de Biologia Molecular e Celular, University of Porto, 4150–180 Porto, Portugal

{dagger} To whom correspondence should be addressed. E-mail: iezzoni{at}msu.edu

Received 20 March 2008; Revised 22 May 2008 Accepted 27 May 2008


    Abstract
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Supplementary data
 References
 
Tetraploid sour cherry (Prunus cerasus L.) exhibits gametophytic self-incompatibility (GSI) whereby the specificity of self-pollen rejection is controlled by alleles of the stylar and pollen specificity genes, S-RNase and SFB (S haplotype-specific F-box protein gene), respectively. As sour cherry selections can be either self-compatible (SC) or self-incompatible (SI), polyploidy per se does not result in SC. Instead the genotype-dependent loss of SI in sour cherry is due to the accumulation of non-functional S-haplotypes. The presence of two or more non-functional S-haplotypes within sour cherry 2x pollen renders that pollen SC. Two new S-haplotypes from sour cherry, S33 and S34, that are presumed to be contributed by the P. fruticosa species parent, the complete S-RNase and SFB sequences of a third S-haplotype, S35, plus the presence of two previously identified sweet cherry S-haplotypes, S14 and S16 are described here. Genetic segregation data demonstrated that the S16-, S33-, S34-, and S35-haplotypes present in sour cherry are fully functional. This result is consistent with our previous finding that ‘hetero-allelic’ pollen is incompatible in sour cherry. Phylogenetic analyses of the SFB and S-RNase sequences from available Prunus species reveal that the relationships among S-haplotypes show no correspondence to known organismal relationships at any taxonomic level within Prunus, indicating that polymorphisms at the S-locus have been maintained throughout the evolution of the genus. Furthermore, the phylogenetic relationships among SFB sequences are generally incongruent with those among S-RNase sequences for the same S-haplotypes. Hypotheses compatible with these results are discussed.

Key words: Prunus cerasus, self-incompatibility, SFB, S-RNase


    Introduction
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Supplementary data
 References
 
Gametophytic self-incompatibility (GSI) is a common genetic mechanism that promotes outcrossing in flowering plants (de Nettancourt, 2001). In GSI, self-incompatibility (SI) is determined by a single multi-allelic locus, called the S-locus, which contains a minimum of two genes, one (stylar-S) controlling stylar specificity and the other (pollen-S) controlling pollen specificity of the SI reaction. The stylar-S gene in three plant families, the Solanaceae, Plantaginaceae, and Rosaceae encodes a ribonuclease (S-RNase; Anderson et al., 1986; McClure et al., 1989; Sassa et al., 1992; Xue et al., 1996), which is expressed in the pistil and specifically degrades the RNA of incompatible pollen (McClure et al., 1990). The pollen-S gene encodes an F-box protein named S-locus F-box protein (SLF) in Antirrhinum (Lai et al., 2002), Petunia inflata (Sijacic et al., 2004), and Prunus mume (Entani et al., 2003), or S haplotype-specific F-box protein (SFB) in Prunus dulcis, P. mume, P. avium, P. spinosa, and P. cerasus (Ushijima et al., 2003; Yamane et al., 2003b; Ikeda et al., 2004a; Nunes et al., 2006). Despite having similar or even identical names in Solanaceae and Prunus the pollen gene is not orthologous (Wheeler and Newbigin, 2007).

Within Prunus (Rosaceae), cherry represents a natural diploid–tetraploid series with the tetraploid sour cherry (P. cerasus) arising through hybridization between sweet cherry (P. avium) and the tetraploid ground cherry (P. fruticosa) (Olden and Nybom, 1968). Like sweet cherry, sour cherry exhibits an S-RNase-based GSI system (Yamane et al., 2001; Hauck et al., 2002; Tobutt et al., 2004), however, in contrast to sweet cherry, natural sour cherry selections include both SI and self-compatible (SC) types (Redalen, 1984; Lansari and Iezzoni, 1990). This genotype-dependent loss of SI in sour cherry indicates that genetic changes, and not polyploidy per se, cause the breakdown of SI. The genetic switch from SI to SC in sour cherry results from the accumulation of non-functional S-haplotypes according to the ‘one-allele-match model’ (Hauck et al., 2006b). In this model, the match between a functional pollen-S in the 2x pollen and its cognate functional S-RNase in the style, results in an incompatible reaction. A similar reaction would occur regardless of whether the pollen contained a single functional pollen-S gene or two different pollen-S genes. The absence of any functional match results in a compatible reaction. Thus for successful self-fertilization, the 2x pollen must contain two non-functional S-haplotypes. Recently, Huang et al. (2008) reported competitive interaction in a SC selection of tetraploid Prunus pseudocerasus, raising the possibility that the SC mechanism between these two tetraploid Prunus species could be different. However, although the data in Huang et al. (2008) is consistent with hetero-allelic pollen being SC, homo-allelic pollen (e.g. S1S1, S5S5, or S7S7) was not shown to be successful in a compatible cross and unsuccessful in an incompatible cross. Therefore, it is possible that the SC in P. pseudocerasus could be caused by mutations in other genes critical for the SI reaction.

Six S-haplotypes present in sweet cherry (S1, S4, S6, S9, S12, and S13) have been shown to be present in sour cherry as well. However, three of these S-haplotypes (S1, S6, and S13) also have non-functional variants in sour cherry that have lost pollen and/or stylar function (Yamane et al., 2003a; Hauck et al., 2006a, b; Tsukamoto et al., 2006). Loss of function in these non-functional S-haplotypes was due to structural alterations of the S-RNase, SFB or S-RNase upstream sequences. Sour cherry also possesses S-haplotypes (S26, S36a, S36b, S36b2, and S36b3) that were presumably derived from the other species parent, P. fruticosa, as these S-haplotypes have not been identified in sweet cherry (Hauck et al., 2006b; T Tsukamoto et al., unpublished data). The extensive sour cherry germplasm collection at Michigan State University (Iezzoni, 2005) provides an excellent resource for the identification of previously undiscovered S-haplotypes for which information about their functionality would aid in the breeding of SC types. Therefore a germplasm survey was undertaken to search for novel S-haplotypes in sour cherry and to determine the functionality of these S-haplotypes. Here two new S-haplotypes identified in sour cherry are described that are presumed to be contributed by the P. fruticosa species parent. The first complete S-RNase and SFB sequences for a third S-haplotype that was previously identified in sour cherry are also reported. Furthermore, genetic segregation data presented demonstrates that all three of these sour cherry S-haplotypes are fully functional. In addition, the S14- and S16-haplotypes have been identified in sour cherry. Phylogenetic analyses of S-RNase and SFB sequences are used to investigate patterns of evolution of S-haplotypes in Prunus.


    Materials and methods
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Supplementary data
 References
 
Plant material
One sweet cherry cultivar ‘Hedelfingen’ (S3S5) and 17 sour cherry cultivars were used: ‘Cigány 59’, ‘Crisana’, ‘Englaise Timpurii’, ‘Erdi Botermo’, ‘Erdi Jubileum’, ‘Erdi Nagygyumolcsu’, ‘Meteor’, ‘Montmorency’, ‘Pandy 38’, ‘Pandy 114’, ‘Rheinische Schattenmorelle’, ‘Surefire’, ‘Tamaris’, ‘Tarina’, ‘Tschernokorka’, ‘Újfehértói furtos’, and ‘MSU III 18 (12)’. These cultivars were grown at the Michigan State University Clarksville Horticultural Experimental Station, Clarksville, MI, USA and a commercial orchard in Suttons Bay, MI, USA. To determine the inheritance of the S-haplotypes, self-pollinated progeny of ‘Meteor’, ‘Montmorency’ and ‘Tamaris’ plus cross-pollinated progeny between ‘Újfehértói furtos’x‘Surefire’ were obtained.

DNA isolation
Young unfolded leaves were collected in the spring, frozen in liquid nitrogen, lyophilized, and stored at –20 °C. Genomic DNA was isolated from lyophilized leaves according to the method of Ikeda et al. (2004b). Extracted leaf DNA was treated with RNase A (Roche, Mannheim, Germany). To genotype the self-pollinated progeny of ‘Meteor’, ‘Montmorency’, and ‘Tamaris’, DNA was extracted from mature seed from which the testa was removed using the procedure of Hauck et al. (2006b). Extracted seed DNA was treated with RNase A (Roche).

PCR amplification
Total DNA was isolated from the cherry selections and used as template DNA for PCR. PCR procedures were identical to those used by Tao et al. (1999). The S-RNase gene specific primer set, Pru-C2 and PCE-R (Tao et al., 1999; Yamane et al., 2001) that correspond to the previously identified C2 and C3 conserved regions (Ushijima et al., 1998), respectively, were used. This primer pair can differentiate among most S-RNase alleles based on polymorphisms in the length of the second intron in the Prunus S-RNase. However, this primer pair cannot amplify S35-RNase. EM-PC2consFD and EM-PC5consRD (Sutherland et al., 2004) were used to amplify S35-RNase. In addition, the primer pair of Pru-C2 and PCE-R cannot differentiate S36a from S36b, S36b2 and S36b3. Instead, these variants of the S36-haplotype were differentiated using the following primer pairs for detection of the S36a-haplotype [PcS36ab-F (5'-GCTAGCCAACCACTTTTACG-3’) and PcS36a-spR (5'-GAAACCCACATGATACAAACTG-3’)] and detection of the S36b-, S36b2-, and S36b3-haplotypes [PcS36ab-F and PcS36b-spR (5'-ATACATTGTAGGCCAGTCTGTG-3’)]. PCR products were run on 2% agarose gels and the DNA bands were visualized by ethidium bromide staining. Furthermore, PCR products were purified using the QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany), ligated into the pGEM-T Easy vector (Promega, Madison, WI, USA), and transformed Escherichia coli JM109 (Promega). Plasmid DNA was prepared using Wizard Plus Minipreps DNA Purification Kit (Promega) and their sequences were determined as described below.

Construction and screening of genomic libraries
Fosmid libraries were constructed using the Copy Control Fosmid library production kit (Epicentre Technologies, Madison, WI, USA). Fosmid libraries from ‘Meteor’, ‘Montmorency’, and ‘Tamaris’ were constructed to clone the S-RNase and SFB alleles from the S-haplotypes being investigated. The fosmid libraries were screened at 60 °C with a mixture of DIG-dUTP-labelled S6-RNase and SFB6 probes, as previously described (Ushijima et al., 2001). The DIG-labelled S6-RNase and SFB6 probes were obtained by PCR-labelling using the PCR DIG Probe Synthesis Kit (Roche) with Pru-C2 and PCE-R primers (Tao et al., 1999; Yamane et al., 2001), and SFB-C1F (Ikeda et al., 2004a) and SFB-C4R primers (Yamane et al., 2003c). Fosmid DNAs from positive clones were prepared using the Wizard Plus Minipreps DNA Purification Kit (Promega). The S-haplotype of each positive fosmid clone was determined by PCR with the S-RNase consensus primer pair (Pru-C2 and PCE-R). Positive clones were also analysed by PCR with the SFB consensus primer pair (SFB-C1F and SFB-C2R) (Ikeda et al., 2004a) to check if the clone has SFB. S-RNase and SFB allele-specific primer pairs were also used to identify the S-haplotypes. For the S35-haplotype in ‘Montmorency’ (S6S13mS35S36a), positive clones were obtained by using a DIG-dUTP-labelled S35-RNase probe. Then more positive clones obtained by hybridization with a mixture of DIG-dUTP-labelled SFB6 and SFB36a probes were subjected to PCR with the SFB consensus primer (SFB-C1F and SFB-C2R) to confirm that the clones contain SFB. Then positive clones were analysed by PCR with SFB allele-specific (SFB6, SFB13, and SFB36a/SFB36b) primer pairs: PaSFB6-F and PaSFB6-R (Ikeda et al., 2005), DdeS13-F (Tsukamoto et al., 2008) and SFB13-spR (Tsukamoto et al., 2006), and PcSFB36ab-F (5'-GGCGGTCGATCCTGATGAC-3') and PcSFB36ab-R (5'-TGTCCGATAACAGCTCCGG-3'), respectively. The positive clones from which a fragment could be amplified with the primer pair SFB-C1F and SFB-C2R, but not amplified with SFB6-, SFB13-, and SFB36a/SFB36b-specific primers, were determined to contain SFB35 and one positive clone was sequenced.

DNA sequencing
DNA sequencing was carried out by using ABI PRISM 3100 Genetic Analyser at the Michigan State University Research Technology Support Facility. The plasmid clones were sequenced by using SP6 and T7 primers. The fosmid clones were sequenced by primer walking using Pru-T2, Pru-C2, Pru-C2R, PCE-F, PCE-R, Pru-C4R, and Pru-C5 (Tao et al., 1999; Yamane et al., 2001; Tsukamoto et al., 2006) for S-RNase, and SFB-C1F, SFB-C2R, SFB-C5F, and FB3R (Ikeda et al., 2005) for SFB. The EM-PC5consRD primer (Sutherland et al., 2004) was also used to sequence S35-RNase.

Phylogenetic and variability analysis
For phylogenetic analyses, available nucleotide sequences for S-RNase and SFB from Prunus armeniaca, P. avium, P. cerasus, P. domestica, P. dulcis, P. mume, P. salicina, P. spinosa, and P. tenella S-haplotypes were assembled (see Fig. 7). Only sequences whose predicted amino acid sequences covered at least 75% of the average length of available complete sequences of S-RNase (227 amino acids) and SFB (375 aminio acids) from the two cherry species (P. avium and P. cerasus) were included. Nucleotide sequences were translated into deduced amino acid sequences using DAMBE (Xia and Xie, 2001); the amino acid sequences were then aligned using ClustalX (Thompson et al., 1997) and adjusted manually. Finally, nucleotide sequences were aligned to the amino acid alignments with DAMBE. Analyses of DNA variability were performed using DnaSP 4.1 (Rozas et al., 2003). Phylogenetic analyses of the aligned amino acid and nucleotide sequences based on maximum parsimony were implemented in PAUP* (Swofford, 2002) with heuristic searches using the TBR branch-swapping algorithm and 1000 random taxon addition replicates and maxtrees allowed automatic increases as necessary. Relative support for clades was assessed using 1000 bootstrap replicates with 10 random taxon addition replicates per bootstrap replicate and maxtrees set at 100.


Figure 7
View larger version (29K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 7. Phylogenetic trees based on phylogenetic analysis of nucleotide sequences of S-RNase (left) and SFB (right) alleles for species of Prunus. Numbers above, below, or adjacent to branches are bootstrap support values greater than 50%; asterisks indicate branches that collapsed in the strict consensus tree from each analysis. The positions of PcS33, PcS34, and PcS35 are indicated, respectively, by open boxes, lightly shaded boxes, and darkly shaded boxes. Left: Single most parsimonious tree (l=2366, ci excluding autapomorphies=0.3271, ri=0.5263) based on S-RNase alleles. The nucleotide sequences for three novel P. cerasus (PcS-RNase) alleles presented in this study were aligned with 12 S-RNase alleles from P. avium (PaS1-RNase, AB028153 [GenBank] ; PaS2-RNase, AB010304 [GenBank] ; PaS3-RNase, AB010306 [GenBank] ; PaS4-RNase, AB028154 [GenBank] ; PaS5-RNase, AJ298314 [GenBank] ; PaS6-RNase, AB010305 [GenBank] ; PaS7-RNase, EU035974 [GenBank] ; PaS12-RNase, AY259115 [GenBank] ; PaS13-RNase, DQ385842 [GenBank] ; PaS23-RNase, AY259114 [GenBank] ; PaS24-RNase, AY259112 [GenBank] ; PaS25-RNase, AY259113 [GenBank] ); one from P. cerasus (PcS26-RNase, EU035975 [GenBank] ); three from P. domestica (PdoS5-RNase, AM746946 [GenBank] ; PdoS6-RNase, AM746947 [GenBank] ; PdoS9-RNase, AM746948 [GenBank] ); nine from P. dulcis (PdSa-RNase, AB026836 [GenBank] ; PdSb-RNase, AB011469 [GenBank] ; PdSc-RNase, AB011470 [GenBank] ; PdSd-RNase, AB011471 [GenBank] ; PdSk-RNase, AB252409 [GenBank] ; PdSm-RNase, DQ099895 [GenBank] ; PdSn-RNase, DQ093825 [GenBank] ; PdS11-RNase, AM231660 [GenBank] ; PdS12-RNase, AM746949 [GenBank] ); two from P. mume (PmS1-RNase, AB101438 [GenBank] ; PmS7-RNase, AB101439 [GenBank] ); four from P. armeniaca (ParS1-RNase, AY587561 [GenBank] ; ParS2-RNase, AY587562 [GenBank] ; ParS4-RNase, AY587564 [GenBank] ; ParS17-RNase, EU516388 [GenBank] ); three from P. cerasifera (PcsfS3-RNase, AM746943 [GenBank] ; PcsfS9-RNase, AM746944 [GenBank] ; PcsfS10-RNase, AM746945 [GenBank] ); ten from P. salacina (PsSa-RNase, AB252411 [GenBank] ; PsSb-RNase, AB252413 [GenBank] ; PsSc-RNase, AB084102 [GenBank] ; PsSd-RNase, AB084103 [GenBank] ; PsSe-RNase, AB280693 [GenBank] ; PsSf-RNase, DQ512911 [GenBank] ; PsSg-RNase, AM746950 [GenBank] ; PsSh-RNase, DQ512914 [GenBank] ; PsS7-RNase, AY781290 [GenBank] ; PsS8-RNase, DQ512913 [GenBank] ); seven from P. spinosa (PspS8-RNase, DQ677587 [GenBank] ; PspS9-RNase, DQ677588 [GenBank] ; PspS10-RNase, DQ677589 [GenBank] ; PspS12-RNase, DQ677590 [GenBank] ; PspS3-1-RNase, DQ677584 [GenBank] ; PspS3-2-RNase, DQ677585 [GenBank] ; PspS7-1-RNase, DQ677586 [GenBank] ); and nine from P. tenella (PtS1-RNase, DQ983373 [GenBank] ; PtS2-RNase, DQ983374 [GenBank] ; PtS3-RNase, DQ983375 [GenBank] ; PtS4-RNase, DQ983363 [GenBank] ; PtS5-RNase, DQ983364 [GenBank] ; PtS6-RNase, DQ983365 [GenBank] ; PtS7-RNase, DQ983366 [GenBank] ; PtS8-RNase, DQ983367 [GenBank] ; PtS9-RNase, DQ983370 [GenBank] ). Data set contained 63 taxa and 744 characters, of which 220 were constant, 129 were variable but uninformative, and 395 were parsimony-informative. Right: One of 65 most parsimonious trees (l = 2968, ci excluding autapomorphies = 0.3913, ri = 0.4860) from phylogenetic analysis of nucleotide sequences of SFB alleles. The nucleotide sequences for three novel P. cerasus (PcSFB) alleles presented in this study were aligned with coding sequences for 12 SFB alleles from P. avium (PaSFB1, AY805048 [GenBank] ; PaSFB2, AB111519 [GenBank] ; PaSFB3, AB096857 [GenBank] ; PaSFB4, AB111521 [GenBank] ; PaSFB5, AB111520 [GenBank] ; PaSFB6, AB096858 [GenBank] ; PaSFB7, EU035976 [GenBank] ; PaSFB9, DQ422809 [GenBank] ; PaSFB10, AY805053 [GenBank] ; PaSFB12, AY805054 [GenBank] ; PaSFB13, DQ385844 [GenBank] ; PaSFB16, AY805056 [GenBank] ); two from P. cerasus (PcSFB1, DQ827715 [GenBank] ; PcSFB26, EU035977 [GenBank] ); three from P. domestica (PdoSFB5, AM746955 [GenBank] ; PdoSFB6, AM746956 [GenBank] ; PdoSFB9, AM746957 [GenBank] ); seven from P. dulcis (PdSFBa, AB092966 [GenBank] ; PdSFBb, AB092967 [GenBank] ; PdSFBc, AB079776 [GenBank] ; PdSFBd, AB081648 [GenBank] ; PdSFBk, AB252408 [GenBank] ; PdSFB11, EF061758 [GenBank] ; PdSFB12, AM746959 [GenBank] ); three from P. mume (PmSFB1, AB101440 [GenBank] ; PmSFB7, AB101441 [GenBank] ; PmSFB9, AB092645 [GenBank] ); four from P. armeniaca (ParSFB1, AY587563 [GenBank] ; ParSFB2, AY587562 [GenBank] ; ParSFB4, AY587565 [GenBank] ; ParSFB17, EU516388 [GenBank] ); three from P. cerasifera (PcsfSFB3, AM746952 [GenBank] ; PcsfSFB9, AM746953 [GenBank] ; PcsfSFB10, AM746954 [GenBank] ); nine from P. salicina (PsSFBa, AB252410 [GenBank] ; PsSFBb, AB252412 [GenBank] ; PsSFBc, DQ849084 [GenBank] ; PsSFBd, AM746962 [GenBank] ; PsSFBe, AB280794 [GenBank] ; PsSFBf, DQ849089 [GenBank] ; PsSFBg, AM746963 [GenBank] ; PsSFBh, DQ849118 [GenBank] ; PsSFB7, DQ849085 [GenBank] ); seven from P. spinosa (PspSFB8, DQ677587 [GenBank] ; PspSFB9, DQ677588 [GenBank] ; PspSFB10, DQ677589 [GenBank] ; PspSFB12, DQ677598 [GenBank] ; PspSFB3-1, DQ677616 [GenBank] ; PspSFB3-2, DQ677615 [GenBank] ; PspSFB7-1, DQ677595 [GenBank] ), and one from P. tenella (PtSFB8, DQ983369 [GenBank] ). Data set contained 54 taxa and 1161 characters, of which 337 were constant, 230 were variable but uninformative, and 594 were parsimony-informative.

 
Several approaches were used to test for significant incongruence between the topologies of phylogenetic trees supported by the SFB and S-RNase data sets, respectively, and all of these approaches were used for both the amino acid and the nucleotide sequences for each gene. First, a combined data set was constructed including only the 49 haplotypes for which sequences of both determinants were available. The partition homogeneity test, implemented in PAUP* with 1000 test replicates and heuristic searches using the TBR branch-swapping algorithm and 10 random taxon addition replicates per test replicate with maxtrees set to 100, was used to test for significant conflict between the S-RNase and SFB partitions within each data set. Second, each partition was analysed separately and all of the most parsimonious trees were saved to a single tree file. The Kishino-Hasegawa (K-H), Templeton, and winning-sites tests were implemented in PAUP* to test whether or not the topologies produced by the two data partitions were significantly different. Third, because bootstrap analyses showed that many relationships were only weakly supported by each of the two data partitions, constraint trees were constructed in which only groups supported with bootstrap values of 80% or more by each of the data partitions were resolved. Four such constraint trees were constructed; one each for SFB and S-RNase for the combined amino acid data set and one each for SFB and S-RNase for the combined nucleotide data set. Each data partition was analysed without constraints and the best trees were saved to a file. Then each data partition was analysed under the constraint corresponding to the well-supported groups from the other partition and the best trees were saved to the same file. The K-H, Templeton, and winning-sites tests were implemented in PAUP* to test whether or not the constrained trees were significantly longer than the unconstrained trees. In order to test whether or not the phylogenetic relationships among SFB and S-RNase alleles were significantly different from the phylogenetic relationships among species of Prunus based on other evidence, each data set was re-analysed with various topological constraints enforced (see Results). For each data set, all of the most parsimonious trees (MPT) from unconstrained and constrained analyses were saved to a single tree file. The K-H, Templeton, and winning-sites tests were implemented in PAUP* to test whether or not the constrained trees were significantly longer than the unconstrained trees.


    Results
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Supplementary data
 References
 
Identification and cloning of novel S-haplotype genes (S-RNase and SFB) from sour cherry
Genomic PCR with the S-RNase consensus primer pair Pru-C2 and PCE-R (Tao et al., 1999; Yamane et al., 2001) revealed three amplification products not previously characterized in three sour cherry cultivars (Fig. 1A). Fragments of ~ 480 bp, ~ 420 bp, and ~ 850 bp were identified in ‘Englaise Timpurii’, ‘Meteor’, and ‘Tamaris’, respectively.


Figure 1
View larger version (34K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 1. PCR amplification for S-RNase alleles of 17 sour cherry selections. (A) Genomic DNA was amplified by PCR with Pru-C2 (Tao et al., 1999) and PCE-R (Yamane et al., 2001) primer set. (B) Genomic DNA was amplified by PCR with EM-PC2consFD and EM-PC5consRD (Sutherland et al., 2004) primer set. PCR products were separated on 2% agarose gels and detected with ethidium bromide staining. The colour of black and white is inverted in this image. The asterisks indicate the band of PCR product of S35-RNase. M, 123 bp DNA ladder (Invitrogen, Carlsbad, CA, USA). Lane abbreviations are: C59, ‘Cigány 59’; Cri, ‘Crisana’; ET, ‘Englaise Timpurii’; EB, ‘Erdi Botermo’; EJ, ‘Erdi Jubileum’; EN, ‘Erdi Nagygyumolcsu’; Met, ‘Meteor’; Mon, ‘Montmorency’; P38, ‘Pandy 38’; P114, ‘Pandy 114’; RS, ‘Rheinische Schattenmorelle’; Sur, ‘Surefire’; Tam, ‘Tamaris’; Tar, ‘Tarina’; Tsc, ‘Tschernokorka’; UF, ‘Újfehértói furtos’; III 18 (12), ‘MSU III 18 (12)’. (This figure is available in colour at JXB online).

 
The nucleotide sequence of the 481 bp S-RNase PCR product from ‘Englaise Timpurii’ revealed high similarity (99.6%) with the P. avium S14-RNase partial sequence (GenBank accession no. AJ635277 [GenBank] ; Sonneveld et al., 2003) and the P. avium S23-RNase complete sequence (GenBank accession no. AY259114 [GenBank] ; Wünsch and Hormaza, 2004; the two sequences are identical). The ‘Englaise Timpurii’ 481 bp S-RNase nucleotide sequence differed from the P. avium S14-/S23-RNase by only two base pairs within the second intron (data not shown). Therefore, the 481 bp PCR S-RNase product from ‘Englaise Timpurii’ was considered to be the S14-RNase, as the P. avium S14- and S23-RNases likely code the same specificity.

BlastN of the 424 bp partial nucleotide sequence of the S-RNase PCR product from ‘Meteor’ revealed high homology (98.3%) with the partial S10-RNase of Japanese apricot (P. mume) sequence (GenBank accession no. DQ011150 [GenBank] ; sequence upstream of the C2 conserved region and downstream of the RC4 conserved region is lacking). This novel S-haplotype was named S33 to follow the previously named functional cherry S-haplotypes (S1 to S7, S9, S10, S12 to S14, S16 to S32) (for a review see Vaughan et al., 2008).

Two S33 clones were obtained by screening a ‘Meteor’ fosmid library. Both clones contained the S-haplotype genes, S33-RNase and SFB33, and the complete nucleotide sequences of the S33-RNase (GenBank accession no. EU054325 [GenBank] ) and SFB33 (GenBank accession no. EU054328 [GenBank] ) were obtained by sequencing the genes from these two clones. The predicted amino acid sequence of the S33-RNase consisted of 238 residues and was aligned with that of functional S-RNases present in sour cherry (Fig. 2). When the nucleotide sequences of the coding regions of the sour cherry S33-RNase and P. mume S10-RNase were compared, there were four synonymous and one non-synonymous differences. In the second intron, there were four nucleotide differences plus two indels (148 bp and 1 bp long, respectively; see Supplementary Fig. 1 at JXB online).


Figure 2
View larger version (55K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 2. Amino acid sequence alignment of three novel S-RNases obtained from sour cherry and that of other functional S-RNases from sour cherry. The alignment was generated by DNASIS version 3.5 (Hitachi Software Engineering Co. Ltd., Tokyo, Japan). Gaps are marked by dashes. Conserved amino acids are shown on a darkened background. The five conserved regions, C1, C2, C3, RC4, and C5 (Ushijima et al., 1998) are marked with solid boxes, and the hypervariable region, RHV (Ushijima et al., 1998) reported in the rosaceous S-RNases, is marked with a dotted box. The positions and directions of the four consensus primers used in genomic PCR are indicated by arrows.

 
SFB33 consisted of 376 amino acid residues that exhibited the characteristic F-box and variability patterns of previously identified SFBs (Fig. 3). S33-RNase and SFB33 specific primer pairs were designed (Table 1) that amplified the allele specific fragments of 819 bp and 860 bp, respectively (Fig. 4A, B).


Figure 3
View larger version (74K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 3. Amino acid sequence alignment of three novel SFBs obtained from sour cherry and that of other functional SFBs from sour cherry. The alignment was generated by DNASIS version 3.5 (Hitachi Software Engineering Co. Ltd., Tokyo, Japan). Gaps are marked by dashes. Conserved amino acids are shown on a darkened background. The locations of the F-box motif, V1, V2, HVa, and HVb (Ikeda et al., 2004a), and Vn (Nunes et al., 2006) are indicated by solid boxes.

 

View this table:
[in this window]
[in a new window]

 
Table 1. Preferred PCR primer pairs to amplify S-RNase and SFB alleles of S33-, S34-, and S35-haplotypes identified in sour cherry

 

Figure 4
View larger version (51K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 4. PCR amplification with S33-, S34-, and S35-allele specific primer pair for S-RNase and SFB in 17 sour cherry selections. PCR products were separated on 2% agarose gel and detected with ethidium bromide staining. M, 123 bp DNA ladder (Invitrogen, Carlsbad, CA, USA). Lane abbreviations are: C59, ‘Cigány 59’; Cri, ‘Crisana’; ET, ‘Englaise Timpurii’; EB, ‘Erdi Botermo’; EJ, ‘Erdi Jubileum’; EN, ‘Erdi Nagygyumolcsu’; Met, ‘Meteor’; Mon, ‘Montmorency’; P38, ‘Pandy 38’; P114, ‘Pandy 114’; RS, ‘Rheinische Schattenmorelle’; Sur, ‘Surefire’; Tam, ‘Tamaris’; Tar, ‘Tarina’; Tsc, ‘Tschernokorka’; UF, ‘Újfehértói furtos’; III 18 (12), ‘MSU III 18 (12)’.

 
BlastN with the 868 bp partial nucleotide sequence of the S-RNase PCR product from ‘Tamaris’ revealed high homology with the S1-RNase obtained from sweet cherry (GenBank accession no. AB031815 [GenBank] ), the P. tenella S8-RNase (GenBank accession no. DQ983367 [GenBank] ), and partial sequences of the P. dulcis S11-RNase (GenBank accession no. AM231660 [GenBank] ) and the P. domestica S5-RNase (GenBank accession no. AM746946 [GenBank] ). This novel S-haplotype was named S34 (Fig. 5; Table 2). Sour cherry S34-RNase falls in this previously described group of S-RNases where the sweet cherry S1-RNase was found to be identical to the P. tenella S8-RNase and to differ from the P. dulcis S11-RNase by just one amino acid (Surbanovski et al., 2007).


Figure 5
View larger version (21K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 5. The deduced amino acid sequence alignment for the Prunus cerasus S34-RNase (PcS34-RNase), P. domestica S5-RNase (PdoS5-RNase), P. avium S1-RNase (PaS1-RNase), P. dulcis S11-RNase (PdS11-RNase), and P. tenella S8-RNase (PtS8-RNase). The asterisks indicate the five amino acid residues that are different between P. cerasus S34-RNase and P. domestica S5-RNase. Conserved nucleotides are shown on a darkened background. The five conserved regions, C1, C2, C3, RC4, and C5 (Ushijima et al., 1998) are marked with solid boxes, and the hypervariable region, RHV (Ushijima et al., 1998) reported in the rosaceous S-RNases, is marked with a dotted box.

 

View this table:
[in this window]
[in a new window]

 
Table 2. Identities of the derived amino acid sequences among the cherry S-locus genes whose complete sequence are available

 
One S34 clone was obtained by screening a ‘Tamaris’ fosmid library. This clone contained both of the S34-RNase and SFB34 permitting the complete nucleotide sequencing of the S34-RNase (GenBank accession no. EU054326 [GenBank] ) and SFB34 (GenBank accession no. EU054329 [GenBank] ). The predicted amino acid sequence of S34-RNase consisted of 226 residues and was aligned with that of functional S-RNases obtained from sour cherry cultivars (Fig. 2). A comparison of the S34-RNase with a partial sequence of the P. domestica S5-RNase (PdoS5-RNase lacks 30 and 27 amino acids at the N- and C-terminals, respectively) revealed high homology (97%) and only five amino acid differences (Fig. 5). When PcS34-RNase was compared with P. avium S1-RNase, 15 amino acid differences were identified (Fig. 5). Nevertheless, the sweet cherry S1-RNase specific primer pair (Sonneveld et al., 2001) amplified a product for the S34-RNase. However, the amplification product from the S34-RNase was the expected size of ~850 bp instead of 615 bp as for the P. avium S1-RNase (data not shown). The second intron of the S34-RNase is 238 bp longer than that of the P. avium S1-RNase (GenBank accession nos AB031815 [GenBank] and AB028153 [GenBank] ; see Supplementary Fig. 2 at JXB online).

SFB34 was composed of 376 amino acid residues that exhibited the characteristic variability patterns of previously identified SFBs (Fig. 3). Due to the high homology of the S34-RNase with the P. avium S1-RNase, an alignment was performed with their respective SFBs. The predicted amino acid sequence of P. cerasus SFB34 was aligned with three cherry SFB1 sequences: (i) the P. cerasus SFB1 obtained from ‘Pandy 114’ (GenBank accession no. DQ827715 [GenBank] ), (ii) the P. avium SFB1 cloned from ‘Skeena’ (GenBank accession no. AY805048 [GenBank] ), and (iii) the P. avium SFB1 sequence obtained from ‘Seneca’ (GenBank accession no. AB111518 [GenBank] ) and P. dulcis SFB11, and P. tenella SFB8 (Fig. 6). The P. cerasus ‘Pandy 114’ SFB1 and P. avium cv. ‘Skeena’ SFB1 sequences are identical, however the ‘Skeena’ sequence is not complete as 10 amino acid residues at the C-terminal are not determined. The amino acid sequence of the ‘Seneca’ SFB1 is slightly different from that of SFB1 of ‘Pandy 114’ and ‘Skeena’ (Fig. 6). Unlike their S-RNases, the predicted amino acid sequences for the sour cherry SFB34 and SFB1 differed by 69 amino acids (Fig. 6) and shared only 81.6% identity (Table 2). Similarly, the amino acid sequence of the P. cerasus SFB34 had only 81.0% homology with that of the identical P. tenella SFB8 and P. dulcis SFB11. By comparison, the P. cerasus SFB34 and the partial sequence of the P. domestica SFB5 (GenBank accession number. AM746955 [GenBank] ; PdoSFB5 is lacking 12 and 33 amino acid residues at the N- and C-terminals, respectively) differed at 15 amino acid residues (95.5% homology), compared with 69 different residues between SFB34 and SFB1 (Fig. 6).


Figure 6
View larger version (49K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 6. The deduced amino acid sequence alignment for the Prunus cerasus SFB34 (PcSFB34), P. domestica SFB5 (PdoSFB5), P. avium SFB1 obtained from ‘Skeena’ (PaSFB1-Skn), P. avium SFB1 obtained from ‘Seneca’ (PaSFB1-Snc), P. cerasus SFB1 obtained from ‘Pandy 114’ (PcSFB1), P. tenella SFB8 (PtSFB8), and P. dulcis SFB11 (PdSFB11). The asterisks indicate the 15 amino acid residues that are different between P. cerasus SFB34 and P. domestica SFB5. Conserved nucleotides are shown on a darkened background. The locations of the F-box motif, V1, V2, HVa, and HVb (Ikeda et al., 2004a), and Vn (Nunes et al., 2006) are indicated by solid boxes.

 
To compare the S34- and S1-haplotypes further, a partial sequence of the intergenic region between S-RNase and SFB of P. cerasus S34 was obtained. However, the intergenic sequence for the S34-haplotype was extremely divergent compared with that from the P. avium S1, P. dulcis S11, and P. tenella S8 and did not permit a conclusive alignment (see Supplementary Fig. 3 at JXB online).

To aid in the confirmation of the S34-haplotype in future germplasm surveys, S34-RNase and SFB34 specific primer pairs were designed (Table 1) that amplified allele specific fragments of expected size (898 bp and 714 bp, respectively) (Fig. 4C, D). The S1-RNase specific primer pair (Sonneveld et al., 2001) amplified fragments for the S1-RNase (~615 bp) and S34-RNase (~850 bp) (data not shown), but the S34-RNase specific primer pair only amplified a product from the S34-RNase, not the S1-RNase (Fig. 4C).

To identify all the S-haplotypes in ‘Tamaris’, genomic PCR was performed with all available S-RNase specific primer pairs. A fragment of ~ 680 bp was amplified with the S16-RNase specific primer pair (Sonneveld et al., 2003). A comparison of the nucleotide sequence from the S16-RNase PCR product from ‘Tamaris’ with the P. avium S16-RNase revealed only one base pair difference in the second intron. Therefore ‘Tamaris’ was considered to have the S16-haplotype.

Identification and cloning of the S35-RNase and SFB35
Previously an S-haplotype survey in sour cherry cultivars was carried out using RFLP analysis, and S-RNase based PCR was also carried out using the consensus S-RNase gene-specific primer set, Pru-C2 and PCE-R (Tao et al., 1999; Yamane et al., 2001). Using these two approaches, only three different S-haplotypes were characterized in many cultivars, including the landrace sour cherry cultivar ‘Pandy’ (syn. ‘Köröser’, ‘Crisana’) and ‘Montmorency’. Yet, genetic studies with many of these cultivars led to the conclusion that none of the three S-haplotypes identified in each cultivar was present in a double dose (Hauck et al., 2006b). For example, the cultivar ‘Újfehértói furtos’ has the S-haplotypes S1', S4, and S36b of which S1' and S36b are non-functional. Therefore, it is possible that ‘Újfehértói furtos’ could contain two copies for either of these two non-functional S-haplotypes. However, our genetic segregation data was consistent in rejecting this hypothesis (Hauck et al., 2006a, b). Therefore, it is postulated that these cultivars contained a fourth allele and termed it Snull as it was not possible to resolve this fourth allele on a Southern blot with either S-RNase or SFB probes.

Boskovic et al. (2006) reported the presence of a different S-RNase, SD, in ‘Köröser’, ‘Montmorency’, and ‘Bruine Waalse’ using a different S-RNase consensus primer pair (EM-PC2consFD and EM-PC5consRD). Therefore, this alternate primer pair was used to test the possibility that SD might be the Snull allele. Using this primer pair, a ~ 530 bp amplification product was amplified in eight selections (‘Crisana’, ‘Erdi Botermo’, ‘Montmorency’, ‘Pandy 38’, ‘Pandy 114’, ‘Surefire’, ‘Tschernokorka’, and ‘Újfehértói furtos’; Fig. 1B), that had a similar size to the SD-RNase (Boskovic et al., 2006). The 527 bp S-RNase PCR products amplified with the EM-PC2consFD and EM-PC5consRD primers from ‘Crisana’, ‘Pandy 114’, and ‘Újfehértói furtos’ were cloned and sequenced and the nucleotide sequences were identical. Unfortunately, these nucleotide sequences could not be compared with that of SD-RNase since the latter is not available in GenBank. Nevertheless, a comparison using the amino acid sequence available in Fig. 2 of Boskovic et al. (2006) revealed that our sequence differed from SD-RNase in two regions (see Supplementary Fig. 4 at JXB online). One region is just before the first intron. The SD-RNase has two additional glycine residues not present in our S-RNase sequence. The other region is the RHV region just after the second intron. The reasons for these discrepancies are not known. However, as the validity of our full length S-RNase sequences were verified multiple times from fosmid clones from three different genotypes (see below), and genetically verified in inheritance studies (see below), it has been tentatively postulated that our sequence indeed represents the previously identified SD-haplotype. This S-haplotype is named S35.

To clone the S35-RNase and SFB35, a fosmid library of ‘Montmorency’ (S6S13mS35S36a) was first screened using a probe of DIG-dUTP-labelled with a 527 bp S35-RNase fragment amplified by EM-PC2consFD and EM-PC5consRD from ‘Crisana’. Sixteen positive clones were obtained. Among them, two clones (Mon37 and Mon46) were shown to contain the S35-RNase since an ~530 bp fragment was amplified by PCR with the S-RNase consensus primer pair (EM-PC2consFD and EM-PC5consRD) but no fragment was obtained with the other S-RNase consensus primer pair (Pru-C2 and PCE-R). Unfortunately these two clones did not contain the SFB35 since the SFB band was not amplified by PCR with the SFB consensus primer pair (SFB-C1F and SFB-C2R). Next the ‘Montmorency’ fosmid library was screened again with a mixture of DIG-dUTP-labelled SFB6 and SFB36a probes at lower stringency (55 °C) and obtained 25 positive clones. Among them, six clones (Mon64, Mon71, Mon74, Mon76, Mon122, and Mon132) were considered to contain the novel SFB since these six clones amplified a SFB fragment by PCR with the SFB consensus primer (SFB-C1F and SFB-C2R) but not with allele specific primer pairs for SFB6, SFB13, and SFB36a/SFB36b. This novel SFB was named SFB35. These six clones were also analysed by PCR with the S-RNase primer pair EM-PC2consFD and EM-PC5consRD, but no fragment was amplified from any of the six clones. Therefore these six clones were considered to contain SFB35 but not to contain S35-RNase. This suggests that the intergenic distance between the S35-RNase and SFB35 is larger than that observed in the majority of cherry S-haplotypes (380 bp to 38 kb) identified to date (Ikeda et al., 2005) as the average insert size of our fosmid clones was ~ 40 kb. Mon37 and Mon64 were sequenced to determine the complete nucleotide sequences of the S35-RNase (GenBank accession no. EU054327 [GenBank] ) and SFB35 (GenBank accession no. EU054330 [GenBank] ), respectively. The S35-RNase, which consists of 232 amino acid residues (Fig. 2), is very different from other S-RNases identified in cherry. The second intron of S35-RNase is extremely short, consisting of only 82 bp (see Supplementary Fig. 4 at JXB online). The S-RNase consensus primer pair (Pru-C2 and PCE-R) could not amplify S35-RNase because the Pru-C2 primer was designed based on the amino acid sequence LWPSNYSN and the S35-RNase has LWPSNYSK. The PCE-R primer was designed based on the amino acid sequence EXEWNK, but the deduced amino acid sequence for S35-RNase is GREWNK.

SFB35 is composed of 371 amino acid residues that exhibited the characteristic variability patterns of previously identified SFBs (Fig. 3). Using these sequences, S35-RNase and SFB35 specific primers were designed (Table 1). The S35-RNase specific primer pair amplified a fragment of the expected size (435 bp) in all eight selections from which the 527 bp PCR product was amplified with the EM-PC2consFD and EM-PC5consRD primer pair (Fig. 4E). The SFB35 specific primer pair also amplified fragments of the expected size (557 bp) in all eight selections that had a 527 bp fragment following amplification using the EM-PC2consFD and EM-PC5consRD primer pair (Fig. 4F). It was further confirmed that the ~560 bp fragment was amplified by the SFB35 specific primer pair in all six fosmid clones (Mon64, Mon71, Mon74, Mon76, Mon122, and Mon132) (data not shown).

S-haplotype functionality
The functionality of the S-haplotypes identified in sour cherry was tested using populations derived from self- and cross-pollination. In the parents of these populations the S-haplotype being tested would segregate in both the eggs and pollen grains in an expected 1:1 ratio assuming that these gametes are viable and equally probable of occurring in a successful gamete. Therefore, if the S-haplotype is functional in the style and pollen, the pollen that contains that S-haplotype would be incompatible and the ratio of that S-haplotype would be 1:1 to represent its expected ratio in the eggs. However, if the S-haplotype being tested has lost either pollen or stylar function, the pollen carrying that S-haplotype would be compatible and that S-haplotype would be expected to be present in the progeny in a 3:1 ratio representing the additional contribution of the pollen S-haplotype. Using this strategy, multiple functional and non-functional S-haplotypes were previously identified in sour cherry (Hauck et al., 2006b).

To determine if the S33-haplotype in ‘Meteor’ was fully functional, the self-pollinated progeny of ‘Meteor’ were genotyped for their S-haplotypes. The progeny segregation ratio for S33 fit the expected 1:1 ratio and rejected the 3:1 ratio indicating that S33 is a functional S-haplotype (Table 3). Self-pollinated progeny of ‘Tamaris’ also segregated according to a 1:1 rejecting the 3:1 ratio for the S34-haplotype indicating that this S-haplotype is functional (Table 3). As ‘Tamaris’ is the first report of the presence of the sweet cherry S16-haplotype in sour cherry, the functionality of this S-haplotype was also tested. S16 segregated according to a 1:1 ratio indicating that the S16-haplotype is also functional.


View this table:
[in this window]
[in a new window]

 
Table 3. Progeny segregation of the S16-, S33-, S34- and S35-haplotypes to test the functionality of each S-haplotype

 
S35 was identified in several cultivars with diverse geographic origin; therefore, three different cultivars with diverse origins were used to test the functionality of S35. Segregation of S35 in the self-pollinated progeny of ‘Montmorency’ fit a 1:1 ratio and rejected the 3:1 ratio indicating that S35 is a functional S-haplotype in this cultivar (Table 3). As S35 is also present in ‘Újfehértói furtos’ and ‘Surefire’ it allowed us to test whether the S35 containing pollen from ‘Surefire’ is compatible in ‘Újfehértói furtos’ styles. The segregation of S35 in the progeny from this cross segregated according to the 1:1 ratio, rejecting the 3:1 ratio confirming that S35 is a functional S-haplotype.

‘Montmorency’, ‘Újfehértói furtos’, and ‘Surefire’ had previously been determined to have the Snull-haplotype as our genetic data predicted the presence of a fourth S-haplotype (Hauck et al., 2006b). To test the hypothesis that Snull is the newly identified S35, the progeny from self-pollinated ‘Montmorency’ and ‘Újfehértói furtos’x‘Surefire’ were S-genotyped with the addition of the S35. In all cases S35 co-segregated with the previous prediction of Snull (data not presented) confirming that S35 is the fourth S-haplotype in these selections.

Phylogenetic analyses
Phylogenetic analyses of amino acid and nucleotide sequences of S-RNase and SFB sequences from eight species of Prunus (Fig. 7) produced trees in which the closest relatives of many alleles were alleles from other species. This pattern, first reported for Solanaceae species, was named trans-specific evolution by Richman et al. (1996). In Prunus this pattern has been described before (see references in Vieira et al., 2008), but in contrast with the observation made for the Solanaceae species, in this genus trans-specific evolution cannot be taken as evidence for the very old age of alleles (Vieira et al., 2008). Although this pattern has been known for more than a decade, no proper phylogenetic analyses have so far been performed to show conclusively that this pattern is not simply due to lack of phylogenetic resolution. In Prunus, bootstrap support values are generally weak. Nevertheless, in no case were all the alleles of any one gene from any one species supported as monophyletic.

Phylogenetic analyses of nucleotide sequences of both S-RNase and SFB data resolved the P. cerasus S35 along with the P. dulcis Sa as highly divergent from the remaining S-haplotypes (Fig. 7). Therefore, S35 was identified as not only divergent from other cherry S-haplotypes but also divergent from other Prunus S-haplotypes. For the P. cerasus S33 and P. cerasus S34, as well as for the alleles from several other S-haplotypes, the resolved relationships were different in the two data sets (Fig. 7). All data sets strongly supported the sister relationship between P. cerasus S34 and P. domestica S5. Both the amino acid and the nucleotide sequence data resolved the S-RNases of the last two haplotypes as sister to a clade including the P. avium S1-RNase, the P. dulcis S11-RNase, and the P. tenella S8-RNase. The SFB data did not, however, show this relationship, and the position of P. cerasus SFB34 plus P. domestica SFB5 was very weakly supported. The P. cerasus S33-RNase was weakly supported as sister to P. mume S1-RNase and P. spinosa S10-RNase; however, P. cerasus SFB33 was resolved, again with very weak support, as sister to the clade of P. avium SFB3, P. cerasus SFB34, and P. domestica SFB5. The relationships resolved by the amino acid sequence data for both genes (not shown) differed in some details from those resolved by the nucleotide sequence data, especially for the more weakly supported relationships.

Results of the Partition Homogeneity test revealed strong differences in phylogenetic signal between the S-RNase and the SFB nucleotide data sets as well as between the amino acid data sets for the two genes (P=0.001 for all tests). Tests were carried out for significant differences in the tree topologies between the S-RNase and SFB data for both nucleotides and amino acids using the Templeton, K–H, and winning-sites tests, as described above. For all cases in which data from one gene were constrained to give topologies identical to those found to be optimal for the other gene, the resulting trees were significantly longer than unconstrained trees (P <0.0001 for all tests). In contrast, when data from one gene were constrained to resolve only the groups supported with bootstrap values of 80% or greater by the other gene, the results were not significant for any data set or any test.

Because resolutions of relationships among the Prunus S-RNase and SFB alleles were generally weak, the possibility was tested that topologies consistent with one or more constraints concordant with our understanding of organismal relationships based on other data (Bortiri et al., 2001) were not significantly worse than the best (most parsimonious) topologies resolved for the two data sets containing just the Prunus sequences. Constraint trees with three distinct topologies were used to test these hypotheses. In the first topology, the ‘species’ constraint, all representatives of each species were constrained to form a monophyletic group, but no restrictions were placed on the relationships among species. In the second topology, the ‘subgenera’ constraint, each of several groups of 2–4 species (P. avium/P. cerasus, P. armeniaca/P. mume, P. dulcis/P. tenella, and P. cerasifera/P. domestica/P. salicina/P. spinosa) was constrained as monophyletic, but no other restrictions were placed on the relationships within or among species. These groups have been resolved as closely related in recent studies and they also correspond to previously recognized subgenera (Rehder, 1940). In the third topology, the ‘largeclades’ constraint, each of two major lineages resolved in recent studies, one including the two cherry species (P. avium and P. cerasus) and one including the rest, was constrained as monophyletic, but no other restrictions were imposed. For both S-RNase and SFB alleles, the shortest unconstrained trees were significantly better (P <0.0001 in all cases) than any tree consistent with any of these constraints, as assessed by the K–H, Templeton, and winning-sites tests.

Variability analyses
When comparing the S-RNase and SFB genes from P. tenella S8- and P. avium S1-haplotypes, Surbanovski et al. (2007), observed no amino acid differences between the S-RNases but 12 amino acid differences between the SFBs. In order to address the generality of this observation, the per site non-synonymous (Ka) and synonymous (Ks) values for pairs of closely related S-haplotypes (those showing a Ks value smaller than 0.1 for either the S-RNase or the SFB) was determined. There are 20 such S-haplotype pairs (Table 4). There is no tendency for an increased rate of non-synonymous or synonymous mutation at one of the genes compared with the other (Non-parametric Sign test, P >0.05 in both cases).


View this table:
[in this window]
[in a new window]

 
Table 4. Synonymous (Ks) and non-synonymous (Ka) per site divergence rates for the S-RNase and SFB genes from closely related S-haplotype pairs

 
Since there is no evidence for an increased rate of mutation at the SFB gene when compared with the S-RNase gene, it is conceivable that a rare recombination event could be responsible for the discrepancy observed by Surbanovski et al. (2007). Using 22 Prunus S-haplotypes, Nunes et al. (2006) noticed that the history of the two genes is positively correlated, although not highly correlated. When using all pair-wise Ks values, the Pearson correlation coefficient was +0.628 (n=231; P <0.01), and when using all pair-wise Ka values, the Pearson correlation coefficient was +0.759 (n=231; P <0.01). This observation is suggestive of recombination at one or both genes. Using the same methodology but a larger data set (n=48), the correlation coefficient was now even lower. When using all pair-wise Ks values, the Pearson correlation coefficient was +0.346 (n=1128; P <0.01) and when using all pair-wise Ka values, the Pearson correlation coefficient was +0.673 (n=1128; P <0.01). Therefore recombination cannot be ruled out as one conceivable reason for the discrepancy observed by Surbanovski et al. (2007). Nevertheless, it should be noted that for both the S-RNase and SFB gene there is no evidence for the clustering of synonymous or non-synonymous mutations, as expected under the hypothesis of a recent recombination event affecting only a portion of the SFB or S-RNase gene.


    Discussion
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Supplementary data
 References
 
Sour cherry is a segmental allotetraploid (Beaver and Iezzoni, 1993) arising from the inter-mating of the diploid sweet cherry and tetraploid ground cherry. Previous studies have shown that six S-haplotypes present in sweet cherry (S1, S4, S6, S9, S12, and S13) are present in sour cherry as well. In this study, it is shown that the sweet cherry S14- and S16-haplotypes are also found in sour cherry, bringing to eight the total number of S-haplotypes shared by the two species. The identification in this study of three S-haplotypes not previously identified in sour cherry and not found in sweet cherry suggests that these S-haplotypes were contributed by the ground cherry (P. fruticosa) parent. To address this question, an ongoing survey of the P. fruticosa present in the Michigan State University cherry germplasm collection has resulted in the identification of multiple accessions possessing the S33- and S35-haplotypes. However, no individual possessing the S34-haplotype has been identified (T Tsukamoto et al., unpublished data).

The finding that the S35-haplotype is functional has significant implications regarding the genetic control of SI and SC in sour cherry. The ‘one allele match model’ states that for any sour cherry cultivar to be SI it must possess a minimum of three functional S-haplotypes (Hauck et al., 2006b). This is because all pollen carrying at least one matching functional S-haplotype is predicted to be incompatible. The landrace cultivars ‘Pandy 38’, ‘Pandy 114’ (syn. ‘Crisana’, ‘Köröser’), and ‘Tschernokorka’ are well documented examples of SI sour cherry cultivars (Redalen, 1984; Lansari and Iezzoni, 1990). Both these cultivars possess three functional S-haplotypes, thereby supporting the current hypothesis for the genetic control of SI and SC in sour cherry (Yamane et al., 2001; Hauck et al., 2006b). The three functional S-haplotypes in ‘Pandy 38’ and ‘Pandy 114’ are S1, S4, and S35 while the three functional S-haplotypes in ‘Tschernokorka’ are S9, S13, and S35. This also provides additional evidence that competitive interaction, for example, the compatibility of heteroallelic pollen that is associated with the breakdown of GSI in the Solanaceae, does not occur in tetraploid sour cherry as pollen containing two functional S-haplotypes is incompatible. Our result is in contrast to Huang et al. (2008) who reported a set of observations compatible with the ‘competitive interaction model’ in a SC selection of tetraploid Prunus pseudocerasus. However, as not all predictions of the model were tested, further data are needed to understand the generality of the ‘one allele match model’ in Prunus.

It is possible that the S33-RNase identified in ‘Meteor’ has the same specificity as the P. mume S10-RNase (subgenus Prunus). Confirmation will require the complete sequence of the P. mume S10-RNase. However, if these genes do indeed code for identical proteins, this would represent another example of ancestral specificities shared between Prunus subgenera. These findings are consistent with S-haplotype divergence predating speciation in Prunus (Nunes et al., 2006; Vieira et al., 2008). Prunus species from the same subgenus share, on average, a higher percentage of ancestral specificities than Prunus species from different subgenera (Vieira et al., 2008). Nevertheless, identical amino acid sequences, as in this case (the comparison of the S-RNase from P. cerasus S33 and P. mume S10), were only found between P. tenella S8-RNase and P. avium S1-RNase (Surbanovski et al., 2007), and P. tenella S3-RNase and P. salicina Sf-RNase (although the latter two sequences are partial; Vieira et al., 2008). The SFB of P. tenella S8- and P. avium S1-haplotypes differ, however, by 12 amino acids scattered along the gene (Surbanovski et al., 2007). Such comparison is not possible for P. cerasus S33- and P. mume S10-SFB alleles since the latter sequence is not available. Moreover, the SFB sequences of the P. tenella S3- and P. salicina Sf-haplotypes are not available.

The inferred rate of synonymous and non-synonymous mutation for closely related haplotypes is not higher at the SFB than at the S-RNase gene. Therefore, recombination may be the cause of the pattern observed by Surbanovski et al. (2007), although we failed to identify where the putative recombination event occurred.

Phylogenetic analyses provide strong support for the placement of the P. cerasus S35 and P. dulcis Sa-RNases and SFBs as a sister group to all the other alleles examined. These specificities are among the oldest ones, and are 15–20 million years old (Vieira et al., 2008). The significant conflict between the S-RNase and SFB data sets, as revealed by the partition homogeneity and tree topology tests, suggest real differences in the underlying phylogenetic histories of the S-RNase and SFB genes (Fig. 7). These results echo those of Nunes et al. (2006). Two observations suggest that homoplasy within one or both data sets, which could result from rapid evolution and/or intragenic recombination involving either or both genes, may be the primary cause of the lack of phylogenetic congruence. First, in phylogenies from both genes, there is generally weak resolution of relationships among alleles, especially for the deeper branches in the trees (Fig. 7). Second, those relationships that are strongly supported tend either to be strongly supported by sequences from both genes (e.g. the sister relationship of P. cerasus S34 and P. domestica S5) or they are strongly supported by one data set and weakly resolved in the other (e.g. the position of P. spinosa S7-1), as reflected by the non-significant results of the topology tests in which each data set was constrained to resolve only groups with 80% or better bootstrap support by the other gene. The fact that the relationships that are strongly supported by either data set, and especially those supported by both data sets, appear to involve recent divergences, while the deeper branches are generally weakly supported, suggests that whatever factors are causing phylogenetic incongruence between the two data sets are also acting to limit our ability to reliably construct more distant relationships among alleles of either gene alone. The lack of congruent topologies of relationship for S-RNase and SFB alleles does not, in any case, preclude the possibility that the two genes have coevolved to maintain self-incompatibility. Indeed, as expected, since amino acid sites responsible for specificity determination at the S-RNase and SFB must coevolve, the much higher correlation coefficient for non-synonymous divergence (+0.673) than for synonymous divergence (+0.346) suggests otherwise. As long as the history of amino acid sites important for specificity determination at the S-RNase and SFB genes is tightly linked, the S-haplotype will be functional. For instance, historical intragenic gene conversion may have preferentially occurred in regions where there are no amino acid sites important for specificity determination. It remains, nevertheless, to be determined how much intragenic recombination needs to be argued in order to account for the observed synonymous and non-synonymous correlation coefficients as well as for the overall topology incongruencies. It should be noted that features compatible with recombination have been previously reported at the S-RNase and SFB genes, although the evidence is still not unequivocal (Vieira et al., 2003; Nunes et al., 2006; Ortega et al., 2006).

The overall amounts of total change and the relative distribution of that change on internal and terminal branches are similar. In both cases, internal branches are generally short compared to many of the terminal branches, suggesting that cladogenic diversification in each gene may have occurred in a common ancestral species, followed by lineage sorting and recombination (not necessarily in that order). In agreement with this view, our phylogenetic analyses of S-RNase and SFB sequences produced gene trees in which alleles from a single species were not resolved as monophyletic. Such results are expected for genes involved in self-incompatibility among closely related species (Lu, 2001; Vieira et al., 2008). Given sufficient divergence time, however, the genes within a particular evolutionary lineage should eventually coalesce. It is known, for example, that within Rosaceae, the S-RNase genes of tribe Pyreae (Potter et al., 2007) form a clade distinct from those of Prunus (Igic and Kohn, 2001). In an attempt to ascertain the phylogenetic magnitude of the lack of correspondence between our gene trees and the presumed phylogeny for these species, a constraint tree approach was used (see Results). Constraining the topologies of the trees in any way to correspond to current understanding of organismal relationships resulted in significantly longer tree lengths. This was true at the level of species, groups of 2–4 closely related species, and two major clades that have been resolved in recent phylogenetic studies of the genus. Thus, our sampling of five of the roughly 200 species of Prunus, which includes representatives of two of the three major lineages within the genus resolved in recent phylogenetic analyses, indicate that neither the SFB nor the S-RNase locus has attained coalescence at any level below that of the entire genus. Sampling of additional species of Prunus, and of species of related genera in Rosaceae, will be required to test this hypothesis and to determine the precise level at which coalescence has occurred.


    Supplementary data
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Supplementary data
 References
 
The following supplementary data for this article are available at JXB online.

Fig. S1. Nucleotide sequence and the deduced amino acid sequence alignments for the Prunus cerasus S33-RNase and Prunus mume S10-RNase.

Fig. S2. Nucleotide sequence alignment for the Prunus cerasus S34-RNase, P. avium S1-RNase, P. dulcis S11-RNase, and P. tenella S8-RNase.

Fig. S3. Comparison between downstream sequence of Prunus cerasus S34-RNase and intergenic sequence between S-RNase and SFB of P. avium S1, P. dulcis S11, and P. tenella S8.

Fig. S4. Nucleotide sequence and the deduced amino acid sequence alignment for the Prunus cerasus S35-RNase and SD-RNase.


    Acknowledgements
 
This work was supported by a grant from the USDA Cooperative State Research, Education and Extension Service – National Research Initiative – Plant Genome Program Grant no. 2004-01543.


    Footnotes
 
* Present address: Department of Plant Sciences, University of Arizona, Tucson, AZ85721, USA. Back


    Abbreviations
 
GSI, gametophytic self-incompatibility; SI, self-incompatible; SC, self-compatible; SFB, S haplotype-specific F-box protein.


    References
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Supplementary data
 References
 
Anderson MA, Cornish EC, Wau S-L, et al. Cloning of cDNA for a stylar glycoprotein associated with expression of self-incompatibility in Nicotiana alata. Nature (1986) 321:38–44.[CrossRef][Web of Science]

Beaver J, Iezzoni AF. Allozyme inheritance in tetraploid sour cherry (Prunus cerasus L.). Journal of the American Society for Horticultural Science (1993) 118:873–877.[Abstract/Free Full Text]

Bortiri E, Oh S, Jiang J, Baggett S, Granger A, Weeks C, Buckingham M, Potter D, Parfitt D. Phylogeny and systematics of Prunus (Rosaceae) as determined by sequence analysis of ITS and the chloroplast trnL-trnF spacer DNA. Systematic Botany (2001) 26:797–807.[Web of Science]

Boskovic RI, Wolfram B, Tobutt KR, Cerovic R, Sonneveld T. Inheritance and interaction of incompatibility alleles in the tetraploid sour cherry. Theoretical and Applied Genetics (2006) 112:315–326.[CrossRef][Web of Science][Medline]

de Nettancourt D. Incompatibility and incongruity in wild and cultivated plants (2001) Berlin: Springer.

Entani T, Iwano M, Shiba H, Che F-S, Isogai A, Takayama S. Comparative analysis of the self-incompatibility (S-) locus region of Prunus mume: identification of a pollen-expressed F-box gene with allelic diversity. Genes to Cells (2003) 8:203–213.[Abstract]

Hauck NR, Yamane H, Tao R, Iezzoni AF. Self-compatibility and incompatibility in tetraploid sour cherry (Prunus cerasus L.). Sexual Plant Reproduction (2002) 15:39–46.[CrossRef][Web of Science]

Hauck NR, Ikeda K, Tao R, Iezzoni AF. The mutated S1-haplotype in sour cherry has an altered S-haplotype specific F-box protein gene. Journal of Heredity (2006a) 97:514–520.[Abstract/Free Full Text]

Hauck NR, Yamane H, Tao R, Iezzoni AF. Accumulation of non-functional S-haplotypes results in the breakdown of gametophytic self-incompatibility in tetraploid Prunus. Genetics (2006b) 172:1191–1198.

Huang S-X, Wu H-Q, Li Y-R, Wu J, Zhang S-J, Heng W, Zhang S-L. Competitive interaction between two functional S-haplotypes confer self-compatibility on tetraploid Chinese cherry (Prunus pseudocerasus Lindl.CV. Nanjing Chuisi). Plant Cell Reports (2008) 27:1075–1085.[CrossRef][Web of Science][Medline]

Iezzoni AF. Acquiring cherry germplasm from Central and Eastern Europe. HortScience (2005) 40:304–308.[Abstract/Free Full Text]

Igic B, Kohn JR. Evolutionary relationships among self-incompatibility RNases. Proceedings of the National Academy of Sciences, USA (2001) 98:13167–13171.[Abstract/Free Full Text]

Ikeda K, Igic B, Ushijima K, Yamane H, Hauck NR, Nakano R, Sassa H, Iezzoni AF, Kohn JR, Tao R. Primary structure futures of the S haplotype-specific F-box protein, SFB, in Prunus. Sexual Plant Reproduction (2004a) 16:235–243.[CrossRef][Web of Science]

Ikeda K, Ushijima K, Yamane H, Tao R, Hauck NR, Sebolt AM, Iezzoni AF. Linkage and physical distances between the S-haplotype S-RNase and SFB genes in sweet cherry. Sexual Plant Reproduction (2005) 17:289–296.[CrossRef][Web of Science]

Ikeda K, Watari A, Ushijima K, Yamane H, Hauck NR, Iezzoni AF, Tao R. Molecular markers for the self-compatible S4'-haplotype, a pollen-part mutant in sweet cherry (Prunus avium L.). Journal of the American Society for Horticultural Science (2004b) 129:724–728.[Abstract/Free Full Text]

Lai Z, Ma W, Han B, Liang L, Zhang Y, Hong G, Xue Y. An F-box gene linked to the self-incompatibility (S) locus of Antirrhinum is expressed specifically in pollen and tapetum. Plant Molecular Biolology (2002) 50:29–42.[CrossRef]

Lansari A, Iezzoni A. A preliminary analysis of self-incompatibility in sour cherry. HortScience (1990) 25:1636–1638.[Abstract/Free Full Text]

Lu Y. Roles of lineage sorting and phylogenetic relationship in the genetic diversity at the self-incompatibility locus of Solanaceae. Heredity (2001) 86:195–205.[Web of Science][Medline]

McClure BA, Haring V, Ebert PR, Anderson MA, Simpson RJ, Sakiyama F, Clarke AE. Style self-incompatibility gene products of Nicotiana alata are ribonucleases. Nature (1989) 342:955–957.[CrossRef][Web of Science][Medline]

McClure BA, Gray JE, Anderson MA, Clarke AE. Self-incompatibility in Nicotiana alata involves degradation of pollen rRNA. Nature (1990) 347:757–760.[CrossRef][Web of Science]

Nunes MDS, Santos RAM, Ferreira SM, Vieira J, Vieira C. Variability patterns and positively selected sites at the gametophytic self-incompatibility pollen SFB gene in a wild self-incompatible Prunus spinosa (Rosaceae) population. New Phytologist (2006) 172:577–587.[CrossRef][Web of Science][Medline]

Olden EJ, Nybom N. On the origin of Prunus cerasus L. Hereditas (1968) 59:327–345.[Web of Science]

Ortega E, Boskovic RI, Sargent DJ, Tobutt KR. Analysis of S-RNase alleles of almond (Prunus dulcis): characterization of new sequences, resolution of synonyms and evidence of intragenic recombination. Molecular and Genetic Genomics (2006) 276:413–426.[CrossRef]

Potter D, Eriksson T, Evans RC, et al. Phylogeny and classification of Rosaceae. Plant Systematics and Evolution (2007) 266:5–43.[CrossRef][Web of Science]

Redalen G. Fertility in sour cherries. Gartenbauwissenschaft (1984) 49:212–217.[Web of Science]

Richman AD, Uyenoyama MK, Kohn JR. Allelic diversity and gene genealogy at the self-incompatibility locus in the Solanaceae. Science (1996) 273:1212–1216.[Abstract]

Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics (2003) 19:2496–2497.[Abstract/Free Full Text]

Sassa H, Hiranao H, Ikehashi H. Self-incompatibility related RNases in styles of Japanese pear (Pyrus serotina Rehd.). Plant and Cell Physiology (1992) 33:811–814.[Abstract/Free Full Text]

Sijacic P, Wang X, Skirpan AL, Wang Y, Dowd PE, McCubbin AG, Huang S, Kao T-h. Identification of the pollen determinant of S-RNase-mediated self-incompatibility. Nature (2004) 429:302–305.[CrossRef][Web of Science][Medline]

Sonneveld T, Robbins TP, Boskovic R, Tobutt KR. Cloning of six cherry self-incompatibility alleles and development of allele specific PCR detection. Theoretical and Applied Genetics (2001) 102:1046–1055.[CrossRef][Web of Science]

Sonneveld T, Tobutt KR, Robbins TP. Allele-specific PCR detection of sweet cherry self-incompatibility (S) alleles S1 to S16 using consensus and allele-specific primers. Theoretical and Applied Genetics (2003) 107:1059–1070.[CrossRef][Web of Science][Medline]

Sutherland BG, Robbins TP, Tobutt KR. Primers amplifying a range of Prunus S-alleles. Plant Breeding (2004) 123:582–584.[CrossRef][Web of Science]

Surbanovski N, Tobutt KR, Konstantinovic M, Maksimovic V, Sargent DJ, Stevanovic V, Ortega I, Boskovic RI. Self-incompatibility of Prunus tenella and evidence that reproductively isolated species of Prunus have different SFB alleles coupled with an identical S-RNase allele. The Plant Journal (2007) 50:723–734.[CrossRef][Web of Science][Medline]

Swofford DL. PAUP*: phylogenetic analysis using parsimony (and other methods), version 4.0b10 (2002) Sunderland: Sinauer.

Tao R, Yamane H, Sugiura A, Murayama H, Sassa H, Mori H. 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 (1999) 124:224–233.[Abstract/Free Full Text]

Thompson J, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The ClustalX window interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research (1997) 25:4876–4882.[Abstract/Free Full Text]

Tobutt KR, Boskovic R, Cerovic R, Sonneveld T, Ruzic D. Identification of incompatibility alleles in the tetraploid species sour cherry. Theoretical and Applied Genetics (2004) 108:775–785.[CrossRef][Web of Science][Medline]

Tsukamoto T, Hauck NR, Tao R, Jiang N, Iezzoni AF. Molecular characterization of three non-functional S-haplotypes in sour cherry (Prunus cerasus). Plant Molecular Biology (2006) 62:371–383.[CrossRef][Web of Science][Medline]

Tsukamoto T, Tao R, Iezzoni AF. PCR markers for mutated S-haplotypes enable discrimination between self-incompatible and self-compatible sour cherry selections. Molecular Breeding (2008) 21:67–80.[CrossRef][Web of Science]

Ushijima K, Sassa H, Dandekar AM, Gradziel TM, Tao R, Hirano H. Structural and transcriptional analysis of the self-incompatibility locus of almond: identification of a pollen-expressed F-box gene with haplotype-specific polymorphism. The Plant Cell (2003) 15:771–781.[Abstract/Free Full Text]

Ushijima K, Sassa H, Tao R, Yamane H, Dandekar AM, Gradziel TM, Hirano H. Cloning and characterization of cDNAs encoding S-RNases in almond (Prunus dulcis): primary structure features and sequence diversity of the S-RNases in Rosaceae. Molecular and General Genetics (1998) 260:261–268.[CrossRef][Web of Science]

Ushijima K, Sassa H, Kusaba M, Tao R, Gradziel TM, Dandekar AM, Hirano H. Characterization of the S-locus region of almond (Prunus dulcis): analysis of a somaclonal mutant and a cosmid contig for an S-haplotype. Genetics (2001) 158:379–386.[Abstract/Free Full Text]

Vaughan SP, Boskovic RI, Gisbert-Climent Russell K, Tobutt KR. Characterization of novel S-alleles from cherry (Prunus avium L.). Tree Genetics and Genomes (2008) 4:531–541.[CrossRef]

Vieira CP, Charlesworth D, Vieira J. Evidence for rare recombination at the gametophytic self-incompatibility locus. Heredity (2003) 91:262–267.[CrossRef][Web of Science][Medline]

Vieira J, Fonseca NA, Santos RAM, Habu T, Tao R, Vieira CP. The number, age, sharing and relatedness of S-locus specificities in Prunus. Genetics Research (2008) 90:17–26.[Medline]

Wheeler D, Newbigin E. Expression of 10 S-class SLF-like genes in Nicotiana alata pollen and its implications for understanding the pollen factor of the S locus. Genetics (2007) 177:2171–2180.

Wünsch A, Hormaza JI. Cloning and characterization of genomic DNA sequences of four self-incompatibility alleles in sweet cherry (Prunus avium L.). Theoretical and Applied Genetics (2004) 108:299–305.[CrossRef][Web of Science][Medline]

Xia X, Xie Z. DAMBE: Software package for data analysis in molecular biology and evolution. Journal of Heredity (2001) 92:371–373.[Abstract/Free Full Text]

Xue Y, Carpenter R, Dickinson HG, Coen ES. Origin of allelic diversity in Antirrhinum S locus RNases. The Plant Cell (1996) 8:805–814.[Abstract]

Yamane H, Ikeda K, Hauck NR, Iezzoni AF, Tao R. Self-incompatibility (S) locus region of the mutated S6-haplotypes of sour cherry (Prunus cerasus) contains a functional pollen S-allele and a non-functional pistil S-allele. Journal of Experimental Bototany (2003a) 54:2431–2434.[CrossRef]

Yamane H, Ikeda K, Ushijima K, Sassa H, Tao R. A pollen-expressed gene for a novel protein with an F-box motif that is very tightly linked to a gene for S-RNase in two species off cherry, Prunus cerasus and P. avium. Plant and Cell Physiology (2003b) 44:764–769.[Abstract/Free Full Text]

Yamane H, Tao R, Sugiura A, Hauck NR, Iezzoni AF. Identification and characterization of S-RNases in tetraploid sour cherry (Prunus cerasus). Journal of the American Society for Horticultural Science (2001) 126:661–667.[Abstract/Free Full Text]

Yamane H, Ushijima K, Sassa H, Tao R. The use of S-haplotype-specific F-box protein gene, SFB, as a molecular marker for S-haplotypes and self-compatibility in Japanese apricot (Prunus mume). Theoretical and Applied Genetics (2003c) 107:1357–1361.[CrossRef][Web of Science][Medline]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Material
Right arrowOA All Versions of this Article:
59/11/3169    most recent
ern172v1
Right arrow E-letters: Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when E-letters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Tsukamoto, T.
Right arrow Articles by Iezzoni, A. F.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Tsukamoto, T.
Right arrow Articles by Iezzoni, A. F.
Agricola
Right arrow Articles by Tsukamoto, T.
Right arrow Articles by Iezzoni, A. F.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?