Skip to main content

Unlocking the story of resistance to Zymoseptoria tritici in Tunisian old durum wheat germplasm based on population structure analysis



Septoria tritici blotch (STB) remains a significant obstacle to durum wheat cultivation on a global scale. This disease remains a challenge for farmers, researchers, and breeders, who are collectively dedicated to reduce its damage and improve wheat resistance. Tunisian durum wheat landraces have been recognized as valuable genetic ressources that exhibit resistance to biotic and abiotic stresses and therefore play a crucial role in breeding program aimed at creating new wheat varieties resistant to fungal diseases as STB, as well as adapted to climate change constraints.


A total of 366 local durum wheat accessions were assessed for resistance to two virulent Tunisian isolates of Zymoseptoria tritici Tun06 and TM220 under field conditions. Population structure analysis of the durum wheat accessions, performed with 286 polymorphic SNPs (PIC > 0.3) covering the entire genome, identified three genetic subpopulations (GS1, GS2 and GS3) with 22% of admixed genotypes. Interestingly, all of the resistant genotypes were among GS2 or admixed with GS2.


This study revealed the population structure and the genetic distribution of the resistance to Z. tritici in the Tunisian durum wheat landraces. Accessions grouping pattern reflected the geographical origins of the landraces. We suggested that GS2 accessions were mostly derived from eastern Mediterranean populations, unlike GS1 and GS3 that originated from the west. Resistant GS2 accessions belonged to landraces Taganrog, Sbei glabre, Richi, Mekki, Badri, Jneh Khotifa and Azizi. Furthermore, we suggested that admixture contributed to transmit STB resistance from GS2 resistant landraces to initially susceptible landraces such as Mahmoudi (GS1), but also resulted in the loss of resistance in the case of GS2 susceptible Azizi and Jneh Khotifa accessions.

Peer Review reports


Durum wheat (Triticum. turgidum L. ssp. durum) is a tetraploid species, originated in the Fertile Crescent about 10,000 BP. It evolved from the domestication of hulled tetraploid wheat subspecies as emmer (Triticum turgidum L. ssp. dicoccum) in the eastern Mediterranean, notably at the mountains of the Fertile Crescent (Iran, Turkey, Syria and Jordan) and at the Tigris and Euphrates basin [1,2,3,4,5,6]. Thereafter, the geographical expansion of durum wheat had closely followed human migration [4], where Phoenicians, Greeks and Romans contributed crucially in the pervasion and adoption of durum wheat cultivation around the Mediterranean Basin [5]. During their migration, domesticated wheat populations underwent strong natural and human selection processes, after which they have adapted specifically to local environments and developed to become landraces [7, 8]. Howbeit, durum wheat followed two dispersal pathways from its area of origin in the Mediterranean Basin, over the north side via Southern Europe (Turkey, Greece, and Italy) and over the south side via North Africa [1, 9]. The spread of durum wheat populations occurred by land through the Balkans and by maritime route through the Mediterranean Sea [5]. Along these pathways, wheat landraces have traced a complex history of dissemination, adaptation and genetic differentiation in time and space [10]. Therefore, the eastern–western dispersal of the Mediterranean durum landraces favored their divergence into genetically distinct groups. In fact, Moragues et al. [9] classified a collection of 63 durum wheat landraces from the Mediterranean basin in a north dispersal group (including south European and south-western Asian landraces) and a south dispersal group (including landraces from North Africa and the Iberian Peninsula). Soriano et al. [11] also structured a collection of 172 durum wheat landraces from 21 Mediterranean countries into four genetic populations related to their geographical origin, namely eastern Mediterranean, eastern Balkans and Turkey, western Balkans and Egypt, and western Mediterranean.

Durum wheat cultivation history in North Africa and particularly Tunisia, involved the intervention of Phoenicians importing wheat from Lebanon to Carthage, along with the development of Carthage trade maritime activity in the Mediterranean Sea favoring seed exchanges between Tunisia and the Mediterranean countries. North African landraces were also introduced and diffused by Romans, who greatly influenced durum wheat cultivation in this area by setting up modernized irrigation systems [5, 12, 13]. Recently, Robanna et al. [13] studied the structure of six Tunisian durum wheat landraces, and reported their genetic similarity with landraces from North African countries and landraces from Greece, Italy, and Lebanon. Durum wheat was prevalent and well established in North Africa in the classical times [5]. Accordingly, North Africa and Abyssinian regions are considered as secondary centers of diversity for durum wheat [14, 15]. Ren et al. [16] also suggested that North Africa should be considered as a microcenter of wheat diversity. Tunisia, being part of the secondary center of diversity for durum wheat, has a rich collection of old local durum wheat landraces [14]. To date, around 40 old durum wheat landraces are known in Tunisia that were morphologically characterized and classified by Bœuf [14] and Dghais et al. [17]. Several studies emphasised high levels of genetic and agro-morphological diversity [13, 18,19,20,21,22], phenological features [23,24,25,26] and resistance to biotic and abiotic stresses i.e. drought, heat and fungal diseases [21, 27,28,29] within the Tunisian old durum wheat germplasm.

Recently, Ben Krima et al. [30] showed a complex structure of 14 Tunisian durum wheat populations that was not entirely related to their geographic origin and variety name. However, Ouaja et al. [22] identified a strong correlation between the genetic structure of 11 Tunisian durum wheat landraces and their morphological characterisation and nomenclature. Therefore, various interacting factors were reported that have influenced the structure and evolutionary dynamics of durum wheat in Tunisia and the Maghreb region overall, among which complex selection trajectory, the significance of variety names, the occurrence of heterogeneous mixtures within populations, local adaptation, local and regional exchanges between farmers and, loss and misidentification [22, 30, 31]. After the Green revolution, old durum wheat landraces were mainly grown and managed by smallholder farmers under low-input traditional agrosystems in the marginal areas of Mediterranean region, notably in southern Europe and North Africa [4, 11] as they were progressively abandoned from the early 1970s and replaced by improved genetically uniform modern varieties/cultivars [32, 33]. In Tunisia, durum wheat landraces are still cultivated by low-input farmers, in northern and central mountainous areas, under traditional farming systems. These landraces, transmitted by farmers from one generation to the next, are designated by a variety name linked to a historical origin, regional location and specific phenotypic characteristics [22, 30].

A high diversity of Tunisian durum accessions has been observed using morphological descriptors and biochemical markers [22, 34,35,36,37]. The genetic diversity of Tunisian durum germplasm was also investigated using different molecular markers such as AFLP and SSR markers, which allowed us to study the genetic variation among and within Tunisian landraces and modern cultivars [20, 22, 38]. Nowdays new high-throughput genotyping technologies such as single nucleotide polymorphism (SNP) arrays or genotyping-by-sequencing (GBS) become a procedure of choice. In fact, based on genotyping, several studies were conducted to analyse genetic diversity and the genetic structure of durum wheat landraces and modern cultivars in the mediteraneen basin and Ethiopia [13, 39,40,41].

Septoria Tritici Blotch (STB) caused by the fungus Zymoseptoria tritici (Z. tritici) (Desm.) (formerly Mycosphaerella graminicola) is currently considered among the most damaging and worldwide distributed fungal disease on cultivated wheat [42, 43]. The appearance of Z. tritici as a host-specialized wheat pathogen occurred about 10,500 years ago via host tracking throughout the time of wheat domestication [44,45,46]. Comparative genomic analysis between Z. tritici and its close relatives highlighted strong adaptive evolution of Z. tritici in relation to specialization on wheat [47]. Howbeit, Stukenbrock et al. [44] demonstrated that wheat-adapted Z. tritici was derived from an ancestral population infecting wild grasses in the Middle East and that the domestication of an agricultural crop was concomitantly accompanied by the domestication of a fungal pathogen. Accordingly, the Fertile Crescent is considered a hotspot of Z. tritici genetic diversity [44, 48]. Therefore, wheat landraces and their wild relatives from the Fertile Crescent, having co-evolved for a long time with Z. tritici, must harbor the greatest diversity for resistance to STB [49]. Moreover, a host species specialization was highlighted in Z. tritici populations to either bread or durum wheat [50,51,52] making bread wheat-derived Z. tritici strains not suitable to decipher Stb genes in durum wheat, as studies into the genetic basis of STB resistance were entirely based on the Z. tritici-bread wheat pathosystem [53,54,55,56,57,58].

In Tunisia and under suitable environmental conditions for infection, STB causes considerable yield losses up to 50–60% [59, 60]. The introduction of the modern cultivar ‘Karim’ in 1980 displaced the cultivation of landraces, resulting in a reduction of the genetic diversity (genetic erosion) and therefore, enhancing the susceptibility to STB [61, 62]. Durum wheat landraces, characterized by a substantial level of genetic diversity [11, 20, 22], represent the main sources of resistance to Z. tritici to be incorporated into breeding programs for a sustainable STB disease control and management [21, 63, 64]. The present study consisted on (a) analyzing the genetic structure of 366 Tunisian durum wheat accessions belonging to 13 landraces using 286 single nucloetoide polymorphic markers (SNPs), (b) evaluating the resistance of these accessions to two durum-wheat derived Z. tritici isolates (Tun06 and TM220) and, (c) relating genotyping and disease screening data of the Tunisian durum wheat accessions to describe the life history of STB resistance in local durum wheat landraces.


Genetic structure of durum wheat landraces

The genetic structure of the durum wheat accessions was determined using the Bayesian clustering model implemented in STRUCTURE [65]. The maximum likelihood (LnP (K)) and delta K (ΔK) methods [66]) showed that the most likely number of genetic subpopulations (K) was 3 (Fig. 1, a and b). The inferred population structure at K = 3 showed that 78% of the genotypes have a membership coefficient higher than 0.7 to one of the genetic subpopulations (GS1, GS2 and GS3), the rest were admixed. GS1, GS2, GS3 and admixed genotypes represented 20%, 45%, 13% and 22% of the entire collection, respectively. Among admixed genotypes, 53%, 27% and 20% were admixed between G1-G2, G1-G3 and G2-G3, respectively (Fig. 1c, Table S1).

Fig. 1
figure 1

Genetic structure analysis of 366 Tunisian durum wheat accessions genotyped with 286 SNP markers: (a) Plot of mean posterior probability (ln P(D)) values per cluster (K); (b) delta-K analysis of Ln P(D), for K ranging from 1 to 20; (c) Membership coefficient bar plot displaying genetic structure at K = 3 from STRUCTURE [65]. Each genotype is represented by a vertical line

PCAs were performed using the 286 SNPs on the 366 genotypes (Fig. 2). Axes 1 and 2 of the PCAs accounted for 24.33% and 16.54% of the total genetic variation, respectively. Figure 2a showed that PCA grouping corroborated the genetic structure inferred by STRUCTURE, pointing a clear differentiation between GS1, GS2 and GS3. Admixed genotypes were essentially distributed between GS1 and GS2 and between GS2 and GS3, reflecting ongoing hybridization and allele exchanges between these groups. Pairwise Fst values also showed considerable genetic differentiations between GS1, GS2 and GS3 (Table 1a). The highest Fst value (0.751) was observed between GS1 and GS3. Fst values between GS1 and GS2 and between GS2 and GS3 were both around 0.400. Similarly, Nm indices between GS1 and GS2 and between GS2 and GS3 were close; while Nm between GS1 and GS3 was the lowest (0.078) reflecting an almost absence of genetic exchange between these two populations. The AMOVA (Table 2) was consistent with pairwise Fst and Nm analysis revealing that the genetic variation between subpopulations (61%) was higher than the variation within subpopulation (39%).

Fig. 2
figure 2

Principal component analysis plot of 366 Tunisian durum wheat accessions belonging to 13 landraces genotyped with 286 SNP markers under GenAlEx [101], color-coded by (a) genetic structure (GS1, GS2, GS3 and admixed genotypes) as identified by STRUCTURE [58] (for K = 3); (b) landraces denomination; and (c) resistance class to Z. tritici isolate Tun06 and (d) isolate TM220

Table 1a Pairwise Fst values (under the diagonal) and gene flow (Nm) (above the diagonal) calculated under GenAlEx 6.501 [101] between the genetic subpopulations as defined by STRUCTURE [65] (at K = 3)*
Table 1b Pairwise Fst values (under the diagonal) and gene flow (Nm) (above the diagonal) calculated under GenAlEx 6.501 [101] between landraces*
Table 2 Analysis of molecular variance (AMOVA) of the Tunisian durum wheat accessions performed under GenAlEx 6.501 [101], by genetic subpopulation as defined by STRUCTURE [65] (at K = 3) and by landrace *

GS1 was solely composed by the landrace Mahmoudi. GS2 was composed of the landraces Azizi, Badri, Biskri, Jneh Khotifa, Mekki, Richi, Sbei glabre and Taganrog. GS3 was composed by the landrace Biada and two modern cultivars (Karim and Razzek). Twenty accessions of Mahmoudi, 11 of Jneh Khotifa, 8 of Richi, one Azizi and one Roussia accession were admixed between GS1 and GS2. Admixed genotypes between GS2 and GS3 included one accession of Azizi, 4 of Biada, 2 of Biskri, 5 of Souri, 5 of Roussia and 2 modern cultivars (Om Rabia andSalim). Finally, one accession of Mahmoudi, 2 of Biada, 8 of Bidi, 6 of Souri, 2 of Roussia and 2 modern cultivars (Nasr and Khiar) were admixed between GS1 and GS3 (Fig. 2b, Table S1).

In addition, population structure was investigated from K = 13 (corresponding to the number of landraces studied) to K = 20 (Table S2). At K = 13, landraces Azizi, Badri, Biada, Bidi, Jneh Khotifa, Mekki, Richi and Sbei glabre were assigned to separate genetic subpopulations. However, landraces Souri and Roussia were grouped in the same genetic subpopulation. Landraces Mahmoudi and Biskri were both divided into two genetic subpopulations. Taganrog landrace was entirely composed of admixed genotypes (between Sbei glabre, Jneh Khotifa, Azizi and Badri). The STRUCTURE assignment that matched exactly the landraces denomination was obtained for K = 15. Indeed, at K = 15, all 13 landraces were discriminated and assigned to genetically distinct subpopulations, except for Mahmoudi and Biskri, which were both divided into two genetic subpopulations as for K = 13. In fact, Mahmoudi subpopulations corresponded to two morphological types of Mahmoudi namely Mahmoudi-122 and Mahmoudi-986. Mahmoudi-122 had larger grain size and a relaxed spike comparing to Mahmoudi-986. Biskri subpopulations corresponded to two morphological types of Biskri namely Biskri-Ac1 and Biskri-glabre (Figure S1). Pairwise Fst values based on landrace varied from 0.020 to 1 and most all the durum wheat landraces were genetically differentiated (Table 1b). Nonetheless, the lowest Fst value (0.020) was observed between Biada and the modern varieties, followed by the Fst value between Taganrog and Sbei glabre (0.184). These varieties shared almost the same agro-morphological characteristics. The lowest gene flow values were recorded between the modern varieties and all landraces, except for Biada. In addition, AMOVA (Table 2) showed that the genetic variation among landraces (81%) was higher than the variation within landraces (19%).

The Unweighted Pair Group Method with Arithmetic Average (UPGMA) tree, generated with 286 SNPs data of 366 genotypes, differenciated three subclusters that were mainly in agreement with the genetic grouping defined by STRUCTURE (Fig. 3). Subcluster « SC-I » grouped all Mahmoudi GS1 and admixed genotypes, jointly with all GS2 and admixed genotypes of Jneh Khotifa, 2 Azizi admix, 3 Biada admix, one Roussia admix, 8 Richi admix and one Richi GS2. Subcluster « SC-II » included all Azizi, Badri, Biskri, Mekki, Sbei glabre, Taganrog and Richi GS2 genotypes, along with 5 Souri admix and 3 Mahmoudi admix. Subcluster « SC-III » grouped all GS3 genotypes notably Biada, together with 8 admixed genotypes of Bidi, 6 of Roussia, 6 of Souri and 2 of Mahmoudi.

Fig. 3
figure 3

Unweighted Pair Group Method with Arithmetic Average (UPGMA) tree of 366 Tunisian durum wheat accessions genotyped with 286 SNP markers. Genotype names are labelled as listed in Table S3 and are color-coded by resistance class to Z. tritici isolate Tun06 . Branches are color-coded according to STRUCTURE results at K = 3 [65]. Three subclusters were identified and labelled SC-I, SC-II and SC-III.

