Next Article in Journal
New Sources, Differentiation, and Therapeutic Uses of Mesenchymal Stem Cells 2.0
Next Article in Special Issue
Seedlessness Trait and Genome Editing—A Review
Previous Article in Journal
Relationship among Low T3 Levels, Type 3 Deiodinase, Oxidative Stress, and Mortality in Sepsis and Septic Shock: Defining Patient Outcomes
Previous Article in Special Issue
Towards Marker-Assisted Breeding for Black Rot Bunch Resistance: Identification of a Major QTL in the Grapevine Cultivar ‘Merzling’
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 4
10.3390/ijms24043932
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

S-Locus Genotyping in Japanese Plum by High Throughput Sequencing Using a Synthetic S-Loci Reference Sequence

by
Afif Hedhly
1,*,
María Engracia Guerra
2,
Jerome Grimplet
1,3 and
Javier Rodrigo
1,3,*
1
Departamento de Ciencia Vegetal, Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA), Avda Montañana 930, 50059 Zaragoza, Spain
2
Área de Fruticultura Mediterránea, CICYTEX-Centro de Investigación ‘Finca La Orden-Valdesequera’, A-V, KM 372, Guadajira, 06187 Badajoz, Spain
3
Instituto Agroalimentario de Aragón-IA2, CITA-Universidad de Zaragoza, 50013 Zaragoza, Spain
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(4), 3932; https://doi.org/10.3390/ijms24043932
Submission received: 31 December 2022 / Revised: 9 February 2023 / Accepted: 10 February 2023 / Published: 15 February 2023
(This article belongs to the Special Issue Advances in Research for Fruit Crop Breeding and Genetics)
Browse Figures
Versions Notes

Abstract

:
Self-incompatibility in Prunus species is governed by a single locus consisting of two highly multi-allelic and tightly linked genes, one coding for an F-box protein—i.e., SFB in Prunus- controlling the pollen specificity and one coding for an S-RNase gene controlling the pistil specificity. Genotyping the allelic combination in a fruit tree species is an essential procedure both for cross-based breeding and for establishing pollination requirements. Gel-based PCR techniques using primer pairs designed from conserved regions and spanning polymorphic intronic regions are traditionally used for this task. However, with the great advance of massive sequencing techniques and the lowering of sequencing costs, new genotyping-by-sequencing procedures are emerging. The alignment of resequenced individuals to reference genomes, commonly used for polymorphism detection, yields little or no coverage in the S-locus region due to high polymorphism between different alleles within the same species, and cannot be used for this purpose. Using the available sequences of Japanese plum S-loci concatenated in a rosary-like structure as synthetic reference sequence, we describe a procedure to accurately genotype resequenced individuals that allowed the analysis of the S-genotype in 88 Japanese plum cultivars, 74 of them are reported for the first time. In addition to unraveling two new S-alleles from published reference genomes, we identified at least two S-alleles in 74 cultivars. According to their S-allele composition, they were assigned to 22 incompatibility groups, including nine new incompatibility groups reported here for the first time (XXVII-XXXV).

1. Introduction

In flowering plants, Gametophytic Self-Incompatibility (GSI) is the most widespread and ancient post-pollination mechanism of non-self and self-recognition that inhibits the growth of self-pollen tubes in the style, preventing self-fertilization [1,2]. Self-incompatibility allows maintaining genetic diversity and avoiding the potentially negative effects of inbreeding [3]. It is controlled by a single genetic locus subjected to strong negative frequency selection, i.e., rare alleles have a mating advantage, [4] that would lead to an increase in the number and diversity of S-alleles [5]. The S-locus consists of two tightly linked genes that determine both pollen and pistil specificity [6]. In Prunus sp., as in other species of the Rosaceae, pollen allele specificity is determined by an F-box protein (SFB) that is expressed in the pollen grain [7], and pistil allele specificity is determined by a ribonuclease (S-RNase) that is expressed in styles [8,9], and inhibits the growth of self-pollen tubes [10]. The genotyping of S-locus in both natural and domesticated plant species has confirmed the high allelic diversity and sequence polymorphism [11,12].
In fruit trees, S-allele genotyping is of great practical and scientific importance, since it allows the adequate design of crosses in breeding programs and the correct choice of pollinators for self-incompatible cultivars in new orchards [13]. In Japanese plum (Prunus salicina and hybrids of P. salicina with other Prunus species), most cultivars are self-incompatible and require cross-pollination to set fruit. Therefore, knowledge of the S-genotype of each cultivar is essential for the design of new orchards with an adequate proportion of inter-compatible cultivars and to solve problems of low fruit set related to lack of pollination [14,15]. It also allows the characterization of allelic series in natural populations to advance our understanding of gene flow and the evolution of mating barriers in plants, and their ecological and evolutionary implications [16,17].
In fruit tree species, self-incompatibility and inter-incompatibility relationships between cultivars were initially determined on the basis of controlled crosses followed by characterization of fruit set under orchard conditions [18]. Later, incompatibility was determined more precisely by microscopic observation of pollen tube growth through the style in self- and cross-pollinated flowers [19], first in the field and later in semi-in vivo cut flowers in the laboratory (e.g., [20,21,22]). These procedures allowed the characterization of cultivar self-(in)compatibility and the establishment of incompatibility groups by using an arbitrary nomenclature system of S-alleles. Towards the end of the 20th century, the identification of proteins (S-RNases and SFBs) associated with S-alleles allowed the identification of S-genotypes by biochemical methods [23,24]. Likewise, the cloning and sequencing of the genes encoded by the S-locus allowed obtaining the first S-allele sequences upon which different PCR techniques evolved (reviewed in [25,26]). In Prunus and other genera belonging to the Rosaceae, degenerate primers designed from conserved exonic regions within the S-RNase and SFB genes and spanning polymorphic intronic regions allowed trans-specific S-genotyping using only a few primer pairs [27]. In the absence of intronic length polymorphism, other PCR-based techniques were developed, such as the single strand conformation polymorphism (SSCP) technique [28]. Nowadays, the tremendous advances in whole genome sequencing are opening new avenues for less time- and resource-consuming S-genotyping procedures. However, reference genomes, which are currently available for many economically important fruit tree species, cannot be used to accurately align and resequence the S-genotype of individuals due to the extreme polymorphism that characterizes the S-locus. To circumvent this pitfall, several approaches have recently emerged, including methods based on bar-coded amplicons [29,30,31], separate alignment on several reference genomes of the same species [32], and de novo assemblies of either the whole genome, followed by BLAST and phylogenetic analyses [31], or the S-allele itself [33,34]. In apple (Malus domestica), de Franceschi et al. managed to S-genotype 63 resequenced cultivars using 32 apple S-RNase sequences concatenated in a string-like structure as reference sequence [33].
In this work, building on the previously described procedure used in M. domestica [33], we present an improved methodology to use information from the entire range of known S-loci (genomic regions including S-RNases and SFBs), rather than just one of the two genes that form the locus, for successful S-genotyping 88 individual shotgun sequences of Japanese plum cultivars. To circumvent the problem of unspecific high depth of coverage in intronic regions, which can lead to misidentification of alleles, we developed a simple but quantitative three-step validation that relies not only on the depth but also on the breadth of coverage. This methodology allowed us to accurately ascertain the S-genotype of most cultivars, enabling accurate S-genotyping of any individual (i.e., the four alleles) as long as it has an already known S-allele, and to identify individuals with supernumerary S-alleles, or with likely novel and unknown alleles.

2. Results and Discussion

2.1. An S-Locus Based Synthetic Sequence

