ORIGINAL RESEARCH article

Front. Plant Sci., 15 September 2025

Sec. Plant Breeding

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1645145

Harnessing chloroplast SSRs to decipher genetic diversity in underutilized Allium species

  • 1. Indian Council of Agricultural Research (ICAR)-Directorate of Onion and Garlic Research, Pune, Maharashtra, India

  • 2. Mahatma Phule Krishi Vidyapeeth, Rahuri, Maharashtra, India

  • 3. Indian Council of Agricultural Research (ICAR)-NBPGR Regional Station, Bhowali, Uttarakhand, India

  • 4. Indian Council of Agricultural Research (ICAR)-IARI Regional Station, Pune, Maharashtra, India

  • 5. Indian Council of Agricultural Research (ICAR)-Indian Agriculture Research Institute, Hazaribagh, Jharkhand, India

Abstract

Alliums, including vital crops such as onion, garlic, chives, bunching onion, and leek, are globally prized for their culinary applications and medicinal attributes. However, their genetic improvement remains constrained by large genome size, high heterozygosity, and limited characterization of genetic resources. To bridge this gap, we developed chloroplast simple sequence repeat (cp-SSR) markers, which are particularly suitable for population genetics studies because of their maternal inheritance, low recombination rates, and high variability. Leveraging the chloroplast genome of Allium fistulosum, we identified 22 cp-SSR loci, with tetranucleotides being the most prevalent, followed by di-, tri-, and pentanucleotides. Screening 96 underutilized Allium accessions using polymorphic cp-SSR markers revealed 89.2% polymorphism, indicating substantial genetic diversity. The polymorphism information content (PIC) ranged from 0.00 to 0.66 (average 0.20), confirming the utility of these markers in diversity assessments. The population structure analysis revealed three distinct genetic clusters, whereas phylogenetic analysis categorized the accessions into six major clades, mirroring their evolutionary divergence. Fixation index (FST) analysis showed high genetic differentiation (mean FST = 0.6) among accessions. These findings underscore the significance of cp-SSRs in revealing genetic structure and diversity across underutilized Allium species. This work lays a crucial foundation for integrating chloroplast markers with nuclear genomic and omics tools to drive the development of resilient, high-value cultivars suited to future agricultural challenges.

Introduction

The genus Allium, comprising approximately 750 underutilized species (), is the largest within the monocot group and includes numerous economically important plants. The key species include onion (Allium cepa L.), garlic (Allium sativum L.), chives (Allium schoenoprasum L.), leek (Allium porrum L.), and bunching onions (Allium fistulosum L.) (). In addition to these well-known crops, the genus also encompasses lesser-known species, such as A. altaicum, A. ramosum, A. chinense, and A. tuberosum etc. which are of interest for their potential contributions to biodiversity research. In addition to their agricultural value, Allium species hold significant ecological importance, with some listed in the Red Book of countries such as Mongolia, Russia, and China due to concerns over their conservation status ().

Underutilized Allium species exhibit significant chromosomal diversity, ranging from diploid to highly polyploid forms. This remarkable cytogenetic variability reflects the genus’ complex evolutionary history and dynamic genomic architecture, with origins tracing back to regions of Asia and Europe. Over time, these lesser-studied species have adapted to a wide array of ecological niches, resulting in a rich spectrum of phenotypic traits and specialized adaptations. Investigating their chromosomal profiles and evolutionary trajectories not only enhances our understanding of Allium genomics but also supports crop improvement efforts and biodiversity conservation. Notably, underexploited Allium populations represent untapped reservoirs of genetic diversity, offering valuable traits for sustainable breeding and long-term genetic resource management.

To date, the genetic diversity of onions and related Allium species has been examined using various molecular markers, including RAPDs (; ), ISSRs (; Sudha et al., 2019; , ), and combinations of RAPD and ISSR (Sudha et al., 2019), RAPD and PCR-RFLP (), and RAPD and SSR (). Additional marker systems such as RFLPs (), AFLPs (Van Heusden et al., 2000; Simó et al., 2014), TRAP (; ), SRAP combined with ISSR (), ILP (; ), and SRAP (; ) have also been widely applied.

More recently, SSR and SNP markers have been extensively employed to assess genetic variation in Allium fistulosum (Yamashita et al., 2010), Allium mongolicum (), and Allium cepa (; ; ; ; ; ; ; ), among others. Studies incorporating STS and SNPs (; Scholten et al., 2016; ; Villano et al., 2019; ; ; ; ; ; ; ) further highlight the growing utility of these high-resolution markers in understanding genetic variation. Among these, SSR markers have emerged as a preferred tool due to their high polymorphism, co-dominant inheritance, reproducibility, and cross-species transferability (Son et al., 2012; Wang and Zhang, 2022). These characteristics make SSRs effective for evaluating genetic diversity and population structure in plant species.

In addition to nuclear markers, the chloroplast genome characterized by maternal inheritance offers a powerful system for elucidating evolutionary relationships, phylogeography, and population genetics within and across Allium species. As a type of SSR marker, chloroplast simple sequence repeats (cp-SSRs) are particularly advantageous due to their high mutability, conservation, variability, co-dominant inheritance, and organelle-specific transmission (). In Allium, cp-SSR markers have successfully revealed genetic diversity, population differentiation, and evolutionary relationships among closely related taxa (). They have also been used to detect historical demographic events, such as bottlenecks and genetic drift, which are critical for understanding population dynamics (). Strikingly, cp-SSRs contributed to onion breeding programs by facilitating targeted genetic analyses and facilitating the development of conservation strategies (Sharma et al., 2020).

Cp-SSR markers offer a non-destructive and efficient means for detecting subtle genetic variations in the chloroplast genome (). Only a few studies have reported the use of cp-SSR markers in Allium species, such as Allium cepa L., Allium sativum L., and Allium paradoxum (M. Bieb.) (). The present study aims to investigate chloroplast genetic divergence, heterozygosity, allelic diversity, population structure, and genetic relatedness across 96 underutilized Allium species using 22 cp-SSR markers. By doing so, this research seeks to provide a foundation for the strategic use of these genotypes in future breeding programs and to guide the conservation of these critical plant resources in the face of ongoing environmental challenges.

Materials and methods

Plant material and DNA extraction

In the present study, a total of 96 underutilized Allium (A.) accessions (Table 1) were randomly collected from their primary regions of distribution across India. The collected samples were planted during the regular crop growing season at an experimental field site at the Indian Council of Agricultural Research-Directorate of Onion and Garlic Research (ICAR-DOGR) in Rajgurunagar, Pune, Maharashtra, India, which is located at geographic locations (18°52’0”N, 73°54’0”E; 645 m above sea level). Young leaf tissues from ten individuals per accession (96 underutilized Allium species) were collected randomly for genomic DNA isolation. Total genomic DNA was isolated from these samples via the CTAB method as described by . Leaf tissues were homogenized in liquid nitrogen and incubated for an hour at 65°C in 1 ml of CTAB buffer, which contained 4% polyvinylpyrrolidone (PVP), 0.5% β-mercaptoethanol, 1.4 M NaCl, 100 mM Tris-HCl, 20 mM EDTA, and 2% cetyl trimethylammonium bromide. The quantity and quality of the extracted DNA were assessed by electrophoresis on a 0.8% agarose gel, using lambda HindIII marker (Thermo Fisher Scientific) used as a reference.

Table 1

