Evaluating agromorphological traits of Greek wheat landraces and exploring their potential for bread and pasta making based on seed physical properties

: Physical and morphological characteristics of Greek wheat landrace and cultivar seeds were evaluated, aiming to assess their bread-and pasta-making quality. Furthermore, the plant agromorphological traits of wheat landraces and cultivars were measured and correlated with seed physical properties. Fifteen Triticum spp. accessions, out of which four cultivars (two Triticum aestivum and two T. durum ), nine landraces (two T. durum , ﬁve T. aestivum , one T. dicoccum and one T. polonicum ) and two T. durum of unknown status were studied. Seventeen morphological plant and seed traits were measured based on UPOV descriptors. Ear emergence was earlier in cultivars than in landraces. The plant weight of the landraces was, in some cases, almost twice the cultivars’. Seed ﬁrmness (62.96–194.85N) was positively correlated with thousand kernel weight and volume (0.840 and 0.791, P < 0.05, respectively). Based on the agromorphological traits, cluster analysis clearly separated the species and the cultivars from the landraces, and the unknown status accessions were grouped with the cultivars. Cluster analysis with all measured traits set the wheat accessions in the following distinct groups: (1) T. durum cultivars including the unknown status accessions, (2) T. aestivum cultivars, and (3) T. aestivum and T. durum landraces. T. polonicum and T. dicoccum each formed a separate group. Seed physical properties of the analyzed Greek landraces indicated their suitability for bread and/or pasta making.


Introduction
Many researchers in the past have tried to define the meaning of a landrace, which is also known as "a local, farmer's, traditional variety or a local population" (Zeven, 1998).Recently, a more general definition for landrace was suggested: "a cultivated, genetically heterogeneous variety that has evolved in a certain ecogeographical area and is therefore adapted to the edaphic and climatic conditions and to its traditional management and uses" (Casañas et al, 2017).
Landraces have been the main crop genetic resources used since the emergence of agriculture (Purugganan, 2019).In the last century, with agricultural intensification, several landraces have been replaced with highyielding cultivars (hereafter the term cultivar means scientifically improved cultivated variety).The main differences between landraces and cultivars are that the latter are the results of formal breeding research (when compared to farmers' mass selection for landraces) and they are provided to farmers and growers through formal seed systems (when compared to traditionally infor-mal channels, such as exchanges and local markets for landraces).Landrace cultivation has been almost abandoned and many landraces have been lost, leading to severe genetic erosion (Robbana et al, 2021), which is the loss of crop variation due to the modernization of agriculture (Van De Wouw et al, 2010).
In the case of wheat, a study conducted by Haudry et al (2007), concluded that the estimated initial biodiversity decreased by 84% for durum wheat and by 69% for bread wheat, although such estimations could be controversial (Khoury et al, 2022).Genetic erosion of Triticum aestivum landraces in Ethiopia was found to be 66.7% in 2016 (Gelelcha et al, 2023).
In Greece, only 3% of the total land surface cultivated with durum wheat was planted with landraces in 1984, and 1.3% in 1991 (Kyzeridis et al, 1995).On an expedition encompassing 36 villages on Lefkada island, from 2010 to 2012, researchers were able to collect only three wheat accessions, two of which were landraces (Thomas et al, 2013).Although the genetic erosion of wheat landraces was shown by collecting expeditions carried out in Greece, they can still be found cultivated on farms (Thanopoulos et al, 2021).The study of landrace properties is crucial both for local communities, which would benefit from the added value of their products, and their use in breeding programmes.Thanopoulos et al (2021), who collected landraces in Arcadia, Greece, estimated wheat genetic erosion between 1942 and 2021 to be 45.7%.Moreover, recent research revealed a 97% genetic erosion of winter wheat for the years 1955-2015 in Central Europe (Cseh et al, 2021).
Cereal cultivars are high-yielding, but the majority of them is less resilient to abiotic and biotic stresses compared to landraces (Newton et al, 2010).Reduced adaptation to organic, low-input production farming systems and poorer genetic diversity were observed in Tunisian improved durum wheat cultivars, in contrast to landraces (Ayadi et al, 2020).Although further research is required, there is also evidence that ancient wheat types, such as emmer wheat, provide health benefits compared to modern wheat cultivars (Shewry, 2018) and could be a potential source of antioxidants, carotenoids and phenolics (Newton et al, 2010).For such ancient species, like T. dicoccum, and neglected and underutilized species, like T. polonicum, (Bieńkowska et al, 2020;Carvalho et al, 2019) more research to investigate their properties is needed.
Currently, climate change, extensive loss of genetic diversity, increased vulnerability of cultivars to abiotic and biotic stresses and consumer awareness of the value of local products, have raised interest in wheat landraces (Brush, 1995;Moghaddam et al, 1997;Saleh, 2020;Tomás et al, 2020).Wheat landraces could play a crucial role in mitigating the negative effects of climate change due to their ability to perform sufficiently under high temperatures (Ulukan, 2021) and exhibit higher yield stability under water stress (Karamanos et al, 2017).Recent research focused on the collection of landraces, their traits and cultivation revival as well as the promotion of sustainable agricultural production (Newton et al, 2010).It also highlighted the characteristics associated with their utilization in the production of high-quality products for consumers (López-Fernández et al, 2021;Ruisi et al, 2021).
The analysis of the agronomic performance of wheat landraces revealed significant differences in plant height, seed colour and biological cycle (Gharib et al, 2021).A study of morphological characteristics of old Sicilian landraces showed high variability among durum wheat accessions (Sciacca et al, 2014).Similar results were also obtained for days to emergence, tillering, booting, flowering and physiological maturity, plant height, and thousand kernel weight (TKW) in 23 Moroccan durum wheat accessions (Zarkti et al, 2010).In a study of 63 durum wheat landraces, grain weight was the most significant criterion for improving their adaptation to the northern parts of the Mediterranean basin (Moragues et al, 2006).
Several studies compared grain characteristics of landraces and other variety types.In one study, ten obsolete varieties (varieties that were developed by systematic breeding efforts in the past and now have been replaced by newer and better-performing cultivars) (Maxted et al, 2011) and two wheat landraces from Bosnia and Herzegovina were evaluated for grain length, width, thickness, volume, surface area, TKW and protein content, and most of the traits presented significant interactions among them.Obsolete wheat genotypes 'Brkulja' and ' Šianka' showed the most promising grain characteristics (Kondić et al, 2020).Kernel morphology is important for manufacturing final food products, which requires specific grain characteristics (Dholakia et al, 2008).Machine-vision techniques (image analysis) were used as a useful tool for the classification of varieties based on morphology, colour and seed texture of different species (Smykalova et al, 2011(Smykalova et al, , 2013;;Chaugule and Mali, 2016) and wheat seed quality assessment (Venora et al, 2009).Other studies evaluated quality parameters, such as starch and protein content of wheat flour obtained from landraces in comparison to cultivars.Boukid et al (2018) detected higher values of protein and gliadin in Tunisian durum wheat landraces, while modern genotypes had significantly higher values of total starch.Zheng et al (2012) studied 477 Chinese bread wheat landraces (T.aestivum L.) and concluded that there was a significant correlation between starch and protein content.Mediterranean durum wheat landraces presented a higher amount of protein content but lower gluten strength than cultivars (Roselló et al, 2018).In addition, a lower content of water-extractable arabinoxylans, the main non-starch polysaccharides in cereals, was found in Turkish landraces (T.monococcum, T. durum and T. aestivum) than in cultivars (T.aestivum) (Cetiner et al, 2020).The comparison of different aromatic profiles and the characteristics of the processed products showed high diversity in Sicilian durum wheat landraces in contrast to the uniformity of cultivars (Ruisi et al, 2021).Recently, some studies focused on the physical or functional characteristics of final products made from wheat landraces.For example, some soft wheat landraces could be used for biscuit production, since their grain hardness did not significantly reduce biscuit quality (Igrejas et al, 2002).Spanish landraces showed lower bread-making quality properties because of their low-quality protein compared to cultivars (Gómez et al, 2009).High variability of Mediterranean landraces could be used for improving gluten quality (Nazco et al, 2014).
Even though many papers separately examine either plant agromorphological traits or seed physicomechanical properties, there is limited published work that studies holistically the above traits with both agronomic and food technology approaches.For instance, in a comparative study of five cultivars and two landraces, the latter were significantly taller, flowered later, and had lower yields and gluten index but higher protein content in comparison to wheat cultivars (Preiti et al, 2022).This field needs more data to critically evaluate this issue.The purpose of the present study was to characterize wheat landraces and cultivars according to morphological, agronomical, seed quality and technological traits and investigate their potential for bread and/or pasta making.Specifically, the research aimed to: (1) study the agromorphological characteristics and seed physical properties, (2) compare T. durum, T. aestivum, T. dicoccum, and T. polonicum wheat landraces with cultivars, (3) classify all the accessions (cultivars and landraces) based on agromorphological traits, separately or on the whole set of traits (both agromorphological and seed physical properties) and (4) identify the wheat landraces more suitable for breadand pasta-making purposes.

Plant material
The plant material consisted of 15 Triticum spp.accessions, namely two cultivars and five landraces of T. aestivum, two cultivars, two landraces and two accessions of unknown status of T. durum, one landrace of T. polonicum and one landrace of T. dicoccum.The latter was collected in the early 20th century in Greece, it was conserved in USDA-GRIN (Code PI 94682) and repatriated recently.Landraces were collected from different regions in Greece and donated to the Plant Breeding and Biometry Laboratory of the Agricultural University of Athens (AUA) (Table 1 and Figure 1).The Institute of Plant Breeding and Genetic Resources (IPGRB) donated four cultivars, recommended for bread making ('Yekora' and 'Elisavet' belonging to T. aestivum) and pasta making ('Mexicali 81' and 'Elpida', of T. durum) which were used as controls.1).

Experimental design
The experiment was conducted at an experimental field of the Agricultural University of Athens (location: 37 • 59' 06.8" N, 23 • 42' 24.7" E, altitude 24m a.s.l.) during the winter growing season of 2019-2020.The meteorological data are presented in Figure 2. The values of mean air temperature and precipitation for the Athens region during the years 2015-2020 were provided by the Hellenic National Meteorological Service (HNMS, 2023).A randomized complete block design with three replications was used.Each replication consisted of 20 plants of each accession sown in rows spaced 0.2m apart.Sowing distances among plants in each row was 20cm.The sowing date was mid-November, a typical period for sowing in central Greece.
Manipulation of sowing time is one of the most effective techniques for the production of high-quality wheat kernels (Butkovskaya and Kozulina, 2021).In the present study, for comparison reasons, plant sowing was carried out simultaneously for all accessions, without considering the best sowing time for each accession.

Measurement of agromorphological traits
During vegetative and reproductive stages, time of ear emergence (days), total number of tillers, ear length excluding awns (cm), ear length including awns (cm), number of spikelets in the spike of the first tiller, plant length (measured from the base of the plant to the tip of the highest awn, cm), stem length (measured from plant base to ear base, cm) and whole plant weight (physically dried by sun, in g/plant) were measured.In addition, the qualitative variables ear colour, ear shape in profile, growth habit and awn colour, were scored.Agromorphological traits were evaluated on all the plants per plot to calculate the mean value for the quantitative traits and the median for the qualitative traits.The measurements of the following traits: plant  length, growth habit, ear emergence, ear length, ear colour, ear shape, and awn colour were carried out using UPOV descriptors (UPOV, 2012(UPOV, , 2017)).

Physical properties of seeds
Before measurements, shrunken, broken kernels, and other impurities were removed.The weight of 1,000 randomly selected kernels of wheat was measured for thousand kernel weight (TKW).The same 1,000 selected kernels were then filled into a volumetric cylinder to determine their volume.Tap volume was measured after gently tapping the cylinder 100 times or until there was no further decrease in the sample level.Bulk density (ρ b ) and tap density (ρ t ) were calculated as weight of sample per volume and per tap volume of sample, respectively (g/ml) (Chaloulos et al, 2021).Determination of flow properties is essential in processes such as transportation, mixing and storage (Fitzpatrick et al, 2004).The Carr index (CI) was used to calculate the flowability of the seeds (Carr, 1965): The firmness of seed (SF), defined as the force required to crush the seed, was estimated using an Instron (Universal Testing Machine, Model 3343, Nordwood, MA, USA) equipped with a 1kN load cell.Fifty intact seeds of each accession were compressed to 50% of their initial width with a 4cm diameter probe and a speed of 1mm/min.
Seed colour was measured using the colourimeters 3nh High-Quality Spectrophotometer NS800S (Shenzhen 3nh Technology, China), according to CIE-L*a*b* uniform colour space, where L* indicates lightness, a* indicates hue on a green (−) to red (+) axis, and b* indicates hue on a blue (−) to yellow (+) axis.White index (WI) was also calculated according to the following equation:

Image analysis-shape factor measurements
One hundred seeds of each accession were scanned twice using a flatbed scanner (HP scan jet 4370, Hewlett Packard, USA).Images of the seeds were captured for further analysis using an image analysis software (Image Pro Plus 7, Media Cybernetics, USA).An indicative image of the scanned seeds is presented in Figure 3A.Area (mm 2 ), aspect, average optical density (density/intensity mean), mean, maximum and minimum diameter (mm), perimeter (mm), roundness, size length (mm) and width (mm) were determined.Detailed descriptions and definitions of shape factors are presented in Figure 3B.The measurements were done in triplicates.

Statistical analysis
Statistical analysis of the data was performed with Statgraphics Centurion XV (Statgraphics, Rockville, MD, USA).For qualitative variables (plant growth habit, colour of ear, ear shape in profile and awn colour) frequency percentages (Supplemental Figure 1) and median values were calculated (Supplemental Table 1).
Friedman Test was used to test the null hypothesis that the mean ranks of the samples groups are the same (Supplemental Table 1).
Quantitative trait data were subjected to analysis of variance after testing the assumptions of normality and homoscedasticity with Levene's and Bartlett's tests respectively (Supplemental Table 2).The comparisons of the means were performed using the Tukey HSD criteria with a level of significance of α = 0.05.
The correlation between quantitative parameters was assessed by Pearson's correlation test (significance level at α = 0.05).The correlation between qualitative parameters (nonparametric measures) was assessed by Spearman rank correlation test (significance level at α = 0.05).To correlate quantitative and qualitative variables the Spearman test was also used (Supplemental Table 3).Cluster analysis was used to group the quantitative observations with similar characteristics using nearest neighbour (single linkage) and distance metric of squared Euclidean on standardized data.

Agromorphological traits
Mean values of quantitative agromorphological traits for all studied wheat accessions are presented in Table 2. Time of ear emergence ranged from 107 to 138 days (Yekora -T.aestivum, cultivar, and Zoulitsa -T.aestivum, landrace, respectively).All cultivars, except Elisavet -T.aestivum, were earlier than landraces.Similarly, Frankin et al (2021) reported cultivars to be on average 12% earlier than landraces.On the other hand, the landraces Zoulitsa (T.aestivum), Kopaida (T.dicoccum) and Leventis (T.polonicum) had the latest heading time of all accessions.The T. aestivum landraces, Asprositi Kozani (W12), Ntopio (W6) and Asprositi Kalavrita (W11) exhibited early maturation compared to other landraces.Consequently, these landraces have the potential to be incorporated into breeding programmes focused on selecting for early maturation traits (Javaid et al, 2005).Time of ear emergence was positively correlated with total number of tillers, plant length (from the base of the plant to the tip of the highest awn) and stem length (0.584, 0.674 and 0.710, P < 0.05, respectively) (Table 3).A positive correlation between days of heading and number of fertile tillers was also recorded in a study with 64 bread wheat genotypes in West Shewa and found that days to heading were positively correlated with the number of productive tillers (Mecha et al, 2017).In the study of Bilgrami et al (2020), where bread wheats were used, the total tiller number was also significantly correlated with days from sowing to heading (0.22, P < 0.05).Although landraces had more tillers than cultivars, the difference between landraces and specific cultivars was not statistically significant and the number of tillers ranged from 4.18 (Mavragani Lemnos -T.durum, unknown status, W13) to 10.38 (Grilos -T.aestivum, landrace, W8).
In most cases, the landraces had higher total plant weight than the cultivars.The weight of the landraces ranged from 28.97g/plant (Asprositi Kalavrita -T.aestivum, landrace, W11) to 73.01g/plant (Asprositi Kozani -T.durum, landrace, W12), while cultivars' weight ranged from 17.55g/plant (Yekora -T.aestivum, cultivar, W1) to 25.33g/plant (Elpida -T.durum, cultivar, W3).This could be attributed to the higher biomass produced by landraces, a characteristic desirable by farmers who need straw for their animals (Thanopoulos et al, 2021).Plant weight was positively and moderately correlated with the total number of tillers and ear length including awns (0.570 and 0.526, P < 0.05, respectively) and highly correlated with the number of spikelets, plant length and stem length (0.644, 0.796 and 0.759, P < 0.01, respectively).The number of spikelets per spike of the first tiller did not differ significantly between accessions.
Cultivars were significantly shorter than the landraces, which is due to the introduction of Rht (Reduced height) genes in modern wheat varieties.These dwarfing genes have been used in wheat breeding to develop cultivars with short stature (Acquaah, 2012).Similar results were recorded in another study, where Italian landraces were found to have taller culms compared to cultivars (Preiti et al, 2022).The plant length for the cultivars was between 83.88cm (Yekora -T.aestivum, cultivar, W1) and 103.87cm (Mexicali 81 -T.durum, cultivar, W4), while for the landraces the plant length was between 118cm (Kopaida -T.dicoccum, landrace, W10) and 159.04cm (Aspratheri -T.durum landrace, W14).Mavragani accessions, W7 and W13, of unknown breeding status, were short, with 87.07cm and 93.18cm, respectively.These values did not differ significantly from cultivar plant length (W1 -W4) (Table 2).Low plant height is a desirable trait for commercial varieties, as it is correlated with good standing ability and resistance to lodging (Würschum et al, 2017).The weak point of wheat landraces is the lodging, but under climatic change and high temperatures, taller plants could demonstrate higher yield than short ones (Jatayev et al, 2020).
Plant growth habit was significantly correlated with weight per plant, total number of tillers, plant length and stem length (base to ear) (0.525, 0.623, 0.538,  0.529, respectively P < 0.05) (Supplemental Table 3).A higher number of tillers leads to an increased inclination of the plant.Cluster analysis based on the quantitative agromorphological traits grouped all cultivars with two Mavragani accessions of unknown status (W7 and W13) in cluster 1, while all landraces were included in a distinct cluster (cluster 2) (Figure 4).The similarity of Mavragani durum wheat accessions from Skyros and Lemnos, with wheat cultivars indicates that these accessions could be old cultivars grown and reproduced for years by farmers, who considered them of local origin (Thomas et al, 2013).Overall, the quantitative agromorphological traits based on UPOV descriptors were sufficient to separate wheat accessions by their breeding status, namely landraces or cultivars.

Physical properties of seeds
TKW is of great interest (Javaid et al, 2005) as it is an important index for the prediction of flour extraction rate (Posner, 2009).Higher values of TKW indicate a higher percentage of endosperm leading to higher flour yield (Wiersma et al, 2001).TKW values of the studied accessions are presented in Figure 5. TKW of W9 (T.polonicum) was significantly higher (80.45g) than all other accessions.W10 (T.dicoccum) presented the lowest value of TKW (26.92g), consistent with previous research, where the TKW of 38 emmer wheat accessions, collected in several European countries, ranged from 22.9g to 42.6g (Mondini et al, 2013).T. dicoccum presented significantly lower values of TKW compared to T. durum, 32.02-36.12g and 45.17-46.31g, respectively (Shoormij et al, 2022).TKW in the rest of the accessions, either cultivars or landraces, ranged from 34.08g (Elisavet -T.aestivum, cultivar, W2) to 67.16g (Mavragani Lemnos -T.durum, unknown status, W13) and were in accordance with Travlos et al (2012) who studied five Greek bread wheat landraces (30-50.7g).Specifically, in the present study, W5 (Zoulitsa) had an average TKW of 43.66g, while in Travlos et al (2012), TKW in Zoulitsa ranged from 36g to 44g, indicating that climatic conditions can affect TKW but also the origin of the genetic material (Papadakis, 1929).Moreover, it is worth noting that in the present study, T. aestivum in most cases presented significantly lower values of TKW compared to T. durum accessions.Landraces had similar or higher values of TKW compared to cultivars of the same species, indicating that landraces could give high flour yield.However, results from previous research are contrasting as no significant difference between landraces' and cultivars' TKW was reported for Turkish (Cetiner et al, 2020) and Sicilian wheats (Ruisi et al, 2021).On the contrary, Spanish bread wheat landraces showed higher TKW compared to a reference set (López-Fernández et al, 2021).In another research, similar TKW values for T. aestivum landraces and cultivars were detected, while T. durum cultivars, compared to landraces, had significantly higher TKW values (Frankin et al, 2021).
Due to their high heritability grain perimeter, area and volume are important indirect indices for grain yield improvement (Abdipour et al, 2016).Seed volume ranged from 48ml to 113ml (Elisavet -T.aestivum, cultivar, W2 and Leventis -T.polonicum, landrace, W9, respectively) (Table 4).It is worth noting that W9 (T.polonicum) had significantly the highest values of TKW, volume and seed firmness (SF).Although Papadakis (1929) mentioned that T. polonicum did not provide any agronomic value, our findings show that this species has bread-making value.
There are many methods to determine the required force to evaluate grain mechanical properties.In the present study, a rapid measurement of the mechanical properties of a single grain was used, which provided valuable results (Shewry and Hey, 2016).SF ranged from 62.96N to 194.85N (W10, Kopaida -T.dicoccum, landrace and W9, T. polonicum, landrace, respectively) and was correlated positively with seed weight and volume (0.840 and 0.791, P < 0.05, respectively) (Table 5).High values of weight and hardness indicate healthy seeds.Seed texture determines, to a great extent, milling properties and end uses.Softer seeds require more energy to mill and produce finer flours suitable for biscuit production, while flour from harder seeds is more appropriate for bread making (Liu, 2008).As W10 (Kopaida -T.dicoccum, landrace) had the lowest SF values, it could be suitable for cake and biscuit production.A recent research concluded that in general old cultivars of T. aestivum were softer than modern cultivars (Cetiner et al, 2020).In the present research, this statement was not fully confirmed, as there was no clear trend.However, it should be mentioned that fewer old and modern varieties were used in the present work compared to the earlier study.
Bulk density could be used as a criterion for yield improvement, as it was correlated positively with grain yield (Karimizadeh et al, 2012).In the present study, density did not differ significantly between accessions except for W10 (Kopaida -T.dicoccum, landrace), having the lowest value (0.55g/ml).
In general, colour parameters did not differ significantly among accessions and were not correlated with the other physical properties (Table 5).The exceptions were the green/red hue (a*), which correlated negatively with volume (-0.516,P < 0.05), indicating that seeds were redder as their volume increased, and lightness (L*), blue/yellow hue (b*) and white index (WI), which were correlated positively with density (0.649, 0.719 and 0.565, P < 0.05, respectively).Increased density led to increased lightness, yellowness and whiteness of the seeds.W10 (Kopaida -T.dicoccum, landrace) was the only accession with a significantly lower value of b*   among all other accessions.Colour determines to a great extent the overall consumer acceptance of food (Mandala and Protonotariou, 2021).Wheat kernel colour is related to different milling and baking attributes of the final products (Dowell, 1998).Bright yellow or amber colour of durum wheat results in high-quality pasta products (Cole et al, 1991), thus sample W10 (Kopaida -T.dicoccum, landrace) probably is the least suitable for pasta making among the studied accessions.Probably, W9 (Leventis -T.polonicum, landrace) and W13 (Mavra-gani Lemnos -T.durum, unknown status) could be the most suitable for pasta making considering their high values of b*, TKW and seed volume.However, further investigation is needed.Mechanical properties and sensory evaluation of produced pasta could verify our claim.
The determination of kernel morphology is of great importance for bread making, as kernel shape and uniformity may influence the milling quality (Campbell et al, 1999).Image analysis could play a key role in distinguishing different wheat varieties by defining the product traceability for wheat landraces (Khoshroo et al, 2014;Grillo et al, 2017).In the present research, mean seed diameter ranged between 4.17mm (W2, Elisavet -T.aestivum, cultivar) and 5.66mm (W9, Leventis -T.polonicum, landrace).Particle size distribution curves differed between accessions, indicating that some landraces had a more uneven seed population than cultivars (Figure 6).Curve of W12 (Asprositi Kozani -T.durum, landrace) presented an almost bimodal distribution, indicating the presence of two distinguished populations of seeds, with different sizes, which could be explained by the fact that sometimes landraces consisted of two or more populations with different morphological characteristics (Papadakis, 1929).This could be linked with seed heteromorphism, characteristic of wild populations.A similar curve was presented by W13 (Mavragani Lemnos -T.durum, unknown status), where the bimodal distribution was less evident, but also presented a shoulder on the left (Figure 6).In contrast, W3 (Elpida -T.durum, cultivar) presented the narrowest unimodal distribution curve, with the peak around 5.5mm.All the other accessions presented frequency curves similar to W2 (Elisavet -T.aestivum, cultivar) or W9 (Leventis -T.polonicum, landrace).The main characteristic of these curves is that distribution is broad, due to the uneven seed populations.W10 (Kopaida -T.dicoccum, landrace) presented a similar curve to W2 (Elisavet -T.aestivum, cultivar), despite the fact that they belonged to different species.
Bread wheat cultivars (W1, Yekora -T.aestivum and W2, Elisavet -T.aestivum) had significantly more roundness and shorter seeds compared to pasta wheat cultivars (W3, Elpida -T.durum, and W4, Mexicali 81 -T.durum), as referred to in the literature (Campbell et al, 1999).Perimeter was significantly lower for W2, Elisavet -T.aestivum, cultivar (15.71mm) and significantly higher for W9, Leventis -T.polonicum, landrace (23.62mm), compared to all other accessions (Supplemental Table 4).Similar results were observed for seed length, with W9 having the longest seeds (9.3mm) and W2 the shortest (5.96mm).T. polonicum (represented by W9) had characteristically long seed (Percival, 1921).Regarding the width of the seeds, the differences were not significant in all the cases and the values ranged from 2.72mm to 3.94mm (W10 and W13, respectively).Grain width could be used as a predictive index for mean grain weight determination (Haghshenas et al, 2022).However, in the present study, the most valuable parameter was the mean diameter.

Correlations among all parameters
In the present study, TKW was significantly positively correlated with area of seed, mean diameter, SF, perimeter, length, width and volume of the seeds (0.871, 0.887, 0840, 0.786, 0.677, 0.829, and 0.975, P < 0.05, respectively, Supplemental Table 5).There is previous research which also associated kernel weight with its width and length (Campbell et al, 1999;Kondić et al, 2020), but some of them could not associate these parameters when they studied only bread wheats (Schuler et al, 1995).Moreover, Javaid et al (2005) found that TKW was positively correlated with plant height in Pakistani bread wheat landraces, and Moghaddam et al (1997) found a significant correlation between TKW and number of tillers in Irani bread wheat landraces, which was not observed in the present study, suggesting possibly different patterns.
There are several studies evaluating wheat kernels by compression tests (Ponce-García et al, 2016).However, there is a research gap on the correlation between wheat seed mechanical properties with plantagromorphological and seed-morphological traits.In the present study, SF was significantly correlated with plant weight, and seed area, mean diameter, perimeter, length, width, volume and TKW (0.626, 0.711, 0.709, 0.688, 0.615, 0.570, and 0.840, P < 0.05, respectively).Thus, SF could be used as an additional tool to predict all the above parameters.Fitting a simple linear regression model to describe the relationship between SF and TKW explained 70.58% of firmness variability (P < 0.05) (Supplemental Figure 2).

Grouping of accessions
Cluster analysis, using all quantitative data, allowed us to group the wheat accessions in accordance with their species and breeding status (Figure 7).W9 (Leventis -T.polonicum, landrace) and W10 (Kopaida -T.dicoccum, landrace) created a distinct clade each, as anticipated, considering that they were different species.The control cultivars of T. durum species (W3 and W4) created a distinct cluster, along with the T. durum accessions of unknown breeding status W7 and W13 (Mavragani -T.durum, unknown status from Skyros and Lemnos, respectively).Interestingly, except for W11 (Asprositi Kalavrita) all landraces, whether T. aestivum or T. durum, formed a separate cluster, related to the T. aestivum cultivars.
It is worth noting that with the addition of physical properties in the cluster analysis, the grouping of Triticum spp.accessions changed compared to Figure 4, where only observations according to agromorphological traits were used.The cluster of landraces was closer to T. aestivum cultivars, indicating that landracesregardless of the species they belonged to -could be more suitable for bread making, considering their physical properties, but also the main preference and use (mainly for bread) by farmers in Greece (Douma et al, 2016).

Conclusion
The agromorphological traits recorded in this study proved to be sufficient to highlight the differences between wheat species, landraces and cultivars, grouping them into two different clusters.TKW did not differ significantly between landraces and cultivars, indicating that landraces could also give a high flour yield.Increased values of TKW and seed volume indicated higher values of seed firmness.Cluster analysis of all data separated the wheat accessions into groups accord- ing to their attributes, indicating that Greek wheat landraces, either T. durum or T. aestivum could be probably used for bread making.Regarding the W10 (T.dicoccum, landrace) and W9 (T.polonicum, landrace) wheat accessions, their incorporation in both bread and pasta making would be of great interest, as they presented distinct characteristics.The evaluation of both agromorphological and physical properties of wheat seeds provided valuable information that could be used to distinguish wheat varieties for pasta and bread making based on quality characteristics.

Figure 1 .
Figure 1.Map of Greece with collection sites of wheat landraces characterized in this study (for code explanations see Table1).

Figure 2 .
Figure 2. Mean air temperature ( • C) and total monthly precipitation (mm) for Athens during the growing period (November 2019 to July 2020) and the average of the last five years (2015-2020).Data provided by HNMS (2023).

Figure 3 .
Figure 3. A) Indicative image from scanned seeds processed with Image Pro Plus 7 analysis software and B) shape factors definitions and schematic representations.

Table 3 .
Correlation coefficients between quantitative agromorphological traits of the studied wheat accessions.The P-values are in brackets.Values in bold denote significance level at

Table 5 .
Correlation coefficients between seed physical properties of all the studied wheat accessions.The P-values are in brackets.a*, hue on a green (−) to red (+) axis; b*, hue on a blue (−) to yellow (+) axis; CI, Carr Index; SF, Seed Firmness; L*, lightness; TKW, Thousand Kernel Weight; and WI, White Index.Values in bold denote significance level at P < 0.05.