Since the first two Japanese plum S-RNase alleles were cloned in 1999 from the Japanese cultivar ‘Sordum’ [35], a further 53 S-RNases and 36 SFBs sequences—corresponding to 39 different loci—have been published or submitted directly to the NCBI database (Table 1 and references therein) representing Chinese, Japanese, North American, Israeli, Hungarian, and Tunisian germplasms. As this is a locus involved in mate recognition, a high degree of diversification and high number of alleles was suggested and confirmed in both wild and domesticated plant species [5,36,37]. Further genotyping of germplasm from China—where the species is believed to have originated—will undoubtedly reveal more polymorphism. To confirm the identity of all S-alleles and detect possible duplicates, i.e., redundant or inconsistent-, we performed phylogenetic and pairwise similarity analyses at the nucleotide and amino acid levels for all alleles (data not shown). We detected some redundancies, such as that the S1, S3, S4, S5, and S6 alleles were actually duplicates of the Sa, Sk, Sc, Se, and Sf alleles, respectively (Table 1). There was one case of two different S-alleles receiving the same name (S23-RNase), and both sequences were conserved. Despite the high similarity between S10 and S32 (99.3% identity at the nucleotide level), we could not fully determine whether they are the same allele. There was only 96.0% identity at the amino acid level, differing at the 3 amino acids immediately downstream of the conserved RC4 region. In addition, both arrangements are shared with other S-alleles. Both alleles were also conserved to build the synthetic sequence.
To complete our search for all available information and, if possible, to retrieve full S-loci sequences, we performed a BLASTN search against the three available reference genomes with either short and conserved regions or full S-RNase and SFB sequences. In addition, we searched for highly similar candidate genes that fulfill a tail-to-tail orientation as in all Prunus S-loci, being able to retrieve three full S-loci. The S-RNase and SFB from cultivar ‘Zhongli nº6′ turned out to be 100% identical to Sb-RNase and SFBb of P. salicina, respectively (Figure 1A,B). All other alleles did not appear to match any sequenced allele and were therefore considered novel alleles (named Ssany and Szgl). The alleles recovered from the ‘Sanyueli’ reference genome showed lower than expected trans-specific similarity within Prunus, whereas the S-RNase and SFB recovered from the ‘ZGL’ reference genome were highly similar to S06 S-RNase from P. speciosa and SFB5 from P. tenella, respectively. High similarity between alleles of different species compared to those of the same species—called trans-specific polymorphism—is a common phenomenon within GSI [12,38,39,40]. Sequence analysis of both genes showed typical structural arrangements and size ranges as in other P. salicina alleles [35,41,42,43,44,45,46]. Alignment of the three S-RNases revealed that they have the hypothetical signal peptide sequence (and the first intron at the junction between the signal peptide and the mature protein), the three conserved regions, and the single hyper-variable region (Figure 1C) typical of ribonucleases of the T2 family [24,39]. Alignment of the three SFB alleles indicated that they have the single F-box domain, the two variable regions, and the two hypervariable regions (Figure 1D) typical of Prunus SFB genes [47].
To build a synthetic reference sequence, we retained the three full S-loci sequences, plus 1000 bp upstream and downstream the locus to increase specificity, and the longest genomic sequence of each S-locus gene pair and concatenated them into a rosary-like structure with 1000 bp N spacers (Figure 1E). To have a reference sequence in single strand-like structure and mimic the structure of the S-locus, reverse complements of S-RNase sequences were concatenated in tail-to-tail orientation with SFB alleles within each locus. In the absence of SFB alleles, the S-RNases sequences were also reverse complemented. Alignment of shotgun sequences on available reference genomes, as used for other applications, does not produce a useful alignment at the S-locus region due to the extreme polymorphism that characterizes this locus [33,34]. The technical advance of the method is based on the use of the entire S-locus allelic series of the species—i.e., all S-RNases and SFBs—to construct a synthetic reference sequence on which the resequenced genotypes are aligned. Compared to the previously described methodology based on the use of only the S-RNase sequence [33], our genotyping procedure is based on the parallel calling of the four SFB and S-RNase alleles that can be found in a diploid individual, and should therefore reduce allele misidentification and result in more robust genotyping.
Table
Table 1. Prunus salicina SFB and S-RNase allele sequences deposited at GenBank. Sequences used to build the concatenated artificial sequence are indicated.
Table 1. Prunus salicina SFB and S-RNase allele sequences deposited at GenBank. Sequences used to build the concatenated artificial sequence are indicated.
AlleleCultivar
(S-locus Genotype)		Nucleotide SequenceProtein SequenceCoding
Sequence (CDS)		Reference
IDLengthIDLength
SFBaBurmosa (SaSb)AM746961.1992 bpCAN90151.1331 aapartial[40]
SFBaOzarkpremier (SaSf)DQ849087.1977 bpABI15333.1325 aapartial[45]
SFBa  1Sordum (SaSb)AB252410.11131 bpBAF42763.1376 aacomplete[48]
Sa-RNase  1Sordum (SaSb)AB252411.11277 bpBAF42764.1226 aacomplete[48]
S1-RNaseRed-beaut (SaSb)AF433649.1576 bpAAP97311.195 aapartial[49]
SFBbBlack GoldenKJ396620.11131 bpAHX39360.1376 aacomplete[50]
SFBbGaixiandali (SbSd) DQ849088.1978 bpABI15334.1326 aapartial[45]
SFBbHamra BedriKJ396618.11131 bpAHX39358.1376 aacomplete[50]
SFBbSanta Rosa (ScSe)KJ396607.1 21131 bpAHX39347.1376 aacomplete[50]
SFBb  1Sordum (SaSb)AB252412.11131 bpBAF42765.1376 aacomplete[48]
Sb-RNase  1Sordum (SaSb) AB252413.12332 bpBAF42766.1221 aacomplete[48]
Sb-RNaseUnknownDQ646488.12406 bpABG36934.1159 aapartial[51]
SFBcAin TorkiaKJ396613.11128 bpAHX39353.1375 aacomplete[50]
SFBcBedri1 (SeSh)KJ396608.11128 bpAHX39348.1375 aacomplete[43,50,52]
SFBcMeiguili (ScSe)DQ849084.11128 bpABI15330.1375 aacomplete[47]
SFBc  1Santa Rosa (ScSe)AB280792.12124 bpBAF91847.1375 aacomplete[53]
Sc-RNaseOishiwasesumomo (ScSd)AB084144.11781 bpBAC20940.1177 aapartial[41]
Sc-RNaseSanta Rosa (ScSe)AB280791.12630 bpBAF91846.1230 aacomplete[53]
Sc-RNaseUnknownDQ646489.11920 bpABG36936.1172 aapartial[51]
S4-RNaseRoyal Zee (ScSe)AF432418.11283 bpAAP97308.195 aapartial[49]
SFBd  1Formosa (SbSd)AM746962.1992 bpCAN90152.1331 aapartial[40]
Sd-RNase  1Oishiwasesumomo (ScSd)AB084145.11976 bpBAC20941.1169 aapartial[41]
SFBeAouina Hamra Bedria (SeS?)KJ396612.11131 bpAHX39352.1376 aacomplete[50]
SFBeBlack Diamant/
Black Diamond (SeSh)		KJ396605.11128 bpAHX39345.1375 aacomplete[43,50]
SFBeCidre (SaSe)KJ396611.11128 bpAHX39351.1375 aacomplete[41,50]
SFBeMeiguili (ScSe)DQ849086.11036 bpABI15332.1345 aapartial[45]
SFBeSanta Rosa (ScSe)AB280794.11248 bpBAF91849.1375 aapartial[53]
SFBeStanleyKJ396606.1 31128 bpAHX39346.1375 aacomplete[50]
SFBe  1Unknown DQ646490.11749 bpABG36937.1373 aacomplete[51]
Se-RNase  1Santa Rosa (ScSe)AB280793.12622 bpBAF91848.1239 aacomplete[53]
S5-RNaseRoyal Zee (ScSe)AF433647.11553 bpAAP97309.195 aapartial[49]
SFBf  1Huangpili (S7Sf)DQ849089.1972 bpABI15335.1324 aapartial[45]
SFBfJanhaKJ396610.1996 bpAHX39350.1331 aacomplete[50]
SFBfMeski HamraKJ396614.11086 bpAHX39354.1361 aacomplete[50]
SFBfMeski KahlaKJ396617.11086 bpAHX39357.1361 aacomplete[50]
SFBfUnknownDQ989578.11081 bpABM54900.1360 aapartial[49]
SFBfZaghwenia (ScSf)KJ396615.1993 bpAHX39355.1330 aacomplete[50,52]
Sf-RNase  1White Plum (SfSg)AB084147.11554 bpBAC20943.1132 aapartial[41]
Sf-RNase mRNAHuangpili (S7Sf)DQ512911.1762 bpABF61820.1215 aapartial[54]
S6-RNaseWikson AF433648.11212 bpAAP97310.195 aapartial[49]
SFBgBonnie (SgSh)AM746963.1992 bpCAN90153.1331 aapartial[40]
SFBg  1UnknownDQ989579.11084 bpABM54901.1361 aapartial[51]
Sg-RNase  1Bonnie (SgSh)AM746950.11536 bpCAN90140.1169 aapartial[40]
Sg-RNaseHoney Rosa (SbSg)AB093131.11266 bpBAC75456.179 aapartial[41]
SFBhAin TasstouriaKJ396616.11125 bpAHX39356.1374 aacomplete[50]
SFBhBedri2KJ396619.11125 bpAHX39359.1374 aacomplete[50]
SFBhJabounia SafraKJ396609.11125 bpAHX39349.1374 aacomplete[50]
SFBh  1Nvgelei (ScSh)DQ849118.11131 bpABI15337.1376 aacomplete[45]
SFBhUnspecifiedDQ646491.11253 bpABG36938.1374 aapartial[51]
Sh-RNase  1Kelsey (SfSh)AB084148.11172 bpBAC20944.1175 aapartial[41]
Si-RNase  1Bakemonosumomo (SbSi)AB084149.1968 bpBAC20945.1170 aapartial[41]
Sj-RNase  1Tecumseh (SfSj)AB093132.12346 bpBAC75457.1173 aapartial[41]
SFBk  1Wickson (SfSk)DQ992485.11083 bpABM54902.1361 aapartial[44,51]
Sk-RNaseFriar (ShSk)DQ790372.1374 bpABH07013.163 aapartial[55]
Sk-RNase  1Starkgold (SgSk)AB093133.11035 bpBAC75458.1187 aapartial[41,56]
Sk-RNaseWickson (SfSk)EU113311.1632 bpABW86860.1149 aapartial[44,57]
S3-RNaseUnspecifiedAF432417.1467 bpAAP97307.195 aapartial[58]
Sl-RNase  1Combination (SgSl)AB0931341533 bpBAC75459.1188 aapartial[41]
Sm-RNase  1Botan (SaSm)AB093135.1930 bpBAC75460.1171 aapartial[41]
Sn-RNase  1Superior (SaSn)AB093136843 bpBAC75461.1121 aapartial[41,56]
St-RNase  1Karari (SbSt)AB573636.1525 bpBAJ13374.1175 aapartial[59]
SFB7  1Huangpili (S7Sf)DQ849085975 bpABI15331.1325 aapartial[45]
S7-RNase  1Huangpili (S7Sf)AY781290.1369 bpAAV34703.194 aapartial[60]
S7-RNase mRNAHuangpili (S7Sf)DQ512912.1761 bpABF61821.1217 aapartial[51]
S8-RNase  182-1-109AY902455.1450 bpAAW83158.194 aapartial[60]
S8-RNase mRNA82-1-109DQ512913.1750 bpABF61822.1213 aapartial[51]
S9-RNase  1Tai Ping Guo AY996051.1815 bpAAY15124.194 aapartial[60]
SFB10  1Mali (S10S14)DQ849090.11041 bpABI15336.1347 aapartial[45]
S10-RNase  1Hong Men LingDQ003310.1880 bpAAX98270.194 aapartial[60]
S11-RNase  1ZhushaliDQ512908.1670 bpABF61817.1119 aapartial[51]
S12-RNase  1Kongquedan (S12S13)DQ512909.1197 bpABF61818.165 aapartial[51]
S13-RNase  1Kongquedan (S12S13)DQ512910.1191 bpABF61819.163 aapartial[51]
S14-RNase  1Mali (S10S14)EF177345.1194 bpABM53475.137 aapartial[45]
S15-RNase  1Kuandiandali (S15S16)EF177346.1571 bpABM53476.139 aapartial[61]
S16-RNase  1Kuandiandali (S15S16)EF177347.1460 bpABM53477.139 aapartial[61]
S17-RNase  1Xinjiangsanhao (S17S19)EF177348.1261 bpABM53478.139 aapartial[61]
S18-RNase  1Danchengxingmei (S18S20)EF177349.1413 bpABM53479.138 aapartial[61]
S19-RNase  1Xinjiangsanhao (S17S19)EU113259.1690 bpABW88922.155 aapartial[62]
S20-RNase  1Danchengxingmei (S18S20)EU113260.11925 bpABW88923.155 aapartial[62]
S21-RNase  1Zhenzhuli (S21S22)EU113261.11601 bpABW88924.1116 aapartial[62]
S22-RNase  1Zhenzhuli (S21S22)EU113262.1853 bpABW88925.1117 aapartial[62]
S23-RNase  1Dameigui (S23S24)EU113263.1750 bpABW88926.154 aapartial[62]
S23-RNase  1GuangdonghongliFJ377732.1 4597 bpACI88749.164 aapartial[63]
S24-RNase  1Dameigui (S23S24)EU113264.11155 bpABW88927.1117 aapartial[62]
S25-RNase  1UnspecifiedEU113265.1920 bpABW88928.1116 aapartial[62]
S26-RNase  1UnspecifiedEU113266.11492 bpABW88929.1117 aapartial[62]
S27-RNase  1UnspecifiedEU113267.1657 bpABW88930.154 aapartial[62]
S32-RNase  1GualiGU574195.1 5898 bpADD20972.1100 aapartial[64]
S41-RNase  1QingzhouliGU968758.1719 bpADF42685.1121 aapartial[65]
1 Used for genotyping. 2 BLASTN sequence is 100% identical to SFBb but the cultivar is likely wrong; ‘Santa Rosa’ was genotyped elsewhere and here as ScSe. 3 BLASTN sequence is 100% identical to SFBe but the cultivar name is reported in the literature as Prunus domestica. 4 A duplicate case: same name given to different alleles. 5 Highly similar to S10-Rnases at the nucleotide level.

