Skip to main content

Genetic variation and structure of maize populations from Saoura and Gourara oasis in Algerian Sahara



The ability of maize populations/landraces to tolerate drastically extreme environments over the past four centuries in Algeria leads to characterize these genetic resources for germplasm management as well as the identification of the best landraces useful for genetic improvement. Thus, the aim of the present work was a fingerprinting of an Algerian maize collection (47 landraces) from Saharan oasis by using 24 agro-morphological traits and18 Simple Sequence Repeats to evaluate genetic diversity and population structure.


Phenotypic traits showed large significant variation in which earliness, plant size, ear and kernel features and crop yield appeared the most discriminant traits among landraces by using principal component analysis (PCA). One hundred ninety-seven different alleles were detected with a high mean number of allele per locus (10.9). The selected SSR were highly informative with PIC values > 0.65 as well as an overall genetic diversity (0.47) highlighting a broad genetic variability in the analyzed landraces. Genetic structure analysis revealed a high genetic differentiation among the 47 maize landraces with an overall Fst value (0.33). Cluster analysis for morphological traits as well as for SSR markers grouped the 47 Algerian populations regardless their geographic origin.


Maize from Algerian desert harbors a wide genetic diversity offering a source of novel/unique alleles useful for maize breeding programs to face the ongoing and future major challenge, the climate changes.


Maize is one of the most diverse crop in the world [1], characterized by a high degree of genetic variability due to an extended selection process before its spread from Central America to the other regions of the world [2, 3].

The introduction of maize in African agro-systems goes back to the sixteenth century, as previously reported [4]. In Algeria, maize cultivation was reported in several regions (Kabylie, Tell and Saharan oasis) by general Duval in 1856 [5], but in the Algerian desert since the sixteenth century [6]. Maize is still appreciated by the Arabo-Berbers of Saharan oasis in South Algeria where several names were attributed to this cereal: Mastoura, Ktania, Safra, Abidiya, Dora (in Arabic), Tifsi or Engafouli (in Tamacheq) [6]. Over this time, maize adapted to the drastic soil and climatic conditions of Algerian desert (dry subtropical climate), where the farmers maintained the useful high genetic diversity through traditional low-pressure selection [7]. More interestingly, crop evolution and maize divergent selection took place and probably still takes place in the traditional farming systems [8]. Therefore, Saharan oasis traditional plant populations or landraces could represent a unique and valuable germplasm.

Although phenotypic variations represent an ideal tool to determine, the agronomic crops performance as well as being inexpensive and comfortable [9], they suffer from low degree of polymorphism and they are affected by environment, thereby determining a limited heritability in germplasm collections [10]. By contrast, DNA-based molecular markers, as microsatellite (Simple Sequence Repeat, SSR), are so far the most used for large scale maize genetic diversity analysis and population structure being highly polymorphic, codominant, abundant in the genome and highly reproducible [11,12,13,14,15,16,17,18]. However, morphological and molecular characterizations are both essential and complementary to study and manage genetic diversity in maize germplasm as well as other crops.

Although maize germplasm from Algerian desert was recently studied by phenotypic and genetic analysis [7, 19], to date the genetic maize pool available in Algeria is not completely evaluated. In particular, Djemel et al. [7] reported a preliminary approach using 10 open-pollinated populations based only on agronomic evaluation, whereas Aci et al. [19] showed genetic analysis of a limited number (15) of Algerian maize populations and therefore not representative of Algerian germplasm. The aim of the present paper was to fingerprint a large collection of Algerian maize landraces (47) from Saharan oasis by both phenotypic traits and molecular markers (18 SSRs) to evaluate genetic diversity and populations structure useful to select genotypes harboring traits of agronomic interest for future breeding programs.


Morphological diversity

ANOVA showed a significant block effect for 4 out of 24 traits (DA, KT, 1000 KW and HMC). The seven checks in each replicate significantly differed for many traits except EMG, EV, ASI, ED, EL, HMC and KYP, whereas landraces significantly differed for 13 out of 24 traits including DS, DA, EH, PLH, ERN, ED, CD, RD, KL,KW, EW, 1000 KW and HMC (Additional file 1: Table S4).

The basic statistics of quantitative traits were calculated among all landraces (Table 1). According to the coefficient of variation (CV%), Algerian maize landraces displayed a wide phenotypic variation for ASI (125%), EH (43.88%) PLH (28.50%), KYP (36.15%) and EW (35.43%). By contrast, the lowest variations were observed for kernel features (KL, KW, KT and K%), phenological traits (DA and DS) and ED. The remaining traits showed intermediate levels of variation. At early developmental stages rather all the landraces showed a significant EV, where ENR appeared the most vigorous (Table 1). The earliest landrace was EHA with 69.5 and 68 days after sowing to silking and anthesis, respectively, whereas the latest was KMA (98.52 and 97 days after sowing to silking and anthesis, respectively). A short mean ASI (anthesis -silking interval) was scored (1.22 days), which ranged from 2.52 to 4.47 in BEC and BML, respectively (Table 1), and negative values were detected for seven landraces (DHT, BEC, KTA, DDL, ONA, TIF and AGU). The mean values of PLH and EH was 127.63 and 43.22 cm, respectively; KMA showed the highest (260 cm and 85.70 cm, respectively), while the lowest PLH and EH (79.41 cm and 17.40 cm) were observed in landraces AGU and MNS, respectively (Table 1). NEP ranged from 1.79 (YAK) to 3.19 (AGU). The mean KL was 11.05 mm, the longest value was observed in BSA2 (11.13 mm), while TKR recorded the shortest ones (7.84 mm). BHT and AOR2 showed the highest and lowest values in KW (10.18 and 7.71 mm, respectively), whereas BYA and AOR2 recorded the slimmest and the thickest kernel (KT range from 5.28 to 4.10 mm, respectively). The average KYP was 1.02 Mg ha−1ranging from 1.18 (KEK) to7.04 (KMA). However, BTH showed the highest 1000 KW (323.6 g) whereas TBN recorded the lowest value (139.95 g) (Table 1).

