Skip to main content

Towards a molecular taxonomic key of the Aurantioideae subfamily using chloroplastic SNP diagnostic markers of the main clades genotyped by competitive allele-specific PCR



Chloroplast DNA is a primary source of molecular variations for phylogenetic analysis of photosynthetic eukaryotes. However, the sequencing and analysis of multiple chloroplastic regions is difficult to apply to large collections or large samples of natural populations. The objective of our work was to demonstrate that a molecular taxonomic key based on easy, scalable and low-cost genotyping method should be developed from a set of Single Nucleotide Polymorphisms (SNPs) diagnostic of well-established clades. It was applied to the Aurantioideae subfamily, the largest group of the Rutaceae family that includes the cultivated citrus species.


The publicly available nucleotide sequences of eight plastid genomic regions were compared for 79 accessions of the Aurantioideae subfamily to search for SNPs revealing taxonomic differentiation at the inter-tribe, inter-subtribe, inter-genus and interspecific levels. Diagnostic SNPs (DSNPs) were found for 46 of the 54 clade levels analysed. Forty DSNPs were selected to develop KASPar markers and their taxonomic value was tested by genotyping 108 accessions of the Aurantioideae subfamily. Twenty-seven markers diagnostic of 24 clades were validated and they displayed a very high rate of transferability in the Aurantioideae subfamily (only 1.2 % of missing data on average). The UPGMA from the validated markers produced a cladistic organisation that was highly coherent with the previous phylogenetic analysis based on the sequence data of the eight plasmid regions. In particular, the monophyletic origin of the “true citrus” genera plus Oxanthera was validated. However, some clarification remains necessary regarding the organisation of the other wild species of the Citreae tribe.


We validated the concept that with well-established clades, DSNPs can be selected and efficiently transformed into competitive allele-specific PCR markers (KASPar method) allowing cost-effective highly efficient cladistic analysis in large collections at subfamily level. The robustness of this genotyping method is an additional decisive advantage for network collaborative research. The availability of WGS data for the main “true citrus” species should soon make it possible to develop a set of DSNP markers allowing very fine resolution of this very important horticultural group.


Chloroplast DNA is a primary source of molecular variations for phylogenetic analysis of photosynthetic eukaryotes. The chloroplast genome has a simple and stable genetic structure, it is haploid, there are no (or very rare) recombination, it is generally uniparentally transmitted. It is known to be highly conserved in both gene order and gene content [1] with substitution rate much lower than nuclear DNA [2]. Therefore, universal primers can be used to amplify targeted sequences for molecular diversity analysis. Several fragments of coding regions, introns, and intergenic spacers, including atpB, atpB-rbcL, matK, ndhF, rbcL, rpl16, rps4-trnS, rps16, trnH-psbA, trnL-F, trnS-G, etc., have been used for phylogenetic reconstructions at various taxonomic levels [38]. Some of these regions such as matK, rbcL, and trnH-psbA have been relied upon heavily for development of candidate markers for plant DNA barcoding [9]. The matK gene is one of the most used sequence so far, because it is useful for identification at family, genus, and even species levels and trnH-psbA is the most variable region in the chloroplast genome across a wide range of groups [9]. Different kinds of molecular markers of chloroplast polymorphisms have been developed such as Restriction Fragment Length Polymorphisms (RFLPs), Cleaved Amplified Polymorphism Sequence (CAPS), Single Sequence Repeats (SSRs). More recently complete chloroplast genome sequences and whole genome resequencing data allowed very fine phylogenetic analysis [10]. Targeted or whole genome sequencing data offer the opportunity to identify polymorphisms diagnostic of well-defined clades. Our hypothesis is that the high level of conservation of the chloroplast sequence and the identification of clade-diagnostic Single Nucleotide Polymorphisms (SNPs) should allow developing molecular taxonomic keys at subfamily or family level, based on an easy, scalable and low-cost method of SNP marker genotyping. These diagnostic set of markers should then be applied in large germplasm collections. In this work we tested this hypothesis within the Aurantioideae subfamily (Rutaceae).

The Rutaceae family comprises seven subfamilies [11], of which the Aurantioideae is the largest group [12]. According to the Swingle and Reece classification [13] largely adopted by the citrus research community, the Aurantioideae subfamily comprises 33 genera and 210 species native to Africa, Australia, North and South America and Asia. They grow in varied climates from equatorial, hot-humid to cool maritime conditions. It is subdivided into two tribes: the Clauseneae with five genera and the Citreae with 28 genera. Swingle and Reece [13] believed that the Clauseneae tribe contained the more primitive genera of the subfamily. According to their classification, the Clauseneae tribe included three subtribes (Micromelinae, Clauseninae and Merrilliinae) and the Citreae tribe also includes three subtribes (Triphasilinae, Balsamocitrinae and Citrinae). The Citrinae is the most important subtribe of Citreae and comprise three groups: the “primitive citrus fruit trees”, the “near citrus fruit trees” and the “true citrus fruit trees”. The third group is formed by the Citrus genus and its closely related genera: Fortunella Swingle, Microcitrus Swingle, Eremocitrus Swingle, Clymenia Swingle and Poncirus Raf [13] [14]. Within the edible Citrus, molecular studies [1523] identified four basic species: C. medica L. (citrons), C. maxima (Burm.) Osbeck (pummelos), C. reticulata Blanco (mandarins) and C. micrantha Wester (papeda). All other cultivated species, such as C. sinensis (sweet oranges), C. aurantium (sour oranges), C. paradisi (grapefruits), C. lemon (lemons) and C. aurantifolia (limes) result from these four ancestral taxa by reticulate evolution.

However, tribal and subtribal classifications of the Aurantioideae subfamily have barely been debated. Tanaka [24] grouped the subfamily into eight tribes and eight subtribes including 28 genera, while Engler [11] grouped all members of the subfamily into a single tribe, Aurantieae, subdivided into two subtribes, the Hesperethusinae (16 genera) and the Citrinae (13 genera). More recently, Mabberley [25] fused Eremocitrus, Fortunella, and Microcitrus with Citrus, and suggested that Poncirus should be so treated as well. These three classifications, as well as Swingle and Reece’s [13], are based on morphological traits generally influenced by environments and require strong human expertize for adequate classification inferences. In recent decades several molecular studies have been carried out to clarify the phylogenetic relationships within the Aurantioideae subfamily. The chloroplast genome of angiosperms is characterised by its evolutionary conservatism, the relative abundance of plant tissue, uniparental inheritance and small size (135 to 160 kb) [26]. Chloroplastic DNA (cpDNA) can avoid the problems of the complicated phylogenetic relationships within the Aurantioideae subfamily, and especially between Citrus and its close relatives, because of hybridization, apomixis and polyploidy [16]. Preliminary studies on the Aurantioideae subfamily using chloroplast DNA sequence data were based on two plastid genes, (the atpB-rbcL intergenic spacer and the rps16 intron) analysed on 15 of the 33 genera from the Aurantioideae subfamily [27]. The polymorphisms provided were not informative enough for resolution or sufficient support for phylogenetic inferences. Araùjo et al. [28] published a phylogenetic study on the tribe Citreae, based on partial sequences from the trnL-F region and some morphological characters. Morton et al. [29] used the rps16 and trnL-trnF introns from 24 genera. The results obtained by Samuel et al. [27], Araùjo et al. [28] and Morton et al. [29] were not fully congruent and did not have enough resolution to address tribal and subtribal delimitations.

Broader studies by Morton [30], and particularly by Bayer et al. [31], testing nine cpDNA gene regions confirmed the monophyly of the Aurantioideae subfamily and suggested none-monophyly for several subtribes. The study by Bayer et al. [31] resulted in a revision of the Swingle and Reece [13] classification and, still today, the more conclusive phylogenetic analysis of the whole Aurantioideae subfamily members.

Several studies tried to clarify the phylogeny within the “true citrus” group. For Bayer et al. [31] the circumscription of Citrus genus was broadened by including seven other closely related genera of the orange subfamily: Clymenia, Fortunella, Poncirus, Microcitrus, Eremocitrus, Oxanthera Montrouz., and Feroniella Swingle. The Citrus genus (as defined by Swingle and Reece [13]) was not monophyletic, with C. medica included in a clade with the Australian and New Caledonian genera. However, another study of the “true citrus fruit trees” group based on three cpDNA fragments (trnL-trnF, psbH-petB and trnS-trnG) confirmed its monophyly [32], but the group divided into six genera as previously proposed by Swingle and Reece [13]. A more definitive picture of the phylogeny of the true citrus group was recently provided by the analysis of complete chloroplast sequence polymorphisms [10] derived from whole genome resequencing data mapped on the sweet orange reference chloroplast genome [33]. It revealed three main clades: the first joining the citron with the Australian species, a second one associating the pummelos with C. micrantha and a third one the mandarins with C. ichangensis, a papeda species. Poncirus and Fortunella appeared as independent units.

Numerous studies analysed the maternal phylogeny within the Citrus genus using different tools. The oldest study was based on plastid Restriction Fragment Length Polymorphism (RFLPs) [26, 34, 35]. Other studies were based on chloroplast gene or spacer sequences such as matK sequences [36], trnL-trnF spacer [37, 38], rbcL-ORF106, trnL-trnF and trnF-trnVr region sequences [39] or psbH-petB, trnL-trnF, rbcL genes [40] and trnS-trnG cpDNA regions [41]. Chloroplastic Simple Sequence Repeats (CpSSRs) were also useful for differentiating the ancestral citrus taxa and for identifying the maternal phylogeny of secondary species [4245]. Cleaved Amplified Polymorphic Sequence (CAPS) markers were also successfully developed [39, 46]. All these studies contributed to establishing the differentiation between the four ancestral taxa of the cultivated citrus (C. reticulata, C. maxima, C. medica, C. micrantha) and their contribution in the maternal phylogeny of the secondary species.

Two complete chloroplastic genome sequences have been published; the first for C. sinensis osbeck L (sweet oranges) [33] and more recently for C. aurantifolia (Omani lime cultivar) [47]. They revealed molecular polymorphisms (Insertion/ Deletion : Indels, Single Nucleotide Polymorphism: SNPs and SSRs) between the two species [47]. Based on whole genome resequencing data for 34 genotypes of the true citrus mapped on the chloroplast genome, Carbonell-Caballero et al. [10] identified 14.5 SNV/kb.

As mentioned before the objective of this work was to develop a molecular taxonomic key of the Aurantioideae subfamily based on an easy and low-cost SNP marker method. The competitive allele-specific PCR method allows efficient, scalable, easy and affordable SNP genotyping. It was therefore chosen for this study. The publicly available nucleotide sequence of eight plastid genomic regions (GeneBank; [31]) were compared for 79 accessions of the Aurantioideae subfamily to search for SNPs revealing taxonomic differentiation at inter-tribe, inter-subtribe, inter-genus and interspecific levels. KASPar markers were developed from the selected SNPs and used to genotype 108 accessions of the Aurantioideae subfamily to validate their taxonomic value, their transferability to the whole subfamily and to specify the maternal phylogeny of some Citrus species and cultivars.


Chloroplastic sequence selection

