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

HAL
open science

# Barley systematics and taxonomy foreseen by seed morphometric variation

Angèle Jeanty, Laurent Bouby, Vincent Bonhomme, François Balfourier, Clément Debiton, Camille Dham, Sarah Ivorra, Jérôme Ros, Allowen Evin

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

## To cite this version:

Angèle Jeanty, Laurent Bouby, Vincent Bonhomme, François Balfourier, Clément Debiton, et al. Barley systematics and taxonomy foreseen by seed morphometric variation. PLoS ONE, 2023, 10.1371/journal.pone.0285195. [https://doi.org/10.1371/journal.pone.0285195](https://doi.org/10.1371/journal.pone.0285195). [hal-04141471v2](https://hal.science/hal-04141471v2)

HAL Id: hal-04141471

https://hal.science/hal-04141471v2

Submitted on 19 Sep 2023

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
To cite this article

Jeanty A, Bouby L, Bonhomme V, Balfourier F, Debiton C, Dham C, et al. (2023) Barley systematics and taxonomy foreseen by seed morphometric variation. PLoS ONE 18(5): e0285195.
https://doi.org/10.1371/journal.pone.0285195

# Barley systematics and taxonomy foreseen by seed morphometric variation

Angèle Jeanty<fnref n="1" />*, Laurent Bouby<fnref n="1" />, Vincent Bonhomme<fnref n="1" />,<fnref n="2" />, François Balfourier<fnref n="3" />, Clément Debiton<fnref n="3" />, Camille Dham<fnref n="1" />, Sarah Ivorra<fnref n="1" />, Jérôme Ros<fnref n="1" />, Allowen Evin<fnref n="1" />

## Abstract

Since its Neolithic domestication in the Fertile Crescent, barley has spread to all continents and represents a major cereal in many modern agrarian systems. Current barley diversity includes thousands of varieties divided into four main categories corresponding to 2-row and 6-row subspecies and naked and hulled types, each of them with winter and spring varieties. This diversity is associated to different uses and allow cultivation in diverse environments. We used a large dataset of 58 varieties of French origin, (1) to assess the taxonomic signal in barley grain measurements comparing 2-row and 6-row subspecies, and naked and hulled types; (2) to test the impact of the sowing period and interannual variation on the grains size and shape; (3) to investigate the existence of morphological differences between winter and spring types; and finally (4) to contrast the relationship between the morphometric and genetic proximity. Size and shape of 1980 modern barley caryopses were quantified through elliptic Fourier Transforms and traditional size measurements. Our results indicate that barley grains record morphological diversity of the ear (89.3% classification accuracy between 2-row/6-row subspecies; 85.2% between hulled and naked type), sowing time of the grains (from 65.6% to 73.3% within barley groups), and environmental conditions during its cultivation and varietal diversity. This study opens perspectives for studying archaeological barley seeds and tracing the barley diversity and evolution since the Neolithic.

## Keywords

Geometric morphometrics; cereals; Elliptical Fourier Transforms; varietal diversity; *Hordeum*

<footnote n="1">¹ ISEM, CNRS, EPHE, IRD, Université de Montpellier, Montpellier, France</footnote>

<footnote n="2">² Athéna, Lacamp, Roquedur, France</footnote>

<footnote n="3">³ UCA, INRAE, GDEC, Clermont-Ferrand, France</footnote>
**Introduction**

Barley (Hordeum vulgare L.) is one of the most important cereal crops in the world with a global production of 160 million tons every year (http://www.fao.org/). Barley tolerates drier, colder and poorer soils than wheat, which explains its wider geographical distribution in Eurasia [1, 2]. Today, barley grains are mainly used for both animal and human consumption while straw is dedicated to livestock, and starch is used in food production and chemical industry [3-5]. At the global scale, several hundreds of barley varieties are recorded today [6].

Two Hordeum subspecies can be distinguished based on spike morphology and number of fertile spikelets present at each node of the rachis. Six-row barley (Hordeum vulgare subsp. vulgare also called Hordeum hexastichum) has 3 fertile spikelets while in two-row barley (Hordeum vulgare subsp. distichum) only the central spikelet is fertile [1, 2]. Within these two subspecies, hulled and naked types are distinguished on the basis of the adherence or non-adherence of the protective envelopes of the caryopses. Hulled barley is the most widely grown type, mainly for animal feeding and malt production for brewing. Naked barley is more scarcely cultivated nowadays and mainly serves as a human food source [1, 4, 5, 7]. In addition, barley varieties can be divided into winter- or spring-sown varieties [1, 4] that differ mainly in their vernalization requirements [8]. Winter barley has a cycle of about 10 months. Sown in autumn, it needs vernalisation, i.e. period of low temperatures, for its development and closes its cycle before the summer droughts. The need of vernalization ensures that ear development takes place after the risk of frost damage has passed [9, 10]. In spring barleys, flowering is not inhibited (no vernalization) because this stage takes place during the good season [11]. They have a short growth cycle, and are well adapted to northern regions where winter conditions are too harsh for the cultivation of winter varieties. Spring barleys, however, usually provide lower yields than winter barleys.

Barley diversity has been partially studied using molecular markers, and the genotyping of 570 French accessions (784 SNPs spread on the whole genome) shows that the main variation is caused by the difference between the sowing season (winter vs. spring), prior to the number of spike rows (6 vs. 2) [12]. Current barley diversity reflects its past history, which is, at least partially, known through the study of archaeobotanical remains as well as recent genetic studies. Barley is one of the founder crops of the Old World food production [2]. It was domesticated from its wild progenitor Hordeum vulgare subsp. spontaneum around 10.000 years ago in the Fertile Crescent. Only the central spikelet of this ancestor is fertile and produce a kernel. Barley domestication has long been thought to be monophyletic [13], but recent genetic analyses of current varieties support the hypothesis of a polyphyletic and multiregional origin of barley domestication [14-17] and discerned European and Asian routes of barley spread based on Simple Sequence Repeat (SSR) markers [18, 19]. Nevertheless, the biogeographical history of barley can be traced mainly from the macro-remains (grain, spikelet or chaff) found in archaeological sites [2]. They give evidence of the spread of barley in Europe and Asia as early as the Neolithic.

Palaeogenetic analyses gave some insights into the early domestication and spread of barley [18, 20, 21]. However, most of the cereal macroremains are, especially in the Mediterranean area, preserved by charring which is detrimental to ancient DNA preservation and palaeogenetic studies [22-24]. Tracing the ancient history of barley from the study of macro-remains must therefore primarily rely on morphological
characteristics. The ratio of straight and twisted caryopses has long been used to identify the presence of two-row and six-row barley. In addition, the presence of naked or hulled barley types can be inferred from the general aspect of the grain and of its outer surface, naked barley grains being more roundish, especially in cross-section, and having cross-ripples on their surface [e.g. 25]. In the Near-East, morphological traits allow to document the presence of six-row hulled and naked barleys soon after the domestication of the two-row morphotype, which is closer to the wild H. vulgare subsp. spontaneum [2]. However, it is difficult to apply systematically and consistently between studies using qualitative morphological criteria. As a consequence, important issues in the history of barley spread and in the cultivation dynamics of barley types remain unclear or under debate. Hordeum types are not always distinguished in archaeobotanical studies. However, in Western Europe, it is often considered that the Early Neolithic agriculture relied mostly on naked six-row barley [e.g. 26, 27]. Two-row barley is however occasionally identified, in particular by chaff, and hulled six-row barley seems to be predominant in Italy and North-Eastern Spain [28, 29].

For these reasons, it is crucial to investigate quantitative (and therefore more objective) methods for discriminating and identifying barley types from seed morphology.

Several studies have shown the interest of morphometrics to discriminate barley types. Traditional morphometrics, based on the study of the length, width and thickness of archaeological barley grains, has been used to confirm the presence of 2-row barley during Roman times in France [30]. More recently, geometric morphometric analysis of grain's outline shape of a sample of 10 present-day varieties evidenced that 2-row and 6-row types could be differentiated even when the grains were experimentally charred [31]. This discrimination was later confirmed by [32] who also showed that within 2-row and 6-row barley, a selection of British and Scandinavian varieties could be distinguished. This last study also gave evidence that differences between barley types are independent of environmental conditions. In addition, it has been suggested that barley size was related to culinary systems and traditions in prehistoric Asia [33].

The north-western Mediterranean basin host an important barley diversity and France is the 5th largest producer of barley in the world according to FAO data available in 2020, with hundreds of very diverse varieties recorded [34, 35], providing a favorable context to investigate morphometric variation within and among barley types.

Based on the quantification of the size and shape variation of barley grains belonging to a set of 58 French varieties, the present study aims to (1) assess more globally how the morphological variability of caryopses is structured according to the different categories of barley (2-row vs six-row, hulled vs naked, varieties), (2) take into account the impact of interannual variation, (3) investigate the existence of morphological differences between winter and spring types, and (4) explore the relationship between morphometric and genomic proximities.

# Material and methods

## Morphometric data

A total of 1980 barley seeds corresponding to 58 varieties were studied (S1 Table). The varieties are divided into four taxonomical groups (Table 1): two-row naked (N = 13), two-row hulled (N = 23), six-row naked (N = 10) and six-row hulled (N = 12). The varieties are also divided into three different sowing periods:
spring (N = 24), winter (N = 31), and alternative (N = 3) without seasonal preference (Table 1). So-called alternative varieties corresponding to varieties that can be sown in spring or winter.

An accession includes grains from a single variety grown and sampled the same year in the same field. A total of 66 accessions were analysed. This includes 54 varieties collected a single year and 4 varieties repeatedly harvested on three different years (varieties 10004-CFL33 and 10024-ESV are 6-row-hulled-winter barley; varieties 12510-DLG and 12900-CHI are 2-row-hulled-spring barley). These last 12 accessions (4 varieties x 3 years sown) were used to test inter-annual grain size and shape variation.

<table><caption>**Table 1. Number of varieties studied per taxonomic and systematic groups.**</caption><thead><tr><th></th><th>Spring</th><th>Winter</th><th>Alternative</th></tr></thead><tbody><tr><td>2-row hulled</td><td>17</td><td>6</td><td>-</td></tr><tr><td>2-row naked</td><td>7</td><td>6</td><td>-</td></tr><tr><td>6-row hulled</td><td>-</td><td>12</td><td>-</td></tr><tr><td>6-row naked</td><td>-</td><td>7</td><td>3</td></tr></tbody></table>

All the seeds originated from the Biological Resources Centre (BRC) small grain cereals at INRAE in Clermont Ferrand (France) where they were cultivated under the similar growing conditions and stored under strictly controlled conditions.

Once received at the laboratory, the grains were placed for 48 hours in a freezer, followed by 24 hours at 38°C in an oven in order to avoid germination and eliminate pests. Each grain was then manually peeled to remove the husks. Only complete and undeformed grains were selected.

Because the analysis of samples of 50 grains from two varieties revealed that a sample size of 30 grains was sufficient to capture the size and shape variation of a population, 30 barley seeds were analysed per accession. This preliminary analysis was performed on variance estimation using rarefaction curves [36-38] (S2 Fig).

The grains were positioned on plasticine and photographed in their dorsal and lateral views using Olympus SZ-ET microscope and DP26 Olympus Camera. The lateral view documents the shape of the grains taken in their thickness while the ventral view documents the shape of the grains visualised with the furrow facing upwards. The combination of the two views allows to document the shape of the grains from two different points of view, which somehow allows to study the grains in 2.5D, approaching a 3D description. A centimetric scale was included in all pictures. The background of each picture was removed, and the grain was converted to a black mask using Adobe Photoshop CS6 software. The 2D outline coordinates was then extracted using R v. 4.1.1 and the package Momocs v. 1.3.0 (https://github.com/MomX/Momocs/; [39]). The outline coordinates were scaled using two landmark coordinates manually localized 1 cm apart on the original picture using Image J [40]. The length, width, thickness and centroid size of the grains were calculated using the Momocs package (S1 Fig).
## Genetic data

Genetic data were available for 51 of our varieties [12], corresponding to 784 SNPs markers covering the 7 chromosomes of the genome, and were used to compute Sokal and Michener genetic distances [41] between varieties. The genetic distance matrix was compared to the morphometric distance matrix using a Mantel test.

## Statistical analysis

Differences between groups were first explored for size, shape and form (size + shape) separately. Different levels of investigation were explored: between accessions, varieties, systematic and taxonomical groups, between sowing period, and between year of collect for 4 varieties.

For size analysis, differences in length, width, thickness and centroid size of the grains were tested using Kruskall-Wallis rank tests and visualised with boxplots. Pairwise differences were tested using Wilcoxon rank tests.

To analyse shape, outlines coordinates were centered and scaled. Subsequently, the elliptical Fourier transforms (EFT) were calculated. Outline of the grains are decomposed into a series of coefficients of trigonometric functions, the harmonics. The shape of the studied object is reconstructed using the inverse transform. In our case, the lateral view of the barley grain is described by 5 harmonics and the ventral view by 7. The number of harmonics was determined using the harmonic power criterium in the Momocs package [39]. These harmonic coefficients correspond to shape variables and are analysed first using a Principal Component Analysis (PCA) in order to reduce the dimensionality of the data, assess the overall grain shape variation and detect potential outliers. Then, differences between groups were tested using Multivariate Analyses of Variance (MANOVA). Pairwise differences were assessed using pairwise multilevel comparison (vegan & pairwiseAdonis R packages). Subsequently, discriminant analyses were used to separate groups and to accuracies presented were calculated using leave-one-out cross-validation (CVP) and their confidence intervals (MASS R package, [42]). Neighboor Joining (NJ) dissimilarity networks, based on the Mahalanobis distances, were computed. Differences in mean shapes were visualised using the MSHAPES function of the Momocs R package.

## Results

### Overall morphometric variation between varieties

Overall, barley varieties differ in their grain size (length, width, thickness, centroid sizes), shape and form of both their lateral and ventral views (all p < 2.2e-16). The centroid sizes of the ventral and lateral views appeared highly correlated with each other and with grain length (Fig 1A). Consequently, centroid sizes were not analysed further. Conversely, length, width and thickness show significative (all p-value < 2.2e-16) small correlation (Fig 1A) and were analysed separately.

Visualisation of between-varieties size variation (Fig 1B-1D) shows a strong general overlap with however some varieties showing larger or smaller values, neither appear related to the 2-row or 6-row, not to hulled or naked categories.
The between-varieties dissimilarity networks for the lateral and ventral shapes (respectively Fig 1E and 1F) appeared, at least partially, taxonomically structured. The ventral shape network (Fig 1F) is mainly structured by a 2-row vs 6-row varieties opposition with only few varieties not clustering within their categories, while the lateral shape network (Fig 1E) appeared less structured. When the lateral and ventral shapes were combined (Fig 1G), the network showed a clearer pattern, with a clear structuring, first between the 2-row and 6-row varieties, then between naked and hulled varieties within the two main groups. Only three exceptions can be noted corresponding to three 2-row varieties being clustered with the 6-row varieties.

<figure><img src="image_4.png" /><figcaption>**Fig 1. Overall morphometric variation between varieties.** Correlation between the size indices of length, width, thickness and centroid sizes (CS) of the ventral (VV) and lateral (VL) views of the grains (A). Boxplot of length (B) width (C) and thickness (D) showing variation between varieties. Dissimilarity network between varieties for the lateral (E) and ventral (F) views of the grain, and the shape of the two views combined (G). The colors differenciate the hulled (red) and naked (orange) 2-row varieties and hulled (purple) and naked (blue) 6-row varieties.</figcaption></figure>

## Interannual variability

Four varieties were sampled for three different years, but not necessarily in the same years, which limits direct comparisons. The varieties show different pattern of interannual size variation with one variety (10004-CFL33) showing no variation, two varieties (10024-ESV and 12900-CHI) showing interannual differences for all comparisons, and one (12510-DLG) showing significant differences only for grain thickness (Table 2 and Fig 2A–2C). The boxplots and dissimilarity network evidenced the interannual variation in grain size (Fig 2A–2C) and shape (Fig 2D, both views combined) to be weaker than variation between varieties.
<table><caption>Table 2. Interannual variation in size and shape for four varieties.</caption><thead><tr><th rowspan="2">Variety</th><th colspan="2">Length</th><th colspan="2">Width</th><th colspan="2">Thickness</th><th colspan="2">Lateral shape</th><th colspan="2">Ventral shape</th></tr><tr><th>χ2</th><th>p</th><th>χ2</th><th>p</th><th>χ2</th><th>p</th><th>F</th><th>p</th><th>F</th><th>p</th></tr></thead><tbody><tr><td>10004-CFL33</td><td>2.27</td><td>0.32</td><td>5.85</td><td>0.054</td><td>2.26</td><td>0.32</td><td>2.49</td><td><strong>0.041</strong></td><td>1.64</td><td>0.114</td></tr><tr><td>10024-ESV</td><td>27.94</td><td><strong>8.6e-07</strong></td><td>21.7</td><td><strong>1.9e-05</strong></td><td>33.72</td><td><strong>4.78e-08</strong></td><td>5.85</td><td><strong>0.001</strong></td><td>7.36</td><td><strong>0.001</strong></td></tr><tr><td>12510-DLG</td><td>4.45</td><td>0.11</td><td>6.68</td><td><strong>0.035</strong></td><td>11.48</td><td><strong>0.0032</strong></td><td>5.99</td><td><strong>0.001</strong></td><td>6.004</td><td><strong>0.001</strong></td></tr><tr><td>12900-CHI</td><td>13.1</td><td><strong>0.001</strong></td><td>9.08</td><td><strong>0.011</strong></td><td>8.48</td><td><strong>0.014</strong></td><td>10.09</td><td><strong>0.001</strong></td><td>13.5</td><td><strong>0.001</strong></td></tr></tbody></table>

Results of Kruskal Wallis tests for size (X2 and p-value) and MANOVA for shape (F and p-values). P-values are still significant after adjustment for multiple comparisons.

<figure><img src="image_5.png" /><figcaption>**Fig 2.** Interannual variation in grain morphometrics. Boxplot of the length (A), width (B) and thickness (C) of the grain for visualising size differences between the sampled years for four varieties. Dissimilarity network based on shape (two views combined) between the four varieties and their three sampled years (D).</figcaption></figure>

## Genetic information and grain morphometrics

Genetic and morphometric datasets revealed no correlation whether it be for the lateral (p = 0.889) and dorsal (p = 0.326) shape of the grains or their length, width, and thickness (p > 0.005).

## Morphometric differences between systematic and taxonomic groups

All pairwise comparisons of the size indices between the four main categories appeared significant (all p < 2.2e-16), except length difference between 6-row hulled and naked types, and thickness difference between the 2-row and 6-row hulled types (Table 3, Fig 3A-3C). Size measurements greatly varied between sowing seasons with hulled spring varieties showing larger measurements than winter varieties when the differences were significant (Fig 3A-3C) and winter varieties showing larger measurements than winter varieties in naked types. For 6-row naked, winter varieties showed larger measurements than varieties characterized as "alternative".
<table><caption>Table 3. Pairwise comparisons between taxa in grain length, width and thickness.</caption><thead><tr><th rowspan="2"></th><th colspan="3">2-row hulled</th><th colspan="3">2-row naked</th><th colspan="3">6-row hulled</th></tr><tr><th>Length</th><th>Width</th><th>Thickness</th><th>Length</th><th>Width</th><th>Thickness</th><th>Length</th><th>Width</th><th>Thickness</th></tr></thead><tbody><tr><td>2-row naked</td><td><strong>1.4e-05</strong></td><td><strong>8.6e-10</strong></td><td><strong>6.2e-11</strong></td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>6-row hulled</td><td><strong>&lt; 2e-16</strong></td><td><strong>&lt; 2e-16</strong></td><td>0.87</td><td><strong>0.014</strong></td><td><strong>&lt; 2e-16</strong></td><td><strong>1.5e-08</strong></td><td>-</td><td>-</td><td>-</td></tr><tr><td>6-row naked</td><td><strong>2.7e-15</strong></td><td><strong>&lt; 2e-16</strong></td><td><strong>&lt; 2e-16</strong></td><td><strong>0.0098</strong></td><td><strong>&lt; 2e-16</strong></td><td><strong>&lt; 2e-16</strong></td><td>0.549</td><td><strong>0.0023</strong></td><td><strong>2.6e-12</strong></td></tr></tbody></table>

Regarding dissimilarity networks, the same patterns were observed with the naked vs hulled dichotomy for the lateral shape (Fig 3D), and the 2-row vs 6-row dichotomy for the ventral shape (Fig 3E). The two sowing periods of the same category clustered together in both networks, with the 6-row (both naked and hulled) categories showing more differences between sowing periods than their 2-row counterparts (Fig 3D and 3E). The network pooling the shape of the two views (not shown) was highly similar to those based on the ventral shape.

<figure><img src="image_6.png" /><figcaption>**Fig 3.** Overall morphometric variation between systematic and taxonomic categories. Boxplot of length (A) width (B) and thickness (C) showing variation between categories. Dissimilarity network between categories for the lateral (D) and ventral (E) views of the grain. The colors differentiate hulled (h, orange) and naked (n, red) 2-row varieties and hulled (h, purple) and naked (n, blue) 6-row varieties.</figcaption></figure>

## Discriminating power of the different morphometric parameters

All varieties were then grouped, first according to the four initial categories (2-row/6-row, hulled/naked) to which was then added the sowing period.

For all comparisons of the four main categories (Fig 4, Table 4), shape and form performed equally, and the mean CVP obtained when the lateral and ventral views of the grains are combined always performed better than their separate analyses for comparing the taxonomic groups (Fig 4A). Results are more contrasted
regarding the sowing season for which the ventral view performs equally or better than the combined analyses (Fig 4B). The three size indices provided always much lower CVP than shape and form analyses. At best, grains are identified to the correct variety with a mean CVP of 48.1% (Confidence Interval (CI): 47.5–48.6%, shape of both views combined).

<figure><img src="image_7.png" /><figcaption>**Fig 4. Systematic and taxonomic signal in grain morphometrics.** Mean Cross-Validation Percentages (CVP) computed from width, thickness and length size indices, as well as from lateral and ventral shape and form of the grain analysed separately and combined contrasting the categories mentioned along the x-axis (A). The ‘4groups’ (4gp) comparison includes the 2-row hulled (2H) and naked (2N), and 6-row naked (6N) categories. The sowing periods were contrasted for each category separately comparing winter vs. spring varieties for the 2H and 2N, and winter vs. alternative for the 6N (B). All mean CVP values and their confidence intervals can be found in (Table 5).</figcaption></figure>

<table><caption>**Table 4. Differences between the 4 groups of barley in size and shape between.**</caption><thead><tr><th></th><th>χ<sup>2</sup></th><th>p</th><th>CVP</th></tr></thead><tbody><tr><td>Length</td><td>105.19</td><td><strong>&lt; 2.2e-16</strong></td><td>30.3% (CI: 29.2–31.9%)</td></tr><tr><td>Width</td><td>373.05</td><td><strong>&lt; 2.2e-16</strong></td><td>39.0% (CI: 37.9–40.3%)</td></tr><tr><td>Thickness</td><td>155.78</td><td><strong>&lt; 2.2e-16</strong></td><td>35.3% (CI: 34.3–36.4%)</td></tr><tr><td>Lateral shape</td><td>51.328</td><td><strong>&lt; 2.2e-16</strong></td><td>67.8% (CI: 66.3–69.2%)</td></tr><tr><td>Ventral shape</td><td>33.169</td><td><strong>&lt; 2.2e-16</strong></td><td>64.8% (CI: 63.3–66.0%)</td></tr><tr><td>Lat. & vent. shapes</td><td>36.185</td><td><strong>&lt; 2.2e-16</strong></td><td>76.7% (CI: 75.7–77.6%)</td></tr></tbody></table>

Results of the t-test of size, MANOVA for shape. Cross-Validation Percentages (CVP) with confidence intervals (CI). The CVP and CI are reported in Fig 5.

The ventral shape performed better in discriminating the 2-row vs 6-row types, while the lateral shape performed better for hulled vs naked discrimination. Visualisation of mean shape differences (Fig 5) revealed shorter and wider grains for naked compared to hulled barley in ventral view, but also an apex mismatch in lateral view. The differences between 2-row and 6-row barley are mainly in the lower part of the grain, with 2-row barley being shorter, wider and rounder than 6-row barley in ventral view. In the lateral view, there is also a shift in the apex of the grain and the 2-row barleys are more rounded in the furrow of the grain.
<figure><img src="image_8.png" /><figcaption>Fig 5. Meanshape differences between hulled and naked varieties (A) and between 2-row and 6-row varieties (B).</figcaption></figure>

Outlines indicate mean shapes for each of barley types: hulled (A) and 2-row (B) in blue, naked (A) and 6-row (B) in purple.
Warmer colours indicate the largest morphological differences between barley features.

The grains can be attributed to the 2-row or 6-row category with a mean CVP of 89.3% (CI: 88.8–89.9%, form of both views combined) and to the naked or hulled category with a mean CVP of 85.3 (CI: 84.6–86.1%, shape of both views combined).

When the four categories are considered simultaneously, the lateral and ventral views performed equally, less efficiently than their combined analysis that allow attributing the grain correctly with a mean CVP of 76.7% (CI: 75.7–77.6%, shape of both views combined).

Within each of the four main categories, the sowing period can be identified with relatively high CVP ranging from 65.1% (CI: 62.2–68.6%, shape of both views combined) for two-row hulled barleys to 77.6% (CI: 73.3–81.1%, shape of both views combined) for six-row naked barleys. The sowing season of hulled barleys was more discriminated than that of naked barleys. While both views performed relatively equally for hulled barleys, the ventral view performed better for naked barleys for discriminating the sowing seasons.

Spring, winter or alternative barley

For each of the four groups, spring varieties have longer, wider and thicker grains than winter varieties (Table 5, Fig 3). The only exception is 2-row naked, showing no length differences between spring and winter varieties.

<table><caption>Table 5. Taxonomic variation in size and shape for each of the 4 groups of barley.</caption><thead><tr><th rowspan="2"></th><th colspan="3">Length</th><th colspan="3">Width</th><th colspan="3">Thickness</th><th colspan="3">Lateral shape</th><th colspan="3">Ventral shape</th><th rowspan="2">Shape of the 2 views</th></tr><tr><th>χ2</th><th>p</th><th>CVP</th><th>χ2</th><th>p</th><th>CVP</th><th>χ2</th><th>p</th><th>CVP</th><th>F</th><th>p</th><th>CVP</th><th>F</th><th>p</th><th>CVP</th></tr></thead><tbody><tr><td>2-row hulled</td><td>15.48</td><td><strong>8.0e-5</strong></td><td>58.9%</td><td>3.63</td><td>0.056</td><td>53.9%</td><td>12.38</td><td><strong>0.0004</strong></td><td>58.7%</td><td>7.20</td><td><strong>&lt;2.2e-16</strong></td><td>66.6%</td><td>6.85</td><td><strong>&lt;2.2e-16</strong></td><td>67.7%</td><td>65.1%</td></tr><tr><td>2-row naked</td><td>0.20</td><td>0.65</td><td>47.3%</td><td>8.34</td><td><strong>0.004</strong></td><td>56.4%</td><td>5.89</td><td><strong>0.015</strong></td><td>52.8%</td><td>2.27</td><td><strong>0.002</strong></td><td>58.3%</td><td>4.31</td><td><strong>2.0e-11</strong></td><td>65.4%</td><td>66.7%</td></tr><tr><td>6-row naked</td><td>29.66</td><td><strong>5.0e-8</strong></td><td>62.8%</td><td>9.57</td><td><strong>0.002</strong></td><td>58.3%</td><td>12.88</td><td><strong>0.0003</strong></td><td>59.5%</td><td>6.07</td><td><strong>2.0e-13</strong></td><td>74.2%</td><td>7.95</td><td><strong>&lt;2.2e-16</strong></td><td>74.9%</td><td>77.6%</td></tr></tbody></table>

Results of Kruskal Wallis tests for size (X2, p-value and CVP) and MANOVA for shape (F, p-values and CVP).
# Discussion

## Morphometric methods

Geometric morphometric analysis of barley grains belonging to 58 French varieties show that 6-row/2-row, hulled/naked types and even varieties differ in their grain size, shape and form. As expected, form and shape variables (lateral, ventral outlines and the combination of the two) are more discriminant than size variables (length, width, thickness, lateral and ventral centroid size).

Barley varieties differ in the size, shape and form of the grains. At best, 48.1% of the grains are correctly re-assigned to the correct variety using a combination of lateral and ventral shapes. Our results are coherent to previous studies, [31] obtained 53.7% of correct cross-classification but using only 10 varieties and a different methodological approach based on sliding semi-landmarks on the ventral view of the grains. In 2019, [32] obtained 61.5% of correct classification based on the study of 54 landraces using an EFT analysis of the dorsal outline only. It appeared therefore possible that outline analysis, using the EFT method, better discriminate barley varieties than Procrustes approaches, though a quantitative comparison of the methods is still lacking.

The combination of the lateral and ventral grain outlines also provides a clear accuracy gain to analyses based on a single view. Here, the combination of lateral and ventral shapes increases the correct reclassification rate by about 5% compared to lateral shape alone. On the other hand, the use of a size variable in combination with lateral and ventral shapes (form) does not significantly improve the performance of the model. The better performance of the analyses based on the combination of lateral and ventral shapes is confirmed at all levels (6-rows/2-rows, hulled/naked).

A possible explanation for the non-correlation between genetics and morphometrics may be the low number of SNPs used. To find grain measurements causative SNPs, a Genome-Wide Association Study (GWAS) would have to be performed with hundreds of thousands of SNPs [e.g. 43].

## Inter-annual variation and influence of environmental conditions in grain morphology

Our results show that grain size, shape and form vary from one year to another. This interannual variability appeared specific to each variety and no common variation emerged. That being said, interannual variation in grain shape is here shown to be lower than varietal variation. This conclusion echoes with the findings of [32]. In their study, 54 accessions were sown in two different environments and they show that the effect of growing conditions and grain morphology was minor compared to genetic factors. Evidencing that grain morphology is more strongly determined by genetics than by environmental conditions is an important milestone with the prospect of applying these models established on modern accessions to identify barley types to which belong archaeological grains. The role of genetics in the control of grain morphology is supported by the identification of QTL involved in grain form determination [44, 45] that do not seem associated with environmental variability.

Environmental conditions interact with genetics to influence barley grain size and weight [46]. The environmental effects vary according to seasonality and barley cycle [47]. The conditions during pre-anthesis (flowering) period directly influence grain weight and size through determining assimilates, while
conditions during post-anthesis period influence cellular division, grain filling or deposits of starch grains. Grain width and thickness are controlled by cell division and grain length by the elongation within the developing endosperm, a process that ceases 20-25 days after initiation of flowering [47]. According to [48], barley grain length is weakly affected by the environment. On the other hand, width and thickness vary more strongly according to environmental factors [47, 49]. Our results showed differences in length, width and thickness for some varieties sown in several different years.

## Patterns of grain morphological variability in relation to spike anatomy

More broadly, the dissimilarity networks calculated on the lateral and ventral grain shapes combined indicates that our dataset is firstly structured by the number of spike rows (2-rows vs. 6 rows), then according to hulled/naked types. The discrimination of 6-row/2-row varieties is mainly determined by the ventral shape while the discrimination of hulled/naked varieties is determined by the lateral outline. This structuration is readily explainable by the constraints imposed by the anatomy of the spike on the morphology of the grains. These constraints are very different for the 6 and 2-row types. In 2-row barley only the central spikelet is fertile at each node of the rachis whereas the 6-row barley has three fertile spikelets at each node. This difference in the number of fertile grains attached to a rachis node induces a twisting of the lateral grains in the 6-row type, which will have less place to develop in the ear and thus will have a different shape from the untwisted central grains [31, 48]. A mutation on the vrs1 gene of the chromosome 2 is responsible for these differences in the fertility of lateral spikelets [50-52].

Morphological differences between hulled and naked barley grains are related to husks adherence to the grain after ripening, tight in hulled barleys and much looser in naked barley varieties [31, 48]. These differences, here quantified and visualized, are well known to archaeobotanists and especially the fact that naked barley grains are rounder than grains from hulled varieties [25]. Here we can see that the naked grains are more particularly rounded on the lateral view.

Cross-validation percentages are coherent with the observed structuration, with the accuracy for discriminating 2-row from 6-row barley being higher for the ventral shape, whereas the CVP for discriminating naked from hulled barley is higher for the lateral shape. Additionally, the combination of lateral and ventral views allows a slightly better discrimination of 6-row/2-row types than hulled/naked types. Grains can be attributed to 2-row or 6-row categories with a mean CVP of 89.3%, and to naked or hulled categories with a CVP of 85.3%, and, combining the four categories, to a CVP of 76.7%. The CVP values for 6-row/2-row types are close to those obtained in previous studies by [31] (CVP = 91%) and [32] (CVP = 87.6%). In their study, [32] only included the central grains of 6-row varieties to compare them to 2-row grains. The strong discrimination they obtained show that even the morphology of central, "untwisted", grains is constrained by the pressure of lateral spikelets in 6-row varieties. Here we decided to consider the whole population of 6-row grains, "twisted" and "untwisted", in order to develop approaches that can be used to identify archaeological barley grains which cannot always be easily sorted according to "twisted" and "untwisted" categories. In this regard it is important to note that geometric morphometrics gives encouraging results to distinguish the four types 6-row-hulled, 6-row-naked, 2-row-hulled and 2-row-naked when uncharred grains are studied.
## Morphological variation and sowing period

Contrary to genotypic data, which evidence a structuration primarily dependent on the sowing period [12], grain morphometric variation appeared less structured between winter and spring varieties. This is consistent with the fact that spike anatomy does not impose different morphological constraints between spring barley and winter barley, unlike the differences between 6-row/2-row and hulled/naked types. A weak morphological structuration nevertheless exists and is strong enough to identify winter/spring types when looking at the combined lateral-ventral shapes of the grain. Within each of the 6-row/2-row vs hulled/naked groups, relatively high CVPs, from 65.1% to 77.6%, are obtained for the distinction between winter/spring types.

Several hypotheses could be suggested to explain morphological differences between winter/spring types, since size and shape differences could be linked to genetic, physiologic or environmental influence during life cycle.

A first hypothesis could be genetic with the association between grain form QTLs [11, 44, 45, 53, 54] and loci related to phenological traits. Flowering time genes are classified into at least three families [55]: photoperiod genes (e.g. Ppd-H1), vernalization genes (e.g. Vrn-H1, Vrn-H2 and Vrn-H3, sgh1, sgh2, sgh3) and earliness per se (eps) genes, the last controlling flowering independently from photoperiod and temperature (e.g. Sdw1 for semi-dwarfing genes). Barley size QTLs are reported to be associated to Ppd-H1 locus [47, 56], eps2 locus [47, 56, 57], swd1 locus, which are linked to late maturity, reduced plant height, increased tillers number and biomass production [47].

The second hypothesis to explain the differences between spring and winter barley grains morphology could be ecological and linked to resources trade-offs. The plant has limited resources for its growth until maturity, which necessarily induces trade-offs in the allocation of these resources to the different plant sinks. For example, tillers formation could limit the resources available for grain filling, which could limit the size of the grains even though their number per plant would be greater. Tillering is the production of multiple stems (tillers) starting from the initial single seedling. This ensures the formation of dense tufts and multiple ears [8]. Tillering is influenced by genetic variation [8, 58-60] but also by environmental variation during pre-anthesis phase [58, 61, 62]. Tillering should be higher for winter barley than spring barley according to [8]. Thus, spring varieties that have less tillers than winter varieties could provide more resources for seed filling, which would explain why spring grains are wider and thicker than winter ones. Another explanation to differences between spring and winter barley can be variation in environmental conditions during barley growth cycle. Several critical growth periods, as tillering, flowering, filling and ripening, are sensitive to rainfall [63], drought [64, 65], available soil water [66]. The impact of environmental variability could be increased in our dataset as in the BRC of Clermont-Ferrand varieties were not always sown in the season they were supposed to be, i.e. spring varieties were occasionally sown in December because droughts during the spring sowing period are more frequent and winters are less harsh due to climate change. Unfortunately, this data was not always recorded, it would therefore be important to further investigate this question based on a detailed record of sowing time and environmental conditions along the life-cycle, and to grow the same accession at different periods.
**Conclusion**

This study of 1980 grains from 58 modern varieties demonstrated the possibility of determining barley
characteristics (2 rows/6 rows, naked/hulled, spring/winter) using morphometric analysis of caryopses.
Despite inter-annual variability, the characteristics related to varietal differences allow varieties to be
distinguished with a cross-validation percentage of 48.1%. The higher identification percentages for the
distinction between 2 rows and 6 rows (CVP = 89.3%), between naked and hulled barley (CVP = 85.3%) and
between sowing periods (between 65.6% and 77.4%) are promising for documenting the characteristics of
archaeological barley. Indeed, barley is found charred in archaeological contexts, preventing the study of
grains taxonomy using genetics. Further studies should include charring experiment of the modern
diversity in order to build a reference collection of known characteristics (such as 2-row/6-row or
naked/hulled types) directly comparable with archaeological charred grains. I would be then possible to
explore the diachronic evolution of barley and the factors shaping its diversity over time.
Supporting information

<table><caption>S1 Table. Barley accessions information as varieties, subspecies, hulled or naked barley, sowing period, agronomic selection number of grains sampled per accession.</caption><thead><tr><th>Variety</th><th>Code</th><th>ERGE number</th><th>Harvest year</th><th>2-row / 6-row</th><th>Hulled / Naked</th><th>Winter / Spring / Alternative</th><th>Genetic data available</th></tr></thead><tbody><tr><td>ALIZE</td><td>ALI-12971</td><td>12971</td><td>2006</td><td>2-row</td><td>hulled</td><td>spring</td><td>x</td></tr><tr><td>CELESTE</td><td>CEL-34435</td><td>34435</td><td>2010</td><td>6-row</td><td> naked</td><td>alternative</td><td></td></tr><tr><td>CFL113~NUE</td><td>CFL113-10737</td><td>10737</td><td>2008</td><td>2-row</td><td> naked</td><td>spring</td><td>x</td></tr><tr><td>CFL194-3-2-3</td><td>CFL194-10750</td><td>10750</td><td>2009</td><td>2-row</td><td> naked</td><td>spring</td><td>x</td></tr><tr><td>CFL33</td><td>CFL33-10004</td><td>10004</td><td>2007</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>CFL33</td><td>CFL33-10004</td><td>10004</td><td>2008</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>CFL33</td><td>CFL33-10004</td><td>10004</td><td>2017</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>CFL56-4-4-4</td><td>CFL56-10776</td><td>10776</td><td>2012</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>CFL78-44</td><td>CFL78-10006</td><td>10006</td><td>2008</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>CFL79-560~NUE</td><td>CFL79-12105</td><td>12105</td><td>2018</td><td>2-row</td><td> naked</td><td>spring</td><td>x</td></tr><tr><td>CFL98-460</td><td>CFL98-22663</td><td>22663</td><td>2010</td><td>2-row</td><td> naked</td><td>winter</td><td>x</td></tr><tr><td>CFL98-461</td><td>CFL98-22664</td><td>22664</td><td>2006</td><td>2-row</td><td> naked</td><td>winter</td><td>x</td></tr><tr><td>CHIMERE</td><td>CHI-12900</td><td>12900</td><td>2006</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>CHIMERE</td><td>CHI-12900</td><td>12900</td><td>2018</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>CHIMERE</td><td>CHI-12900</td><td>12900</td><td>2019</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>CHOPINE</td><td>CHO-24628</td><td>24628</td><td>2006</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>CINNA</td><td>CIN-12184</td><td>12184</td><td>2010</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>CONDE</td><td>CON-10366</td><td>10366</td><td>2019</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>CYTRIS</td><td>CYT-12155</td><td>12155</td><td>2012</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>DE LONGUEVILLE</td><td>DLG-12510</td><td>12510</td><td>2006</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>DE LONGUEVILLE</td><td>DLG-12510</td><td>12510</td><td>2018</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>DE LONGUEVILLE</td><td>DLG-12510</td><td>12510</td><td>2019</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>E71 OISANS</td><td>E710-12112</td><td>12112</td><td>2011</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>EMERAUDE~moment</td><td>EME-10212</td><td>10212</td><td>2016</td><td>2-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>ESCOURGEON DE CHAMPAGNE</td><td>ESC-10027</td><td>10027</td><td>2012</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>ESCOURGEON DE VENDEE</td><td>ESV-9941</td><td>9941</td><td>2009</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>ESCOURGEON NOIR</td><td>ESCN-25601</td><td>25601</td><td>2004</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>ESCOURGEON NOUVEAU</td><td>ESN-10024</td><td>10024</td><td>2016</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>ESCOURGEON NOUVEAU</td><td>ESN-10024</td><td>10024</td><td>2017</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>ESCOURGEON NOUVEAU</td><td>ESN-10024</td><td>10024</td><td>2020</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>FAUVETTE</td><td>FAU-10381</td><td>10381</td><td>2019</td><td>2-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>HATIF DU DOUAISIS</td><td>HDOU-9939</td><td>9939</td><td>2015</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>HDE8133-1</td><td>HDE8133-10177</td><td>10177</td><td>2015</td><td>6-row</td><td> naked</td><td>winter</td><td></td></tr><tr><td>HDE-8133-2</td><td>HDE8133-25490</td><td>25490</td><td>2004</td><td>6-row</td><td> naked</td><td>winter</td><td></td></tr><tr><td>HDE-8133-8</td><td>HDE8133-25492</td><td>25492</td><td>2016</td><td>6-row</td><td> naked</td><td>winter</td><td></td></tr><tr><td>HDE-833-12</td><td>HDE833-25493</td><td>25493</td><td>2007</td><td>6-row</td><td> naked</td><td>winter</td><td></td></tr><tr><td>HEXAGONALE DE WAGNONVILLE</td><td>HEW-9735</td><td>9735</td><td>2007</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>JOHANNA</td><td>JOH-11116</td><td>11116</td><td>2006</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>KT3</td><td>KT3-11140</td><td>11140</td><td>2013</td><td>2-row</td><td> naked</td><td>spring</td><td>x</td></tr><tr><td>KTR28</td><td>KTR28-11141</td><td>11141</td><td>2013</td><td>2-row</td><td> naked</td><td>spring</td><td>x</td></tr><tr><td>KTR3R-177</td><td>KTR3R-11142</td><td>11142</td><td>2012</td><td>2-row</td><td> naked</td><td>spring</td><td>x</td></tr><tr><td>KTR4</td><td>KTR4-11143</td><td>11143</td><td>2010</td><td>2-row</td><td> naked</td><td>spring</td><td>x</td></tr><tr><td>LA.VENOURS.H11</td><td>VEH11-10459</td><td>10459</td><td>2020</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>LE PUY N1</td><td>PUY-11157</td><td>11157</td><td>2016</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>LM3904</td><td>LM3904-10180</td><td>10180</td><td>2012</td><td>2-row</td><td> naked</td><td>winter</td><td>x</td></tr><tr><td>LM-5-18</td><td>LM5-28031</td><td>28031</td><td>2009</td><td>2-row</td><td> naked</td><td>winter</td><td>x</td></tr><tr><td>LM5-3~NUE</td><td>LM5-9843</td><td>9843</td><td>2012</td><td>2-row</td><td> naked</td><td>winter</td><td>x</td></tr><tr><td>LM911</td><td>LM911-10179</td><td>10179</td><td>2006</td><td>2-row</td><td> naked</td><td>winter</td><td>x</td></tr><tr><td>N60AV-1-4-1</td><td>N60AV-13965</td><td>13965</td><td>2014</td><td>6-row</td><td> naked</td><td>alternative</td><td></td></tr><tr><td>NATASHA</td><td>NAT-12127</td><td>12127</td><td>2006</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>NCA 301</td><td>NCA301-11238</td><td>11238</td><td>2009</td><td>6-row</td><td> naked</td><td>winter</td><td>x</td></tr><tr><td>NUE DE TOULOUSE</td><td>NDT-28360</td><td>28360</td><td>2005</td><td>6-row</td><td> naked</td><td>winter</td><td>x</td></tr><tr><td>ORGE ALBERT</td><td>OALB-34414</td><td>34414</td><td>2010</td><td>6-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>ORGE DU TARN</td><td>ODTA-12589</td><td>12589</td><td>2010</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>OROSE</td><td>ORO-20050</td><td>20050</td><td>2016</td><td>2-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>OSIRIS~VERN</td><td>OSV-23033</td><td>23033</td><td>2016</td><td>2-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>PACIFIC</td><td>PAC-13512</td><td>13512</td><td>2007</td><td>2-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>PETULA</td><td>PET-10227</td><td>10227</td><td>2016</td><td>2-row</td><td> hulled</td><td>winter</td><td>x</td></tr><tr><td>RENE GUILLEMART</td><td>RGUI-11568</td><td>11568</td><td>2010</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>SYLVIANA</td><td>SYL-11636</td><td>11636</td><td>2010</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>TRAPMAR-4-1-1</td><td>TRA-13962</td><td>13962</td><td>2014</td><td>6-row</td><td> naked</td><td>alternative</td><td></td></tr><tr><td>TRESOR DU VELAY</td><td>TDV-11658</td><td>11658</td><td>2009</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>TROPHEE</td><td>TRO-12143</td><td>12143</td><td>2015</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>TURQUOISE</td><td>TUR-11685</td><td>11685</td><td>2014</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>VAHINE</td><td>VAH-12922</td><td>12922</td><td>2006</td><td>2-row</td><td> hulled</td><td>spring</td><td>x</td></tr><tr><td>VIOLETTE DE GUIMALAYE</td><td>VDG-12673</td><td>12673</td><td>2006</td><td>6-row</td><td> naked</td><td>winter</td><td>x</td></tr></tbody></table>
<figure><img src="image_9.png" /><figcaption>S1 Fig. Morphometric protocol applied to barley seeds.</figcaption></figure>

1. From photography to binary mask.

2. Size measurements and outlines coordinates calculated.

<figure><img src="image_10.png" /><figcaption>S2 Fig. Rarefaction curves of barley seeds in X and percentage of population variation recorded in Y.</figcaption></figure>

**Acknowledgments**

We gratefully thank the Biological Resource Center-INRAE Clermont Ferrand for providing the studied
grains.

**Author Contributions**

**Conceptualization:** Laurent Bouby, Vincent Bonhomme, Allowen Evin.

**Data curation:** Camille Dham.

**Formal analysis:** Angèle Jeanty, Allowen Evin.

**Funding acquisition:** Allowen Evin.

**Methodology:** Angèle Jeanty, Sarah Ivorra, Allowen Evin.

**Project administration:** Allowen Evin.
**Resources:** François Balfourier, Clément Debiton, Camille Dham.

**Software:** Allowen Evin.

**Supervision:** Laurent Bouby, Jérôme Ros, Allowen Evin.

**Validation:** Laurent Bouby, Sarah Ivorra, Jérôme Ros, Allowen Evin.

**Visualization:** Vincent Bonhomme, Sarah Ivorra, Allowen Evin.

**Writing – original draft:** Angèle Jeanty, Allowen Evin.

**Writing – review & editing:** Laurent Bouby, Vincent Bonhomme, François Balfourier, Clément Debiton, Sarah Ivorra, Jérôme Ros, Allowen Evin.

## References

1. Bothmer R von, editor. Diversity in barley: Hordeum vulgare. 1st ed. Amsterdam; Boston: Elsevier; 2003.

2. Zohary D, Hopf M, Weiss E. Domestication of plants in the Old World: the origin and spread of domesticated plants in Southwest Asia, Europe, and the Mediterranean Basin. 4th ed. Oxford: Oxford University Press; 2012.

3. Newman. A Brief History of Barley Foods. CFW. 2006 [cited 23 Jun 2022].

4. Carena MJ, editor. Cereals. New York, NY: Springer US; 2009.

5. Riehl S. Barley in Archaeology and Early History. Oxford Research Encyclopedia of Environmental Science. Oxford University Press; 2019.

6. Knupffer H, Terentyeva I, Hammer K, Kovaleva O, Sato K. Ecogeographical diversity–a Vavilovian approach. Developments in Plant Genetics and Breeding. Elsevier; 2003. pp. 53–76.

7. Lister DL, Jones H, Oliveira HR, Petrie CA, Liu X, Cockram J, et al. Barley heads east: Genetic analyses reveal routes of spread through diverse Eurasian landscapes. Sun G, editor. PLoS ONE. 2018; 13: e0196652. https://doi.org/10.1371/journal.pone.0196652 PMID: 30020920

8. del Moral MBG, del Moral LFG. Tiller production and survival in relation to grain yield in winter and spring barley. Field Crops Research. 1995; 44: 85–93. https://doi.org/10.1016/0378-4290(95)00072-0

9. Ellis RP, Russell G. Plant development and grain yield in spring and winter barley. The Journal of Agricultural Science. 1984; 102: 85–95. https://doi.org/10.1017/S0021859600041496

10. Turner A, Beales J, Faure S, Dunford RP, Laurie DA. The Pseudo-Response Regulator Ppd-H1 Provides Adaptation to Photoperiod in Barley. Science. 2005; 310: 1031–1034. https://doi.org/10.1126/science.1117619 PMID: 16284181

11. Pourkheirandish M, Komatsuda T. The Importance of Barley Genetics and Domestication in a Global Perspective. Annals of Botany. 2007; 100: 999–1008. https://doi.org/10.1093/aob/mcm139 PMID: 17761690

12. Balfourier F, Bardy L, Deloche M, Exbrayat FE-V, Luguin A, Capon S, et al. COLNATOR—Caractérisation de la collection nationale de ressources génétiques d'orge. Innovations Agronomiques. 2021; 84: 145–154.

13. Badr A, M K, Sch R, Rabey HE, Effgen S, Ibrahim HH, et al. On the Origin and Domestication History of Barley (Hordeum vulgare). Molecular Biology and Evolution. 2000; 17: 499–510. https://doi.org/10.1093/oxfordjournals.molbev.a026330 PMID: 10742042
14. Morrell PL, Clegg MT. Genetic evidence for a second domestication of barley (Hordeum vulgare) east of the Fertile Crescent. Proc Natl Acad Sci USA. 2007; 104: 3289-3294. https://doi.org/10.1073/pnas.0611377104 PMID: 17360640

15. Allaby RG. Barley domestication: the end of a central dogma? Genome Biol. 2015; 16: 176. https://doi.org/10.1186/s13059-015-0743-9 PMID: 26306561

16. Poets AM, Fang Z, Clegg MT, Morrell PL. Barley landraces are characterized by geographically heterogeneous genomic origins. Genome Biology. 2015; 16: 173. https://doi.org/10.1186/s13059-015-0712-3 PMID: 26293830

17. Pankin A, Altmuller J, Becker C, von Korff M. Targeted resequencing reveals genomic signatures of barley domestication. New Phytol. 2018; 218: 1247-1259. https://doi.org/10.1111/nph.15077 PMID: 29528492

18. Jones G, Charles MP, Jones MK, Colledge S, Leigh FJ, Lister DA, et al. DNA evidence for multiple introductions of barley into Europe following dispersed domestications in Western Asia. Antiquity. 2013; 87: 701-713. https://doi.org/10.1017/S0003598X00049401

19. Lister DL, Jones MK. Is naked barley an eastern or a western crop? The combined evidence of archaeobotany and genetics. Veget Hist Archaeobot. 2013; 22: 439-446. https://doi.org/10.1007/s00334-012-0376-9

20. Mascher M, Schuenemann VJ, Davidovich U, Marom N, Himmelbach A, Hubner S, et al. Genomic analysis of 6,000-year-old cultivated grain illuminates the domestication history of barley. Nature Genetics. 2016; 48: 1089. https://doi.org/10.1038/ng.3611 PMID: 27428749

21. Liu X, Lister DL, Zhao Z, Petrie CA, Zeng X, Jones PJ, et al. Journey to the east: Diverse routes and variable flowering times for wheat and barley en route to prehistoric China. Chiang T-Y, editor. PLoS ONE. 2017; 12: e0187405. https://doi.org/10.1371/journal.pone.0187405 PMID: 29095896

22. Palmer SA, Moore JD, Clapham AJ, Rose P, Allaby RG. Archaeogenetic Evidence of Ancient Nubian Barley Evolution from Six to Two-Row Indicates Local Adaptation. Fuller DQ, editor. PLoS ONE. 2009; 4: e6301. https://doi.org/10.1371/journal.pone.0006301 PMID: 19623249

23. Hagenblad J, Morales J, Leino MW, Rodríguez-Rodríguez AC. Farmer fidelity in the Canary Islands revealed by ancient DNA from prehistoric seeds. Journal of Archaeological Science. 2017; 78: 78-87. https://doi.org/10.1016/j.jas.2016.12.001

24. Lundstrom M, Forsberg NEG, Heimdahl J, Hagenblad J, Leino MW. Genetic analyses of Scandinavian desiccated, charred and waterlogged remains of barley (Hordeum vulgare L.). Journal of Archaeological Science: Reports. 2018; 22: 11-20. https://doi.org/10.1016/j.jasrep.2018.09.006

25. Jacomet S, Greig J. Identification of cereal remains from archaeological sites. 2006 [cited 23 Jun 2022]. https://www.semanticscholar.org/paper/Identification-of-cereal-remains-from-sites-Jacomet-Greig/8b12ab151c0322124fd6279895b1f5ecb918cad1

26. Bakels CC. The Western European Loess Belt: Agrarian History, 5300 BC—AD 1000. Springer Science & Business Media; 2009.

27. Bouby L, Marinval P, Durand F, Figueiral I, Briois F, Martzluff M, et al. Early Neolithic (ca. 5850-4500 cal BC) agricultural diffusion in the Western Mediterranean: An update of archaeobotanical data in SW France.
Biehl PF, editor. PLoS ONE. 2020; 15: e0230731. https://doi.org/10.1371/journal.pone.0230731 PMID: 32240184

28. Rottoli M, Pessina A. Neolithic agriculture in Italy: an update of archaeobotanical data with particular emphasis on northern settlements. The Origins and Spread of Domestic Plants in Southwest Asia and Europe. 2007; 141–153.

29. Antolin F. Local, intensive and diverse?: Early farmers and plant economy in the North-East of the Iberian Peninsula (5500–2300 cal BC). Barkhuis; 2016. https://doi.org/10.2307/j.ctt20p56hs

30. Bouby L. L'orge à deux rangs (Hordeum Distichum) dans l'agriculture gallo-romaine: données archéobotaniques. Bulletin analytique d'histoire romaine. 2001.

31. Ros J, Evin A, Bouby L, Ruas M-P. Geometric morphometric analysis of grain shape and the identification of two-rowed barley (Hordeum vulgare subsp. distichum L.) in southern France. Journal of Archaeological Science. 2014; 41: 568–575. https://doi.org/10.1016/j.jas.2013.09.015

32. Wallace M, Bonhomme V, Russell J, Stillman E, George TS, Ramsay L, et al. Searching for the Origins of Bere Barley: a Geometric Morphometric Approach to Cereal Landrace Recognition in Archaeology. J Archaeol Method Theory. 2019; 26: 1125–1142. https://doi.org/10.1007/s10816-018-9402-2

33. Ritchey MM, Sun Y, Motuzaite Matuzeviciute G, Shoda S, Pokharia AK, Spate M, et al. The wind that shakes the barley: the role of East Asian cuisines on barley grain size. World Archaeology. 2021; 53: 287–304. https://doi.org/10.1080/00438243.2022.2030792

34. Jestin L, Le Blanc A, Le Gouis J. Status and evolvement of the barley germplasm collections in France in 2000. 2000. https://hal.archives-ouvertes.fr/hal-01189629

35. Didier A, Bardy L, Boulat E, Koenig J, Lapierre A, Exbrayat-Vinson F, et al. The PGR Networks in France: Collaboration of users and the genetic resources centre on small grain cereals. Genetika / Russian Journal of Genetics. 2012; 16: 52–60.

36. Sanders HL. Marine Benthic Diversity: A Comparative Study. The American Naturalist. 1968; 102: 243–282.

37. Hurlbert SH. The Nonconcept of Species Diversity: A Critique and Alternative Parameters. Ecology. 1971; 52: 577–586. https://doi.org/10.2307/1934145 PMID: 28973811

38. Simberloff D. Properties of the Rarefaction Diversity Measurement. The American Naturalist. 1972; 106: 414–418. https://doi.org/10.1086/282781

39. Bonhomme V, Picq S, Gaucherel C, Claude J. Momocs: Outline Analysis Using R. Journal of Statistical Software. 2014; 56: 1–24. https://doi.org/10.18637/jss.v056.i13

40. Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012; 9: 671–675. https://doi.org/10.1038/nmeth.2089 PMID: 22930834

41. Michener CD, Sokal RR. A Quantitative Approach to a Problem in Classification. Evolution. 1957; 11: 130. https://doi.org/10.2307/2406046

42. Venables WN, Ripley BD. Modern Applied Statistics with S. New York, NY: Springer New York; 2002.

43. Sharma R, Draicchio F, Bull H, Herzig P, Maurer A, Pillen K, et al. Genome-wide association of yield traits in a nested association mapping population of barley reveals new gene diversity for future breeding. Journal of Experimental Botany. 2018; 69: 3811–3822. https://doi.org/10.1093/jxb/ery178 PMID: 29767798
44. Wang J, Wu X, Yue W, Zhao C, Yang J, Zhou M. Identification of QTL for barley grain size. PeerJ. 2021; 9: e11287. https://doi.org/10.7717/peerj.11287 PMID: 33986999

45. Wang Q, Sun G, Ren X, Du B, Cheng Y, Wang Y, et al. Dissecting the Genetic Basis of Grain Size and Weight in Barley (Hordeum vulgare L.) by QTL and Comparative Genetic Analyses. Frontiers in Plant Science. 2019; 10. Available: https://www.frontiersin.org/article/10.3389/fpls.2019.00469 PMID: 31105718

46. Vasilescu L, Petcu E, Sîrbu A, Petcu E, Lazăr C. Grain morphometry analysis of Romanian winter barley cultivars registered during 1959–2019 period. Romanian agricultural research. 2021; 9.

47. Coventry SJ, Barr AR, Eglinton JK, McDonald GK. The determinants and genome locations influencing grain weight and size in barley (Hordeum vulgare L.). Aust J Agric Res. 2003; 54: 1103. https://doi.org/10.1071/AR02194

48. Reid DA. Morphology and Anatomy of the Barley Plant. In: Rasmusson DC, editor. Agronomy Monographs. Madison, WI, USA: American Society of Agronomy, Crop Science Society of America, Soil Science Society of America; 1985. pp. 73–101.

49. Fox G, Kelly A, Poulsen D, Inkerman A, Henry R. Selecting for increased barley grain size. Journal of Cereal Science. 2006; 43: 198–208. https://doi.org/10.1016/j.jcs.2005.08.004

50. Komatsuda T, Pourkheirandish M, He C, Azhaguvel P, Kanamori H, Perovic D, et al. Six-rowed barley originated from a mutation in a homeodomain-leucine zipper I-class homeobox gene. Proceedings of the National Academy of Sciences. 2007; 104: 1424–1429. https://doi.org/10.1073/pnas.0608580104 PMID: 17220272

51. Sakuma S, Pourkheirandish M, Hensel G, Kumlehn J, Stein N, Tagiri A, et al. Divergence of expression pattern contributed to neofunctionalization of duplicated HD-Zip I transcription factor in barley. New Phytologist. 2013; 197: 939–948. https://doi.org/10.1111/nph.12068 PMID: 23293955

52. Zwirek M, Waugh R, McKim SM. Interaction between row-type genes in barley controls meristem determinacy and reveals novel routes to improved grain. New Phytol. 2019; 221: 1950–1965. https://doi.org/10.1111/nph.15548 PMID: 30339269

53. Du B, Liu L, Wang Q, Sun G, Ren X, Li C, et al. Identification of QTL underlying the leaf length and area of different leaves in barley. Sci Rep. 2019; 9: 4431. https://doi.org/10.1038/s41598-019-40703-6 PMID: 30872632

54. Walker CK, Ford R, Muñoz-Amatriain M, Panozzo JF. The detection of QTLs in barley associated with endosperm hardness, grain density, grain size and malting quality using rapid phenotyping tools. Theor Appl Genet. 2013; 126: 2533–2551. https://doi.org/10.1007/s00122-013-2153-2 PMID: 23884598

55. Wiegmann M, Maurer A, Pham A, March TJ, Al-Abdallat A, Thomas WTB, et al. Barley yield formation under abiotic stress depends on the interplay between flowering time genes and environmental cues. Sci Rep. 2019; 9: 6397. https://doi.org/10.1038/s41598-019-42673-1 PMID: 31024028

56. Laurie DA, Pratchett N, Bezant JH, Snape JW. Genetic analysis of a photoperiod response gene on the short arm of chromosome 2(2H) of Hordeum vulgare (barley). Heredity. 1994; 72: 619–627. https://doi.org/10.1038/hdy.1994.85

57. Boyd WJR, Li CD, Grime CR, Cakir M, Potipibool S, Kaveeta L, et al. Conventional and molecular genetic analysis of factors contributing to variation in the timing of heading among spring barley (Hordeum vulgare
L.) genotypes grown over a mild winter growing season. Aust J Agric Res. 2003; 54: 1277-1301.
https://doi.org/10.1071/ar03014

58. Aspinall D. The control of tillering in the barley plant. 1961; 13.

59. Cannell RQ. The tillering pattern in barley varieties: I. Production, survival and contribution to yield by component tillers. J Agric Sci. 1969; 72: 405-422. https://doi.org/10.1017/S0021859600024837

60. Hussien A, Tavakol E, Horner DS, Muñoz-Amatriain M, Muehlbauer GJ, Rossini L. Genetics of Tillering in Rice and Barley. The Plant Genome. 2014; 7. https://doi.org/10.3835/plantgenome2013.10.0032

61. Alqudah AM, Koppolu R, Wolde GM, Graner A, Schnurbusch T. The Genetic Architecture of Barley Plant Stature. Front Genet. 2016; 7. https://doi.org/10.3389/fgene.2016.00117 PMID: 27446200

62. Masoni A, Pampana S, Arduini I. Barley Response to Waterlogging Duration at Tillering. Crop Science. 2016; 56: 2722-2730. https://doi.org/10.2135/cropsci2016.02.0106

63. Fernandez-Getino AP, Santin-Montanya MI, Zambrana E, de Andres EF, Tenorio JL. The response of barley to rainfall and temperature in different tillage and crop rotation systems in semi-arid conditions: Response of barley to rainfall and temperature in semi-arid conditions. Ann Appl Biol. 2015; 166: 143-153.
https://doi.org/10.1111/aab.12172

64. Abebe TD, Bauer AM, Leon J. Morphological diversity of Ethiopian barleys (Hordeum vulgare L.) in relation to geographic regions and altitudes. Hereditas. 2010; 147: 154-164.
https://doi.org/10.1111/j.1601-5223.2010.02173.x PMID: 20887601

65. Djuric N, Kandic V, Dodig D, Zorić M, Nikolic A, Surlan-Momirovic G, et al. Grain filling parameters of two- and six-rowed barley genotypes in terminal drought conditions. Italian Journal of Agrometeorology. 2018; 2018. https://doi.org/10.19199/2018.2.2038-5625.005

66. Vahamidis P, Stefopoulou A, Kotoulas V, Voloudakis D, Dercas N, Economou G. A further insight into the environmental factors determining potential grain size in malt barley under Mediterranean conditions. European Journal of Agronomy. 2021; 122: 126184. https://doi.org/10.1016/j.eja.2020.126184