Table 1 Means and variation of agro-morphological traits in the 47 Algerian maize landraces

The correlations among agro-morphological traits were included in Additional file 2: Table S5. The highest significant positive correlation was between DS and DA (r = 0.97) followed by PLH and EH (r = 0.92). Interestingly, KYP is highly correlated with PLH and EH, DA and DS, ERN and NKR, while others significant correlations between traits were observed (Additional file 2: Table S5). A limited number of highly significant negative correlations was found, such as KT with DS and DA (r = − 0.51 and − 0.46, respectively), ERN and KW (r = − 0.39), K% and 1000 KW (r = − 0.62) (Additional file 2: Table S5).

PCA based on agro-morphological traits displayed three principal components that contributed for 65% to overall phenotypic variability among the Saharan maize landraces (Additional file 3: Table S6). In particular, the PC1 explained 44.29% of total variation, where ear features, earliness, PLH and EH as well as KYP appeared those mainly involved. The PC2accounted for 11.35% of total variation and the traits involved were 1000 KW, kernel sizes and NEP (Additional file 3: Table S6). Furthermore, EV, ASI and KP were traits contributing to the PC3 (9.35% of the total phenotypic variation).The first two PCs (55.64%) were able to distinguish the maize landraces according to their major contributing traits (Additional file 4: Figure S1). Genetic dissimilarity was calculated from agro-morphologic data by cluster analysis based on Euclidean distance and Ward’s method and displayed as heatmap too. Cluster analysis grouped the 47 landraces into three main clusters with distinguished genetic profiles (Fig. 1). The first enclosed the 15 early maturing, short stature and lowest yield landraces. The second cluster grouped 20 landraces with great kernels and high 1000 KW but with the lowest NEP. The third and last cluster was mainly composed of highly vigorous landraces with long ASI and the highest KYP.

Fig. 1

Agro-morphological traits heatmap of the 47 Algerian maize landraces based on Euclidean distances

Genetic diversity

SSR summary statistics and genetic indexes among 47Saharan maize landraces are reported in Table 2. All 18 SSR loci, evenly distributed throughout maize genome, were polymorphic. A high number of alleles (197) were detected in bulked DNA samples; allele numbers observed for each locus ranged from 6 (phi127) to 20 (umc1222) with an average of 10.94 alleles per locus (Table 2). The mean effective number of alleles (Ae), the total gene diversity (H) and Shannon’s information index (I), commonly used to assess population genetic diversity, were also estimated (Ae =1.92, H = 0.46 and I = 0.84) exhibiting the highest values for umc1335 (Ae = 3.00, H = 0.67 and I = 1.45) and the lowest for phi127 (Ae = 1.33, H = 0.25 and I = 0.47) (Table 2). Furthermore, the polymorphism information content (PIC) ranged from 0.86 (umc1335) to 0.31 (phi127) with an whole average value of 0.62 (Table 2). Total genetic differentiation Fst was moderate (0.33), while an excess of homozygotes was observed in our Algerian maize gene pool. A set of 18unique alleles were detected at 14 different SSR loci in the analyzed landraces, and the highest number was recorded in umc1222 (3 out of 20 alleles) as well as in umc1403 (2 out of 7). Thirteen landraces out of the 47 revealed unique alleles (Additional file 5: Table S7); the highest number was detected in DHT (3).

Table 2 Genetic diversity summary of statistics for 18 SSR loci across the 47 Algerian maize landraces

Average of different alleles (Na) was the highest in YAK (4.53) and the lowest in EID (2.63), and gene diversity (H) of a single landrace ranged from 0.38 (EID) to 0.60 (YAK). (Additional file 5: Table S7).

Finally, the analysis of the molecular variance (AMOVA) revealed that genetic diversity within populations contributes to 70% of overall genetic diversity (Additional file 6: Table S8).

Cluster analysis and genetic structure

Cluster analysis based on SSR data was performed using the genetic distances of Nei and Li [20] and UPGMA algorithm (Fig. 2). Forty-seven landraces were classified in three large clusters (from A to C) and two landraces as out-groups. In particular, cluster A enclosed the highest number of landraces (29), clusters B and C encompassed 14 and 2 landraces, respectively, while the landraces GHT and OAU were the out groups. The highest genetic distance was recorded between cluster C and A, while the lowest was observed between clusters A and B and the landraces OAU and GHT. The average gene diversity was the highest within cluster A and it was higher than total genetic diversity (Additional file 7: Table S9).

Fig. 2

Dendrogram of the 47 Algerian maize landraces based on SSR data according to UPGMA method. Different symbols indicate clusters (♦ = A, □ = B, ▲ = C, ■ = D,  = E)

The correlation coefficient (r = − 0.15, P < 0.001) for the two clustering matrices tasted by Mantel’s test revealed a highly significant but negative correlation between the agro-morphological and SSR information.

The differences among landraces were further highlighted by STRUCTURE analysis. The optimum number of genetic groups (K) within our collection was determined K = 2 based on ΔK peaks (Fig. 3). Considering the admixture coefficient (Q) ≥ 0.9 as the assignment probability of each population to a group, it was possible to assign 32 and 13 populations to group 1 (red) and 2 (green), respectively (Fig. 3). Two populations, BML and TKK, exhibited Q values of 0.686 and 0.88, respectively, thus showing an admixture genetic structure. Finally, allele frequency divergence between populations was 0.0479, while mean genetic distance within population was 0.6575 and 0.6012 for group 1 and 2, respectively.

Fig. 3

Genetic structure of the 47 maize landraces as inferred by STRUCTURE based on 18 SSR data set


Landraces represent an odd genetic resource for crops with a high genetic variation useful for modern breeding programs [17]. Approximately only 5% of maize genetic diversity is reported to be in commercial use [21], while the remaining germplasm have been neglected due to undesirable traits, such as low yield. In maize, landraces are low yielding varieties, but their large genetic diversity ensure crop adaptation to different stresses appearing adequate to meet the ongoing and future major challenge, the climate changes [22].

Therefore, the conservation of landraces as genetic resources and the management of such gene pools are important by ensuring suitable agronomic and genetic knowledge [23]. In the present study, we approached a genetic diversity assessment of 47 maize landraces from Saoura and Gourara Oasis (Algeria) using both agro-morphological and molecular data.

Morphological diversity

Algerian maize is a promising source of favorable alleles for stress tolerance, but it has been neglected until recently [7], thus the main purpose of the current research was to overview the agro-morphological potential of Saharan maize landraces based on quantitative traits. The estimation of variance revealed significant variability among landraces, suggesting a high phenotypic diversity within the Algerian maize collection. The maize populations varied greatly for plant size, ASI, kernel yield, and 10 ears weight. KMA was the latest, the tallest and the most productive landrace exhibiting the highest grain yield. Indeed, the late maturing maize landraces characterized by short ASI were also the most productive [24] At the phenological stages, all the landraces showed a normal flowering behavior (ASI = 1.22 days) except DHT, BEC, KTA, DDL, ONA, TIF and AGU which showed negative ASI values as a consequence of silk emergence prior to pollen emission. This feature could be a beneficial trait under drought conditions [25] causing delay in silk emergence, increasing the synchronicity between male and female flowering as well as resulting into higher grain yield. In this respect, negative and short ASI as secondary trait employed for breeding programs for drought tolerance was previously reported [26] Morphological traits were subjected to the environment influence as well as natural and human selection [27]. The high variability among flowering behavior (days to silking and anthesis as well as ASI) and crop yield, which are related to maize adaptation to the environment [28], suggested that Saharan landraces were adapted to a large range of environments supported also by alternative agronomic practices in different growing maize regions. Thus, our results are in agreement with a previous report on Algerian maize landraces, which pointed out that even though the landraces originated from subtropical regions, they showed adaptation to both sub-humid and dry continental environments [7]. Significant correlations among agro-morphological traits were identified, representing a useful tool for directing breeding programs [29]. However, environment plays also a pivotal role, affecting concurrently the traits in the same or opposite direction [30]. Such forceful correlation between agro-morphological traits in maize was already reported [15, 31].

The principal component analysis (PCA) was performed to cluster landraces based on discriminant morphologic traits. According to Clifford and Stephenson [32], the first three principal components included the overall phenotypic variation, taken into account discriminant traits able to distinguish the landraces. Among them, plant size, days to flowering, yield, ear as well as kernel features, and 1000 kernel weight were highly significant for clustering maize landraces. Our results are also in agreement with Gouesnard et al. [33], Beyene et al. [34] and Hartings et al. [35], which reported flowering, plant and ear size as the most important traits for grouping French, Ethiopian and Italian landraces, respectively. Further, Sharma et al. [36] also reported that kernel yield and 1000 kernel weight as highly informative traits in maize differentiations. CA as well as PCA analysis were not able to distinguish maize populations from the same geographic origin being dispersed in different clusters. Similar results were already reported by Hartings et al. [35] and Sharma et al. [36], which observed a lack of relationship between clustering and landraces geographic origin. Finally, it is noteworthy that landraces from dry, of poor soil and highly radiated environments origin showed adaptation traits as shortness and early flowering, considered important sources of gene diversity for developing drought tolerant genotypes [31, 34].

Genetic diversity

The first step for crop improvement strategies consist on the assessment of genetic variation and the relationships among landraces/varieties/populations. SSR markers represent an efficient tool in providing direct and reliable information for a helpful management and conservation of germplasm collections. The genetic diversity of 47 Algerian maize landraces was analyzed in bulk fingerprinting method, leading to detect high levels of polymorphism (10.94 alleles per locus). Higher results were reported by Wasala and Prasanna [37], which found 13.1 alleles per locus (across 42 SSR loci) studying 48 Indian landraces. A less number of alleles per locus compared to our study were detected in maize landraces from Switzerland, Ghana and Turkey with 8, 7.3 and 6.21 alleles per locus, respectively [15, 38, 39]. In a previous report, 5.8 alleles per locus were detected with a 6% polyacrylamide gel electrophoresis using the same SSR set of the present work [19]. The present study strengthened a wide genetic variation in our Algerian maize landraces from Sahara. As expected, we proved the higher efficiency of fluorescent dye labeled SSR markers, with a higher allele resolution and detection using DNA-sequencer, compared to the method adopted by Aci et al. [19]. The PIC average value (0.62) revealed an allelic variation in SSR loci higher than in the US germplasm (0.59), Swiss (0.52) and Ghana (0.50) maize landraces [11, 38, 39]. PIC is considered as the better parameter for measuring genetic diversity than the number of alleles, taking into account the relative frequencies of each allele [40]. This supports the informativeness of the chosen SSR to assess the Algerian maize genetic diversity.

Almost 50% of the SSR used in our study was highly informative (PIC > 0.65) and the overall gene diversity (0.46) was higher compared to the previous report on 15 Algerian landraces (0.40) [19] and comparable to those found in Ghanaian landraces [39] and Spanish populations [41].

Regardless the abundance of SSR in maize genome and their high polymorphism [12, 42], the dinucleotide motif SSR were reported to overestimate the genetic diversity and the number of alleles per locus [43], because of their higher mutation rate compared to oligonucleotide motif SSR [44]. In agreement, 7 out of 9 most informative SSR are dinucleotide repeats (umc1222 and umc1335), showing a high number of alleles (20 and 19, respectively). Furthermore, plant populations harboring unique alleles, due to SSR high mutation rate [45], but also to natural selection of specific alleles along the adaptation to environments [46], are reported to be potential for crop improvement as source of new alleles [47]. It is noteworthy that 18 unique alleles were identified in our maize landraces (almost 10% of total), confirming this percentage of unique alleles previously detected by Aci et al. [19]. Choukan et al. [48] observed 22% of unique alleles on 36 maize inbred lines from Iranby using 42 SSR, while Wasala and Prasanna [37] led to the identification of 31% unique alleles. The present study indicated a moderate differentiation among Saharan maize landraces related to such specific allele fixation in the core set. The total genetic differentiation (Fst = 0.33) was closed to that reveled in the first report across Algerian maize landraces (Fst = 0.31) [19]. Similar Fst values (ranged from 0.33 to 0.35) were previously reported in maize landraces from India, Switzerland and Northern America [38, 45, 49].

The population structure based on the Bayesian admixture-model is in agreement with PCA and cluster analyses as well as the geographic origin of landraces. This result appears in agree with previous studies [42, 50], which reported a lack of clustering based on phenotypes, environmental adaptation, grain color or type, and maturity.

Even if maize is an allogamous plant species, in agreement with our results, a higher genetic variation within population than among them was yet described in other studies [51]. Da Silva et al. [51] considered that SSR loci revealed higher heterozygosity within population, which produced low differentiation among them. As a result, a wide variation is engendered useful for crop improvement [52]. Cluster analysis revealed a significant genetic diversity within and among landraces at both phenotypic and genetic levels, but failed to cluster the landraces according to their geographical origin as previously observed by Cömertpay et al. [15] and Noldin et al. [53].

A negative correlation between morphological and genetic distance matrices was revealed by the Mantel test, in agreement with previous report on ryegrass varieties by Roldan-Ruiz et al. [54]. This suggests that phenotypic and genetic variances are due to independent genome regions [55], being molecular markers selectively neutral and not always mirrored the diversity of expressed traits [56]. These results justify the use of both molecular and phenotypic markers, resulting complementary for assessing maize genetic diversity.


The Algerian maize landraces hold a significant genetic diversity suggesting their potential adaptation to the extreme Sahara environment, according to genetic parameters and the clustering based on both molecular and agro-morphological markers. Four centuries of growing in the Sahara could lead to a widening of gene pool diversity in maize germplasm, due to genetic recombination, mutations, natural and human selection, as stated for plant populations by Hartl and Clark [57]. Natural selection could gathered favorable genes for tolerance to abiotic stress. Noteworthy, in the present study the best promising tolerant landraces to abiotic stress were identified with the aim to be included in breeding programs for the development of novel water use efficient and drought tolerant maize genotypes.


Plant material and agro-morphological descriptors

The maize germplasm used in this study consisted of forty-seven populations/landraces (a landrace of out-pollinated plant species can be defined “population”) collected during 2009 and 2010 in different Oasis of South-Western Algeria (from 27°52′ to 31°37’ N latitude and 0°17′ to 2°13’ W longitude). In particular, six and thirty-eight populations/landraces were collected from Bechar and Adrar oasis, respectively; whereas three populations/landraces were collected from Gherdaia (1) and Saida (2) (see Additional file 8: Table S1, also for the acronym of each population/landrace).

The trial was carried out at the experimental station of the National School of Agronomy of Algiers (36° 47’ N, 2° 03′E altitude 32 m) in the sub humid North of Algeria during the spring-2013 under field conditions. The experimental design was an Augmented design with three blocks [58] and seven checks of which 2 maize populations from Iowa State University, USA (BS17 and BSL(S)C6), one population from University of Guelph, Canada (Longfellow), and 4 from Spain (Norteño, Rastrojero and Tremesino populations from the Spanish Ministry of Agriculture, and Tuy from Misiòn Biològica de Galicia, Pontevedra) were repeated in each block.

Each experimental plot consisted of one row, 30 hills per row and one grain per hill. Rows were spaced 0.70 m apart and hills were 0.20 m. The data were recorded on ten random plants for twenty-four agro-morphological traits using the standard descriptors suggested by the International Board for Plant Genetic Resources [59] (Additional file 9: Table S2). They include plant emergence (EMG,%), number of days to reach 50% of field germination (T50, days), early vigor (EV), days to silking (DS, days), days to anthesis (DA, days), anthesis-silking interval (ASI, days), number of leaves (NL) and ears per plant (NEP), plant (PLH, cm)and ear height (EH, cm), ear row number (ERN), number of kernels per row (NKR), ear length (EL, cm) and diameter (ED, cm), cob (CD, cm) and rachis diameter (RD, cm), kernel length (KL, cm), width (KW, cm) and thickness (KT, cm), kernel proportion (K%), weight of 10 ears (EW, kg), 1000 kernel weight (1000 KW, g), moisture content at harvest (HMC, %), and kernel yield per plot (KYP, Mgha− 1).

DNA extraction

Maize seeds of each landrace were surface sterilized for 20 min in 20% (v/v) sodium hypochlorite solution, rinsed with deionized water and transferred for germination in Petri dish (diameter 9 cm), on filter paper with 0.5 mM CaSO4 at 26 °C in dark condition for one week.

The DNA pooled-sampling strategy was employed [60, 61] and the genomic DNA was isolated from 15 seedlings [37] for each population using DNeasy Plant Mini Kit (Qiagen, Milano, Italy), according to manufacturer’s protocol and its quality and quantification was assayed by BioPhotometer D30 (Eppendorf, Hamburg, Germany).

SSR analysis

