# Molecular Genotyping (SSR) and Agronomic Phenotyping for Utilization of Durum Wheat (Triticum durum Desf.) Ex Situ Collection from Southern Italy: A Combined Approach Including Pedigreed Varieties

Stefania Marzario, Giuseppina Logozzo, Jacques David, Pierluigi Spagnoletti Zeuli, Tania Gioia

► To cite this version:

Stefania Marzario, Giuseppina Logozzo, Jacques David, Pierluigi Spagnoletti Zeuli, Tania Gioia. *Molecular Genotyping (SSR) and Agronomic Phenotyping for Utilization of Durum Wheat (Triticum durum Desf.) Ex Situ Collection from Southern Italy: A Combined Approach Including Pedigreed Varieties*. Genes, 2018, 9 (10), 10.3390/genes9100465. [DOI: 10.26268/hal-02626831](https://wwwHAL.inrae.fr/hal-02626831)

HAL Id: hal-02626831

https://hal.inrae.fr/hal-02626831

Submitted on 26 May 2020

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>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Molecular Genotyping (SSR) and Agronomic Phenotyping for Utilization of Durum Wheat (*Triticum durum* Desf.) Ex Situ Collection from Southern Italy: A Combined Approach Including Pedigreed Varieties

Stefania Marzario¹, Giuseppina Logozzo ¹,*, Jacques L. David ², Pierluigi Spagnoletti Zeuli ¹
and Tania Gioia ¹ ID

¹ Scuola di Scienze Agrarie, Forestali, Alimentari ed Ambientali, Università degli Studi della Basilicata, viale dell'Ateneo Lucano 10, 85100 Potenza, Italy; stefy.85s@libero.it (S.M.); pierluigi.spagnoletti@unibas.it (P.S.Z.); tania.gioia@unibas.it (T.G.)

² Montpellier SupAgro, UMR AGAP Amélioration Génétique et Adaptation des Plantes méditerranéennes et Tropicales, 34060 Montpellier, France; jacques.david@supagro.fr

* Correspondence: giuseppina.logozzo@unibas.it; Tel.: +39-0971205533

Received: 16 August 2018; Accepted: 18 September 2018; Published: 20 September 2018

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

check for
updates

**Abstract:** In South Italy durum wheat (*Triticum durum* Desf.) has a long-time tradition of growing and breeding. Accessions collected and now preserved ex situ are a valuable genetic resource, but their effective use in agriculture and breeding programs remains very low. In this study, a small number (44) of simple sequence repeats (SSR) molecular markers were used to detect pattern of diversity for 136 accessions collected in South Italy over time, to identify the genepool of origin, and establish similarities with 28 Italian varieties with known pedigree grown in Italy over the same time-period. Phenotyping was conducted for 12 morphophysiological characters of agronomic interest. Based on discriminant analysis of principal components (DAPC) and STRUCTURE analysis six groups were identified, the assignment of varieties reflected the genetic basis and breeding strategies involved in their development. Some “old” varieties grown today are the result of evolution through natural hybridization and conservative pure line selection. A small number of molecular markers and little phenotyping coupled with powerful statistical analysis and comparison to pedigreed varieties can provide enough information on the genetic structure of durum wheat germplasm for a quick screening of the germplasm collection able to identify accessions for breeding or introduction in low input agriculture.

**Keywords:** *Triticum durum* (Desf.); germplasm; simple sequence repeats; morphophysiological traits; genetic diversity

## 1. Introduction

As part of a worldwide effort large collections of plant genetic resources have been established to meet future needs for agriculture and breeding and maintained ex situ in genebanks [1] but, due to the lack of genetic information, their effective use has been fairly small [2].

Durum wheat (*Triticum turgidum* L. ssp. *durum* Desf., 2n = 4x = 28; AABB genomes) is one of the most widely cultivated cereal species worldwide. It represents a major renewable resource for food, feed, and industrial raw materials [3] and is traditionally grown around the Mediterranean, mainly in Italy, Spain, France, Greece, West Asian, and North African countries [4]. The recent history of durum wheat has been characterized by the replacement of landraces and the improvement of modern
cultivars by inbred varieties in which dwarfing genes were introduced (second part of the 20th century) and adapted to favorable lands where high input agriculture is practiced [5,6]. So, the accessions collected might include germplasm from different genepools: cultivated varieties (cultivars) in current use; newly developed varieties; obsolete cultivars, primitive cultivars (landraces); but also, special genetic stocks (including elite and current breeders' line and mutants).

Studies on genetic diversity within ex situ germplasm collections are needed since their effective utilization for breeding relies on the ability to identify accessions with a given characteristic or possibly with the desired gene. The passport data, when available, are useful information since they relate the accession to where and when it was obtained/collected, thus providing hints on the more likely genetic characteristics or sometimes the gene pool of origin. Phenotypic characterization data of agronomic interest both for qualitative and quantitative traits (e.g., spike shape, spike density, seed color, days to heading, plant height, and spike length) are very useful but for the characterization of large collections a huge amount of money and human resources are needed and must deal with little polymorphism, low heritability, and genotype by environment interaction [1]. In addition, morphological traits selected for scoring might end up assessing diversity for a small, not representative amount of the whole genome.

The advances in molecular genotyping provide tools easier to handle the study variation in germplasm collections [7]. DNA molecular markers have been used to assess the amount of variation and describe the genetic relationships among accessions [8]. For their high polymorphism, codominance and locus specificity, simple sequence repeats or microsatellite (SSRs) markers have proved to be highly efficient molecular tools for the characterization of durum wheat germplasm collections [9–14].

To date most studies on Italian durum germplasm have analyzed collections including old and new elite varieties for morphophysiological and qualitative traits [15,16], and the use of molecular markers was focused on temporal trends of diversity [8,17,18], relatedness among genotypes [10], genetic structure [19], also in comparison to *Triticum turgidum* L. subspecies [13].

In a recent study [18] a panel of 370 durum wheat genotype including 35 Italian genotypes (29 elite varieties and six landraces) where genotyped using 500 simple sequence repeats (SNPs) markers. All varieties were obtained from crosses and not by pure line selection, and, except for one (Capeiti 8), were all "modern" released after 1974 when Creso, the first Italian successful dwarf variety, was released.

In Italy durum wheat has a long-time tradition of growing and breeding and accessions collected in Southern Italy, now preserved ex situ, are a valuable genetic resource. A large number of Italian and some north African landraces were grown in Italy since the end 1800 and after 1920 [20] were gradually substituted by the pure line varieties selected from landraces from Italy, but also from Syria-Palestinia and north Africa (Lybia) [21]. Almost all the genetic material was from *Mediterraneum typicum* and *Syriacum typicum* groups [22,23] where Syriacum with erect leaves and short head, did show a shorter plant, earlier heading, higher tillering, than Mediterraneum.

New varieties were then obtained from crosses between lines from Mediterraneum (Senatore Cappelli) and Syriacum (Aziziah, Eiti, Sinai, and Tripolino) groups (1950–1965), from the introgression from tetraploid and hexaploid species and mutagenesis (1965–1973), the introduction of Norin10 dwarfing genes (1970–1980) and of new germplasm mostly from CIMMYT (Centro Internacional de Mejoramiento de maiz y trigo, Messico) and from recombination of the whole available tetraploid gene pool (1980–1990), followed [24–28].

Luckily, as the result of extensive germplasm collecting expedition programs that did start in 1948, durum wheat germplasm from inland South Italy was collected; also in the most remote areas and preserved ex situ collections [21]. Today these landraces/old varieties, along with their breeding potential, are of interest for cultivation in marginal areas and in organic or low-input agricultural systems since the typical local products pay a premium price and ex situ collection can provide this germplasm. The detailed passport data that are available (e.g., accession name, collection site, geographical origin, and collecting year) only in few cases include a "variety" name. Almost no information is available to identify the gene pool of origin and length of time for evolution at the
collection site or about their genetic characteristics, which are important because of the relevant role of durum wheat breeding in Italy.

In this study a strategy combining molecular genotyping, phenotyping, and the use of pedigreed varieties is proposed to characterize a durum wheat ex situ germplasm collection: (i) a small number of SSR molecular markers were used to detect pattern of diversity for 136 accessions collected in South Italy over time; (ii) to identify the genepool of origin of these ex situ accessions the same molecular markers were used to establish similarities with 28 Italian varieties with known pedigree; and (iii) phenotyping was then conducted for 12 morphophysiological characters of agronomic interest and the structure of diversity compared to identify the most diverse accessions for breeding purposes or for use in low input agriculture.

The results show that a small number of molecular markers and little phenotyping coupled with powerful statistical analysis and comparison to pedigreed varieties can provide enough information on the genetic structure of durum wheat germplasm for a quick screening of the ex situ germplasm collection, in order to identify accessions for breeding or introduction in low input agriculture.

## 2. Methods

### 2.1. Plant Material

The analyzed durum wheat germplasm included 136 accessions collected in Southern Italy from 1947 to 2003, now preserved ex situ in genebanks, and of 28 elite varieties cultivated in the past decades in Italy and representative of the Italian breeding programs.

Seed samples for all accessions were obtained by the Institute of Plant Genetics and Crop Plant Research (IPK), Germany; Centro per la Salvaguardia delle Risorse Genetiche Vegetali "Pierino Iannelli", Università degli Studi della Basilicata (Unibas), Italy; and the Istituto di Bioscienze e Biorisorse del Consiglio Nazionale delle Ricerche (IBBR-CNR), Italy. Passport information (e.g., accession name, collsite, geographical origin, and collecting year) was available for all the accessions. In only a few cases the accession was also identified with a "variety" name. Luckily the accessions of eight very old "varieties" have been collected in 1950 from an institution well-known for durum wheat breeding (Stazione Sperimentale di Granicoltura, Catania). But in most cases, none of the information available was useful to distinguish "landraces" from varieties or to identify the gene pool of origin and length of time for evolution at the collection site. The seeds of 28 varieties were obtained from the Centro di Ricerca per la Cerealcoltura-Foggia (CREA-CER), Italy. Year of release, registered pedigree, and breeder for each elite variety were recorded. A complete list of the accessions analyzed is available in Supplementary Materials file 1: Tables S1 and S2.

For the purpose of the study, accessions were subdivided into groups according to the year in which the collecting expeditions took place: Group1 (number of accessions = 37) comprised of accessions collected in the period 1947–1950; Group2 (number of accessions = 77) included accessions collected from 1973 to 1982; and Group3 (number of accessions = 22) consisted of accessions collected during the period 1983 to 2003.

According to the 28 varieties' release year, two temporal groups were identified: (i) "old and intermediate" group (1915–1975) including "old" genotypes selected from indigenous and exotic landraces, and "intermediate" genotypes selected from crosses or mutagenesis involving old materials (n = 9); (ii) "modern" including genotypes released after the 1980s, mostly selected from crosses between CYMMIT (International Maize and Wheat Improvement Center) breeding lines and materials belonging to the former set of varieties (n = 19).

### 2.2. Molecular Characterization

In order to estimate molecular diversity, the germplasm collection was characterized using 44 SSRs markers. All the nuclear SSRs markers were selected on the basis of their known genetic locations to give a uniform coverage for all seven linkage groups as well as their temperature of annealing (Tan).
and degree of polymorphism [29–32]. More information about the SSR markers used, including the designation, chromosome location, primer pair sequences, the repeat motif, and the allele size range of the amplified loci are available in Supplementary Materials file 1: Table S3.

## 2.3. Genomic DNA Extraction and Genotyping

Genomic DNA was extracted from leaf tissues at the tillering stage of greenhouse-grown plants using the CTAB (Cetyltrimethylammonium bromide) method [33] with minor modifications. The purity and concentration of DNA was determined spectrophotometrically at 260 and 280 nm by using Nano Drop® Spectrophotometer ND 1000 (NanoDrop Technologies, Wilmington, DE, USA). The DNA samples were diluted to a final concentration of 20 ng/µL with TE (Tris-EDTA) for the following SSRs analysis.

Multiplex PCR (two to five SSRs primer pairs with the same annealing temperature in one PCR reaction) and multiple sample loading (loading more than one PCR reaction in each well) were used. The reverse primer for each locus was labeled with one of the three fluorophores (6FAM, NED, or HEX) and PCR amplification was performed with a final volume of 20 µL containing 30 ng of genomic DNA, 0.2 mM of each dNTP, 2 nM MgCl₂, 1–4 pM of each primer, and 0.5 U of Taq DNA polymerase. PCR was carried out as follows: after 5 min at 94 °C, 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 a final extension step of 30 min at 72 °C. The PCR products were detected by capillary electrophoresis using an ABI PRISM 3130xl Genetic Analyzer (Applied Biosystems, Foster City, CA, USA) and analyzed using the GENEMAPPER V3.7 genotyping software (Applied Biosystems).

## 2.4. Morphological Characterization

Three plants for each of the 136 accessions from ex situ genebanks and 28 varieties were grown in the greenhouse at the Università degli Studi della Basilicata (Potenza, Italy) during the growing season of 2012 to 2013. Heading date (expressed as days after 10 April) was recorded when the culms showed emerging spikes and plant height (excluding spike) was measured at harvest time. Individual plants were harvested manually in July 2013 and ten quantitative morphological traits were scored according to descriptors for wheat defined by IBPGR (International Board for Plant Genetic Resources) [34]: spike length (cm), spikelets number/spike (n), fertile spikelets number/spike (n), kernel length (mm), kernel width (mm), kernels number/spike (n), kernels weight/spike (g), 1000 kernels weight (g), kernels number/spikelet (n), and kernels number/fertile spikelets (n). Three characters were calculated: (i) thousand kernels weight (g) as (kernels weight per spike/kernel number per spike) × 1000; (ii) kernels number per spikelet as kernel number per spike/spikelet number per spike; (iii) kernels number per fertile spikelet as kernel number per spike/fertile spikelet number per spike.

## 2.5. Statistical Analysis

Summary statistics for each SSRs locus, including the number of alleles detected and the gene diversity or unbiased expected heterozygosity ($H_e$ [35]) were calculated using the software Genetix V4.05 [36].

The level of variation in different germplasm types (landraces and elite varieties) and within temporal groups of accessions was evaluated in terms of number of alleles per locus ($N_a$), gene diversity, or unbiased expected heterozygosity ($H_e$ [34]) and allelic richness ($R_s$ [37]). All these indices were calculated using FSTAT 2.9.3.2 software [38]. For computing $R_s$, FSTAT implements the rarefaction method to trim the unequal accessions number to the same standardized sample size, a number equal to the smallest sample across the population. Overall differences in the expected heterozygosity and allelic richness between temporal groups of accessions were assessed for significance using FPTestR software [39]. Population structure within the germplasm collection was examined by first applying the Discriminant analysis of principal components (DAPC) [40], a new multivariate method designed to identify and describe clusters of genetically related individuals. This approach allows
extracting rich information from genetic data, providing assignment of individuals to groups, a
visual assessment of between-population differentiation, and contribution of individual alleles to
population structuring. Discriminant analysis of principal components does not assume unrelatedness,
and therefore, potentially closely related individuals can be included in the analysis. The method relies
on allele data transformation using principal component analysis (PCA) as a prior step to discriminant
analysis (DA). Discriminant analysis of principal components was performed using the adegenet
package [41] in R (R development Core Team, 2009). 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 or increased only subtly. We then
applied the *dapc* function to describe the relationship between the inferred groups. This function
constructs synthetic variables and discriminant functions (DFs) that maximize variation between
while minimizing variation within groups and computes coordinates along these functions for each
individual. In order to obtain reliable group membership probabilities and to avoid overfitting, we
retained only the four first principal components (PCs) from the preliminary data transformation step
(indicated to be the optimal number based on the *optim.a.score* function).

The assignment of each accession to the different clusters was examined taking into account the
germplasm type (accessions or varieties) and the temporal group to which each genotype belonged.

The model-based approach implemented in the software package STRUCTURE [42] was then
further applied to infer population structure. Initially, the number of subgroups (K) was set from
two to ten. Twenty independent simulations were performed for each K setting using the admixture
model, with each simulation set to a 5000 burn-in period and 50,000 Markov chain Monte Carlo
(MCMC) repetitions. To determine the optimal number of clusters [43], STRUCTURE HARVESTER [44],
was used for calculate the Delta K statistical test in combination with the likelihoods (posterior
probabilities) of each preset K. Results from simulations with the highest likelihood within each
number of different K simulations were chosen to assign accessions to populations. Accessions with a
population membership coefficient of less than 0.7 were identified as potential hybrids.

Pairwise FST metrics [45] were also calculated in Genetix V4.05 [36] to estimate the divergence
between the clusters identified by DAPC analysis. The value of FST vary from zero to one; when
FST = 0 the subpopulations are identical, while when FST = 1 they are completely differentiated in
relation to the fixation of different alleles in the subpopulations.

Differences of trait means were checked by a nested analysis of variance (ANOVA) based on
the following sources of variation, between clusters, between germplasm types (landraces and elite
varieties) within clusters, and between accessions among germplasm types within clusters.

Principal component analysis (PCA) was performed to determine the overall morphological traits
distinctiveness, and to investigate the relationships between the traits. Principal component analysis
was performed using the software JMP version 8.0 [46].

3. Results

3.1. *Overall Levels of Genetic Diversity*

Forty-four SSR markers broadly distributed between the A and B durum wheat genomes were
analyzed in the present study. All SSR markers showed polymorphism and a total of 242 alleles were
detected across the 164 durum wheat genotypes. The average N<sub>a</sub> per SSR was 5.50, ranging from
two alleles (*Xgwm357*, *Xgwm165-4A*, *Xgwm415*, *Xgwm408*, and *Xgwm169*) to 14 alleles (*Xgpw2302*)
(Supplementary Materials file 1: Table S3). Gene diversity or H<sub>e</sub> across the total accessions was 0.607,
and ranged from a low 0.105 (*Xgwm374*) to a high 0.832 (*Xgwm6*) (Supplementary Materials file 1:
Table S3).

When the whole germplasm collection was subdivided into two groups on the basis of the
improvement status of each accession (accessions vs. varieties), a higher mean N_a per locus was
observed in the accessions (N<sub>a</sub> = 5.39) compared to the varieties (N<sub>a</sub> = 3.86) (Table 1). FPTest [39] among accessions and varieties indicated significant differences in allelic counts (R<sub>s</sub> = 4.562 and R<sub>s</sub> = 3.830, *p* < 0.01).

In order to evaluate how genetic diversity changed over time different groups of accessions were identified and compared. When considering the accessions, a slight, but not significant, increase of diversity was observed from Group1 (H<sub>e</sub> = 0.579, R<sub>s</sub> = 3.715) to Group2 (H<sub>e</sub> = 0.603, R<sub>s</sub> = 3.882), followed by a decrease in variation in Group3 (H<sub>e</sub> = 0.531, R<sub>s</sub> = 3.303). Nevertheless, when considering the varieties, a loss of genetic diversity was observed from the beginning of the 20th century to nowadays germplasm. Indeed, the highest H<sub>e</sub> and R<sub>s</sub> were observed in "old and intermediate" genotypes (H<sub>e</sub> = 0.604, R<sub>s</sub> = 3.434) as compared to "modern" varieties (H<sub>e</sub> = 0.534, R<sub>s</sub> = 3.011).

<table><caption>**Table 1.** Genetic diversity estimates for 136 ex situ accessions per collection year and 28 durum wheat varieties per year of release.</caption><thead><tr><th>Germplasm Group</th><th>n</th><th>N<sub>a</sub></th><th>H<sub>e</sub></th><th>R<sub>s</sub></th></tr></thead><tbody><tr><td>Accessions</td><td>136</td><td>5.39</td><td>0.601</td><td>4.562</td></tr><tr><td>Group1 (1947–1950)</td><td>37</td><td>4.34</td><td>0.579</td><td>3.715</td></tr><tr><td>Group2 (1973–1982)</td><td>77</td><td>5.00</td><td>0.603</td><td>3.882</td></tr><tr><td>Group3 (1983–2003)</td><td>22</td><td>3.57</td><td>0.531c</td><td>3.303</td></tr><tr><td>Varieties</td><td>28</td><td>3.86</td><td>0.576</td><td>3.830</td></tr><tr><td>Old and intermediate varieties (1915–1973)</td><td>9</td><td>3.45</td><td>0.604</td><td>3.434</td></tr><tr><td>Modern varieties (1974–2007)</td><td>19</td><td>3.09</td><td>0.534</td><td>3.011</td></tr></tbody></table>

*n* = number of samples; N<sub>a</sub> = mean number of alleles per locus; H<sub>e</sub> = expected heterozygosity; R<sub>s</sub> = allelic richness computed for a minimum sample size of nine diploid individuals per group.

## 3.2. Genetic Structure of the Wheat Collection

The population structure of the germplasm collection, including both accessions and pedigreed varieties, was first assessed using a DAPC analysis, a method that identifies and describes clusters of genetically related individuals within a dataset [40]. The BIC analysis used to identify the optimal number of clusters clearly identified six genetic clusters (K = 6, Supplementary Materials file 1: Figure S1).

To describe the relationships among the six clusters, a scatterplot of the first two principal components of the DAPC, which accounted for 28% of the total variance, is presented in Figure 1. The distribution of the ex situ accessions and the varieties among the six clusters is reported in Table 2. Germplasm ex situ accessions were not separated from varieties, but varieties representing a homogeneous genepool (related pedigrees) were included in the same cluster. Three clusters (C3, C4, and C6) included germplasm accessions and only one variety in each one; C5 and C1 included mostly germplasm accessions, with only two (C5) and four (C1) varieties. A large number of the analyzed varieties (n = 19) did cluster in C2 with 10 germplasm accessions (Table 2). The distribution of diversity over collection time shows how the pattern of substitution of old varieties by newly released varieties did occur. Additionally, accessions collected from (1947–1950) were never included in C1 and C2, accessions collected in (1973–1982) were present in each of the DAPC clusters and genotypes collected from 1983 to 2003 (third temporal group) were assigned to all but one cluster (C5) (Table 2). Clustering revealed that the choice of collecting sites was very careful at avoiding the most recently introduced varieties. Out of 136 accessions only 10 did cluster in C2 that included most "modern" varieties and 23 in C1.

The assignment of the varieties largely reflected the genetic basis and breeding strategies applied for their development. Cappelli, selected from the North African population "Jean Retifah" and considered one of the most relevant ancestors of the modern durum wheat gene pool, did cluster in C3. The same cluster also included two accessions (TRI3873 and TRI3882) collected in 1950 in Calabria that, according to the passport data information, were also named, "Senatore Cappelli" and the accession named "Bidi" that, according to De Cillis (1927), was from a variety from Tunisia, was not different from the one that originated Cappelli. Russello, a variety selected from a Sicilian local population, was the only variety included in C4 where six out of the ten accessions were collected in Sicily in 1950, including
accession TRI3857 also named "Russello", TRI3851 named "Gioia", and TRI3855 named "Regina", which are old Sicilian varieties. It should be noted that, according to molecular data, accessions and varieties with the same "variety" name that cluster in the same group do not have the same genotypes thus cannot be considered "duplicates". Two old varieties, namely Aziziah, selected from the exotic landraces "Near East", and Timilia, selected from Sicilian landraces, did cluster in C5 with accession TRI14111, named "Triminia", collected in Campania in 1982. Capeiti 8, Appulo, Simeto, and Ciccio, that could be traced back to a common origin on the bases of pedigree information (Cappelli × Eiti 6), were assigned to C1 which also included accession TRI16687, named "Trinakria", a variety that was obtained using Capeiti 8 as parent. Except for the variety Isa, that did cluster to C6 with accessions TRI3518 named "Cotrone", TRI3519 named "Inglesa", TRI3520 named "Scorsonera", and TRI3853 named "Bianmaia", which are also old varieties, all the remaining 19 varieties were assigned to C2. These genotypes were the most recently released varieties and were related, according to their pedigree information, to Creso, and to the innovative semidwarf materials developed at CYMMIT. Accession TRI13079 named "Creso" and TRI16565 named "Capeiti" were also included in C2.

<table><caption>**Table 2.** Number of accessions per collection year and varieties per year of release in six clusters obtained by discriminant analysis of principal components (DAPC) for 136 ex situ accessions and 28 durum wheat varieties.</caption><thead><tr><th>Germplasm</th><th>Year</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>Total</th></tr></thead><tbody><tr><td>Accession</td><td>Collection</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Group1 (1947–1950)</td><td>1947–1950</td><td></td><td></td><td>16</td><td>7</td><td>2</td><td>12</td><td>37</td></tr><tr><td>Group2 (1973–1982)</td><td>1973–1982</td><td>17</td><td>5</td><td>19</td><td>1</td><td>17</td><td>18</td><td>77</td></tr><tr><td>Group3 (1983–2003)</td><td>1983–2007</td><td>6</td><td>5</td><td>6</td><td>2</td><td>-</td><td>3</td><td>22</td></tr><tr><td>Total</td><td></td><td>23</td><td>10</td><td>41</td><td>10</td><td>19</td><td>33</td><td>136</td></tr><tr><td>Varieties</td><td>Release</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Old and intermediate</td><td>1915–1973</td><td>2</td><td>3</td><td>1</td><td>1</td><td>2</td><td>-</td><td>9</td></tr><tr><td>Modern</td><td>1974–2007</td><td>2</td><td>16</td><td>-</td><td>-</td><td>-</td><td>1</td><td>19</td></tr><tr><td>Total</td><td></td><td>4</td><td>19</td><td>1</td><td>1</td><td>2</td><td>1</td><td>28</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Scatterplot of the first two principal components of the DAPC applied on 136 ex situ accessions and 28 durum wheat varieties. Minimum spanning tree connects the six groups. Varieties (*) and named accessions are reported. Numbers and colors identify the clusters. LD: loadings.</figcaption></figure>
To further describe the population structure of the whole germplasm collection, a Bayesian model-based clustering method implemented in the STRUCTURE program was performed. The number of subpopulations (K) was identified based on maximum likelihood and Delta K values [43]. The Delta K test suggested that our sample was made up of four main genetic groups (K = 4), with the next largest peak found at seven clusters (K = 7) (Supplementary Materials file 1 Figure S2). At K = 4, four clusters were clearly separated: the first cluster (K1; blue color in Figure 2) was composed of 17 accessions and three old-intermediate varieties; the second cluster (K2, pink color in Figure 2) included 16 varieties and only six accessions; the third cluster (K3; green color in Figure 2) included a total of 36 accessions and no varieties; and the fourth cluster (K4; red color in Figure 2) was composed of 58 accessions and only two old-intermediate.

The STRUCTURE bar graphic (Figure 2) also provides information on the level of admixture in the study sample. At K = 4, 138 genotypes out of 164 (84%) were assigned to one or another group with more than 70% posterior probability (Figure 2). The remaining 26 not assigned accessions (16%) were assumed to have a mixed ancestry.

Individual assignments provided by STRUCTURE differed only slightly from those by DAPC: one of the STRUCTURE clusters (K3) was split into three DAPC groups, while the other three clusters where consistent with the DAPC assignments (Figure 3).

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Population structure analysis for 136 accessions and 28 durum wheat varieties based on SSRs analysis. K-values of subpopulations are shown to right and accessions and varieties grouped according to the six clusters identified by DAPC are given below. Each individual is represented by a vertical line, and cluster assignments are indicated by color. Individuals are considered assigned to a cluster if their posterior probability in that cluster is at least 0.7. Varieties are indicated with a star.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Comparison of individual assignment between DAPC and STRUCTURE analyses for 136 ex situ accessions and 28 durum wheat varieties based on 44 SSRs. Individuals are considered assigned to a cluster if their posterior probability in that cluster is at least 0.7.</figcaption></figure>

### 3.3. Genetic Differentiation among DAPC Clusters

Genetic differentiation among the six DAPC clusters was investigated by computing pairwise F<sub>ST</sub> values. All pairwise F<sub>ST</sub> values were significant and ranged from 0.122 (C5–C6) and 0.502 (C3–C4), as shown in Table 3.

<table><caption>**Table 3.** Differentiation indices (F<sub>ST</sub>) between all pairwise combination of clusters identified by DAPC performed on 136 durum wheat ex situ accessions and 28 varieties.</caption><thead><tr><th>DAPC Cluster</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th></tr></thead><tbody><tr><td>C1</td><td></td><td>0.225 **</td><td>0.284 **</td><td>0.415 **</td><td>0.317 **</td><td>0.220 **</td></tr><tr><td>C2</td><td></td><td></td><td>0.256 **</td><td>0.276 **</td><td>0.218 **</td><td>0.145 **</td></tr><tr><td>C3</td><td></td><td></td><td></td><td>0.502 **</td><td>0.458 **</td><td>0.327 **</td></tr><tr><td>C4</td><td></td><td></td><td></td><td></td><td>0.251 **</td><td>0.158 **</td></tr><tr><td>C5</td><td></td><td></td><td></td><td></td><td></td><td>0.122 **</td></tr><tr><td>C6</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

To get some insight into the underlying causes of the differentiation among clusters, we inspected the associated allele loadings. DAPC is perfectly adapted to this task, as it finds linear combination of alleles (the discriminant functions) which best separate the clusters. Alleles with the largest contribution to this discrimination are the most different among groups. Therefore, the plot of allele contribution (Figure 4) could be useful for a graphical assessment of alleles of major interest.

The locus *Xgwm408*, located on the long arm of chromosome 5B, mostly contributed to the first principal component with allele 180 and allele 200. The first allele discriminated C4 and C5, which was not represented, from the other clusters, whereas the second allele was not recorded in accessions attributed to C1 and C3. Allele 170 of locus *Xgwm282*, located on the long arm of chromosome 7A, was detected in all DAPC clusters except for C1 and C3. The locus *Xgwm234*, on the short arm of chromosome 5B, contributed to the first principal component with allele 224 detected only in accessions belonging to C1, C2, C3, and C6. The same discrimination pattern was observed for allele 175 of the locus *Xgpw2302*, located on the short arm of chromosome 4A.

<footnote>(** significant for *p* < 0.01).</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Contribution of alleles to the first (a) and the second (b) principal components of DAPC applied on 136 accessions and 28 durum wheat varieties. The height of each bar is proportional to the contribution of the corresponding allele to the first and second principal components of the analysis, respectively. Only alleles whose contribution was above a threshold (grey horizontal line) are indicated for the sake of clarity.</figcaption></figure>
An additional locus mapped on the short arm of chromosome 4A, *Xgpw2239*, contributed to the first principal component with two alleles, namely allele 196 and allele 198. The first allele was recorded in all clusters except for C4, while allele 198 was detected only in C3, C4, C5, and C6. The allele 219 identified at locus *Xgwm285* on the short arm of chromosome 3B was useful to discriminate C1, C2, and C3 from the remaining clusters identified by DAPC analysis, in which this allele was not recorded. Finally, allele 96 at locus *Xgwm413* (short arm of chromosome 1B), allele 101 identified at locus *Xgpw7148* (long arm of chromosome 3B), and allele 224 for locus *Xgpw4004* (long arm of chromosome 5A), showed similar behavior in discriminating DAPC clusters. Indeed, these alleles were all detected in C1, C2, C3, and C6.

When considering the second principal component, four loci showed the largest contribution to the discrimination among clusters. Regarding the locus *Xgwm193*, located on the short arm of chromosome 6B, two alleles showed a high contribution; allele 170 was recorded in all DAPC clusters except for C5, whereas allele 172 was identified in accessions assigned to all clusters except for C4. A similar level of discrimination was observed for an additional locus mapped on the short arm of chromosome 6B, namely *Xgwm518*. Allele 120 was observed in all clusters except for C4; allele 170 was useful to discriminate C5, whose accessions were not characterized by this allele, from the remaining clusters. The locus *Xgpw2103*, located on the short arm of chromosome 7A, contributed with allele 232 and allele 234. The first allele was not recorded in C3, whereas the second allele was not detected in C1, C4, and C5. Finally, allele 102 identified at locus *Xgwm413* (short arm of chromosome 1B), characterized genotypes attributed to C1, C2, C5, and C6.

### 3.4. Morphological Characterization

To identify the relationship between molecular SSR diversity and phenotypic characters of agronomic interest, a greenhouse phenotyping was conducted. Twelve morphophysiological characters were scored on three plants for each accession and variety.

The differences among the means of 12 characters were tested by a nested analysis of variance (ANOVA): between clusters, between germplasm type (accessions and varieties), and between accessions among germplasm type (Supplementary Materials file 1: Table S4). Differences among clusters were highly significant for all characters except for kernel width and kernels number/spike. For all remaining sources of variation means were always significantly different. For all the traits the summary statistics for each group (accessions and varieties) included in the analyzed germplasm are reported in Supplementary Materials file 1: Table S5. Ex situ accessions showed a higher mean heading date (38.87 days) than varieties (36.17 days) and were taller (mean plant height = 101.69 cm) than varieties (mean plant height = 72.09 cm). In accessions the spike was longer (19.53 cm) and had a higher number of spikelets (n = 21.73) than in varieties (18.06 cm and n = 17.40, respectively). Similarly, the number of fertile spikelets per spike was higher in accessions (n = 17.50) when compared to the varieties (n = 15.69). By contrast, kernel number per spikelet was lower in accessions (n = 1.68) than in varieties (n = 2.21). Additionally, higher values were recorded for both kernel weight/spike and 1000 kernels weight in accessions (1.95 g and 52.26 g, respectively) in comparison to the varieties (1.81 g and 47.42 g, respectively).

### 3.5. Morphological Differentiation among Discriminant Analysis of Principal Components Clusters

The six genetic structure clusters revealed by DAPC were further investigated to detect pattern of diversity for morphological traits of agronomic interest. A nested analysis of variance (ANOVA) was run. Means were significantly different (*p* < 0.001) between clusters for all the traits except for kernel width and kernels number/spike. Additionally, means were highly significant (*p* < 0.001) for all the traits between germplasm types within clusters and between accessions among germplasm types within cluster (Table 4).

The mean and the range of variation for the six clusters grouped according to DAPC are shown in Supplementary Materials file 1: Figure S3 along with the overall mean. The means and the ranges can
be used to identify phenotypically divergent sources for traits of interest in breeding programs. For all traits within the six clusters groups different ranges were expressed. For some traits (e.g., heading date and fertile spikelets number per spike) the ranges expressed within cluster 2 (made of 19 varieties) were narrower than in other clusters, suggesting a low degree of diversity. Sampling from different clusters is expected to increase allelic diversity and selection response in breeding programs.

<table><caption>**Table 4.** Analysis of variance (ANOVA) of 136 durum wheat ex situ accessions and 28 varieties for 12 morphophysiological plant characters comparing the six clusters identified by DAPC, the germplasm types (accessions and varieties) within clusters, and the accessions within type group and cluster.</caption><thead><tr><th rowspan="2">Traits</th><th colspan="3">Source of Variations</th><th rowspan="2">Error</th></tr><tr><th>Between Clusters</th><th>Between Germplasm Types (Clusters)</th><th>Between Accessions (Germplasm Type × Clusters)</th></tr></thead><tbody><tr><td>Heading date (from 10 April)</td><td>268.1 ***</td><td>186.2 ***</td><td>126.4 ***</td><td>8.7</td></tr><tr><td>Plant height (cm)</td><td>13,786.8 ***</td><td>394.1 ***</td><td>519.5 ***</td><td>95.1</td></tr><tr><td>Spike length (cm)</td><td>98.7 ***</td><td>28.9 ***</td><td>9.4 ***</td><td>2.2</td></tr><tr><td>Spikelets number/spike (n)</td><td>244.2 ***</td><td>47.6 ***</td><td>16.9 ***</td><td>1.8</td></tr><tr><td>Fertile spikelets number/spike (n)</td><td>81.9 ***</td><td>25.2 ***</td><td>20.3 ***</td><td>6.3</td></tr><tr><td>Kernel length (mm)</td><td>1.7 ***</td><td>0.6 ***</td><td>0.7 ***</td><td>0.1</td></tr><tr><td>Kernel width (mm)</td><td>0.06 ns</td><td>0.3 ***</td><td>0.4 ***</td><td>0.05</td></tr><tr><td>Kernels number/spike (n)</td><td>263.0ns</td><td>39.1 **</td><td>182.4 ***</td><td>65.4</td></tr><tr><td>Kernels weight/spike (g)</td><td>1.7 ***</td><td>1.1 ***</td><td>1.0 ***</td><td>0.2</td></tr><tr><td>1000 kernels weight (g)</td><td>370.7 ***</td><td>569.9 ***</td><td>323.9 ***</td><td>46.5</td></tr><tr><td>Kernels number/spikelet (n)</td><td>2.6 ***</td><td>0.7 ***</td><td>0.3 ***</td><td>0.1</td></tr><tr><td>Kernels number/fertile spikelet (n)</td><td>1.2 ***</td><td>0.5 ***</td><td>0.2 ***</td><td>0.1</td></tr></tbody></table>

Principal components analysis (PCA) analysis on all traits revealed four principal components (PC1–PC4) explaining 90% of total variance (Table 5). The first principal component (PC1) explained 42.7% of total variance and was mostly associated (higher loadings) with kernels weight/spike, kernels number/spike, 1000 kernels weight, kernel width, kernel length, and heading date, 26.8% of additional variation projected by PC2 came primarily from plant height and spikelets number/spike. PC3 and PC4 added, respectively, another 14.9% and 4.5% of the observed traits variability: these components accounted mainly for kernels number/spike, 1000 kernels weight, and kernel length and width (Table 5). A plot of genotype distribution for the first two principal components is presented in Figure 5.

<table><caption>**Table 5.** Correlations between four principal components variables and 12 original morphophysiological plant characters for six groups obtained from 136 durum wheat "ex situ" accession and 28 varieties.</caption><thead><tr><th>Principal Components</th><th>PC1</th><th>PC2</th><th>PC3</th><th>PC4</th></tr></thead><tbody><tr><td>Kernels weight/spike</td><td><strong>0.42</strong></td><td>-0.04</td><td>-0.01</td><td>0.11</td></tr><tr><td>Heading date</td><td><strong>-0.31</strong></td><td>0.27</td><td>0.21</td><td>0.37</td></tr><tr><td>Plant height</td><td>0.14</td><td><strong>0.46</strong></td><td>0.08</td><td>-0.14</td></tr><tr><td>Spikelets number/spike</td><td>0.23</td><td><strong>0.43</strong></td><td>-0.03</td><td>-0.02</td></tr><tr><td>Kernels number/fertile spikelets</td><td>0.22</td><td><strong>-0.42</strong></td><td>0.25</td><td>0.04</td></tr><tr><td>Kernels number/spikelet</td><td>0.19</td><td><strong>-0.41</strong></td><td>0.38</td><td>0.09</td></tr><tr><td>Kernels number/spike</td><td><strong>0.37</strong></td><td>-0.05</td><td><strong>0.39</strong></td><td>0.07</td></tr><tr><td>Spike length</td><td>0.17</td><td>0.29</td><td><strong>0.34</strong></td><td>-0.25</td></tr><tr><td>Fertile spikelets number/spike</td><td>0.29</td><td>0.29</td><td><strong>0.29</strong></td><td>0.03</td></tr><tr><td>1000 kernels weight</td><td><strong>0.33</strong></td><td>-0.03</td><td><strong>-0.44</strong></td><td>0.06</td></tr><tr><td>Kernel length</td><td><strong>0.30</strong></td><td>0.09</td><td>-0.28</td><td><strong>0.72</strong></td></tr><tr><td>Kernel width</td><td><strong>0.34</strong></td><td>-0.08</td><td>-0.35</td><td><strong>-0.48</strong></td></tr><tr><td>Eigenvalues</td><td></td><td></td><td></td><td></td></tr><tr><td>Variation explained (%)</td><td>42.7</td><td>26.8</td><td>14.9</td><td>4.5</td></tr><tr><td>Cumulative proportion of total variance</td><td>42.7</td><td>69.5</td><td>84.4</td><td>88.9</td></tr></tbody></table>

In bold, variable loading scores with the greatest loads on each component.

<footnote>(**, *** significant for *p* < 0.01; *p* < 0.001. ns = not significant).</footnote>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Principal component analysis (PCA) and cumulative proportion of total variance for the first and second PC for 136 accessions and 28 durum wheat varieties. Different colors are used for the six groups uncovered by the DAPC. Arrows indicate the loadings for each trait along the first two components.</figcaption></figure>

Principal components analysis (PCA) analysis based on the measured morpho-agronomic traits was also able to differentiate contrasting wheat groups as identified using the DAPC analysis. Genotypes that did cluster in C1 and C2 with the DAPC analysis were also included in the same group in PCA analysis and, according to the first PCA loadings, were early heading and short plant height, while genotypes in C3, C4, C5, and C6 showed late heading, tall plant height, and low kernel number/fertile spikelet and kernels number/spikelet (Figure 5).

**4. Discussion**

The recent history of durum wheat has been characterized by the replacement of landraces at first with selected pure line varieties and after with modern cultivars obtained by crosses between lines that did also include lines from mutagenesis, exotic germplasm, and different wheat species. Luckily ex situ collections preserve precious genetic resources collected over time that include germplasm from different genepools. Since in most cases the passport data do not include genetic information, their effective use in agriculture and breeding programs is very low.

The identification of useful genotypes is largely improved if the genetic structure of the germplasm collection is known. In the present study, we combined molecular genotyping and plant phenotyping with the comparison to well-known pedigreed varieties to describe the genetic structure of durum wheat germplasm collected in Southern Italy and preserved ex situ.

4.1. *Diversity of South Italian Durum Wheat Germplasm*

Overall, gene diversity for SSR markers was 0.607, and within the range of previous studies in durum wheat [9–11,47]. Gene diversity and allelic richness within ex situ accessions were higher than within pedigreed varieties. This indicates a substantial loss of genetic diversity over time [48] that might involve alleles valuable for plant improvement and future demands of producers and
consumers. Moreover, the gene diversity within accessions and within varieties is representative of different genepools. In fact, most modern varieties, that are the results of hybridization and also include newly introduced germplasm cluster in group C2 (n = 16) or in C1 (n = 2), while most accessions cluster in four different groups (C3, C4, C5, and C6) that include "landraces" or varieties from pure line selection largely representative of the Italian-Mediterranean genepool. In ex situ accessions a small increase of diversity was observed from Group 1 (n = 37 accessions collected from 1947 to 1950) to Group 2 (n = 77 accessions collected from 1973 to 1982). Group 2 accessions, which were included in all clusters, were "old" and "new" germplasm grown at the same time in the same area at the time when intense varietal substitution was in progress. As expected, the amount of collected diversity decreased in the following time (n = 22 accessions collected from 1983 to 2003) when most obsolete varieties had been already replaced. Clustering also showed that the choice of collecting sites was very careful to avoid the most recently introduced varieties. Out of 136 accessions, only 10 did cluster in C2 that included most "modern" varieties, and 23 in C1 which includes varieties derived from hybridization.

It has been argued that the selection pressure applied in breeding programs may have reduced the level of genetic diversity in durum wheat germplasm [49]; in our study modern varieties (1974–2007) did show a degree of genetic diversity lower than in the "old and intermediate" group (1915–1973).

Our findings on the elite germplasm agree with previous reports by Medini et al. [50], Reif et al. [48], and Figliuolo et al. [11], however, Maccaferri et al. [10] did show a progressive increase of the genetic basis in the elite durum wheat germplasm. Moreover, Martos et al. [19] and Laidò et al. [13] observed that the overall molecular diversity of durum wheat remained quite constant throughout genetic improvement occurring during the 20th century. The reduction of diversity might be explained by the "Green Revolution", which was characterized by breeding semidwarf varieties with a higher yielding potential due to an increased harvest index and better lodging tolerance, particularly under high fertilizer and water inputs [6]. These new high-yielding semidwarf modern varieties were based on a limited number of founder genotypes and rapidly dominated the wheat germplasm base [48,51]. This assumption may suggest that a high proportion of SSRs markers used in this study may be associated with chromosomal regions selected during breeding programs. These could be chromosomal regions harboring some agronomic relevant trait loci (such as those determining semidwarf habit or yield) that resulted in more uniform and stable modern varieties.

The now available high-density consensus maps of durum wheat coupled with sequencing information provide an advanced tool [52,53] for durum wheat fingerprinting. Indeed, the increasing knowledge of polymorphic marker distribution in the durum wheat genome will facilitate the selection of SSR makers useful to detect pattern of diversity associated to different gene pools and will improve the use of ex situ collections in breeding programs designed to exploit opportunities offered by the wider wheat gene pool.

## 4.2. *Analysis of Population Structure*

Discriminant Analysis of Principal Components and STRUCTURE analysis were very informative and complementary about the structure of the durum wheat ex situ germplasm collection and were useful to identify the most diverse accessions for utilization in low input agriculture and in breeding programs. For instance, accessions were included in all clusters, whereas modern varieties from hybridization, including in their pedigree the innovative semidwarf CIMMYT materials, were assigned to only two clusters (C1 and C2). These findings are consistent with the repeated use of a few founder genotypes that played a relevant role in the creation of the genetic basis of modern genetic pools, in turn becoming progenitors of new elite varieties and completely replacing traditional varieties or landraces [5,54]. The continued use of these genotypes made the gene pool smaller for all of the durum wheat elite varieties and resulted in the loss of genetic diversity [13]. Low levels of genetic diversity of Italian cultivars developed by Italian breeders were also observed by Kabbaj et al. [18]. The low level of genetic diversity can be explained as the combined result of frequent hybridization of a reduced
number of founders and the strong selection pressure for the same trait needed for Italian growing conditions and the requirements of the pasta industry [18].

New variability is needed to face the challenges of modern low input agriculture. Accessions from the ex situ collection that are genetically distant, as those included in different clusters, could be useful for this purpose. Pedigreed varieties were useful to identify the gene pool of origin of the accessions that did cluster in different groups. It should be noted that the varietal purity was higher for those released more recently, while genotypes associated with older releases might be different from accessions with the same name. In some cases, Saragolla, a brand new variety that has been released in 2004 with an old name, or variety and accessions with the same name as "Russello", might be selected as pure lines from the same "landrace" (Russello 329-S.97-S.G.7) [21], or the variety Senatore Cappelli, that now is very popular in organic agriculture, is genetically distinct from the two ex situ accessions with the same name, that cluster in the same group. The accession variety name must be also carefully verified: the Capeiti 8 variety is included in C1 as expected according to the pedigree, while the accession named "Capeiti", which was probably misnamed, did cluster in C2 and is genetically admixed. Also, the Trinakria variety (C2) and "Trinakria" accession (C1) were included in different clusters, but in this case, the mixed genotype was the variety. The Tripolino that we analyzed was probably not the "old" variety since it did cluster in C2 with most modern varieties while the synonymous variety Azizia did cluster in C5. Moreover, the STRUCTURE analysis revealed that it is actually an admixed genotype, thus suggesting that seeds of "pedigreed" varieties of "true" sources must be carefully tested for varietal purity. Two accessions and one variety with name Senatore Cappelli cluster in C3. The two accessions look fairly similar genetically and are also similar to Bidi that originated from the same "landrace", while the now grown Senatore Cappelli variety is genetically distinct and shows some degree of admixture with the old germplasm from Sicily and with C1 and C2 where the most advanced varieties are included. It seems that evolution through natural hybridization took place for Senatore Cappelli and that the pure line, recently reselected for distribution to farmers, is actually different from the old one. These results were also confirmed by phenotypic data that show Senatore Cappelli variety is different from both accessions and from Bidi: the plant is shorter and has later heading, longer head, less fertile spikelets, smaller kernels, and fewer kernels per spike.

In a recent study [18], a panel of 370 durum wheat genotypes, including 35 Italian genotypes (29 elite variety and 6 landraces), were genotyped using SNP markers. The Italian genotypes were attributed to seven subclusters. Three "landraces" did cluster with elite varieties indicating that these accessions were probably not "landraces" but cultivars that were not properly identified during the collecting missions [18]. For the other three Italian "landraces" no information about their origin was available in passport data so a direct comparison with the accessions used in this study was not possible. Out of twenty-nine Italian elite varieties studied by Kabbai et al. [18] only five were common to our study (Capeiti 8, Creso, Applo, Claudio, and Svevo). Capeiti 8, an "old" Italian variety from a cross between Syriacum × Mediterraneum groups (Eiti 6 × Cappelli), did cluster with lines and cultivars from the International Center for Agricultural Research in the Dry Areas (ICARDA) breeding program that include the old cultivar "Om Rabi" in their pedigree [18]. Also in our study was Capeiti 8, which did cluster with "old" genotypes selected from indigenous and exotic landraces. The remaining four varieties were included in three different clusters in Kabbai et al. [18], while in our study they did cluster in the same group with the "modern" dwarf genotypes released after 1974 that were selected from crosses that also included CYMMIT breeding lines.

Discriminant analysis of principal components can be also used to identify the alleles with the largest contributions to the discriminant functions, as an approach to detect putative patterns among the genes determining the group differentiation [40]. A plot of SSR allele contributions was used to identify alleles of major interest and, substantially, most of them were detected at loci associated to quantitative trait locus (QTLs) for phenological traits, as well as to resistance to wheat diseases and to grain qualitative parameters. Indeed, when analyzing the alleles' contributions to the first principal component, the gene *Vrn-B1*, affecting the vernalization response, was found to be tightly linked to the
marker *Xgwm408* [55,56]. Noteworthy, the genes of the *Vrn* system are the most important to determine the rate of generative development in wheat and, thereby, its flowering time [55–57]. Additionally, a QTL for Fusarium head blight (FHB) resistance was mapped proximal to marker *Xgwm282* [58], whereas the locus *Xgwm234* was found to be significantly associated with dough strength [59,60] and the marker *Xgwm285* was located at only 2.1 cM from the *tsn2* gene, conditioning resistance to race 3 of the fungus *Pyrenophora tritici-repentis* (Died.), causal agent of tan spot in durum wheat and inducing necrosis [61]. Since the inspection of alleles for the abovementioned markers showed that C1 and C3 were highly differentiated from C4 and C5, this could be an effective criterion for selection.

Similarly, when considering the alleles' contribution to the second principal component, we noticed that the marker *Xgwm193* encompasses the peak of the *QGpc.ndsu.6Bb* QTL for grain protein content. Thus, it could be useful to use this SSRs marker to select for the high-GPC allele in Marker-Assisted Selection programs [62,63]. Finally, the marker *Xgwm518* was found to be linked to the *QLr.caas-6BS.2* QTL for resistance to leaf rust caused by *Puccinia triticina* Westend [64] and the marker *Xgwm413* to be associated to a QTL for grain weight on the short arm of chromosome 1B [65]. These three markers differentiated C3 and C4 from C5, allowing to highlight the potential of the DAPC method to go beyond mere group delimitation and to identify accessions potentially useful in breeding programs.

Genetic groups identified by DAPC were also differentiated on the basis of morphological traits, as confirmed by PCA. Moreover, accessions and varieties assigned to the same cluster, although genetically related on the base of molecular markers, exhibited different morphological traits. These results could be helpful to select, especially from C3 and C4, early-heading landraces well adapted to the Southern Italian semi-arid conditions, where earliness allows farmers to better tackle the drought season [16]. Within each of the six DAPC clusters accessions plants were taller than varieties. Plant height could be a preventive measure to face weed competition, which remains one of the main problems in organic wheat crops since no herbicides are allowed [66]. Indeed, durum wheat landraces, characterized by tall plants and early vigor, could provide early groundcover which is vital to compete for light, water or nutrients. This trait represents a competitive advantage over early emerging weeds and their suppression [67]. Finally, higher values for yield components in ex situ accessions in C3, C4, and C5 could be due to higher N-use efficiency and their ability to produce high yields at low soil-N availability typical in organic and low-input agricultural systems [68–70]. Therefore, since wheat landraces have evolved mostly in environments with low nutrient availability [71], they represent a source of variation for selection of varieties adapted to cropping systems with low fertilizer input.

**5. Conclusions**

Both the morphological traits and SSR markers used in the present study were equally appropriate to provide a first overview of the genetic diversity levels and of the population structure within the whole durum wheat germplasm collection from Southern Italy.

Evolution through natural hybridization is suggested for old varieties recently reintroduced in agriculture and subjected to careful conservative selection by pure line.

Despite similar level of diversity for molecular markers in ex situ accessions and varieties DAPC did show that they cluster in different groups, thus representing different genepools capable to provide different sources of genes.

The analysis of diversity based only on molecular markers must be supported by passport data and or some phenotyping or knowledge in crops. In crops where plant breeding played a relevant role, social factors that promote seed exchange might be as significant as ecological factors in determining the distribution of genetic diversity and variation in such "ex situ" collections might be heavily structured.

The germplasm evaluation highlighted the existence of a broad genetic base in ex situ accessions and a narrowing of diversity in varieties due to breeding activities. The investigation of population structure suggested the genetic potential of landraces for the detection of unexplored sources of variation and allowed to identify groups of accessions differentiated both at molecular level and for
morphological, phenological, and qualitative traits not directly evaluated and potentially useful in
sustainable production systems.

The generated knowledge about the levels of diversity and population structure could be an
important contribution for parent selection in durum wheat breeding programs and for germplasm
management and conservation.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2073-4425/9/10/465/s1,
Figure S1: Bayesian Information Content (BIC) for increasing values of the number of clusters. The chosen number
of clusters was K = 6, Figure S2: STRUCTURE estimation of the number of subpopulations for K ranging from
2 to 10 by mean likelihoods (A) and Delta K values (ΔK) (B), Figure S3: Means and ranges of variation for
12 morphophysiological plant characters for six DAPC group clusters obtained from 136 durum wheat accession
(in black) and 28 varieties (in red), Table S1: List of durum wheat accessions used in the study. Seed source,
collecting date, geographical origin, and DAPC cluster attribution of each accession are reported, Table S2: List
of durum wheat varieties used in the study. Year of release, pedigree, breeder, and seed source of each variety
are reported, Table S3: Primers sequence and characteristics of the 44 SSR markers used to study diversity in a
durum wheat germplasm collection from Southern Italy, Table S4: ANOVA analysis for the morphological traits
measured on the wheat germplasm collection. Variation: between germplasm type (accessions and varieties)
and between accessions among germplasm type (*, **, *** significant for *p* < 0.05; *p* < 0.01; *p* < 0.001, Table S5:
Summary statistics with mean values, coefficient of variation (CV), minimum (Min), and maximum (Max) for
morpho-quantitative traits measured on 136 accessions and 28 varieties.

**Author Contributions:** S.M., T.G., and G.L. performed research; S.M. and G.L. analyzed data; J.L.D. provided materials; T.G., G.L., and P.S.Z. designed the research; and S.M. and T.G. wrote the paper. All authors read and approved the final manuscript.

**Funding:** This research was funded by SIR-MIUR grant number [RBSI14SJ81] and by the Basilicata Region (C.F. 80002950766) project "Azioni integrate per la salvaguardia e la valorizzazione di specie vegetali autoctone a rischio di estinzione" (Rural Development Programme 2007-2013, founded by Basilicata Region PSR- misura 214 azione 5, D.G.R. n. 1096/2012).

**Conflicts of Interest:** The authors declare no conflicts of interest.

**References**

1. Odong, T.L.; Jansen, J.; van Eeuwijk, F.; van Hintum, T.J.L. Quality of core collections for effective utilisation of genetic resources review, discussion and interpretation. *Theor. Appl. Genet.* **2013**, *126*, 289–305. [CrossRef] [PubMed]

2. Sharma, S.; Deo Upadhyaya, H.; Varshney, R.K.; Gowda, C.L.L. Pre-breeding for diversification of primary gene pool and genetic enhancement of grain legumes. *Front. Plant Sci.* **2013**, *4*, 309. [CrossRef] [PubMed]

3. Paux, E.; Sourdille, P.; Mackay, I.; Feuillet, C. Sequence-Based marker development in wheat: Advances and applications to breeding. *Biotechnol. Adv.* **2012**, *30*, 1071–1088. [CrossRef] [PubMed]

4. Ren, J.; Sun, D.; Chen, L.; You, F.M.; Wang, J.; Peng, Y.; Nevo, E.; Sun, D.; Luo, M.C.; Peng, J. Genetic diversity revealed by single nucleotide polymorphism markers in a worldwide germplasm collection of durum wheat. *Int. J. Mol. Sci.* **2013**, *14*, 7061–7088. [CrossRef] [PubMed]

5. Thuillet, A.C.; Bataillon, T.; Poirier, S.; Santoni, S.; David, J.L. Estimation of long-term effective population sizes through the history of durum wheat using microsatellite data. *Genetics* **2005**, *169*, 1589–1599. [CrossRef] [PubMed]

6. Evenson, R.E.; Gollin, D. Assessing the impact of the Green revolution. *Science* **2003**, *300*, 758–762. [CrossRef] [PubMed]

7. Newton, A.C.; Akar, T.; Baresel, J.P.; Bebeli, P.J.; Bettencourt, E.; Bladenopoulos, K.V.; Czembor, J.H.; Fasoula, D.A.; Katsiotis, A.; Koutis, K.; et al. Cereal landraces for sustainable agriculture. A review. *Agron. Sustain. Dev.* **2010**, *30*, 237–269. [CrossRef]

8. Hoisington, D.; Khairallah, M.; Reeves, T.; Ribault, J.M.; Skovmand, B.; Taba, S.; Warburton, M. Plant genetic resources: What can they contribute toward increased crop productivity? *Proc. Natl. Acad. Sci. USA* **1999**, *96*, 5937–5943. [CrossRef] [PubMed]

9. Moragues, M.; Moralejo, M.; Sorrells, M.E.; Royo, C. Dispersal of durum wheat [Triticum turgidum L. ssp. turgidum convar. durum (Desf.) MacKey] landraces across the Mediterranean basin assessed by AFLPs and microsatellites. *Genet. Resour. Crop Evol.* **2007**, *54*, 1133–1144. [CrossRef]
10. Maccaferri, M.; Sanguineti, M.C.; Donini, P.; Tuberosa, R. Microsatellite analysis reveals a progressive widening of the genetic basis in the elite durum wheat germplasm. *Theor. Appl. Genet.* **2003**, *107*, 783–797.
[CrossRef] [PubMed]

11. Figliuolo, G.; Mazzeo, M.; Greco, I. Temporal variation of diversity in Italian durum wheat germplasm. *Genet. Resour. Crop Evol.* **2007**, *54*, 615–626. [CrossRef]

12. Ruiz, M.; Giraldo, P.; Royo, C.; Villegas, D.; Aranzana, M.J.; Carrillo, J.M. Diversity and genetic structure of a collection of Spanish durum wheat landraces. *Crop Sci.* **2012**, *52*, 2262–2275. [CrossRef]

13. Laidò, G.; Mangini, G.; Taranto, F.; Gadaleta, A.; Blanco, A.; Cattivelli, L.; Marone, D.; Mastrangelo, A.M.; 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*, e67280. [CrossRef] [PubMed]

14. Sahri, A.; Chentoufi, L.; Arbaoui, M.; Ardisson, M.; Belqadi, L.; Birouk, A.; Roumet, P.; Muller, M.H. Towards a comprehensive characterization of durum wheat landraces in Moroccan traditional agrosystems: Analysing genetic diversity in the light of geography, farmer's taxonomy and tetraploid wheat domestication history. *BMC Evol. Biol.* **2014**, *14*, 264. [CrossRef] [PubMed]

15. Motzo, R.; Giunta, F. The effect of breeding on the phenology of Italian durum wheats: From landraces to modern cultivars. *Eur. J. Agron.* **2007**, *26*, 462–470. [CrossRef]

16. De Vita, P.; Li Destri Nicosia, O.; Nigro, F.; Platani, C.; Riefolo, C.; Di Fonzo, N.; Cattivelli, L. Breeding progress in morphophysiological, agronomical and qualitative traits of durum wheat cultivars released in Italy during the 20th century. *Eur. J. Agron.* **2007**, *26*, 39–53. [CrossRef]

17. Giunta, F.; Motzo, R.; Pruneddu, G. Trends since 1900 in the yield potential of Italian-bred durum wheat cultivars. *Eur. J. Agron.* **2007**, *27*, 12–24. [CrossRef]

18. Kabbaj, H.; Sall, A.T.; Al-Abdallat, A.; Geleta, M.; Amri, A.; Filali-Maltouf, A.; Belkadi, B.; Ortiz, R.; Bassi, F.M. Genetic diversity within a global panel of durum wheat (*Triticum durum*) landraces and modern germplasm reveals the history of alleles exchange. *Front. Plant Sci.* **2017**, *8*, 1277. [CrossRef] [PubMed]

19. Martos, V.; Royo, C.; Rharrabti, Y.; Garcia del Moral, L.F. Using AFLPs to determine phylogenetic relationships and genetic erosion in durum wheat cultivars released in Italy and Spain throughout the 20th century. *Field Crops Res.* **2005**, *91*, 107–116. [CrossRef]

20. De Cillis, E. *I Grani d'Italia*; Tipografia della Camera dei Deputati: Roma, Italy, 1927; p. 173.

21. Deidda, M.; Motzo, R.; Giunta, F.; Fois, S. *Evoluzione Varietale e Qualità in Frumento Duro (Triticum turgidum subsp. Durum)*: Dalle Vecchie Popolazioni Alle Attuali Cultivar; Dipartimento di Scienze agronomiche e Genetica vegetale agraria, Facoltà di Agraria, Università degli Studi di Sassari: Sassari, Italy, 2001; p. 3.

22. Grignac, P. Contribution à l'étude de *Triticum durum* Desf. Thése Doctorat, Université de Toulouse, Toulouse, France, 1965; p. 151.

23. Bozzini, A. Genetics and breeding of durum wheats. *Genet. Agraria* **1970**, *24*, 145–193.

24. Porceddu, E.; De Pace, C.; Spagnoletti Zeuli, P.L.; Scarascia Mugnozza, G.T. Genetic variability in segregating populations of durum wheat (*Triticum durum*). *Genet. Agraria* XXXIV **1980**, *34*, 182–183.

25. Blanco, A.; Scarascia-Mugnozza, G.T. Messapia. A new early and productive variety of durum wheat. *L'Inf. Agr. XLV* **1989**, *39*, 61–63.

26. Blanco, A.; Simeone, R.; Resta, P. The addition of *Dasypirum villosum* (L.) Candargy chromosome in durum wheat (*Triticum durum* Desf.). *Theor. Appl. Genet.* **1987**, *74*, 328–333. [CrossRef] [PubMed]

27. Blanco, A.; Giorgi, B.; Perrino, P.; Simeone, R. Genetic resources and breeding of durum wheat quality. *Agricoltura e Ricerca* **1991**, *144*, 41–58.

28. Hammer, K.; Laghetti, G.; Cifarelli, S.; Perrino, P. Collection of plant genetic resources in Italy 1989. *Die Kulturpflanze* **1990**, *38*, 311–323. [CrossRef]

29. Somers, D.J.; Isaac, P.; Edwards, K. A high-density microsatellite consensus map for bread wheat (*Triticum aestivum* L.). *Theor. Appl. Genet.* **2004**, *109*, 1105–1114. [CrossRef] [PubMed]

30. Korzun, V.; Börner, A.; Worland, A.J.; Law, C.N.; Röder, M.S. Application of microsatellite markers to distinguish inter-varietal chromosome substitution lines of wheat (*Triticum aestivum* L.). *Euphytica* **1997**, *95*, 149–155. [CrossRef]

31. Röder, M.S.; Korzun, V.; Wendehake, K.; Plaschke, J.; Tixer, M.H.; Leroy, P.; Ganal, M.W. A microsatellite map of wheat. *Genetics* **1998**, *149*, 2007–2023. [PubMed]
32. Sourdille, P.; Gandon, B.; Chiquet, V.; Nicot, N.; Somers, D.; Murigneux, A.; Bernard, M. Wheat Genoplante SSR Mapping Data Release: A New Set of Markers and Comprehensive Genetic and Physical Mapping Data. *Graigenes*. 2004. Available online: http://wheat.pw.usda.gov/GG2/index.shtml (accessed on 15 March 2018).

33. Doyle, J.J.; Doyle, J.L. A rapid DNA isolation procedure from small quantities of fresh leaf tissue. *Phytochem. Bull*. **1987**, *19*, 11–15.

34. International Board for Plant Genetic Resources (IPGRI). *Descriptors for Wheat (Revised)*; International Plant Genetic Resources Institute: Rome, Italy, 1985.

35. Nei, M. Estimation of average heterozygosity and genetic distance from a small number of individuals. *Genetics* **1978**, *89*, 583–590. [PubMed]

36. Belkhir, K.; Borsa, P.; Chikhi, L.; Raufaste, N.; Bonhomme, F. GENETIX 4.02, Logiciel Sous Windows™ pour la Génétique des Populations; Laboratoire Génome, Populations, Interactions, CNRS UMR 5000, Université de Montpellier II: Montpellier, France, 2001.

37. El Mousadik, A.; Petit, R.J. 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. [CrossRef] [PubMed]

38. Goudet, J. FSTAT, a Program to Estimate and Test Gene Diversities and Fixation Indices (Version 2.9.3). 2002. Available online: http://www2.unil.ch/popgen/softwares/fstat.htm (accessed on 15 March 2018).

39. Yang, M.H.; Fu, Y.B. An R version of FPTest for testing differences in allelic count. *Conserv. Genet. Resour.* **2017**, *9*, 313–318. [CrossRef]

40. 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. [CrossRef] [PubMed]

41. Jombart, T. adegenet: A R package for the multivariate analysis of genetic markers. *Bioinformatics* **2008**, *24*, 1403–1405. [CrossRef] [PubMed]

42. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. *Genetics* **2000**, *155*, 945–959. [CrossRef] [PubMed]

43. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. *Mol. Ecol.* **2005**, *14*, 2611–2620. [CrossRef] [PubMed]

44. Earl, D.A.; von Holdt, B.M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. *Conserv. Genet. Resour.* **2012**, *4*, 359–361. [CrossRef]

45. Weir, B.; Cockerham, C.C. Estimating F statistics for the analysis of population structure. *Evolution* **1984**, *38*, 1358–1370. [CrossRef] [PubMed]

46. Institute Inc. SAS. SAS/STAT®9.2. User's Guide, 2nd ed.; SAS Institute Inc.: Cary, NC, USA, 2009.

47. Al-Khanjari, S.; Hammer, K.; Buerkert, A.; Röder, M.S. Molecular diversity of Omani wheat revealed by microsatellites: I. Tetraploid landraces. *Genet. Resour. Crop Evol.* **2007**, *54*, 1291–1300. [CrossRef]

48. Reif, J.C.; Zhang, P.; Dreisigacker, S.; Warburton, M.L.; van Ginkel, M.; Hoisington, D.; Bohn, M.; Melchinger, A.E. Wheat genetic diversity trends during domestication and breeding. *Theor. Appl. Genet.* **2005**, *110*, 859–864. [CrossRef] [PubMed]

49. Vellvè, R. The decline of diversity in European agriculture. *Ecologist* **1993**, *23*, 64–69.

50. Medini, M.; Hamza, S.; Rebai, A.; Baum, M. Analysis of genetic diversity in Tunisian durum wheat cultivars and related wild species by SSR and AFLP markers. *Genet. Resour. Crop Evol.* **2005**, *52*, 21–31. [CrossRef]

51. Dalrymple, D.G. *Development and Spread of High-Yielding Wheat Varieties in Developing Countries*, 7th ed.; US Agency for International Development: Washington, DC, USA, 1986.

52. Maccaferri, M.; Cane', M.A.; Sanguineti, M.C.; Salvi, S.; Colalongo, M.C.; Massi, A.; Clarke, F.; Knox, R.; Pozniak, C.J.; Clarke, J.M.; et al. A consensus framework map of durum wheat (*Triticum durum* Desf.) suitable for linkage disequilibrium analysis and genome-wide association mapping. *BMC Genom.* **2014**, *15*, 873. [CrossRef] [PubMed]

53. Maccaferri, M.; Ricci, A.; Salvi, S.; Milner, S.G.; Noli, E.; Martelli, P.L.; Casadio, R.; Akhunov, E.; Scalabrin, S.; Vendramin, V.; et al. A high-density, SNP-based consensus map of tetraploid wheat as a bridge to integrate durum and bread wheat genomics and breeding. *Plant Biotechnol. J.* **2015**, *13*, 648–663. [CrossRef] [PubMed]

54. Porceddu, E. Varietal evolution and current problems of breeding winter cereals. *Rivista di Agronomia* **1987**, *21*, 33–54.
55. Tóth, B.; Galiba, G.; Fehěr, E.; Sutka, J.; Snape, J.W. Mapping genes affecting flowering time and frost resistance on chromosome 5B of wheat. *Theor. Appl. Genet.* **2003**, *107*, 509–514. [CrossRef] [PubMed]

56. Snape, J.W.; Butterworth, K.; Whitechurch, E.; Worland, A.J. Waiting for fine times: Genetics of flowering time in wheat. *Euphytica* **2001**, *119*, 185–190. [CrossRef]

57. Leonova, I.; Pestsova, E.; Salina, E.; Efremova, T.; Röder, M.S.; Börner, A.; Fischbeck, G. Mapping of the *Vrn-B1* gene in *Triticum aestivum* using microsatellite marker. *Plant Breed.* **2003**, *122*, 209–212. [CrossRef]

58. Jia, G.; Chen, P.; Qin, G.; Bai, G.; Wang, X.; Wang, S.; Zhou, B.; Zhang, S.; Liu, D. QTLs for Fusarium head blight response in a wheat DH population of Wangshuibai/Alondra's'. *Euphytica* **2005**, *146*, 183–191. [CrossRef]

59. Zanetti, S.; Winzeler, M.; Feuillet, G.; Keller, B.; Messmer, M. Genetic analysis of bread-making quality in wheat and spelt. *Plant Breed.* **2001**, *120*, 13–20. [CrossRef]

60. Arbelbide, M.; Bernardo, R. Mixed-model QTL mapping for kernel hardness and dough strength in bread wheat. *Theor. Appl. Genet.* **2006**, *112*, 885–890. [CrossRef] [PubMed]

61. Singh, P.K.; Gonzalez-Hernandez, J.L.; Mergoum, M.; Ali, S.; Adhikari, T.B.; Kianian, S.F.; Elias, E.M.; Hughes, G.R. Identification and molecular mapping of a gene conferring resistance to *Pyrenophora tritici-repentis* race 3 in tetraploid wheat. *Phytopathology* **2006**, *96*, 885–889. [CrossRef] [PubMed]

62. Khan, I.A.; Procunier, J.D.; Humphreys, D.G.; Tranquilli, G.; Schlatter, A.R.; Marcucci-Poltri, S.; Frohberg, R.; Dubcovsky, J. Development of PCR based markers for a high grain protein content gene from *Triticum turgidum* ssp. *dicoccoides* transferred to bread wheat. *Crop Sci.* **2000**, *40*, 518–524. [CrossRef]

63. Olmos, S.; Distelfeld, A.; Chicaiza, O.; Schlatter, A.R.; Fahima, T.; Echenique, V.; Dubcovsky, J. Precise mapping of a locus affecting grain protein content in durum wheat. *Theor. Appl. Genet.* **2003**, *107*, 1243–1251. [CrossRef] [PubMed]

64. Ren, Y.; Li, Z.; He, Z.; Wu, L.; Bai, B.; Lan, C.; Wang, C.; Zhou, G.; Zhu, H.; Xia, X. QTL mapping of adult-plant resistances to stripe rust and leaf rust in Chinese wheat cultivar Bainong 64. *Theor. Appl. Genet.* **2012**, *125*, 1253–1262. [CrossRef] [PubMed]

65. Mir, R.R.; Kumar, N.; Jaiswal, V.; Girdharwal, N.; Prasad, M.; Balyan, H.S.; Gupta, P.K. Genetic dissection of grain weight in bread wheat through quantitative trait locus interval and association mapping. *Mol. Breed.* **2012**, *29*, 963–972. [CrossRef]

66. Van Bueren, E.L.; Jones, S.S.; Tamm, L.; Murphy, K.M.; Myers, J.R.; Leifert, C.; Messmer, M.M. The need to breed crop varieties suitable for organic farming, using wheat, tomato and broccoli as examples: A review. *NJAS Wagening. J. Life Sci.* **2011**, *58*, 193–205. [CrossRef]

67. Van Bueren, E.L.; Struik, P.C.; Jacobsen, E. Ecological concepts in organic farming and their consequences for an organic crop ideotype. *NJAS Wagening J. Life Sci.* **2002**, *50*, 1–26. [CrossRef]

68. Anderson, W.K.; Seymour, M.; D'Antuono, M.F. Evidence for differences between cultivars in responsiveness of wheat to applied nitrogen. *Aust. J. Agric. Res.* **1991**, *42*, 363–377. [CrossRef]

69. Foulkes, M.J.; Sylvester-Bradley, R.; Scott, R.K. Evidence for differences between winter wheat cultivars in acquisition of soil mineral nitrogen and uptake and utilization of applied fertilizer nitrogen. *J. Agric. Sci.* **1998**, *130*, 29–44. [CrossRef]

70. Dawson, J.C.; Murphy, K.M.; Jones, S.S. Decentralized selection and participatory approaches in plant breeding for low-input systems. *Euphytica* **2008**, *160*, 143–154. [CrossRef]

71. Gioia, T.; Nagel, K.A.; Beleggia, R.; Fragasso, M.; Ficco, D.B.M.; Pieruschka, R.; De Vita, P.; Fiorani, F.; Papa, R. Impact of domestication on the phenotypic architecture of durum wheat under contrasting nitrogen fertilization. *J. Exp. Bot.* **2015**, *66*, 5519–5530. [CrossRef] [PubMed]

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