2.2. A Powerful S-Allele Calling Procedure

Raw individual shotgun sequences from 88 Japanese plum cultivar samples (83 genomic, four exon targeted-capture, and one transcriptome libraries (Table S1 and references therein) were aligned to the synthetic sequence using Bowtie2, following the end-to-end default sensitive mode (Figure 2). To increase the stringency of the alignment for higher S-allele specificity, reads with a mapping quality lower than 20, reads shorter than 20bp in length, and reads with more than two mismatches with the reference sequence were filtered out.
To accurately genotype the S-locus of all samples, we followed a strategy that combined a visual inspection of the alignment with the implementation of a quantitative analysis. First, we visually inspected all alignments and established a preliminary S-genotype based on gene coverage. Second, we used different tools that report per-base depth and per-gene breadth of the coverage. The per-base depth was plotted as raw data (area graph) and fitted as probability density distribution. The per-gene breadth of the coverage was plotted as bar plots of the number of reads per gene and the proportion of the gene covered by at least one read base (see Figure S1 for all genotypes analyzed). Whereas each plot showed a bit of information about the likely S-genotype, careful examination of all figures followed by quantitative analysis was necessary to contrast the preliminary visual inspection and to call a positive S-allele. The probability density plot was especially informative about the actual likelihood to observe a given peak. The size of the area below the curve that a given gene covers is in fact and approximation of the probability of observing that peak considering that the whole area under the curve should sum up to one (see Table S2 for all genotypes analyzed). The proportion of the gene covered by at least one read base was especially relevant for discarding those unspecific alignments with high coverage but limited breadth, which are characteristic of intronic and intergenic regions with repetitive and transposable elements that are frequently found throughout the genome. We considered adding the rest of the genome to the synthetic sequence, but this approach significantly increased the duration of the alignment procedure with similar results to our three-steps validation procedure. Positive S-allele calls were based on the highest gene breadth of coverage and high likelihood in congruence with the visual inspection of the alignment.
Of the 88 shotgun samples, we were able to unambiguously genotype two S-alleles in 74 cultivars, only one S-allele in 11 cultivars, and none of the S-alleles in four cultivars (Table 2). The low depth and breadth of coverage in the S-locus regions extracted from the reference genomes (Sb, Ssany and Szgl, see Figure S1) for most cultivars not carrying those three alleles illustrate the high degree of polymorphism at the S-locus, and how reference genomes cannot be used for genotyping this locus. Our methodology clearly improved the approach used in M. domestica based solely on per base depth of coverage, which is highly sensitive to non-specific coverage of intronic regions [33]. Although this limitation was successfully circumvented in the previous work by setting the per-base depth of coverage only for exonic regions, this might lead to misidentification of alleles in species with lower polymorphism than apple in transcribed regions. Our approach, which includes intronic regions while solving the problem of non-specific alignment, allowed us to reduce these identifications errors.
Ijms 24 03932 g001 550
Figure 1. Full S-loci sequences retrieved from published reference genomes and synthetic reference sequence structure. (A,B) Pairwise trans-specific identities of the 6 retrieved S-RNase and SFB alleles. (C) Comparison of the amino-acids predicted from the coding regions of the three S-RNases. Signal peptide, conserved regions, and hypervariable regions (as defined by [66]) are highlighted. (D) Comparison of the amino acids predicted from the coding regions of the three SFBs. F-box, variable, and hypervariable regions (as defined by [47]) are highlighted. The amino acids of the S-RNases and SFBs were aligned using Muscle5 [67]. (E) Structure of the synthetic reference sequence used for aligning resequenced individuals. Genes (green), CDSs (blue), upstream, intergenic, and downstream regions (orange), and 1000 bp spacers (grey) and their orientation are illustrated. Sizes of the boxes are proportional to sequence length. Accession numbers of the sequences of each allele are given in Table 1. Full S-loci regions for Sb, Ssany, and Szgl, including 1000 bp upstream and downstream the locus and the intergenic regions, were retrieved from the threes published reference genomes of ‘Zhongly’ (Chr6:14567566:14573591), ‘Sanyueli’ (Chr6:33823425:33831243) and ZGL (Chr2:2766997:2772560), respectively.
Figure 1. Full S-loci sequences retrieved from published reference genomes and synthetic reference sequence structure. (A,B) Pairwise trans-specific identities of the 6 retrieved S-RNase and SFB alleles. (C) Comparison of the amino-acids predicted from the coding regions of the three S-RNases. Signal peptide, conserved regions, and hypervariable regions (as defined by [66]) are highlighted. (D) Comparison of the amino acids predicted from the coding regions of the three SFBs. F-box, variable, and hypervariable regions (as defined by [47]) are highlighted. The amino acids of the S-RNases and SFBs were aligned using Muscle5 [67]. (E) Structure of the synthetic reference sequence used for aligning resequenced individuals. Genes (green), CDSs (blue), upstream, intergenic, and downstream regions (orange), and 1000 bp spacers (grey) and their orientation are illustrated. Sizes of the boxes are proportional to sequence length. Accession numbers of the sequences of each allele are given in Table 1. Full S-loci regions for Sb, Ssany, and Szgl, including 1000 bp upstream and downstream the locus and the intergenic regions, were retrieved from the threes published reference genomes of ‘Zhongly’ (Chr6:14567566:14573591), ‘Sanyueli’ (Chr6:33823425:33831243) and ZGL (Chr2:2766997:2772560), respectively.
Ijms 24 03932 g001
All individuals with a putative S10 or S32 allele did not show an exclusive alignment to any of the two loci. Although this result was expected due to the highly similarity between both alleles at the nucleotide level, we finally decided to assign a mixed genotype S10/32 to all individuals with this genotype. De novo assembly or full sequencing of both S-alleles in their respective cultivars will clarify whether they represent the same or different alleles. All unidentified S-alleles were assumed to be novel unsequenced alleles. Furthermore, since the S-RNase and SFB genes are known to be expressed in flower tissues during the anthesis period [7,8,9,10], we demonstrated that the S-genotype of the cultivar ‘Sanyueli’, as established from genomic reads (S10/32Ssany) could also be inferred from flower transcriptome reads (Figure S1, Table S2). The lower likelihood of positive calls and shorter breadths of gene coverage can be explained by the absence of coverage within the intronic regions. In addition, we extended our analysis to test whether shotgun sequences from exome-targeted-capture libraries could be useful for S-genotyping. We chose four libraries from cultivars of previously known genotypes [35,41,42] and aligned them to the synthetic sequence. Although we successfully confirmed the S-genotype of all 4 cultivars (Table 2), lower values of likelihood and breadth of coverage should also reflect the absence of intronic coverage among others (Figure S1, Table S2). In summary, our approach represents a relevant addition to the methodologies previously described in M. domestica [33] and Arabidopsis halleri [34], as it combined the strategy of using the entire S-locus, the high stringency of the alignment of shotgun sequences from different types of libraries, and the quantitative analysis based on the determination of both the probability of positive calls and the percentage of gene coverage.

2.3. S-Allele Genotypes, Methodology Validation, and Incompatibility Groups

The approach developed here allowed for the identification of the S-genotype in 88 Japanese plum cultivars, 74 of them reported for the first time (Table 2). At least two S-alleles were identified in 74 cultivars, which, according to their S-allele composition, were assigned to 22 incompatibility groups, including nine new incompatibility groups reported for the first time (XXVII-XXXV, Table 2). In addition, 18 cultivars were included in group 0, since no other cultivars with the same S-genotype have been reported so far and, therefore, they can be considered as potential universal pollinators [14]. This information could be highly valuable for growers and breeders, since the knowledge of incompatibility relationships is intended to help fruit growers to select compatible pollinizers and breeders to choose parental genotypes.
To validate our methodology, we compared our results with already known S-genotypes of those cultivars for which S-locus sequences are available in the GenBank (Table 1 and Table 2). The results confirmed the S-genotype of 13 cultivars: ‘Sordum’ [35], ‘Oishiwase’ and ‘Taiyo’ [41], ‘Angeleno’ [43], ‘Fortune’ and ‘Friar’ [68], ‘Jinshali’ and ‘Qiuji’ [45], ‘Pingguoli’ and ‘Zuili’ [46], and the two self-compatible cultivars ‘Honeyrosa’ [42] and ‘Santa Rosa’ [68]. However, the S-genotype of ‘Beauty’ (SbSd) and ‘Friar’ (SbSh) differed from previously reported genotypes. ‘Beauty’ was previously genotyped as ScSe [41] and ‘Friar’ as ShSk [55]. Based on a previous report on the S-allele inheritance in ‘Friar’ [68], its most likely S-genotype is SbSh. On the other hand, further research using additional plant material from ‘Beauty’ would be needed to clarify its S-genotype. Although further validation of the remaining cultivars -as more S-locus sequences become available- and clarification of those conflicting results would provide definitive support for the validity of this methodology, confirmation of the S-genotype in 14 cultivars is a strong indication that our methodology is reliable.
In 11 cultivars, only one S-allele could be identified and, therefore, they were not assigned to any incompatibility group and were considered unclassified (Table 2). In four cultivars, no S-allele could be identified, suggesting that at least two novel S-alleles were detected. S-allele de novo assembly complemented with PCR analysis based on the amplification of S-RNase and SFB genes in these 15 unclassified cultivars would be needed to identify the new S-alleles and clarify their S-genotype. Three pairs of cultivars have a highly similar names and the same S-genotype (‘Changlixiangjiaoli’ and ‘Changlixiangbian’, ‘Guofeng2’ and ‘N2 Guofeng’, ‘Cuihongli’ and ‘Cuihong’) so they could be considered synonyms. On the other hand, ‘Li he’ and ‘Lihe’, with highly similar names, showed different S-genotypes, so they are distinct genotypes and could be considered homonyms. The presence of three and four alleles in the cultivars ‘Guiyang’ and ‘Xingyikongxinli’, respectively, means that either they have undergone some S-locus duplication events, or, more likely, they are triploid and tetraploid individuals, respectively.
No information on self-(in)compatibility is available for most of the cultivars tested (all except two self-compatible cultivars ’Honeyrosa’ [42] and ‘Santa Rosa’ [68]). Controlled pollinations and subsequent microscopic analysis of pollen tube growth in self-pollinated flowers would be necessary, since the identification of the S-alleles of a cultivar does not accurately predict its self-(in)compatibility [14].
Japanese plum, which include pure P. salicina cultivars and hybrids of P. salicina with other Prunus species [14,69], is an ideal system to undertake such S-genotyping for several scientific and practical reasons. P. salicina is considered to have originated in south-western China from where it spread to other Chinese territories, then to Japan, from Japan to California, and from there to other parts of the world such as Europe, New Zealand, South America and South Africa [70,71,72]. The highest diversity of S-alleles was observed in cultivars putatively originating from the southern provinces of China (Table 2, Table S1), further supporting the putative origin of the species in this region. Furthermore, no cultivar of putative Chinese origin had the Sa and Sl alleles, suggesting a probable introgression of these alleles from other Prunus species through hybridization. Interestingly, the cultivars ‘Changlixiangjiaoli’ and ‘Changlixiangbian’ from Hebei province in China, genotyped here as SeSk (Table 2 and Table S1), were claimed to belong to P. simonii in a previous study aimed at characterizing the origin and diversity of this species [72]. Whether this S-genotype reflects trans-specific polymorphism or an introgression of one or both alleles in Japanese plum hybrids is unknown. Undoubtedly, further studies on Chinese germplasm aimed at distinguishing pure P. salicina cultivars from interspecific hybrids will reveal more introgressed S-alleles.
Ijms 24 03932 g002 550
Figure 2. Bioinformatics pipeline for S-genotyping by high throughput sequencing. The visual inspection illustrates cultivars ‘Santa Rosa’ and ‘Beauty; the per base depth and per gene depth of the coverage illustrate the cultivar ‘Santa Rosa’ (see Figure S1 for all genotypes analyzed).
Figure 2. Bioinformatics pipeline for S-genotyping by high throughput sequencing. The visual inspection illustrates cultivars ‘Santa Rosa’ and ‘Beauty; the per base depth and per gene depth of the coverage illustrate the cultivar ‘Santa Rosa’ (see Figure S1 for all genotypes analyzed).
Ijms 24 03932 g002
Table
Table 2. Incompatibility groups (I. G.), S-genotype and library type (E: exon-targeted capture, G: genomic, T: transcriptome) of 88 Japanese plum cultivars.
Table 2. Incompatibility groups (I. G.), S-genotype and library type (E: exon-targeted capture, G: genomic, T: transcriptome) of 88 Japanese plum cultivars.
I. G.S-allelesCultivarOrigin 1LibraryReference
ISaSbSordumJapanE[35]
IISbScFortuneCalifornia, USAG[68]
ZaoshengyueguangJapanG
DashizhongshengJapanG
TaiyoJapanE[41]
IIISbSfZhenyuanqiyuexiangGuizhou ChinaG
IVSbShFriarCalifornia, USAG[68]
Friar 2California, USAG[55] (ShSk)
QiujiJapanG[45]
Guofeng7Liaoning, China G
C20SpainG
Zhongly nº 6Northern China G
VSbSiFujianfurongliFujian, ChinaG
VISfShSatsumaCalifornia, USAG
VIIScShFulihongShaanxi, ChinaG
AngelenoCalifornia, USAG[43]
IXSfSgJinshaliYunnan, ChinaG[45]
CehengjixueliGuizhou, ChinaG
XIVSaScC46SpainG
XIXSbSdAbazhoumeiguiliSichuan, ChinaG
Beauty 2California, USAG[41] (ScSe)
XXISeSkWuxiangliHebei, ChinaG
ChanglixiangjiaoliHebei, ChinaG
ChanglixiangbianHebei, ChinaG
XXIIIScSdDashizaoshengJapanG
JingshangJapanG
OishiwaseJapanG[41]
OishiwaseJapanE[41]
XXIVSeS11JinxiqiuliLiaoning, China G
QiyuexiangLiaoning, China G
XXVS15S16PingguoliLiaoning, China G[46]
XXVII 3SbS7HeyuansanhualiGuangdong, ChinaG
QingpihongxinGuangdong, ChinaG
ZaohuangliAnhui, ChinaG
JinganzhushaliJiangxi, ChinaG
JaiqingziJiangsu, ChinaG
XXVII 3SbS8HuahongliYunnan, ChinaG
BraceCalifornia, USAG
XXIX 3SeS10/32ShuiliGuangxi, ChinaG
DahuiliHenan, ChinaG
XiangjiaoliLiaoning, China G
QiuxiangLiaoning, China G
XXX 3SgSkYuanshuaiJapanG
XXXI 3ShS10/32Guofeng2Liaoning, China G
N2 GuofengLiaoning, China G
XXXII 3SiS9ChangkuangsanhualiGuangdong, ChinaG
Sanhua plumGuangdong, ChinaG
XXXIII 3S10/32S11TiankanmaliSichuan, ChinaG
Li HeChinaG
Wushan plumChongging, ChinaG
Wanshuang plumChongging, ChinaG
XXXIV 3S11S20YuhuangliShaanxi, ChinaG
HuangliShandong, ChinaG
LiheHebei, ChinaG
XXXV 3SbScSfGuiyangJapanG
Group 0SaSlHongshouJapanG
SbS11YanzhiliFujian, ChinaG
SdSgFenghualiZhejiang, ChinaG
SeS12HuangguliZhejiang, ChinaG
SeS20DaqingkeShandong, ChinaG
ShS15Guofeng17Liaoning, China G
ShS7ZuiliZhejiang, ChinaG[46]
SkS8BullbankCalifornia, USAG
S7S20PingdingxiangShandong, ChinaG
S8S9WanshuhuanaiFujian, ChinaG
S10/32S15ZhengzhouzaoliHenan, ChinaG
S10/32SsanySanyueli Guangdong, ChinaG
Sanyueli Guangdong, ChinaT
S10/32SzglCuihongliSichuan, ChinaG
Cuihong plumSichuan, ChinaG
S11S16SaozouliGuangxi, ChinaG
SeSkS10/32S11XingyikongxinliGuizhou, ChinaG
Cultivars in which only one S allele could be identified
SbCuipinwannaiFujian, ChinaG
ScGuoliLiaoning, China G
SeTianmumiliZhejiang, ChinaG
SeHongxinliAnhui, ChinaG
ShGuomeiLiaoning, China G
S10/32Yinhong plumSichuan, ChinaG
S10/32Fengtang plumGuizhou, ChinaG
S11TongkeliGuangxi, ChinaG
S19Non-specifiedTibet, ChinaG
S21DamiliGuangdong, ChinaG
S22ZhushaliYunnan, ChinaG
Cultivars in which no S allele could be identified
KuaishiliGuangdong, ChinaG
LushanliJiangxi, ChinaG
Ni Ma Qu Ji LiChinaG
WuyuecuiSichuan, ChinaG
Self-compatible cultivars
SCSbSgHoneyrosaJapanE[42]
SCScSeSanta RosaCalifornia, USAG[68]
1 Presumed origin, and in its absence the sampling location is given. 2 Cultivars described with different S-genotype. 3 New Incompatibility Groups first reported herein.

3. Materials and Methods

3.1. Building the Reference Synthetic S-Loci Sequence

Japanese plum sequences for S-RNase and SFB were retrieved from NCBI based on the literature and on direct database search. In addition, we conducted a BLASTN search on the three available Japanese plum reference genomes of Chinese origin (cultivars ‘Sanyueli’ assembly V2.0 at rosaceae.org [73], ‘Zhongli’ assembly V1.0 at rosaceae.org [74] and ‘ZGL’ assembly GCA_020226455.1 at NCBI [75]) and recovered three full length S-locus regions. The synthetic sequence was built by concatenating individual sequences in a string-like structure. Depending on the sequences available per locus, it contained sequences from three full S-locus regions, individual S-locus sequences in tail-to-tail orientation (SFB is transcribed on the forward strand and S-RNase on the reverse), and individual S-RNase sequences for loci with unsequenced SFB alleles. BLAST searches were performed at NCBI [76] or locally on the online platform useglaxy.eu using the Galaxy BLAST+ blastn v 2.10.1+galaxy0 tool [77,78]. Sequence editing, alignment, and concatenation were performed with Geneious prime 202.2.2.

3.2. Alignment to Reference Sequence

Individual shotgun sequencing reads from Japanese plum genomic, transcriptomic and exon targeted-capture libraries were aligned to the synthetic sequence. Alignment was performed on the online platform usegalaxy.eu [78] following the pipeline outlined in Figure 2. The raw reads were first subjected to quality control (FastQC, [79]), then trimmed (Trimmomatic [80]), and finally aligned to the reference sequence using Bowtie2 [81] with the end-to-end default sensitive mode. Samtools View [82] was used to filter out reads with mapping quality below 20, and Bam Filter [83] was used to filter out reads shorter than 20 bp or having more than 2 mismatches with the reference sequence.

3.3. S-Allele Typing Methodology

To accurately type S-allele genotypes, we used a combination of visual inspection of the alignment using IGV [84] and quantitative analyses using SAMtools [82] and BEDtools coverage [85], followed by graphical analyses using the R package ggplot2 [86]. The likelihood of a positive call was inferred by integrating the probability density distribution of read depth at the gene level, and by determining the proportion of the gene covered by reads. Final S-allele call was based on full congruence between the three approaches. A positive call should have a high probability, the highest proportion of gene coverage, and should be congruent with visual inspection of the alignment.

4. Conclusions

Building on previous S-genotyping procedures [33,34], we developed a modified and robust methodology to S-genotype Japanese plum cultivars using whole genome shotgun sequences. The technical advance of the method is based on the use of the entire S-locus allelic series of the species to construct a synthetic reference sequence on which resequenced genotypes are aligned. S-genotyping is based on parallel calling of the four alleles that form the S-locus, combined with a simple but powerful quantitative methodology. Further work is underway to de-novo-assemble Japanese plum S-alleles and obtain longer sequence for many alleles for which only short sequences are available or discover novel SFB and S-RNases alleles.
S-genotyping cultivars of Chinese origin would ultimately provide valuable information on the original allelic series that can be found in pure P. salicina cultivars and would allow us to identify all those S-alleles, and their respective species, that are the result of the frequent hybridization events in modern Japanese plum hybrid cultivars. It would also enable us to build a comprehensive database of inter-compatibility relationships for all available germplasm, which would be of great value for breeders and growers. The approach developed here could be applied to any other plant species—with a GSI system—provided that several S-loci have been sequenced. A study is underway to apply this S-genotyping-by-sequencing procedure to several economically important Prunus fruit tree species such as Prunus armeniaca, Prunus avium and Prunus dulcis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24043932/s1, References [87,88,89,90,91,92,93,94,95] are cited in the Supplementary Materials.

Author Contributions

A.H., M.E.G., J.G. and J.R. conceived the study; A.H. performed the analyses a prepared the final data and figures. A.H., M.E.G., J.G. and J.R. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by MCIN/AEI and by “ERDF A way to make Europe MCIN/AEI/10.13039/501100011033” (project PID2020-115473RR-I00), Junta de Extremadura—Fondo Europeo de Desarrollo Regional, FEDER (Grupo de Investigación AGA021; RBFruit), and Gobierno de Aragón—European Social Fund, European Union (Grupo Consolidado A12_17R).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. de Nettancourt, D. Incompatibility in Angiosperms; Springer: Berlin/Heidelberg, Germany, 1977. [Google Scholar]
  2. Steinbachs, J.E.; Holsinger, K.E. S-RNase-medited Gametophytic Self-Incompatibility is Ancestral in Eudicots. Mol. Biol. Evol. 2002, 19, 825–829. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. de Nettancourt, D. Incompatibility and Incongruity in Wild and Cultivated Plants; Springer: Berlin/Heidelberg, Germany, 2001. [Google Scholar]
  4. Uyenoyama, M.K. Evolutionary Dynamics of Self-Incompatibility Alleles in Brassica. Genetics 2000, 156, 351–359. [Google Scholar] [CrossRef] [PubMed]
  5. Wright, S. The distribution of self-sterility alleles in populations. Genetics 1939, 24, 538. [Google Scholar] [CrossRef]
  6. Wu, J.; Gu, C.; Khan, M.A.; Gao, Y.; Wang, C.; Korban, S.S.; Zhang, S. Molecular Determinants and Mechanisms of Gametophytic Self-Incompatibility in Fruit Trees of Rosaceae. CRC Crit. Rev. Plant Sci. 2013, 32, 53–68. [Google Scholar] [CrossRef]
  7. Ushijima, K.; Sassa, H.; Dandekar, A.M.; Gradziel, T.M.; 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. Plant Cell 2003, 15, 771–781. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Tao, R.; Yamane, H.; Sassa, H.; Mori, H.; Gradziel, T.M.; Dandekar, A.M.; Sugiura, A. Identification of Stylar RNases Associated with Gametophytic Self-Incompatibility in Almond (Prunus dulcis). Plant Cell Physiol. 1997, 38, 304–311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. 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. J. Am. Soc. Hortic. Sci. 1999, 124, 224–233. [Google Scholar] [CrossRef] [Green Version]
  10. McCubbin, A.G.; Kao, T. Molecular Recognition and Response in Pollen and Pistil Interactions. Annu. Rev. Cell Dev. Biol. 2000, 16, 333–364. [Google Scholar] [CrossRef]
  11. Castric, V.; Vekemans, X. Invited Review: Plant self-incompatibility in natural populations: A critical assessment of recent theoretical and empirical advances. Mol. Ecol. 2004, 13, 2873–2889. [Google Scholar] [CrossRef]
  12. Halász, J.; Molnár, A.B.; Ilhan, G.; Ercisli, S.; Hegedüs, A. Identification and Molecular Analysis of Putative Self-Incompatibility Ribonuclease Alleles in an Extreme Polyploid Species, Prunus laurocerasus L. Front. Plant Sci. 2021, 12, 715414. [Google Scholar] [CrossRef]
  13. Herrera, S.; Lora, J.; Hormaza, J.I.; Rodrigo, J. Pollination Management in Stone Fruit Crops. In Production Technology of Stone Fruits; Mir, M.M., Iqbal, U., Mir, S.A., Eds.; Springer: Singapore, 2021; pp. 75–102. [Google Scholar]
  14. Guerra, M.E.; Rodrigo, J. Japanese plum pollination: A review. Sci. Hortic. 2015, 197, 674–686. [Google Scholar] [CrossRef]
  15. Guerra, M.E.; Guerrero, B.I.; Casadomet, C.; Rodrigo, J. Self-(in)compatibility, S-RNase allele identification, and selection of pollinizers in new Japanese plum-type cultivars. Sci. Hortic. 2020, 261, 109022. [Google Scholar] [CrossRef]
  16. Busch, J.W.; Schoen, D.J. The evolution of self-incompatibility when mates are limiting. Trends Plant Sci. 2008, 1s, 128–136. [Google Scholar] [CrossRef] [PubMed]
  17. Pickup, M.; Brandvain, Y.; Fraïsse, C.; Yakimowski, S.; Barton, N.H.; Dixit, T.; Lexer, C.; Cereghetti, E.; Field, D.L. Mating system variation in hybrid zones: Facilitation, barriers and asymmetries to gene flow. New Phytol. 2019, 224, 1035–1047. [Google Scholar] [CrossRef] [PubMed]
  18. Williams, R.R. Techniques used in fruit-set experiments. In Towards Regulated Cropping; Williams, R.R., Wilson, D., Eds.; Grower Books: London, UK, 1970; pp. 57–61. [Google Scholar]
  19. Herrero, M.; Dickinson, H.G. Pollen tube growth following compatible and incompatible intraspecific pollinations in Petunia hybrida. Planta 1980, 148, 217–221. [Google Scholar] [CrossRef] [PubMed]
  20. Rodrigo, J.; Herrero, M. Evaluation of pollination as the cause of erratic fruit set in apricot ‘Moniqui’. J. Hortic. Sci. 1996, 71, 801–805. [Google Scholar] [CrossRef]
  21. Sanzol, J.; Herrero, M. Identification of self-incompatibility alleles in pear cultivars (Pyrus communis L.). Euphytica 2002, 128, 325–331. [Google Scholar] [CrossRef]
  22. Guerra, M.E.; Wünsch, A.; López-Corrales, M.; Rodrigo, J. Lack of Fruit Set Caused by Ovule Degeneration in Japanese plum. J. Am. Soc. Hortic. Sci. 2011, 136, 375–381. [Google Scholar] [CrossRef] [Green Version]
  23. Bošković, R.; Tobutt, K.R. Correlation of Stylar Ribonuclease Zymograms with Incompatibility Alleles in Sweet Cherry. Euphytica 1996, 90, 245–250. [Google Scholar] [CrossRef]
  24. Sassa, H.; Nishio, T.; Kowyama, Y.; Hirano, H.; Koba, T.; Ikehashi, H. Self-incompatibility (S) alleles of the Rosaceae encode members of a distinct class of the T2/S ribonuclease superfamily. Mol. Gen. Genet. 1996, 250, 547–557. [Google Scholar]
  25. Halász, J.; Hegedüs, A. A critical evaluation of methods used for S-genotyping: From trees to DNA level. Int. J. Hortic. Sci. 2006, 12, 19–29. [Google Scholar] [CrossRef]
  26. Yamane, H.; Tao, R. Molecular Basis of Self-(in)compatibility and Current Status of S-genotyping in Rosaceous Fruit Trees. J. Jpn. Soc. Hortic. Sci. 2009, 78, 137–157. [Google Scholar] [CrossRef] [Green Version]
  27. Sutherland, B.G.; Robbins, T.P.; Tobutt, K.R. Primers amplifying a range of Prunus S-alleles. Plant Breed. 2004, 123, 582–584. [Google Scholar] [CrossRef]
  28. Delorme, V.; Gaude, T.; Heizmann, P.; Dumas, C. Use of Inmunochemical and SSCP analyses to test homozygosity at the S Locus of Brassica oleracea genotypes. Mol. Breed. 1995, 1, 237–244. [Google Scholar] [CrossRef]
  29. Jorgensen, R.A. Plant Genetics and Genomics: Crops and Models; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  30. Mable, B.K.; Hagmann, J.; Kim, S.T.; Adam, A.; Kilbride, E.; Weigel, D.; Stift, M. What causes mating system shifts in plants? Arabidopsis lyrata as a case study. Heredity 2017, 118, 52–63. [Google Scholar]
  31. Mable, B.K.; Brysting, A.K.; Jørgensen, M.H.; Carbonell, A.K.Z.; Kiefer, C.; Ruiz-Duarte, P.; Lagesen, K.; Koch, M.A. Adding Complexity to Complexity: Gene Family Evolution in Polyploids. Front. Ecol. Evol. 2018, 6, 114. [Google Scholar] [CrossRef] [Green Version]
  32. Tsuchimatsu, T.; Goubet, P.M.; Gallina, S.; Holl, A.C.; Fobis-Loisy, I.; Bergès, H.; Marande, W.; Prat, E.; Meng, D.; Long, Q.; et al. Patterns of Polymorphism at the Self-Incompatibility Locus in 1083 Arabidopsis thaliana Genomes. Mol. Biol. Evol. 2017, 34, 1878–1889. [Google Scholar] [CrossRef] [Green Version]
  33. De Franceschi, P.; Bianco, L.; Cestaro, A.; Dondini, L.; Velasco, R. Characterization of 25 full-length S-RNase alleles, including flanking regions, from a pool of resequenced apple cultivars. Plant Mol. Biol. 2018, 97, 279–296. [Google Scholar] [CrossRef] [Green Version]
  34. Genete, M.; Castric, V.; Vekemans, X. Genotyping and De Novo Discovery of Allelic Variants at the Brassicaceae Self-Incompatibility Locus from Short-Read Sequencing Data. Mol. Biol. Evol. 2020, 37, 1193–1201. [Google Scholar] [CrossRef]
  35. Yamane, H.; Tao, R.; Sigiura, A. Identification an cDNA Cloning for S-RNases in Self-incompatible Japanese Plum (Prunus salicina Lindl. cv. Sordum). Plant Biotechnol. 1999, 16, 389–396. [Google Scholar] [CrossRef] [Green Version]
  36. Lawrence, M.J. Population Genetics of the Homomorphic Self-incompatibility Polymorphisms in Flowering Plants. Ann. Bot. 2000, 85, 221–226. [Google Scholar] [CrossRef] [Green Version]
  37. Igic, B.; Bohs, L.; Kohn, J.R. Historical inferences from the self-incompatibility locus. New Phytol. 2004, 161, 97–105. [Google Scholar] [CrossRef]
  38. Richman, A.D.; Kohn, J.R. Learning from rejection: The evolutionary biology of single-locus incompatibility. Trends Ecol. Evol. 1996, 11, 497–502. [Google Scholar] [CrossRef] [PubMed]
  39. Ioerger, T.R.; Clark, A.G.; Kao, T. Polymorphism at the self-incompatibility locus in Solanaceae predates speciation. Proc. Natl. Acad. Sci. USA 1990, 87, 9732–9735. [Google Scholar] [CrossRef] [Green Version]
  40. Sutherland, B.G.; Tobutt, K.R.; Robbins, T.P. Trans-specific S-RNase and SFB alleles in Prunus self-incompatibility haplotypes. Mol. Genet. Genom. 2008, 279, 95–106. [Google Scholar] [CrossRef]
  41. Beppu, K.; Yamane, H.; Yaegaki, H.; Yamaguchi, M.; Kataoka, I.; Ryutaro, T. Diversity of S-RNase genes and S-haplotypes in Japanese plum (Prunus salicina Lindl.). J. Hortic. Sci. Biotechnol. 2002, 77, 658–664. [Google Scholar] [CrossRef]
  42. Beppu, K.; Takemoto, Y.; Yamane, H.; Yaegaki, H.; Yamaguchi, M.; Kataoka, I.; Tao, R. Determination of S-haplotypes of Japansese plum (Prunus salicina Lindl.) cultivars by PCR and cross-pollination test. J. Hortic. Sci. Biotechnol. 2003, 78, 315–318. [Google Scholar] [CrossRef]
  43. Sapir, G.; Stern, R.A.; Eisikowitch, D.; Goldway, M. Cloning of four new Japanese plum S-alleles and determination of the compatibility between cultivars by PCR analysis. J. Hortic. Sci. Biotechnol. 2004, 79, 223–227. [Google Scholar] [CrossRef]
  44. Sapir, G.; Stern, R.A.; Goldway, M.; Shafir, S. SFBs of Japanese Plum (Prunus salicina): Cloning Seven Alleles and Determining Their Linkage to the S-RNase Gene. HortScience 2007, 42, 1509–1512. [Google Scholar] [CrossRef] [Green Version]
  45. Zhang, S.L.; Huang, S.X.; Kitashiba, H.; Nishio, T. Identification of S-haplotype-specific F-box gene in Japanese plum (Prunus salicina Lindl.). Sex. Plant Reprod. 2007, 20, 1–8. [Google Scholar] [CrossRef]
  46. Zhang, S.J.; Huang, S.X.; Heng, W.; Wu, H.Q.; Wu, J.; Zhang, S.L. Identification of S-genotypes in 17 Chinese cultivars of Japanese plum (Prunus salicina Lindl.) and molecular characterisation of 13 novel S-alleles. J. Hortic. Sci. Biotechnol. 2008, 83, 635–640. [Google Scholar] [CrossRef]
  47. Ikeda, K.; Igic, B.; Ushijima, K.; Yamane, H.; Hauck, N.R.; Nakano, R.; Sassa, H.; Iezzoni, A.F.; Kohn, J.; Tao, R. Primary structural features of the S haplotype-specific F-box protein, SFB, in Prunus. Sex. Plant Reprod. 2004, 15, 235–243. [Google Scholar] [CrossRef]
  48. Tao, R.; Watari, A.; Hanada, T.; Habu, T.; Yaegaki, H.; Yamaguchi, M.; Yamane, H. Self-compatible peach (Prunus persica) has mutant versions of the S-haplotypes found in self-incompatible Prunus species. Plant Mol. Biol. 2007, 63, 109–123. [Google Scholar] [CrossRef] [PubMed]
  49. Sapir, G.; Goldway, M. Prunus Salicina S-RNases. Nucleotide Sequences AF432418.1, AF433647.1, AF433648.1, AF433649.1. Available online: https://www.ncbi.nlm.nih.gov/nuccore (accessed on 28 December 2022).
  50. Abdallah, D.; Baraket, G.; Ben Tamarzizt, H.; Ben Mustapha, S.; Salhi Hannachi, A. Identification, Evolutionary Patterns and Intragenic Recombination of the Gametophytic Self Incompatibility Pollen Gene (SFB) in Tunisian Prunus Species (Rosaceae). Plant Mol. Biol. Report. 2016, 34, 339–352. [Google Scholar]
  51. Sapir, G.; Goldway, M. Identification of F-box Gene in Prunus Salicina. Nucleotide Sequences AF432417.1, AF433649.1, DQ646488.1, DQ646489.1, DQ646490.1, DQ646491.1, DQ989578.1, DQ989579.1, DQ992485.1. Available online: https://www.ncbi.nlm.nih.gov/nuccore (accessed on 28 December 2022).
  52. Abdallah, D.; Baraket, G.; Perez, V.; Ben Mustapha, S.; Salhi-Hannachi, A.; Hormaza, J.I. Analysis of self-incompatibility and genetic diversity in diploid and hexaploid plum genotypes. Front. Plant Sci. 2019, 10, 896. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Watari, A.; Hanada, T.; Yamane, H.; Esumi, T.; Tao, R.; Yaegaki, H.; Yamaguchi, M.; Beppu, K.; Kataoka, I. A Low Transcriptional Level of se-RNase in the Se-haplotype Confers Self-compatibility in Japanese Plum. J. Am. Soc. Hortic. Sci. 2007, 132, 396–406. [Google Scholar] [CrossRef] [Green Version]
  54. Zhang, S.L.; Huang, S.X.; Zhang, S.J.; Wu, J.; Wu, H.Q.; Heng, W.; Zhang, Y.Y. Prunus Salicina S-RNase Genes. Nucleotide Sequences DQ512908.1, DQ512909.1, DQ512910.1, DQ512911.1, DQ512912.1, DQ512913.1. Available online: https://www.ncbi.nlm.nih.gov/nuccore (accessed on 28 December 2022).
  55. Halász, J.; Hegedus, A.; Szabó, Z.; Nyéki, J.; Pedryc, A. DNA-based S-genotyping of Japanese Plum and Pluot Cultivars to Clarify Incompatibility Relationships. HortScience 2007, 42, 46–50. [Google Scholar] [CrossRef] [Green Version]
  56. Beppu, K.; Syogase, K.; Yamane, H.; Tao, R.; Kataoka, I. Inheritance of self-compatibility conferred by the Se-haplotype of Japanese plum and development of Se-RNase gene-specific PCR primers. J. Hortic. Sci. Biotechnol. 2010, 85, 215–218. [Google Scholar] [CrossRef]
  57. Sapir, G.; Goldway, M.; Stern, R.A.; Shafir, S. S-RNase Based S-genotyping of Japanese Plum (Prunus salicina Lindl.) Cultivars. Nucleotide Sequence EU113311.1. Available online: https://www.ncbi.nlm.nih.gov/nuccore (accessed on 28 December 2022).
  58. Sapir, G.; Goldway, M. Prunus salicina S-RNases. Nucleotide Sequence AF432417.1. Available online: https://www.ncbi.nlm.nih.gov/nuccore (accessed on 28 December 2022).
  59. Beppu, K.; Kumai, M.; Yamane, H.; Tao, R.; Kataoka, I. Molecular and genetic analyses of the S-haplotype of the self-compatible Japanese plum (Prunus salicina Lindl.) ‘Methley’. J. Hortic. Sci. Biotechnol. 2012, 87, 493–498. [Google Scholar] [CrossRef]
  60. Huang, S.X.; Zhang, S.L.; Zhang, Y.Y.; Wu, J.; Heng, W.; Zhao, C.P. Identification of New S Genotype in Prunus salicina. Nucleotide Sequences AY781290.1, AY902455.1, AY996051.1, DQ003310.1. Available online: https://www.ncbi.nlm.nih.gov/nuccore (accessed on 28 December 2022).
  61. Zhang, S.L.; Zhang, S.J.; Huang, S.X.; Wu, H.Q.; Wu, J. Identification of Self-Incompatibility Genotypes of Japanese Plum. Nucleotide Sequences EF177346.1, EF177347.1, EF177348.1, EF177349.1. Available online: https://www.ncbi.nlm.nih.gov/nuccore (accessed on 28 December 2022).
  62. Zhang, S.L.; Zhang, S.J.; Huang, S.X.; Heng, W.; Wu, H.Q.; Wu, J. Identification of Self-Incompatibility Genotypes of Japanese plum by S-Allele-Specific PCR Analysis. Nucleotide sequences EU113259.1, EU113260.1, EU113261.1, EU113262.1, EU113263.1, EU113264.1, EU113265.1, EU113266.1, EU113267.1. Available online: https://www.ncbi.nlm.nih.gov/nuccore (accessed on 28 December 2022).
  63. Wuyun, T.N.; Li, H.G. Isolation of Six SFB New Genes in the Prunus armeniaca (China Apricot) Nucleotide Sequence FJ377732.1. Available online: https://www.ncbi.nlm.nih.gov/nuccore (accessed on 28 December 2022).
  64. Wu Yun, T.N.; Li, H.G. Isolation of Six New S-RNase Genes in the Prunus armeniaca. Nucleotide Sequence GU574195.1. Available online: https://www.ncbi.nlm.nih.gov/nuccore (accessed on 28 December 2022).
  65. Na, T.; Yun, W.; Li, H.G. Isolation of Six New S-RNase Genes in the Prunus salicina. Nucleotide Sequence GU968758.1. Available online: https://www.ncbi.nlm.nih.gov/nuccore (accessed on 28 December 2022).
  66. Ushijima, K.; Sassa, H.; Tao, R.; Yamane, H.; Dandekar, A.M.; Gradziel, T.M.; Hirano, H. Cloning and characterization of cDNAs encoding S-RNases from almond (Prunus dulcis): Primary structural features and sequence diversity of the S-RNases in Rosaceae. Mol. Gen. Genet. 1998, 260, 261–268. [Google Scholar] [CrossRef]
  67. Edgar, R.C. Muscle5: High-accuracy alignment ensembles enable unbiased assessments of sequence homology and phylogeny. Nat. Commun. 2022, 13, 6968. [Google Scholar] [CrossRef] [PubMed]
  68. Guerra, M.E.; Rodrigo, J.; López-Corrales, M.; Wünsch, A. S-RNase genotyping and incompatibility group assignment by PCR and pollination experiments in Japanese plum. Plant Breed. 2009, 128, 304–311. [Google Scholar] [CrossRef]
  69. Okie, W.R.; Hancock, J.F. Plums. In Temperate Fruit Crop Breeding, Hancock, J.F., Ed.; Springer: Dordrecht, Germany, 2008; pp. 337–357. [Google Scholar]
  70. Faust, M.; Surányi, D. Origin and Dissemination of Plums. In Horticultural Reviews; Janick, J., Ed.; John Wiley & Sons Inc.: New York, NY, USA, 1999; pp. 179–231. [Google Scholar]
  71. Okie, W.R.; Weinberger, J.H. Plums. In Fruit Breeding, Volume I: Tree and Tropical Fruits; Janick, J., Moore, J.N., Eds.; John Wiley and Sons, Inc.: New York, NY, USA, 1996; Volume 1, pp. 559–608. [Google Scholar]
  72. Wei, X.W.; Shen, F.; Zhang, Q.; Liu, N.; Zhang, Y.; Xu, M.; Liu, S.; Zhang, Y.; Ma, X.; Liu, W. Genetic diversity analysis of Chinese plum (Prunus salicina L.) based on whole-genome resequencing. Tree Genet. Genomes 2021, 17, 26. [Google Scholar] [CrossRef]
  73. Liu, C.; Feng, C.; Peng, W.; Hao, J.; Wang, J.; Pan, J.; He, Y. Chromosome-level draft genome of a diploid plum (Prunus salicina). GigaScience 2020, 9, giaa13. [Google Scholar] [CrossRef] [PubMed]
  74. Huang, Z.; Shen, F.; Chen, Y.; Cao, K.; Wang, L. Chromosome-scale genome assembly and population genomics provide insights into the adaptation, domestication, and flavonoid metabolism of Chinese plum. Plant J. 2021, 108, 1174–1192. [Google Scholar] [CrossRef] [PubMed]
  75. Tan, Q.; Li, S.; Zhang, Y.; Chen, M.; Wen, B.; Jiang, S.; Chen, X.; Fu, X.; Li, D.; Wu, H.; et al. Chromosome-level genome assemblies of five Prunus species and genome-wide association studies for key agronomic traits in peach. Hortic. Res. 2021, 8, 213. [Google Scholar] [CrossRef]
  76. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  77. Cock, P.J.A.; Chilton, J.M.; Grüning, B.; Johnson, J.E.; Soranzo, N. NCBI BLAST+ integrated into Galaxy. GigaScience 2015, 4, 39. [Google Scholar] [CrossRef] [Green Version]
  78. Afgan, E.; Baker, D.; Batut, B.; van den Beek, M.; Bouvier, D.; Cech, M.; Chilton, J.; Clements, D.; Coraor, N.; Grüning, B.A.; et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses. Nucleic Acids Res. 2018, 48, W395–W402. [Google Scholar]
  79. Andrews, S. FastQC: A Quality Control Tool for High Througput Sequence Data. Available online: http://www.bioinformataics.babraham.ac.uk/proyects/fastqc (accessed on 28 December 2022).
  80. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  81. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Danecek, P.; Bonfield, J.K.; Liddle, J.; Marshall, J.; Ohan, V.; Pollard, M.O.; Whitwham, A.; Keane, T.; McCarthy, S.A.; Davies, R.M.; et al. Twelve years of SAMtools and BCFtools. GigaScience 2021, 10, giab008. [Google Scholar] [CrossRef]
  83. Breese, M.R.; Liu, Y. NGSUtils: A software suite for analyzing and manipulating next-generation sequencing datasets. Bioinformatics 2013, 29, 494–496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Robinson, J.T.; Thorvaldsdóttir, H.; Winckler, W.; Guttman, M.; Lander, E.S.; Getz, G.; Merisov, J.P. Integrative genomics viewer. Nat. Biotech. 2011, 29, 24–26. [Google Scholar] [CrossRef]
  85. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  87. Cao, K.; Li, Y.; Deng, C.H.; Gardiner, S.E.; Zhu, G.; Fang, W.; Chen, C.; Wang, X.; Wang, L. Comparative population genomics identified genomic regions and candidate genes associated with fruit domestication traits in peach. Plant Biotechnol. J. 2019, 17, 1954–1970. [Google Scholar] [CrossRef] [Green Version]
  88. Fang, B.; Xu, Y.; Yu, J. The complete chloroplast genome sequence of Prunus salicina ‘Wushan plum’. Mitochondrial DNA B. 2021, 26, 1243–1244. [Google Scholar] [CrossRef]
  89. Wang, X.; Liu, S.; Zuo, H.; Zheng, W.; Zhang, S.; Huang, Y.; Pingcuo, G.; Ying, H.; Zhao, F.; Li, Y.; et al. Genomic basis of high-altitude adaptation in Tibetan Prunus fruit trees. Curr. Biol. 2021, 31, 3848–3860. [Google Scholar] [CrossRef]
  90. Cao, K.; Peng, Z.; Zhao, X.; Li, Y.; Liu, K.; Arús, P.; Fang, W.; Chen, C.; Wang, X.; Wu, J.; et al. Chromosome level genome assemblies of four wild peach species provide insights into genome evolution and genetic basis of stress resistance. BMC Biol. 2022, 139. [Google Scholar] [CrossRef]
  91. Xu, Y.; Fang, B.; Wang, Y.; Liu, J.; Liu, C.; Yu, J. Phylogenomic analysis and development of molecular markers for the determination of twelve plum cultivars (Prunus, Rosaceae). BMC Genomics 2022, 23, 745. [Google Scholar] [CrossRef]
  92. Fiol, A.; Garcí a-Gómez, B.E.; Jurado-Ruiz, F.; Alexiou, K.; Howard, W.; Aranzana, M.J. Characterization of Japanese Plum (Prunus salicina) PsMYB10 Alleles Reveals Structural Variation and Polymorphisms Correlating with Fruit Skin Color. Frontiers Plant Sci. 2021, 12, 655267. [Google Scholar] [CrossRef] [PubMed]
  93. Fiol, A.; Jurado-Ruiz, F.; López-Girona, E.; Aranzana, M.J. An efficient CRISPR-Cas9 enrichment sequencing strategy for characterizing complex and highly duplicated genomic regions. A case study in the Prunus salicina LG3 MYB10 genes cluster. Plant Methods 2022, 18, 105. [Google Scholar] [CrossRef] [PubMed]
  94. Fang, Z.-Z.; Wang, K.-L.; Dai, H.; Zhou, D.-R.; Jiang, C.-C.; Espley, R.V.; Deng, C.; Lin, Y.-J.; Pan, S.-L.; Ye, X.-F. The genome of low-chill Chinese plum “Sanyueli” (Prunus salicina Lindl.) provides insights into the regulation of the chilling requirement of flower buds. Mol. Ecol. Resour. 2022, 22, 1919–1938. [Google Scholar] [CrossRef] [PubMed]
  95. Numaguchi, K.; Akagi, T.; Kitamura, Y.; Ishikawa, R.; Ishii, T. Interspecific introgression and natural selection in the evolution of Japanese apricot (Prunus mume). The Plant J. 2020, 104, 1551–1567. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hedhly, A.; Guerra, M.E.; Grimplet, J.; Rodrigo, J. S-Locus Genotyping in Japanese Plum by High Throughput Sequencing Using a Synthetic S-Loci Reference Sequence. Int. J. Mol. Sci. 2023, 24, 3932. https://doi.org/10.3390/ijms24043932

AMA Style

Hedhly A, Guerra ME, Grimplet J, Rodrigo J. S-Locus Genotyping in Japanese Plum by High Throughput Sequencing Using a Synthetic S-Loci Reference Sequence. International Journal of Molecular Sciences. 2023; 24(4):3932. https://doi.org/10.3390/ijms24043932

Chicago/Turabian Style

Hedhly, Afif, María Engracia Guerra, Jerome Grimplet, and Javier Rodrigo. 2023. "S-Locus Genotyping in Japanese Plum by High Throughput Sequencing Using a Synthetic S-Loci Reference Sequence" International Journal of Molecular Sciences 24, no. 4: 3932. https://doi.org/10.3390/ijms24043932

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop