Self-(in)compatibility in apricot germplasm is controlled by two major loci, S and M.

BACKGROUND: Apricot (Prunus armeniaca L.) exhibits a gametophytic self-incompatibility (GSI) system and it is mostly considered as a self-incompatible species though numerous self-compatible exceptions occur. These are mainly linked to the mutated S C-haplotype carrying an insertion in the S-locus F-box gene that leads to a truncated protein. However, two S-locus unlinked pollen-part mutations (PPMs) termed m and m' have also been reported to confer self-compatibility (SC) in the apricot cultivars 'Canino' and 'Katy', respectively. This work was aimed to explore whether other additional mutations might explain SC in apricot as well.
RESULTS: A set of 67 cultivars/accessions with different geographic origins were analyzed by PCR-screening of the S- and M-loci genotypes, contrasting results with the available phenotype data. Up to 20 S-alleles, including 3 new ones, were detected and sequence analysis revealed interesting synonymies and homonymies in particular with S-alleles found in Chinese cultivars. Haplotype analysis performed by genotyping and determining linkage-phases of 7 SSR markers, showed that the m and m' PPMs are linked to the same m 0-haplotype. Results indicate that m 0-haplotype is tightly associated with SC in apricot germplasm being quite frequent in Europe and North-America. However, its prevalence is lower than that for S C in terms of frequency and geographic distribution. Structures of 34 additional M-haplotypes were inferred and analyzed to depict phylogenetic relationships and M 1-2 was found to be the closest haplotype to m 0. Genotyping results showed that four cultivars classified as self-compatible do not have neither the S C- nor the m 0-haplotype.
CONCLUSIONS: According to apricot germplasm S-genotyping, a loss of genetic diversity affecting the S-locus has been produced probably due to crop dissemination. Genotyping and phenotyping data support that self-(in)compatibility in apricot relies mainly on the S- but also on the M-locus. Regarding this latter, we have shown that the m 0-haplotype associated with SC is shared by 'Canino', 'Katy' and many other cultivars. Its origin is still unknown but phylogenetic analysis supports that m 0 arose later in time than S C from a widely distributed M-haplotype. Lastly, other mutants putatively carrying new mutations conferring SC have also been identified deserving future research.

Background

Gametophytic self-incompatibility (GSI) is a system widely distributed in the plant kingdom [1] that prevents self-fertilization favoring outcrossing [2]. GSI specific recognition is under the control of a multi-allelic locus, termed S-locus, containing at least two linked genes: a pistil expressed S-RNase [3] and the pollen expressed S-locus F-box [4-6]. S-locus F-box proteins are thought to be components of E3 ubiquitin ligase complexes that recognize non-self S-RNases promoting their ubiquitination and degradation by the 26S proteasome proteolytic pathway [7, 8]. Recently, the collaborative model proposed in Solanaceae (Petunia) suggests that several F-box proteins are necessary to recruit S-RNases for degradation [9, 10]. This system seems to be extended to other plant families exhibiting GSI such as Rosaceae and particularly the Maloideae subfamily [10, 11]. However, Prunus does not seem to follow this model since knock-out of the S-locus F-box (SFB) gene leads to self-compatibility (SC) in contrast with the observations in Solanaceae. Reasons behind this behavior have been speculated for a long time [12, 13] but only recently evidence have been provided supporting S-locus F-box like-2 protein as a ‘general inhibitor’ that detoxifies S-RNases unspecifically, unless affected by SFB [14].

In general, self-incompatibility (SI) trait predominates in stone fruits (Prunus genus) in accordance with the high degree of heterozygosity showed by many of these species. However, different ‘degrees’ of SC have been detected in this genus ranging from sweet cherry, almost strictly self-incompatible (with a few exceptions), or apricot, mainly self-incompatible (SC is restricted to European cultivars), to complete SC in peach [15]. SC is mostly related to mutations in the S-locus genes and mutations affecting both S-RNases and SFBs have been detected in many Prunus species [13, 16]. In apricot (Prunus armeniaca L.) the S C-allele known to confer SC has been well characterized showing that a 358 bp insertion in the SFB gene leads to a putative truncated protein lacking the two essential 3′-hypervariable domains HVa and HVb [17]. Furthermore, the origin and dissemination of S C has also been studied, identifying the non-mutated ancestor S 8-allele and detecting its presence in accessions from different geographic areas [18]. In fact, most S-genotyped self-compatible apricot cultivars have been shown to carry the S C-allele [18-22].

Along with the S-locus specific products, other S-locus unlinked factors are also necessary for the GSI system to work. These factors known as ‘modifiers’ were firstly identified in Solanaceae [23]. Nevertheless, genetic evidence supporting modifiers have also been accumulated in other species including Prunus spp. [17, 24]. In apricot, pollen-part mutations (PPMs) conferring SC by putatively affecting modifiers have been identified in the Spanish local cultivar ‘Canino’ and the North-American one ‘Katy’. Both PPMs were mapped at the distal end of chr.3 within the so-called M- and M’-loci, respectively [25, 26]. Given their distinct geographic origins and the significant genetic distance estimated between both cultivars, it was initially hypothesized that these two PPMs could have arisen independently but affecting the same gene [26].

In this context, we thought that additional S-locus unlinked mutations (and not just restricted to the M-locus) might also be present in apricot germplasm. Thus, in this study we have genotyped both the S- and M-loci in a wide set of apricot cultivars and accessions from different geographic origins, by using SFB and S-RNase introns and M-locus linked SSR markers. These data were combined with those of phenotyping to dissect SC causes and distribution in apricot germplasm.

Results

Self-incompatibility versus self-compatibility in apricot

A set of 67 apricot cultivars and accessions (hereinafter referred only as accessions) representing a wide range of geographic origins has been analyzed. Information about pedigree was only available for a few of them but it was quite useful to support genotyping data. Blooming time and self-(in)compatibility phenotypes were scored according to literature reports. Thus, apricot accessions were classified into 5 gross classes (early, mid-early, mid, mid-late and late) regarding the first phenotype and as self-incompatible or self-compatible regarding the second (Table 1).

Table 1

Country of origin, pedigree, self-(in)compatibility and blooming time phenotypes of the analyzed apricot accessions

	Cultivar	Sourceb	Geog. areac	Country of origind	Pedigreee	SI/SCe	Blooming timee
1	Alba	CAPA	WE	Italy	Unknown	?	Midf
2	ASPa	FM	WE	Spain (V)	Unknown (isolated tree) [25]	?	Midf
3	Aurora	CAPA	NA	USA	RR17–62 × NJA-13 [26]	SI [73]	Mid [74]
4	Bebecou	MAGRAMA	SE/NAf	Greece	Unknown [41]	SC [73]	Mid-early [41]
5	Bergeron	MAGRAMA	WE	France	Unknown [41]	SC [73]	Mid-late [41, 47]
6	Budapest	St. Istvan U	EE	Hungary	‘Nancy’ × (‘Acme’, ‘Hungarian Best’, ‘Kései Rózsa’) [16]	SC [16]	Mid-late [75]
7	Búlida	MAGRAMA	SE/NAf	Spain (M)	Unknown [41]	SC [73]	Mid [41]
8	Canino	IVIA	WE	Spain (V)	Unknown [41]	SC [73]	Mid [41]
9	Canino 9–7	IVIA	WE	Spain (V)	Clonal selection from Canino [27]	?	Mid [27]
10	Canino 14–4	IVIA	WE	Spain (V)	Clonal selection from Canino [27]	?	Mid [27]
11	Canino 14–6	IVIA	WE	Spain (V)	Clonal selection from Canino [27]	?	Mid [27]
12	Castlebrite	CEBAS	NA	USA	o.p. (Perfection × Castleton) [76]	SC [73]	Mid-latef
13	Castleton	CEBAS	NA	USA	Perfection × Newcastle [29]	SC [76]	Mid [41]
14	Ceglédi óriás	St. Istvan U	EE	Hungary	Unknown (local selection) [16]	SI [32]	Mid [35]
15	Colorao	CEBAS	SE/NAf	Spain (M)	Unknown [41]	↗ sterile [73]	Mid [41]
16	Corbató	MAGRAMA	WE	Spain (V)	Unknown [25]	?	Mid-late [74]
17	Cow-1	CAPA	WE	France	INRA	?	Mid [35]
18	Cow-2	CAPA	WE	France	INRA	?	?
19	Cristalí	FM	WE	Spain (V)	Unknown [25]	?	Mid [74]
20	Currot	IVIA	WE	Spain (V)	Unknown [25]	SC [73]	Early [74]
21	Dulcinea	CAPA	WE	Italy	Unknown (Toscana variety)	SC [73]	Mid [35]
22	Effect	St. Istvan U	EE	Ukraine	Krupnolodnyi o.p. [16]	SC [16]	Late [77]
23	Ezzine	CAPA	SE/NAf	Tunisia	INRAT	?	Earlyf
24	Fergani	St. Istvan U	EE	Former USSR	Unknown	?	?
25	Galta Rojaa	CAPA	WE	Spain (V)	Unknown [25]	SC [73]	Mid-early [41]
26	GVV	FM	WE	Spain (V)	Unknown (isolated tree) [25]	?	?
27	Gandía	FM	WE	Spain (V)	Unknown [25]	?	Mid-early [74]
28	Gavatxet	FM	WE	Spain (V)	Unknown [25]	?	Mid [74]
29	Ginesta	IVIA	WE	Spain (V)	Unknown [25]	SCe	Mid-early [74]
30	Goldrich	IVIA	NA	USA	Sunglo × Perfection [76]	SI [73]	Mid [41, 47]
31	Gönci Magyar	St. Istvan U	EE	Hungary	Clone/hybrid? of Hungarian Best [16]	SC [16]	Mid-late [75]
32	Harcota	IVIA	NA	Canada	[(Geneva × Naramata) × Morden 604] × NJA1(Phelps × Perfection) [76]	SI [73]	Mid [41]
33	Hargrand	St. Istvan U	NA	Canada	V51092 [(Reliable × o.p.) × o.p.] × NJA1 (Phelps × Perfection) [76]	SI [73]	Latef
34	Harlayne	IVIA	NA	Canada	V51092 [(Reliable × o.p.) × o.p.] × Sunglo [76]	SC [73]	Mid-late [41, 77]
35	Henderson	IVIA	NA	USA	Unknown [76]	SC?	Mid-late [41, 77]
36	Katy	IVIA	NA	USA	Zaiger’s genetics (USA) [35]	SC [35]	Earlyf
37	Kech-pshar	St. Istvan U	EE	Uzbekistan	Unknown (local selection) [16]	?	?
38	K. Pozdniia	St. Istvan U	EE	Ukraine	Unknown (chance seedling) [16]	SC [16]	Mid-late [75]
39	Lambertin-1	CEBAS	NA	USA	[Perfection × (Royal x Blush)] o.p. × (Perfection × Royal) [35]	SI [73]	Mid [35]
40	Lito	IVIA	SE/NAf	Greece	SEO × Tyrinthos [32]	SC [32]	Mid-latef
41	Manrí	FM	WE	Spain (V)	Unknown [25]	?	Mid-early [74]
42	Mari de Cenad	St. Istvan U	EE	Romania	Unknown [16]	SC [16]	?
43	Mariem	CAPA	SE/NAf	Tunisia	2nd generation (Bergeron × Ouardi) × (Carraut × Crossa-Raynaud) [35]	?	Mid [35]
44	Martinet	FM	WE	Spain (V)	Unknown [25]	?	Mid [74]
45	Mitger	IVIA	WE	Spain (V)	Unknown [25]	SCe	Mid [35]
46	Moniquí	CEBAS	SE/NAf	Spain (M)	Unknown [41]	SI [73]	Mid [41]
47	Ninfa	CAPA	WE	Italy	Ouardi × Tyrinthos [35]	SC [73]	Mid-early [47]
48	Orange Red	IVIA	NA	USA	Lasgerdi Mashhad × NJA2 [35]	SI [73]	Mid [41, 47]
49	Ouardia	IVIA	SE/NAf	Tunisia	Canino × Hamidi [41]	SI [73]	Mid-early [41]
50	Palabras	IVIA	WE	Spain (V)	Unknown [25]	SCe	Mid-early [74]
51	Palau	IVIA	WE	Spain (V)	Unknown [25]	SCe	Mid-early [74]
52	Patterson	CEBAS	NA	USA	F2 seedling Perfection × unknown [76]	SC [73]	Mid-late [76]
53	Perla	CAPA	WE	Italy	ANFIC (Italy)	?	Mid-early [35]
54	Porticia	FM	WE	Italy	Unknown [41]	SC [73]	Mid-early [41]
55	Rojo de Carlet	IVIA	WE	Spain (V)	Unknown [25]	SCe	Mid [74]
56	Rózsakajszi	St. Istvan U	EE	Hungary	Local selection Nagykórös [16]	SC [16]	Mid-late [75]
57	Sayeba	CEBAS	SE/NAf	Tunisia	Canino × Hamidi [41]	SC [73]	Mid-early [41]
58	SEOa	IVIA	NA	USA	Unknown [76]	SI [73]	Mid-late [41]
59	Shalaha	St. Istvan U	EE	Armenia	Unknown	SC/SI? [73]	Late [77]
60	Stella	IVIA	NA	USA	Unknown [41]	SI [73]	Late [41, 77]
61	S. Mammuta	St. Istvan U	EE	Hungary	Unknown (local selection)	SI [73]	Mid [75]
62	Tadeo	IVIA	WE	Spain (V)	Unknown [25]	SC [73]	Mid-late [74]
63	Tirynthos	IVIA	SE/NAf	Greece	Unknown [41]	SC [73]	Early [41]
64	Trevatt	CEBAS	NA	Australia	Unknown [35]	SC [32]	Mid-late [41]
65	Veecot	IVIA	NA	Canada	o.p. from Reliable [76]	SI [73]	Mid [41]
66	Velázquez	MAGRAMA	SE/NAf	Spain (M)	Unknown	SI [73]	Latef
67	Xirivello	CAPA	WE	Spain (V)	Unknown [25]	?	Mid-late [74]

aSynonyms and acronyms: ASP ‘Albaricoquero Sin Polen’; Galta Roja ‘Galta Roja de Mitger’; GVV ‘Galta Vermella Valenciana’; Ouardi ‘Priana’; Portici ‘Portici-6’; Sayeb ‘Beliana’; SEO ‘Stark Early Orange’ and Shalah ‘Erevani’; K. Pozdnii ‘Konservnyi Pozdnii’; S. Mammut ‘Szegedi Mammut’

bMaterial sources (see Methods)

cMain geographic areas defined in this study: Western Europe (WE); Southern Europe/North Africa (SE/NAf); Eastern Europe (EE) and North America (NA)

dAccessions from Spain come from two different regions: Valencia (V) and Murcia (M). The first ones were assigned to WE and the second ones to SE/NAf

eReferences for pedigree are indicated when available. Reports for SI/SC and blooming time phenotypes are also indicated marking with (?) when no consistent data were found: García et al. [41]; Halázs et al. [27]; Della Strada et al. [76]; Halázs et al. [18]; Orero et al. [35]; Brooks and Olmo [32]; Russell [73]; Syrgianidis and Mainou [74]; Burgos et al. [47]; Badenes et al. [75]; Massai [77]; Nyujtó et al. [78]; Mehlenbacher et al. [36]

fIVIA’s own data (see footnotes in Table 3)

Table 3

S- and M-haplotypes assigned to the apricot accessions analyzed in this study

	Cultivar	SI/SC	S-genotd	M-genote		Cultivar	SI/SC	S-genot	M-genot