Sl. no.Name of the speciesCodeSl. no.Sample_nameCodeSl. no.Sample_nameCode
1Alliun cepa L.
Var. Bhima Super
AceB.Super33Allium macranthum L. NMK 3233AmacNMK323365Allium hookeri L.AhooNG3155
2Allium sativum L.
Var. Bhima Purple
AsaB.Purple34Allium macranthum L. NMK 3232AmacNMK323266Allium hookeri L.AhooNMK3235
3Allium altaicum L.
EC 328485 Pall
AalEC328485P35Allium macranthum L. NMK 3229AmacNMK322967Allium hookeri L.AhooNG3156
4Allium altaicum L.
Pall CGN 14769
AalCGN14769P36Allium macranthum L.Amacranthum68Allium hookeri L.A.hookeri
5Allium cepa var aggregatum
3 Meitai Tilou
Aceaggr3MT37Allium tuberosum L.AtubBKCGN1574969Allium auriculatum L.A.auriculatum
6Allium cepa var aggregatum
4 Eshing Tilou
Aceaggr4ET38Allium tuberosum L.AtubCGN1641870Allium albidum L.A.albidum
7Allium cepa var aggregatum
5 Manipur
Aceaggr5M39Allium tuberosum L.AtubESKCGN1641271Allium oreoprasum L.A.oreoprasum
8Allium cepa × Allium fistulosum Beltsville BunchingAce × AfisB40Allium tuberosum L.AtubKazhakistan158772Allium ramosum L.A.ramosum
9Allium fistulosum L. NIC 23426AAfisNIC23426A41Allium tuberosum L.AtubZimmu73Allium fasciculatum L.A.fasciculatum
10Allium fistulosum L. NIC 20221AfisNIC2022142Allium tuberosum L.AtubMKG2474Allium viviparum L.A.viviparum
11Allium fistulosum L. EC 321643.1AfisEC321643.143Allium tuberosum L.AtubMKG321475Allium stracheyi L.A.stracheyi
12Allium fistulosum L. EC 321643.2AfisEC321643.244Allium tuberosum L.AtubNMK321476Allium negianum L.A.negianum
13Allium fistulosum L. AKO-1 (China)AfisAKO-1.China45Allium tuberosum L.AtubIC35352477Allium consanquianum L.A.consanguineum
14Allium fistulosum L.
BG Autumn Nepthane
AfiBGAN46Allium tuberosum L.AtubNG318378Allium roylei L.A.roylei
15Allium fistulosum L.A.fistulosum47Allium tuberosum L.AtubNMK321979Allium proliferum L.A.proliferum
16Allium fistulosum L. GeorgianAfisGeorgien48Allium tuberosum L.Atub.OP80Allium wallichii L.A.wallichii
17Allium fistulosum L. (China) All 647AfisAll647.China49Allium tuberosum L.AtubNMK323181Allium barsczewskii L.AbarMK95
18Allium macranthum L. NMK 3240AmacNMK324050Allium tuberosum L.AtubNMK322882Allium senescense L.AsenEC328503
19Allium chinense
Cholang White RAK100
AChiRAK10051Allium tuberosum L.AtubNMK320783Allium schoenoprasum L.AschNMK12
20Allium chinense L. NMK 3247AChiNMK324752Allium tuberosum L.AtubEC60748384Allium schoenoprasum L.Aschoenoprasum
21Allium chinense L.A.chinense53Allium tuberosum L.Atuberosum85Allium schoenoprasum L.AschhNR6NGB14774
22Allium macranthum L. NMK 3216AmacNMK321654Allium prszewalskianum L.AprzMMK12086Allium schoenoprasum L.AschNRNGB597
23Allium macranthum L. NMK 3242AmacNMK324255Allium prszewalskianum L.AprzMMK11987Allium ascalonicum L.AascalC353523
24Allium macranthum L. NMK 3248AmacNMK324856Allium prszewalskianum L.AprzMMK12188Allium ascalonicum L.AascalC99923
25Allium macranthum L. NMK 125AmacNMK12557Allium chinense L.AChiNMK316589Allium porrum L.AporlC632238
26Allium macranthum L. NMK 3227AmacNMK322758Allium chinense L.AchiNMK323690Allium porrum L.Aporl353526
27Allium macranthum L. NMK 3246AmacNMK324659Allium chinenes L.AchiMMK13191Allium altaicum L.AalCGN23934P
28Allium macranthum L. NMK 3245AmacNMK324560Allium ascalonicum L.A.ascalonicum92Allium altaicum
CGN 14769
AalCGN14769P
29Allium macranthum L. NMK 3244AmacNMK324461Allium fragrance
EC 383446
AfraEC38344693Allium altaicum L.
CGN 14771
AalCGN14771
30Allium macranthum L. NMK 3243AmacNMK324362Allium angulosum
EC 328486
AangEC32848694Allium fragrance L.AfraEC383447
31Allium macranthum L. NMK 3238AmacNMK323863Allium ampeloprasum NMK 3211AampNMK321195Allium fistulosum L.Afis.OP
32Allium macranthum L. NMK 3237AmacNMK323764Allium carolinianum L.AcarMMK13596Allium schoenoprasum L. NR6 NGB 147745AschNR6NGB14775

Description of the samples used for the characterization of chloroplast derived simple sequence repeat markers (cp-SSR).

Chloroplast SSR marker development