Eighteen SSR loci were employed for genotyping the forty-seven landraces [19] ( (Additional file 10: Table S3). The PCR amplification was carried out with a Thermal Cycler 2720 (Applied Biosystem, Thermo Fisher Scientific Inc.) using specific temperature of annealing (Ta) for each primers pair. The reaction was carried out in 20 μl volume containing 20 ng DNA for each bulk, 1 U Taq DNA polymerase (Thermo Fisher Scientific Inc.), 0.32 μM reverse primer, 0.16 μM forward and 0.16 μM fluorescence (FAM) labeled universal primer M13 (− 21) as reported by [62]. PCR was programmed at 95 °C (5 min), 30 cycles 94 °C (30 s), Ta (45 s) and 72 °C (45 s) followed by 8 cycles 94 °C (30 s), 53 °C (45 s) and 72 °C (45 s), and final extension at 72 °C for 10 min. PCR products were separated by capillary electrophoresis and genotyped with an ABI PRISM 3500 Genetic Analyzer (Applied Biosystem). Allelic data were exported and selected from sequencer using Gene Mapper v.5 software (Applied Biosystem, Thermo Fisher Scientific Inc). Finally, to filter raw data, allele calling and estimate allele frequencies for each bulk and SSR locus the “FreqsR” software was used. In addition, the conversion of allele frequencies into allele sizes in individuals was carried out by “F-to-L” software [12, 16]. Both software were run on R platform (R Development Core Team 2008).

Statistical analysis

An augmented design was adopted to analyze the landraces for morphological traits [58]. Means of 47 landraces were adjusted by blocking effects of the replicated checks using ACBD-R software released by CYMMYT (International Maize and Wheat Improvement Center). Data were analyzed to find significant differences among landraces for each trait. Standardized means to remove the effects of different scales of measurements was adopted in the Principal Component Analysis (PCA) and Cluster Analysis (CA) As suggested by Clifford and Stephenson [32] and Guei et al. [63], the first three components were used for characterizing and differentiating the landraces. For cluster analysis, data were analyzed to determine Euclidean distance based on paired group method to determine dissimilar groups among landraces. Analyses were performed using R software (R Development Core Team 2008).

To determine genetic parameters such as number of alleles (A), average of different alleles (Na), effective number of allele (e), Nei’s gene diversity index (H), Shannon diversity index (I), forboth SSR loci and population GenAlex software version 6.3 was utilized [64]. Private alleles were determinate using GDA software (, whereas the Polymorphic information content (PIC) for each SSR locus was also estimated using Cervus v. 3.0.7 software (Copyright Tristan Marshal, Field Genetic, Ltd). The Analysis of Molecular Variance (AMOVA) was also performed to analyze genetic variation among and within individuals [65] by Arlequin software ( testing Fst by 9999 random permutations. A dendrogram to define differences among populations was constructed based on Nei and Li [20] pair-wise distances matrix and the unweighted pair group method of arithmetic clustering algorithm (UPGMA) [66] using MEGA v. 6 software [67].

In addition, model-based (Bayesian) clustering was performed to evaluate genetic relationship among individuals and population structure by using Software package STRUCTURE [68]. The program was set up and run as reported in Mercati et al. [69]. Then, the criterion (ΔK) of Evanno et al. [70] was used to determine the most probable K value, in order to compensate for overestimation of subgroup number by STRUCTURE. Finally, as reported by Wang et al. [71], lines with membership probabilities ≥0.90 were assigned to the corresponding subgroups and lines with membership < 0.90 were assigned to a mixed subgroup.



Analysis of Molecular Variance


International Board for Plant Genetic Resources


inbreeding coefficient within subpopulations relative to the total population (genetic differentiation)


Nei’s gene diversity index


principal component analysis


Polymerase chain reaction


Polymorphic information content


Simple sequence repeat


unweighted pair group method using arithmetic averages


  1. 1.

    Whitt SR, Wilson M, Tenaillon MI, Gaut BS, Buckler ES. Genetic diversity and selection in the maize starch pathway. Proc Natl Acad Sci U S A. 2002;99:12959–62.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. 2.

    Matsuoka Y, Vigouroux Y, Goodman MM, Sanchez Garcia J, Buckler E, Doebley J. A single domestication for maize shown by multilocus microsatellite genotyping. Proc Natl Acad Sci U S A. 2002;99:6080–4.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. 3.

    Doebley JF. The genetics of maize evolution. Ann Rev Genet. 2004;38:37–59.

    Article  PubMed  CAS  Google Scholar 

  4. 4.

    Revilla P, Soengas P, Cartea ME, Malvar RA, Ordás A. Isozyme variability among European maize populations and the introduction of maize in Europe. Maydica. 2003;48:141–52.

    Google Scholar 

  5. 5.

    Laumont P, Laby H. Le maïs et sa culture en Algérie. Doc et Rens agricoles. 1950; Bulletin n. 155, Alger.

    Google Scholar 

  6. 6.

    Chevalier A. Les Productions végétales du Sahara et de ses confins Nord et Sud. Passé-Présent-Avenir. Revue de Botanique Appliquée & D’Agriculture Tropicale. 1932;12:669–924. bulletin 133-134

    Article  Google Scholar 

  7. 7.

    Djemel A, Revilla P, Hanifi-Mekliche L, Malvar RA, Alvarez A, Khelifi L. Maize (Zea mays L.) from the Saharan oasis: adaptation to temperate areas and agronomic performance. Genet Res Crop Evol. 2011;59:1493–504.

    Article  Google Scholar 

  8. 8.

    Pressoir G, Berthaud J. Population structure and strong divergent selection shape phenotypic diversification in maize landraces. Heredity. 2004;92:95–101.

    Article  PubMed  CAS  Google Scholar 

  9. 9.

    Franco J, Crossa J, Ribaut JM, Betran J, Warburton ML, Khairallah M. A method for combining molecular markers and phenotypic attributes for classifying plant genotypes. Theor Appl Genet. 2001;103:944–52.

    Article  Google Scholar 

  10. 10.

    Smith JSC, Smith OS. Fingerprinting crop varieties. Adv Agron. 1992;47:85–140.

    Article  CAS  Google Scholar 

  11. 11.

    Senior ML, Murphy JP, Goodman MM, Stuber CW. Utility of SSRs for determining genetic similarities and relationships in maize using agarose gel system. Crop Sci. 1998;38:1088–98.

    Article  Google Scholar 

  12. 12.

    Warburton ML, Xianchun X, Crossa J, Franco J, Melchinger AE, Frisch M, et al. Genetic characterization of CIMMYT inbred maize lines and open pollinated populations using large scale fingerprinting methods. Crop Sci. 2002;42:1832–40.

    Article  Google Scholar 

  13. 13.

    Reif JC, Hamrit S, Heckenberger M, Schipprack W, Maurer PH, Bohn M, et al. Genetic structure and diversity of European flint maize populations determined with SSR analysis of individuals and bulks. Theor Appl Genet. 2005;111:906–13.

    Article  PubMed  CAS  Google Scholar 

  14. 14.

    Vigouroux Y, Glaubitz JC, Matsuoka Y, Goodman MM, Sánchez JG, Doebley J. Population structure and genetic diversity of New World maize races assessed by DNA microsatellites. Am J Botany. 2008;95:1240–53.

    Article  Google Scholar 

  15. 15.

    Cömertpay G, Baloch FS, Kilian B, Ülger AC, Hözkan H. Diversity assessment of Turkish maize landraces based on fluorescent labeled SSR markers. Plant Mol Biol Rep. 2012;30:261–74.

    Article  Google Scholar 

  16. 16.

    Dubreuil P, Warburton M, Chastanet M, Hoisington D, Charcosset A. More on the introduction of temperate maize into Europe: large-scale bulk SSR genotyping and new historical elements. Maydica. 2006;51:281–91.

    Google Scholar 

  17. 17.

    Warburton ML, Reif J, Frisch M, Bohn M, Bedoya C, Xia X. Genetic diversity in CIMMYT non-temperate maize germplasm: landraces, open pollinated varieties, and inbred lines. Crop Sci. 2008;48:617–24.

    Article  Google Scholar 

  18. 18.

    Liu Z, Guo R, Zaho J, Cai Y, Wang F, Cao M, et al. Genetic diversity of two important groups of maize landraces with the same name in China revealed by M13 tailed-primers SSRs. Agr Sci China. 2009;8:15–23.

    Article  CAS  Google Scholar 

  19. 19.

    Aci MM, Revilla P, Morsli A, Djemel A, Belalia N, Kadri Y, et al. Genetic diversity in Algerian maize (Zea mays L.) landraces using SSR markers. Maydica. 2013;58:304–10.

    Google Scholar 

  20. 20.

    Nei M, Li W. Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc Natl Acad Sci U S A. 1979;76:5269–73.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. 21.

    Carena MJ, Yang J, Caffarel C, Mergoum M, Hallauer AR. Do different production environments justify separate maize breeding programs? Euphytica. 2009;169:141–50.

    Article  Google Scholar 

  22. 22.

    Ignjatovic-Micic D, Ristic D, Babic V, Andjelkovic V. Vancetovic J. A simple SSR analysis for genetic diversity estimation of maize landraces. Genetika. 2015;47:53–62.

    Article  Google Scholar 

  23. 23.

    Zeven AC. Landraces: a review of definition and classification. Euphytica. 1998;104:127–39.

    Article  Google Scholar 

  24. 24.

    Bolaños J, Edmeades GO. The importance of the anthesis-silking interval in breeding for drought tolerance in tropical maize. Field Crops Res. 1996;48:65–80.

    Article  Google Scholar 

  25. 25.

    Ribaut JM, Hoisington DA, Deutsch JA, Jiang C, Gonzalez de Leon D. Identification of quantitative trait loci under drought condition in tropical maize. 1. Flowering parameters and the anthesis-silking interval. Theor Appl. Genetics. 1996;92:905–14.

    CAS  Google Scholar 

  26. 26.

    Ngugi K, Cheserek J, Muchira C, Chemining’wa G. Anthesis to Silking interval usefulness in developing drought tolerant maize. J Renew Agr. 2013;1:84–90.

    Google Scholar 

  27. 27.

    Asfaw Z. Relationship between spike morphology, hordeins and altitude within Ethiopian barley, Hordeum vulgare L. (Poaceae). Hereditas. 1989;110:203–9.

    Article  Google Scholar 

  28. 28.

    San Vicente FM, Hallauer AR. Mass selection for adaptation in Antigua maize (Zea mays L.) composite. J Iowa Acad Sci. 1993;100:9–12.

    Google Scholar 

  29. 29.

    Dewey DR, Lu KH. A correlation and path coefficient analysis of components of crested wheat grass seed production. Agron J. 1959;51:515–8.

    Article  Google Scholar 

  30. 30.

    Yücel C, Hizli H, Firincioglu HK, Cil A, Anlarsal AE. Forage yield stability of common vetch (Vicia sativa L.) genotypes in the Cukurova and GAP regions of Turkey. Turk J Agric For. 2009;33:119–25.

    Google Scholar 

  31. 31.

    Iqbal J, Shinwari ZK, Rabbani MA. Maize (Zea mays L.) germplasm agro-morphological characterization based on descriptive, cluster and principal component analysis. Pak J Bot. 2015;47:255–64.

    CAS  Google Scholar 

  32. 32.

    Clifford HT, Stephenson W. An introduction to numerical classification. New York: Academic Press; 1975.

    Google Scholar 

  33. 33.

    Gouesnard B, Dallard J, Panouillé A, Boyat A. Classification of French maize populations based on morphological traits. Agronomie. 1997;17:491–8.

    Article  Google Scholar 

  34. 34.

    Beyene Y, Botha AM, Myburg AA. A comparative study of molecular and morphological methods of describing genetic relationships in traditional Ethiopian highland maize. Afr J Biotech. 2005;4:586–95.

    CAS  Google Scholar 

  35. 35.

    Hartings H, Berardo N, Mazzinelli GF, Valoti P, Verderio A, Motto M. Assessment of genetic diversity and relationship among maize (Zea mays L.) Italian landraces by morphological traits and AFLP profiling. Theor Appl Genet. 2008;117:831–42.

    Article  PubMed  CAS  Google Scholar 

  36. 36.

    Sharma L, Prasanna BM, Ramesh B. Analysis of phenotypic and microsatellite-based diversity of maize landraces in India, especially from the north east Himalayan region. Genetica. 2010;138:619–31.

    Article  PubMed  CAS  Google Scholar 

  37. 37.

    Wasala SW, Prasanna BM. Microsatellite markers-based diversity and population genetic analysis of selected lowland and mid-altitude maize landraces accessions of India. J Plant Biochem Biot. 2012;22:392–400.

    Article  CAS  Google Scholar 

  38. 38.

    Freitag N, Schneider D, Mir C, Stamp P, Hund A, Messmer R. Swiss maize (Zea mays L.) landraces. Their genetic diversity and distinctiveness in a global comparison. Maydica. 2012;57:226–35.

    Google Scholar 

  39. 39.

    Oppong A, Bedoya CA, Ewool MB, Asante MD, Thompson RN, Adu-Dapaah H, et al. Bulk genetic characterization of Ghanaian maize landraces using microsatellite markers. Maydica. 2014;59:1–8.

    Google Scholar 

  40. 40.

    Laborda PR, Oliveira KM, Garcia AAF, Paterniani M, De Souza AP. Tropical maize germplasm: what can we say about its genetic diversity in the light of molecular markers? Theor Appl Genet. 2005;111:1288–99.

    Article  PubMed  CAS  Google Scholar 

  41. 41.

    Romay MC, Butrón A, Ordás A, Revilla P, Ordás B. Effects of recurrent selection on the genetic structure of two broad base Spanish maize synthetics. Crop Sci. 2012;52:1493–502.

    Article  Google Scholar 

  42. 42.

    Xia X, Reif J, Hoisington D, Melchinger A, Frisch M, Warburton ML. Genetic diversity among CIMMYT maize inbred lines investigated with SSR markers: I. Lowland tropical maize. Crop Sci. 2004;44:2230–7.

    Article  Google Scholar 

  43. 43.

    Liu K, Goodman M, Muse S, Smith JS, Buckler E, Doebley J. Genetic structure and diversity among maize inbred lines as inferred from DNA microsatellites. Genetics. 2003;165:2117–28.

    PubMed  PubMed Central  CAS  Google Scholar 

  44. 44.

    Bachtrog D, Agis M, Imhof M, Schlötterer C. Microsatellite variability differs between dinucleotide repeat motifs - evidence from Drosophila melanogaster. Mol Biol Evol. 2000;17:1277–85.

    Article  PubMed  CAS  Google Scholar 

  45. 45.

    Vigouroux Y, Jaqueth JS, Matsuoka Y, Smith OS, Beavis WD, Smith JSC, Doebley J. Rate and pattern of mutation at microsatellite loci in maize. Mol Biol Evol. 2002;19:1251–60.

    Article  PubMed  CAS  Google Scholar 

  46. 46.

    Beyene Y, Botha AM, Myburg A. Genetic diversity among traditional Ethiopian highland maize accessions assessed by simple sequences repeat (SSR) markers. Genet Res Crop Evol. 2006;53:1579–88.

    Article  CAS  Google Scholar 

  47. 47.

    Singode A, Prasanna BM. Analysis of genetic diversity in the north eastern Himalayan (NEH) maize landraces of India using microsatellite markers. J Plant Biochem Biot. 2010;19:33–41.

    Article  CAS  Google Scholar 

  48. 48.

    Choukan R, Hossinzade A, Ghanadha MR, Warburton ML, Talei A, Mohammadi SA. Use of SSR data to determine relationships and potential heterotic groupings within medium to late maturing Iranian maize inbred lines. Field Crops Res. 2006;95:212–22.

    Article  Google Scholar 

  49. 49.

    Gonzalez-Ugalde WG. Genetic characterization of northern American pointed Rice popcorns, a group commercially flints and flours maize (Zea mays L. ssp. mays) with isozyme, SSR and morphological markers. Dissertation: Iowa State University; 1997.

  50. 50.

    Xia XC, Reif JC, Melchinger AE, Frisch M, Hoisington DA, Beck D. Genetic diversity among CIMMYT maize inbred lines investigated with SSR markers: II. Subtropical, tropical mid altitude, and highland maize inbred lines and their relationships with elite U.S. and European maize. Crop Sci. 2005;45:2573–82.

    Article  CAS  Google Scholar 

  51. 51.

    Da Silva TA, Cantagalli LB, Saavedra J, Lopes AD, Mangolin CA, da Silva Machado M, et al. Population structure and genetic diversity of Brazilian popcorn germplasm inferred by microsatellite markers. Electron J Biotechnol. 2015;18:181–7.

    Article  Google Scholar 

  52. 52.

    Wietholter P, Cruz de Melo Sereno MJ, de Freitas Terra T, Delmar dos Anjos e Silva S, Barbosa Neto JF. Genetic variability in corn landraces from Southern Brazil. Maydica. 2008;53:151–9.

    Google Scholar 

  53. 53.

    Noldin O, Revilla P, Ordás B. Genetic diversity of the floury race of maize AvatiMorotî from the Guaraní tribe in Paraguay. Span J Agr Res. 2016;14:e0707.

    Article  Google Scholar 

  54. 54.

    Roldan-Ruiz I, van Eeuwijk FA, Gilliland TJ, Dubreuil P, Dillmann C, Lallemand J, et al. A comparative study of molecular and morphological methods of describing relationships between perennial ryegrass (Loliumperenne L) varieties. Theor Appl Genet. 2001;103:1138–50.

    Article  CAS  Google Scholar 

  55. 55.

    Navarro C, Cavers S, Pappinen A, Tigerstedt P, Lowe A, Merila J. Contrasting quantitative traits and neutral genetic markers for genetic resource assessment of Mesoamerican Cedrelaodorata. Silvae Genet. 2005;54:281–92.

    Article  Google Scholar 

  56. 56.

    Li F, Gan S, Wang Q, Zhao X, Huang S, Li M, et al. RAPD and morphological diversity among four populations of the tropical tree species Paramichelia aillonii (Pierre) Hu in China. Forest Ecol Manag. 2008;255:1793–801.

    Article  Google Scholar 

  57. 57.

    Hartl DL, Clark AG. Principles of population genetics. Sunderland: Sinauer Associates, Inc; 1997.

    Google Scholar 

  58. 58.

    Federer WT. Augmented designs with one-way elimination of heterogeneity. Biometrics. 1961;17:447–73.

    Article  Google Scholar 

  59. 59.

    CIMMYT, IBPGR. Descriptors of maize. CIMMYT/IBPGR. Rome; 1991.

  60. 60.

    Rebourg C, Gouesnard B, Charcosset A. Large scale molecular analysis of traditional European maize populations. Relationships with morphological variation. Heredity. 2001;86:574–87.

    Article  PubMed  CAS  Google Scholar 

  61. 61.

    Rebourg C, Chastanet M, Gouesnard B, Welcker C, Dubreuil P, Charcosset A. Maize introduction into Europe: the history reviewed in the light of molecular data. Theor Appl Genet. 2003;106:895–903.

    Article  PubMed  CAS  Google Scholar 

  62. 62.

    Schuelke M. An economic method for the fluorescent labelling of PCR fragments. Nat Biotechnol. 2000;18:233–4.

    Article  PubMed  CAS  Google Scholar 

  63. 63.

    Guei RG, Sanni KA, Abamu FJ, Fawole I. Genetic diversity of rice (Oryza sativa L.). Agron Africaine. 2005;5:17–28.

    Google Scholar 

  64. 64.

    Peakall R, Smouse PE. GenAlEx6: genetic analysis in excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6:288–95.

    Article  Google Scholar 

  65. 65.

    Excoffier LL, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992;131:479–91.

    PubMed  PubMed Central  CAS  Google Scholar 

  66. 66.

    Sneath PHA, Sokal RR. Numerical taxonomy: the principles and practice of numerical classification. San Francisco: W.H. Freeman and Company; 1973.

    Google Scholar 

  67. 67.

    Kumar S, Stecher G, Tamura K. MEGA/: Molecular evolutionary genetics analysis version 7.0 for bigger dataset. Mol Biol Evol. 2015;33:1870–4.

    Article  CAS  Google Scholar 

  68. 68.

    Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.

    PubMed  PubMed Central  CAS  Google Scholar 

  69. 69.

    Mercati F, Riccardi P, Leebens-Mack J, Abenavoli MR, Falavigna A, Sunseri F. Single nucleotide polymorphism isolated from a novel EST dataset in garden asparagus (Asparagus officinalis L.). Plant Sci. 2013;203–204:115–23.

    Article  PubMed  CAS  Google Scholar 

  70. 70.

    Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  71. 71.

    Wang R, Yu Y, Zhao J, Shi Y, Song Y, Wang T, et al. Population structure and linkage disequilibrium of a mini core set of maize inbred lines in China. Theor App Genet. 2008;117:1141–53.

    Article  CAS  Google Scholar 

Download references


The authors wish to acknowledge the Ministry of the high education and scientific research from Algeria that financed the project and the Università Mediterranea di Reggio Calabria, Italy for supporting some technical and financial aspects for the realization of this work.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding authors on reasonable request.

Author information




All the authors conceived and designed the experiments; MMA, AL, AM, and MA carried out the experiments; MMA, AL and FS performed the statistical analysis; MMA, AL, LK and FS commented the results and drafted the manuscript. All the authors read and approved the final manuscript.

Corresponding authors

Correspondence to Lakhdar Khelifi or Francesco Sunseri.

Ethics declarations

Ethics approval and consent to participate

The authors declare that the experiments of this study comply with the current laws. The populations used in the present research belong to a crop species not at risk of extinction. The maize populations utilized in the present study are available at LRGB, École Nationale Supérieure Agronomique (ENSA), Avenue Pasteur, Hassan Badi, El Harrach-Algiers 16,200, Algeria.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S4. Augmented ANOVA (mean squares) for 24 agro-morphological traits in 47 Algerian maize landraces and 7 checks. (DOCX 17 kb)

Additional file 2:

Table S5. Correlation matrix of 24 agro-morphological traits used to characterize 47 Algerian maize landraces. (XLSX 16 kb)

Additional file 3:

Table S6. Eigenvalues, variances and coefficients associated with first three principal components. (DOCX 15 kb)

Additional file 4:

Figure S1. Principal component analysis of the 47 Algerian maize landraces based on 24 agro-morphological traits. (JPEG 152 kb)

Additional file 5:

Table S7. Landraces summary statistics based on SSR analysis. (DOCX 14 kb)

Additional file 6:

Table S8. Analysis of Molecular Variance. (DOCX 12 kb)

Additional file 7:

Table S9. Partition of the gene diversity among clusters. (DOCX 12 kb)

Additional file 8:

Table S1. List of the 47 Algerian maize landraces. (DOCX 22 kb)

Additional file 9:

Table S2. Agro-morphological traits recorded in the 47 Algerian maize landraces. (DOCX 15 kb)

Additional file 10:

Table S3. Primers used to detect SSR marker. (DOCX 14 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

Aci, M.M., Lupini, A., Mauceri, A. et al. Genetic variation and structure of maize populations from Saoura and Gourara oasis in Algerian Sahara. BMC Genet 19, 51 (2018).

Download citation


  • Zea mays L.
  • Landraces
  • Genetic resources
  • Genetic diversity
  • Phenotypic variation
  • SSR markers