The sequences of eight chloroplastic regions used by Bayer et al. [31] for their phylogenetic study of the Aurantioideae subfamily, were obtained from the National Center of Biotechnology Information (NCBI reference, can be found in (Additional file 1): atpB-coding region, rbcL-atpB spacer, rps16 spacer, trnL-F region, rps4-trnT spacer, matK-5trnK spacer, psbM-trnD GUC spacer and trnG intron. For each region, the sequences of seventy-nine accessions of the Aurantioideae subfamily were used. Ruta graveolens L. from the Rutaceae family was chosen as the outgroup. All genera of the Aurantioideae subfamily except Limnocitrus were sampled. The Clauseneae tribe was represented by three subtribes: Clauseninae (seven accessions), Merrilliinae (one accession) and Micromelinae (one accession). The Citreae tribe was represented by three subtribes: Triphasilinae (ten accessions), Balsamocitrinae (seven accessions) and Citrinae (fifty-three accessions). For the Citrinae subtribe we tried to select all known species of all genera (Table 1 and Table 2).

Table 1 Number of species/genus and plant accession/genus used for in silico SNP mining and KASPar analysis; classification according to Swingle and Reece [3]
Table 2 Species and number of plant accessions/species of the true citrus used for in silico SNP mining and KASPar analysis; classification according to Swingle and Reece [3]

SNP identification and selection

Using BioEdit software [48], the sequences were aligned for each gene to the C. sinensis reference chloroplast genome sequence [33] available in [GenBank: NC008334]. The accessions were classed according to the consensus trees generated by Bayer et al. [31] to select SNPs between clades on different taxonomic levels. A SNP was considered as diagnostic of a clade when all accessions within the clade were identical and different from all the accessions outside the clade.

SNP loci were further selected for marker development using broad marker design criteria to allow utilization with several SNP chemistries. These criteria were: availability of DNA sequence 50 bp upstream and downstream from the SNP site and 15 bp around the target SNP with no additional molecular polymorphism.

KASPar marker development and application

Plant material

One hundred and eight genotypes of the Aurantioideae subfamily, including five subtribes, namely the Balsamocitrinae, Citrinae, Triphasiinae (from the Citreae tribe) and Clauseninae and Merrilliinae (from the Clauseneae tribe) were sampled for the SNP markers selected. In this study, we adopted the Swingle and Reece subdivision [13]. The Clauseninae subtribe included three genera: Clausena (three accessions), Glycosmis (one accession) and Murraya (two accessions). The Merrilliinae subtribe was represented by the only genus Merrillia with one accession. The Balsamocitrinae subtribe was represented by seven genera: Swinglea (one accession), Aegle (one accession), Afraegle (two accessions), Aeglopsis (one accession), Balsamocitrus (one accession), Feronia (one accession), and Feroniella (one accession). The Triphasilinae subtribe included five genera: Oxanthera, Pamburus, Paramignya, Triphasia, and Wenzelia. One accession represented each genus. Within the Citrinae subtribe three groups were represented: the primitive citrus fruits, the near citrus fruits, and the true citrus fruits. The near citrus fruit group included two genera, Citropsis (four accessions) and Atalantia (four accessions), the primitive citrus fruit group included the genera Severinia (two accessions), Pleiospermium (one accession) and Hesperethusa (one accession), and the true citrus fruit group was represented by six genera, namely Fortunella (six accessions), Eremocitrus (one accession), Poncirus (two accessions), Clymenia (one accession), Microcitrus (eight accessions) and Citrus (fifty-eight accessions) (Table 1, Table 2 and detailed in (Additional file 2).

DNA extraction

High molecular weight genomic DNA was extracted from leaf samples using the DNeasy Plant Mini Kit (Qiagen S.A.; Madrid, Spain) according to the manufacturer’s instructions.

KASPar™ Genotyping

The KASPar™ Genotyping System from LGC Genomics® is a competitive allele-specific dual FRET based assay for SNP genotyping. It uses two Fluorescent Resonance Energy Transfert (FRET) cassettes where fluorometric dye, HEX or FAM, is conjugated to primer but quenched via resonance energy transfer. Sample DNA is amplified with a thermal cycler using allele-specific primers, leading to the separation of fluorometric dye and quencher when the FRET cassette primer is hybridized with DNA [49]. The primers were designed by LGC Genomics® based on the SNP locus-flanking sequence (approximately 50 nucleotides either side of the SNP). Detailed information for all SNP markers can be found in Additional file 3. The fluorescence signals of PCR products were measured with Fluostar Omega (BMG) and genotype calling was done with KlusterCaller software (LGC Genomics).

SNP marker data analysis

Theoretical Clade differentiation

The GST parameter [50] was used to estimate the efficiency of each developed marker to effectively differentiate the clade identified by Bayer et al. [31]. GST estimations were computed using Excel software considering two subpopulations: (1) the accessions theoretically included in the clade in question, (Ci) and (2) a theoretical population of all other accessions (C-i). The analysis was performed from the estimated allele frequency of each group considering the same population size to estimate the frequency of the whole population (Tot) frequency.

$$ GS{T}_{\mathrm{C}\mathrm{ladei}}=\frac{H{e}_{\mathrm{Tot}}-\frac{\left( He{}_{\mathrm{C}\mathrm{i}}+ He{}_{\mathrm{C}\hbox{-} \mathrm{i}}\right)}{2}}{H{e}_{\mathrm{Tot}}} $$

Where He is the genetic diversity within population (He = 1-Σ pi 2, where pi is the frequency of a given allele in the considered population).

Polymorphism Information Content (PIC) was estimated for each marker as PIC = 1-Σ (Pi) 2

Clade analysis: a cladistic analysis was performed from the KASPar SNP data. Unweighted pair group method with the arithmetic average (UPGMA) was performed with Mega version 6 software [51] from a matrix of Euclidian distances between each pair of accessions.

Results and discussion

Clade-specific SNP identification

The available sequences of 79 Aurantioideae accessions for 8 chloroplastic genome regions were aligned and organized according to the Bayer et al. [31] classification (see Additional file 1 for the nomenclature of the Bayer Clades). They were also aligned to the reference chloroplastic genomes of sweet orange [33] in order to locate the selected SNPs on that reference genome. This allowed easy visual identification of SNPs totally differentiating the different clades from the rest of the accessions. Additional file 4 gives an example of such diagnostic SNP (DSNP) mining in the region rps4-trnT spacer for the Bayer clade Q corresponding to the true citrus of Swingle and Reece [13] plus the Oxanthera and Feroniella genera. In position 49301 of the reference chloroplastic genome a SNP (A/G) totally differentiated the accessions of clade Q (A) from all other Aurantioideae accessions (G). The alignments of the height analysed chloroplast regions can be found in the Additional file 5. Of the 54 clade levels analysed, 8 did not display any DSNPs. (one of them displaying a SNP shared by the clade GG1 and the species C. aurantifolia and Hesperethusa crenulata). The majority (28) displayed between one and 3 DSNPs, 15 clades displayed between 4 and 8 DSNPs and only three displayed 12 or more DSNPS. A total of 166 DSNPs was identified (Fig. 1 and Tables 3).

Fig. 1

Distribution of the number of DSNPs for the 54 considered clades from Bayer et al. [21]

Table 3 Number of diagnostic SNPs of the Bayer’s Clade encountered from the in silico mining

For two of the clades with 12 or more DSNPs (Poncirus trifoliata and Clymenia polyandra), these high values could be partially explained by the fact that these taxa were represented by only one accession and that we therefore missed a major part of intra-taxa diversity. However, it also argues in favour of considerable differentiation of these genera within the “true citrus” clade. The third clade with a high DSNP value was the Glycosmis clade with accessions of three species, confirming the conclusion of Bayer et al. [31] and Samuel et al. [27] regarding the strong monophyletism of Glycosmis. The identification of 8 and 7 DSNPs for the clades F (Merrillia and Murraya paniculata), K (Aegle, Aeglopsis, Afraegle and Balsamocitrus) and Q (true citrus plus Feroniella and Oxanthera) validated the robustness of these clades and the monophyletic character of each one. For clade F, it confirmed the close relationship of Merrillia and Murraya paniculata proposed by Bayer et al. [31]. The K clade was defined by Bayer et al. as the true Balsamocitrinae clade corresponding to the bael-fruit group of Swingle and Reece [13]. The monophyletic origin of the Q clade including the sequences of the true Citrus [13], Feroniella and Oxanthera, is also clearly confirmed by these numerous DSNPs. Also with 7 DSNPs, C. macroptera appeared clearly differentiated from other papeda such as C. micrantha and C. celebica, even though they were joined in the DD clade validated with 2 DSNPs. The 6 DSNPs for the OO (Oxanthera) and V (M. australasica; M. Garrawayi) clades revealed substantial differentiation within the Australian/New Caledonian clade. With 5 DSNPs the C. medica/C. indica clade (QQ) was clearly validated. Only one DSNP was found for the R clade joining the C. indica/C. medica clade and the Australian/New Caledonian one (T). It revealed the distended relationship between the two groups, even though the R clade was also inferred from Whole Genome Sequence (WGS) data [10]. The identification of 5 DSNPs for clade E and clade G strongly supported the monophyletic origins of Bergera/Clausena/Micromelum/Glycosmis on the one hand and all the other considered genera on the other hand. Also with 5 DSNPs C. ichangensis appeared clearly differentiated from the other papeda. With 4 DSNPs the monophyletic origins of clade C (Clausena species), D (Bergera and Clausena), L (Merope/Monanthocitrus/Wenzelia), and U (Australian citrus) were also clearly validated. The absence of DSNPs for several clades identified by Bayer et al. [31] corresponded to branching with a low maximum parsimony bootstrap in their consensus phylogenetic tree analysis and highlighted the weakness of these clades. It concerned the monophyletic origin of (i) clade I, Triphasia trifolia, and the strong clade Merope/Monanthocitrus/Wenzelia, (ii) clade O (Feronia and Atalantia) and (true citrus plus 0xanthera and Feroniella). It therefore appears that the resolution of the Citreae tribe into sub-groups remains confused and data from whole chloroplast genome sequences will probably be necessary to definitively conclude on the organisation of this tribe. The monophyletic origin of the subsequent groups is also questioned: (iii) clade T: Clymenia polyandra and the Australian/New Caledonian clade, (iv) clade Y joining the clades BB (mandarins) and Z (pummelos and papedas). The other branching with no DSNPs concerned more recent differentiation within closely related accessions (L1; GG1; MM; NN).

Of the 166 DSNPs identified, 40 SNPs diagnostic of 33 clades (Additional file 3) were selected. For each clade the main selection criterion was the absence of additional polymorphisms (SNPs or Indels) in close vicinity (ideally in the 30 pb flanking the SNP each side; but at least in the 15 pb for clades when the previous condition was not encountered).

KASPar analysis

Primers for the 40 selected SNPs were designed by LGC Genomics® based on the SNP locus-flanking sequence (approximately 50 nucleotides either side of the SNP (Additional file 3). Eleven of these markers provided an inconsistent or no signal and were not useful: CP_3707 (clade DD), CP_10192 (clade U), CP_32282 (clade N), CP_5927 (clade EE), CP_51032 (clade II), CP_32684 (clade Z), CP_57022 (clade A), CP_58121 (clade J), CP_51540 (clade JJ), CP_57256 (clade L), CP_3996 (clade Q). The CP_57686 marker, theoretically diagnostic of Ichang papeda, provided a good signal but was not polymorphic so it was discarded. Moreover, CP_57742, which was theoretically diagnostic of Clymenia, differentiated Atalantia ceylanica but not Clymenia from all the other accessions and was also discarded.

A genetic analysis was performed with the 27 remaining markers diagnostic of 24 Bayer et al. [31] clades, for 108 accessions (Table 4). Figure 2 illustrates the clear differentiation obtained with the CP_49301 marker between the true citrus accessions of Swingle and Reece, plus Oxanthera and Clymenia, and all the other genera. The rate of missing data for these markers was very low, ranging from 0 to 2.8 % with an average of 1.2 %. This situation is very different from the one encountered for SNP KASPar nuclear markers mined from a panel of Citrus accessions where the frequency of missing data was higher for the citrus relatives and increased with taxonomic distances within the Aurantioideae subfamily [20]. This showed that selection based on the absence of additional polymorphism in close vicinity to the selected SNPs was efficient. This very good transferability also resulted from the much lower rate of chloroplast sequence evolution compared to nuclear evolution. The PIC varied from 0.02 to 0.5 in relation to the taxonomic level revealed by each marker, but also to our sample with an overrepresentation of the true citrus compared to the other genera.

Table 4 KASPar marker results
Fig. 2

An example of KASPar analysis (marker CP_49301) distinguishing the accessions of the Bayer’s Clade Q. Q clade represent true Citrus of Swingle and Reece [3] plus Oxanthera. Red dots: SNP “A”, Blue dots: SNP “G”, Black dot: blank data, and Pink dot: missing data

The GST analysis (the clade theoretically diagnosed by the SNP considered as one population and the rest of the sample as a second one) globally validated the cladistic diagnostic value of the selected markers with an average value of 0.98. Twenty of the markers totally fitted (GST = 1). For the marker of the CC clade (CP_4159) the Yuzu (C. junos) was different from the mandarins, while it was included in the Bayer et al. analysis [31]. Our results are in agreement with the results of Abkenar et al. [52] and Yamamoto et al. [39] who found the Yuzu differentiated from mandarins by Chloroplastic CAPS analysis. Yuzu was found to be closely associated with Ichang papeda in several studies [37, 39]. Therefore, except for Yuzu which still has a debated status, this marker appears as a perfect univocal diagnostic marker of the mandarin group. For the M clade the two selected markers displayed one outgroup accession with the same allele as clade M (Clausena excavata for CP_3829 and Severinia disticha for CP_49405). CP_58062 for clade X (Fortunella Sp and C. halimii) also displayed one accession outside the clade (Atalantia ceylanica), sharing the clade allele. For CP_10342 (clade AA; Fortunella sp.) the diagnostic allele was shared with two accessions of clade D (Clausena anisata and Murraya koenigii) in addition to the 6 Fortunella species analysed. For M. koenigii it confirmed the results of the SNP mining from the Bayer et al. [31] sequence alignments. For CP_3982 (clade P; Limonia and Atalantia) the diagnostic allele was shared with two accessions of clade N (Hesperethusa and Pleiospermium). For CP_3747 (clade N) displaying the lower GST value (0.73) one accession of the theoretical clade shared the out-clade allele while Severinia buxifolia shared the clade allele. The results of the last two markers confirmed the weakness of the resolution of the Citrinae clade into consistent subclades.

An UPGMA analysis was performed with the 27 SNP markers (Fig. 3). The representation was highly coherent with the clade of Bayer et al. [31].

Fig. 3

Classification of 108 accessions of the Aurantioideae subfamily using 27 SNP Bayer’s Clade diagnostic markers; the different colors correspond to different clade levels according to the Bayer et al. classification (2009)

The six genera of the true citrus of Swingle and Reece [13] were included in the Q clade as well as the representative of the genus Oxanthera as proposed by Bayer et al. [31]. However, the representative of the genus Feroniella (F. oblata) was found in the M’ clade, closely associated with representatives of the Balsamocitrinae (as in the Swingle and Reece classification [13]), primitive and near citrus fruits, and is not in the Q clade as concluded by Bayer et al. [31]. Therefore, the genetic conformity of the accession of Feroniella used in Bayer et al. [31], or in our study, can be questioned.

Within the Q clade the general organisation globally tallied with 15 of the Bayer clades encountered (with only the position of Yuzu and Ichang papeda in the pummelo clade being different). The Poncirus clade (RR) was clearly differentiated. The mandarin clade (JJ) was largely extended to 14 species of the Tanaka classification (C. depressa, C. tangerina, C. daoxianensis, C. clementina, C. kinokuni, C. erythrosa, C. suhuiensis, C. reticulata, C. reshni, C. nobilis, C. paratangerina, C. sunki C. deliciosa, C. unshiu); the supposed natural tangor Ortanique and Murcott also share this clade as concluded by Garcia-Lor et al. [43] from chloroplastic SSRs and mitochondrial indels. Three acid citrus (Rangpur lime, Rough lemon and Volkameriana) were also associated with this clade in agreement with Bayer et al. [31] for Rangpur lime and Rough lemon. For these three acid citrus used as rootstocks, we confirmed the conclusions of Froelicher et al. [18] with mitochondrial indels, and of Curk et al. [45] derived from a large combined cytoplasmic and nuclear analysis demonstrating that these three limes and lemons result from direct hybridization between the mandarin and citron gene pools. With the available data from SNP mining, we were not able to develop markers of the differentiation between sweet and sour mandarin observed with mitochondrial data [18], chloroplastic SSRs [43] and WGS data [10]. The recent availability of these WGS data will certainly make it possible to develop SNP KASPar markers for this differentiation. The pummelo clade (EE) included all the analysed pummelos, sour oranges, sweet oranges, and grapefruits, Meyer and Lisbon lemons, Ichang papeda and Yuzu. The pummelo maternal phylogeny of sour and sweet oranges, and grapefruit was already described by several authors [10, 1618, 3941, 53]. Nicolosi et al. [16] proposed that C. lemon should be a hybrid between sour orange and citron and it was definitively demonstrated by Curk et al. [45], as well as the pummelo maternal phylogeny of Meyer lemon and Bergamot. The development of two additional SNP markers diagnostic of the differentiation (i) between sour oranges and pummelos, already revealed by chloroplastic SSRs [45] and WGS data [10] and (ii) between sweet oranges pummelo and sour orange revealed by WGS [10] would allow a very fine definition of this clade with few SNP markers. The FF clade (Mexican lime) included also the Combava (C. hystrix), C. micrantha, C. macrophylla and C. amblycarpa. These results are in agreement with those of Froelicher et al. [18], Garcia Lor et al. [43] and Curk et al. [45] from mitochondrial indels and chloroplastic SSRs. The six analysed Fortunella species were identical (clade AA) and were grouped in clade X with C. halimii according to Bayer et al. [31]. These findings confirm the studies of Roose [54], Scora [55], the isozymatic data of Herrero et al. [56], the chemotaxonomy study of Ogawa et al. [57], along with cpRFLP data [58] and nuclear gene data [59]. This genetic relationship should be related to the occurrence of Fortunella and C. halimii in the Malaysian Peninsula [58]. All citrons were grouped in clade QQ which was joined in clade R with all the Australian and New Caledonian citrus. This relationship between C. medica and the Australian/New Caledonian citrus was first revealed by Pfeil et al. [60] and Bayer et al. [31] and recently confirmed by a chloroplastic WGS analysis [10]. Within the T clade of the Australian and New Caledonian citrus, the organisation displayed by Bayer et al. [31] was mostly conserved within clades V and W.

The M’ clade corresponded to the M clade of Bayer et al. [31] but it did not include the Q clade of the true Citrus and Oxanthera. This may have been due to the selection of two DSNPs of clade Q. M is the clade where few DSNPs of the Bayer sub-clade were found and where the within genetic organization was not in agreement with Bayer et al. [31] nor with Swingle and Reece [13]. More molecular data would be necessary to have good resolution for this group. Similarly, Penjor et al. [36] and Morton et al. [29] did not support the distinction in the Citrinae sub-tribe between “primitive citrus fruit trees” and “near citrus fruit trees” made by Swingle and Reece [13].

The main clades for the rest of the Citreae and the Clauseneae is respected compared with the Bayer et al. [31] classification with (i) the K clade corresponding to Afraegle, Aeglopsis, Balsamocitrus and Aegle; (ii) the I clade joining Wenzelia and Triphasia trifolia; (iii) the J clade with Pamburus and Paramignya; (iv) the F clade with Merrillia and Murraya paniculata and finally (v) the Clauseneae clade (E) with Murraya koenigii, Clausena Sps and Glycosmis pentaphylla. An in-depth discussion on the meaning of this cladogram and its relations with previous data can be found in Bayer et al. [31]. To date their work remains the most in-depth at this taxonomic level.

Usefulness and complementarity of KASPar DSNPs analysis compared with others molecular markers

Although cpDNA analysis for phylogenetic studies is especially based on noncoding regions as well as gene sequence data analysis [37, 38, 41, 52, 61], several studies analysed Aurantioideae maternal phylogeny using different molecular markers such as RFLPs [26, 34, 35, 38, 52, 53, 58], cpSSRs [4244], Single Strand Conformation Polymorphisms (SSCPs) [62, 63] and CAPS [39, 46]. The complete sequence chloroplast genome of C. sinensis [33] led on to in-depth studies such as the use of WGS to elucidate citrus phylogenetic relationships and divergence time estimations [10], and opened up the way for SNP marker development over the whole chloroplast genome. SNPs are the most abundant type of sequence variation in eukaryotic genomes [64, 65] and are considered to be the simplest, the ultimate and the smallest unit of inheritance. Chloroplast SNPs where identified in Populus species and a CAPS approach was developed to reveal polymorphisms between seven species [66]. The main disadvantages of the CAPS approach are low throughput and the relatively high cost of genotyping. A large set of soybean chloroplast SNP markers was identified, selected and successfully included in a genotyping array [67], as has been done for potato, Solanum tuberosum [68]. However, the use of genotyping arrays remains costly and does not allow flexibility for adaptation to specific research questions needing only a limited specific subset of markers. In our lab the genotyping of one sample for one diagnostic SNP with the KASPar method cost around 0.10 euros per plant sample. For the first time in plant taxonomy, we have shown proof of the concept that a well-selected set of SNPs at subfamily level combined with an efficient competitive allele-specific PCR method (KASPar technology) opens up the way for implementing a highly effective, simple and low-cost molecular taxonomic key.

For recently diverged groups such as the “true citrus” that include the cultivated forms, the recent release of WG chloroplast data for 34 accessions and the corresponding phylogenetic analysis [10] testify that diagnostic polymorphisms of the differentiation between the different genera and all Swingle and Reece Citrus species could be found. Therefore, regarding the “true citrus”, for the maternal phylogeny point of view, the set of marker may be complemented and the developed approach will provide a very efficient cladistics molecular key at interspecific level. It is clear, however, that the chloroplast phylogeny reveal only a part of the evolutive history of a gene pool due its mono-parental heredity and should not be enough resolutive for very recently diverged group due to its low rate of variation. Therefore chloroplast and nuclear data are totally complementary to clarify phylogenomic structure and evolutionary stories of recent divergent groups particularly when reticulation events were frequents. It is the case for cultivated Citrus arising from hybridisation between four ancestral taxa and displaying nuclear interspecific admixture [20, 22]. Diagnostic nuclear markers of these four taxa coupled with cytoplasmic information were powerful to decipher the origin of citrus secondary species [23, 45] and WGS data provided a detailed picture of the nuclear interspecific mosaic phylogenomic structures of some cultivars [69]. Next Generation Sequencing (NGS) coupled with a reduction of genome complexity using methods such as Genotyping By Sequencing (GBS) or Restriction site Associated DNA Sequencing RAD- seq appears as a good option to provide phylogenomic pictures over the whole nuclear genome complementary to the information given by CpDNA DSNPs.


CpDNA analysis is a powerful approach for phylogenetic studies on a wide taxonomic level. Monoparental inheritance is particularly useful in a gene pool with reticulate evolution, such as cultivated citrus, to infer the maternal phylogeny. Several methods have been developed for cpDNA analysis, such as CAPS and SSR markers, amplicon sequencing and more recently the use of WGS data mapped in the chloroplast reference genome. In this work, for the first time, we propose using SNPs that are diagnostic of the different clades to establish a molecular taxonomic key of the Aurantioideae. The concept was to focus the molecular analysis on targeted polymorphisms with high taxonomic value and to analyse those polymorphisms with a simple, cost effective, but highly robust genotyping method having wide transferability across the Aurantioideae subfamily. Diagnostic-Clade SNP mining was performed in silico from publicly available cpDNA sequences considering the classification and clades identified by Bayer et al. [21] as the reference. We then tested the efficiency of the competitive allele-specific PCR method (KASPar) to achieve our objectives of wide transferability. From the 40 DSNPs selected from the in silico mining analysis, 27 KASPar markers, diagnostic of 24 Bayer clades, were successfully developed with a very high rate of transferability in the Aurantioideae subfamily (only 1.2 % of missing data on average). Twenty-one key clade markers, univocally diagnostic of 19 clades, were identified. We have demonstrated proof of the concept that with well-established clades, DSNPs can be selected and efficiently transformed into competitive allele-specific PCR markers allowing cost-effective, highly efficient cladistic analysis in large collections at subfamily level. The robustness of the method is a decisive advantage compared with other kinds of cpDNA markers for network collaborative research. The application is considerably easier (laboratory work and data analysis) than a comparison of targeted sequences of amplicons. The availability of WGS data for the main true citrus species may soon enable the development of a set of DSNP markers allowing very fine resolution.



Cleaved amplified polymorphic sequence


ChloroPlastic DNA


Deoxyribo nucleic acid


Diagnostic SNP


Fluorescent resonance energy transfert




National center of biotechnology information


Polymerase chain reaction


Polymorphism information content


Restriction fragment length polymorphism


Single nucleotide polymorphism


Single strand conformation polymporphisms


Simple sequence repeats


Unweighted pair group method with the arithmetic average


Whole genome sequence


  1. 1.

    Raubeson LA, Jansen RK. Chloroplast genomes of plants. In: Henry RJ, editor. Diversity and evolution of plants—genotypic and phenotypic variation in higher plants. Wallingford: CABI Publishing; 2005. p. 45–68.

    Google Scholar 

  2. 2.

    Wolfe KH, Li WH, Sharp PM. Rates of nucleotide substitution vary greatly among plant mitochondrial, chloroplast, and nuclear DNAs. Proc Natl Acad Sci U S A. 1987;84:9054–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Gao X, Zhu YP, Wu BC, Zhao YM, Chen JQ, Hang YY. Phylogeny of Dioscorea sect. Stenophora based on chloroplast matK, rbcL and trnL-F sequences. J Syst Evol. 2008;46:315–21.

    Google Scholar 

  4. 4.

    Li JH. Phylogeny of Catalpa (Bignoniaceae) inferred from sequences of chloroplast ndhF and nuclear ribosomal DNA. J Syst Evol. 2008;46:341–8.

    Google Scholar 

  5. 5.

    Wilson CA. Phylogenetic relationships among the recognized series in Iris Section Limniris. Syst Bot. 2009;34:277–84.

    Article  Google Scholar 

  6. 6.

    Peterson PM, Romaschenko K, Johnson G. A classification of the Chloridoideae (Poaceae) based on multi-gene phylogenetic trees. Mol Phylogenet Evol. 2010;55:580–98.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Hilu KW, Black C, Diouf D, Burleigh JG. Phylogenetic signal in matK vs. trnK: A case study in early diverging eudicots (angiosperms). Mol Phylogenet Evol. 2008;48:1120–30.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Kim KJ, Jansen RK. NdhF sequence evolution and the major clades in the sunflower Family. Proc Natl Acad Sci U S A. 1995;92:10379–83.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Dong W, Liu J, Yu J, Wang L, Zhou S. Highly Variable Chloroplast Markers for Evaluating Plant Phylogeny at Low Taxonomic Levels and for DNA Barcoding. PLoS ONE. 2012;7(4):e35071. doi:10.1371/journal.pone.0035071.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Carbonell-Caballero J, Alonso R, Ibanez V, Terol J, Talon M, Dopazo J. A phylogenetic analysis of 35 chloroplast genomes elucidates the relationships between wild and domestic species within the genus Citrus. Mol Biol Evol. 2015;32:2015–35.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Engler A. Rutaceae. In: Engler A, Prantl K, editors. Die Naturlichen Pflanzenfamilian? Germany: Engelmann Leipzig; 1931. p. 187–359.

    Google Scholar 

  12. 12.

    Thorne RF. The classification and geography of the flowering plant: Dicotyledons of the class Angiospermae (subclasses Magnoliidae, Ranunculidae, Caryphyllidae, Dilleniidae, Rosidae, Asteridae, and Lamiidae). Bot Rev. 2000;66:441–647.

    Article  Google Scholar 

  13. 13.

    Swingle WT, Reece PC. The botany of Citrus and its wild relatives in the orange subfamily. In: Reuther W, Webber HJ, Batchelor DL, editors. The Citrus industry. Berkeley: University of California; 1967. p. 190–340.

    Google Scholar 

  14. 14.

    Krueger RR, Navarro L. Citrus germplasm resources. In: Khan I, editor. Citrus genetics, breeding and biotechnology. Wallington: CABI Publishing; 2007. p. 45–140.

    Google Scholar 

  15. 15.

    Federici CT, Fang DQ, Scora RW, Roose ML. Phylogenetic relationships within the genus Citrus Rutaceae and related genera as revealed by RFLP and RAPD analysis. Theor Appl Genet. 1998;96:812–22.

    CAS  Article  Google Scholar 

  16. 16.

    Nicolosi E, Deng ZN, Gentile A, La Malfa S, Continella G, Tribulato E. Citrus phylogeny and genetic origin of important species as investigated by molecular markers. Theor Appl Genet. 2000;100:1155–66.

    CAS  Article  Google Scholar 

  17. 17.

    Barkley NA, Roose ML, Krueger RR, Federici CT. Assessing genetic diversity and population structure in a Citrus germplasm collection utilizing simple sequence repeat markers SSRs. Theor Appl Genet. 2006;112:1519–31.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Froelicher Y, Mouhaya W, Bassene JB, Costantino G, Kamiri M. New universal mitochondrial PCR markers reveal new information on maternal citrus phylogeny. Tree Genet Genomes. 2011;7:49–61.

    Article  Google Scholar 

  19. 19.

    Garcia-Lor A, Luro F, Navarro L, Ollitrault P. Comparative use of Indel and SSR markers in deciphering the interspecific structure of cultivated Citrus genetic diversity: a perspective for genetic association studies. Mol Genet Genomics. 2012;2871:77–94.

    Article  Google Scholar 

  20. 20.

    Garcia-Lor A, Curk F, Snoussi-Trifa H, Morillon R, Ancillo G, Luro F, et al. A nuclear phylogenetic analysis: SNPs, indels and SSRs deliver new insights into the relationships in the ‘true citrus fruit trees’ group (Citrinae, Rutaceae) and the origin of cultivated species. Ann Bot-London. 2013;111:1–19.

    CAS  Article  Google Scholar 

  21. 21.

    Ollitrault P, Terol J, Garcia-Lor A, Aurélie B, Aurélie C, Yann F, et al. SNP mining in C.clementina BAC end sequences; transferability in the Citrus genus (Rutaceae), phylogenetic inferences and perspectives for genetic mapping. BMC Genomics. 2012;13:13.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Curk F, Ancillo G, Garcia-Lor A, Luro F, Perrier X, Jacquemoud-Collet JP, et al. Next generation haplotyping to decipher nuclear genomic interspecific admixture in Citrus species: analysis of chromosome 2. BMC Genet. 2014;15:152.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Curk F, Ancillo G, Ollitrault F, Perrier X, Jacquemoud-Collet JP, Garcia-Lor A, et al. Nuclear species-diagnostic snp markers mined from 454 amplicon sequencing reveal admixture genomic structure of modern citrus varieties. PLoS ONE. 2015;10(5):e0125628.

    Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Tanaka T. The taxonomy and nomenclature of Rutaceae Aurantioideae. Blumea. 1936;2:101–10.

    Google Scholar 

  25. 25.

    Mabberley DJ. Australian Citreae with notes on other Aurantioideae (Rutaceae). Telopea. 1998;7:333–44.

    Article  Google Scholar 

  26. 26.

    Olmstead RG, Palmer JD. Chloroplast DNA systematics: a review of methods and data analysis. Am J Bot. 1994;81:1205–24.

    CAS  Article  Google Scholar 

  27. 27.

    Samuel R, Ehrenddorfer F, Chase MW, Greger H. Phylogenetic analyses of Aurantoideae (Rutaceae) based on non-coding plastid DNA sequences and phytochemical features. Plant Biol. 2001;3:77–87.

    CAS  Article  Google Scholar 

  28. 28.

    Araujo EF, Queiroz LP, Machado MA. What is Citrus? Taxonomic implications from a study of cp-DNA evolution in the tribe Citreae Rutaceae subfamily Aurantioideae. Org Divers Evol. 2003;3:55–62.

    Article  Google Scholar 

  29. 29.

    Morton C, Grant M, Blackmore S. Phylogenetic relationships of the Aurantioideae inferred from chloroplast DNA sequence data. Am J Bot. 2003;90:1463–9.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Morton C. Phylogenetic relationships of the Aurantioideae Rutaceae based on the nuclear ribosomal DNA ITS region and three noncoding chloroplast DNA regions, atpB-rblc spacer, rps16, and trnL-trnF. Org Divers Evol. 2008;9:52–68.

    Article  Google Scholar 

  31. 31.

    Bayer R, Mabberley D, Morton CM, Cathy H, Sharma I, Pfeil BE, et al. A molecular phylogeny of the orange subfamily Rutaceae: Aurantioideae using nine cpDNA sequences. Am J Bot. 2009;96:668–85.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Lu Z, Zhou Z, Xie R. Molecular phylogeny of the true Citrus fruit trees group Aurantioideae, Rutaceae as inferred from chloroplast DNA sequences. Agr Sc China. 2011;10:49–57.

    CAS  Article  Google Scholar 

  33. 33.

    Bausher M, Singh N, Lee S, Jansen R, Daniell H. The complete chloroplast genome sequence of Citrus sinensis (L.) Osbeck var ‘Ridge Pineapple’: organization and phylogenetic relationships to other angiosperms. BMC Plant Biol. 2006;6:21.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Green RM, Vardi A, Galun E. The plastome of Citrus. Physical map, variation among Citrus cultivars and species and comparison with related genera. Theor Appl Genet. 1986;72:170–7.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Yamamoto M, Kobayashi S, Nakamura Y, Yamada Y. Phylogenic relationships of citrus revealed by RFLP analysis of mitochondrial and chloroplast DNA. Jpn J Breed. 1993;43:355–65.

    CAS  Article  Google Scholar 

  36. 36.

    Penjor T, Yamamoto M, Uehara M, Ide M, Matsumoto N, Matsumoto R, Nagano Y. Phylogenetic Relationships of Citrus and Its Relatives Based on MatK Gene Sequences. Plos One. 2013;8, e62574.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Li Y, Cheng Y, Tao N, Deng X. Phylogenetic Analysis of Mandarin Landraces, Wild Mandarins, and Related Species in China Using Nuclear LEAFY Second Intron and Plastid trnL-trnF Sequence. J Am Soc Hortic Sci. 2007;1326:796–806.

    Google Scholar 

  38. 38.

    Jena SN, Kumar S, Nair KN. Molecular phylogeny in Indian Citrus L. Rutaceae inferred through PCR-RFLP and trnL-trnF sequence data of chloroplast DNA. Sci Hortic-Amsterdam. 2009;199:403–16.

    Article  Google Scholar 

  39. 39.

    Yamamoto M, Tsuchimochi Y, Ninomiya T, Koga T, Kitajima A, Yamasaki A, et al. Diversity of Chloroplast DNA in Various Mandarins Citrus spp. and Other Citrus Demonstrated by CAPS Analysis. J Jpn Soc Hortic Sci. 2013;2:106–13.

    Article  Google Scholar 

  40. 40.

    Penjor T, Anai T, Nagano T, Matsumoto R, Yamamoto M. Phylogenetic relationships of Citrus and it relatives based on rbcL gene sequences. Tree Genet Genomes. 2010;6:931–9.

    Article  Google Scholar 

  41. 41.

    Li X, Xie R, Lu Z, Zhou Z. The Origin of Cultivated Citrus as Inferred from Internal Transcribed Spacer and Chloroplast DNA Sequence and Amplified Fragment Length Polymorphism Fingerprints. J Am Soc Hortic Sci. 2010;1354:341–50.

    Google Scholar 

  42. 42.

    Deng Z, La Malfa S, Xie Y, Xiong X, Gentile A. Identification and evolution of chloroplast uni- and trinucleotide sequence repeats in citrus. Sci Hortic-Amsterdam. 2007;111:186–92.

    CAS  Article  Google Scholar 

  43. 43.

    Garcia-Lor A, Luro F, Ollitrault P, Navarro L. Genetic diversity and population structure analysis of mandarin germplasm by nuclear, chloroplastic and mitochondrial markers. Tree Genet Genomes. 2015;11:123.

    Article  Google Scholar 

  44. 44.

    Hang L, Xiaoming Y, Lianshu Z, Hualin Y, Lijun C, Xiuxin D. Parentage analysis of natural Citrus hybrid ‘Zhelong Zhoupigan’ based on nuclear and chloroplast SSR markers. Sci Hortic-Amsterdam. 2015;186:24–30.

    Article  Google Scholar 

  45. 45.

    Curk F, Ollitrault F, Garcia-Lor A, Luro F, Navarro L, Ollitrault P. Phylogenetic origin of limes and lemons revealed by cytoplasmic and nuclear markers. Ann Bot-London. 2016;117(4):565–83.

    Article  Google Scholar 

  46. 46.

    Lotfy S, Luro F, Carreel F, Froelicher Y, Rist D, Ollitrault P. Application of cleaved amplified polymorphic sequence method for analysis of cytoplasmic genome among Aurantioideae intergeneric somatic hybrids. J Am Soc Hortic Sci. 2003;128:225–30.

    CAS  Google Scholar 

  47. 47.

    Su HJ, AHogenhout SA, Al-Sadi AM, Kuo CH. Complete chloroplast Genome Sequence of Omani Lime (Citrus aurantiifolia) and Comapative Analysis within the Rosids. Plos one. 2014;9:e113049.

    Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleid Acids Symp. 1999;41:95–8.

    CAS  Google Scholar 

  49. 49.

    Cuppen E. Genotyping by Allele-Specific Amplification (KASPar). CSH Protoc. 2007:pdb.prot4841

  50. 50.

    Nei M. Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci. 1973;70:3321–3.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Abkenar AA, Mardi M, Matsumoto R. Comparative analysis of organelle DNAs in citrus to clarify the origin of yuzu (citrus junos sieb. Ex tanaka) and its relationship with other acid citrus. Acta Hortic. 2011;892:6.

    Google Scholar 

  53. 53.

    Abkenar AA, Isshiki S, Tashino Y. Phylogenetic relationships in the true Citrus fruit trees revealed by PCR-RFLP analysis of cp DNA. Sci Hortic-Amsterdam. 2004;102:233–42.

    CAS  Article  Google Scholar 

  54. 54.

    Roose ML. Isozymes and DNA restriction fragment polymorphism in citrus breeding and systematics. Proc Int Soc Citricult. 1998;1:155–65.

    Google Scholar 

  55. 55.

    Scora RW. Biochemistry, taxonomy and evolution of modern cultivated Citrus. Proc Int Soc Citricult. 1988;1:277–89.

    Google Scholar 

  56. 56.

    Herrero R, Asìns MJ, Carbonell EA, Navarro L. Genetic diversity in the orange subfamily Aurantioideae. I. Intraspecies and intragenus genetic variability. Theor Appl Genet. 1996;92:599–609.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Ogawa K, Kawasakib A, Omurab M, Yoshidab T, Ikomab Y, Yanob M. 2001. 3′,5′-Di- C-β glucopyranosylphloretin, a flavonoid characteristic of the genus Fortunella. Phytochemistry. 2001;57:737–42.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Gulsen O, Roose ML. Chloroplast and nuclear genome analysis of the parentage of lemons. J Am Soc Hortic Sci. 2001;126:210–5.

    CAS  Google Scholar 

  59. 59.

    Ramadugu C, Pfeil BE, Keremane ML, Lee RF, Maureira-Butler IJ, Roose ML. A six nuclear gene phylogeny of Citrus (Rutaceae) taking into account hybridization and lineage sorting. PLoS ONE. 2013;8, e68410.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  60. 60.

    Pfeil BE, Crisp MD. The age and biogeography of Citrus and the orange subfamily (Rutaceae: Aurantioideae). Am J Bot. 2008;95:1621–31.

    Article  PubMed  Google Scholar 

  61. 61.

    Jung YH, Kwon HM, Kang SH, Kang JH, Kim SC. Investigation of the phylogenetic relationships within the genus Citrus Rutaceae and related species in Korea using plastid trnL-trnF sequences. Sci Hortic-Amsterdam. 2005;104:179–88.

    CAS  Article  Google Scholar 

  62. 62.

    Cheng Y, De Vicente MC, Meng H, Guo W, Tao N, Deng X. A set of primers for analyzing chloroplast DNA diversity in Citrus and related genera. Tree Physiol. 2005;25:661–72.

    Article  PubMed  Google Scholar 

  63. 63.

    Olivares-Fuster O, Hernandez-Garrido M, Guerri J, Navarro L. Plant somatic hybrid cytoplasmic DNA characterization by single-strand conformation polymorphism. Tree Physiol. 2007;27:785–92.

    CAS  Article  PubMed  Google Scholar 

  64. 64.

    Garg K, Green P, Nickerson DA. Identification of candidate coding region single nucleotide polymorphisms in 165 human genes using assembled expressed sequence tags. Genome Res. 1999;9:1087–92.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Batley J, Barker G, O’Sullivan H, Edwards KJ, Edwards D. Mining for single nucleotide polymorphisms and insertions/deletions in maize expressed sequence tag data. Plant Physiol. 2003;132:84–91.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  66. 66.

    Schroeder H, Fladung M. Differentiation of Populus species by chloroplast SNP markers for barcoding and breeding approaches. In: Vettori C, Fladung M, editors. Journal of Biogeosciences and Forestry. 2014. p. 544–6.

    Google Scholar 

  67. 67.

    Lee YG, Jeong N, Kim JH, Lee K, Kim KH, Pirani A, et al. Development, validation and genetic analysis of a large soybean SNP genotyping array. Plant J. 2015;81:625–36.

    CAS  Article  PubMed  Google Scholar 

  68. 68.

    Vos P, Uitdewilligen J, Voorrips R, Visser R, Van Eck H. Development and analysis of a 20 K SNP array for potato (Solanum tuberosum): an insight into the breeding history. Theor Appl Genet. 2015;128:2387–401.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  69. 69.

    Wu GA, Prochnik S, Jenkins J, Salse J, Hellsten U, Murat F, et al. Sequencing of diverse mandarin, pummelo and orange genomes reveals complex history of admixture during citrus domestication. Nat Biotechnol. 2014;32:656–62.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


We acknowledge the CRB Citrus (Inra/Cirad, France) and USDA-ARS National Clonal Germplasm Repository for Citrus & Dates (Riverside, USA) for plant material samples.


This work was supported by a grant (AGL2011-26490) from the Ministry of Economy and Innovation-Fondo Europeo de Desarrollo Regional (FEDER), Spain, and a grant from Tunisian “Ministère de l’enseignement supérieur et de la recherche scientifique » and ED STVST of El Manar University from Tunis, Tunisia.

Availability of data and materials

The datasets supporting the results of this article are included within the article (and its additional files 1 and 3).

Authors’ contributions

AO performed sequence selection, SNP identification and selection, primer design, genetic analysis of SNP data, molecular genetics studies and drafted the manuscript. FO participated in the molecular marker analysis. GB participated in molecular genetics studies. ASH participated in the coordination of the study and revised the manuscript. LN participated in the coordination of the study and provided several germplasm accessions. PO drew up and designed the study, participated in its coordination and data analysis, and drafted the manuscript. All the authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information



Corresponding author

Correspondence to Patrick Ollitrault.

Additional files

Additional file 1:

Sequence accession (GenBank) used for SNP in silico mining and their classification in the Bayer et al. [21] analysis. (XLSX 20 kb)

Additional file 2:

List of plant material used for KASPar analysis. (XLSX 143 kb)

Additional file 3:

List of the selected SNPs for KASPar marker development with the sequences of surrounding 50 bases each side and the Bayer et al. (2009) clade theoretically differentiated. (XLSX 13 kb)

Additional file 4:

Example of diagnostic SNP mining in the region rps4-trnT in the position 49301 distinguishing the accessions of clade Q (A) from all other Aurantioideae accessions (G). (XLSX 38 kb)

Additional file 5:

Sequence alignment of all Aurantioideae accessions for the eight chloroplast regions. (XLSX 3440 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Oueslati, A., Ollitrault, F., Baraket, G. et al. Towards a molecular taxonomic key of the Aurantioideae subfamily using chloroplastic SNP diagnostic markers of the main clades genotyped by competitive allele-specific PCR. BMC Genet 17, 118 (2016).

Download citation


  • Aurantioideae
  • Phylogeny
  • Clades
  • Chloroplastic SNP markers
  • KASPar™ genotyping