1	Alba	SCb	S 1/S C	M 1–1/M 3	35	Henderson	SC?	S 3/S 17	M 2–0/M 16
2	ASP	↗sterileb,c	S 5/S C	M 4–0/M 5–1	36	Katya	SC/SCb	S 1/S 2	M 3/m 0–0
3	Auroraa	SI	S 1/S 17	M 3/M 11	37	Kech-pshara	?	S 15/S Z	M ?/M ?
4	Bebecou	SC	S 6/S C	M 4–0/M 4–2	38	Konservnyi P. a	SC	S 2/S C	M 4–1/M 8–1
5	Bergerona	SC	S 2/S C	M 4–0/M 8–0	39	Lambertin-1a	SI	S 1/S 2	M 1–0/M 17
6	Budapesta	SC	S 2/S C	M 1–1/M 12	40	Lito	SCb	S 6 / S C	M 6/M 10
7	Búlida	SC	S 5/S C	M 4–1/M 5–1	41	Manrí	?	S C / S C	m 0–0 / m 0–0
8	Caninoa	SC/SCb	S 2/S C	M 1–0/m 0–0	42	Mari de Cenada	SC	S C/S 19	M 8–0/M 18
9	Canino 9–7	?	S 2/S C	M 1–0/m 0–0	43	Mariem	SC?b	S 7/S 20	M 1–0/M 8–2
10	Canino 14–4	SCb	S 2/S C	M 1–0/m 0–0	44	Martinet	?	S 2/S 2	m 0–0 / m 0–1
11	Canino 14–6	SCb	S 2/S C	M 1–0/m 0–0	45	Mitger	SCb	S C / S C	M 5–0/M 5–0
12	Castlebrite	SC	S 2/S 2	M 3/m 0–0	46	Moniquía	SI	S 2/S 6	M 4–1/M 14–0
13	Castleton	SC	S 1/S 2	M 3/m 0–0	47	Ninfa	SC	S 7/S C	M 7–0/M 10
14	Ceglédi óriása	SI	S 8/S 9	M 12/M 13	48	Orange Red	SI	S 6/S 17	M 2–0/M 2–1
15	Coloraoa	↗sterile	S 5/S C	M 4–0/M 4–0	49	Ouardia	SI	S 2/S 7	M 1–0/M 7–0
16	Corbató	SCb	S 2/S 5	M 5–0/m 0–0	50	Palabras	SC/SCb	S C / S C	m 0–0 / m 0–0
17	Cow-1	SCb	S 1/S 31	M 3/m 0–0	51	Palaua	SC/SCb	S C / S C	m 0–0 / m 0–0
18	Cow-2	SIb	S 20/S 31	M 1–0/m 0–0	52	Patterson	SC	S 1/S C	M 1–3/M 3
19	Cristalí	SCb	S 20/S C	M 5–0/m 0–0	53	Perla	SIb	S 2/S 20	M 1–2/M 15–1
20	Currota	SC/SCb	S C / S C	m 0–0 / m 0–0	54	Portici	SC/SCb	S 2/S 20	M 1–4/m 0–0
21	Dulcinea	SC	S 2/S C	M 7–3/M 14–1	55	Rojo Carlet	SCb	S C / S C	M 4–0/M 5–0
22	Effecta	SC	S 8/S C	M 8–0/M 12	56	Rózsakajszia,f	SC	S 2/S C	M 1–0/M 12
23	Ezzine	SCb	S 24/S C	M 1–0/M 7–1	57	Sayeba	SC	S 7/S C	M 1–0/M 7–2
24	Fergani	?	S V/S X	M ?/M ?	58	Shalaha	SC	S 5/S 11	M 8–1/M 19
25	GaltaRoja	SC/SCb	S C / S C	M 4–0/M 5–2	59	SEO	SI	S 6/S 17	M 2–0/M 6
26	GVV	?	S 2/S C	M 4–0/M 5–2	60	Stella	SI	S 6/S 20	M 9/M 9
27	Gandía	?	S C / S C	m 0–0 / m 0–0	61	Szegedi M.	SI	S 8/S 9	M 12/M 13
28	Gavatxet	?	S 20/S C	m 0–0 / m 0–0	62	Tadeo	SC	S 20/S C	M 15–0/m 0–0
29	Ginestaa	SC/SCb	S C / S C	m 0–0 / m 0–0	63	Tirynthos	SC	S C / S C	M 10/M 10
30	Goldricha	SI/SIb	S 1/S 2	M 1–0/M 2–0	64	Trevatt	SC	S 2/S C	M 1–0/m 0–0
31	Gönci Magyara	SC	S 8/S C	M 8–0/M 12	65	Veecot	SI/SIb	S 2/S 20	M 2–0/M 3
32	Harcota	SI/SIb	S 1/S 4	M 1–0/M 2–2	66	Velázquez	SI	S 5/S 20	M 4–0/M 4–1
33	Hargranda	SI	S 1/S 2	M 1–0/M 2–0	67	Xirivello	?	S C / S C	M 4–0/M 4–1
34	Harlayne	SC	S 3/S 20	M 2–0/M 9

aCultivars previously S-genotyped by Halázs et al. [18]; Vilanova et al. [19]; Halázs et al. [20]; Zuriaga et al. [26]; Halázs et al. [27]; Burgos et al. [42]; Alburquerque et al. [79]

bOwn data on SI/SC phenotype obtained in this work (see Table 2). Additionally, SC had been observed in a set of accessions grown under insect-proof screen house at IVIA (‘Rojo de Carlet’, ‘Mitger’, ‘Palabras’, ‘Palau’, ‘Currot’, ‘Ginesta’, ‘Canino 14–4’ and ‘Canino 14–6’) showing moderate fruit-setting (not quantified) across several years

cMale-sterility in the ASP accession was indicated by shrunken pale anthers

d S-allele nomenclature is proposed according to Vilanova et al. [19]; Halázs et al. [27]; Zhang et al. [29]; Wu et al. [30]; Gu et al. [31] and Halázs et al. [28]. S-haplotype associated with SC (S C) is written in bold

e M-haplotypes were named with two digits. The first one corresponds to the M-haplotype ‘main class’ and the second to the subtype. M-haplotype variants associated with SC (m 0–0 and m 0–1) are written in bold. Haplotypes designated by M ? could not be defined

f S-genotype determined for Rózsakajszi (S 2 S C) was not in agreement with that previously reported by Halázs et al. [18] (S C S C). Reasons for this discrepancy are still unknown

Self-(in)compatibility phenotype data were completed with our own results when adult trees were available. In this regard, self-pollination tests were used to determine (8) or to confirm (13) self-(in)compatibility phenotypes but also to obtain progenies useful to search for S-locus unlinked mutations. Data suggest that 5 accessions are self-incompatible (‘Aurora’, ‘Cow-2’, ‘Perla’, ‘Veecot’ and ‘Velázquez’) while the remaining 16 show variable fruit-setting ranging from 0.5% (‘Búlida’) to 55% (‘Ninfa’), being recorded as self-compatibles (Table 2). Some progeny could be finally obtained from all accessions of this latter group except from ‘Búlida’. Those accessions with a progeny large enough were subsequently analyzed for non-S-locus mutations by S-genotyping embryos. The self-compatible controls (‘Canino’ and ‘Katy’), already studied, and the S C homozygote cultivar ‘Galta Roja’ were excluded. In summary, 40 out of the 67 accessions analyzed were classified as self-compatible, 16 as self-incompatible, 2 as male-sterile and the 9 remaining as undetermined (Table 3).

Table 2

Number of bagged flowers, fruit-setting percentage and SI/SC deduced phenotype in self-pollination tests

Cultivar/Accession	Bagged flowers	Fruit-Setting	% Setting	Phenotype	Progenyd
Albaa	355	37	10.4	SC	1
Aurora	350	0	0	SI	--
Bebecou	760	108	14.2	SC	96
Búlida	200	1	0.5	SC?b	--
Canino	412	99	24.0	SC	99
Castlebrite	300	36	12.0	SC	2
Corbatóa	320	52	16.3	SC	44
Cow-1a	200	8	4.0	SC	7
Cow-2a	315	0	0	SI	--
Cristalía	200	24	12.0	SC	13
Dulcinea	850	111	13.1	SC	104
Ezzinea	350	160	45.7	SC	21
Galta Roja	775	106	13.7	SC	106
Katy	731	80	10.9	SC	80
Mariema	450	11	2.4	SC?c	11
Ninfa	400	221	55.3	SC	12
Perlaa	370	0	0	SI	--
Portici	850	63	7.4	SC	59
Tadeo	375	14	3.7	SC	5
Veecot	450	0	0	SI	--
Velázquez	270	0	0	SI	--

aCultivars and accessions of unknown SI/SC phenotype

bIn general, fruit-setting percentages below 2–3% should not be undoubtedly associated to SC since some degree of pollen contamination cannot be fully discarded

c S-genotyping of the progeny suggested some degree of non-self pollination and therefore SC could not be confirmed

dNumber of individual plants obtained from the set fruits

S-genotyping: identification of new S-alleles

All the 67 accessions were S-genotyped using primer pairs amplifying different fragments of the S-haplotype region: the first and second S-RNase introns as well as the 5′-UTR SFB intron (Table 3). S-alleles were determined by comparing intron size patterns with previously S-genotyped accessions and, when needed, sequence analyses supported these assignments. Fragment analyses of the PCR products obtained with 5 primer combinations allowed us to detect up to 20 different S-alleles. Seventeen out of them had already been reported in apricot (S 1, S 2, S 3, S 4, S 5, S 6, S 7, S 8, S 9, S 11, S 15, S 17, S 18, S 19, S 20, S 24 and S C) and in the present work were designed basically according to the nomenclature established by Vilanova et al. [19], Halázs et al. [27] and Halázs [28], except for S 20, corresponding to the one described by Zhang et al. [29] and S 24 reported by Wu et al. [30] and Gu et al. [31]. Another S-allele (present in ‘Cow-1’ and ‘Cow-2’) was found to have a previously non-described intron size pattern being named as S 22, according to its homology with the NCBI GenBank accession HM053569, corresponding to the S 22-RNase from the Chinese apricot cultivar ‘Tuxiangbai’. Lastly, two additional S-alleles, preliminary suggested as S V and S X, were detected only once. The amplification of first S-RNase intron with the single primer pair SR1-F/SR1-R distinguished up to 13 S-alleles ranging from 260 (S 4) to 427 bp (S 17) (see Additional file 1: Table S1). Another primer combination (Pru-T2/SR1-R) allowed to amplify and distinguish two additional S-alleles S 6 and S 24. Exceptions were S 1 and S 7, and S C and S 8 pairs having exactly the same fragment sizes and S 19, S 20 and S X that could not be PCR-amplified with any of the four different primer combinations tested (data not shown). Sizes for the second S-RNase intron were approximately determined by agarose gel electrophoresis ranging from 300 (S 24) to 2800 bp (S C) and allowed to define 15 S-alleles with a single primer pair (Pru-C2/Pru-C6R). Exceptions were S 5 and S 6, and S 8 and S C pairs as well as S 19, S 20 and S V sharing similar fragment sizes, and S 3 that could not be amplified with any of the two primer combinations used. Variability in size of the 5′-UTR SFB intron containing fragments ranged from 189 (S 17) to 210 bp (S 1) and facilitated the identification of 9 S-alleles. Exceptions were S 2 and S 11, S 18 and S 22, and S C and S 8 pairs sharing the same intron sizes and S 4, S 5, S 9, S 19, S 20, S 24, S V and S X which fragments could not be amplified. Altogether, combined data allowed us to distinguish unambiguously 18 out of the 20 S-alleles detected excluding S V and S X (see Additional file 1: Table S1). S C and S 8-alleles could only be distinguished by PCR-amplifying a fragment containing the 358 bp insertion present in SFB C and absent in SFB 8 [17, 18]. As a result, S 8 was only detected in four Hungarian cultivars (‘Ceglédi óriás’, ‘Effect’, ‘Gönci Magyar’ and ‘Szegedi Mammut’) while S C was found in 38 cultivars from many different origins (Table 3). Some already identified S-alleles (S 3, S 20 and S 24) were not present in previously S-genotyped control cultivars (Table 3). For instance, putative S 3 detected in ‘Harlayne’ and ‘Henderson’ has an estimated size for the first S-RNase intron that matches that of the ‘Sunglo’ S 3-allele identified by Vilanova et al. [19] and nucleotide sequences from the three cultivars aligned by CLUSTALW showed 99% identity (data not shown). Moreover, according to the pedigree ‘Sunglo’ is the ‘Harlayne’ male parent [32]. S 20 could not be amplified for the first S-RNase intron but the estimated size as well as the nucleotide sequence (99% identity) (see Additional file 2: Table S2) of the second one matches almost perfectly that of the S 20-allele reported by Zhang et al. [29]. Interestingly, S 19 could only be weakly amplified from ‘Mari de Cenad’ with PruC2/PruC6R primers and the fragment size observed in this work and that previously reported by Halázs et al. [20] is similar to S 20, suggesting that both S-alleles might be the same. Unfortunately, this fragment could not be finally sequenced to confirm this point. S 24 was only detected in ‘Ezzine’ but the estimated sizes for the first and second S-RNase introns, and the nucleotide sequence of the first one (248 bp), matches that S 24-allele reported by Gu et al. [31] (see Additional file 2: Table S2).

M-locus genotyping

The apricot M- and M’-loci were previously identified in cultivars ‘Canino’ and ‘Katy’ as carrying S-locus unlinked PPMs (m and m’, respectively) conferring SC [25, 26]. Up to 35 SSRs were found in the peach syntenic region of the M-locus (defined according to the ‘Canino’ genetic map) but just a five of them could be finally mapped due to several reasons (i.e. multiloci patterns, monomorphism, non- amplification, etc.). Two F1 populations segregating for the PPMs obtained by crossing the self-incompatible cultivar ‘Goldrich’ with ‘Canino’ (‘G × Ca’) and ‘Katy’ (‘G × K’) were used for mapping [25, 26]. Three SSRs were mapped using ‘G × Ca’ (PGS3.71, PGS3.62 and PGS3.96) and two with ‘G × K’ (PGS3.22 and PGS3.23). In this work, linkage phases were determined by using JoinMap 3.0 software [33] and SSR-alleles in coupling were subsequently used to define haplotypes for the M-locus in ‘Goldrich’, ‘Canino’ and ‘Katy’. These 5 SSR markers were also PCR-amplified from M- and M’-loci homozygote individuals derived from the self-pollination of ‘Canino’ [CC-77 (MM) and CC-67 (mm)] and ‘Katy’ [K06–17 (m’m’)] confirming estimated linkage phases and adding information to complete the haplotypes. Lastly, two additional SSRs identified from apricot genomic sequences located at the M-locus (AGS.20 and AGS.30) were also incorporated using the same procedure. SSR heterozygosity was estimated according to Nei [34] (Fig. 1; see Additional file 3: Table S3). As first finding, all SSR-alleles in coupling with the ‘Canino’ PPM m were found to be also in coupling with the ‘Katy’ PPM m’. In other words, both cultivars share the same pollen-part mutated haplotype designated hereinafter by m 0. As a whole, these results lead to define initially four different M-haplotypes (m 0, M 1, M 2 and M 3) according to the genotypes established for ‘Canino’, ‘Katy’ and the reference cultivar ‘Goldrich’ (Fig. 1).

Fig. 1

Apricot M-locus haplotypes structure. The peach syntenic region at the distal end of chr. 3 (black box) comprising the M-locus (grey box) is zoomed twice and shows SSRs (dashed lines) PCR-amplified in apricot cultivars ‘Goldrich’ (G), ‘Canino’ (Ca) and ‘Katy’ (K). SSR positions in peach genome (Kb) and allele sizes (bp) determined in apricot are indicated. White, black, diagonal striped and grey thick lines represent apricot m 0, M 1, M 2 and M 3 haplotypes, respectively. SSR anchoring positions are shown in centimorgans (cM) (boxed numbers) according to the available mapping populations (‘G × Ca’, ‘K × K’ and ‘G × K’)

These seven SSR markers were subsequently PCR-genotyped in the remaining accessions. Structures of up to 34 additional M-haplotypes were then statistically inferred (see Methods for details) while m 0, M 1, M 2 and M 3 were fully confirmed (see Additional file 4: Table S4). For the sake of simplicity and to facilitate the graphical representation of their relative frequencies, M-haplotypes were grouped into ‘main classes’ when they differ in no more than three SSR alleles (resulting in a total of 20 M-haplotype ‘main classes’: m 0 and M 1 to M 19). Thus, M-haplotypes were designed by two sub-indexes, the first one corresponds to the ‘main class’ and the second defines the subtype (see Additional file 4: Table S4). Structures of the M-haplotypes corresponding to the cultivars ‘Kech-pshar’ and ‘Fergani’ could not be inferred but according to the SSR allele sizes all four are distinct (data not shown). Phylogenetic reconstructions from Jaccard’s and Bruvo’s distances show that some ‘main classes’ clustered close together suggesting common ancestry such as, for instance, M 7 and M 8; M 2, M 4 and M 16 or m 0, M 1 and M 13 (Fig. 3). In this latter group, the nearest haplotypes to the majority haplotype m 0–0 were M 1–2, M 1–1 and M 1–0 showing Bruvo’s distance values of 0.18, 0.22 and 0.26, respectively.

Fig. 3

Clustering analysis of apricot M-locus haplotypes based on genetic distances. a Clustering obtained by Neighbor-Joining algorithm using Jaccard’s distance. b Clustering obtained by Neighbor-Joining algorithm using Bruvo’s distance. Colors represent geographic areas where the distinct M-locus haplotypes were detected (see legend)

Mutations conferring self-compatibility in apricot

The S C-haplotype known to confer SC [17] was found in 38 accessions (being homozygous in 7) (Table 3). Thirty out of the 38 were self-compatible, 6 have undetermined phenotype and two were male-sterile. The m 0-haplotype associated with SC by Zuriaga et al. [25, 26] was detected in a total of 19 accessions (being homozygous in 8): 8 were previously shown to be self-compatible, 5 were undetermined and one was classified as self-incompatible (‘Cow-2’) since it did not produce any fruit after self-pollination (Tables 2 and 3). As a whole, all self-compatible accessions analyzed have at least one of these two haplotypes already known to confer SC (S C and/or m 0) except for ‘Harlayne’, ‘Henderson’, ‘Shalah’ and ‘Mariem’. The number of self-compatible haplotypes carried by each accession varied from 0 to 4. A total of 21 accessions do not carry neither S C nor m 0, 26 carry only one self-compatible haplotype (corresponding to S- and M-genotypes S C/x or m 0/x), 13 carry up to two (S C/S C, m 0/m 0 or S C/x-m 0/x), 1 carry three (i.e. Gavatxet S 20/S C-m 0/m 0) and 6 carry four (S C/S C-m 0/m 0) (Table 3).

Segregation of the S-genotypes in self-progenies obtained from two self-compatible cultivars carrying a single copy of the m 0-haplotype (‘Portici’ and ‘Corbató’) supported a mutation outside the S-locus as the cause for the phenotype (Table 4). Progenies from self-compatible cultivars (‘Dulcinea’ and ‘Bebecou’) not carrying the m 0-haplotype were used as controls to confirm this finding (Table 4). Similar results were observed in other accessions (with fewer offspring) carrying (‘Cow-1’ and ‘Cristalí’) and not carrying (‘Ezzine’ and ‘Ninfa’) the m 0-haplotype (data not shown). Moreover, distortion ratios detected in the segregation of SSR markers tightly linked to the M-locus (AGS.20 and PGS3.23) point out that the mutation is located at the M-locus (Table 4).

Table 4

Segregation of S-RNase alleles in controlled self-pollinations

Cultivar S-genotype	Progeny S-genotypes			Totala	Ratio 1:1b	Ratio 1:2:1b	Ratio 1:1c
	A	H	B
Bebecou (S 6/S C)	0 (S 6/S 6)	33 (S 6/S C)	45 (S C/S C)	96 (78)	1.85 (0.17)	53.8 (0)	n.d.d .
Dulcinea (S C/S 2)	36 (S C/S C)	37 (S C/S 2)	0 (S 2/S 2)	104 (73)	0.34 (0.56)	37.9 (1e-8)	n.d.
Corbató (S 5/S 2)	10 (S 2/S 2)	17 (S 5/S 2)	9 (S 5/S 5)	44 (36)	---	0.17 (0.92)	0.17 (0.68) [24]
Portici (S 2/S 20)	16 (S 2/S 2)	24 (S 2/S 20)	4 (S 20/S 20)	59 (44)	---	6.91 (0.03)	0.90 (0.35) [79]

aNumber of individuals from the self-progeny and total analyzed (between parentheses) are indicated

bχ2 and (P) values for S-genotypes expected ratios considering a single mutation unlinked (1:2:1) or linked to the S-locus (1:1)

cχ2 and (P) values for PGS3.23 (in ‘Corbató’ self-progeny) and AGS.20 (in ‘Portici’ self-progeny) genotypes expected ratios considering a single mutation located at the M-locus (1:1). No homozygote for the SSR-allele in repulsion-phase linkage with the m 0-haplotype was detected in either case. [n] indicates number of individuals tested

dn.d. Not determined

Interestingly, S-genotypes segregation found in the ‘Portici’ progeny might also indicate the presence of another mutation affecting the S 2-haplotype, since the number of S 2 S 2 genotypes was unexpectedly high (Table 4). Analysis of genomic DNA fragments containing the complete sequence of S 2-RNase and SFB 2 alleles from ‘Portici’ revealed only one mismatch (A/G) within SFB 2 located at position 1.296. This change leads to a non-synonymous substitution (lysine by arginine) in the hypervariable region HVb.

Lastly, S- and M-genotyping of ‘Canino’ selected clonal sibs (14–4, 14–6 and 9–7) [35] did not reveal differences with ‘Canino’ (Table 3), pointing out that S C and m mutations are shared by all these four accessions.

Geographical distribution of self-(in)compatibility

For the sake of the analysis, cultivars and accessions studied in this work were grouped as belonging to four big geographic areas according to the country of origin, pedigree (Table 1) and data about dissemination [36-38]: North America (NA), Western Europe (WE), Eastern Europe (EE) and Southern Europe/North Africa (SE/NAf). Nonetheless, this classification is obviously arbitrary and exhibits some inconsistencies. This is particularly true for the WE and SE/NAf groups considered separately but clearly interrelated. Attending to geographic criteria, cultivars from Italy, France and Spain (Valencia region) were grouped into the first one while those from Tunisia, Greece and Spain (Murcia region) into the second. Spanish cultivars were divided into two subgroups that might reflect some differences according to [49]. Parentage relations among some cultivars from WE and SE/NAf groups (i.e. ‘Canino’ with ‘Ouardi’ and ‘Sayeb’) underscore the inaccuracy of this classification but, assuming these limitations, it produces more balanced groups without affecting data interpretation. Thus, SC phenotype is distributed across all four groups but frequencies varied from 2/19 (SI/SC) in WE to 9/16 in NA (Fig. 2). Similarly, blooming time types from early/mid-early to mid-late/late are present in all four groups, but frequency ranged from 9/24 in WE to 1/16 in NA (Table 1). In general it could be said that most early/mid-early classified accessions are self-compatible (13/15) but this phenotype is also generally found in many mid-late/late classified accessions (12/17). Therefore, even considering the reduced number of scored accessions and the limited phenotype data, it seems that there is no correlation between SC and early blooming in cultivated apricots.

Fig. 2

S- and M-locus haplotypes distribution according to geographic areas. Apricot accessions analyzed in this study were grouped in four arbitrarily defined geographic areas represented by the bottom map: Western Europe (WE), North America (NA), Southern Europe and North Africa (SE/NAf) and Eastern Europe (EE). Accordingly, relative frequencies for S- and M-haplotypes (pie charts) are shown for each area (clonal sibs from ‘Canino’ were excluded from estimations). For the sake of simplicity M-locus haplotypes are represented in their ‘main classes’. Question marks (?) designate not defined haplotypes. Number of accessions showing self-(in)compatible phenotype are encircled (green color means SC, red SI and blue undetermined phenotype)

Frequencies of the S-alleles also varied between the different regions analyzed. For instance, S 2 and S C are present in all four groups but, on the contrary, S 8 and S 9 are only present in EE, S 3 in NA and S 7 in SE/NAf (Fig. 2). Regarding S-alleles distribution it is also important to take into account synonymies and homonymies. In the first case, BLASTN has revealed that some of the S-alleles initially found and numbered from European accessions have also been detected in Chinese cultivars but named differently. For instance, a 659 bp fragment containing the second intron of the S 17-RNase present in the North-American cultivar ‘Aurora’ (as well as in ‘SEO’, ‘Orange Red’ and ‘Henderson’) shows >99% identity with those from S 44 and S 9-RNases identified from the Chinese cultivars ‘Shailaiyulvke’ and ‘Xinshiji’, respectively, suggesting these three S-alleles to be the same. Similar examples are reported for S 1, S 4, S 9, S 12, S 16 and S 20-RNases. The other way around, BLASTN also points out homonymies. For example, S 9 derived from the Hungarian cultivar ‘Ceglédi óriás’ is different to that S 9 reported from the Chinese cultivar ‘Xinshiji’ (see Additional file 2: Table S2). BLASTN has also detected significant similarities between S 13 reported in accessions from Armenia (including ‘Shalah’), Turkey, Tunisia and Morocco [20-22] and S 5, previously found in Southern Spanish cultivars [19], and among S 12, detected in Turkey and Tunisia [20, 22], S 24 and S 4 (see Additional file 2: Table S2).

Distribution of the 38 inferred M-haplotypes is not uniform either. The M 1 ‘main class’ is present in all four groups, m 0 is only found in WE and NA and others, such as M 12 and M 2, were exclusive for EE and NA groups, respectively (Figs. 2 and 3).

Discussion

S-alleles and the apricot diffusion from Armenian and Chinese centers to Europe

Most apricot S-alleles identified in this work had been previously reported except for S 22, S V and S X [19, 27–30]. Geographical distribution of these S-alleles is heterogeneous including a few widely disseminated and others restricted to certain areas. These results are in agreement with previous works. For instance, S 2 and S C had already been found in accessions from all countries analyzed in this work but also in Turkey [20], Morocco [21], Afghanistan [18], Iran [40], etc. On the contrary, S 1 is only present in NA and in a few accessions from WE. Meanwhile, S 10 to S 14-alleles, suggested to be from Armenian origin by Halázs et al. [20], have been detected in Turkey, Tunisia, Morocco and in a few materials from EE but they were thought to be absent in WE [20–22, 41]. However, S 5 suggested to be the same that S 13 connects (along with S 6) Armenian, Eastern-Turkish and Moroccan accessions [20-22] with Southern-Spanish accessions [19, 42] supporting the apricot Southwest-Mediterranean diffusion route, from the Irano-Caucasian gene pool, proposed by Bourguiba et al. [38].

The heterogeneous distribution of the apricot S-alleles is a consequence of the crop domestication and diffusion and this process has also produced a genetic bottleneck detected in neutral SSR markers [38]. Though not quantified this loss of genetic diversity can also be grossly appreciated at the S-locus. For instance, China is not only considered a center of origin but also a probable domestication center for the species [38] and at least 20 different S-alleles have been found in just 30 Chinese cultivars [29, 30]. This variability is considerably much higher than that detected in European accessions and, interestingly, only 5 of these S-alleles (S 9, S 14, S 17, S 20 and S 24) seem to be represented in apricot cultivars outside China. In addition, other S-alleles detected in far-east region cultivars such as S 15/S 18 in ‘Kech-pshar’ (Uzbekistan), S 10 in ‘Harmat’ (Armenia) [27] and S V/S X in ‘Fergani’ (former USSR) are not found in European accessions. S-alleles can therefore be a useful tool to study apricot dissemination but the establishment of a unified nomenclature would be advisable to favor this goal.

Pollen-part mutated m-haplotype origin and dissemination

The two S-locus unlinked PPMs (m and m’) conferring SC in apricot cultivars ‘Canino’ and ‘Katy’ were found to be located in an overlapping region at the distal end of chr. 3 [26]. This finding suggested that both PPMs might be the same. However, in addition to their different geographic origins, ‘Canino’ and ‘Katy’ exhibited a high value for the Nei’s genetic distance (0.83) estimated on the basis of 85 SSR markers distributed across the whole genome [26]. Thus, initially it was hypothesized that these two PPMs originated independently but might affect the same gene. Nevertheless, a more detailed analysis of a set of SSR markers linked to the M-locus has allowed us to determine that haplotypes in coupling with the PPMs have the same structure in both cultivars. Thus, pollen-part mutated m- and m’-haplotypes previously associated with SC in ‘Canino’ and ‘Katy’ cultivars [25, 26] are now considered to be the same (m 0). Moreover, in this work the m 0-haplotype has been detected in 17 additional accessions (excluding ‘Canino’ clonal sibs) mainly Spanish (12 in total) but also from USA, Australia, France and Italy. Fifteen of them were confirmed as self-compatible (exceptions were ‘Cow-2′ described above as well as ‘Gandía’, ‘Gavatxet’, ‘Manrí’ and ‘Martinet’ with undetermined phenotype). The m 0-haplotype was frequently accompanied by the S C-allele in confirmed self-compatible accessions (7 cases), which might suggest that mutations conferring SC tend to accumulate once the system is broken due to relaxed selective constraints [43]. Following this reasoning and considering the higher prevalence of the S C-allele in confirmed self-compatible accessions (74%) compared with the m 0-haplotype (23%), it might also be speculated that PPM m arose in a self-compatible accession carrying S C. However, against this hypothesis is the fact that the m 0-haplotype alone was also found in 6 confirmed self-compatible accessions and, therefore, further analyses are necessary to clarify this point. This latter group includes ‘Portici’ (S 2/S 20-M 1–4/m 0–0), previously shown to lack the S C-haplotype but which causative mutation was still unknown [44].

Beside the m 0-haplotype, 37 additional M-haplotypes were identified by SSR analysis being grouped in 19 ‘main classes’. Microsatellite haplotype distances were analyzed by two alternative clustering methods to validate results. The first one relies on the proportion of shared alleles assuming independence and ignoring mutational processes which can bias distances, particularly when loci are highly polymorphic. This method was based on the similarity coefficient for binary data developed by Jaccard [45]. The second takes into account the stepwise mutation model considering higher likelihood for small than for large changes in microsatellite repeat number and it is based on the Bruvo’s distance [46]. Results obtained with both methods were equivalent but the second draws a more accurate phylogeny, where ‘main classes’ m 0, M 1 and M 13 group close together suggesting a not too distant common ancestor. Remarkably, the closest haplotypes to m 0–0 (M 1–0, M 1–1 and M 1–2) are widely distributed geographically, while m 0–0 seems to be restricted to NA and WE, and M 13 to EE. The situation is different for the mutated S C-allele present in all geographic areas studied [18-21] whereas the ancestor S 8-allele was only detected in Hungarian and Turkish cultivars [20]. Altogether, these results suggest that the m 0-haplotype arose later in time than S C since its distribution is much more limited. It can also be hypothesized that the last common ancestor for m 0, M 1 and M 13 spread throughout different regions and mutated somewhere in WE. Nevertheless, in general, it has not been observed a significant correlation between geographic distribution and clustering of the M-haplotypes (see Fig. 3), and this may be due to the frequent crosses between genotypes from different groups [39].

SC and genotyping: uncovering accessions carrying unanalyzed mutations

Part of the materials analyzed in this work could not be self-pollinated because adult trees were not available. However, though phenotype could not be directly assessed in these cases, according to the almost perfect association observed between SC and S C/m 0-alleles it should be expected that most, if not all, accessions carrying whatever of these two alleles were self-compatible. However, some incongruities between phenotype and genotype need to be highlighted. For instance, in spite of having the m 0–0 haplotype self-pollination of ‘Cow-2’ (S 20/S 31 -M 1–0/m 0–0) did not produce any fruit. This behavior could be due to fruit-set problems, as exemplified by the self-compatible cultivar ‘Búlida’ [47] which setting was also nearly null. Alternately, it may also be possible that m 0–0 in ‘Cow-2’ does not carry the expected PPM but further analyses are needed to check this point. On the contrary, four cultivars classified as self-compatible do not carry neither the S C nor the m 0–0 haplotypes: ‘Mariem’ (S 7/S 20-M 1–0/M 8–2), ‘Shalah’ (S 5/S 11-M 8–1/M 19), ‘Harlayne’ (S 3/S 20-M 2–0/M 9) and ‘Henderson’ (S 3/S 17-M 2–0/M 16). Again, erroneous phenotyping can lead to discrepancies. In fact, ‘Mariem’ progeny could not be fully confirmed as derived from self-pollination (Table 2), ‘Shalah’ was reported as self-compatible by Burgos et al. [47] but scored as self-incompatible by Halázs et al. [20] and information describing ‘Henderson’ as self-compatible just comes from nursery catalogs. So, phenotype can be questioned in all these three cultivars. However, SC in ‘Harlayne’ is fully confirmed since it has been tested by several independent sources [47, 48]. Regarding genotypes, all the S-alleles carried by these cultivars can be frequently found in self-incompatible accessions and therefore expected to be ‘functional’. Something similar happens with M 1–0 and M 2–0 haplotypes but M 8–1 and M 8–2 belong to the M 8-haplotype ‘main class’ always detected in self-compatible accessions (a total of 7 in this work and mostly East-European). However, no progenies are available and therefore it could not be tested whether M 8 is mutated (non-functional) or not. As for M 16 and M 19, they were detected only once in the set of accessions analyzed. Overall, there are still no evidence supporting the nature of the mutation/s conferring SC in these cases and neither S/M-loci unlinked mutations nor other mutations affecting S- and M-loci can be discarded.

Lastly, genetic analysis suggests the presence of a SNP mutation within SFB 2 HVb region in ‘Portici’ that could also be associated with SC. It could be speculated that a single non-synonymous change within a SFB hypervariable region might alter its specificity, since these domains (strongly hydrophobic and under positive selection) were already suggested to have a role in the specific recognition of S-RNases [49]. Interestingly, this change has also been found in ‘Katy’ SFB 2 where a similar distortion was detected in the self-progeny [26]. However, further research is needed to validate this observation and other possibilities cannot be discarded (i.e. mutations affecting pollen viability).

Forces selecting for self-compatibility in apricot

SC is therefore quite common in apricot but distribution is not uniform across geographic areas. SI is the prevalent phenotype in three out of the four major eco-geographical groups for apricot (centers of origin): Central Asian, Irano-Caucasian and Dzhungar-Zailij and also in the later proposed Chinese group [36]. In fact, studies about Chinese apricots do not report self-compatible cultivars [29, 30, 50]. On the contrary, SC predominates in the European group but some disequilibrium can also be observed. Among the materials analyzed in this work, SI is frequent in commercial North-American cultivars (two thirds) but unusual in West and East-European countries (one fifth). This might point out a non-European SI ancestral donor for North-American cultivars apricots as previously suggested for PPV resistance trait [51].

It is generally accepted that apricot genetic diversity decreases from east to south-west [36] and, in this context, it is questionable whether SC might be one of the causes. Recent works suggest that a substantial part of this loss is independent of the selection impact due to the domestication bottleneck, being more related to apricot diffusion routes [38]. In general, self-fertility would be advantageous under unfavorable pollination conditions including early blooming [52] but this correlation has not been observed in this study. However, we cannot discard that SC was originally selected from early-flowering apricot types. In fact, in cherry (Prunus avium) only a very few cultivars are reported to be naturally self-compatible and these include two exceptionally early-flowering cherries, ‘Cristobalina’ and ‘Kronio’ [53, 54]. Moreover, they were selected independently since each one carry a different mutation conferring SC, ‘Cristobalina’ in a non-S-locus pollen modifier [24] and ‘Kronio’ in the SFB 5 allele [54]. On the other hand, in spite of the benefits linked to SC for growing stone fruits (increased yield, removing of interspersed pollinators, etc.) it was not usually considered as a target trait. Some morphological and agronomic traits are putatively linked to the S- (fruit shape) and M-loci (polycarpel and flower color) in Prunus [55] but it seems unlikely that SC was brought in by linkage drag with these particular traits. Thus, SC in apricot is more likely the indirect result of selective breeding on major traits, such as ecological adaptation [36], and only recently it has become a breeding objective by itself. In any case, SC must have had an effect on reducing diversity since S C S C homozygote cultivars (potentially resulting from self-pollination events) are frequently found in most European countries.

In this work we have confirmed two independent PPMs (S C and m) as the main causes for SC in apricot. According to [16] a total of 29 mutations conferring SC have been identified in Prunus. Most of them affect S-locus (27) and only two correspond to mutated modifiers, the apricot mutated m 0-haplotype and a non S-locus PPM in sweet cherry [17, 24]. Indeed, these two are located at the distal end of chr.3 but it is still unknown whether they affect different genes [26]. Within the supposedly complex GSI mechanism it is noteworthy this low number of mutations affecting modifiers. One possible explanation is that other modifiers participating in the control of the GSI system may be redundant and when mutated do not confer SC. Furthermore, loss of function of some factors might lead to SI as may be predicted for SLFL2 [14]. Lastly, in spite of the high number of accessions already evaluated in all Prunus species, it can also be argued that a more exhaustive screening might reveal novel mutations. Indeed, in the light of the results of this work, it is conceivable that additional mutations conferring SC may remain ‘hidden’ in the large amount of unanalyzed plant material.

Conclusions

It has been shown that apricot S-allele geographical distribution is heterogeneous. In addition, it seems that the loss of genetic diversity, hypothesized for this crop during the dissemination process, affected not only neutral markers but also the S-locus. Indeed, only a few S-alleles present in self-incompatible cultivars from China (the domestication center of the species) can be found in European cultivars. However, S C, one of the most diffused S-alleles and the main cause of SC in apricot has not been detected in China. Apart from the PPM associated to S C, other S-locus unlinked PPMs conferring SC have been found in apricot cultivars ‘Canino’ and ‘Katy’. At odds with the initial hypothesis, we have shown that both are associated with the same m 0-haplotype, in spite of the genetic distance between the two cultivars. Moreover, this haplotype has been surprisingly detected in other Spanish, French, Italian and North-American accessions. Its origin is still a matter of speculation but the phylogenetic analysis supports that m 0 arose later in time than S C from a widely distributed M-haplotype. As a whole, self-(in)compatibility in apricot seems to rely mainly on the S and M loci. However, other self-compatible mutants putatively carrying different mutations have been identified in this work. These results deepen our knowledge on the self-(in)compatibility trait in stone fruits and open new questions that deserve future research.

Methods

Plant Material and self-pollination test

Sixty seven apricot cultivars and accessions from diverse geographic origins were used in this study (Table 1). Most are currently kept at the collection of the Instituto Valenciano de Investigaciones Agrarias (IVIA) in Valencia (Spain). Part of this collection was kindly provided by Frutales Mediterráneo S.A. (FM) company, by the Consellería de Agricultura, Pesca y Alimentación (CAPA) and by the Ministerio de Agricultura, Alimentación y Medio Ambiente of Spain (MAGRAMA). Other samples, already used in previous works [56], were provided by the Departamento de Mejora y Patología Vegetal del CEBAS-CSIC in Murcia (Spain) and by the University of St. Istvan (Budapest, Hungary).

Trees from different cultivars and accessions (Table 2) were tested for SC by self-pollination in the field. Before anthesis, insect-proof bags were put over several branches, containing approximately 200–250 flower buds in total per cultivar ad minimum, to prevent cross pollination. Subsequent fruit set was recorded and fruits collected about 3 months later. Seed-derived embryos were dissected from the rest of the seed tissue and stored at −20 °C.

DNA extraction

Two leaf discs were collected from each accession, frozen in liquid N2 and stored at −80 °C before DNA isolation. Genomic DNA was extracted following the method of [57] with slight modifications: the ratio of fresh leaf tissue/CTAB buffer was reduced to 10 mg/250 μl and nucleic acids were recovered by precipitation with 1× volume of isopropanol. DNA quantification was performed by NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE) and integrity was checked on 1% agarose gel by comparison with lambda DNA (Promega, Madison, WI, USA). Embryo DNA was extracted by incubating for 10 min at 95 °C with 20 μL of TPS (100 mM Tris-HCl, pH 9.5; 1 M KCl; 10 mM EDTA) isolation buffer [58].

S-locus genotyping

Apricot S-alleles were identified by PCR-amplifying fragments comprising first and second introns of the S-RNase as well as the 5′-UTR SFB intron, respectively. Genomic DNA isolated from cultivars and accessions listed in Table 1 was used as PCR template. PCRs were performed in a final volume of 20 μL containing 1 × DreamTaq buffer, 0.2 mM of each dNTP, 20 ng of genomic DNA and 1 U of DreamTaq DNA polymerase (Thermo Fisher Scientific, Waltham, MA, USA). Previously developed primers designed from conserved regions of Prunus armeniaca S-RNase genomic sequences, SRc-F and SRc-R [19, 58] and from Prunus avium S-RNase-cDNA sequences, Pru-T2 and Pru-C2R [59], were used to amplify the first intron (see Additional file 5: Table S5) in all four possible combinations. The amplification was carried out using a temperature profile with an initial denaturing of 95 °C for 2 min; 30 cycles of 95 °C for 45 s, 52 °C for 1 min and 72 °C for 1 min 30 s; and a final extension of 60 °C for 30 min (UNO96, VWR, Radnor, PA, USA). Each PCR was performed by the procedure of [60] using three primers: the specific forward primer with M13(−21) tail at its 5′ end at 0.4 mM, the reverse primer at 0.8 mM, and the universal fluorescent-labeled M13(−21) primer at 0.4 mM. Allele lengths were determined using an ABI Prism 3130 Genetic Analyzer with the aid of GeneMapper software, version 4.0 (Applied Biosystems, Foster City, CA, USA).

The second intron was amplified using two sets of primers designed from Prunus avium S-RNase-cDNA sequences [61]: Pru-C2/Pru-C4R and Pru-C2/Pru-C6R (see Additional file 5: Table S5). PCRs were performed using 0.25 μM of each primer and the program previously described by [62] to amplify long PCR products. PCR products were electrophoresed in 0.8% (w/v) agarose gels using 1 x TBE (89 mM Tris, 89 mM boric acid, and 2 mM EDTA (pH 8.0)) buffer, stained with RedSafe Nucleic acid Staining Solution (iNtRON Biotechnology, Korea) and visualized under UV light. Molecular sizes of amplified fragments were estimated using GeneRuler 100 bp Plus DNA ladder (Thermo Fisher Scientific).

The 5′-UTR SFB intron was also amplified using the degenerate primer pair (F-BOX5′A/F-BOXintronR) developed by Vaughan et al. [63] from sweet cherry sequences (see Additional file 5: Table S5). PCR components, thermo-cycler conditions and detection procedure were identical to that described above for the first S-RNase intron. Two additional primers (RFBc-F/SFBins-R), designed from the consensus sequence of the Prunus SFB alleles [17], were used to amplify the SFB C insertion from genomic DNA of several apricot cultivars in order to distinguish S C and S 8-alleles (see Additional file 5: Table S5).

Seven SSR markers comprised within (or flanking) the M-locus were genotyped: PGS3.71, PGS3.22, PGS3.62, PGS3.23 and PGS3.96 [25, 26] and AGS-20 and AGS-30 (this work) (see Additional file 3: Table S3). SSR amplifications were performed in a GeneAmp® PCR System 9700 thermal cycler (Perkin–Elmer, Freemont, CA, USA) in a final volume of 20 μL, containing 1 × DreamTaq buffer, 0.2 mM of each dNTP; 20 ng of genomic DNA and 1 U of DreamTaq DNA polymerase (Thermo Fisher Scientific). Regarding primers, PCRs were performed as described above by the procedure of [60]. The following temperature profile was used: 94 °C for 2 min, then 35 cycles of 94 °C for 45 s, 50–60 °C for 1 min, and 72 °C for 1 min and 15 s, finishing with 72 °C for 5 min. Allele lengths were determined using an ABI Prism 3130 Genetic Analyzer with the aid of GeneMapper software, version 4.0 (Applied Biosystems).

Sequence analysis of PCR products containing S-RNase introns

PCR products containing the first and second introns of S 20 and S 24 -RNase alleles previously obtained from genomic DNA of ‘Ezzine’ (S C S 24), on one side, and from cultivars ‘Cow-2’ (S 20 S 22), ‘Harlayne’ (S 3 S 20), ‘Cristalí’ (S 20 S C), ‘Gavatxet’ (S 20 S C), ‘Mariem’ (S 20 S 7), ‘Perla’ (S 20 S 2), ‘Portici’ (S 20 S 2), ‘Stella’ (S 20 S 6), ‘Tadeo’ (S 20 S C), ‘Veecot’ (S 20 S 2) and ‘Velázquez’ (S 20 S 5), on the other, were sequenced to check their identities. Similarly, PCR products associated with S 3, S 17 and S 22 -RNase alleles obtained from cultivars ‘Harlayne’ (S 3 S 20), ‘Henderson’ (S 3 S 17), ‘Orange Red’ (S 6 S 17), ‘SEO’ (S 6 S 17), ‘Cow-1’ (S 1 S 22) and ‘Cow-2’ (S 20 S 22) were also sequenced. Primer combinations SRc-F/SRc-R and Pru-C2/Pru-C4R were used in all cases for the first and second introns respectively except for S 24 (Pru-T2/SRc-R) and S 22 (Pru-C2/Pru-C6R). PCR products were electrophoresed in 0.8% or 2% (w/v) agarose gels (second or first intron, respectively) stained with RedSafe Nucleic acid Staining Solution (iNtRON Biotechnology). Molecular sizes of the amplified fragments were estimated using GeneRuler 100 bp Plus DNA ladder (Thermo Fisher Scientific). Fragments were extracted and purified from the agarose gels using the Zymoclean Gel DNA Recovery kit (Zymo Research, Irvine, CA, USA). Sequences were determined automatically using an ABI PRISM 3100 Genetic Analyzer (Applied Biosystems) and the Big Dye Terminator Cycle Sequencing Kit v3.1 (Applied Biosystems) following the manufacturer’s instructions. Sequences were assembled and edited with the Staden package v.2.0 [64]. Multiple alignments were made with CLUSTALW v.1.82 [65] and edited in BioEdit v.7.0.9 [66]. Homology searches were performed against the NCBI Genbank database using the BLASTN program [67].

S2-locus sequence analysis

Specific primers designed from the apricot S 2-haplotype sequence [17] were used to amplify genomic fragments containing the complete S 2-RNase (Sf-Hap2/Sr-Hap2) and SFB 2 coding sequences (FBf-Hap2/FBr-Hap2) (see Additional file 5: Table S5) using ‘Portici’ (S 2 S 20) genomic DNA as template. PCR conditions and methods for isolating and sequencing (using the three additional primers SFBc-F, FBF5 and FBF6) these bands (see Additional file 5: Table S5) were the same reported above for fragments containing the S-RNase second intron.

Clustering analysis

Reference m 0–0 , M 1–0 , M 2–0 and M 3-haplotypes were established using genetic maps from ‘G × Ca’, ‘K × K’ and ‘G × K’ populations through the automatic determination of linkage phases by JoinMap 3.0 [33]. Remaining M-haplotypes were inferred from SSR genotypes by comparing with the references and confirmed by using the EM algorithm [68] implemented in PowerMarker V3.25 software [69]. When available pedigree was used to confirm assignments. Similarities between M-haplotypes were estimated by using Jaccard’s similarity coefficient [45] through the Phylip 3.62 package [70] and Bruvo’s genetic distance [46] through a hand-made phyton script. Bootstrapped data matrices were obtained using the Phyltools 1.32 software [71] and used to test the stability of the neighbor joining trees constructed with the Phylip 3.62 package [70]. HyperTree software [72] was used to visualize the obtained trees.

Fragment sizes (bp) determined for S-RNase and SFB allele introns. Fragment sizes for the first S-RNase and the 5′-UTR SFB introns were exactly determined by an ABI PRISM 3100 Genetic Analyzer and those for the second S-RNase intron were approximately estimated from 0.8% agarose gels. References reporting S-allele molecular sizes for the first time are provided. (DOCX 15 kb)

Main results of homology searches performed by BLASTN against the NCBI Genbank database using Prunus armeniaca S-RNase sequences. Query sizes, Query cover %, Total score, E-value and Identity % are indicated. Putative names were assigned to the S-alleles according to the fragment sizes and sequence analyses. (DOCX 18 kb)

Characteristics of SSR primers developed from peach (PGS) and apricot genome sequences (AGS) located at the M-locus. Repeat motifs, Primer sequences, Scaffold positions, ORF (Prupe), Size range, N° of alleles and Heterozygosity (in the set of accessions analyzed) are indicated. (DOCX 14 kb)

Allele sizes (bp) for the SSR markers comprised into the M-haplotypes. ‘Main class’ as well as the subtype are indicated for every M-haplotype. Absolute frequency in the set of accessions analyzed is also provided. (DOCX 18 kb)

Primers used in this study. Sequences and references are indicated. (DOCX 14 kb)