Diversity analysis (Table 3) of the genetic subpopulations showed that GS2 presented higher diversity indices (I = 0.519; He = 0.353) than GS1 and GS3. The overall gene flow between the genetic subpopulations (Nm = 0.414) reflected a low gene exchange between the three groups. Diversity indices of the accessions based on their landrace nomenclature were the lowest for Badri and Bidi (I = 0.00; He = 0.00) and the highest for Roussia (I = 0.392; He = 0.259). These observations are consistent with pairwise Fst values (Table 1b and AMOVA (Table 2) indicating that landraces were genetically distinct (low intra-population variability).

Table 3 Diversity indexes of 366 Tunisian durum wheat accessions grouped by genetic subpopulations as defined by STRUCTURE [65] (at K = 3) and by landraces

Distribution of the resistance to Zymoseptoria tritici among subpopulations

The 366 durum wheat accessions were evaluated for their resistance to two Tunisian Z. tritici isolates Tun06 and TM220. Based on disease scoring, the landraces were categorized into three classes of response to Z. tritici (R, I and S) as defined by Ouaja et al. [21]. Overall, 20%, 15.8% and 63.8% of the collection were R, I and S genotypes to isolate Tun06, respectively. In addition, 32.8%, 8.9% and 57.9% of the collection were R, I and S genotypes to isolate TM220, respectively.

The ANOVA (Table 4) revealed that the genetic subpopulation effect is highly significant, indicating that the variation in the responses to Z. tritici infection was consistent with the population genetic structure. However, a non-significant effect of isolates on the RAUDPC scores was revealed. Genetic subpopulations x isolates interaction was also non-significant. Accordingly, isolates Tun06 and TM220 did not significantly varied in their severity (RAUDPC scores) towards wheat accessions.

Table 4 Analyses of variance (ANOVA) performed under R 3.3.2 [97] on the Relative Area Under the Disease Progress Curve (RAUDPC) on the Tunisian durum wheat accessions genotyped with 286 SNPs. Main effects of genetic subpopulations, isolates and genetic subpopulations x isolates interaction were investigated

The distribution of R, I and S classes by genetic subpopulation was uneven. However, the distribution of each class of resistance differed slightly between isolates Tun06 and TM220 (Fig. 2.c, d, Table S1). GS1 and GS3 were mainly formed by accessions of class S. Infact 93% and 85% of the accessions were susceptible to Tun06 within GS1 and GS3, respectively. For the isolate TM220, 64% and 34% of the accessions were susceptible within GS1 and GS3 respectively. Interestingly, 44% and 31% of resistant accessions to Tun06 belonged to GS2 and admixed genotypes with GS2, respectively. Similarly, 48% and 26% of the resistant to TM220 were GS2 and admixed genotypes with GS2 respectively. Genotypes with intermediate responses to Tun06 and TM220 isolates were also mainly among GS2 or admixed with GS2. While GS2 comprised less than 35% of susceptible accessions. About 44% of admixed genotypes (mostly GS1-GS2 and GS2-GS3) were susceptible.

The PCA of 366 Tunisian durum wheat accessions sorted by classes of response to isolate Tun06 (Fig. 2c) showed that resistant genotypes were mainly laying among GS2 and admixed genotypes GS1-GS2 and GS2-GS3. Susceptible genotypes were spread over the three subpopulations. UPGMA tree (Fig. 3) showed that the three subclusters grouped genotypes of different resistance classes. All susceptible Mahmoudi GS1 genotypes and mostly resistant admixed genotypes between GS1 and GS2 (Mahmoudi, Jneh Khotifa, Badri and Richi) were in SC-I. SC-II was essentially composed of resistant GS2 landraces (Taganrog, Sbei glabre, Richi, Mekki, Badri and Azizi). However, 100% of Biskri accessions, 54% of Azizi, 7% of Mekki, 5% of Sbei glabre and 45% of Souri admix were within SC-II and showed susceptibility. SC-III is entirely composed of susceptible genotypes of Biada GS3, Biada admix, Bidi admix, Souri admix and Roussia admix.


Old durum wheat germplasm represents a precious genetic heritage. Understanding the genetic and phenotypic structure of old local landraces will help retracing their life history, resistance and durability of use. Durum wheat population structure study also help deciphering new sources of resistance to cope with challenging abiotic and biotic stresses, notably STB, one of the most devastating fungal disease on durum wheat crop. In the present study, we genotyped 366 Tunisian durum wheat accessions belonging to 13 old Tunisian landraces [22], collected from three central and two southern regions in Tunisia, using 286 SNPs derived from a High-density 90 K wheat SNP array [67]. This study revealed the population structure and the genetic distribution of the resistance to Z. tritici in the Tunisian durum wheat landraces.

Identification of a population structure in the tunisian durum wheat landraces related to their potential introduction pathways

The genetic structure of 366 durum wheat accessions was investigated using 286 SNPs. Three major genetic supbpopulations (K = 3: GS1, GS2 and GS3) were obtained under STRUCTURE [65]. At K = 15, we were able to attribute each landrace to a distinct genetic group, with the exception of Mahmoudi and Biskri landraces that were both divided into two groups. This result is in agreement with our previous study [22], where 8 out of 11 Tunisian durum wheat landraces corresponded to distinct genetic groups using 10 SSR markers. AMOVA analysis showed a high genetic variability (61%) between subpopulations, suggesting that the three subpopulations were derived from different gene pools. This is consistent with previous studies of tetraploid wheats [2, 68] and barley [69] showing that local landraces were derived from multiple ancestral populations and had reticulated phylogenetic relationships. However, several other durum wheat landrace studies reported a higher or an equal genetic variation whitin population than among population [30, 70, 71]. Soriano et al. [11] detected only 13% of genetic varibility between four genetic subpopulations that were tracing distinct geographical pattern of the Mediterranean durum wheat germplasm. In this study, GS1 and GS3 were composed essentially with Mahmoudi and Biada landrace, respectively. These landraces have preferred specific agro-morphological traits recognized by the farmers, selected and multiplied over time. The formation of distinct domesticated gene pools were also reported for several other crops such as common bean, which underwent parallel evolution and spread further through the development of landraces with distinct characteristics and specific adaptations [72,73,74,75].

Informations about putative origins and years of introduction in Tunisia of the 13 herein studied landraces, along with their agro-morphological characteristics as described by Ouaja et al. [22], were analyzed to explain their inferred genetic structure. In fact, the grouping pattern of accessions appear to be associated, to some extent, with the geographical pattern of the landraces. Indeed, Slim et al. [20] highlighted a strong North-South stratification in Tunisian durum wheat landraces with the prevalence of modern cultivars in the North versus landraces, grown by small-holder farmers under low-input traditional agrosystems in the marginal areas of the Center and South However, other durum wheat landrace studies showed a genetic clustering irrespective of their geographical origin, underlying the presence of plant material exchanges that could have reduced the genetic differentiation [41]. In this study, the subpopulation GS2 mainly included landraces of North African and East Mediterranean origins. GS2 includes landrace Mekki from Morocco, landrace Taganrog from Southern Russia, landraces Azizi, Jneh Khotifa and Sbei glabre that were considered as local populations and Richi which was reported as foreign [17, 33, 76]. Robbana et al. [13] also reported that Tunisian landraces were genetically associated to North African landraces. However, using a set of DArtSeq markers to describe the genetic diversity of Tunisian landraces, Robbana et al. [13] reported that Jneh Zarzoura, a close relative to Jneh Khotifa [17] clustered distinctly with accessions from Jordan. Moreover, local landraces names have been traditionally selected by natives, generally according to morphological features or locality, and are often consciously used by farmers for management, selection or exchanges [14, 77]. Sahri et al. [77] particularly highlighted the significance of variety name, which have largely influenced the structure and evolutionary dynamics of durum wheat in Morocco. Therefore, nominal analogies were frequently reported between landraces of different Mediterranean regions, probably reflecting trade and migration of the same landraces around the Mediterranean basin. Xynias et al. [6] reported that the Italian cultivar « Senatore Capelli » was selected in 1915 from the local North African landrace « Jean Retifah » which was very prominent and marked the cultivation of durum wheat in Italy. The landrace « Jean Retifah » must be the known Jneh Khotifa in Tunisia [14, 17]. In addition, Soriano et al. [11] studied the structure of durum wheat landraces from 21 Mediterranean countries, using SSR markers, and reported two Italian landraces named « Hymera » and « Aziziah » that were associated to eastern Mediterranean genetic group. Knowing that Tunisia was the former bread basket of the Roman Empire [78, 79], these latters could be the landraces known Hmira and Azizi in Tunisia [14, 17, 20, 22], suggesting though that Azizi landrace may also have an eastern origin. A recent study on 170 durum wheat landraces from 24 Mediterranean countries revealed that more than 23% of Tunisian landraces were from eastern mediterraneen countries [40]. Boeuf [14] mentioned that the landraces of North Africa had dominant characters, specific to Abyssinian wheats, such as red, purple or black spikes, pubescent glumes and red or dark-colored grains. These characters were totally absent and unknown in Europe. Herein it was noted that landraces of GS2, with the exception of Biskri, shared features of spikes and grains [22] corresponding to the Abyssinian wheats as Bœuf [14] reported.

Furthermore, according to Bœuf [14], the geographical expansion and domestication of wheats from Abyssinia enhanced the accumulation of recessive characters, in particular white spikes, hairless glumes and light-colored grains widely adopted in Europe. During crop domestication process, several changes were induced for major morphological, structural and functional traits associated with adaptation and cultivation in order to meet human needs, as reported for the common bean ) [72]. Bœuf [14] mentioned that the whiteness of the spike and light-colored grains were among the most sought after and preferred criteria in wheat by European farmers during commercial trade in North Africa. So far, these traits are characteristics of the landraces Mahmoudi (GS1), Biskri (GS2), Biada and Bidi (GS3) [22]. Landraces Biskri, Biada and Bidi were introduced in Tunisia from Algeria, while Mahmoudi was considered as a local landrace population with various reported origins including Algeria and Italy. Another example reflecting nominal analogy between Mediterranean landraces is the Tunisian landrace Biada and the Spanish landraces « Blancal » and « Blanco de Baleares » meaning, among others, the white wheats [17, 22, 80], thus indicating that such phenotypic characteristic was probably derived from western Mediterranean. Based on these findings/statements, although the majority of Tunisian landraces are North Africa, we suggest that subpopulations GS1 and GS3 were introduced to North Africa and particularly to Tunisia from Europe/western Mediterranean, unlike subpopulation GS2 which may be originated from the Middle East. Ben Krima et al. [30] also agree that the combination of both genetic and agro-morphological approches are essential for retracing the history, origin and dynamic lifestory of Tunisian durum wheat landraces. In the same context, Moragues et al. [9] highlighted two dispersal pathways of the Mediterranean durum wheat landraces which had contributed to the divergence of these landraces into distinct genetic groups following their adaptation to different local environments. A first pathway through the North-East of the Mediterranean basin to Europe, and a second pathway through the South of the Mediterranean basin to North Africa reaching the Iberian Peninsula.

Genetic distribution of the resistance of Z. tritici in the tunisian durum wheat populations

Tunisian durum wheat landraces have been reported to exhibit valuable sources of resistance to STB, useful to include in breeding programs and to develop varieties with durable and broad spectrum of resistance [21, 63, 64]. In the present study, the 366 genotyped durum wheat accessions were also screened for Z. tritici resistance, under field conditions, using two Z. tritici isolates Tun06 and TM220 collected from two Tunisian regions, Bizerte and Manouba, respectively. The analysis of variance showed a non-significant variation between isolates Tun06 and TM220 towards the durum wheat accessions.This result agrees with Ferjaoui et al. [63] findings, detecting only two virulence profils among 55 Tunisian Z. tritici isolates screened at seedling stage.

Overall, 60% of the accessions were susceptible, showing that Tun06 and TM220 isolates were virulent on the majority of the accessions, which reflect the adaptation of Z. tritici virulence to durum wheat landraces in Tunisia. Likewise, Ouaja et al. [21] suggested that the Tunisian Z. tritici isolate Tun06 still preserve virulences against old durum wheat landraces even though they are currently marginally grown in wheat production areas in Tunisia. In addition, several studies reported that Z. tritici undergoes frequent sexual reproduction on durum wheat in Tunisia [81, 82]. In fact, sexual reproduction plays a key role in the evolution of pathogenicity traits, including virulence and aggressiveness [83], enabling the fungus to quickly evolve and circumvent the resistance genes by creating new combinations of alleles and, in combination with the asexual reproduction allowing frequent generation of the new genotypes [49, 82]. Nevertheless, about 27% of the collection was resistant, suggesting that Tunisian durum wheat landraces still carry effective STB resistance genes.

The analysis of variance showed a large and significant variation between GS1, GS2 and GS3 regarding the resistance to isolate Tun06. Subpopulations GS1 (composed of Mahmoudi accessions) and GS3 (Biada accessions) showed higher frequencies of susceptible responses than GS2. The distribution of the resistance within each landrace observed in the UPGMA tree, indicate that the resistance relies on the landrace instead of the genetic structure (K = 3), as GS2 grouped both resistant and susceptible genotypes. These results also indicate that although landraces of GS2 formed an invidivualized genetic group, they may harbor combination of resistance genes that differ in nature, number, structure, chromosomal localisations and type of interaction that still need to be depicted and revealed by a genome wide association study (GWAS).

The susceptibility observed in the GS1, GS3, Biskri, and Azizi of GS2 genotypes can be attributed to their widespread use and commercial share, leading to extensive cultivation, particularly in northern Tunisia. Over time, these landraces have lost their resistance to the disease [14, 17]. In fact, Z. tritici, the causative agent of the disease, is predominantly prevalent in the northern and northwestern regions of Tunisia where sub-humid zones become significant hotspots for STB, exhibiting high levels of infection [59, 60, 84]. As a result of the rapid adaptation of Z. tritici isolates, the landraces extensively grown in the northern zones became susceptible to the disease [49, 56, 61, 85]. Moreover, rapid adaptation of Z. tritici to landraces from western mediterannen origin (GS1 and GS3) could be facilitated by domestication process mainly accompanied by a strong reduction in genetic diversity and/or high levels of gene loss compared to wild ancestors or wild gene pool [86,87,88]. These events reduced the adaptation of cultivated wheat to erratic environmental variations, where wild traits show much greater fitness over domesticated ones [75]. This senario might be consistent with the significant variability of the resistance to Z. tritici observed between subpopulations in the present study, suggesting that resistance (R) genes of the three subpopulations may have evolved divergently and crucially under a combination of environmental and human pressures. Therefore, we hypothesize that major genes conferring specific resistance (as dominant characters) have undergone modifications and alterations by mutations during the geographic expansion and with the wide exchange network of durum wheat landraces among Mediterranean regions, which resulted in the loss of dominant R genes/alleles and the spread of susceptibility as observed within landraces of subpopulations GS1 and GS3. Alternatively, the results suggest that landraces of GS1 and GS3 were initially susceptible to Z. tritici when introduced to Tunisia. Nevertheless, unlike GS1 and GS3, most landraces of GS2 (except for Biskri and Azizi landraces) were resistant because they had recourse to different geographical patterns/pathways as they were local or directly derived from an eastern origin and did not pass via Europe. Western Europe farmers were among the first to create and adopt modern methods of plant breeding and exert wheat genetic improvement, involving direct selection for homogeneous material, thus, indirectly reducing the variability of the genetic sources of resistance or even unintentionally selecting for recessive genes [14, 89, 90]. In contrast, smallholder farmers in North Africa have been preserving the local durum wheat diversity with on-farm conservation practices over generations [14, 17, 22, 77].

In this study, 22% of the accessions were admixed, among which 51% of admixed genotypes between GS2 and GS1 and 24% admixed between GS2 and GS3. Admixtures occurs mainly by gene flows, through the frequent introduction of new genotypes into fields and seed exchange network within and between farmer communities [2, 11, 16]. In fact, Nm between GS1 and GS2 and between GS2 and GS3 were both around 0.4. All the resistant accessions of Mahmoudi were admixed between GS1 and GS2 and all the resistant accessions of Biada, Souri and Roussia were admixed between GS3 and GS2; suggesting that resistant landraces in GS2 were probably the sources transmitting resistance to Z. tritici via admixture. Indeed, resistant Mahmoudi genotypes may have acquired their resistance from Jneh Khotifa, as they appeared phylogenetically close in the UPGMA tree. The UPGMA tree also showed that modern varieties are genetically close to Biada (GS3), indicating that modern varieties have been selected for certain agro-morphological and phenological characteristics of GS3, such as white and short spikes, short plant and precocity. In addition, modern varieties with a genome derived from GS3 or GS1 or admixed, have been selected from a susceptible background to Z. tritici. It is therefore necessary to re-direct the breeding programs towards another genepool presenting resistance to Z. tritici, for example by developping marker-assisted selection to introgress the resistance to Z. tritici from GS2. On the other hand, admixture may have also caused loss of resistance as some admixed Jneh Khotifa and Richi accessions were susceptible to Tun06. This result indicated that admixture between genetically distinct landraces/populations and frequently recurring gene exchanges (or gametic association between gene loci) [91], may have elicited susceptibility within initially resistant accessions throughout an alteration or loss of the resistance genes/alleles.


The present study revealed that Tunisian durum wheat life history of resistance to STB involve the interaction of miscellaneous factors, including the landrace genetic structure and introductory pathways, the local commercial share defining the geographic and regional distribution of the landraces and the occurrence of admixtures within these landraces. In fact, landraces were subjected to genetic differentiation in time and space during their introductory pathways in the Mediterranean area, their adaptation to specific environments and the Human selection pressure/domestication, contributing to their divergence in their resistance spectra. In addition, the regional distribution of the landraces across Tunisia influenced the evolutionary history of resistance genes under different climatic constraints in relation to hotspots of STB disease. Admixtures also represent one of the main driver factors of the resistance to Z. tritici involving old or ongoing exchanges that contributed to the introgression and/or the loss of resistance genes. Our results suggest the presence of diverse resistance sources towards two Z. tritici isolates Tun06 and TM220 in the Tunisian durum wheat landraces, which implies the need for more in-depth research to investigate and characterize the related resistance genes. Continuous identification of new sources of resistance to STB is required for the development of wheat cultivars with sustainable field resistance. Therefore, old local durum wheat landraces can play an important role as donor of resistance genes in breeding programs. In addition, further investigation and study of virulence patterns in Z. tritici populations in Tunisia will provide new insights and understanding of the Z. tritici-durum wheat interactions and trace their co-evolution.


Wheat accessions

A set of 375 durum wheat accessions was used in this study. Accessions were collected by the National Gene Bank (BNG) from four regions in Central (the Sahel and Kairoun) and southern (Gabes and Medenine) Tunisia (Table S3). Accessions were morphologically characterized, identified and classified into 13 landrace-populations namely Mahmoudi, Biada, Bidi, Biskri, Azizi, Badri, Mekki, Jneh Khotifa, Sbei glabre, Taganrog, Richi, Souri and Roussia as described by Ouaja et al. [22]. Information about the origin and year of introduction and selection of these landraces are depicted in Table S4.

Field trials and Z. tritici isolates

Durum wheat landraces were screened for their resistance to Z. tritici at the adult plant stage in the experimental station of CRRGC Beja in northwest Tunisia. Inoculation assays were performed using two durum wheat-derived Z. tritici isolates on separate field experimental plots; notably the well-characterized and virulent reference isolate Tun06 (also reported as TunBz-1) sampled from the Bizerte region of Tunisia in 2006 and the isolate TM220 collected from the Chili landrace cultivated in Lansarine region of Tunisia in 2014 [21, 61, 92, 93]. The virulence patterns of TM220 and Tun06 isolates were previously assessed at seedling stage on 21 old durum wheat accessions set and revealed that TM220 was virulent on 8 accessions comparing to isolate Tun06 which was avirulent (Thierry Marcel personnel com.). Tun06- and TM220 field trials were realized as reported by Ouaja et al. [21] and followed an Augmented Randomized Complete Block Design (ARCBD), including 6 blocks spaced apart of 1 m. Blocks were 1 m width linearly drilled and accessions were sown 20 cm spaced apart. Each block was composed of 70 accessions and 7 checks. The checks included six susceptible to moderately resistant modern durum wheat cultivars “Karim, Khiar, Om Rabia, Salim, Maali and Nasr” [60, 61]. Durum wheat landraces were evaluated for resistance to Tun06 isolate during two consecutive growing seasons (2015–2016 and 2016–2017), and for resistance to isolate TM220 on the growing season of 2016–2017.

Inoculation and screening for resistance to Z. tritici

Tun06 and TM220 isolates were maintained frozen at -80 °C. A preculture of the inoculum was prepared by shaking Z. tritici isolates (at 100 rpm/ 25 °C) for six to seven days, in 100 ml yeast glucose liquid medium (30 g glucose, 10 g yeast per liter demineralized water). The produced spore suspensions of Tun06 and TM220 were subsequently transferred to 500 ml yeast glucose liquid media and were incubated under the aforementioned conditions to provide sufficient inoculum for the field trials. Spores of both isolates were collected after overnight settling in static cultures, concentrated by decanting the supernatant medium, and were adjusted to 106 spore/ml. Accessions in all experimental plots were inoculated twice, at the three-leaf stage (approximately Growing stage 21) and at the stem elongation stage (approximately Growing stage 37) [94] as described by Ouaja et al. [21].

The disease severity was evaluated by estimating pycnidia coveragepercentages which correspond to sporulating area. The same leaf layer (F3-F4) of each accession in all experimental plots was assessed for STB resistance, at three time points for isolate Tun06 and at two time points for isolate TM220. Tun06-disease scores were assessed at 18, 38 and 46 days post the second inoculation (dpi) during 2015–2016 growing season, and at 15, 35 and 53 dpi during 2016–2017 growing season. TM220-disease scores were assessed at 20 and 40 dpi during 2016–2017 growing season. Screening data were used to calculate the Area Under the Disease Progress Curve (AUDPC) and the Relative Area Under the Disease Progress Curve (RAUDPC) for quantitative analyses of the temporal differences in disease progress as detailed by Ouaja et al. [21].

Three classes of response to Z. tritici at the adult plant stage, defined by Ouaja et al. [21], were considered in this study. The resistant class (R) include accessions having pycnidia score < 25%, while the susceptible class (S) include accessions showing > 47% of pycnidia on infected leaves. Accessions with pycnidia score ranging between 25% and 46% constitute the class of intermediate (I). Two and 63 accessions were not scored for STB resistancefor Tun06 and TM220, respectively. The variation of classes of response R, I and S between the isolates Tun06 and TM220 is shown in Table S7.

Genotyping and data analysis

A total of 20,279 polymorphic SNPs were generated by Illumina sequencing 375 Tunisian durum wheat accessions using a High-density 90 K wheat SNP array (iSelect, San Diego, USA) [67], among which 13,998 (~ 70%) were mapped in the consensus durum wheat genetic map [95]. The filtered SNPs had < 5% of missing data and a minor allele frequency of 5% were filtered. From these markers, 286 SNPs covering the entire genome with 12 to 23 SNPs per chromosome of 10 cM apart and with a PIC > 0.3, were selected to perform the analysis in the present study (Table S5). PIC values were calculated by determining the frequency of alleles per locus as Powell et al. [96] (Table S6). Genotypes with > 10% of missing data out of the 286 SNPs were eliminated, and 366 genotypes were included in the rest of the analysis. Based on 286 SNP data generated for 366 accessions, 147 multilocus genotypes (MLG) were identified with R 3.3.2 [97] (Table S1). According to genotype accumulation curves (Figure S2), performed under R 3.3.2 [97], the 286 SNPs (Figure S2a) provide enough discrimination between the indiviuduals than the 13,998 SNPs (Figure S2b) and have a good power of resolution to study the population structure and genetic diversity of the Tunisian durum wheat.

The genotypic data based on the 286 selected SNP markers were used to study the population structure of the Tunisian durum wheat accessions with the model-based clustering algorithm STRUCTURE 2.3.4 [65]. STRUCTURE program was run on 147 MLG with K values between 1 and 20, by applying 10 independent runs for each value of K, 100,000 burnins and 100,000 Markov Chain Monte Carlo (MCMC) repetitions. The optimal number of genetic subpopulations was determined using the mean posterior probability (ln P(D)) value per cluster (K) and the delta-K method of ln P(D) under STRUCTURE harvester [66]. Population structure was investigated at the optimal K as well as at K ranging from 13 to 20 in order to assertain the genetic grouping according to landrace denomination. Individuals presenting a membership coeficient less than 0.7 to any genetic subpopulation were considered as admixed genotypes. To study the phylogenetic relationships between accessions/landraces, an Unweighted Pair Group Method with Arithmetic Average (UPGMA) tree was created using Nei standard genetic distance [98]. A bootstrap analysis was performed with the program Populations 1.2.32 [99] and branch support values were estimated using 1000 bootstrap randomizations. A consensus tree was generated and viewed using TreeView 1.6.6 [100].

In addition, pairwise Fst and gene flow (Nm) coefficients were calculated with GenAlEx 6.501 [101] by genetic subpopulations and by landraces (based on their nomenclature). An analysis of molecular variance (AMOVA) were also carried out using GenAlEx 6.501 [101] in order to investigate the significance of genetic subpopulations differentiation as well as genetic differentiation between and among landraces. A Principal Components Analysis (PCAs) on the genotypic data was performed under GenAlEx 6.501 [101]. This multivariate analysis was completed to visualize the structure of the 366 durum wheat accessions sorted by genetic subpopulations (as defined by STRUCTURE), by landrace denomination [22] and by STB resistance class to Tun06 (R, I and S). Furthermore, genetic diversity indices (Ne, I, He, P) were calculated by genetic subpopulation and by landrace [101].

An analyses of variance (ANOVA) were carried out under R 3.3.2 [97] to determine the effect of genetic subpopulations and isolates on RAUDPC scores. Genetic subpopulations x isolates interaction were also investigated.

Data Availability

The seeds of each accession can be requested from the National Gene Bank of Tunisia GRIN-Global repository ( according to their NGBTUN ID under the following web link The data sets supporting the results of this article are included in this manuscript in its additional information files. Genotyping data are available in an excel sheet ‘Table S5’ with a corresponding SNP code.



National Gene Bank of Tunisia


Septoria tritici blotch


Single Nucloetoide Polymorphic marker


Unweighted Pair Group Method with Arithmetic Average


Genetic subcluster; GS = genetic subpopulation


  1. Mackey J. Wheat: its concept, evolution and taxonomy In: Royo C, Nachit MM, Di Fonzo N, Araus JL, Pfeiffer WH, Slafer GA, editors. Durum wheat breeding: current approaches and future strategies. In.: New York: Food Products Press; 2005.

  2. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Kabbaj H, Sall AT, Al-Abdallat A, Geleta M, Amri A, Filali-Maltouf A, Belkadi B, Ortiz R, Bassi FM. 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).

  4. Baloch FS, Alsaleh A, Shahid MQ, Çiftçi V, de Miera LES, Aasim M, Nadeem MA, Aktaş H, Özkan H, Hatipoğlu R. A whole genome DArTseq and SNP analysis for genetic diversity assessment in durum wheat from central fertile crescent. PLoS ONE. 2017;12(1):e0167821.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Martínez-Moreno F, Solís I, Noguero D, Blanco A, Özberk İ, Nsarellah N, Elias E, Mylonas I, Soriano JM. Durum wheat in the Mediterranean Rim: historical evolution and genetic resources. Genetic Resources and Crop Evolution; 2020.

  6. Xynias IN, Mylonas I, Korpetis EG, Ninou E, Tsaballa A, Avdikos ID, Mavromatis AG. Durum wheat breeding in the Mediterranean region: current status and future prospects. Agronomy. 2020;10(3):432.

    Article  Google Scholar 

  7. De Ron AM, Bebeli PJ, Negri V, Vaz Patto MC, Revilla P. Warm season grain legume landraces from the south of Europe for germplasm conservation and genetic improvement. Front Plant Sci. 2018;9:1524.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Balfourier F, Bouchet S, Robert S, De Oliveira R, Rimbert H, Kitt J, Choulet F, International Wheat Genome Sequençing Consortium, BreadWheat Consortium, Paux E. Worldwide phylogeography and history of wheat genetic diversity. Sci Adv. 2019;5:eaav0536.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Moragues M, Moralejo M, Sorrells ME, 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(5):1133–44.

    Article  CAS  Google Scholar 

  10. Haudry A, Cenci A, Ravel C, Bataillon T, Brunel D, Poncet C, Hochu I, Poirier S, Santoni S, Glémin S et al. Grinding up wheat: a massive loss of nucleotide diversity since domestication. Mol Biol Evol 2007, 24.

  11. Soriano JM, Villegas D, Aranzana MJ, García del Moral LF, Royo C. Genetic structure of modern Durum Wheat Cultivars and Mediterranean Landraces Matches with their agronomic performance. PLoS ONE. 2016;11(8):e0160983.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Essid MY. History of Mediterranean food. MediTERRA 2012 (english). Presses de Sciences Po; 2012: 51–69.

  13. Robbana C, Kehel Z, Naceur B, Sansaloni C, Bassi F, Amri A. Genome-wide genetic diversity and population structure of tunisian durum wheat landraces based on DArTseq technology. Int J Mol Sci. 2019;20(6):1352.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Boeuf F. Le Blé en Tunisie: La plante. Le milieu physico-chimique: Société anonyme de l’impr. rapide de Tunis; 1932.

  15. Vavilov NI. The origin, variation, immunity and breeding of cultivated plants. Volume 72. LWW; 1951.

  16. Ren J, Sun D, Chen L, You FM, 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(4):7061–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Deghaïs M, Kouki M, Gharbi M, El Felah M. Les variétés de céréales cultivées en Tunisie. Ministère de L’Agriculture et des Ressources Hydrauliques Tunis, Tunisie 2007.

  18. Sourour A, Chahine K, Youssef T, Olfa S-A, Hajer S-A. Phenotypic diversity of tunisian durum wheat landraces. Afr Crop Sci J 2010, 18(1).

  19. Amine S, Sourour A, Olfa S-A, Cyrine R, da Silva JAT, Hajer S-A. Morphological diversity of some qualitative traits in tetraploid wheat landrace populations collected in the south of Tunisia. Int J Plant Breed. 2011;5(1):67–70.

    Google Scholar 

  20. Slim A, Piarulli L, Chennaoui Kourda H, Rouaissi M, Robbana C, Chaabane R, Pignone D, Montemurro C, Mangini G. Genetic structure analysis of a collection of tunisian durum wheat germplasm. Int J Mol Sci. 2019;20(13):3362.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Ouaja M, Aouini L, Bahri B, Ferjaoui S, Medini M, Marcel TC, Hamza S. Identification of valuable sources of resistance to Zymoseptoria tritici in the tunisian durum wheat landraces. Eur J Plant Pathol. 2020;156(2):647–61.

    Article  CAS  Google Scholar 

  22. Ouaja M, Bahri BA, Aouini L, Ferjaoui S, Medini M, Marcel TC, Hamza S. Morphological characterization and genetic diversity analysis of tunisian durum wheat (Triticum turgidum var. Durum) accessions. BMC Genomic Data. 2021;22(1):1–17.

    Article  Google Scholar 

  23. Ayadi S, Karmous C, Hammami Z, Tamani N, Trifa Y, Esposito S, Rezgui S. Genetic variability of Nitrogen Use Efficiency components in tunisian improved genotypes and landraces of durum wheat. Agric Sci Res J. 2012;2(11):591–601.

    Google Scholar 

  24. Chamekh Z, Karmous C, Ayadi S, Sahli A, Hammami Z, Fraj MB, Benaissa N, Trifa Y, Slim-Amara H. Stability analysis of yield component traits in 25 durum wheat (Triticum durum Desf.) Genotypes under contrasting irrigation water salinity. Agric Water Manage. 2015;152:1–6.

    Article  Google Scholar 

  25. Babay E, Mnasri SR, Mzid R, Hanana M. Quality selection and genetic diversity of tunisian durum wheat varieties using SSR markers. Bioscience J 2019, 35(4).

  26. Bouacha OD, Rezgui S. Spaghetti quality: comparison between landraces and high yielding tunisian durum wheat varieties. J New Sci. 2019;64(7):4056–60.

    Google Scholar 

  27. Yacoubi I, Nigro D, Sayar R, Masmoudi K, Seo YW, Brini F, Giove S, Mangini G, Giancaspro A, Marcotuli I. New insight into the North-African durum wheat biodiversity: phenotypic variations for adaptive and agronomic traits. Genet Resour Crop Evol. 2020;67(2):445–55.

    Article  CAS  Google Scholar 

  28. Huhn MR. Tetraploid tunisian wheat germplasm as a new source of Fusarium head blight resistance. North Dakota State University; 2008.

  29. Ferjaoui S, Sbei A, Aouadi N, Hamza S. Monogenic inheritance of resistance to septoria tritici blotch in durum wheat ‘Agili’. Int J Plant Breed. 2011;5:17–20.

    Google Scholar 

  30. Ben Krima SB, Slim A, Gelisse S, Kouki H, Nadaud I, Sourdille P, Yahyoui A, M’Barek SB, Suffert F, Marcel TC. Life story of Tunisian durum wheat landraces revealed by their genetic and phenotypic diversity. bioRxiv 2020.

  31. 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, farmers’ taxonomy and tetraploid wheat domestication history. BMC Evol Biol. 2014;14(1):1–18.

    Article  Google Scholar 

  32. Ortiz R, Trethowan R, Ferrara GO, Iwanaga M, Dodds JH, Crouch JH, Crossa J, Braun H-J. High yield potential, shuttle breeding, genetic diversity, and a new international wheat improvement strategy. Euphytica. 2007;157(3):365–84.

    Article  Google Scholar 

  33. Bonjean AP, Angus WJ. The World Wheat Book: a history of wheat breeding. Lavoisier Publishing; 2001.

  34. Ayed S, Slim AH. Distribution and phenotypic variability aspects of some quantitiative traits among durum wheat accessions. Afr Crop Sci Soc. 2008;16(4):219–24.

    Google Scholar 

  35. Ayed S, Karmous C, Trifa Y, Slama-Ayed O, Slim-Amara H. Phenotypic diversity of tunisian durum wheat landraces. Afr Crop Sci Soc. 2010;18(1):35–42.

    Article  Google Scholar 

  36. Slim A, Ayed S, Slama-Ayed O, Robbana C, Jaime A, da Silva T, Slim-Amara H. Morphological diversity of some qualitative traits in tetraploid wheat landrace populations collected in the South of Tunisia. Int J Plant Breed. 2011;5(1):67–70.

    Google Scholar 

  37. Babay E, MnasriI S, Khamassi K, Ouriri I, Hanana M. Morphological and end-use quality characterization of tunisian durum wheat. Agrobiologia. 2019;9(2):1560–7.

    Google Scholar 

  38. 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 Ev. 2005;52:21–31.

    Article  CAS  Google Scholar 

  39. Miazzi MM, Babay E, De Vita P, Montemurro C, Chaabane R, Taranto F, Mangini G. Comparative genetic analysis of Durum Wheat Landraces and Cultivars widespread in Tunisia. Front Plant Sci. 2022;13:939609. PMID: 35909756; PMCID: PMC9326505.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Rufo R, Alvaro F, Royo C, Soriano JM. From landraces to improved cultivars: assessment of genetic diversity and population structure of Mediterranean wheat using SNP markers. PLoS ONE. 2019;14(7):e0219867.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Alemu A, Feyissa T, Letta T, Abeyo B. Genetic diversity and population structure analysis based on the high-density SNP markers in Ethiopian durum wheat (Triticum turgidum ssp. durum):. BMC Genetics 2020, 21: 18. DOI: 1,0.1186/s12863-020-0825.

  42. O’Driscoll A, Kildea S, Doohan F, Spink J, Mullins E. The wheat–Septoria conflict: a new front opening up? Trends Plant Sci. 2014;19(9):602–10.

    Article  PubMed  Google Scholar 

  43. Dean R, Van Kan JAL, Pretorius ZA, Hammond-Kosack KE, Di Pietro A, Spanu PD, Rudd JJ, Dickman M, Kahmann R, Ellis J, et al. The top 10 fungal pathogens in molecular plant pathology. Mol Plant Pathol. 2012;13(4):414–30.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Stukenbrock EH, Banke S, Javan-Nikkhah M, McDonald BA. Origin and domestication of the fungal wheat pathogen Mycosphaerella graminicola via sympatric speciation. Mol Biol Evol. 2007;24(2):398–411.

    Article  CAS  PubMed  Google Scholar 

  45. Stukenbrock EH, McDonald BA. The origins of plant pathogens in agro-ecosystems. Annu Rev Phytopathol. 2008;46:75–100.

    Article  CAS  PubMed  Google Scholar 

  46. Torriani SFF, Stukenbrock EH, Brunner PC, McDonald BA, Croll D. Evidence for extensive recent intron transposition in closely related fungi. 2011, 21:2017–22.

  47. McDonald MC, McDonald BA, Solomon PS. Recent advances in the Zymoseptoria tritici–wheat interaction: insights from pathogenomics. Front Plant Sci 2015, 6.

  48. Stukenbrock EH, Bataillon T, Dutheil JY, Hansen TT, Li R, Zala M. The making of a new pathogen: insights from comparative population genomics of the domesticated wheat pathogen Mycosphaerella graminicola and its wild sister species. Genome Res 2011, 21.

  49. McDonald BA, Mundt CC. How knowledge of pathogen population biology informs management of Septoria tritici blotch. Phytopathology. 2016;106(9):948–55.

    Article  PubMed  Google Scholar 

  50. Kema G, Annone JG, Sayoud R, Van Silfhout CH, Van Ginkel M, De Bree J. Genetic variation for virulence and resistance in the wheat-Mycosphaerella graminicola pathosystem I. interactions between pathogen isolates and host cultivars. Phytopathology. 1996;86(2):200–12.

    Article  Google Scholar 

  51. Kema G, Sayoud R, Annone J, Van Silfhout C. Genetic variation for virulence and resistance in the wheat-Mycosphaerella graminicola pathosystem. II: analysis of interactions between pathogen isolates and host cultivars. Phytopathology. 1996;86(2):213–20.

    Article  Google Scholar 

  52. Kema GH, van Silfhout CH. Genetic variation for virulence and resistance in the wheat-Mycosphaerella graminicola pathosystem III. Comparative seedling and adult plant experiments. Phytopathology. 1997;87(3):266–72.

    Article  CAS  PubMed  Google Scholar 

  53. Ghaffary SMT, Robert O, Laurent V, Lonnet P, Margalé E, van der Lee TA, Visser RG, Kema GH. Genetic analysis of resistance to septoria tritici blotch in the french winter wheat cultivars Balance and Apache. Theor Appl Genet. 2011;123(5):741–54.

    Article  PubMed  Google Scholar 

  54. Brown JK, Chartrain L, Lasserre-Zuber P, Saintenac C. Genetics of resistance to Zymoseptoria tritici and applications to wheat breeding. Fungal Genet Biol. 2015;79:33–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Arraiano LS, Brown JK. Sources of resistance and susceptibility to Septoria tritici blotch of wheat. Mol Plant Pathol 2016.

  56. Zhong Z, Marcel TC, Hartmann FE, Ma X, Plissonneau C, Zala M, Palma-Guerrero J. A small secreted protein in Zymoseptoria tritici is responsible for avirulence on wheat cultivars carrying the Stb6 resistance gene. New Phytol 2017, 214.

  57. Kema GH, Verstappen EC, Waalwijk C. Avirulence in the wheat Septoria tritici leaf blotch fungus Mycosphaerella graminicola is controlled by a single locus. Mol Plant Microbe Interact. 2000;13(12):1375–9.

    Article  CAS  PubMed  Google Scholar 

  58. Kema GH, Gohari AM, Aouini L, Gibriel HA, Ware SB, van Den Bosch F, Manning-Smith R, Alonso-Chavez V, Helps J, M’Barek SB. Stress and sexual reproduction affect the dynamics of the wheat pathogen effector AvrStb6 and strobilurin resistance. Nat Genet. 2018;50(3):375.

    Article  CAS  PubMed  Google Scholar 

  59. Fakhfakh M, Rezgui S, Jarrahi T, Zouani R, Bousalmi A, Sgatni K, Yahyaoui A. Mapping distribution of Septoria allowed targeting management options of the disease. In: 8th International Symposium on Mycosphaerella and Stagonospora Diseases of Cereals, Mexico City: 2011; 2011: 10–14.

  60. Berraies S, Ammar K, Salah Gharbi M, Yahyaoui A, Rezgui S. Quantitative inheritance of resistance to Septoria tritici blotch in durum wheat in Tunisia. Chil J agricultural Res. 2014;74(1):35–40.

    Article  Google Scholar 

  61. Gharbi M, Deghais M, Ben Amar F. Breeding for resistance to Septoria tritici in durum wheat. In: Proceedings of Durum Wheat Conference, Zaragoza 2000: 2000; 2000: 397–401.

  62. Gharbi M, El Felah M. Article de synthèse Les céréales en tunisie: plus d’un siècle de recherche variétale. Ann DE L’INRATT. 2013;86:45.

    Google Scholar 

  63. Ferjaoui S, M’Barek S, Bahri B, Slimane R, Hamza S. Identification of resistance sources to septoria tritici blotch in old tunisian durum wheat germplasm appliad for the analysis of the Zymoseptoria tritici-durum wheat interaction. J Plant Pathol. 2015;97(3):471–81.

    Google Scholar 

  64. Aouini L. Durum wheat and septoria tritici blotch: genes and prospects for breeding. Wageningen University; 2018.

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

  66. Earl DA. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4(2):359–61.

    Article  Google Scholar 

  67. Wang S, Wong D, Forrest K, Allen A, Chao S, Huang BE, Maccaferri M, Salvi S, Milner SG, Cattivelli L. Characterization of polyploid wheat genomic diversity using a high-density 90 000 single nucleotide polymorphism array. Plant Biotechnol J. 2014;12(6):787–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Civáň P, Ivaničová Z, Brown TA. Reticulated origin of domesticated emmer wheat supports a dynamic model for the emergence of agriculture in the fertile crescent. PLoS ONE. 2013;8(11):e81955.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Poets AM, Fang Z, Clegg MT, Morrell PL. Barley landraces are characterized by geographically heterogeneous genomic origins. Genome Biol. 2015;16(1):1–11.

    Article  Google Scholar 

  70. Asmamaw M, Keneni G, Kassahun T. Genetic diversity of Ethiopian durum wheat (Triticum durum Desf.) landrace collections as revealed by SSR Markers: Advances in Crop Science and Technology, 2019, 7(1): 903 1000413. DOI:

  71. Kyratzis A, Nikoloudakis N, Katsiotis A. Genetic variability in landraces populations and the risk to lose genetic variation. The example of landrace ‘Kyperounda’ and its implications for ex situ 1040 conservation. PLoS ONE. 2019;14(10):e0224255.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Papa R, Gepts P. Asymmetry of gene flow and differential geographical structure of molecular diversity in wild and domesticated common bean (Phaseolus vulgaris L.) from Mesoamerica. Theor Appl Genet. 2003;106(2):239–50.

    Article  CAS  PubMed  Google Scholar 

  73. Bitocchi E, Rau D, Benazzo A, Bellucci E, Goretti D, Biagetti E, Panziera A, Laido G, Rodriguez M, Gioia T. High level of nonsynonymous changes in common bean suggests that selection under domestication increased functional diversity at target traits. Front Plant Sci. 2017;7:2005.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Bitocchi E, Rau D, Bellucci E, Rodriguez M, Murgia ML, Gioia T, Santo D, Nanni L, Attene G, Papa R. Beans (Phaseolus ssp.) as a model for understanding crop evolution. Front Plant Sci. 2017;8:722.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Di Vittori V, Bellucci E, Bitocchi E, Rau D, Rodriguez M, Murgia ML, Nanni L, Attene G, Papa R. Domestication and crop history. The common bean genome. Springer; 2017: 21–55.

  76. Deghais M, El Faleh M, Gharbi M, Zarcouna T, Chakroun M. Les acquis de l’amélioration génétique des céréales en Tunisie. In: Annales de l’INRAT: 1999; 1999: 3–30.

  77. 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, farmers’ taxonomy and tetraploid wheat domestication history. BMC Evol Biol. 2014;14(1):264.

    Article  PubMed  PubMed Central  Google Scholar 

  78. Fabricant F. Rome’s glory is now Tunisia’s. New York Times; 1998.

  79. Latiri K, Lhomme J-P, Annabi M, Setter TL. Wheat production in Tunisia: progress, inter-annual variability and relation to rainfall. Eur J Agron. 2010;33(1):33–42.

    Article  Google Scholar 

  80. Moragues M, Zarco-Hernández J, Moralejo MA, Royo C. Genetic diversity of glutenin protein subunits composition in durum wheat landraces [Triticum turgidum ssp. turgidum convar. Durum (desf.) MacKey] from the Mediterranean basin. Genet Resour Crop Evol. 2006;53(5):993–1002.

    Article  CAS  Google Scholar 

  81. Boukef S, McDonald BA, Yahyaoui A, Rezgui S, Brunner PC. Frequency of mutations associated with fungicide resistance and population structure of Mycosphaerella graminicola in Tunisia. Eur J Plant Pathol. 2012;132(1):111–22.

    Article  CAS  Google Scholar 

  82. Hassine M, Siah A, Hellin P, Cadalen T, Halama P, Hilbert J-L, Hamada W, Baraket M, Yahyaoui A, Legrève A. Sexual reproduction of Zymoseptoria tritici on durum wheat in Tunisia revealed by presence of airborne inoculum, fruiting bodies and high levels of genetic diversity. Fungal biology. 2019;123(10):763–72.

    Article  CAS  PubMed  Google Scholar 

  83. Suffert F, Delestre G, Gélisse S. Sexual reproduction in the fungal foliar pathogen Zymoseptoria tritici is driven by antagonistic density dependence mechanisms. Microb Ecol. 2019;77(1):110–23.

    Article  PubMed  Google Scholar 

  84. Chedli RBH, M’Barek SB, Yahyaoui A, Kehel Z, Rezgui S. Occurrence of Septoria tritici blotch (Zymoseptoria tritici) disease on durum wheat, triticale, and bread wheat in Northern Tunisia. Chil J agricultural Res. 2018;78(4):559–68.

    Google Scholar 

  85. Morais D, Gélisse S, Laval V, Sache I, Suffert F. Inferring the origin of primary inoculum of Zymoseptoria tritici from differential adaptation of resident and immigrant populations to wheat cultivars. Eur J Plant Pathol. 2016;145(2):393–404.

    Article  Google Scholar 

  86. Olson MV. When less is more: gene loss as an engine of evolutionary change. Am J Hum Genet. 1999;64(1):18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Gepts P, Papa R. Evolution during domestication. e LS 2001.

  88. Rahman S, Islam S, Yu Z, She M, Nevo E, Ma W. Current progress in understanding and recovering the wheat genes lost in evolution and domestication. Int J Mol Sci. 2020;21(16):5836.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Roussel V, Koenig J, Beckert M, Balfourier F. Molecular diversity in french bread wheat accessions related to temporal trends and breeding programmes. Theor Appl Genet. 2004;108(5):920–30.

    Article  CAS  PubMed  Google Scholar 

  90. Thomas M, Dawson JC, Goldringer I, Bonneuil C. Seed exchanges, a key to analyze crop diversity dynamics in farmer-led on-farm conservation. Genet Resour Crop Evol. 2011;58(3):321–38.

    Article  Google Scholar 

  91. Chakraborty R, Weiss KM. Admixture as a tool for finding linked genes and detecting that difference from allelic association between loci. Proceedings of the National Academy of Sciences 1988, 85(23):9119–9123.

  92. Medini M, Hamza S. Pathotype and molecular characterisation of Mycosphaerella graminicola isolates collected from Tunisia, Algeria and Canada. J Plant Pathol. 2008;90(1):65–73.

    CAS  Google Scholar 

  93. Ferjaoui S, Aouini L, Slimane RB et al. Deciphering resistance to Zymoseptoria tritici in the Tunisian durum wheat landrace accession ‘Agili39’. 2022. BMC Genomics 23, 372 (2022).

  94. Zadoks JC, Chang TT, Konzak CF. A decimal code for the growth stages of cereals. Weed Res. 1974;14(6):415–21.

    Article  Google Scholar 

  95. Maccaferri M, Ricci A, Salvi S, Milner SG, Noli E, Martelli PL, 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–63.

    Article  CAS  PubMed  Google Scholar 

  96. Powell W, Morgante M, Andre C, Hanafey M, Vogel J, Tingey S, Rafalski A. The comparison of RFLP, RAPD, AFLP and SSR (microsatellite) markers for germplasm analysis. Mol Breeding. 1996;2(3):225–38.

    Article  CAS  Google Scholar 

  97. Team RC. R: A language and environment for statistical computing. 2013.

  98. Nei M. Genetic distance between populations. Am Nat. 1972;106(949):283–92.

    Article  Google Scholar 

  99. Langella O. Populations a free population genetic software. Version 1.2. 31. In.; 2002.

  100. Page RD. Tree view: an application to display phylogenetic trees on personal computers. Bioinformatics. 1996;12(4):357–8.

    Article  CAS  Google Scholar 

  101. Peakall R, Smouse PE. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6(1):288–95.

    Article  Google Scholar 

Download references


We are thankful to Shiaoman Ciao from the USDA small grains genotyping laboratory in Fargo, ND for conducting SNP genotyping, and to Aouini Lamia for her help in editing the manuscript.


Field experiments supported by the federated project untitled ‘Identification of durum wheat resistant genotypes to biotic and drougth stress and their valorization for sustainable agriculture’ acronym RESIDUR, supported by IRESA under the Tunisian Ministry of Agriculture.

Genotyping was supported by the International Treaty for Plant Genetic Resources for Food and Agriculture/FAO (FAO/ITPGRFA) and European Union.

M.O. was supported by the PhD fellowship from the Ministry of Higher Education and Scientific Research - Tunisia.

Author information

Authors and Affiliations



M.O. and B.B. conducted the data collection, carried out the analyses and the interpretation of the results. M.O. drafted the manuscript. S.F. contributed to the data collection. M.M. provided the studied durum wheat accessions from the BNG. S.M.U. participated on organizing the genotyping with Illumina 90 K bread wheat array. S.H. contributed in the conception, design and coordination of the study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Sonia Hamza.

Ethics declarations

Ethics approval and consent to participate

The used plant material consist on Triticum durum landraces registered at the Tunisian gene bank (BNG) and provided by the co-author M.M. Experimental research and field studies conducted at the National Institute of Agronomy (INAT) and at the Regional Field Crops Research Center of Beja (CRRGC) including the collection of plant material complies with relevant institutional, national and international guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ouaja, M., Bahri, B.A., Ferjaoui, S. et al. Unlocking the story of resistance to Zymoseptoria tritici in Tunisian old durum wheat germplasm based on population structure analysis. BMC Genomics 24, 328 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: