# Towards a comprehensive characterization of durum wheat landraces in Moroccan traditional agrosystems: analysing genetic diversity in the light of geography, farmers’ taxonomy and tetraploid wheat domestication history

Ali Sahri, Lamyae Chentoufi, Mustapha Arbaoui, Morgane Ardisson, Loubna

Belqadi, Ahmed Birouk, Pierre Roumet, Marie-Helene Muller

► To cite this version:

Ali Sahri, Lamyae Chentoufi, Mustapha Arbaoui, Morgane Ardisson, Loubna Belqadi, et al.. Towards a comprehensive characterization of durum wheat landraces in Moroccan traditional agrosystems: analysing genetic diversity in the light of geography, farmers’ taxonomy and tetraploid wheat domestication history. BMC Evolutionary Biology, 2014, 14, 10.1186/s12862-014-0264-2 . hal-01140157

HAL Id: hal-01140157

https://hal.science/hal-01140157

Submitted on 8 Mar 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Towards a comprehensive characterization of durum wheat landraces in Moroccan traditional agrosystems: analysing genetic diversity in the light of geography, farmers’ taxonomy and tetraploid wheat domestication history

Ali Sahri¹, Lamyae Chentoufi¹, Mustapha Arbaoui¹, Morgane Ardisson², Loubna Belqadi¹, Ahmed Birouk¹, Pierre Roumet² and Marie-Hélène Muller²*

**Abstract**

**Background:** Crop diversity managed by smallholder farmers in traditional agrosystems is the outcome of historical and current processes interacting at various spatial scales, and influenced by factors such as farming practices and environmental pressures. Only recently have studies started to consider the complexity of these processes instead of simply describing diversity for breeding purposes. A first step in that aim is to add multiple references to the collection of genetic data, including the farmers’ varietal taxonomy and practices and the historical background of the crop.

**Results:** On the basis of interview data collected in a previous study, we sampled 166 populations of durum wheat varieties in two traditional Moroccan agrosystems, in the Pre-Rif and Atlas Mountains regions. Using a common garden experiment, we detected a high phenotypic variability on traits indicative of taxonomical position and breeding status, namely spike shape and plant height. Populations often combined modern (short) with traditional-like (tall) statures, and classical durum squared spike shape (5 flowers/spikelet) with flat spike shape (3 flowers/ spikelet) representative of primitive domesticated tetraploid wheat (ssp. dicoccum). By contrast, the genetic diversity assessed using 14 microsatellite markers was relatively limited. When compared to the genetic diversity found in a large collection of tetraploid wheat, it corresponded to free-threshing tetraploid wheat. Within Morocco, the two studied regions differed for both genetic diversity and variety names. Within regions, neither geography nor variety names nor even breeding status constituted strong barriers to gene exchange despite a few significant patterns.

**Conclusions:** This first assessment of morphological and genetic diversity allowed pointing out some important factors that may have influenced the structure and evolutionary dynamics of durum wheat in Morocco: the significance of variety names, the occurrence of mixtures within populations, the relative strength of seed exchange between farmers and local adaptation, as well as the fate of modern varieties once they have been introduced. Further, multidisciplinary studies at different spatial scales are needed to better understand these complex agrosystems of invaluable importance for food security.

**Keywords:** Genetic diversity, Traditional agrosystems, Variety names, Landraces, Durum wheat

<footnote>* Correspondence: Marie-Helene.Muller@supagro.inra.fr
²INRA, UMR 1334, Amélioration Génétique et Adaptation des Plantes méditerranéennes et tropicales (AGAP), 2 place Pierre Viala, F-34060 Montpellier Cedex 1, France
Full list of author information is available at the end of the article</footnote>
## Background

With the Green Revolution, most of breeding effort has been devoted to the development of elite cultivars adapted to favorable lands where high input agriculture is predominantly practiced [1]. The attention dedicated to marginal lands is more recent [2,3]. These areas, characterized by specific and heterogeneous environmental conditions, reduced economic return and/or severe limitations for intensive agricultural use, host a large share (if not the majority in Africa, Asia, Central and South America) of rural populations. There, smallholder farmers practice low-input traditional agriculture and grow a large diversity of crop species [4]. One of the key production factors ensuring the sustainability of these agricultural systems is the set of traditional varieties (landraces) having evolved in these places for generations under a combination of environmental and human pressures [4]. Their high genetic diversity buffers against spatial and temporal heterogeneity and enhances resilience to abiotic and biotic stresses compared to modern varieties [5,6]. These landraces have often been considered as resources for contemporary agriculture; they have been used in breeding programs to enlarge the genetic diversity of modern genetic pools and to improve the adaptation of modern cultivars to abiotic or biotic constraints for instance, [7-9]. Considering the incapacity of modern breeding to provide adapted varieties to traditional agrosystems [10], landraces need to be considered for themselves [11], as the basis for food security. Efforts have to be devoted to their conservation and their use, especially in their place of origin [12].

Studies of landraces often focus on describing their genetic diversity in order to identify valuable and diversified genetic resources; the genetic information is generally only coupled with the geographical origin of the accessions [13,14]. To implement coherent actions for the conservation and use of the diversity present in traditional agrosystems, it is necessary to bridge the gap between the diversity observed within and across fields and the historical and current processes that have been shaping such diversity. These processes involve the interaction of multiple factors: Human practices such as farming practices, seed regeneration for the upcoming agricultural season, seed storage conditions, consumption habits, as well as seed exchanges within and between communities of farmers, generate selection pressures, drift or founder effects [15-17]. Particularly, cultural barriers act jointly with geographical ones on reproductive isolation, limiting or extending gene flow within a region [18,19]. In parallel, environmental conditions such as climate and soil characteristics are likely to strongly affect the evolutionary dynamics of landraces through heterogeneous selection pressures [20,21]. Moreover, extreme climatic events can lead to local extinction and seed renewing affecting

genetic diversity patterns [22]. Finally, the recent introduction of modern varieties might break the existing equilibrium, but can at the same time enrich the local genetic diversity [23].

Disentangling these mechanisms requires as first step integrating additional information to the morphologic and genetic characterization of landraces, in order to consider them in their environmental and social context [24]. For instance, various studies relying on interviews with farmers have been conducted, such interviews appearing essential to properly interpret the genetic diversity [11,25-27]. Together with the characteristics of the farm and the farmers' practices, the name of the variety especially appears as a crucial piece of information, since it is most often consciously used by farmers for management, selection, exchanges and uses [28].

Durum wheat (Triticum turgidum ssp. durum) is a selfing tetraploid cereal. The outcrossing rate is generally considered lower than 5% by breeders, and has been estimated between 1 and 3% in Ethiopian landraces [29]. Durum wheat history began with the domestication of wild emmer (Triticum turgidum ssp. dicoccoides) in the mountains of Fertile Crescent around the 8th millennium BC. Firstly it led to a non-brittle hulled subspecies (T. t. ssp. dicoccum) which was one of the first cereals domesticated [30]. Then, human selection of mutations at multiple loci in particular at Tg – tenacious glume – and Q loci, [31-33] gave rise to traits such as square-shaped spike, soft glumes and non-hulled grains improving threshing efficiency [34,35] and facilitating widespread cultivation. This phenotypic evolution together with hybridization between different forms [36] led to free-threshing subspecies (T. t. ssp. polonicum, T. t. ssp. carthlicum, T. t. ssp. turgidum, T. t. ssp. durum). Hulled and free-threshing forms (from 7th and 3rd millenium BP respectively) played a crucial role in the development of Mediterranean civilizations. They were widely spread out with the early agricultural movements leading to agriculture systems based on tetraploid wheat. Nowadays, elite durum varieties and landraces grown in different environments are of major importance for grain production in the Mediterranean basin [37]. Molecular markers-based studies showed that durum wheat history was associated with a decrease of the level of genetic diversity, from the wild ancestor to the most recent modern varieties [38,39], and that the different free-threshing subspecies, differing mainly in spike traits, were generally not distinguishable genetically [40,41].

In Morocco, tetraploid wheat represents a major crop (20% of cereals surface, www.agriculture.gov.ma) and an important source of staple food. Most of the cultivated varieties in the Great Plains result from breeding programs developed by the Moroccan National Institute
of Agricultural Research (INRA) or by international research centres (ICARDA, CIMMYT). These modern varieties are single-genotype based, early flowering spring (i.e. they do not require low temperature to induce flowering and/or they are insensitive to day length) with a short stature due to dwarfing genes, [42] to avoid lodging in relation with nitrogen fertilization. In contrast, in mountainous regions of Morocco (Rif, Central and High Atlas), the traditional agriculture system remains largely based on traditional durum wheat varieties, recognizable with their tall stature allowing the use of straw and grain. In these marginal areas, a diversity of tetraploid wheat types have been reported historically, including hulled wheat (denoted as subsp. *maroccanum* [43] or more classically as ssp. *dicoccum* [44,45]).

A few studies have described the genetic diversity of Moroccan landraces in terms of allele number and specific alleles [46,47]. Kehel et al. [46] detected a structure of the genetic diversity according to geography, water and temperature regimes; noticeably, landraces from the Rif and from the high continental plateaus from eastern High Atlas (including Saharan slope) were differentiated genetically. But in these studies, no other information than geographical location and large scale agro-ecological conditions was incorporated. Chentoufi et al. [48] recently conducted interviews with farmers from the Pre-Rif and the oases of the Atlas Mountains. They inventoried the varietal names and described their geographical distribution. Among their key results, they evidenced that the seed was in majority produced on farm, and that exchange practices included a combination on (mainly) purchase from local markets, supply from friends and relatives, and rarely, supply of modern varieties from agricultural extension offices. In these regions, modern varieties have indeed been introduced at variable levels and Chentoufi et al. [48] showed that they were integrated into the traditional practices of seed exchange and management. Nevertheless, farmers clearly distinguished traditional from modern varieties, and appreciated the traditional varieties especially for the straw and the quality for food processing.

In the present work, we enriched these interviews with genetic and morphologic data to address three questions: (i) How can we locate the Moroccan cultivated populations in the historical background of tetraploid wheat domestication? Do we find the different evolutionary steps described above (subspecies of *Triticum turgidum* and products of modern breeding)? (ii) How is the Moroccan genetic and morphologic diversity structured relative to two factors: the geographic location at small and large scales, and the farmers' taxonomy (i.e. the varietal names attributed to the cultivated populations). These two factors are potentially important determinants of gene flow through seed exchange or intercrossing. (iii) What are the consequences

of the introduction of modern varieties? Do the traditional and modern genetic pools remain distinct and how can the diversity of traditional varieties be affected on the long term?

To address these questions, we used two data sets. One, denoted as the Moroccan sample, is made of a sample of Moroccan farmers' population, identified by their varietal name and their farmer of origin [48]; the other, denoted as the INRA collection, is a worldwide collection of *Triticum turgidum* conserved at INRA (French National Institute for Agricultural Research) and including wild and domesticated subspecies. We considered three kinds of markers: The spike shape is an indicator of the taxonomical position, the stature of the plants is an indicator of breeding status, and a set of microsatellite markers gives us access to the neutral genetic diversity. The results are presented with respect with these three kinds of markers, and considered together to give insights into our three axes of investigation. For the interpretation, we also made use of the qualitative information gathered by the interviews of Chentoufi et al. [48].

## Results

We sampled 166 durum wheat populations in two Moroccan regions (Moroccan sample). These populations belonged to 30 different varieties that were classified as modern or traditional. Each region was subdivided into zones including geographically close villages (Table 1, Figure 1).

## Morphological diversity in the Moroccan sample

We characterized the height and spike shape of 30 maternal families (denoted as progenies) per population and of two Moroccan modern varieties used as checks.

### Spike shape variability: geographic structuration and relative inconsistency of its distribution across varietal names

Using a combination of spike shape and flow cytometric analysis, each progeny was assigned to one of four spike types: one corresponding to hexaploid bread wheat (*Triticum aestivum*, BW), and three corresponding to a shape indicative of different subspecies of the tetraploid species *Triticum turgidum* (CA for ssp. *carthlicum*, DC for ssp. *dicoccum* and DW for ssp. *durum*). Even if this description is referring to different subspecies of tetraploid wheat, it cannot be considered sufficient to confidently assign a plant to a subspecies. Variability for spike type was observed among and within populations. Based on within-populations frequencies of spike types, different categories of pure and mixed populations were defined (see Methods for details).
<table><caption>Table 1 Summary of the sampling design in Morocco</caption><thead><tr><th>Region</th><th>Zone</th><th>Altitude range (m)</th><th>Number of different varieties collected *</th><th>Number of sampled populations **</th></tr></thead><tbody><tr><td rowspan="5">Pre-Rif</td><td>ZN1</td><td>187-428</td><td>8 (1)</td><td>18 (1)</td></tr><tr><td>ZN2</td><td>191-422</td><td>5 (1)</td><td>10 (5)</td></tr><tr><td>ZN3</td><td>180-650</td><td>12 (3)</td><td>32 (8)</td></tr><tr><td>ZN4</td><td>872</td><td>2</td><td>5</td></tr><tr><td>total</td><td>180-872</td><td>18 (3)</td><td>65 (14)</td></tr><tr><td rowspan="7">Atlas Mountains</td><td>ZS1</td><td>1385-1464</td><td>5 (1)</td><td>15 (1)</td></tr><tr><td>ZS2</td><td>1593-2133</td><td>6</td><td>22</td></tr><tr><td>ZS3</td><td>2160-2387</td><td>5</td><td>27</td></tr><tr><td>ZS4</td><td>1662-1674</td><td>4</td><td>8</td></tr><tr><td>ZS5</td><td>1293-1296</td><td>3</td><td>6</td></tr><tr><td>ZS6</td><td>1492-1894</td><td>2</td><td>23</td></tr><tr><td>total</td><td>1293-2387</td><td>12 (1)</td><td>101 (1)</td></tr></tbody></table>

Altitude was recorded at the village and not at the field level.

Bread wheat was detected in 7 populations, all located in zone ZN1 of the Pre-Rif (Table 2, Figure 1A), constituting in two cases pure bread wheat populations.

Within *Triticum turgidum*, two spike types were predominant overall, DC and DW. In the Pre-Rif, ZN1 was characterized by a large diversity of spike types and a notable frequency of mixed populations (DW-DC and DW-DC-CA). The other zones were dominated by the classical DW type (Table 2, Figure 1A).

The DC and CA spike-types were always rare at the within-population level.

In the Atlas Mountains, spikes were mainly of DW and DC types. Two groups of zones appeared: ZS1, ZS2, ZS4 and ZS5 were characterized by a high proportion of the DW type. ZS3 and ZS6 were dominated by the DC type with similar overall frequencies, but with a prevalence of pure populations in ZS6. A gradient was detected along the Ziz valley with an increase of

<figure><img src="image_2.png" /></figure>

<footnote>* (resp.**) In brackets, number of modern varieties (resp. number of sampled populations belonging to modern varieties).</footnote>
<table><caption>Table 2 Frequencies of spike types and of short progenies within each zone</caption><thead><tr><th rowspan="2">Regions</th><th rowspan="2">Zones</th><th rowspan="2">Number of<br>progenies</th><th colspan="4">Spike types</th><th rowspan="2">Short<br>progenies (S)</th></tr><tr><th colspan="3">T. turgidum</th><th>T. aestivum</th></tr><tr><th></th><th></th><th></th><th>DW</th><th>DC</th><th>CA</th><th>BW</th><th></th></tr></thead><tbody><tr><td rowspan="4">Pre-Rif</td><td>ZN1</td><td>402</td><td>0.53</td><td>0.27</td><td>0.06</td><td>0.14</td><td>0.08</td></tr><tr><td>ZN2</td><td>299</td><td>0.89</td><td>0.03</td><td>0.08</td><td>0</td><td>0.53</td></tr><tr><td>ZN3</td><td>959</td><td>0.95</td><td>0.04</td><td>0.01</td><td>0</td><td>0.38</td></tr><tr><td>ZN4</td><td>149</td><td>0.9</td><td>0.1</td><td>0</td><td>0</td><td>0</td></tr><tr><td rowspan="6">Atlas Mountains</td><td>ZS1</td><td>450</td><td>0.81</td><td>0.19</td><td>0</td><td>0</td><td>0.07</td></tr><tr><td>ZS2</td><td>660</td><td>0.64</td><td>0.36</td><td>0</td><td>0</td><td>0</td></tr><tr><td>ZS3</td><td>810</td><td>0.11</td><td>0.87</td><td>0.02</td><td>0</td><td>0</td></tr><tr><td>ZS4</td><td>240</td><td>0.84</td><td>0.14</td><td>0.02</td><td>0</td><td>0.11</td></tr><tr><td>ZS5</td><td>180</td><td>0.98</td><td>0.02</td><td>0</td><td>0</td><td>0</td></tr><tr><td>ZS6</td><td>690</td><td>0.25</td><td>0.72</td><td>0.03</td><td>0</td><td>0.03</td></tr></tbody></table>

Frequencies have been computed among the progenies sampled in the field and sowed in the common garden.

DC type frequencies from ZS1 to ZS3. Mixed populations prevailed in the intermediate zone ZS2 by comparison with zones ZS1 and ZS3 where pure populations of respectively DW and DC spike-types were more frequent (Table 2, Figure 1A).

Most of the varieties were represented by several populations. We thus attempted to categorize the varieties according to the morphological composition of these populations (Tables 3 and 4). Some varieties were associated to a given

spike type: for instance, modern varieties were of the DW-type, as well as most traditional varieties of the Pre-Rif (e.g. Guemh, Krifla kehla, etc). In the Atlas Mountains, traditional varieties split up in two groups: DC varieties, including pure populations and mixed ones with a majority of DC spike-types, and DW varieties, also including pure populations and mixtures. But these categories reflected tendencies rather than strict correspondences: some varieties assigned to a given category also included outlier populations:

<table><caption>Table 3 Occurrence and morphological characteristics of the varieties cultivated in the Pre-Rif region</caption><thead><tr><th>Variety</th><th>Status</th><th>Sample size</th><th>Zones <sup>a</sup></th><th>Categorization of the variety</th></tr></thead><tbody><tr><td>Guemh</td><td>traditional</td><td>4</td><td>ZN3</td><td>Tall DW-type with impurities for stature</td></tr><tr><td>Guemh beldi</td><td>traditional</td><td>8</td><td><strong>ZN1, ZN2, ZN3</strong></td><td>No assignment. Great heterogeneity among populations</td></tr><tr><td>Guemh khel</td><td>traditional</td><td>1</td><td>ZN3</td><td>Short DW-type with impurities for stature</td></tr><tr><td>Guemh Ihmer</td><td>traditional</td><td>9</td><td><strong>ZN1, ZN2</strong></td><td>No assignment. Great heterogeneity among populations</td></tr><tr><td>Karim</td><td>modern</td><td>11</td><td><strong>ZN1, ZN2, ZN3</strong></td><td>Short DW type with impurities for spike morphology and stature. One outlier (mainly tall mixture DW_DC_CA)</td></tr><tr><td>Krifla beda</td><td>traditional</td><td>3</td><td>ZN4</td><td>Tall DW-type including mixtures with DC.</td></tr><tr><td>Krifla kehla</td><td>traditional</td><td>7</td><td><strong>ZN3, ZN4</strong></td><td>Tall DW-type including impurities for stature</td></tr><tr><td>Lehjaoui</td><td>traditional</td><td>1</td><td>ZN1</td><td>Tall mixture DW_DC</td></tr><tr><td>Le'khel</td><td>traditional</td><td>4</td><td><strong>ZN1, ZN3</strong></td><td>DW-type, with great heterogeneity among populations for stature</td></tr><tr><td>Massa</td><td>modern</td><td>1</td><td>ZN3</td><td>Short DW-type with impurities for stature</td></tr><tr><td>Mezrouba</td><td>traditional</td><td>3</td><td><strong>ZN3</strong></td><td>DW-type including mixtures of statures</td></tr><tr><td>Pedro</td><td>modern</td><td>1</td><td>ZN3</td><td>Short DW type with impurities for stature</td></tr><tr><td>Technique</td><td>traditional</td><td>2</td><td>ZN1, ZN2</td><td>Short DW type including mixtures of statures</td></tr><tr><td>Twinssia</td><td>traditional</td><td>1</td><td>ZN3</td><td>DW type, mixture of statures</td></tr><tr><td>Twinssia kehla</td><td>traditional</td><td>2</td><td><strong>ZN3</strong></td><td>No assignment. Tall and heterogeneity among populations for spike morphology</td></tr><tr><td>Zeriâa</td><td>traditional</td><td>4</td><td><strong>ZN3</strong></td><td>DW type, including mixtures with DC, and great heterogeneity among populations for stature</td></tr><tr><td>Zeriâa kehla</td><td>traditional</td><td>1</td><td>ZN1</td><td>Tall DW-type</td></tr><tr><td>Zeriâa twila</td><td>traditional</td><td>1</td><td>ZN2</td><td>Mixture of statures and DW_DC</td></tr></tbody></table>

Sample size: number of sampled populations.

<footnote><sup>a</sup> zones where populations of the variety have been sampled. In bold when the sample size was more than one in the zone.</footnote>
<table><caption>Table 4 Occurrence and morphological characteristics of the varieties cultivated in the Atlas Mountains region</caption><thead><tr><th>Variety</th><th>Status</th><th>Sample size</th><th>Zones <sup>a</sup></th><th>Categorization of the variety</th></tr></thead><tbody><tr><td>Aberyoun</td><td>traditional</td><td>15</td><td>ZS2, ZS3</td><td>Tall DC-type, including mixtures with DW.</td></tr><tr><td>Chgira Ibida</td><td>traditional</td><td>4</td><td>ZS1</td><td>Tall DW-type, including mixtures with DC and impurities for stature</td></tr><tr><td>Cocorit</td><td>modern</td><td>1</td><td>ZS1</td><td>Short DW-type</td></tr><tr><td>Ifermourgh</td><td>traditional</td><td>28</td><td>ZS2, ZS3, ZS6</td><td>Tall DC-type, including mixtures with DW. One DW outlier</td></tr><tr><td>Ilks</td><td>traditional</td><td>2</td><td>ZS3</td><td>No assignment. Tall mixture DW-DC</td></tr><tr><td>Irden taghezaft</td><td>traditional</td><td>6</td><td>ZS3, ZS4, ZS5</td><td>No assignment. Tall and great heterogeneity among populations for spike type.</td></tr><tr><td>Ibida touila</td><td>traditional</td><td>6</td><td>ZS1</td><td>Tall DW-type, including mixtures with DC and impurities for stature</td></tr><tr><td>Tabakhoucht</td><td>traditional</td><td>3</td><td>ZS2, ZS4</td><td>No assignment. Tall mixture DW-DC. One short outlier.</td></tr><tr><td>Taberyount</td><td>traditional</td><td>6</td><td>ZS2, ZS3, ZS5</td><td>Tall DC-type, including mixtures with DW. One DW outlier (ZS5)</td></tr><tr><td>Tamellalt</td><td>traditional</td><td>9</td><td>ZS1, ZS2, ZS5</td><td>Tall DW-type, including mixtures with DC.</td></tr><tr><td>Toumlilt</td><td>traditional</td><td>20</td><td>ZS1, ZS2, ZS4, ZS6</td><td>Tall DW-type, including mixtures with DC. One short outlier.</td></tr><tr><td>Zerbana</td><td>traditional</td><td>1</td><td>ZS4</td><td>No assignment. Tall mixture DW-DC</td></tr></tbody></table>

Sample size: number of sampled populations.

<sup>a</sup> zones where populations of the variety have been sampled. In bold when the sample size was more than one in the zone.

For instance, one population referred as Ifermourgh, a DC variety, contained only DW spike-type plants.

Some varieties couldn't be categorized, since a wide range of compositions were observed across populations of those varieties (Tables 3 and 4). Namely, most varieties of ZN1 (Guehm beldi and Guemh Ihmer) were represented by pure or mixed populations with various spike types, including sometimes bread wheat.

**Plant stature: the visible impact of modern breeding, even within traditional varieties**

By comparison with the height of the checks, the progenies were classified into one of two statures: tall (T), corresponding to the usual stature of traditional

varieties, and short (S), indicative of modern breeding (Figure 2). As for spike type, pure and mixed populations were both observed for stature.

The tall type T was the most frequent within each region and each zone, except in zone ZN2 in the Pre-Rif (Table 2). The Pre-Rif region, where more modern varieties have been recorded [48] and sampled (Table 1), was characterized by a higher proportion of short or mixed populations. In the Atlas Mountains, the vast majority of the sampled populations were made of tall individuals (Figure 1B).

The correspondence between status of the variety (modern vs traditional) and plant stature (short vs tall) was not strict (Tables 3 and 4). For instance, one population of

<figure><img src="image_3.png" /><figcaption>**Figure 2** Frequency distribution of the average height of the checks and of the progenies. The dashed line delimits the "short" (S) and "tall" (T) statures.</figcaption></figure>
the modern variety Karim was characterized by tall progenies (Table 3), whereas one population of Tabakhoucht, a traditional variety of the Atlas Mountains mostly contained short individuals (Table 4).

To account for the discrepancy between declared status and stature composition of the population, we will denote as modern-like (resp. traditional-like), a population containing a majority of short individuals (resp. tall individuals).

## Patterns of genetic diversity in *Triticum turgidum*: the INRA collection as a background for the Moroccan sample

One seed per Moroccan population (Moroccan sample, 161 individuals) and Moroccan check (2 varieties) was genotyped with 14 microsatellite loci, as well as 495 samples from a large collection of wild and domesticated subspecies of *Triticum turgidum* (INRA collection).

### Levels of genetic diversity

In total, 335 alleles were detected for the 14 microsatellite loci over the 658 individuals. The loci showed highly different levels of diversity, from locus Xgwm257 (5 alleles and $H_e = 0.327$) to Xgwm285 (47 alleles and $H_e = 0.870$ Table 5).

Among the subspecies of *Triticum turgidum* from the INRA collection, the wild form ssp. dicoccoides was the most diverse, followed by the primitive domesticated form ssp. dicoccum. The subspecies ssp. carthlicum and ssp. polonicum were the less diverse, while ssp. turgidum and ssp. durum showed intermediate values (Table 6). Within ssp. durum, we further distinguished a set of 161 varieties registered after 1950 (see Methods). When

<table><caption>**Table 5 Description and summary statistics of the 14 microsatellite loci used in the present study**</caption><thead><tr><th>Locus</th><th>Chromosome location</th><th>Allele size range</th><th>Number of alleles</th><th>H<sub>e</sub></th></tr></thead><tbody><tr><td>Xgpw7577</td><td>1B</td><td>214-235</td><td>9</td><td>0.440</td></tr><tr><td>Xgwm312</td><td>2A</td><td>184-291</td><td>42</td><td>0.785</td></tr><tr><td>Xgwm257</td><td>2B</td><td>180-194</td><td>5</td><td>0.240</td></tr><tr><td>Xgwm374</td><td>2B</td><td>191-232</td><td>19</td><td>0.327</td></tr><tr><td>Xgwm413</td><td>2B</td><td>84-123</td><td>15</td><td>0.850</td></tr><tr><td>Xgwm2</td><td>3A</td><td>211-257</td><td>19</td><td>0.451</td></tr><tr><td>Xgwm285</td><td>3B</td><td>207-336</td><td>47</td><td>0.870</td></tr><tr><td>Xgwm601</td><td>4A</td><td>144-176</td><td>16</td><td>0.763</td></tr><tr><td>Xgwm495</td><td>4B</td><td>161-189</td><td>17</td><td>0.857</td></tr><tr><td>Xgwm234</td><td>5B</td><td>216-262</td><td>26</td><td>0.841</td></tr><tr><td>Xgwm193</td><td>6B</td><td>166-206</td><td>23</td><td>0.859</td></tr><tr><td>Xgpw2103</td><td>7A</td><td>222-300</td><td>39</td><td>0.862</td></tr><tr><td>Xgwm297</td><td>7B</td><td>146-190</td><td>26</td><td>0.847</td></tr><tr><td>Xgwm537</td><td>7B</td><td>200-250</td><td>32</td><td>0.869</td></tr></tbody></table>

$H_e$: expected heterozygosity in the whole sample.

<table><caption>**Table 6 Genetic variability estimates in the subspecies of *Triticum turgidum* from the INRA collection and in the Moroccan sample**</caption><thead><tr><th>Subspecies of *Triticum turgidum* and Moroccan sample</th><th>Sample size</th><th>N<sub>a</sub></th><th>H<sub>e</sub></th><th>R<sub>s</sub></th></tr></thead><tbody><tr><td>ssp. dicoccoides</td><td>53</td><td>16.1</td><td>0.91<sup>a</sup></td><td>13<sup>a</sup></td></tr><tr><td>ssp. dicoccum</td><td>94</td><td>13.8</td><td>0.78<sup>b</sup></td><td>9.78<sup>b</sup></td></tr><tr><td>ssp. carthlicum</td><td>34</td><td>5.1</td><td>0.54<sup>ce</sup></td><td>4.69<sup>ce</sup></td></tr><tr><td>ssp. polonicum</td><td>29</td><td>4.8</td><td>0.50<sup>de</sup></td><td>4.47<sup>de</sup></td></tr><tr><td>ssp. turgidum</td><td>33</td><td>7.8</td><td>0.68<sup>bc</sup></td><td>7.15<sup>bc</sup></td></tr><tr><td>ssp. durum</td><td>252</td><td>8.6</td><td>0.56<sup>cd</sup></td><td>5.59<sup>cd</sup></td></tr><tr><td>Moroccan sample</td><td>161</td><td>6.5</td><td>0.49<sup>de</sup></td><td>4.57<sup>e</sup></td></tr></tbody></table>

$N_a$, mean number of alleles per locus; $H_e$, mean expected heterozygosity and $R_s$, mean allelic richness per locus computed for a minimum sample size of 17 diploid individuals per group.

Values with the same letters are not significantly different (P < 0.01).

Only 161 Moroccan individuals have been included in the final data set, due to the identification of 5 *Triticum aestivum* individuals in the original sample.

computed only over this set, the diversity of ssp. *durum* decreased to $H_e = 0.496$ and $R_s = 4.53$.

Overall, the Moroccan sample exhibited the lowest levels of genetic diversity. The level of genetic diversity was not significantly different between the two sampled regions ($H_e = 0.48$, $R_s = 4.69$ for the Pre-Rif and $H_e = 0.39$, $R_s = 3.45$ for the Atlas Mountains, P > 0.01). Most alleles identified in Morocco were included in the allele pool of the INRA collection: only 4 private alleles were detected among the 91 identified. These four alleles were present in the Pre-Rif samples at loci Xgpw7577, Xgwm312, Xgwm495 and Xgwm537, with a frequency ranging from 0.3% to 1.86%.

### Genetic structure of the INRA collection and its relationship with the Moroccan sample

Our aim was to describe the position of the Moroccan landraces in the background of durum wheat domestication history. For that purpose, we needed first to describe the genetic structure of the INRA collection with respect to its taxonomic and historic subdivisions: The domesticated subspecies of the INRA collection were analyzed using DAPC, a method that identifies clusters of genetically close individuals within a data set [49]. We then projected the Moroccan sample and the checks as supplementary individuals on the obtained clustering.

The choice of the optimal number of clusters (K) was tricky since there was no clear and repeatable value of K for which the Bayesian Information Criteria (BIC) value was lowest. Practically, we chose K = 10, a value from which the BIC decreased or increased only by negligible amounts in many runs of the analysis (Figure 3A). Note that slightly different values of K provided comparable summaries of the data and thus similar qualitative conclusions. Here, we only present the result of one run of DAPC (denoted as the reference run in the following) as
<figure><img src="image_4.png" /></figure>

well as the repeatability of the obtained clustering based on 15 additional runs with K = 10.

Among the 10 clusters identified in the reference run (Table 7, Figure 3B), three corresponded both to a homogeneous composition regarding the taxonomical position and to a relatively high average membership probability: a cluster with all but one individuals of *ssp. carthlicum* (C1), a cluster only made of *ssp. dicoccum* individuals (C10), and a cluster only made of *ssp. durum* (C9). Three more clusters corresponded to a relative homogeneous composition regarding the taxonomical position and/or to a relatively high average membership probability, although to a lower extent: a cluster almost only made of *ssp. dicoccum* individuals (C7), a cluster almost only made of *ssp. durum* (C3), and a more atypical cluster, including accessions of *ssp. polonicum* and *ssp. turgidum* and well distinguished from the others (C8). Finally, four clusters (C2, C4, C5 and C6) were a mixture of two to four different subspecies, and occupied the central part of the graph. This group of four clusters will be denoted as the “central pack” in the following.

Interestingly, the two homogeneous clusters of *ssp. durum* C3 and C9 included almost exclusively varieties registered from the 50ies, whereas all landraces were included in the central pack (data not shown).

Overall, the 15 additional runs gave similar qualitative results. The 6 clusters not belonging to the central pack were easily recognizable in the new clusterings, even if some were subdivided or merged in a few runs. A central pack was always present but with different subdivisions and a different number of clusters among runs. We thus decided to group some clusters into “metagroups” and evaluated the repeatability of these metagroups as the average proportion of individuals of the metagroup of the reference run, that was assigned to a similar metagroup in the new runs, or that clustered altogether to another metagroup in case of merging. We obtained a high repeatability score (>0.899, Table 7).

The Moroccan genotypes projected as supplementary individuals were mainly assigned to the four clusters of the central pack, with a few (11) assigned to the C9
<table><caption>Table 7 Results of DAPC applied to the collection of domesticated *Triticum turgidum*</caption><thead><tr><th>Subspecies</th><th>C1</th><th>C10</th><th>C7</th><th>C2</th><th>C4</th><th>C5</th><th>C6</th><th>C8</th><th>C3</th><th>C9</th></tr></thead><tbody><tr><td><i>carthlicum</i></td><td>33</td><td></td><td></td><td></td><td></td><td></td><td></td><td>1</td><td></td><td></td></tr><tr><td><i>dicoccum</i></td><td></td><td>38</td><td>30</td><td></td><td>18</td><td>2</td><td>5</td><td></td><td>1</td><td></td></tr><tr><td><i>polonicum</i></td><td></td><td></td><td></td><td>10</td><td>10</td><td>2</td><td>4</td><td>3</td><td></td><td></td></tr><tr><td><i>turgidum</i></td><td></td><td></td><td>1</td><td></td><td>15</td><td>2</td><td>7</td><td>8</td><td></td><td></td></tr><tr><td><i>durum</i></td><td></td><td></td><td></td><td>25</td><td>2</td><td>46</td><td>65</td><td></td><td>32</td><td>82</td></tr><tr><td><strong>Total</strong></td><td>33</td><td>38</td><td>31</td><td>35</td><td>45</td><td>52</td><td>81</td><td>12</td><td>33</td><td>82</td></tr><tr><td><strong>average membership probability</strong></td><td>0.999</td><td>1.000</td><td>0.962</td><td>0.959</td><td>0.953</td><td>0.918</td><td>0.961</td><td>0.986</td><td>0.925</td><td>0.967</td></tr><tr><td><strong>metagroups</strong></td><td>CARTH</td><td>DIC1</td><td>DIC2</td><td></td><td>CP</td><td></td><td></td><td>PT</td><td>DURUM</td><td></td></tr><tr><td><strong>repeatability</strong></td><td>0.984</td><td>0.998</td><td>0.899</td><td></td><td>0.974</td><td></td><td></td><td>0.978</td><td>0.956</td><td></td></tr></tbody></table>

Repartition of the accessions of the different subspecies among the 10 clusters identified by DAPC, and quality measures of the assignments to clusters and metagroups. C2, C4, C5, C6 are the 4 genetic clusters belonging to the 'central pack'.

Metagroups have been defined based on the repeatability of the groupings in 15 additional runs (see text). The repeatability corresponds to the average proportion of the individuals from the original metagroup that were assigned together to a similar metagroup in the repetitions.

cluster containing the most recent varieties (Figure 3C). The Moroccan plants assigned to the C9 cluster were mainly from modern varieties or modern-like populations. These assignments contrasted with the high frequency of spike-types indicative of ssp. *dicoccum* (Table 2). The unstability of the strucure of the central pack prevented us to be more precise on the historical position of the Moroccan landraces.

## Geographical patterns of genetic diversity of durum wheat in Morocco

### Levels of genetic diversity

Within regions, the zones differed only slightly for the amount of genetic diversity measured as $H_e$ and $R_s$: we detected 5 significant differences out of 42 tests. In the Pre-Rif Region, only ZN4 had a lower allelic richness compared to the three other zones. In the Atlas Mountains, ZS3 (at the highest elevation) and ZS5 exhibited the lowest level of genetic diversity, whereas ZS2 was the most genetically diverse zone. As expected for a selfing species, genotypes were highly homozygous. A few accessions were heterozygote at one locus at least: they represented 13% of our sample in both regions (Table 8).

### Genetic structure among geographic groups and among varieties

A hierarchical analysis of molecular variance (AMOVA) was performed with zones nested within regions. It detected significant differentiation both between regions ($F_{CT} = 0.146, p < 0.001$) and between zones within regions ($F_{SC} = 0.171, p < 0.001$). Most of the genetic variance was observed within zones (70.7%) whereas the remaining variance was almost equally partitioned among zones within regions (14.7%) and among regions (14.6%). One hundred and twelve different genotypes were identified in the Moroccan data set. Only three were shared between

the two regions: two were identified in modern-like populations and encompassed respectively 2 and 4 different varietal names (different between regions); the third was present in ten Atlas traditional populations (belonging to 4 different varieties) and in one population described as Krifla beda, a traditional variety of the Pre-Rif. We applied DAPC to the Moroccan sample: individuals from different regions were overall assigned to different clusters. Only 9 individuals were assigned to clusters made of individuals from another region (Additional file 1: Figure S1 and Additional file 2: Table S2). Six individuals from the Pre-Rif were assigned to the 4 Atlas-Mountains clusters (5.7% of the individuals of these clusters). Three individuals from the Atlas Mountains were assigned to the 5 Pre-Rif clusters (5.3% of the individuals of these clusters). These exceptions were in accordance with the pattern of

<table><caption>Table 8 Genetic diversity estimates for the Moroccan sample of durum wheat</caption><thead><tr><th>Regions</th><th>Zones</th><th>Sample size</th><th>N<sub>a</sub></th><th>H<sub>e</sub></th><th>R<sub>s</sub></th><th>Nb Het</th></tr></thead><tbody><tr><td rowspan="4">Pre-Rif</td><td>ZN1</td><td>13</td><td>3.5</td><td>0.412<sup>a</sup></td><td>2.75<sup>ab</sup></td><td>3</td></tr><tr><td>ZN2</td><td>10</td><td>3.1</td><td>0.418<sup>a</sup></td><td>2.67<sup>a</sup></td><td>1</td></tr><tr><td>ZN3</td><td>32</td><td>4.2</td><td>0.449<sup>a</sup></td><td>2.89<sup>ab</sup></td><td>3</td></tr><tr><td>ZN4</td><td>5</td><td>2.1</td><td>0.411<sup>a</sup></td><td>2.07<sup>b</sup></td><td>1</td></tr><tr><td rowspan="6">Atlas Mountains</td><td>ZS1</td><td>15</td><td>2.6</td><td>0.367<sup>ab</sup></td><td>2.23<sup>a</sup></td><td>0</td></tr><tr><td>ZS2</td><td>22</td><td>3.3</td><td>0.446<sup>a</sup></td><td>2.53<sup>a</sup></td><td>6</td></tr><tr><td>ZS3</td><td>27</td><td>3</td><td>0.334<sup>b</sup></td><td>2.19<sup>ab</sup></td><td>3</td></tr><tr><td>ZS4</td><td>8</td><td>2.5</td><td>0.426<sup>ab</sup></td><td>2.40<sup>a</sup></td><td>0</td></tr><tr><td>ZS5</td><td>6</td><td>1.8</td><td>0.314<sup>ab</sup></td><td>1.77<sup>b</sup></td><td>1</td></tr><tr><td>ZS6</td><td>23</td><td>3</td><td>0.356<sup>ab</sup></td><td>2.25<sup>ab</sup></td><td>3</td></tr></tbody></table>

$N_{a}$, mean number of alleles per locus; $H_e$, mean expected heterozygosity and $R_s$, mean allelic richness per locus computed for a minimum sample size of 5 diploid individuals per group. Nb Het, number of heterozygote genotypes. Values with the same letters are not significantly different ($P < 0.01$). Pairwise tests were performed within regions.
shared genotypes between regions described above. We can thus confidently go on by analysing separately the two regions.

We computed two kinds of $F_{ST}$ values between zones: one used the microsatellite genotypes of the samples; the other used as a genotype the name of the variety of origin of the sample. In the Pre-Rif, microsatellite-$F_{ST}$ values were roughly higher between more distant zones. All zones were significantly differentiated from each other, except ZN2 and ZN3 which, relative to the others, were also similar for the frequency of both short – modern-like – statures and spike types (Table 2, Table 9). The pattern was similar for $F_{ST}$ values based on variety names, even if the relationships was not stringent.

In the Atlas Mountains, two genetically distinct groups emerged from pairwise microsatellite-$F_{ST}$ comparisons between zones (Table 10). One group included ZS1, ZS2, ZS4 and ZS5 and the other ZS3 and ZS6. Indeed, zones of the first group were all significantly differentiated from zones of the other group. Moreover, no significant differentiation between zones was found within each group. This pattern mirrored the different spike-type frequencies measured in these two groups of zones (Table 2). By contrast, it was not coherent with the $F_{ST}$ values based on variety names. For instance, a relatively high differentiation based on variety names was detected between ZS3 and ZS6. This inconsistency suggests that the varietal taxonomy is a poor indicator of genetic differentiation.

Within each region, we built a neighbor joining tree based on the genetic dissimilarity between multi-locus genotypes. In both regions, there was no strict equivalence between variety name and multilocus genotype (Figures 4 and 5). Although some samples belonging to the same variety could be genetically identical or very close (e.g. Krifla kehla in the Pre-Rif Region, cluster B, Figure 4), some were very distant from each other (e.g. Guemh that appeared in cluster A and C, Figure 4). Reciprocally, some varieties appeared undistinguishable from each other: for instance, individuals belonging to the traditional varieties Ifermourgh and Aberyoun grown in the Atlas Mountains were intermixed (Figure 5, cluster B).

<table><caption>**Table 9 Pairwise differentiation indices between zones of the Pre-Rif**</caption><thead><tr><th>Studied areas</th><th>ZN1</th><th>ZN2</th><th>ZN3</th><th>ZN4</th></tr></thead><tbody><tr><td>ZN1</td><td></td><td><strong>0.165</strong></td><td><strong>0.185</strong></td><td><strong>0.256</strong></td></tr><tr><td>ZN2</td><td><strong>0.136</strong></td><td></td><td>NS</td><td><strong>0.170</strong></td></tr><tr><td>ZN3</td><td><strong>0.149</strong></td><td><strong>0.059</strong></td><td></td><td><strong>0.160</strong></td></tr><tr><td>ZN4</td><td><strong>0.318</strong></td><td><strong>0.281</strong></td><td><strong>0.150</strong></td><td></td></tr></tbody></table>

Above diagonal: pairwise genetic differentiation index ($F_{ST}$). Below diagonal: $F_{ST}$ based on variety names (see Methods). Only values significant at the 0.05 levels are reported. NS: non significant.

<table><caption>**Table 10 Pairwise differentiation indices between zones of the Atlas Mountains**</caption><thead><tr><th>Studied areas</th><th>ZN1</th><th>ZN2</th><th>ZN3</th><th>ZN4</th><th>ZN5</th><th>ZN6</th></tr></thead><tbody><tr><td>ZS1</td><td></td><td>NS</td><td><strong>0.243</strong></td><td>NS</td><td>NS</td><td><strong>0.150</strong></td></tr><tr><td>ZS2</td><td><strong>0.160</strong></td><td></td><td><strong>0.204</strong></td><td>NS</td><td>NS</td><td><strong>0.172</strong></td></tr><tr><td>ZS3</td><td><strong>0.265</strong></td><td><strong>0.205</strong></td><td></td><td><strong>0.215</strong></td><td><strong>0.306</strong></td><td>0.014</td></tr><tr><td>ZS4</td><td><strong>0.176</strong></td><td>NS</td><td><strong>0.257</strong></td><td></td><td>NS</td><td><strong>0.136</strong></td></tr><tr><td>ZS5</td><td><strong>0.184</strong></td><td><strong>0.170</strong></td><td><strong>0.262</strong></td><td>NS</td><td></td><td><strong>0.196</strong></td></tr><tr><td>ZS6</td><td><strong>0.439</strong></td><td><strong>0.301</strong></td><td><strong>0.332</strong></td><td><strong>0.443</strong></td><td><strong>0.530</strong></td><td></td></tr></tbody></table>

Above diagonal: pairwise genetic differentiation index ($F_{ST}$). Below diagonal: $F_{ST}$ based on variety names (see Methods). Only values significant at the 0.05 levels are reported. NS: non significant

High bootstrap values were only found for some terminal clusters, and the deep branching patterns were not well supported. However, some tendencies deserved to be highlighted, as they were partly related either to the status of the varieties or to some of their morphological categorization (Tables 3 and 4). In the Pre-Rif, five clusters could be distinguished (Figure 5). One cluster (E) grouped populations from ZN1, which often included different spike-types or statures. The four other clusters distinguished modern-like populations (namely, made of a majority of short progenies, clusters A and D), from traditional-like populations (namely, made of a majority of tall progenies, clusters B and C). The frequency of short or tall statures in the population was thus more discriminating than the declared status for interpreting the topology of the tree (Tables 3 and 4); for instance, one traditional-like population of the modern variety Karim appeared in cluster C.

In the Atlas Mountains, three clusters could be identified. One cluster included modern varieties or modern-like populations (C). The two others reflected a combination of the zone of origin (ZS3 and ZS6 for cluster B vs. ZS1, ZS2, ZS4 and ZS5 for cluster A) and of the spike-type class of the population (DC-type for cluster B vs DW-type for cluster A) (Figure 5). The correspondence between clusters and zone or spike-type was not strict, but this partitioning remarkably mirrored the one inferred from spike-type frequencies and from pairwise $F_{ST}$ values between zones (Table 2, Table 10).

## Discussion

### Morocco in the light of wheat history: a diversified morphological basis contrasting with a relatively reduced and unspecific neutral diversity

We identified in the Moroccan sample individuals and whole populations as bread wheat, even if these populations had been explicitly designated by farmers as durum wheat. In some countries, farmers consciously grow mixtures of durum and bread wheat [36,45] but such practices have not been mentioned here. We lack information to understand our observation. It could be partly
<figure><img src="image_5.png" /></figure>

explained by the fact that bread wheat was introduced
only a century ago in the area and is not yet well known
by the farmers. One of the consequence of these mixtures
it that they could allow interspecific hybridization, which
has already occurred in wheat history [36], and may lead
to original phenotypes and/or genotypes. The contribution
of interspecific hybridization to traditional varieties
diversity has for instance been evidenced between
Asian (Oryza sativa) and African rice (Oryza glaberrima)
in different African countries i.e. occurrence of varieties
described as 'hybrid', [50]. Interestingly, the ZN1 zone
where bread wheat was identified appeared in majority in
a distinct cluster in the NJ tree (even if no bread
wheat individual was included in the genetic dataset,
cluster E, Figure 4). But this is an insufficient indication to
conclude on interspecific hybridization.

Plants with ancestral spike shape of the dicoccum type
were observed within Moroccan landrace populations
but they were free-threshing at maturity. Moreover, they

did not cluster with the specific ssp. dicoccum clusters of
the collection (C7, C10, Figure 3C). In other countries,
free-threshing dicoccum landraces have been described,
as well as durum landraces partially introgressed with
dicoccum genes resp. in Yemen and Ethiopia, [45].
We have no information on present-day emmer wheat
(ssp. dicoccum) cultivation in Morocco, but it was
documented during the first part of the twentieth
century [45]. Here, the dicoccum-like landraces could
result from recombination between different ancestral
gene pools. Whatever their origins, the reason why
such forms are preferentially cultivated in some zones
(for instance for their local adaptation or for some
particular uses) and whether they relate to more difficult
post-harvest processing deserve further investigation.
According to the interview data [48], even though
farmers identified some varieties as easier to thresh, these
varieties were most often the modern ones. Using direct
measurements, variability has already been evidenced for
<figure><img src="image_6.png" /></figure>

threshing efficiency within a diversified sample of free-threshing durum wheat landraces [35]. Such measurements could be applied to our sample.

As a whole, the high morphological diversity in the Moroccan tetraploid wheat sample did not correspond to a high and/or original genetic diversity relative to the *Triticum turgidum* INRA collection. Moroccan genetic diversity was relatively reduced and similar to that of the group of modern varieties from the INRA collection registered since the 1950ies. Such result is however not surprising considering the restricted geographical range of our sampling. Oliveira et al. [41] also identified North-West Africa as the Mediterranean region with the least genetic diversity for durum wheat. The Moroccan genetic diversity has also been shown to be lower than for Syrian landraces [46]. Moreover, except for ssp. *carthlicum* which is supposed to result from hybridization with bread wheat, [36], the free-threshing subspecies did not correspond to genetically distinguishable entities

also as previously shown, [41] and were combined in a metagroup with no clear or consistent composition (metagroup CP, Table 7, Figure 3B): almost all Moroccan individuals were assigned to that metagroup. More information on the historical origin of Moroccan durum wheat diversity appears then difficult to gather with our sample, considering the complex and still unresolved history of durum wheat across the Mediterranean basin long period of cultivation, repeated introductions, and gene flow between the different origins, [41].

### Geographical and varietal differentiation in Morocco: large-scale patterns but loose local barriers to gene flow and mixing

Chentoufi et al. [48] inventoried the varieties grown by farmers and showed that except for modern varieties, no names were shared between the two study regions. Such differentiation was supported by our genetic results. Indeed, the Pre-Rif and the Atlas Mountains samples
were strongly differentiated from each other. We thus have
two regions that kept their originality, which con-
firmed the large-scale geographic structure detected
by Kehel et al. [46]. Such high differentiation is not
surprising given that these two distant regions display
very different environmental, cultural and agronomic
conditions. This is also a common result in the study
of diversity of traditional varieties, as exemplified by
sorghum [18]. However, in some countries, studies have
shown a lack of geographical structure despite a strong
environmental heterogeneity, a pattern attributed to large-
scale seed exchange and continuous seed introduction
e.g. barley in Ethiopia, [51]. This is apparently not
(yet) the case in Morocco.

Within regions, although variety names were rather
localized [48], only a few patterns emerged from the
genetic and morphologic characterization, as if neither
distance nor different variety names constituted consistent
barriers to genetic exchange. Moreover, 48% of popu-
lations combined different statures and spike types
(40% when considering only spike type); some populations
even had a radically different composition than the other
populations within their varieties. Such inconsistencies
have been evidenced recurrently, at various levels of inves-
tigation e.g. [52]. They illustrate the limits of the use of
variety names to assess genetic and morphological
diversity [53] and to elaborate sampling strategies for
the conservation of genetic resources [28].

Besides this, we can point out some mechanisms
potentially involved in these patterns (and lacks of pattern)
and open tracks of future research to decipher them.

First, the large proportion of mixed populations
(on the base of spike type and stature) suggests the
occurrence of intermixing and genetic exchanges between
varieties with different morphological composition (DC or
DW spike-types, and short or tall statures). These
exchanges might have different sources: unconscious
mixing in threshing areas used by different farmers
and for different varieties, unreliability of seed exchange
networks (e.g. markets providing seed lots that do not cor-
respond to the declared names) [54], limited farmers'
interest for durum wheat cultivation and seed production.
These mixtures create strong opportunities to generate
diversity through cross-fertilization (heterozygote geno-
types were observed) and recombination. One question of
interest is whether these mixtures will homogenize the
pool of traditional varieties, or whether mechanisms such
as human practices for instance choice of less admixed
plant for seed production [55] or divergent environmental
pressures for instance, when populations of a given
landrace are grown in different environments, [21,56]
will maintain some differentiation between them.

Second, the taxonomy of varieties and how farmers
experience this taxonomy is questioned by our results.

Indeed, our morphological and genetic data suggest that
names of some varieties are interchangeable (Table 4
and Figure 5), for instance Ifermourgh and Aberyoun in
the Atlas Mountains. By contrast, some names or groups
of names carry information on genetic distinctiveness: for
instance Krifla Kehla in the Pre-Rif (cluster B, Figure 4)
and the two differentiated groups of traditional varieties in
the Atlas Mountains (i.e. landraces corresponding to DC
and DW spike types). (i) First, if these names exist, there
is probably a reason to which we do not have any access
yet. We need to understand which criteria farmers use to
distinguish and name their traditional varieties. The inter-
views described in Chentoufi et al. [48] only identified
criteria distinguishing modern and traditional varieties.
The use of less oriented protocols should allow the
expression of the plants’ and populations’ characteristics
actually recognized and used by the farmers themselves
[57]. (ii) Second, if there are criteria, the question is then
whether all the farmers still know and use them. Interest-
ingly, the interviews have shown a rapid recent renewal of
varieties at the farm level in the Pre-Rif region, possibly
for climatic reasons (crop failures due to limited rainfall)
[48]. Such phenomenon has also been observed for barley
in this region [26]. Farmers from the Pre-Rif region regu-
larly have to outsource their seeds and often grow differ-
ent varieties over years: this may reduce the knowledge
associated with each variety and allow mixing and errors,
depending on the reliability of the seed exchange
networks. In this region, except for a few groupings, we
could not discriminate traditional varieties based on their
genotypes. In contrast, the two genetically and morpho-
logically well-differentiated groups of varieties in the Atlas
Mountains coexist even when cultivated in the same area
(ZS6). This suggests that here, farmers’ practices maintain
these two groups well-differentiated. Such variable import-
ance given to a crop in different areas, associated with
more or less frequent crop failure and seed renewal has
been evidenced by Samberg et al. [24] for barley in
Ethiopia. Investigations on the practices of selection and
choices, in the field or during seed exchange, are required
to understand these contrasting patterns.

Finally, we have shown that the two well-differentiated
groups of varieties of the Atlas Mountains partially
overlapped with the two geographic groups of zones
(ZS3-ZS6 vs the four other zones). Assessing the relative
effect of local adaptation and gene flow would help
understanding this pattern. For instance, local adaptation
may arise from different selective pressures along the
altitudinal gradient in the Ziz valley [20]. Reciprocally,
informal discussions suggest that farmers from ZS3
regularly renew their seeds from the much lower altitude
zone ZS6 following crop failures in cold years: such
preferential exchanges might explain the differentiation of
these two zones relative to the others.
## The signature of the incorporation of modern varieties into the traditional network

Modern varieties have been introduced in Moroccan agrosystems in the 1970ies, during the Green revolution. Chentoufi et al. [48] showed that they were particularly present in the Pre-Rif region and that they were submitted to the same agricultural practices than the traditional varieties [48]. The present study showed that the cultivation of modern varieties resulted in gene flow between modern and traditional varieties. Indeed, some populations of traditional and modern varieties combined short (indicative of modern breeding) and tall plants (illustrating the double use of traditional varieties, for straw and grain). Moreover, identical, or similar genotypes were shared between varieties of different status.

Such exchanges have occurred repeatedly in cultivated plants whatever their mating system. For instance, introgression from modern into traditional varieties has been evidenced in maize (outcrossing) and barley (selfing) [58,59]. In Mexico, the introduction of modern varieties of maize into traditional agrosystems gave rise to what is called “creolized” varieties [23]. In these cases, landraces often maintained their genetic identity, and an increase in genetic diversity could even be observed at the agrosystem level due to the introduction of new alleles. However, whenever the pressure of modern varieties becomes too strong, i.e. when they become more prevalent than the traditional ones, an overall loss of diversity is expected [60]. In the Pre-Rif region, the surveys performed by Chentoufi et al. [48] revealed that the cultivation of modern varieties increased in the previous years, and that traditional varieties even disappeared from some villages. Together with our results on introgression, this claims for measures dedicated to the conservation of the traditional varieties.

## Conclusions

Our study, together with the interviews reported in Chentoufi et al. [48], provides an essential first assessment of the varietal, phenotypic and genetic diversity present in Moroccan agroecosystems. Despite our limited level of investigation (namely large-scale geographical sampling, with only a single genotyped individual per population), we interpreted our results in the light of data collected from farmers' interview, and identified a number of investigation tracks to arrive at a good comprehension of the processes generating diversity in these traditional agrosystems, and of the factors that might affect their future evolution.

Following these different tracks requires to concentrate investigations at a small scale, potentially different for each question and to integrate contributions from different disciplines [24]: consideration of the socio-cultural context [16], investigation on the seed exchange practices [27,54],

sampling simultaneously intrapopulation and interpopulation diversity for some landraces [19], considering the genetic diversity present on local markets [24,26], investigating more detailed selection, threshing or other practices [55,61], and collecting precise environmental variables on the study sites [21].

Climatic, economic and political contexts have been changing rapidly in traditional agrosystems during the last few decades following for instance the development of communication routes and agricultural policies, as described for Ethiopia in Samberg et al. [27]. Understanding the current dynamic of traditional agrosystems, by reciprocal learning between farmers and scientists [4], is the necessary step for conservation and use of local diversity in the face of these changes.

## Methods

### Study areas and sample collection

Chentoufi et al. [48] surveyed the durum varieties cultivated in two contrasted agro-ecological regions: the Pre-Rif region in the North of Morocco (medium mountains, rain-fed agriculture) and the oases of the Atlas Mountains (marked altitudinal gradient, irrigated agriculture) in the South-East. The Pre-Rif area is located at the interface between plains were agriculture has been strongly modernized and mountains were traditional agriculture still prevails. It was subdivided into four zones on the basis of geography. The oases of the Atlas Mountains remained more isolated from the other parts of the countries due to transport difficulties in mountainous areas; this region subdivided into six zones on the basis of altitudinal range and geography (Figure 1). Each cultivated variety was identified by the name provided by the farmer: 19 different names were recorded in the Pre-Rif from 163 farmers' interviews, and 14 in the Atlas Mountains from 101 farmers' interviews. More details on the study areas and on the varietal survey can be found in Chentoufi et al. [48]. Varieties were classified as traditional or modern by the farmers. Because a given variety can be grown in different farms, we denote by “population” the particular seed lot of a variety actually sown by a given farmer. Farmers don't grow multiple varieties in the same field: a field then corresponds to a population.

Based on these results, we followed two different sampling strategies in the two regions. In the Pre-Rif region, we adopted a sampling strategy allowing the collection of as many as possible of the wheat varieties. In total, 65 durum wheat populations were sampled. These populations corresponded to 18 different varieties: 17 identified during the initial survey of Chentoufi et al. [48] and one additional variety identified during our collecting trip. In the oases of the Atlas Mountains, we sampled one population from each of the surveyed farmers. In total, we
sampled 101 wheat populations in this region, corresponding to 12 different varieties (Table 1, Figure 1). No names were shared between the two regions.

For each of the 166 wheat populations, 30 spikes from 30 randomly chosen plants were harvested in the farmer's field.

## Morphological characterization

During the 2011/2012 growing season, 20 seeds from each individual spike were sown in an ear-row design (1 m length and 0.3 m between rows) at the experimental station of the IAV Hassan 2 (Rabat, N 33°58'35", W 6°51'59", 50 m above sea level). Two Moroccan durum wheat modern varieties (Karim and Marzak, seeds provided by Moroccan National Institute for Agricultural Research) were used as checks to account for potential environmental variation at the field station: 74 rows of Karim and 39 of Marzak were sown across the experiment (20 seeds per row).

In the following, we will denote as a progeny the material family corresponding to the seeds from a given spike.

We recorded two key morphological traits that are informative on durum wheat history, namely spike shape and plant height. We classified progenies into four categories of spike: (i) dicoccum-like spikes, i.e. flat ear with only three flowers per spikelet (denoted as DC); this corresponds to the shape expected for plants carrying the *q* allele of the *Q* gene [34] but note that these spikes have not systematically a hulled phenotype at maturity; (ii) lax and flat spikes with long rachis and elongated rachilla internodes. Because this morphology can be indicative of *T. turgidum* ssp. *carthlicum* or of bread wheat *T. aestivum*, species status was checked by flow cytometric analysis for all progenies of the populations where it was observed protocol adapted from [62]. Hexaploid progenies were classified as bread wheat (denoted as BW) and tetraploid ones as carthlicum type (denoted as CA); and (iii) classical square-shaped durum wheat spike type (compact spike with 5 flowers per spikelet, denoted as DW, corresponding to the shape of *Q* allele-carrying plants [34]). These four categories (DC, BW, CA and DW) will be denoted as spike type in the following.

Average plant height was recorded for each progeny by averaging 3 measures taken across the row (Figure 2). For the checks, it ranged from 55 to 90 cm; this range of values was considered as representative of modern varieties where dwarfing genes have been introduced [42]. In the Moroccan sample, the distribution was clearly bimodal, the first mode matching the checks' distribution. The upper limit of the checks' distribution was then used to assign each progeny to one of two plant statures: short (S, average height less or equal to 90 cm, possibly carrying dwarfing genes), or tall (T, average height larger than 90 cm).

## DNA extraction and genotyping procedure

For each population composed of thirty collected spikes, we randomly picked one spike. One seed per spike and one for each Moroccan check were sown for DNA extraction.

Genomic DNA was extracted from 100 mg of fresh leaf tissue of three weeks old seedlings. The extraction was performed following a protocol adapted from Tai and Tanksley (1991). DNA quality was checked on agarose gel and the DNA concentration was estimated using spectrofluorometry. Five samples of the Pre-Rif belonged to progenies identified as hexaploid wheat (BW type, see above) and were excluded from the genetic analysis. Thus, 163 Moroccan individuals were included in the final data set: 161 populations (Moroccan sample) and two checks.

In addition, we analyzed a set of 495 DNA samples of genotypes from a large collection representing different subspecies of *Triticum turgidum* conserved at INRA (INRA collection). It included 53 accessions of the wild ancestor *T. turgidum* ssp. *dicoccoides*, 94 accessions of ssp. *dicoccum*, 34 of ssp. *carthlicum*, 29 of ssp. *polonicum*, 33 of ssp. *turgidum* and 252 of ssp. *durum* subsets of this sample have been analyzed in [39,63]. Available information on the ssp. *durum* accessions allowed distinguishing among them 69 landraces from Lebanon, Jordan, Syria and Turkey and 161 varieties of durum wheat registered after 1950 (122 from France and 39 from other countries).

Fourteen nuclear SSR markers used in previous assessments of durum wheat diversity and mapping to a single, genome specific, locus sequences available at http://wheat.pw.usda.gov/GG2/index.shtml, [38,61] were selected. Polymerase chain reaction (PCR) amplification was performed in a total volume of 20 µl containing 50 ng of genomic DNA, 2 nM MgCl₂, 0.2 mM of each dNTP, 0.5 units of Taq DNA polymerase and 1-4 pM of each primer. The forward primer was labeled with one of the three fluorophores (6FAM, NED or HEX). PCR was carried out as follows: after 5 min at 94°, 35 cycles were performed with 30 s at 94°C, 30 s at either 55 or 60°C (depending on the locus), and 1 min at 72°C, followed by final extension step of 30 min at 72°C. Amplified products were detected on an ABI 3130xl Genetic Analyser (Applied Biosystems, Foster City, CA, USA) and analyzed using the GENEMAPPER V3.7 software (Applied Biosystems).

## Morphological data analysis

Frequencies of spike types and of plant statures among progenies were computed for each population and each zone.

Then, frequencies within populations were used to assign each population to a class. First, for each of the two morphological traits, we denoted as pure populations the populations where all sampled progenies were of the same
type (i.e. spike type or stature) with the allowance of one outlier, and as mixed populations the ones combined different types, with at least two progenies of each. For spike types, we defined five classes: (i) BW, corresponding to populations containing bread wheat (*T. aestivum*) individuals. This class included pure populations of bread wheat and mixed populations with a majority of bread wheat progenies associated with *T. turgidum*, (ii) DW and (iii) DC, corresponding to pure populations of respectively DW and DC spike type, (iv) DW\_DC, corresponding to mixed populations combining DC and DW spike types, and (v) DC-DW-CA, corresponding to mixed populations including these three spike types. Note that the CA type was always rare within populations.

For plant stature, we defined five classes of populations: S and T, corresponding respectively to pure populations of short or tall progenies, and three kinds of mixed populations ($S < T, T = S$, and $S > T$, when respectively the frequency of short progenies was lower than 33%, between 33 and 66%, and higher than 66%).

## Genetic diversity and population structure

The mean number of alleles per locus ($N_a$), mean allelic richness $Rs$, [64], and expected heterozygosity $He$, [65] were computed using FSTAT 2.9.3.2 software [66] in different groups of the dataset, namely subspecies of *T. turgidum*, whole Moroccan sample, regions within Morocco, and zones within regions. Differences in the amount of genetic diversity between groups were tested by Wilcoxon signed-rank tests, comparing values for the same loci in different groups, using R (R Development Core Team 2008). For these tests, we used the significance threshold $P = 0.01$.

Population structure within the INRA collection of *T. turgidum* was examined by applying the discriminant analysis of principal components DAPC, [49], a multivariate method designed to identify and describe clusters of genetically related individuals. The method relies on allele data transformation using principal component analysis as a prior step to discriminant analysis. DAPC was performed using the adegenet package [67] in R (R development Core Team). DAPC is particularly useful for identifying population structure without assuming a population model that would not be supported in our data set.

We first ran DAPC on the INRA collection, including only the domesticated subspecies of *T. turgidum*. The optimal number of clusters was determined using the *find.clusters* function which implements successive K-means clustering. The rate of decrease of the Bayesian information criterion (BIC) was visually examined, and the number of clusters was determined as the value of *K* above which BIC values decreased only subtly. In order to obtain reliable group membership

probabilities, the number of Principal Components (PC) retained to compute the discriminant functions was determined using the function *optim.a.score*, which examines the trade-off between power of discrimination and over-fitting. Here we retained 11 to 13 PC depending on the runs. We then projected the Moroccan genotypes as supplementary individuals onto the discriminant functions and examined their assignment to the different clusters. A DAPC analysis was run on the Moroccan sample alone, using the same process as described for the INRA collection.

Within Morocco, a hierarchical analysis of molecular variance was performed with zones nested within regions, using ARLEQUIN v 3.5 [68]. The significance of the fixation indices and of the variance components was assessed using 10 000 permutations.

Genetic differentiation between zones within region ($F_{ST}$) was estimated using the method of Weir and Cockerham [69], with the software Genetix V4.05 [70]. To quantify how much zones differed for the varieties they host, we also computed $F_{ST}$ estimates based on variety names: Namely, each individual was described by its variety name, which was considered as equivalent to a homozygous genotype at a single locus. These $F_{ST}$ values were compared with the classical $F_{ST}$ values based on microsatellite genotypes. It is a way to assess the discrepancy between the information provided by the genotype and by the name of the variety of origin of a sample.

In order to depict the genetic relationships between populations, DARwin version 5.0 was used to compute a genetic dissimilarity matrix using the simple matching index [71]. The dissimilarity matrix was then subjected to cluster analysis using the weighted neighbour-joining method to obtain a dendrogram for the Pre-Rif accessions and the Atlas Mountains populations. The uncertainty of the dendrogram structure was assessed with 1000 bootstraps.

## Additional files

**Additional file 1: Figure S1.** Results of DAPC applied to the Moroccan sample.A. Bayesian Information Criteria (BIC) for increasing values of the number of clusters. The chosen number of clusters was $K = 9$. B. Scatterplot of the first two principal components of the DAPC on the sample of Moroccan populations. Individuals are represented by symbols according to their region of origin. Numbers, colours and inertia ellipses identify the clusters. The bottom-right inset shows the eigenvalues of the 6 first principal discriminant functions.

**Additional file 2: Table S1.** Results of DAPC applied to the Moroccan Sample Repartition of the Moroccan populations from the different regions among the 9 clusters identified by DAPC, and quality measures of the assignments to clusters.

## Competing interests

The authors declare that they have no competing interest.
## Authors' contributions

AS conducted the data collection, carried out the analyses and the interpretation of the results and drafted the manuscript. LC conducted the data collection and contributed to the interpretation of the results. MArb contributed to the design and coordination of the study, to data analysis and interpretation of the results. MArb conducted the data collection. LB and AB participated in the conception, design and coordination of the study, and contributed to the interpretation of the results. PR participated in the conception, design and coordination of the study, contributed to data analysis and interpretation of the results and drafted the manuscript. MHM participated in the conception, design and coordination of the study, carried out the data analysis and interpretation of the results and drafted the manuscript. All authors read and approved the final manuscript.

## Authors' information

LC and AS are MSc at the Institut Agronomique et Vétérinaire Hassan II (Rabat, Morocco); this study forms part of their PhD project. MArb, LB and AB are PhD, teachers and researchers at the Institut Agronomique et Vétérinaire Hassan II (Rabat, Morocco). MArd is lab technician and PR and MHM are PhD, researchers at the Institut National de la Recherche Agronomique (UMR Amélioration Génétique et Adaptation des Plantes méditerranéennes & tropicales, Montpellier, France).

## Acknowledgements

We thank Veronique Viader for flow cytometric analysis and Hélène Fréville and two anonymous referees for a critical review of the manuscript. This work was supported by Agropolis Resource Centre for Crop Conservation, Adaptation and Diversity, a flagship project of Agropolis Foundation, and by the Programme de Recherche Agronomique pour le Développement (project from Campus France).

## Author details

¹Département de Production, Protection et Biotechnologies Végétales, Institut Agronomique et Vétérinaire Hassan II, B.P. 6202, Rabat-Instituts, Rabat, Morocco. ²INRA, UMR 1334, Amélioration Génétique et Adaptation des Plantes méditerranéennes et tropicales (AGAP), 2 place Pierre Viala, F-34060 Montpellier Cedex 1, France.

Received: 10 July 2014 Accepted: 11 December 2014

Published online: 21 December 2014

## References

1. Evenson RE, Gollin D: Assessing the impact of the Green revolution. *Science* 2003, **300**:758-762.

2. Jarvis DI, Padoch C, Cooper HD: Biodiversity, agriculture, and ecosystem services. In *Managing biodiversity in agricultural ecosystems*. Edited by Jarvis DI, Padoch C, Cooper D. New York, USA: Columbia University Press; 2007:1-10.

3. Kang SK, Post WM, Nichols JA, Wang D, West TO, Bandaru V, Izaurralde RC: Marginal lands: concept, assessment and management. *J Agric Sci* 2013, **5**(5):129-139.

4. Altieri MA: Linking ecologists and traditional farmers in the search for sustainable agriculture. *Front Ecol Environ* 2004, **2**(1):35-42.

5. Lin BB: Resilience in agriculture through crop diversification: adaptive management for environmental change. *Bioscience* 2011, **61**(3):183-193.

6. Newton AC, Akar T, Baresel JP, Bebeli PJ, Bettencourt E, Bladenopoulos KV, Czembor JH, Fasoula DA, Katsiotis A, Koutis K, Koutsika-Sotiriou M, Kovacs G, Larsson H, Pinheiro Carvalho de MAA, Rubiales D, Russell J, Dos Santos TMM, Vaz Patto MC: Cereal landraces for sustainable agriculture a review. *Agron Sustain Dev* 2010, **30**:237-269.

7. Bradshaw JE, Ramsay G: Utilisation of the commonwealth potato collection in potato breeding. *Euphytica* 2005, **146**:9-19.

8. Reynolds M, Dreccer F, Trethowan R: Drought-adaptive traits derived from wheat wild relatives and landraces. *J Exp Bot* 2007, **58**(2):177-186.

9. Sharma S, Upadhyaya HD, Varshney RK, Gowda CLL: Pre-breeding for diversification of primary gene pool and genetic enhancement of grain legumes. *Front Plant Sci* 2013, **4**:309.

10. Abay F, Bjornstad A: Specific adaptation of barley varieties in different locations in Ethiopia. *Euphytica* 2009, **167**:181-195.

11. Jarvis DI, Brown AHD, Cuong PH, Collado-Panduro L, Latournerie-Moreno L, Gyawali S, Tanto T, Sawadogo M, Mar I, Sadiki M, Hue NT, Arias-Reyes L,

Balma D, Bajracharya J, Castillo F, Rijal D, Belqadi L, Rana R, Saidi S, Ouedraogo JT, Zangre R, Rhrib K, Chavez JL, Schoen D, Sthapit B, De Santis P, Fadda C, Hodgkin T: A global perspective of the richness and evenness of traditional crop-variety diversity maintained by farming community. *Proc Natl Acad Sci U S A* 2008, **105**(14):5326-5331.

12. Jarvis DI, Hodgkin T, Sthapit BR, Fadda C, Lopez-Noriega I: An heuristic framework for identifying multiple ways of supporting the conservation and use of traditional crop varieties within the agricultural production system. *Crit Rev Plant Sci* 2011, **30**(1-2):125-176.

13. Cai X, Fan J, Jiang Z, Basso B, Sala F, Spada A, Grassi F, Lu B-R: The puzzle of Italian rice origin and evolution: determining genetic divergence and affinity of rice germplasm from Italy and Asia. *PLoS One* 2013, **8**(11):e80351.

14. Odong TL, Jansen J, van Eeuwijk F, van Hintum TJL: Quality of core collections for effective utilisation of genetic resources review, discussion and interpretation. *Theor Appl Genet* 2013, **126**:289-305.

15. Dyer GA, Taylor JE: A crop population perspective on maize seed system in Mexico. *Proc Natl Acad Sci U S A* 2008, **105**(2):470-475.

16. Leclerc C, Coppens d'Eeckenbrugge G: Social organization of crop genetic diversity. The G × E × S interaction model. *Diversity* 2012, **4**:1-32.

17. Delétre M, McKey DB, Hodkinson TR: Marriage exchanges, seed exchanges, and the dynamics of manioc diversity. *Proc Natl Acad Sci U S A* 2011, **108**(45):18249-18254.

18. Deu M, Sagnard F, Chantereau J, Catalayud C, Hérault D, Mariac C, Pham J-L, Vigouroux Y, Kapran I, Traore PS, Mamadou A, Gerard B, Ndjunga J, Bezançon G: Niger-wide assessment of in situ sorghum genetic diversity with microsatellite markers. *Theor Appl Genet* 2008, **116**:903-913.

19. Pusadee T, Jamjod S, Chiang Y-C, Rerkasem B, Schaal BA: Genetic structure and isolation by distance in a landrace of Thai rice. *Proc Natl Acad Sci U S A* 2009, **106**(33):13880-13885.

20. Mercer KL, Martinez-Vasquez A, Perales HR: Asymmetrical local adaptation of maize landraces along an altitudinal gradient. *Evol Appl* 2008, **1**:489-500.

21. Tiranti B, Negri V: Selective microenvironmental effects play a role in shaping genetic diversity and structure in a *Phaseolus vulgaris* L. landrace: implications for on-farm conservation. *Mol Ecol* 2007, **16**:4942-4955.

22. Mercer KL, Perales HR: Evolutionary response of landraces to climate change in centers of crop diversity. *Evol Appl* 2010, **3**:480-493.

23. Van Heerwaarden J, Hellin J, Visser RF, Van Eeuwijk FA: Estimating maize genetic erosion in modernized smallholder agriculture. *Theor Appl Genet* 2009, **119**:875-888.

24. Samberg LH, Fishman L, Allendorf FW: Population genetic structure in a social landscape: barley in a traditional Ethiopian agricultural system. *Evol Appl* 2013, **6**:1133-1145.

25. Bajracharya J, Rana RB, Gauchan D, Sthapit BR, Jarvis DI, Witcombe JR: Rice landrace diversity in Nepal. Socio-economic and ecological factors determining rice landrace diversity in three agro-ecozones of Nepal based on farm surveys. *Genet Resour Crop Evol* 2010, **57**:1013-1022.

26. Jensen HR, Belqadi L, De Santis P, Sadiki M, Jarvis DI, Schoen DJ: A case study of seed exchange networks and gene flow for barley (*Hordeum vulgare* subsp. *vulgare*) in Morocco. *Genet Resour Crop Evol* 2013, **60**:1119-1138.

27. Samberg LH, Shennan C, Zavaleta E: Farmer seed exchange and crop diversity in a changing agricultural landscape in the southern highlands of Ethiopia. *Hum Ecol* 2013, **41**:477-485.

28. Sadiki M, Jarvis DI, Rijal D, Bajracharya J, Hue NN, Camachi TC, Burgos-May LA, Sawadogo M, Balma D, Lope D, Arias L, Mar I, Karamura D, Williams D, Chavez-Servia JL, Sthapit B, Rao VR: Variety names: an entry point to crop genetic diversity and distribution in agroecosystems? In *Managing biodiversity in agricultural ecosystems*. Edited by Jarvis DI, Padoch C, Cooper D. New York, USA: Columbia University Press; 2007:34-76.

29. Tsegaye S: Estimation of outcrossing rate in landraces of tetraploid wheat (*Triticum turgidum* L.). *Plant Breed* 1996, **115**(3):195-197.

30. Nesbitt M, Samuel D, Heller J: From staple crop to extinction? The archeology and history of the hulled wheat. In *First International Workshop on Hulled Wheats*. Edited by Padulosi S, Hammer K. Castelvecchio Pascoli Tuscany, Italy: IPGRI; 1996:41-100.

31. Peleg Z, Fahima T, Korol AB, Abbo S, Saranga Y: Genetic analysis of wheat domestication and evolution under domestication. *J Exp Bot* 2011, **62**(14):5051-5061.

32. Simonetti MC, Bellomo MP, Laghetti G, Perrino P, Simeone R, Blanco A: Quantitative trait loci influencing free-threshing habit in tetraploid wheats. *Genet Resour Crop Evol* 1999, **46**:267-271.
33. Simons KJ, Fellers JP, Trick HN, Zhang Z, Tai Y-S, Gill BS, Faris JD: Molecular characterization of the major wheat domestication Gene Q. *Genetics* 2006, **172**:547–555.

34. Peng JH, Sun D, Nevo E: Domestication evolution, genetics and genomics in wheat. *Mol Breed* 2011, **28**:281–301.

35. Tzarfati R, Saranga Y, Barak V, Gopher A, Korol A, Abbo S: Threshing efficiency as an incentive for rapid domestication of emmer wheat. *Ann Bot* 2013, **112**:829–837.

36. Matsuoka Y: Evolution of polyploid *Triticum* wheats under cultivation: The role of domestication, natural hybridization and allopolyploid speciation in their diversification. *Plant Cell Physiol* 2011, **52**(5):750–764.

37. Araus JL, Ferrio JP, Buxo R, Voltas J: The historical perspective of dryland agriculture: lessons learned from 10 000 years of wheat cultivation. *J Exp Bot* 2007, **58**(2):131–145.

38. Haudry A, Cenci A, Ravel C, Bataillon T, Brunel D, Poncet C, Hochu I, Poirier S, Santoni S, Glémin S, David J: Grinding up wheat: a massive loss of nucleotide diversity since domestication. *Mol Biol Evol* 2007, **24**(7):1506–1517.

39. Thuillet A-C, Bataillon T, Poirier S, Santoni S, David JL: Estimation of long-term effective population size through the history of durum wheat using microsatellite data. *Genetics* 2005, **169**:1589–1599.

40. Laido G, Mangini G, Taranto F, Gadaleta A, Blanco A, Cattivelli L, Marone D, Mastrangelo AM, Papa R, De Vita P: Genetic diversity and population structure of tetraploid wheats (*Triticum turgidum* L.) estimated by SSR, DArT and pedigree data. *PLOS One* 2013, **8**(6):e67280.

41. Oliveira HR, Campana MG, Jones H, Hunt HV, Leigh F, Redhouse DI, Lister DL, Jones MK: Tetraploid wheat landraces in the Mediterranean basin: taxonomy, evolution and genetic diversity. *PLoS One* 2012, **7**(5):e37063.

42. Hedden P: The genes of the green revolution. *Trends Genet* 2002, **19**(1):5–9.

43. Szabo AT, Hammer K: Notes on the taxonomy of farro: *Triticum monococcum*, *T. dicoccon* and *T. spelta*. In *Hulled wheats, promoting the conservation and use of underutilized and neglected crops*. Edited by Padulosi S, Hammer K, Heller J. Rome: IPGRi; 1996:2–40.

44. Miège M: Sur la présence au Maroc de *Triticum dicoccum* Sch. *Bull Soc Sci Nat Maroc* 1925, **5**(3):98–109.

45. Zaharieva M, Ayana NG, Al Hakimi A, Misra SC, Monneveux P: Cultivated emmer wheat (*Triticum dicoccon* Schrank), an old crop with promising future: a review. *Genet Resour Crop Evol* 2010, **57**:937–962.

46. Kehel Z, Garcia-Ferrer A, Nachit MM: Using bayesian and eigen approaches to study spatial genetic structure of Moroccan and Syrian durum wheat landraces. *Am J Mol Biol* 2013, **3**:17–31.

47. Teklu Y, Hammer K, Röder M: Simple sequence repeats marker polymorphism in emmer wheat (*Triticum dicoccon* Schrank): analysis of genetic diversity and differentiation. *Genet Resour Crop Evol* 2007, **54**:543–554.

48. Chentoufi L, Sahri A, Arbaoui M, Belqadi L, Birouk A, Roumet P, Muller M-H: Anchoring durum wheat diversity in the reality of traditional agricultural systems: varieties, seed management, and farmers’ perception in two Moroccan regions. *J Ethnobiol Ethnomed* 2014, **10**:

49. Jombart T, Devillard S, Balloux F: Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. *BMC Genet* 2010, **11**:94.

50. Mokuwa A, Nuijten E, Okry F, Teeken B, Maat H, Richards P, Struik PC: Processes underpinning development and maintenance of diversity in rice in West Africa: evidence from combining morphological and molecular markers. *PLoS One* 2014, **9**(1):e85953.

51. Abebe TD, Léon J: Spatial and temporal genetic analyses of Ethiopian barley (*Hordeum vulgare* L.) landraces reveal the absence of a distinct population structure. *Genet Resour Crop Evol* 2013, **60**:1547–1558.

52. Chakanda R, Van Treuren R, Visser B, van den Berg R: Analysis of genetic diversity in farmers’ rice varieties in Sierra Leone using morphological and AFLP markers. *Genet Resour Crop Evol* 2013, **60**:1237–1250.

53. Bonneuil C, Goffaux R, Bonnin I, Montalent P, Hamon C, Balfourier F, Goldringer I: A new integrative indicator to assess crop genetic diversity. *Ecol Indic* 2012, **23**:280–289.

54. Pautasso M, Aistara G, Barnaud A, Caillon S, Clouvel P, Coomes OT, Delêtre M, Demeulenaere E, De Santis P, Döring T, Eloy L, Emperaire L, Grarine E, Goldringer I, Jarvis DI, Joly HI, Leclerc C, Louafi S, Martin P, Massol F, McGuire S, McKey D, Padoch C, Soler C, Thomas M, Tramontini S:

Seed exchange networks for agrobiodiversity conservation. A review. *Agron Sustain Dev* 2013, **33**:151–175.

55. Barnaud A, Trigueros G, McKey D, Joly H: High outcrossing rates in fields with mixed sorghum landraces: how are landraces maintained? *Hed征 2008*, **101**:445–452.

56. Fraser JA, Alves-Pereira A, Junqueira AB, Peroni N, Clement CR: Convergent adaptations: bitter manioc cultivation systems in fertile anthropogenic dark earths and floodplain soils in central Amazonia. *PLoS One* 2012, **7**(8):e43636.

57. Labeyrie V, Deu M, Barnaud A, Catalayud C, Buiron M, Wambugu P, Manel S, Glaszmann JC, Leclerc C: Influence of Ethnolinguistic diversity on the sorghum genetic patterns in subsistence farming systems in eastern Kenya. *PLoS One* 2014, **9**(3):e92178.

58. Bellucci E, Bitocchi E, Rau D, Nanni L, Ferradini N, Giardini A, Rodriguez M, Attene G, Papa R: Population structure of barley landrace populations and gene-flow with modern varieties. *PLoS One* 2013, **8**(12):e83891.

59. Bitocchi E, Bellucci E, Giardini A, Rau D, Rodriguez M, Biagetti E, Santilocchi R, Spagnoletti Zeuli PL, Gioia T, Logozzo G, Attene G, Nanni L, Papa R: Molecular analysis of the parallel domestication of the common bean (*Phaseolus vulgaris*) in Mesoamerica and the Andes. *New Phytol* 2013, **197**:300–313.

60. Steele KA, Gyawali S, Joshi KD, Shrestha P, Sthapit BR, Witcombe JR: Has the introduction of modern rice varieties changed rice genetic diversity in a high-altitude region of Nepal? *Field Crop Res* 2009, **113**:24–30.

61. Elias M, Lenoir H, McKey D: Propagule quantity and quality in traditional Makushi farming of cassava (*Manihot esculenta*): a case study for understanding domestication and evolution of vegetatively propagated crops. *Genet Resour Crop Evol* 2007, **54**:99–115.

62. Dolezel J, Binarova P, Lucretti S: Analysis of nuclear DNA content in plant cells by flow cytometry. *Biol Plant* 1989, **31**:113–120.

63. Thuillet A-C, Roumet P, Bordat A, Tollon C, Poux G, Santoni S, David JL: Dépérdition ancienne et restauration moderne de la diversité génétique du blé. *Actes du BRG* 2008, **7**:405–422.

64. El Mousadik A, Petit RJ: High level of genetic differentiation for allelic richness among populations of the argan tree [*Argania spinosa* (L.) Skeels] endemic to Morocco. *Theor Appl Genet* 1996, **92**:832–839.

65. Nei M: *Molecular Evolutionary Genetics*. New York: Columbia Univeristy Press; 1987.

66. Goudet J: FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). Available from http://www2.unil.ch/popgen/software/fstat.htm.

67. Jombart T: *adegenet*: a R package for the multivariate analysis of genetic markers. *Bioinformatics* 2008, **24**(11):1403–1405.

68. Excoffier L, Laval LG, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analysis. *Evol Bioinformatics Online* 2005, **1**:47–50.

69. Weir BS: Genetic data analysis II. Sunderland: Sinauer Associates Inc.; 1996.

70. Belkir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F: GENETIX 4.02, logiciel sous Windows™ pour la génétique des populations. In: Montpellier, France: Laboratoire Génome, Populations, Interactions, CNRS UMR 5000, Université de Montpellier II; 2001.

71. Perrier X, Flori A, Bonnot F: Data analysis methods. In *Genetic diversity of cultivated tropical plants*. Edited by Hamon P, Seguin M, Perrier X, Glaszmann JC. Montpellier, France: Enfield, Science Publishers; 2003:43–76.

Submit your next manuscript to BioMed Central and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central