The simple sequence repeats (SSR) loci within the Allium fistulosum chloroplast genome (Voucher No. PRJNA927338; NCBI Reference ID: NC_040222.1; ) were identified using the MISA tool (http://misaweb.ipk-gatersleben.de/). The cp-SSR motif analyzed consisted of repeat units ranging from di- to hexanucleotides, meeting the minimum repeat thresholds set by MISA. Specifically, six motifs contained dinucleotide repeats, four contained trinucleotides, and three included tetra-, penta-, or hexanucleotide repeats. Mononucleotide repeats were excluded from further analysis. Both perfect and compound SSRs were detected via the MISA pipeline, with compound repeats defined as SSRs interrupted by non-repeat sequences of up to 100 bp. Primer pairs flanking the cp-SSR loci were designed using the BatchPrimer3 v1.0 online tool (https://probes.pw.usda.gov/batchprimer3; You et al., 2008). The primer design parameters included primer lengths of 22–27 nucleotides, amplicon sizes of 100 to 300 bp, melting temperature ranging from 48°C to 55°C, and GC content between 40% and 70%, with an optimal GC content of 50% (Table 2).

Table 2

cp-SSR IDPrimer nameSSR typeSSRSizeStartEndPosition
Afi01cp-SSR01p5(TAAAA)31537633777trnK-UUU
Afi02cp-SSR02p4(ATAA)31242954306trnK- rps16
Afi03cp-SSR03p3(TTA)41271427153psbK- psbI
Afi04cp-SSR04p4(ATTT)3121981019821rpoC1-rpoC2
Afi05cp-SSR05p5(ATTGA)3152902329037pet9- psbM
Afi06cp-SSR06p2(AT)6123197931990trnT-GGU- psbD
Afi07cp-SSR07p2(TA)6123510135112trnS-UGA- psbZ
Afi08cp-SSR08p4(TTTC)3124278942800ycf3 Intron
Afi09cp-SSR09p2(TA)8164556845583rps4- trnT-UGU
Afi10cp-SSR10p5(TATAA)3155819858212accD- psaI
Afi11cp-SSR11p4(AATG)3126027860289cemA
Afi12cp-SSR12p4(TAAA)3126189161902petA- psbJ
Afi13cp-SSR13p2(TA)6126575565766psaJ- rpl33
Afi14cp-SSR14p3(TCT)4126900069011ClpP Intron
Afi15cp-SSR15p4(TAAA)3126983169842ClpP Intron
Afi16cp-SSR16p4(TTTA)3127268172692psbN- psbH
Afi17cp-SSR17p4(GGAT)3127475974770petB- PetD
Afi18cp-SSR18p4(ATTG)312111468111479ndhF- rpl32
Afi19cp-SSR19p3(AAT)412115962115973psaC- ndhE
Afi20cp-SSR20p2(TA)714116305116318psaC- ndhE

Distributions of cp-SSR motifs observed in the chloroplast genome of Allium fistulosum.

p2, Dinucleotide; p3, Trinucleotide; p4, Tetranucleotide; p5, Pentanucleotide.

cp-SSR marker analysis

A total of 22 cp-SSR primer pairs were selected and synthesized by Eurofins Genomics (Eurofins, India). After an initial run with the newly developed primer pairs, 20 cp-SSRs exhibiting high resolution, stability, and significant polymorphism were selected for further analysis. The cp-SSR amplification was carried out in a 20 µl reaction volume, which included 2 µl of 10X reaction buffer and 50 ng of template DNA per reaction. The PCR reaction mixture consisted of 50 ng of genomic DNA (1µl), 1.5 mM MgCl2, 0.2 mM of each dNTP, 0.2 µM of each primer (forward and reverse), and 5 U of Taq DNA polymerase. PCR amplification was performed using a Bio-Rad iCycler thermal cycler. The cycling conditions included an initial denaturation at 94°C for 4 minutes, followed by 35 cycles of denaturation at 94°C for 1 minute, annealing at the optimized temperature specific to each primer (as listed in Table 3), and extension at 72°C for 40 seconds. A final extension was carried out at 72°C for 10 minutes. The PCR products were analyzed via gel electrophoresis on a 3.2% agarose gel. Bands were visualized with a 1 kb Plus DNA ladder (Thermo Fisher Scientific) as a reference and documented using a gel documentation system.

Table 3

Sl no.Primer namePrimer F/RPrimer sequence 5’ to 3’No. of basesTm (°C)Allele size range
1cp-SSR01cp-SSR01-FCCAGAATTAGAGCCGTAGAGC2150.22160-240
2cp-SSR01-RCCACGACTGATCCTGAAAGG20
3cp-SSR02cp-SSR02-FTGGCAAACCCATAATTTGAA2048.2200
4cp-SSR02-RTGTGCCAATCCAACACAAAT20
5cp-SSR03cp-SSR03-FTCCTCGTTCTGACCTTCCAG2049.8200
6cp-SSR03-RTGTTGACATAGTGCCCCAAA20
7cp-SSR04cp-SSR04-FATAAACCCGACTTCCCAAGG2052.280-110
8cp-SSR04-RGAAGCCATACAGGGGTTTTG20
9cp-SSR05cp-SSR05-FTCAGCGCAATCATTTCATTT2055.1200-220
10cp-SSR05-RTCGCACTTATTGCTACTGCAC21
11cp-SSR06cp-SSR06-FTGATTTTCTTGTTAATGGACGC2253.4200-220
12cp-SSR06-RTGCATTGCTGAAACAAAACAA21
13cp-SSR07cp-SSR07-FTGTAGAAACCTCCCGGATTG2054.2100
14cp-SSR07-RATTCGGACATGGAGTCGAAG20
15cp-SSR08cp-SSR08-FATCGTTGCTTTGAACGATGC2055.1200-210
16cp-SSR08-RTATTTCCGGGCATTAGAACG20
17cp-SSR09cp-SSR09-FAAACAAAGCAAAGCGAAATCT2152.3180-220
18cp-SSR09-RCCATTTCTACAAACGTTGAGTCAC24
19cp-SSR10cp-SSR10-FTGGGTTGTCATACATATTCGTG2254.2200-260
20cp-SSR10-RTGTCATAGAACGGGTACCTCAA22
21cp-SSR11cp-SSR11-FATCGCGTATCTCCTTCGCT1951.2200
22cp-SSR11-RCCTATCCACGAGTCTGCCAT20
23cp-SSR12cp-SSR12-FTGCTTTTTCTCTTGTTCACCA2150.6100
24cp-SSR12-RCTCAATGAATGACTCCCCTCT21
25cp-SSR13cp-SSR13-FAAATGAAATACTGGAAAGAATAATTGA2751.4220- 240
26cp-SSR13-RACCCTTAGCCATGAACCTCC20
27cp-SSR14cp-SSR14-FATTCAATATGGCGAAGGCAT2051.4220
28cp-SSR14-RGATCCTTCATTCTGGTCGGA20
29cp-SSR15cp-SSR15-FTCAATTCGTTTCATGTCTCCA2152.4200
30cp-SSR15-RTGGAGTATCCAGGCTCTGCT20
31cp-SSR16cp-SSR16-FATGGCGACTAAGGTTGCTGT2055.2180
32cp-SSR16-RCTCAACGGTTTGTGTAGCCA20
33cp-SSR17cp-SSR17-FATCCTATCGGGAAGGAACAA2054.3180
34cp-SSR17-RGCATGGCCCAATCAATAGTT20
35cp-SSR18cp-SSR18-FGCAAAGAAAAAGTAAGAAAAGTAAGCA2753.880 -280
36cp-SSR18-RTTCCGTTAGCTAAGAAAAGGAACT24
37cp-SSR19cp-SSR19-FGCAGGCTCGTACACATTGAG2053.5100
38cp-SSR19-RTTGGTCCCTTCTGATGAACA20
39cp-SSR20cp-SSR20-FTCACGTTGACTAATGGATCGTC2254.7120- 260
40cp-SSR20-RTGGGCTAGCTATTGTTTCGTC21

Details of 20 chloroplast SSR markers, including sequences, annealing temperatures, and allele size.

Scoring and data analysis

To ensure the accuracy of the results, each pair of primers was used for PCR amplification and electrophoresis twice, and only the cp-SSR markers with high definition and good stability were scored. The scoring was performed on the basis of the absence (0) or presence (1) of each band for all the isolates in each primer. Genetic variation at each locus was characterized in terms of the number of alleles, and the PIC value was calculated. The binary matrix was subjected to Jaccard similarity coefficient analysis using NTSYS-pc version 2.02i (), and the unweighted pair group method with arithmetic mean (UPGMA) clustering map based on Nei’s genetic distance was constructed using MEGA X (). Principal component analysis (PCA) and a cluster matrix were plotted on the basis of correlation distance and average genetic linkage via the web tool Clust-Vis (https://biit.cs.ut.ee/clustvis/) (). Genetic diversity parameters, including minor allele frequency (MAF), observed number of alleles (Na), observed heterozygosity (Ho), expected heterozygosity (He), and PIC, were calculated by GenAlEx 6.51 software (). GenAlEx 6.51 was used to calculate the fixation index (FST), which measures the proportional increase in homozygosity. FST values range from 0 (no differentiation) to 1 (complete differentiation) (Wright, 1984). The population genetic structure was analyzed using Bayesian clustering methods via STRUCTURE 2.3.4 software. The number of populations (K) was tested sequentially from 1-10. Each run included a burn-in phase of 50,000 steps, followed by 200,000 Markov chain Monte Carlo (MCMC) iterations, which enhanced the reliability of clustering, as suggested by . The optimal K value was determined via the average lnP(K) and StructureSelector (https://lmme.ac.cn/StructureSelector/), revealing a significant peak in the ΔK values for the most suitable population grouping.

Results

Characterization of the developed cp-SSR markers

In the current study, a total of 22 cp-SSR marker pairs were identified from the Allium fistulosum chloroplast genome. Among the identified markers, tetranucleotide motifs were the most abundant (45.45%), followed by dinucleotide (27.27%), trinucleotide (13.63%), and pentanucleotide (13.63%) motifs (Table 2 and Figure 1). Notably, hexanucleotide repeats were completely absent from the chloroplast genome of A. fistulosum. The most frequently occurring motif was TA (22.73%), followed by the TAAA motif (9.09%). All other motifs were evenly distributed. Mononucleotide repeats, primarily A/T-rich, were excluded from further analysis due to their low polymorphic potential and the risk of sequencing errors caused by homopolymer runs. The average repeat lengths for the di-, tri-, tetra-, and pentanucleotide cp-SSRs were 13, 12, 12, and 15 base pairs, respectively. Among the 22 designed primer pairs, 20 (90.91%) successfully amplified clear and reproducible bands during PCR screening with underutilized Allium species (Table 3). These markers showed high polymorphism and stability indicating their potential suitability for assessing genetic diversity in Allium germplasm.

Figure 1

cp-SSR marker analysis

The analysis of the developed cp-SSR markers revealed a high level of polymorphism and genetic diversity among the underutilized Allium species (Supplementary Figure 1). A total of 37 allelic bands were detected using 20 polymorphic cp-SSR markers. Among these, 89.2% of the amplified alleles were polymorphic, indicating the hypervariable nature of the cp-SSR loci and the broad genetic variation in the underutilized Alliums. The observed allele sizes ranged from 80 to 280 bp (Table 3), with allele frequencies ranging from 0.00 to 0.99, demonstrating the effectiveness of these markers in capturing intra- and interspecific genetic variation. The polymorphism information content (PIC) values of the cp-SSR markers ranged from 0.00 to 0.66, with an average of 0.20 (Table 4). Notably, most of the markers exhibited more than 80% polymorphism, indicating their high utility for diversity and population genetic studies. These highly informative loci can serve as valuable molecular tools in future genetic analyses of Allium species.

Table 4

Primer nameAllele sizeNo. of bandsTotal no. of bands presentFreq AFreq BPICAvg. PIC%polymorphism
cp-SSR01160520.0210.9790.0410.04821.250
18010.0100.9900.021
200930.9690.0310.061
22050.0520.9480.099
24010.0100.9900.021
cp-SSR022001410.4270.5730.4890.48942.708
cp-SSR032001940.9790.0210.0410.04197.917
cp-SSR04802750.7810.2190.3420.28881.771
110130.1350.8650.234
cp-SSR052002920.9580.0420.0800.06093.042
22020.0210.9790.041
cp-SSR062002560.5830.4170.4860.28329.167
22040.0420.9580.080
cp-SSR07100110.0001.0000.0000.000
cp-SSR082002890.9270.0730.1350.08890.042
21020.0210.9790.041
cp-SSR091803110.1150.8850.2030.19191.389
200790.8230.1770.291
22040.0420.9580.080
cp-SSR102004360.3750.6250.4690.25362.748
220170.1770.8230.291
24080.0830.9170.153
26050.0520.9480.099
cp-SSR112001890.9270.0730.1350.13592.708
cp-SSR121001840.8750.1250.2190.21987.500
cp-SSR132202950.9900.0100.0210.02195.521
24010.0100.9900.021
cp-SSR142201870.9060.0940.1700.17090.625
cp-SSR15200190.0940.9060.1700.1709.375
cp-SSR161801800.8330.1670.2780.27783.333
cp-SSR171801700.7290.2710.3950.39572.917
cp-SSR1880360.0630.9380.1170.14713.250
26010.0100.9900.021
280180.1880.8130.305
cp-SSR19160100.0001.0000.0000.0000.000
cp-SSR202603750.7810.2190.3420.66095.389
240190.1980.8020.317
12040.0420.9580.080
Total3713680.196

Statistical analyses of genetic polymorphisms in 20 cp-SSR primer pairs.

Data scoring and analysis

The genetic diversity analysis of 20 cp-SSR markers across 96 underutilized Allium species along and their respective accessions provided a comprehensive view of the genetic structure and variability within the genus. The key diversity indices assessed included the number of observed alleles (Na), effective number of alleles (Ne), Shannon’s information index (I), observed heterozygosity (Ho), expected heterozygosity (He), and unbiased expected heterozygosity (uHe) (Table 5). The number of alleles (Na) ranged from 1 to 4, with markers such as SSR09 and SSR10 showing the highest diversity. The effective number of alleles (Ne) ranged from 1.00 to 2.488. Shannon’s index (I) varied from 0.0 for nonpolymorphic markers to 1.1 for SSR10, indicating high intra-accession diversity for that marker.

Table 5

Primer nameNNaNeIHoHeuHeFST
cp-SSR019331.0670.1570.0650.0630.0630.0
cp-SSR024111.0000.0000.0000.0000.000
cp-SSR039411.0000.0000.0000.0000.000
cp-SSR048821.3370.4190.0000.2520.2531.0
cp-SSR059421.0430.1030.0000.0420.0421.0
cp-SSR066021.1230.2220.0170.1100.1110.8
cp-SSR088921.0230.0620.0220.0220.0220.0
cp-SSR099441.6060.7570.0000.3770.3791.0
cp-SSR106342.4881.0970.0480.5980.6030.9
cp-SSR 118911.0000.0000.0000.0000.000
cp-SSR128411.0000.0000.0000.0000.000
cp-SSR 139521.0110.0330.0110.0100.011-0.1
cp-SSR 148711.0000.0000.0000.0000.000
cp-SSR15911.0000.0000.0000.0000.000
cp-SSR168011.0000.0000.0000.0000.000
cp-SSR177011.0000.0000.0000.0000.000
cp-SSR182431.7070.7000.0420.4140.4230.9
cp-SSR207721.3080.3980.2210.2360.2370.1

Genetic diversity analysis of cp-SSR markers across 96 wild Allium species.

The observed heterozygosity (Ho) was predominantly 0.000 for most of the markers, suggesting low heterozygosity levels, whereas cp-SSR10 exhibited a slightly greater value of 0.048. The expected heterozygosity (He) ranged from 0.000 to 0.598, and the unbiased expected heterozygosity (uHe) ranged from 0.000 to 0.603, both of which were highest for cp-SSR10. Markers such as cp-SSR09, cp-SSR10, and cp-SSR18 exhibited high genetic diversity, whereas cp-SSR02, cp-SSR03, cp-SSR11, and cp-SSR15 were monomorphic with no diversity (Na and Ne = 1.000; I, Ho, and He = 0.000). The fixation index (FST) ranged from 0.0 to 1.0 among the accessions within the six major clades, with a mean FST value of 0.6. Cross-transferability analysis among 30 Allium genotypes (Table 6) revealed that the transferability percentage of cp-SSR alleles ranged from 55% to 90%. The highest transferability was observed in A. viviparum (90%), whereas A. altaicum Pall CGN 14769 and A. altaicum CGN 1477 presented the lowest transferability (55%). Other important species, such as A. hookeri and A. fragrance, presented intermediate transfer percentages of 73.75% and 72.5%, respectively.

Table 6

Sl no.Name of the genotypeTransferable alleles% of transferability
1A. cepa var. Bhima Super12.0060.00
2A. sativum var. Bhima Purple14.0070.00
3A. altaicum EC 328485 Pall11.557.50
4A. cepa var aggr 3 Meitai Tilou12.0060.00
5A. fistulosum12.4562.27
6A. chinensis NMK 323613.3366.66
7A. macranthum NMK 32412.5662.81
8A. tuberosum Bawang Kuchaai CGN 1574912.6163.05
9A. prszewalskianum MMK 1214.0070.00
10A. fragrance EC 38344614.5072.50
11A. angulosum EC 32848613.6668.33
12A. hookeri NG 315514.7573.75
13A. auriculatum14.0070.00
14A. albidium16.0080.00
15A. oreoprasum13.0065.00
16A. ramosum15.0075.00
17A. fasciculatum15.0075.00
18A. viviparum18.0090.00
19A. stracheyi16.0080.00
20A. negianum17.0085.00
21A. consanguineum16.0080.00
22A. roylei15.0075.00
23A. proliferum16.0080.00
24A. barsczewskii13.0065.00
25A. senescens EC 3285-313.0065.00
26A. ascalonicum MMK 1313.0065.00
27A. porrum13.0065.00
28A. altaicum pall CGN 1476911.0055.00
29A. altaicum CGN 147711.0055.00
30A. schoenoprasum NR 6 NGB 1477412.0060.00

Cross-transferability observed in different underutilized Alliums.

Genetic relationships among underutilized and cultivated Alliums

Chloroplast microsatellite markers were utilized to assess genetic relationships among 96 underutilized Allium species through neighbor-joining (NJ) cluster analysis. The dendrogram (Figure 2) grouped the accessions into six distinct clusters (I–VI), each representing varying degrees of genetic relatedness. Cluster I was the largest, consisting of 38 accessions primarily representing Allium tuberosum and closely related taxa such as A. ramosum, A. viviparum, A. wallichi, A. consanguineum, A. albidium, and A. auriculatum. These accessions presented have high genetic similarity, likely due to shared ancestry, ecological adaptation, and geographical proximity. Cluster II included eight underutilized accessions, such as the Allium chinense MMK131, Allium chinense NMK3247, Allium chinense (RAK100), Allium chinense NMK3165, A. chinense, A. chinense NMK3236 and A. fragrans genotypes, which form a genetically cohesive group on the basis of habitat and cytoplasmic traits. Cluster III comprised 23 accessions dominated by A. fistulosum, along with related species such as A. ascalonicum, A. stracheyi, A. porrum, A. roylei etc. reflecting significant diversity and wide geographical origins. Cluster IV consisted of six cultivated Allium accessions, including landraces and hybrids (Allium cepa var. aggregatum 3, 4 and 5, Allium cepa var. Bhima Super, Allium sativum var. Bhima Purple, and Allium cepa × Allium fistulosum Beltsville Bunching, showing limited diversity due to breeding bottlenecks. Cluster V included 17 underutilized accessions, such as A. macranthum and its relatives, as well as A. porrum adapted to high-altitude environments, reflecting substantial genetic divergence. Cluster VI was the smallest, comprising four A. hookeri accessions showing a distinct genetic lineage. Notably, Clusters III and V presented the highest levels of intracluster genetic diversity, whereas Clusters I and IV were relatively homogeneous.

Figure 2

Population structure analysis

The population structure was analyzed using STRUCTURE v2.3.4 software based on data from 20 cp-SSR markers. The analysis was conducted for K-values ranging from 1 to 10. As K increased, the log probability of the data [lnP(K)] also increased (Figure 3A), and the optimal number of clusters was determined using the ΔK method (). A clear peak at K = 3 was observed (Figure 3B), indicating the most likely number of genetic clusters. Accordingly, the accessions were grouped into three distinct sub-populations: pop1, pop2, and pop3 (Figure 3C). The mean intracluster genetic distances for these three populations were 0.2208, 0.1664, and 0.1699, respectively, whereas the average allele-frequency divergence among populations was 0.1100. The alpha mean value was 0.051, and the proportion of membership for each cluster was estimated at 0.313, 0.185, and 0.502, respectively. These results reveal a moderate level of genetic structure and highlight substantial within-population diversity among underutilized Allium species.

Figure 3

Principal component analysis and heatmap analysis

PCA analysis revealed that the first two principal components, PC1 and PC2, accounted for 13.7% and 9.5% of the total genetic variation, respectively (Figure 4). The distribution of accessions across the PCA biplot indicated the presence of two major genetic clusters. The first cluster displayed broader dispersion, primarily in the negative PC1 axis, suggesting greater genetic diversity. In contrast, the second cluster appeared more compact and tightly grouped, indicating higher genetic similarity among its members. The overlap observed between the two clusters suggested the presence of genetic admixture among certain accessions. In accordance with the PCA results, the heatmap analysis (Figure 5) presented a color-coded visualization of pairwise genetic similarity, with red shading indicating high similarity (values closer to 1.0) and gray shading denoting low similarity (values closer to -1.0). The accompanying hierarchical clustering dendrogram revealed distinct clusters of genetically similar accessions, whereas vertical patterns across the heatmap highlighted conserved genetic markers.

Figure 4

Figure 5

Discussion

Chloroplast SSR markers have proven to be highly informative in studies of plant genetic diversity, phylogenetics, and population structure because of their uniparental inheritance, low recombination rates, and conserved genomic context. The central focus of this study was the development and characterization of chloroplast simple sequence repeat (cp-SSR) markers in Allium fistulosum, with subsequent amplification testing across 96 underutilized Allium species. A total of 22 novel cp-SSR loci were identified and validated, revealing distinctive patterns in repeat motif distribution that contribute to our understanding of cp-SSR evolution and utility within the Allium genus. Our findings revealed that tetranucleotide repeats were the most abundant (45.45%), followed by dinucleotide repeats (27.27%). Notably, mononucleotide repeats were deliberately excluded due to their lower informativeness and susceptibility to sequencing errors and polymerase slippage, a strategy aligned with best practices in cp-SSR marker development (). This exclusion also avoided overrepresentation of polyA/polyT stretches, which are highly abundant but offer limited polymorphic information. Furthermore, hexanucleotide motifs were absent, consistent with previous cp-SSR studies where longer motifs are generally rare in chloroplast genomes.

Our results partly align with prior studies in Allium species. For instance, reported the development of 22 cp-SSRs in Allium cepa and Allium sativum, and 15 cp-SSRs in Allium paradoxum. In contrast to our current findings, their work found that tri-nucleotide repeats were the most frequent motif type (50%), suggesting potential interspecific variation in repeat motif composition within the genus Allium. This difference could be attributed to genomic structural variation or differing evolutionary pressures across species. When compared to broader angiosperm studies, our cp-SSR motif composition is consistent with patterns observed in other taxa. For example, identified 139 cp-SSR loci across 11 tree peony plastomes, while Shukla et al. (2018) reported 21–25 cp-SSRs in various Vigna species (V. angularis, V. radiata, and V. unguiculata), and found 92 SSR loci across six Cupressaceae plastomes. In these studies, di- and tetranucleotide repeats also predominated, underscoring a conserved pattern of SSR distribution in chloroplast genomes across plant lineages.

Interestingly, our findings contrast with those of in Physalis angulata, where mononucleotide repeats were the most abundant (68.24%), followed by tetranucleotides (12.28%). The abundance of mononucleotide motifs in that study likely reflects a different analytical approach that included these repeats, which, while common, are typically avoided in marker development due to their lower polymorphic potential. Overall, the distribution of cp-SSR motif types in Allium fistulosum reflects both conserved and species-specific patterns observed across plant taxa. Our deliberate methodological choices such as excluding mononucleotide repeats support the development of highly informative, polymorphic, and stable markers, which are essential for downstream applications such as population genetics, phylogenetic reconstruction, and genetic diversity studies in Allium and related genera. This work contributes to the growing genomic toolkit for Allium research and supports future efforts in conservation and breeding of underutilized species. The high success rate of amplification (90.91%) and clear electrophoretic profiles of these cp-SSR markers demonstrate their robustness and reliability for genetic analysis. Similar success has been reported in other species, such as tree peonies, where 19 out of 21 cp-SSR markers amplified strongly (). The cp-SSR markers developed here complement existing nuclear SSRs by capturing maternal lineage information, thereby enhancing the resolution of genetic diversity studies in Allium species. Moreover, these markers are expected to support broader applications in phylogeography, conservation genetics, and breeding. As more chloroplast genome sequences become available, the cross-transferability and expansion of cp-SSR marker sets will continue to facilitate species-specific and cross-species analyses (; ). This study lays the groundwork for future genomic research and supports the strategic use of cp-SSR markers in Allium crop improvement and biodiversity assessment programs.

The cp-SSR marker analysis confirmed the effectiveness of the developed markers in revealing polymorphisms and genetic variation within Allium germplasm. The high rate of polymorphism (89.2%) and broad allele size range reflect the utility of these markers for studying genetic structure and relationships, especially in underutilized populations. The average PIC value (0.20) aligns with previous findings by , who reported PIC values ranging from 0.007 to 0.427 in Allium germplasm via chloroplast-derived SSR markers. Although slightly lower than in studies by (PIC: 0.24–0.98; avg. 0.608), (avg. 0.64), and (up to 0.7), the moderate PIC values in this study may be attributed to differences in genome source (chloroplast vs. nuclear SSRs), marker selection criteria, and the genetic backgrounds of the tested accessions.

Similar studies have reported varying PIC values, with and reporting averages of 0.45 and 0.51, respectively. reported 53 alleles using 145 SSR markers, with PIC values ranging from 0.219 to 0.715 and an average of 3.54 alleles per locus. reported PIC values ranging from 0.00 to 0.89 with 60 primers, detecting 54 alleles across 19 primers, with an average of 2.84 alleles per locus. These comparisons highlight the impact of population structure, genomic origin, and SSR motif type on marker informativeness. Highly polymorphic markers such as cp-SSR3 and cp-SSR14, with more than 80% polymorphism, offer strong potential for use in genetic mapping and diversity studies. Similar findings by , where 92 polymorphic loci were identified using 19 SSR markers across 95 Allium accessions, underscore the importance of SSRs in understanding genetic variation.

further demonstrated the utility of chloroplast-derived SSR markers in A. cepa and A. sativum, reporting heterozygosity values ranging from 0.009 to 0.540 and PIC values ranging from 0.007 to 0.427. These markers offer a reliable platform for evaluating genetic relationships between underutilized and cultivated Allium species. The observed genetic diversity in underutilized Allium accessions holds critical value for crop improvement. Traits such as disease resistance, yield enhancement, and abiotic stress tolerance can be introgressed from underutilized relatives into cultivated backgrounds. Therefore, the conservation and characterization of underutilized Allium germplasm remains essential for the resilience and sustainability of breeding programs. Recent advances underscore the complementary role of cp-SSRs markers alongside genomic tools in exploring Allium genetic diversity and evolutionary history (; ; ). The integration of cp-SSR data with nuclear SSR and genome-wide SNP dataset will further enrich our understanding of the genetic makeup of Allium species and support targeted breeding and conservation strategies. The observed number of alleles (Na) and effective number of alleles (Ne) support the existence of moderate polymorphism across the cp-SSR markers used. Markers such as cp-SSR09 and cp-SSR10, which exhibited relatively higher number of alleles and diversity indices, are particularly useful in revealing genetic differences among Allium accessions. In contrast, monomorphic markers such as cp-SSR02, cp-SSR03, and cp-SSR11 are likely associated with conserved regions of the chloroplast genome, offering limited intraspecies resolution but potential value for interspecific or phylogenetic studies (; ).

Generally, low Ho values align with the uniparental (mostly maternal) inheritance and haploid nature of the chloroplast genome, as well as the self-pollinating behaviour of many Allium species (). Nevertheless, the high He and I values of cp-SSR10 demonstrate its potential for distinguishing diverse genotypes and tracking lineage relationships. These findings are consistent with previous cp-SSR studies in Allium, where allele numbers typically ranged between 2 and 5 (), and in other genera, such as Ziziphus (). Compared with nuclear SSRs and EST-SSRs, cp-SSRs tend to be less polymorphic, because they are located in more conserved regions of the genome (; ; ). However, their high cross-species transferability and evolutionary stability make them ideal for phylogenetic studies and for characterizing maternal lineages.

The cross-transferability results indicate a broad genetic base within the genus Allium. High transferability rates in species such as A. viviparum and A. negianum suggest their close genetic relationships with other members of the genus and their potential utility in breeding programs. In contrast, lower transferability in accessions such as A. altaicum indicates possible genomic divergence or evolutionary distance. These patterns of allele sharing and divergence can be exploited for both germplasm conservation and introgression breeding strategies aimed at enhancing stress tolerance or other desirable traits. Overall, the cp-SSR markers were effective in evaluating genetic diversity, understanding evolutionary relationships, and identifying candidate accessions for conservation and breeding. The combination of highly polymorphic and conserved markers allows a dual-purpose application: detailed intraspecies diversity analysis and broader phylogenetic studies across underutilized and cultivated Allium species.

Owing to their uniparental inheritance and lack of recombination, chloroplast SSR markers are well-suited for studying genetic relationships, evolutionary history, and domestication processes in plants (Sharma et al., 2020a). In the present study, these markers effectively distinguished underutilized and cultivated Allium species into six well-defined genetic clusters. The genetic homogeneity of Cluster I highlights the close relatedness among A. tuberosum and its allies, likely due to shared ancestry and cultivation across similar ecological regions, which aligns with the findings of , who reported similar clustering patterns. Cluster II, composed of various underutilized accessions, including A. chinense and A. fragrans, reflects cytoplasmic similarity and ecological coherence, corroborating earlier findings by . The presence of both cultivated and underutilized species in Cluster III emphasizes its potential as a genetic bridge, with A. fistulosum and related species offering valuable traits such as disease resistance and abiotic stress tolerance. Cluster IV, comprising cultivated A. cepa genotypes, exhibited reduced diversity, a consequence of domestication and selective breeding, which is consistent with domestication bottlenecks observed by . Cluster V presented the greatest genetic divergence, harboring underutilized species adapted to niche environments such as high altitudes, in agreement with findings from and , highlighting their value as reservoirs of unique alleles. Cluster VI, comprising A. hookeri, presented species-specific genetic uniformity, likely due to restricted geographical distribution and domestication. Despite their narrow diversity, the members of this cluster possess unique traits that are valuable for region-specific applications. The high genetic variation observed in Clusters III and V underscores the evolutionary potential of underutilized species, reaffirming the importance of integrating underutilized relatives into breeding programs to increase stress tolerance, disease resistance, and adaptability in cultivated Allium species.

Understanding population structure is critical for effective germplasm conservation, trait mapping, and breeding applications. The identification of three distinct genetic clusters among the 96 underutilized Allium accessions aligns with earlier findings by , who also reported three subpopulations using cp-SSR markers in Allium accessions. The moderate allele–frequency divergence (0.11) and varying intracluster distances observed in this study suggest both shared ancestry and independent evolutionary trajectories among the populations. The relatively high proportion of membership in pop3 (50.2%) suggests a broad and genetically diverse group, potentially encompassing accessions with mixed ancestry. In contrast, pop2, with a lower proportion (18.5%), may represent a more genetically uniform or isolated subset. The findings also correlate with those of , who identified two groups in 96 underutilized Allium accessions through STRUCTURE analysis and five groups via discriminant analysis, illustrating how methodology and marker type influence the resolution of population structure. reported four clusters using nuclear SSRs, highlighting differences attributable to marker origin (chloroplast vs. nuclear). Similarly, reported population differentiation in Allium landraces based on accession type rather than phenotypic traits. These collective observations underscore the utility of cp-SSR markers in deciphering maternal lineage and cytoplasmic diversity, which are particularly important for breeding strategies involving cytoplasmic male sterility or other organelle-linked traits. Thus, marker-based population structure analysis not only facilitates an understanding of genetic diversity but also provides a valuable framework for selecting parental lines and managing Allium germplasm effectively.

Together, the PCA and heatmap analyses provided a nuanced understanding of the genetic diversity and structure within the underutilized Allium accessions. The separation along PC1 likely reflects deep evolutionary divergence, whereas PC2 captures more recent or subtle genetic differentiation. The broader dispersion observed in one cluster indicates high intragroup variability, potentially representing genetically diverse underutilized relatives with adaptive significance. Moreover, the tighter grouping of the second cluster suggests a subset of accessions with conserved genomic features, possibly shaped by shared ancestry or ecological adaptation. These observations are consistent with common patterns in plant population genetics, where variable levels of diversity are often observed within and among clusters (Spanoghe et al., 2020). Notably, the PCA-based clustering results corresponded well with the three subpopulations identified through STRUCTURE analysis, supporting the robustness and complementary nature of both methods. The overlapping zones in the PCA further corroborate previous studies reporting gene flow and admixture among Allium species (Xiong et al., 2022). The heatmap visualization further reinforced the PCA outcomes by graphically representing the levels of genetic similarity and divergence, with the dendrogram effectively grouping genetically close accessions. Similar integrative approaches have proven valuable in deciphering population structure and evolutionary relationships in Allium and other crop species (; Saina et al., 2023). Overall, these findings not only validate the genetic groupings but also emphasize the utility of multivariate and hierarchical clustering tools in germplasm characterization, aiding in the selection of genetically diverse and elite accessions for breeding and conservation programs.

Conclusion

In conclusion, this study has significantly advanced our understanding of the genetic diversity of underutilized Allium species through the development and application of novel chloroplast SSR markers. The identification of 22 cp-SSR motifs from the A. fistulosum chloroplast genome, with 20 markers exhibiting high polymorphism and stability, provides a robust toolkit for genetic analysis in Allium. The high level of polymorphism (89.2%) observed across 96 underutilized Allium species underscores the effectiveness of these markers in capturing genetic variation. The population structure analysis revealed three distinct genetic clusters, complemented by phylogenetic grouping into six major clusters, which offers valuable insights into the evolutionary relationships and genetic differentiation within the genus. These findings have important implications for Allium conservation strategies and breeding programs, highlighting the potential of underutilized germplasm as a reservoir of genetic diversity for crop improvement. This study demonstrates the utility of cp-SSR markers in revealing the complex genetic tapestry of Allium species, paving the way for precision-guided conservation efforts and the development of improved cultivars with enhanced traits such as disease resistance and stress tolerance.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Author contributions

YK: Software, Methodology, Writing – review & editing, Supervision, Writing – original draft, Investigation, Conceptualization, Funding acquisition, Formal Analysis, Visualization, Data curation, Resources, Project administration, Validation. PM: Formal Analysis, Writing – original draft, Methodology, Investigation. AC: Writing – review & editing, Formal Analysis. KR: Writing – review & editing, Resources, Formal Analysis. SS: Software, Data curation, Formal Analysis, Methodology, Writing – original draft. MJ: Writing – original draft, Methodology, Software, Data curation. AP: Methodology, Writing – original draft, Software, Data curation. VH: Writing – original draft. AR: Methodology, Formal Analysis, Data curation, Writing – review & editing. SM: Writing – review & editing, Formal Analysis, Methodology, Data curation. AK: Validation, Methodology, Writing – original draft, Software, Formal Analysis. HB: Investigation, Software, Data curation, Writing – original draft, Project administration, Methodology. AG: Writing – review & editing, Methodology, Data curation. RK: Resources, Writing – review & editing, Software, Formal Analysis. KP: Methodology, Software, Data curation, Writing – review & editing. VM: Writing – review & editing, Funding acquisition, Resources, Project administration.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Acknowledgments

The authors thank the Indian Council of Agricultural Research- National Bureau of Plant Genetic Resources, New Delhi and its regional station, Bhowali, Uttarakhand, India for sharing underutilized species of Allium used in the research.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1645145/full#supplementary-material

References

  • 1

    AnandhanS.NairA.KumkarD. S.GopalJ. (2015). Retrotransposon based TRAP marker displays diversity among onion (Allium cepa L.) genotypes. Scientia Hortic.190, 123127. doi: 10.1016/j.scienta.2015.04.024

  • 2

    AnwarG. M.HelmeyR. K.MoustafaY. M. M. (2017). Genetic diversity among some species of the genus Allium L. using SSR and ISSR markers. Egyptian J. Genet. Cytology.46, 221233. doi: 10.21608/ejgc.2018.9516

  • 3

    ArifinN. S.OzakiY.OkuboH. (2000). Genetic diversity in Indonesian shallot (Allium cepa var. ascalonicum) and Allium × wakegi revealed by RAPD markers and origin of A. × wakegi identified by RFLP analyses of amplified chloroplast genes. Euphytica111, 2331. doi: 10.1023/A:1003732403402

  • 4

    BaldwinS.Pither JoyceM.WrightK.ChenL.McCallumJ. (2012). Development of robust genomic simple sequence repeat markers for estimation of genetic diversity within and among bulb onion (Allium cepa L.) populations. Mol. Breeding.30, 14011411. doi: 10.1007/s11032-012-9727-6

  • 5

    BrahimiA.HajjajH.MazouzH. (2024). Genetic characterization of white onion genotypes (Allium cepa L.) using ISSR markers. Ecol. Genet. Genomics33, 100287. doi: 10.1016/j.egg.2024.100287

  • 6

    BrahimiA.LandschootS.BekaertB.HajjiL.HajjajH.AudenaertK.et al. (2022). Exploring the genetic and phenotypic diversity within and between onion (Allium cepa L.) ecotypes in Morocco. J. Genet. Eng. Biotechnol.20, 96. doi: 10.1186/s43141-022-00381-w

  • 7

    ChalbiA.Chikh-RouhouH.MezghaniN.SlimA.FayosO.Bel-KadhiM. S.et al. (2023). Genetic diversity analysis of onion (Allium cepa L.) from the arid region of Tunisia using phenotypic traits and SSR markers. Horticulturae9, 1098. doi: 10.3390/horticulturae9101098

  • 8

    ChandS. K.NandaS.JoshiR. K. (2018). Genetics and molecular mapping of a novel purple blotch-resistant gene ApR1 in onion (Allium cepa L.) using STS and SSR markers. Mol. Breed.38, 109. doi: 10.1007/s11032-018-0864-4

  • 9

    DamonS. J.HaveyM. J. (2014). Quantitative trait loci controlling amounts and types of epicuticular waxes in onion. J. Am. Soc. Hortic. Sci.139, 597602. doi: 10.21273/JASHS.139.5.597

  • 10

    dos SantosM. D.RagassiC. F.FonsecaM. E.BuzarA. G.OliveiraV. R.de MeloP. C.et al. (2012). Genetic diversity of tropical-adapted onion germplasm assessed by RAPD markers. Horticultura Bras.30, 112118.

  • 11

    DuangjitJ.BohanecB.ChanA. P.TownC. D.HaveyM. J. (2013). Transcriptome sequencing to produce SNP-based genetic maps of onion. TAG. Theor. Appl. Genet.126, 20932101. doi: 10.1007/s00122-013-2121-x

  • 12

    EbertD.PeakallR. O. D. (2009). Chloroplast simple sequence repeats (cp-SSRs): technical resources and recommendations for expanding cp-SSR discovery and applications to a wide array of plant species. Mol. Ecol. Resources.9, 673690. doi: 10.1111/j.1755-0998.2008.02319.x

  • 13

    EvannoG.RegnautR.GoudetJ. (2005). Detecting the number of clusters of individuals using the software structure: a simulation study. Mol. Ecol.14, 26112620. doi: 10.1111/j.1365294x.2005.02553.x

  • 14

    FengS.JiaoK.ZhangZ.YangS.GaoY.JinY.et al. (2023). Development of chloroplast microsatellite markers and evaluation of genetic diversity and population structure of cutleaf groundcherry (Physalis angulata L.) in China. Plants12, 1755. doi: 10.3390/plants12091755

  • 15

    FritschR. M.FriesenN. (2002). “Evolution, domestication and taxonomy,” in Allium crop science: recent advances (CABI, Wallingford UK), 530. doi: 10.1079/9780851995106.0005

  • 16

    GowdT. Y. M.DeoC.ManjunathagowdaD. C.MahajanV.DuttaR.BhutiaN. D.et al. (2024). Deciphering genetic diversity phylogeny and assembly of Allium species through micro satellite markers on nuclear DNA. Heliyon10, e31650. doi: 10.1016/j.heliyon.2024.e31650

  • 17

    GowdT. Y. M.ManjunathagowdaD. C.MahajanV.DuttaR.BhutiaN. D.SinghB.et al. (2023). Deployment of intron length polymorphic (ILP) markers in dissipating diversity of Allium species. South Afr. J. Bot.160, 157165.

  • 18

    GuoQ.GuoL.LiY.YangH.HuX.SongC.et al. (2022). Development and characterization of microsatellite markers based on the chloroplast genome of tree peony. Genes13, 1543. doi: 10.3390/genes13091543

  • 19

    GuptaA. J.AnandhanS.KadS. K.SinghM. (2020). Confirmation of hybridity in DOGR hybrids of onion (Allium cepa L.) using SSR markers. Vegetable Sci.47, 183188. doi: 10.61180/vegsci.2020.v47.i2.03

  • 20

    HabsatouB.BroumP.BellangerL.CrouzillatD. (2024). Genetic diversity of red onion using microsatellite (SSRs) markers. Scholars Acad. J. Biosci.8, 278284. doi: 10.36347/sajb.2024.v12i08.007

  • 21

    HanciF.GökçeA. F. (2016). Molecular characterization of Turkish onion germplasm using SSR markers. Czech J. Genet. Plant Breeding.52, 7176. doi: 10.17221/162/2015-CJGPB

  • 22

    HancıF.PaşazadeE. (2025). A comparison of efficiency parameters of SRAP and ISSR markers in revealing variation in Allium germplasm. Horticulturae11, 294. doi: 10.3390/horticulturae11030294

  • 23

    HuX.HuJ.ZhangY.JiangS.YuQ. (2022). Genetic structure and demographic history of Allium mongolicum based on SSR markers. Plant Systematics Evol.308, 12. doi: 10.1007/s00606-021-01802-y

  • 24

    HuangL. S.SunY. Q.JinY.GaoQ.HuX. G.GaoF. L.et al. (2018). Development of high transferability cp-SSR markers for individual identification and genetic investigation in Cupressaceae species. Ecol. Evol.8, 49674977. doi: 10.1002/ece3.4053

  • 25

    HuangJ.YangX.ZhangC.YinX.LiuS.LiX. (2015). Development of chloroplast microsatellite markers and analysis of chloroplast diversity in Chinese jujube (Ziziphus jujuba Mill.) and wild jujube (Ziziphus acidojujuba Mill.). PloS One10, e0134519. doi: 10.1371/journal.pone.0134519

  • 26

    HuoY.GaoL.LiuB.YangY.KongS.SunY.et al. (2019). Complete chloroplast genome sequences of four Allium species: comparative and phylogenetic analyses. Sci. Rep.9, 12250. doi: 10.1038/s41598-019-48708-x

  • 27

    JayaswallK.BhandawatA.SharmaH.YadavV. K.MahajanV.SinghM. (2019). Characterization of Allium germplasms for conservation and sustainable management using SSR markers. Indian J. Tradit Knowl.18, 193199.

  • 28

    JayaswallK.SagarR.JayaswalD.KumarA.SinghS. P.SethR.et al. (2024). Development of Allium cepa potential intron polymorphism markers for molecular breeding of Alliums. South Afr. J. Bot.164, 209220. doi: 10.1016/j.sajb.2023.11.050

  • 29

    JayaswallK.SharmaH.BhandawatA.SagarR.JayaswalD.KumarA.et al. (2022). Chloroplast derived SSRs reveals genetic relationships in domesticated alliums and wild relatives. Genet. Resour Crop Evol.69, 363372. doi: 10.1007/s10722-021-01235-z

  • 30

    JayaswallK.SharmaH.JayaswalD.SagarR.BhandawatA.KumarA.et al. (2023). Development of chloroplast derived SSR markers for genus Allium and their characterization in the allies for genetic improvement of Alliums. South Afr. J. Botany.162, 304313. doi: 10.1016/j.sajb.2023.09.021

  • 31

    KaricL.GolzardiM.GlamoclijaS. (2018). Genetic diversity assessment of Allium cepa L. cultivars from Bosnia and Herzegovina using SSR makers. Genet. Mol. Res.17, gmr16039870. doi: 10.4238/gmr16039870

  • 32

    KhadeY. P.SalunkheS. R.ManjunathagowdaD. C.SinhasaneS. R.GowdT. Y. M.MahajanV.et al. (2022). Molecular characterization of short-day onion genotypes by intron length polymorphic (ILP) markers. Genet. Resour Crop Evol.69, 20772086. doi: 10.1007/s10722-022-01398-3

  • 33

    KhadeY. P.SinhasaneS. R.MainkarP.RaiK. M.SalunkheS.SinghP. R.et al. (2024). Exploring sequence-related amplified polymorphism (SRAP) markers for assessing genetic diversity in onion (Allium cepa L.) genotypes and their wild relatives. Genet. Resour. Crop Evol.72, 399415. doi: 10.1007/s10722-024-01978-5

  • 34

    KharA.LawandeK. E.NegiK. S. (2011). Microsatellite marker based analysis of genetic diversity in short day tropical Indian onion and cross amplification in related Allium spp. Genet. Resour. Crop evolution.58, 741752. doi: 10.1007/s10722-010-9616-y

  • 35

    KhosaJ. S.DhattA. S.NegiK. S.KharA. (2014). Utility of simple sequence repeat (SSR) markers to realize worth of germplasm in genus Allium. Indian J. Plant Genet. Resour.27, 238245. doi: 10.5958/0976-1926.2014.00020.5

  • 36

    KhosaJ. S.McCallumJ.DhattA. S.MacknightR. C. (2015). Enhancing onion breeding using molecular tools. Plant Breeding.135, 920. doi: 10.1111/pbr.12330

  • 37

    KishaT. J.CramerC. S. (2011). Determining redundancy of short-day onion accessions in a germplasm collection using microsatellite and targeted region amplified polymorphic markers. J. Am. Soc. Hortic. Sci.136, 129134. doi: 10.21273/JASHS.136.2.129

  • 38

    KumarS.StecherG.LiM.KnyazC.TamuraK. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. evolution.35, 15471549. doi: 10.1093/molbev/msy096

  • 39

    KuttyM.Veere GowdaR.AnandL. (2006). Analysis of genetic diversity among Indian short-day onion (Allium cepa L.) cultivars using RAPD markers. J. Hortic. Sci. Biotechnol.81, 774778. doi: 10.1080/14620316.2006.11512136

  • 40

    LabateJ. A.GlaubitzJ. C.HaveyM. J. (2020). Genotyping by sequencing for SNP marker development in onion. Genome63, 607613. doi: 10.1139/gen-2020-0011

  • 41

    LiC.ZhengY.HuangP. (2020). Molecular markers from the chloroplast genome of rose provide a complementary tool for variety discrimination and profiling. Sci. Rep.10, 12188. doi: 10.1038/s41598-020-68092-1

  • 42

    LyngkhoiF.SainiN.GaikwadA. B.ThirunavukkarasuN.VermaP.SilvarC.et al. (2021). Genetic diversity and population structure in onion (Allium cepa L.) accessions based on morphological and molecular approaches. Physiol. Mol. Biol. Plants.27, 25172532. doi: 10.1007/s12298-021-01101-3

  • 43

    MahajanV.JakseJ.HaveyM. J.LawandeK. E. (2009). Genetic fingerprinting of onion cultivars using SSR markers. Indian J. Horticulture66, 6268.

  • 44

    MallorC.Arnedo-AndresM. S.Garces-ClaverA. (2014). Assessing the genetic diversity of Spanish Allium cepa landraces for onion breeding using microsatellite markers. Scientia Horticulturae.170, 2431. doi: 10.1016/j.scienta.2014.02.040

  • 45

    MandakhU.BattserenM.GanbatD.AyangaT.AdiyaZ.BorjigidaiA.et al. (2020). Folk nomenclature of plants in Cistanche deserticola-associated community in South Gobi, Mongolia. Plant Diversity.42, 434442. doi: 10.1016/j.pld.2020.09.008

  • 46

    MansourE. M.ElballaM. M.El HusseinA. A.AbdallaA. W. H.GadirI. K. A.AbboA. S.et al. (2020). Assessment of genetic variability among onion (Allium cepa L.) cultivars using RAPD and SRAP molecular markers. Univ. Khartoum J. Agric. Sci.28, 124. doi: 10.53332/uofkjas.v28i.177

  • 47

    McCallumJ.BaldwinS.ShigyoM.DengY.van HeusdenS.PitherJoyceM.et al. (2012). Allium Map-A comparative genomics resource for cultivated Allium vegetables. BMC Genomics13, 168. doi: 10.1186/1471-2164-13-168

  • 48

    McCallumJ.LeiteD.PitherJoyceM.HaveyM. J. (2001). Expressed sequence markers for genetic analysis of bulb onion (Allium cepa L.). Theor. Appl. Genet.103, 979991. doi: 10.1007/s001220100630

  • 49

    McCauleyD. E. (1995). The use of chloroplast DNA polymorphism in studies of gene flow in plants. Trends Ecol. evolution.10, 198202. doi: 10.1016/S0169-5347(00)89052-7

  • 50

    MetsaluT.ViloJ. (2015). ClustVis: a web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap. Nucleic Acids Res.43, W566W570. doi: 10.1093/nar/gkv468

  • 51

    MohapatraP. P.SeleimanM. F.MandalR.PramanikK.MaityT. K.TarafdarJ.et al. (2023). Efficiency of RAPD and SSR markers in assessing genetic diversity in summer onion (Allium cepa L.) genotypes. Notulae Botanicae Horti Agrobotanici ClujNapoca51, 13369. doi: 10.15835/nbha51313369

  • 52

    MonteverdeE.GalvanG. A.SperanzaP. (2015). Genetic diversification of local onion populations under different production systems in Uruguay. Plant Genet. Resour.13, 238246. doi: 10.1017/S1479262114000963

  • 53

    MurrayM. G.ThompsonW. F. (1980). Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res.8, 43214326. doi: 10.1093/nar/8.19.4321

  • 54

    NandaS.ChandS. K.MandalP.TripathyP.JoshiR. K. (2016). Identification of novel source of resistance and differential response of Allium genotypes to purple blotch pathogen, Alternaria porri (Ellis) Ciferri. Plant Pathol. J.32, 519. doi: 10.5423/PPJ.OA.02.2016.0034

  • 55

    NouriN.DevineniN.WereV.KhanbilvardiR. (2021). Explaining the trends and variability in the United States tornado records using climate teleconnections and shifts in observational practices. Sci. Rep.11, 1741. doi: 10.1038/s41598-021-81143-5

  • 56

    OmelchenkoD. O.KrinitsinaA. A.BelenikinM. S.KonorovE. A.KuptsovS. V.LogachevaM. D.et al. (2020). Complete plastome sequencing of Allium paradoxum reveals unusual rearrangements and the loss of the ndh genes as compared to Allium ursinum and other onions. Gene726, 144154. doi: 10.5423/PPJ.OA.02.2016.0034

  • 57

    PeakallR.SmouseP. E. (2012). GenAlEx 6.5: Genetic Analysis in Excel. Population genetic software for teaching and research—An update. Bioinf. Appl. Note.28, 25372539. doi: 10.1093/bioinformatics/bts460

  • 58

    PowellW.MorganteM.AndreC.McNicolJ. W.MachrayG. C.DoyleJ. J.et al. (1995). Hypervariable microsatellites provide a general source of polymorphic DNA markers for the chloroplast genome. Curr. Biol.5, 10231029. doi: 10.1016/S0960-9822(95)00206-5

  • 59

    RicciardiL.MazzeoR.MarcotrigianoA. R.RainaldiG.IovienoP.ZonnoV.et al. (2020). Assessment of genetic diversity of the “Acquaviva red onion” (Allium cepa L.) apulian landrace. Plants9, 260. doi: 10.3390/plants9020260

  • 60

    RiveraA.MallorC.Garcés-ClaverA.García-UlloaA.PomarF.SilvarC. (2016). Assessing the genetic diversity in onion (Allium cepa L.) landraces from northwest Spain and comparison with the European variability. New Z. J. Crop Hortic. Science.44, 103120. doi: 10.1080/01140671.2016.1150308

  • 61

    RohlfF. J. (1998). “NTSYSpc,” in Numerical taxonomy and multivariate analysis system. Version 2.0 (User Guide: Exeter Software, Setauket, New York), 31.

  • 62

    SahooJ.MahantyB.MishraR.JoshiR. K. (2023). Development of SNP markers linked to purple blotch resistance for marker-assisted selection in onion (Allium cepa L.) breeding. 3 Biotech.13, 137. doi: 10.1007/s13205-023-03562-7

  • 63

    SainaJ. K.LiZ. Z.NgaregaB. K.GituruR. W.ChenJ. M.LiaoY. Y. (2023). Exploring the genetic diversity and population structure of Ailanthus altissima using chloroplast and nuclear microsatellite DNA markers across its native range. Front. Plant Science.14. doi: 10.3389/fpls.2023.1197137

  • 64

    ScholtenO. E.van AauwenM. P. W.ShahinA.HendrickxP. M.KeizerL. C. P.BurgerK.et al. (2016). SNP-markers in Allium species to facilitate introgression breeding in onion. BMC Plant Biol.16, 187. doi: 10.1186/s12870-016-0879-0

  • 65

    SharmaH.BhandawatA.RawatS. (2020a). Cross-transferability of SSR markers developed in Rhododendron species of Himalaya. Mol. Biol. Rep.47, 63996406. doi: 10.1007/s11033-020-05606-0

  • 66

    SharmaH.HyvonenJ.PoczaiP. (2020). Development of chloroplast microsatellite markers for giant ragweed (Ambrosia trifida). Appl. Plant Sci.8, e11313. doi: 10.1002/aps3.11313

  • 67

    ShuklaN.KuntalH.ShankerA.SharmaS. N. (2018). Mining and analysis of simple sequence repeats in the chloroplast genomes of genus Vigna. Biotechnol. Res. Innovation.2, 918. doi: 10.1016/j.biori.2018.08.001

  • 68

    SimóJ.PascualL.CanizaresJ.CasanasF. (2014). Spanish onion landraces (Allium cepa L.) as sources of germplasm for breeding calcots: a morphological and molecular survey. Euphytica195, 287300. doi: 10.1007/s10681-013-0995-y

  • 69

    SonJ. H.ParkK. C.LeeS. I.KimJ. H.KimN. S. (2012). Species relationships among Allium species by ISSR analysis. Horticulture Environment Biotechnol.53, 256262. doi: 10.1007/s13580-012-0130-3

  • 70

    SpanogheM. C.MariqueT.RivièreJ.MoulinM.DekuijperC.NirshaA.et al. (2020). Genetic patterns recognition in crop species using self-organizing map: The example of the highly heterozygous autotetraploid potato (Solanum tuberosum L.). Genet. Resour. Crop Evolution.67, 947966. doi: 10.1007/s10722-020-00894-8

  • 71

    SudhaG. S.RameshP.SekharA. C.KrishnaT. S.BramhachariP. V.RiazunnisaK. (2019). Genetic diversity analysis of selected onion (Allium cepa L.) germplasm using specific RAPD and ISSR polymorphism markers. Biocatalysis Agric. Biotechnol.17, 110118. doi: 10.1016/j.bcab.2018.11.007

  • 72

    Van HeusdenA. W.Van OoijenJ. W.Vrielinkvan GinkelR.VerbeekW. H. J.WietsmaW. A.KikC. (2000). A genetic map of an interspecific cross in Allium based on amplified fragment length polymorphism (AFLPTM) markers. Theor. Appl. Genet.100, 118126. doi: 10.1007/s001220050017

  • 73

    VillanoC.EspositoS.CarucciF.IorizzoM.FruscianteL.CarputoD.et al. (2019). High-throughput genotyping in onion reveals structure of genetic diversity and informative SNPs useful for molecular breeding. Mol. Breed.39, 111. doi: 10.1007/s11032-018-0912-0

  • 74

    WangY.ZhangY. (2022). Advances in molecular breeding of vegetable crops. Horticulturae8, 821. doi: 10.3390/horticulturae8090821

  • 75

    WrightS. (1984). Evolution and the genetics of populations, volume 4: variability within and among natural populations (Chicago: University of Chicago Press).

  • 76

    XiongY.XiongY.ShuX.YuQ.LeiX.LiD.et al. (2022). Molecular phylogeography and intraspecific divergences in siberian wildrye (Elymus sibiricus L.) wild populations in China, inferred from chloroplast DNA sequence and cp-SSR markers. Front. Plant Science.13. doi: 10.3389/fpls.2022.862759

  • 77

    YamashitaK.TsukazakiH.KojimaA.OharaT.WakoT. (2010). Inheritance mode of male sterility in bunching onion (Allium fistulosum L.) accessions. Euphytica173, 357367. doi: 10.1007/s10681-009-0101-7

  • 78

    YouF. M.HuoN.GuY. Q.LuoM. C.MaY.HaneD.et al. (2008). BatchPrimer3: a high throughput web application for PCR and sequencing primer design. BMC Bioinf.9, 253. doi: 10.1186/1471-2105-9-253

Summary

Keywords

cp-SSR, Allium fistulosum, underutilized species, population structure, cross transferability

Citation

Khade YP, Mainkar P, Chandanshive A, Rai KM, Sinhasane SR, Jadhav M, Patil A, Hembade VL, Radhakrishna A, More SJ, Khar A, Bhandari HR, Gupta AJ, Kale RB, Prakash K and Mahajan V (2025) Harnessing chloroplast SSRs to decipher genetic diversity in underutilized Allium species. Front. Plant Sci. 16:1645145. doi: 10.3389/fpls.2025.1645145

Received

12 June 2025

Accepted

18 August 2025

Published

15 September 2025

Volume

16 - 2025

Edited by

Mohan Lal, North East Institute of Science and Technology (CSIR), India

Reviewed by

Shivani Rohilla, Forest Research Institute, India

Ebubekir Paşazade, Orta Karadeniz Geçit Kuşağı Tarımsal Araştırma Enstitüsü Müdürlüğü, Türkiye

Updates

Copyright

*Correspondence: Yogesh P Khade, ; Hem Raj Bhandari, ; Vijay Mahajan,

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics