# Eco-evo-devo implications and archaeobiological perspectives of trait covariance in fruits of wild and domesticated grapevines

Vincent Bonhomme, Sandrine Picq, Sarah Ivorra, Allowen Evin, Thierry Pastor, Roberto Bacilieri, Thierry Lacombe, Isabel Figueiral, Jean-Frédéric Terral, Laurent Bouby

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

## To cite this version:

Vincent Bonhomme, Sandrine Picq, Sarah Ivorra, Allowen Evin, Thierry Pastor, et al.. Eco-evo-devo implications and archaeobiological perspectives of trait covariance in fruits of wild and domesticated grapevines. PLoS ONE, 2020, 15 (11), e0239863. [https://doi.org/10.1371/journal.pone.0239863](https://doi.org/10.1371/journal.pone.0239863)

HAL Id: hal-03050703

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

Submitted on 10 Dec 2020

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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Eco-evo-devo implications and archaeobiological perspectives of trait covariance in fruits of wild and domesticated grapevines

Vincent Bonhomme<sup>1,2,*</sup>, Sandrine Picq<sup>1,2,3*</sup>, Sarah Ivorra<sup>1,2</sup>, Allowen Evin<sup>1,2</sup>, Thierry Pastor<sup>1,2</sup>, Roberto Bacilieri<sup>4</sup>, Thierry Lacombe<sup>4</sup>, Isabel Figueiral<sup>1,2</sup>, Jean-Frédéric Terral<sup>1,2‡</sup>, Laurent Bouby<sup>1,2‡</sup>

1 ISEM, Univ Montpellier, CNRS, EPHE, IRD, Montpellier, France, 2 Équipe « Dynamique de la biodiversité, anthropo-écologie », CC065 Montpellier Cedex 5, France, 3 Laurentian Forestry Centre, Natural Resources Canada, Québec, Canada, 4 UMR AGAP, Univ Montpellier, CIRAD, INRA, Montpellier SupAgro, Equipe « Diversité, Adaptation et Amélioration de la Vigne », Montpellier, France

* These authors contributed equally to this work.

‡ These authors also contributed equally to this work.

* bonhomme.vincent@gmail.com

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

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

## OPEN ACCESS

**Citation:** Bonhomme V, Picq S, Ivorra S, Evin A, Pastor T, Bacilieri R, et al. (2020) Eco-evo-devo implications and archaeobiological perspectives of trait covariance in fruits of wild and domesticated grapevines. PLoS ONE 15(11): e0239863. https://doi.org/10.1371/journal.pone.0239863

**Editor:** Mohar Singh, National Bureau of Plant Genetic Resources, INDIA

**Received:** April 22, 2020

**Accepted:** September 15, 2020

**Published:** November 6, 2020

**Peer Review History:** PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.pone.0239863

**Copyright:** © 2020 Bonhomme et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The full tables analysed are available (and citable) as: Bonhomme V, Picq S, Ivorra S, Pastor T. BerryPip dataset. doi:10.6084/m9 figshare.12696602.v2.

## Abstract

The phenotypic changes that occurred during the domestication and diversification of grapevine are well known, particularly changes in seed morphology, but the functional causes and consequences behind these variations are poorly understood. Wild and domesticate grapes differ, among others, in the form of their pips: wild grapes produce roundish pips with short stalks and cultivated varieties have more elongated pips with longer stalks. Such variations of form are of first importance for archaeobotany since the pip form is, most often, the only remaining information in archaeological settings. This study aims to enlighten archaeobotanical record and grapevine pip development by better understanding how size and shape (co)variates between pip and berry in both wild and domesticated Vitis vinifera. The covariance of berry size, number of seeds per berry ("piposity"), pip size and pip shape were explored on 49 grapevine accessions sampled among Euro-Mediterranean traditional cultivars and wild grapevines. We show that for wild grapevine, the higher the piposity, the bigger the berry and the more elongated the pip. For both wild and domesticated grapevine, the longer is the pip, the more it has a "domesticated" shape. Consequences for archaeobotanical studies are tested and discussed, and these covariances allowed the inference of berry dimensions from archaeological pips from a Southern France Roman site. This systematic exploration sheds light on new aspects of pip-berry relationship, in both size and shape, on grapevine eco-evo-devo changes during domestication, and invites to explore further the functional ecology of grapevine pip and berry and notably the impact of cultivation practices and human selection on grapevine morphology.
**Funding:** ANR project "Vignes et vins en France du Néolithique au Moyen Âge. Approche intégrée en archéosciences" (PI: Laurent Bouby) https://anr.fr/Projet-ANR-16-CE27-0013.

**Competing interests:** The authors have declared that no competing interests exist.

# Introduction

Grapevine (*Vitis vinifera* L.) is one of the most cultivated fruit species in the world [1], and has held a central economic and cultural role since ancient times, particularly in the Mediterranean area [2, 3]. The berries of grapevine are primarily used in wine production, but can be consumed fresh or dried (i.e. table grape). The wild progenitor of grapevine, *Vitis vinifera* subsp. *sylvestris*, was first domesticated in the South Caucasian area [4], which has yielded the oldest wine making evidence [5], dated to early Neolithic period (~8000 BP). The existence of other domestication centres has also been argued [6, 7]. Since the early times of domestication, grapevine varieties (or cultivars or "cépages") of *Vitis vinifera* subsp. *vinifera* have been selected and propagated; today there are several thousand varieties, identified by ampelography (i.e. grape morphology) and molecular markers [8, 9]. *V. vinifera* subspecies differ in their reproductive biology, and other phenotypic changes following domestication include larger bunches, larger berries, higher diversity in berry shape and skin colour, and higher sugar content [10, 11].

The quantitative morphological description of archaeobotanical material has brought major insights into the intertwined relationships between humans and domesticated plants [12–20], including grapevine [21–27]. So far, molecular approaches on ancient grapevine have yielded limited information on domestication [28, 29], with the study of ancient DNA hindered by its poor preservation in charred archaeobotanical material (but see [30]).

Wild and domesticate grape seeds differ in their form (size plus shape): wild grapes produce roundish pips with short stalks and cultivated varieties produce more elongated pips with longer stalks [31]. Such form variations have been identified on archaeological grapevine pips [22, 24, 25]. Archaeological material is often charred which can cause domesticated pips to appear more similar to wild pips [32], yet experimental charring has demonstrated the robustness of identification [26, 33].

The functional causes and consequences behind the form variation of grapevine pip are poorly understood. If size, shape, taste and colour of berries are traits that have been selected by humans, pip shape was likely not a direct target of selective pressures but may possibly be affected by: the berry size; the number of pips per berry; the growing environment and cultivation practices; the domestication status and the variety for domesticated grapevine; and developmental stochasticity. For instance, previous works suggested that pip size and the number of pips per berry are positively correlated to berry size [8, 34, 35] particularly for wild grapevines [24].

To what extent the form changes observed in archaeological pips imply changes in the form of berries? How this could be affected by the cultivation and domestication of wild individuals?

This paper scrutinizes how the form of berries and pips they contain covariate. A dataset of domesticated and wild contemporary grapevines allowed to compare patterns of covariation between wild and domestic *Vitis vinifera* compartments. This article is divided into four questions: i) how does size (co)vary between pips and berries, and depending on the number of pips?; ii) same question for shape; iii) how much pip shape depends on berry size, number of pips per berry, status, accession, and which practical consequences for archaeobotanical studies?; iv) can we infer the berry dimensions from archaeological pips?

# Results

## Preliminary analyses on modern material

The average piposity is equivalent between domesticated and wild accessions (mean±sd: domesticated = 2.01±0.891, wild = 2.1±0.968; GLM with Poisson error: df = 1468, z = 1.234,
P = 0.217 –Fig 1). The distribution of piposity, however, does differ (Fisher's exact test: P = 0.004) due to a higher proportion of 4-pips berries in wild grapevines (with 4-pips berries removed: P = 0.6272). No difference was observed between cultivated (domesticated and wild) accessions and those collected from wild (GLM with Poisson error: df = 1468, z = -0.676, P = 0.499).

## Covariation between pip and berry size in relation to the number of pips

**Wild vs. domesticated.** All berries and pips measurements were overall smaller for wild accessions compared to their domesticated counterparts (Wilcoxon one-tailed rank tests: all $P<10^{-8}$ –S1 Fig). Differences between wild and domesticated varied among the pip dimensions ($pip_{LengthStalk} > pip_{PositionChalaza} > pip_{Length} > pip_{Thickness} > pip_{Breadth}$).

Overall, the higher the piposity, the lower the contrast between domesticated and wild (Fig 2A). Pip dimensions of wild grapevines increase more substantially with piposity than their domesticated counterpart decrease. In wild grapevines, larger berries have more and larger pips (Fig 2A). No differences in berry dimensions/mass along increasing piposity were found for the cultivated grapevines, excepted between 1- and 2-pip for $pip_{Thickness}$ (Wilcoxon rank test: P = 0.006).

**Table vs. wine cultivars.** Table grapes have higher dimensions than wine varieties (Fig 2B). With increasing piposity, table varieties have bigger berries which is not the case for wine varieties. For pips, the only difference between low and high piposity were found for wine varieties and for $pip_{Breadth}$ and $pip_{Thickness}$ ($P<10^{-16}$).

**Wild grown in collection vs. wild in natura.** Wild grapevine pips and berries are bigger when in cultivation than their counterparts growing *in natura* (Fig 2C). Otherwise, trends of all measured variables are similar along increasing piposity. The berry mass ratios, relatively to wild collected *in natura*, were on average, 6.4 for wine varieties, 15.6 for table ones and 1.8 for cultivated wild.

Bivariate comparisons (S2 Fig) indicate positive correlations between all measurements. The total $pip_{Length}$ is the most consistent variable, between domesticated and wild grapevines: indeed, only the correlation with the $pip_{LengthStalk}$ show a significant interaction. Inversely, the correlations implying $pip_{LengthStalk}$ always show a significant interaction. For pips dimensions, the best correlations were found between $pip_{Length}$ and $pip_{PositionChalaza}$ (adj. $r^2$ = 0.8) among those with non-significant interactions, and between $pip_{LengthStalk}$ and $pip_{PositionChalaza}$ (adj. $r^2_{wild}$ = 0.615, adj. $r^2_{domesticated}$ = 0.717) among those with significant interactions. Compared to pips dimensions, correlations between berry dimensions were much better and the three possible interactions were all significant.

## Covariation between pip and berry shape in relation to the number of pips

The PCA shows that the first two PCs (Fig 3) gathered 69% of the total shape variation, and higher rank components levelled off (PC1 = 43.0%; PC2 = 26.2%; PC3 = 6.7%; PC4 = 3.9%), only the first two PCs were used as synthetic shape variables. Shape differences between wild and domesticated grapevines are mostly captured on PC1 yet scores on both PC1 (Wilcoxon rank tests, $P<10^{-16}$) and on PC2 ($P<10^{-16}$) were found different. PC1 represents how prominent is the stalk and how round is the pip; PC2 represents the circularity, a more global length/width ratio of pips, for the two views.

Regarding shape versus pip dimensions, $pip_{Length}$, correlated to all other measurements, is itself correlated with position on PC1. Two regressions were justified (Analysis of covariance: df = 1, F = 362.7, $P<10^{-16}$); their slope were identical (df = 1, F = 0.037, P = 0.848) but their intercept differed between wild and domesticated. These two regressions were significant yet
<figure><img src="image_5.png" /><figcaption>Fig 1. Distribution of the number of pips per berry for wild and domesticated grapevines.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0239863.g001

$r^2$ were low (wild: P<10<sup>-16</sup>, adj. $r^2$ = 0.195; domesticated: P<10<sup>-16</sup>, adj. $r^2$ = 0.240 - Fig 4A). When PC1 and PC2 are considered jointly, two regressions were not justified (P = 0.04) and the $r^2$ was lower (P = 0.04, adj. $r^2$ = 0.181 - Fig 4A). The longer the pip is, the more “domesticate” it looks, particularly in terms of stalk prominence.

As concerns shape versus piposity, the latter is associated with shape changes on PC1 between wild and domesticated both overall (see above) and within levels (Wilcoxon rank tests, all P<10<sup>-10</sup> - Fig 4B). Within domesticated accessions, differences were never significant. Within wild accessions, differences were not found between pairs of successive piposity levels but those between 1-3, 2-3 and 2-4 (all with P<10<sup>-8</sup> -not shown). For PC2, differences observed between domesticated and wild vanished for high piposity (1-pip: P<10<sup>-12</sup>; 2-pips:
<figure><img src="image_6.png" /><figcaption>**Fig 2. Comparisons of (logged) lengths and (log cubic-rooted) mass measurements.** Each row represents a different comparison: a) wild and domesticated grapevines, b) table and wine varieties for domesticated accessions only, and c) cultivated and collected from wild for wild grapevine only. For each measurement, boxplots are displayed for each piposity level. Differences are tested using multivariate analyses of covariances, and differences of $P<10^{-5}$, are indicated by stars in the facet title (interaction), on the right (overall difference) and above each piposity level (difference within a given piposity).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0239863.g002

$P<10^{-9}$; 3-pips: P = 0.016; 4-pips: P = 0.035 - Fig 4B). No differences within wild/domesticated and between successive piposities were found significant.

Mean shapes (Fig 5) illustrate these results. The mean absolute difference (MD) confirms that larger changes between extreme piposities are observed within wild grapevines (particularly for cultivated ones) and reveals that most of these changes affect the dorsal side of the pips (Fig 5).

## Pip shape and size in relation to status, accession and piposity; consequences for archaeobotanical inference

The respective contributions of berry height, accession and piposity on the shape of pips (Fig 6) show that the accession is the factor affecting the most the pip shape. The accession factor has a higher impact on domesticated grapevine than on wild, and on cultivated wild accessions than on those collected *in natura*. By contrast, its contributions for wine and table domesticated varieties were similar. Here again, piposity and berry height both affect the pip shape of wild accessions, but have a limited (piposity) and very limited (berry height) contribution for domesticated accessions.

Classification accuracies were compared using different training data and on different subsets (Fig 7). When different piposity levels were pooled, mirroring archaeobotanical admixtures, classification was very good at the status level (Fig 7A). Size + shape performed better (95%), than shape (93.7%) and size (92.5%) alone. When these models were evaluated on
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Principal component analysis on the joint matrices of Fourier coefficients obtained for the two views. The first two principal component gathered 69% of the total variance. The component of shape variation they capture are illustrated with reconstructed shapes at each extreme of their range. Colour of markers and convex hulls indicate pips from wild (green) and domesticated (blue) grapevines.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0239863.g003

piposity subsets, they all have an accuracy above 91%, except for 4-pips berries. As expected, accuracies were lower at the accession level (Fig 7B) and when piposity levels were pooled, size
<figure><img src="image_8.png" /><figcaption>**Fig 4.** a) Regressions PC1 and PC2 versus pip length for domesticated pips (blue) and cultivation (green); b) boxplots for each piposity level for domesticated pips (blue) and cultivation (green). Differences are tested using Wilcoxon rank tests, and differences of $P<10^{-5}$, are indicated by brackets above graphs.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0239863.g004

+ shape (89.8%) outperformed shape alone (81.3%) and size alone (46.3%). The same model ranking was observed on piposity subset, except for 4-pips berries. Overall, accuracies obtained were much higher than chance alone.

**Application to archaeological pips: Can we infer the dimensions of the (vanished) berry dimensions from the (recovered) pips?**

On modern material, we used the size of pips to predict berry heights and diameters. Both regressions show a significant interaction of the domestication status (berryDiameter: df = 1, F = 8369, $P<10^{-16}$; berryDiameter: df = 1, F = 7730, $P<10^{-16}$), and two regressions for berry
<figure><img src="image_9.png" /></figure>
**Fig 5.** Mean shapes calculated for all pips sampled from berry with 1 (blue) or 4 pips (red), for different subsets (in rows) and for the dorsal and lateral views (in columns). Between the two views, an index of shape differences between these two extreme piposity levels (0 = identical shapes; unit = as much differences as between wild and domesticated average shapes).

https://doi.org/10.1371/journal.pone.0239863.g005

diameter and two others for its height were obtained (Fig 8). All were significant (all P<10<sup>-16</sup>) yet the adjusted r<sup>2</sup> were low (berryDiameter adj. r<sup>2</sup><sub>wild</sub> = 0.585, adj. r<sup>2</sup><sub>domesticated</sub> = 0.491; berry-Height adj. r<sup>2</sup><sub>wild</sub> = 0.615, r<sup>2</sup><sub>domesticated</sub> = 0.511). Final models all used pipLength, pipThickness, and at least one PC. Table 1). On unlogged berry diameter and height, the relative deviations were obtained (S3 Fig). Mean relative deviation per accession for berryDiameter ranged from -12.9% to +10.3% for wild, and from -22.9% to +17.7% for domesticated; for berryHeight they ranged from -13.0% to +13.1% for wild, and from -29.4% to +29.4% for domesticated. The average predictions were all centred (on zero) ±1.6%.

Then, these four models were applied on the archaeobotanical material after being classified at the wild/domesticated level using LDA. 46 pips (22%) were classified with a posterior probability <0.8 and were filtered out. Among the remaining pips, 114 (72%) were classified as domesticated and 45 (28%) as wild. When compared to their modern analogues (Fig 9), the length of “domesticated pips” were closer to those of wine varieties than table varieties; the lengths of “wild pips” were intermediate between wild accessions collected in their habitat and those cultivated. For archaeological pips identified as domesticated, both inferred berry height and diameter were intermediate between wine and table modern varieties yet closer to wine ones. Similarly, for wild archaeological material, inferred berry height and diameter were intermediate between wild collected in their habitat and those grown in collection.

## Discussion

This study opens new fronts in our understanding of *Vitis vinifera* phenotypic changes under domestication and helps disentangle the interplay of the number of pips per berry, berry dimensions, domestication, pip shape, varietal diversity and cultivation practices in both wild and domesticated grapevines. We discuss implications for *Vitis vinifera* eco-evo-devo and perspectives for archaeobotanical studies for which a possible application is proposed.

## Patterns of covariation between the form of the pip, the form of the berry and the piposity

For grapevine, the theoretical maximum number of pips per berry is four, yet one was observed with five. Such abnormal piposity has already been reported [8, 36]. More than 70% berries had only one or two pips which is in accordance with previous publication [35]. There were no differences in piposity neither between domesticated or wild (Fig 1), nor between cultivated wild and those collected in natura.

Wild pips and berries are smaller than their domesticated counterparts; those from wine varieties are smaller than those of table varieties; and those of wild grapevines collected in natura are also smaller than those from cultivated wild individuals. This study details the effect of piposity on the pip form reported by previous studies [10, 24, 34, 35]. Among vertebrate dispersed plants, the reward (the fruit pulp mass) associated with a given seed mass is commensurate with work required to move it, and is expected to scale relatively [37]. For wild grapevine the berry and pip dimensions are thus expected to be constrained by their dispersers and by a general trade-off between pip size and number [38].

For all but wine varieties, the higher the piposity the longer the pip and the bigger the berry in which they develop (Fig 2). For these groups, it seems that more numerous pips are not
<figure><img src="image_10.png" /><figcaption>Fig 6. Relative contribution of berry height, accession and number of pips per berry (coloured bars) onto the shape of pips for different subsets.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0239863.g006

limited by space or nutrients but rather contribute the development of bigger berries. The two stages of berry development are well known [39, 40]. The first, prior to anthesis, is a period of rapid berry growth mostly due to cell division. After anthesis, berry growth is largely due to cell enlargement and it has been suggested that pip growth may also increases cell mitosis in the developing berry [41]. Auxins, cytokinins and gibberelins, upregulated shortly after fertilization in grapevine ovaries, are likely to trigger berry growth by cell expansion [35].

The absence of positive (or even negative) correlations between piposity, pip and berry dimensions for wine varieties remains unclear. For these varieties, the regulation, if any, may be at the bunch or stock scale, whether it has been selected (for example to concentrate sugars, aromas and flavours) or it is a by-product of another trait under selection. Since table varieties are larger than wine varieties, the berry dimensions of the latter cannot be argued to have reached a developmental limit.
<figure><img src="image_11.png" /><figcaption>**Fig 7.** Classification accuracy (LDA leave-one-out) at the a) status and b) accessions levels. Models are trained (and evaluated) on the admixture of pips, then evaluated on piposity subsets. Different combinations of training data are used (Fourier coefficients of shape, lengths/mass measurements, both). Lines provide a random baseline and summarise 10000 permutations: solid line correspond to mean accuracy; dashed line to the maximal values obtained.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0239863.g007

Finally, bivariate correlations concerning berry dimensions and mass are the strongest observed. This indicate robust allometries between berry size and mass, in other words that berry largely remain ellipsoid in shape, independently of their dimensions.

<figure><img src="image_12.png" /><figcaption>**Fig 8.** Regressions for berry dimensions from pip dimensions, obtained on modern material: Predicted versus actual (logged) berry dimensions at the domestication status level (all accessions). Columns are for berry diameter and height, respectively.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0239863.g008
<table><caption>Table 1. Estimates for pip lengths used to infer berry dimensions.</caption><thead><tr><th rowspan="2">log Predicted:</th><th rowspan="2">subset</th><th rowspan="2">Intercept</th><th colspan="6">Estimates of predictors</th><th rowspan="2">PC1</th><th rowspan="2">PC2</th></tr><tr><th>log pip<sub>Length</sub></th><th>log pip<sub>LengthStalk</sub></th><th>log pip<sub>PositionChalaza</sub></th><th>log pip<sub>Breadth</sub></th><th>log pip<sub>Thickness</sub></th></tr></thead><tbody><tr><td rowspan="2">berry<sub>Diameter</sub></td><td>wild</td><td>0.689</td><td>1.136</td><td>-</td><td>-</td><td>-</td><td>-0.438</td><td>-0.223</td><td>0.209</td></tr><tr><td>domesticated</td><td>1.298</td><td>0.478</td><td>-</td><td>0.149</td><td>-0.150</td><td>0.520</td><td>-</td><td>0.430</td></tr><tr><td rowspan="2">berry<sub>Height</sub></td><td>wild</td><td>0.838</td><td>0.732</td><td>-</td><td>0.216</td><td>0.171</td><td>-0.318</td><td>-0.524</td><td>0.258</td></tr><tr><td>domesticated</td><td>0.760</td><td>1.312</td><td>0.279</td><td>-</td><td>-0.783</td><td>0.584</td><td>1.021</td><td>-</td></tr></tbody></table>

Variables were all logged; so that *berry<sub>Diameter</sub>* (in mm and for wild) can be obtained with exp[0.65163 +1.19914×log(pip<sub>Length</sub>) +0.10617×log(pip<sub>PositionChalaza</sub>) -0.13263×log(pip<sub>Breadth</sub>) -0.45449×log(pip<sub>Thickness</sub>)].

https://doi.org/10.1371/journal.pone.0239863.t001

## Morphometrics and domestication as a wedge into grapevine eco-evo-devo

For grapevine and domesticated plants in general, domestication results in a change of desirable phenotypic patterns but also releases many “natural” constraints such as

<figure><img src="image_13.png" /><figcaption>Fig 9. Distribution of pip lengths (observed for all), berry height and diameter inferred from pip dimensions for the archaeological material from Sauvian–La Lesse. Rows distinguish domesticated and wild grapes since separate regressions were required.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0239863.g009
dispersion [42]. Cultivation practices such as pruning may explain why wild individuals grown in collection have bigger berries for higher piposity: the number of bunches is reduced, leading to larger pips. Cultivation also reduces growth constraints such as competition for water and light, self-supporting and climbing costs, those related to dispersers, etc.

Evidence of plastic and canalized phenotypic expression may be fuel for further eco-evo-
devo studies. The latter brings a conceptual and experimental framework that relies on
environmentally mediated regulatory systems to better understand ecological and evolutionary
changes [43]. Here, the norm of reaction of the pip size and shape, along increasing piposity
and berry dimensions, is clearly different at the three investigated levels: between wild and
domesticated, between wine and table, between wild individuals grown in collection and those
collected in their natural habitat.

Consequences for archaeological inference

Some pip lengths differences between wild and domesticated appear more “robust” to increas-
ing piposity, notably pipLengthStalk and pipPositionChalaza (Fig 2A; and S1 Fig). The use of lengths
in discriminating between wild and domesticated types has long been used, including when
archaeological material is charred [21, 32, 33]. Shape is confirmed here to refine such identifi-
cation [22, 24]. We show here that shape is also affected by piposity for wild accessions (Fig 4).
Most of these changes affect the dorsal side of wild pips, particularly when piposity is high.
There is as much as ~70% difference, between 1- and 4-pips for wild accessions, than those
observed between wild and domesticated types. Differences in extreme piposity are even larger
than this “domestication gap”, for the cultivated wild. This does not answer the question
whether past vineyards cultivated “true” wild grapevines or “weakly” domesticated forms [24,
25] but it nevertheless points out how piposity and cultivation practices may contribute to this
confusion by enhancing the continuum of pip forms.

Pip shape being largely used in archaeobotany, it was crucial to point out which factors con-
tribute to its variability and if they could preclude archaeobotanical. We show that the main
factor associated to pip differences was, by far, the accession and it was even more important
for domesticated accessions. In other words, accession effect appears stronger than domestica-
tion (Fig 6). Comparatively, berry dimensions and piposity poorly contributed to shape differ-
ences. This confirms the usefulness and robustness of shape to identify shape archetypes. It
may also indicate that domestication favours pip shape diversification whether this results
from genetical linkage or drift.

Classification accuracies at the status level were all good (Fig 7). Shape was superior to size
but, when considered jointly, the accuracies was improved. Whenever possible, both should
thus be included. Concerning the accession level, results evidenced even more clearly these
conclusions. Finally, we show that piposity is very unlikely to affect archaeobotanical identifi-
cation either at the status or accession levels.

An application on archaeological material: Inferring berry size from pip

Berry is very likely home to the most selected traits, from the beginnings of domestication to
varietal breeding and diversification times. Unfortunately, such fleshy parts are usually too
degraded, or even absent, to be quantified in archaeological material. The only route is thus
inference based on modern material. Here, multiple regressions on pip dimensions show that
berry diameter (Fig 8A) and height (Fig 8B) were centred and in the ±25% range.

Our archaeological application used material from Sauvian—La Lesse, a Roman farming
establishment involved in wine production, were an admixture of wild and domesticated type
is attested [44]. Berry dimensions inferred from pips are intermediate between the wild grow-
ing *in natura* and those cultivated (Fig 9). This may suggest that wild, or weakly domesticated,
individuals were cultivated in Roman vineyards. The berry dimensions inferred for domesti-
cated varieties were closer to modern wine varieties than to table ones. This is congruent with
the wine production attested at this period and in this region [44, 45].

**Conclusion**

The main finding of this exploration of berry and pip form covariation is that for wild grape-
vine, the higher the piposity, the bigger is the berry and the longer is the pip. For both wild and
domesticated, the longer is the pip, the more its shape looks like “domesticated”. Further stud-
ies and dedicated experimental plans will help clarify the contribution of cultivation practices
and, more generally on the interplay of genetic, functional and evolutionary changes that
occurred in *Vitis vinifera* between the pip, its reproductive unit, and the berry, its dispersal
reward and the main target of its domestication and varietal improvement.

**Material and methods**

**Statistical environment**

Statistical analyses were performed using the R 3.6.2 environment [46], the package Momocs 1.3.0 for everything morphometrics [47] and the tidyverse 1.2.1 packages for data manipulation and most graphics [48]. Alpha significance level was chosen equal to 10⁻⁵ all along analyses. This level both ensure marked differences for subsequent archaeobotanical application and an overall alpha level below 0.05 when repeated tests were done.

**Nomenclature**

*Status* designates compartment (domesticated vs. wild); *accession* designates the variety/individual for domesticated/wild grapevine; synecdochically, a domesticated/wild pip/berry refers to the accession they were collected from; *cultivation* designs whether wild individuals were grown in field collection or sampled *in natura*; *form* is used when *shape* and *size* are used in combination; “piposity” is short for “given a pip, the number of pips in the berry where it was sampled”.

**Modern and archaeological material**

The modern reference material included 49 accessions (30 domesticated and 19 wild) from
Euro-Mediterranean traditional cultivars and wild grapevines (Table 2). Fourteen wild grape-
vines were collected at ripeness in their habitat, and five were cultivated in the French central
ampelographic collection (INRA, Vassal-Montpellier Grapevine Biological Resources Center;
https://www6.montpellier.inra.fr/vassal), along with the domesticated accessions, which filters
out most of environmental confounding effects such as different soils, climate or cultivation
practices. Of the domesticated accessions, 21 were wine varieties and 9 table varieties. For each
accession, 30 normally developed berries have been haphazardly collected from a single, fully
ripe bunch.

Archaeobotanical material comes from two wells at the Roman farm of Sauvian—La Lesse,
extensively described elsewhere (US3022, US3063, US3171 and US3183 in [44]). These archae-
ological layers were dated to 2025–1725 BP based on pottery and coins. The waterlogged con-
ditions ensured very good preservation of the pips (N = 205).
<table><caption>Table 2. Accessions used in this study.</caption><thead><tr><th>id</th><th>accession</th><th>status</th><th>dest./cult.</th><th>pip_length</th><th>pip_lengthSalk</th><th>pipPositionChaza</th><th>pipThickness</th><th>berryHeight</th><th>berryDiameter</th><th>berryMass</th></tr></thead><tbody><tr><td>cAlvarB</td><td>Alvarinho</td><td>dom.</td><td>wine</td><td>5.18 ± 0.28</td><td>1.35 ± 0.19</td><td>2.66 ± 0.23</td><td>3.5 ± 0.21</td><td>2.5 ± 0.15</td><td>10.6 ± 0.81</td><td>10.93 ± 0.98</td><td>0.82 ± 0.21</td></tr><tr><td>cBarbE</td><td>Barbera</td><td>dom.</td><td>wine</td><td>6.73 ± 0.33</td><td>1.95 ± 0.14</td><td>3.58 ± 0.2</td><td>3.67 ± 0.23</td><td>2.7 ± 0.14</td><td>16.38 ± 1.07</td><td>15 ± 1.07</td><td>2.44 ± 0.49</td></tr><tr><td>cCabSaN</td><td>Cabernet-Sauvignon</td><td>dom.</td><td>wine</td><td>5.93 ± 0.34</td><td>1.78 ± 0.16</td><td>3.09 ± 0.2</td><td>3.57 ± 0.22</td><td>2.65 ± 0.17</td><td>14.05 ± 0.89</td><td>13.47 ± 0.84</td><td>1.57 ± 0.27</td></tr><tr><td>cCarigN</td><td>Carignan</td><td>dom.</td><td>wine</td><td>6.85 ± 0.26</td><td>2.14 ± 0.21</td><td>3.85 ± 0.24</td><td>3.54 ± 0.22</td><td>2.73 ± 0.09</td><td>18.1 ± 0.83</td><td>16.85 ± 0.84</td><td>3.16 ± 0.41</td></tr><tr><td>cChardB</td><td>Chardonnay</td><td>dom.</td><td>wine</td><td>5.28 ± 0.41</td><td>1.38 ± 0.2</td><td>3.06 ± 0.28</td><td>3.72 ± 0.28</td><td>2.78 ± 0.27</td><td>13.42 ± 1.24</td><td>12.9 ± 1.16</td><td>1.69 ± 0.41</td></tr><tr><td>cChevKN</td><td>Chevka</td><td>dom.</td><td>wine</td><td>5.97 ± 0.33</td><td>1.48 ± 0.13</td><td>3.05 ± 0.19</td><td>3.73 ± 0.19</td><td>2.53 ± 0.1</td><td>18.08 ± 1.17</td><td>16.51 ± 1.11</td><td>3.09 ± 0.58</td></tr><tr><td>cDebinB</td><td>Debina</td><td>dom.</td><td>wine</td><td>7.01 ± 0.48</td><td>1.73 ± 0.21</td><td>3.7 ± 0.3</td><td>4.21 ± 0.28</td><td>3.41 ± 0.19</td><td>19.61 ± 1.9</td><td>17.67 ± 1.58</td><td>3.89 ± 0.87</td></tr><tr><td>cFesAlB</td><td>Feteasca alba</td><td>dom.</td><td>wine</td><td>5.19 ± 0.23</td><td>1.4 ± 0.13</td><td>2.65 ± 0.19</td><td>3.46 ± 0.19</td><td>2.5 ± 0.13</td><td>12.49 ± 0.79</td><td>12.99 ± 0.69</td><td>1.62 ± 0.26</td></tr><tr><td>cGaidoB</td><td>Gaidouria</td><td>dom.</td><td>wine</td><td>5.55 ± 0.37</td><td>1.61 ± 0.19</td><td>2.92 ± 0.23</td><td>3.8 ± 0.23</td><td>2.94 ± 0.14</td><td>16.51 ± 1.36</td><td>17.1 ± 1.3</td><td>3.46 ± 0.73</td></tr><tr><td>cErenaN</td><td>Grenache</td><td>dom.</td><td>wine</td><td>5.32 ± 0.29</td><td>1.46 ± 0.1</td><td>2.84 ± 0.17</td><td>3.38 ± 0.2</td><td>2.45 ± 0.14</td><td>15.91 ± 0.78</td><td>15.12 ± 0.66</td><td>2.3 ± 0.27</td></tr><tr><td>cKypreN</td><td>Kypreiko</td><td>dom.</td><td>wine</td><td>6.79 ± 0.39</td><td>2.05 ± 0.23</td><td>3.81 ± 0.26</td><td>3.99 ± 0.27</td><td>3.07 ± 0.12</td><td>21.38 ± 1.35</td><td>17.98 ± 1.16</td><td>4.26 ± 0.79</td></tr><tr><td>cMavriN</td><td>Mavrud</td><td>dom.</td><td>wine</td><td>6.9 ± 0.26</td><td>2.35 ± 0.18</td><td>3.93 ± 0.2</td><td>4.3 ± 0.19</td><td>3.19 ± 0.14</td><td>14.91 ± 0.9</td><td>14.53 ± 0.84</td><td>2.04 ± 0.31</td></tr><tr><td>cMerloN</td><td>Merlot</td><td>dom.</td><td>wine</td><td>6.02 ± 0.32</td><td>1.5 ± 0.15</td><td>3.01 ± 0.24</td><td>3.77 ± 0.24</td><td>2.75 ± 0.18</td><td>12.62 ± 1.04</td><td>12.83 ± 1.16</td><td>1.41 ± 0.31</td></tr><tr><td>cMesFrB</td><td>Meslier Saint François</td><td>dom.</td><td>wine</td><td>5.51 ± 0.59</td><td>1.53 ± 0.2</td><td>3.05 ± 0.34</td><td>3.27 ± 0.31</td><td>2.61 ± 0.29</td><td>16.59 ± 1.64</td><td>15.54 ± 1.55</td><td>2.86 ± 0.77</td></tr><tr><td>cMourvèN</td><td>Mourvèdre</td><td>dom.</td><td>wine</td><td>6.17 ± 0.41</td><td>1.51 ± 0.18</td><td>3.11 ± 0.27</td><td>3.73 ± 0.28</td><td>2.89 ± 0.2</td><td>15.71 ± 1.04</td><td>15.27 ± 1.16</td><td>2.28 ± 0.49</td></tr><tr><td>cMusPGRs</td><td>Muscat à petits grains roses</td><td>dom.</td><td>wine</td><td>5.23 ± 0.3</td><td>1.33 ± 0.19</td><td>2.95 ± 0.23</td><td>3.12 ± 0.18</td><td>2.39 ± 0.13</td><td>12.93 ± 1.16</td><td>13.82 ± 0.9</td><td>1.97 ± 0.28</td></tr><tr><td>cMusPGN</td><td>Muscat noir à petits grains</td><td>dom.</td><td>wine</td><td>7.31 ± 0.29</td><td>2.27 ± 0.19</td><td>4.18 ± 0.26</td><td>4.2 ± 0.25</td><td>3.3 ± 0.16</td><td>20.07 ± 1.3</td><td>18.99 ± 1.39</td><td></td></tr><tr><td>cPinNoN</td><td>Pinot noir</td><td>dom.</td><td>wine</td><td>6.35 ± 0.31</td><td>1.77 ± 0.21</td><td>3.54 ± 0.2</td><td>3.93 ± 0.4</td><td>2.78 ± 0.19</td><td>13.87 ± 1.05</td><td>13.42 ± 1.19</td><td>1.83 ± 0.39</td></tr><tr><td>cRoussB</td><td>Roussanne</td><td>dom.</td><td>wine</td><td>6.57 ± 0.23</td><td>1.46 ± 0.15</td><td>3.39 ± 0.19</td><td>4.16 ± 0.23</td><td>3.09 ± 0.2</td><td>15.37 ± 0.77</td><td>14.91 ± 0.83</td><td>2.13 ± 0.32</td></tr><tr><td>cSauviB</td><td>Sauvignon</td><td>dom.</td><td>wine</td><td>6.53 ± 0.39</td><td>1.91 ± 0.19</td><td>3.49 ± 0.28</td><td>3.74 ± 0.18</td><td>2.68 ± 0.16</td><td>15.46 ± 1.07</td><td>13.45 ± 0.87</td><td>1.78 ± 0.3</td></tr><tr><td>cSyrahN</td><td>Syrah</td><td>dom.</td><td>wine</td><td>5.39 ± 0.35</td><td>1.64 ± 0.18</td><td>3.14 ± 0.28</td><td>3.21 ± 0.17</td><td>2.54 ± 0.15</td><td>15.37 ± 0.91</td><td>13.51 ± 0.85</td><td>1.99 ± 0.29</td></tr><tr><td>cAinBoB</td><td>Ain el Bouma</td><td>dom.</td><td>table</td><td>5.56 ± 0.41</td><td>1.9 ± 0.2</td><td>3.42 ± 0.23</td><td>3.53 ± 0.29</td><td>2.63 ± 0.14</td><td>21.11 ± 2.06</td><td>17.96 ± 1.8</td><td>4.51 ± 1.1</td></tr><tr><td>cChaBlB</td><td>Chaouch blanc</td><td>dom.</td><td>table</td><td>6.95 ± 0.26</td><td>1.74 ± 0.16</td><td>3.86 ± 0.31</td><td>4.08 ± 0.25</td><td>3.44 ± 0.21</td><td>24.67 ± 1.85</td><td>22.55 ± 1.82</td><td>7.76 ± 1.59</td></tr><tr><td>cFertTAR</td><td>Ferral tamara</td><td>dom.</td><td>table</td><td>6.79 ± 0.25</td><td>1.95 ± 0.18</td><td>3.52 ± 0.17</td><td>4.21 ± 0.18</td><td>3.29 ± 0.16</td><td>23.45 ± 1.22</td><td>19.23 ± 1.33</td><td>5.79 ± 0.83</td></tr><tr><td>cHadarB</td><td>Hadari</td><td>dom.</td><td>table</td><td>7.42 ± 0.38</td><td>1.89 ± 0.19</td><td>3.49 ± 0.28</td><td>4.2 ± 0.19</td><td>3.24 ± 0.15</td><td>21.54 ± 1.2</td><td>19.37 ± 1.03</td><td>4.73 ± 0.7</td></tr><tr><td>cHuniS</td><td>Hunisa</td><td>dom.</td><td>table</td><td>8.32 ± 0.39</td><td>2.27 ± 0.23</td><td>4.41 ± 0.27</td><td>4.94 ± 0.25</td><td>3.4 ± 0.2</td><td>30.28 ± 2.53</td><td>23.92 ± 1.68</td><td>9.87 ± 1.48</td></tr><tr><td>cKatPaRs</td><td>Kara papigi</td><td>dom.</td><td>table</td><td>6.58 ± 0.36</td><td>1.92 ± 0.2</td><td>3.21 ± 0.29</td><td>3.65 ± 0.22</td><td>2.8 ± 0.18</td><td>18.45 ± 1.67</td><td>16.29 ± 1.33</td><td>3 ± 0.73</td></tr><tr><td>cKraTZN</td><td>Kravi tzitzi</td><td>dom.</td><td>table</td><td>7 ± 0.31</td><td>1.92 ± 0.21</td><td>3.51 ± 0.22</td><td>3.96 ± 0.17</td><td>3.09 ± 0.14</td><td>28.02 ± 2.11</td><td>19.37 ± 1.43</td><td>6.43 ± 1.19</td></tr><tr><td>cPemGeR</td><td>Pembe Gemre</td><td>dom.</td><td>table</td><td>6.98 ± 0.33</td><td>2.11 ± 0.18</td><td>3.93 ± 0.26</td><td>3.84 ± 0.18</td><td>2.91 ± 0.09</td><td>16.46 ± 0.83</td><td>16.16 ± 1.02</td><td>3.07 ± 0.46</td></tr><tr><td>cSlivaN</td><td>Sliva</td><td>dom.</td><td>table</td><td>6.98 ± 0.47</td><td>1.95 ± 0.29</td><td>3.67 ± 0.35</td><td>4.02 ± 0.2</td><td>3.23 ± 0.11</td><td>23.41 ± 1.65</td><td>19.8 ± 1.74</td><td>5.59 ± 1.28</td></tr><tr><td>wCamSa4</td><td>Camp Saure</td><td>wild</td><td>in natura</td><td>5.23 ± 0.26</td><td>1.04 ± 0.15</td><td>2.46 ± 0.17</td><td>3.38 ± 0.17</td><td>2.5 ± 0.2</td><td>8.29 ± 0.83</td><td>8.15 ± 1.02</td><td></td></tr><tr><td>wChala7</td><td>Chalabre</td><td>wild</td><td>in natura</td><td>4.83 ± 0.28</td><td>0.92 ± 0.14</td><td>2.29 ± 0.17</td><td>3.59 ± 0.16</td><td>2.52 ± 0.25</td><td>7.68 ± 0.75</td><td>8.21 ± 0.92</td><td></td></tr><tr><td>wCoBab2</td><td>Col de la Babourade</td><td>wild</td><td>in natura</td><td>4.8 ± 0.56</td><td>0.91 ± 0.22</td><td>2.41 ± 0.31</td><td>3.24 ± 0.32</td><td>2.31 ± 0.27</td><td>8.26 ± 0.96</td><td>8.56 ± 0.93</td><td></td></tr><tr><td>wEscall13</td><td>L'Escal (13)</td><td>wild</td><td>in natura</td><td>5.51 ± 0.38</td><td>1.07 ± 0.18</td><td>2.84 ± 0.22</td><td>3.58 ± 0.2</td><td>2.76 ± 0.2</td><td>9.44 ± 0.97</td><td>9.14 ± 1.06</td><td></td></tr><tr><td>wEscall4BIs</td><td>L'Escal (14)</td><td>wild</td><td>in natura</td><td>4.84 ± 0.32</td><td>0.82 ± 0.12</td><td>2.24 ± 0.18</td><td>3.32 ± 0.15</td><td>2.4 ± 0.17</td><td>7.72 ± 0.81</td><td>7.81 ± 0.95</td><td>0.32 ± 0.1</td></tr><tr><td>wEscall16</td><td>L'Escal (16)</td><td>wild</td><td>in natura</td><td>5.34 ± 0.28</td><td>1.18 ± 0.14</td><td>2.69 ± 0.17</td><td>3.6 ± 0.26</td><td>2.56 ± 0.23</td><td>7.74 ± 0.7</td><td>7.76 ± 0.92</td><td>0.35 ± 0.09</td></tr><tr><td>wEscall17</td><td>L'Escal (17)</td><td>wild</td><td>in natura</td><td>4.94 ± 0.27</td><td>1.02 ± 0.16</td><td>2.38 ± 0.16</td><td>3.61 ± 0.18</td><td>2.77 ± 0.19</td><td>7.37 ± 0.68</td><td>7.41 ± 0.85</td><td>0.33 ± 0.09</td></tr><tr><td>wEscall18</td><td>L'Escal (18)</td><td>wild</td><td>in natura</td><td>5.17 ± 0.46</td><td>0.72 ± 0.16</td><td>2.34 ± 0.24</td><td>3.95 ± 0.37</td><td>2.78 ± 0.33</td><td>7.45 ± 0.78</td><td>7.84 ± 0.8</td><td></td></tr><tr><td>wEscall20</td><td>L'Escal (20)</td><td>wild</td><td>in natura</td><td>5.05 ± 0.24</td><td>0.89 ± 0.14</td><td>2.39 ± 0.17</td><td>3.48 ± 0.2</td><td>2.53 ± 0.21</td><td>8.05 ± 0.67</td><td>8.23 ± 0.87</td><td>0.37 ± 0.1</td></tr><tr><td>wCalme10</td><td>La Calmette (10)</td><td>wild</td><td>in natura</td><td>5.28 ± 0.31</td><td>0.89 ± 0.15</td><td>2.58 ± 0.23</td><td>3.7 ± 0.24</td><td>2.77 ± 0.27</td><td>7.98 ± 0.74</td><td>8.14 ± 0.82</td><td></td></tr><tr><td>wCalme11</td><td>La Calmette (11)</td><td>wild</td><td>in natura</td><td>4.93 ± 0.42</td><td>1.15 ± 0.25</td><td>2.66 ± 0.31</td><td>3.33 ± 0.24</td><td>2.59 ± 0.24</td><td>7.66 ± 1.09</td><td>7.57 ± 1.14</td><td></td></tr><tr><td>wPSL13</td><td>Pic Saint Loup (13)</td><td>wild</td><td>in natura</td><td>5.35 ± 0.29</td><td>0.98 ± 0.12</td><td>2.85 ± 0.21</td><td>3.75 ± 0.22</td><td>2.8 ± 0.3</td><td>8.75 ± 0.63</td><td>8.55 ± 0.73</td><td>0.42 ± 0.1</td></tr><tr><td>wPSLH</td><td>Pic Saint Loup (H)</td><td>wild</td><td>in natura</td><td>5.07 ± 0.44</td><td>0.75 ± 0.19</td><td>2.21 ± 0.27</td><td>3.78 ± 0.29</td><td>2.69 ± 0.26</td><td>7.67 ± 0.76</td><td>8.14 ± 0.88</td><td></td></tr><tr><td>wRivel1</td><td>Rivel</td><td>wild</td><td>in natura</td><td>5.83 ± 0.55</td><td>0.97 ± 0.27</td><td>2.77 ± 0.28</td><td>4.08 ± 0.28</td><td>3.24 ± 0.26</td><td>7.8 ± 0.74</td><td>7.88 ± 0.95</td><td>0.38 ± 0.11</td></tr></tbody></table>

(Continued)
<table><caption>Table 2. (Continued)</caption><thead><tr><th>id</th><th>accession</th><th>status</th><th>dest./cult.</th><th>pip<br>length</th><th>pip<br>lengthStalk</th><th>pip<br>PositionChalaza</th><th>pip<br>Breadth</th><th>pip<br>Thickness</th><th>berry<br>Height</th><th>berry<br>Diameter</th><th>berry<br>Mass</th></tr></thead><tbody><tr><td>cKets27</td><td>Ile de Ketch (27)</td><td>wild</td><td>cultivated</td><td>5.68 ± 0.27</td><td>1.08 ± 0.15</td><td>2.63 ± 0.2</td><td>3.86 ± 0.27</td><td>2.84 ± 0.19</td><td>8.88 ± 0.72</td><td>8.75 ± 0.71</td><td>0.59 ± 0.11</td></tr><tr><td>cPalmA</td><td>Palma</td><td>wild</td><td>cultivated</td><td>5.78 ± 0.32</td><td>1.51 ± 0.19</td><td>3.06 ± 0.2</td><td>3.72 ± 0.14</td><td>2.66 ± 0.23</td><td>11.27 ± 1.3</td><td>10.78 ± 1.44</td><td>0.91 ± 0.33</td></tr><tr><td>cPSL13</td><td>Pic Saint Loup (13)</td><td>wild</td><td>cultivated</td><td>5.8 ± 0.34</td><td>0.91 ± 0.19</td><td>2.74 ± 0.25</td><td>4.12 ± 0.2</td><td>3.14 ± 0.25</td><td>8.79 ± 0.73</td><td>8.72 ± 0.96</td><td>0.52 ± 0.14</td></tr><tr><td>cPSL5</td><td>Pic Saint Loup (5)</td><td>wild</td><td>cultivated</td><td>6.05 ± 0.17</td><td>1.14 ± 0.09</td><td>2.82 ± 0.13</td><td>3.9 ± 0.18</td><td>2.79 ± 0.13</td><td>10.67 ± 0.57</td><td>10.9 ± 0.65</td><td>0.77 ± 0.12</td></tr><tr><td>wLambrN</td><td>wLambrN</td><td>wild</td><td>cultivated</td><td>5.48 ± 0.3</td><td>1.06 ± 0.1</td><td>2.71 ± 0.2</td><td>3.85 ± 0.18</td><td>2.94 ± 0.25</td><td>9.12 ± 0.84</td><td>8.53 ± 0.96</td><td>0.54 ± 0.15</td></tr></tbody></table>

D: domesticated; W: wild; Wn: wine grape; Tb: table grape. For domesticated grapevines, names correspond to the variety names. Dimensions are reported with mean±sd and given in mm, except for *berry<sub>mass</sub>* which is expressed in g.

https://doi.org/10.1371/journal.pone.0239863.t002
## Traditional measurements

On modern material, the berry diameter (berryDiameter), height (berryHeight) and mass (berryMass) were obtained before dissection (Table 2). Mass was not available for 9 accessions that were removed from further analyses involving mass. Then, the number of pips (hereafter “piposity”) was recorded and one pip was randomly chosen. A single berry from the variety “Kravi tzitzi” was found with 5 pips and was discarded from further analyses. The final dataset thus consisted of 1469 pips (48 accessions × 30 pips + 1×29).

All pips, were photographed in dorsal and lateral views by the same operator (TP) using an Olympus SZ-ET stereomicroscope and an Olympus DP camera. On each pip, five length measurements were recorded by the same operator (LB) using ImageJ [49] Table 2, Fig 10): total length (pipLength), length of stalk (pipLengthStalk), position of the chalaza (pipPositionChalaza), breadth (pipBreadth) and thickness (pipThickness). All length measurements were log-transformed to focus on relative changes and minimize size differences; the mass was log cubic-root transformed for the same reason [24].

As preliminary analyses on modern material, differences between average piposity were tested using generalized linear model with Poisson error; differences in their distributions were tested using two-sided Fisher's exact tests.

## Testing the covariation between pip and berry size in relation to the number of pips

On modern material, three sets of differences in pips and berries measurements were tested using multivariate analysis of covariance: i) the interaction between status and piposity; ii) if the latter was significant, we also tested differences between status for a given piposity level; iii) whether the average piposity differs between status. These three possible sets of differences were tested between different subsets: domesticated and wild accessions; wine and table varieties for domesticated accessions; cultivated wild individuals and those collected *in natura*. Piposity was then discarded and sets were compared using Wilcoxon rank tests.

Bivariate comparisons were explored between the domesticated and wild accessions (discarding piposity), and tested with an analysis of covariance. When the domestication status was significant, separate regressions were tested and, if significant, the adjusted r² was obtained.

## Testing the covariation between pip and berry shapes in relation to the piposity

For pips, shape data were extracted from the dorsal and lateral outlines. 2D coordinates were extracted from photographs, centred, scaled, aligned along their longer axis and normalized for the position of their first point. These preliminary steps removed positional, size, rotation and phasing differences between outlines before elliptical Fourier transforms (EFT). The latter were performed on the dorsal and lateral views separately, and the number of harmonics was chosen to gather 99% of the total harmonic power (8 for both views). This generated 64 coefficients for each pip (2 views × 8 harmonics per view × 4 coefficients per harmonic).

To explore the overall variability of shape, a principal component analysis (PCA) was calculated on the full matrix of coefficients. The first two PCs (see Results) were used as synthetic shape variables (Fig 3). To test the effect of piposity and pip dimension on pip shape, the same approach than for length measurements using PC1 and PC2 as the response variables. To test the relation between shape and pip length (only pipLength was used), analyses of covariance first tested if separate regressions were justified. Then Wilcoxon tests were used to test for shape differences between and within piposity levels.
<figure><img src="image_14.png" /><figcaption>Fig 10. Dorsal and lateral views of a grapevine pip, here from the *Vitis vinifera* subsp. *sylvestris* wild" individual "Pic Saint-Loup 13", with indications of morphometric measurements: pip$_{length}$ (total length), pip$_{stalk}$ (length of the stalk), pip$_{chalaza}$ (position of the chalaza), pip$_{width}$ and pip$_{thickness}$.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0239863.g010

To visualize shape differences between extreme piposity levels (1 and 4), mean shapes for the dorsal and lateral views were calculated on the matrix of coefficients. These differences were quantified with the mean absolute difference (MD) between each sets of Fourier coefficients. To make these differences meaningful, they were divided by the mean difference of Fourier coefficients between cultivated and wild accessions with all piposity levels pooled. For each subset, MD was calculated as: $\frac{|\text{coefficients}_{\text{subset, 1-pip}} - \text{coefficients}_{\text{subset, 4-pips}}|}{|\text{coefficients}_{\text{domesticated, all pips}} - \text{coefficients}_{\text{wild, all pips}}|}$. For example, a MD equals to 0 indicates no difference between pips with a piposity of 1 or 4; a MD greater than unit indicates more differences relatively to differences that exist between domesticated and wild individuals.

## Pip shape and size in relation to status, accession and piposity

To quantify the respective contribution of berry dimensions, accession and piposity onto pip shape, a multivariate analysis of variance used the following model: all Fourier coefficients ~ *berryHeight* + *accessions* + *piposity* within accession. Since it was highly correlated to other
berry measurements (see Results), only $berry_{Height}$ was used to describe berry dimensions. The contribution of each variable is the ratio of its sum of squares over the total sum of squares (including residuals). Again, this is tested on the different subsets of interest (Fig 6).

Linear discriminant analyses (LDA) were used to evaluate whether piposity could preclude status and accession classification accuracies. Different combinations of predictors (sizes; shape; sizes + shape) were evaluated to benchmark their performance to classify the pips to their correct status and accession (Fig 7).

Given a combination of (status, accession) × (sizes+shape, sizes, shape), a leave-one-out cross-validation was used to assess classification accuracies, evaluated on all pips, to mirror archaeological admixtures where piposity is unknown (Fig 7). To cope with unbalanced group structures, we calculated a baseline for each subset that estimates the mean and maximum accuracy one can obtain by chance, using $10^4$ permutations (see [50]). If the accuracy observed is higher than the maximum value obtained using permutations, the LDA can be considered to perform better than chance alone, with an estimated alpha below $10^{-4}$.

## Predicting the dimensions of the archaeological berry dimensions

Separate multivariate regressions were calculated on the modern material, for berry height and diameter, using the five length measurements on pips. As predictor variable, we used length measurements (for dimensions) and the first two principal components (for shape). The difference between domesticated and wild grapevines regressions was first tested using an analysis of covariance: two regressions (one for cultivated, one for wild) were obtained for the berry height and two others for its diameter (Fig 8). These four regressions were fitted using stepwise regression with backward elimination based on the AIC, and started with full models: $berry_{Height/Diameter}$ for wild/ domesticated ~ $pip_{Length}$ + $pip_{LengthStalk}$ + $pip_{PositionChalaza}$ + $pip_{Breadth}$ + $pip_{Thickness}$ + PC1 + PC2, all but PCs were log-transformed). Then, archaeological pips were classified into domesticated or wild using an LDA trained using the same variables but of modern pips. Pips assigned to wild/ domesticated with a posterior probability <0.8 were filtered out. Finally, the berry height and diameter of this archaeological material were inferred using the corresponding models (Fig 9).

## Supporting information

**S1 Fig. Comparisons of (logged) lengths and (cubic-rooted) mass measurements.** On rows are displayed different subsets: a) wild and domesticated grapevines, b) for domesticated accessions, table and wine varieties and, c) for wild accessions, those collected *in natura* and others cultivated as domesticated varieties. Different piposity levels are pooled (see Fig 3 for the detail). Differences are tested using Wilcoxon rank tests and all of them have a $P<10^{-5}$.
(TIF)

**S2 Fig. Bivariate pairwise plot between (logged) lengths and (cubic-rooted) mass measurements.** For the sake of readability, only the wild versus domesticated status are displayed using different colours (green for wild; blue for domesticated). If two regressions are justified, then they are shown using the corresponding colours; otherwise a single regression line is showed in black. Then, for each regression, the correlations are tested and, if significant, the adjusted $R^2$ is displayed on the regression lines.
(TIF)

**S3 Fig. Predictions obtained from regressions for berry dimensions from pip dimensions on modern material.** The relative deviation, at the accession level, and for unlogged measurements. Columns are for berry diameter and height, respectively.
(TIF)
**Acknowledgments**

We thank the OSU-OREME (https://oreme.org/) that helped to the constitution of the wild grape pip collection. We are grateful to the INRA Vassal-Montpellier grapevine collection (Marseillan-Plage, France) that provided all the pips from cultivated varieties and cultivated wild grapes. We warmly acknowledge Michael Wallace for his help with English.

**Author Contributions**

**Conceptualization:** Vincent Bonhomme, Sandrine Picq, Jean-Frédéric Terral, Laurent Bouby.

**Data curation:** Vincent Bonhomme, Sarah Ivorra, Thierry Pastor, Laurent Bouby.

**Formal analysis:** Vincent Bonhomme, Allowen Evin.

**Funding acquisition:** Laurent Bouby.

**Investigation:** Sandrine Picq, Isabel Figueiral, Jean-Frédéric Terral.

**Methodology:** Vincent Bonhomme, Sandrine Picq, Sarah Ivorra, Allowen Evin, Jean-Frédéric Terral, Laurent Bouby.

**Project administration:** Laurent Bouby.

**Resources:** Isabel Figueiral.

**Software:** Vincent Bonhomme.

**Supervision:** Roberto Bacilieri, Jean-Frédéric Terral, Laurent Bouby.

**Validation:** Thierry Lacombe, Jean-Frédéric Terral.

**Visualization:** Vincent Bonhomme.

**Writing – original draft:** Vincent Bonhomme.

**Writing – review & editing:** Vincent Bonhomme, Sandrine Picq, Sarah Ivorra, Allowen Evin, Roberto Bacilieri, Thierry Lacombe, Jean-Frédéric Terral, Laurent Bouby.

**References**

1. FAO. FAOSTAT. 2017.

2. Brun J-P. Le vin et l'huile dans la Méditerranée antique: Viticulture, oléiculture et procédés de fabrication. Editions Errance; 2003.

3. McGovern PE. Ancient Wine The Search for the Origins of Viniculture. Princeton University Press; 2007.

4. Imazio S, Maghradze D, De Lorenzis G, Bacilieri R, Laucou V, This P, et al. From the cradle of grapevine domestication: molecular overview and description of Georgian grapevine (Vitis vinifera L.) germplasm. Tree Genetics & Genomes. 2013; 9: 641–658. https://doi.org/10.1007/s11295-013-0597-9

5. McGovern PE, Jalabadze M, Batiuk S, Callahan MP, Smith KE, Hall GR, et al. Early Neolithic wine of Georgia in the South Caucasus. Proceedings of the National Academy of Sciences. 2017; 114: E10309–E10318. https://doi.org/10.1073/pnas.1714728114 PMID: 29133421

6. Grassi F, Labra M, Imazio S, Spada A, Sgorbati S, Scienza A, et al. Evidence of a secondary grapevine domestication centre detected by SSR analysis. 2003; 1315–1320. https://doi.org/10.1007/s00122-003-1321-1 PMID: 13679993

7. Arroyo-Garcia R, Ruiz-Garcia L, Bolling L, Ocete R, Lopez MA, Arnold C, et al. Multiple origins of cultivated grapevine (Vitis vinifera L. ssp. sativa) based on chloroplast DNA polymorphisms. Molecular Ecology. 2006; 15: 3707–3714. https://doi.org/10.1111/j.1365-294X.2006.03049.x PMID: 17032268

8. Galet P. Cépages et vignobles de France. Dehan, editor. Montpellier; 1988.

9. This P, Jung A, Boccacci P, Borrego J, Botta R, Costantini L, et al. Development of a standard set of microsatellite reference alleles for identification of grape cultivars. Theoretical and Applied Genetics. 2004; 1448–1458. https://doi.org/10.1007/s00122-004-1760-3 PMID: 15565426
10. Olmo H. The Origins and Ancient History of Wine. In: McGovern P, editor. The Origins and Ancient History of Wine. Amsterdam: Gordon and Breach; 1995. pp. 31–43.

11. This P, Lacombe T, Thomas M. Historical origins and genetic diversity of wine grapes. Trends in Genetics. 2006; 22: 511–519. https://doi.org/10.1016/j.tig.2006.07.008 PMID: 16872714

12. Zohary D, Hopf M. Domestication of plants in the world. The origin and spread of cultivated plants in west Asia, Europe and the Nile valley. Clarendon Press, Oxford; 1993. https://doi.org/10.1016/0042-6989(93)90127-i PMID: 8249327

13. Jones G, Valamoti S, Charles M. Early crop diversity: A ‘new’ glume wheat from northern Greece. Vegetation History and Archaeobotany. 2000; 9: 133–146. https://doi.org/10.1007/BF01299798

14. Terral J-F, Alonso N, Chatti N, Capdevila RB i, Fabre L, Fiorentino G, et al. Historical biogeography of olive domestication (Olea europaea L.) as revealed by geometrical morphometry applied to biological and archaeological material. Journal of Biogeography. 2004; 31: 63–77. https://doi.org/10.1046/j.0305-0270.2003.01019.x

15. Terral J-F, Newton C, Ivorra S, Gros-Balthazard M, de Morais CT, Picq S, et al. Insights into the historical biogeography of the date palm (Phoenix dactylifera L.) using geometric morphometry of modern and ancient seeds. Journal of Biogeography. 2012; 39: 929–941. https://doi.org/10.1111/j.1365-2699.2011.02649.x

16. Willcox G. Measuring grain size and identifying Near Eastern cereal domestication: evidence from the Euphrates valley. Journal of Archaeological Science. 2004; 31: 145–150. https://doi.org/10.1016/j.jas.2003.07.003

17. Fuller DQ. Contrasting Patterns in Crop Domestication and Domestication Rates: Recent Archaeobotanical Insights from the Old World. Annals of Botany. 2007; 100: 903–924. https://doi.org/10.1093/aob/mcm048 PMID: 17495986

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

19. Bonhomme V, Forster E, Wallace M, Stillman E, Charles M, Jones G. The first shoots of a modern morphometrics approach to the origins of agriculture. Web Ecology. 2015; 16: 1–2. https://doi.org/10.5194/we-16-1-2016

20. 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. Journal of Archaeological Method and Theory. 2018. https://doi.org/10.1007/s10816-018-9402-2

21. Mangafa M, Kotsakis K. A New Method for the Identification of Wild and Cultivated Charred Grape Seeds. 1996; 409–418.

22. Terral J-F, Tabard E, Bouby L, Ivorra S, Pastor T, Figueiral I, et al. Evolution and history of grapevine (Vitis vinifera) under domestication: new morphometric perspectives to understand seed domestication syndrome and reveal origins of ancient European cultivars. Annals of Botany. 2010; 105: 443–455. https://doi.org/10.1093/aob/mcp298 PMID: 20034966

23. Orrù M, Grillo O, Lovicu G, Venora G, Bacchetta G. Morphological characterisation of Vitis vinifera L. seeds by image analysis and comparison with archaeological remains. Vegetation History and Archaeobotany. 2013; 22: 231–242. https://doi.org/10.1007/s00334-012-0362-2

24. Bouby L, Figueiral I, Bouchette A, Rovira N, Ivorra S, Lacombe T, et al. Bioarchaeological Insights into the Process of Domestication of Grapevine (Vitis vinifera L.) during Roman Times in Southern France. Degryse P, editor. PLoS ONE. 2013; 8: e63195. https://doi.org/10.1371/journal.pone.0063195 PMID: 23690998

25. Pagnoux C, Bouby L, Ivorra S, Petit C, Valamoti S-M, Pastor T, et al. Inferring the agrobiodiversity of Vitis vinifera L. (grapevine) in ancient Greece by comparative shape analysis of archaeological and modern seeds. Vegetation History and Archaeobotany. 2015; 24: 75–84. https://doi.org/10.1007/s00334-014-0482-y

26. Ucchesu M, Orrù M, Grillo O, Venora G, Paglietti G, Ardu A, et al. Predictive method for correct identification of archaeological charred grape seeds: Support for advances in knowledge of grape domestication process. PLoS ONE. 2016; 11: 1–18. https://doi.org/10.1371/journal.pone.0149814 PMID: 26901361

27. Karasik A, Rahimi O, David M, Weiss E, Drori E. Development of a 3D seed morphological tool for grapevine variety identification, and its comparison with SSR analysis. Scientific Reports. 2018; 8: 6545. https://doi.org/10.1038/s41598-018-24738-9 PMID: 29695830

28. Wales N, Ramos Madrigal J, Cappellini E, Carmona Baez A, Samaniego Castruita JA, Romero-Navarro JA, et al. The limits and potential of paleogenomic techniques for reconstructing grapevine domestication. Journal of Archaeological Science. 2016; 72: 57–70. https://doi.org/10.1016/j.jas.2016.05.014
29. Bacilieri R, Bouby L, Figueiral I, Schaal C, Terral J-F, Breton C, et al. Potential of combining morphometry and ancient DNA information to investigate grapevine domestication. Vegetation History and Archaeobotany. 2017; 26: 345–356. https://doi.org/10.1007/s00334-016-0597-4

30. Ramos-Madrigal J, Runge AKW, Bouby L, Lacombe T, Samaniego Castruita JA, Adam-Blondon A-F, et al. Palaeogenomic insights into the origins of French grapevine diversity. Nature Plants. 2019; 5: 595–603. https://doi.org/10.1038/s41477-019-0437-5 PMID: 31182840

31. Levadoux L. Les populations sauvages et cultivées des Vitis vinifera L. Annales de l'amélioration des plantes. 1956; 6: 59–118.

32. Smith H, Jones G. Experiments on the effects of charring on cultivated grape seeds. Journal of Archaeological Science. 1990; 17: 317–327. https://doi.org/10.1016/0305-4403(90)90026-2

33. Bouby L, Bonhomme V, Ivorra S, Pastor T, Rovira N, Tillier M, et al. Back from burn out: are experimentally charred grapevine pips too distorted to be characterized using morphometrics? Archaeological and Anthropological Sciences. 2018; 10: 943–954. https://doi.org/10.1007/s12520-016-0425-x

34. Negrul A. Evolucija razmera semjan i jagod u vinograda. Izvao TSHA. 1960; 167–176.

35. Houel C, Martin-Magniette M-L, Nicolas SD, Lacombe T, Le Cunff L, Franck D, et al. Genetic variability of berry size in the grapevine (Vitis vinifera L.). Australian Journal of Grape and Wine Research. 2013; 19: 208–220. https://doi.org/10.1111/ajgw.12021

36. Bonnet A. L'hybridité chez la vigne. Etude de la graine. Annales de l'Ecole Nationale d'Agriculture de Montpellier. 1912; 2: 5–134.

37. Edwards W. Plants reward seed dispersers in proportion to their effort: The relationship between pulp mass and seed mass in vertebrate dispersed plants. Evolutionary Ecology. 2006; 20: 365–376. https://doi.org/10.1007/s10682-006-0006-z

38. Leishman MR. Does the seed size/number trade-off model determine plant community structure? An assessment of the model mechanisms and their generality. Oikos. 2001; 93: 294–302. https://doi.org/10.1034/j.1600-0706.2001.930212.x

39. Dokoozlian N. Grape berry growth and development. Raisin Production Manual. Oakland: University of California, Agricultural and Natural Resources; 2000. pp. 30–37.

40. Ristic R, Iland PG. Relationships between seed and berry development of Vitis Vinifera L. cv Shiraz: Developmental changes in seed morphology and phenolic composition. Australian Journal of Grape and Wine Research. 2005; 11: 43–58. https://doi.org/10.1111/j.1755-0238.2005.tb00278.x

41. Ojeda H, Deloire A, Carbonneau A, Ageorges A, Romieu C. Berry development of grapevines: Relations between the growth of berries and their DNA content indicate cell multiplication and enlargement. Vitis. 1999; 38: 145–150.

42. Larson G, Piperno DR, Allaby RG, Purugganan MD, Andersson L, Arroyo-Kalin M, et al. Current perspectives and the future of domestication studies. Proceedings of the National Academy of Sciences. 2014; 111: 6139–6146. https://doi.org/10.1073/pnas.1323964111 PMID: 24757054

43. Sultan SE. Development in context: the timely emergence of eco-devo. Trends in Ecology & Evolution. 2007; 22: 575–582. https://doi.org/10.1016/j.tree.2007.06.014 PMID: 18029053

44. Figueiral I, Pomarèdes H, Court-Picon M, Bouby L, Tardy C, Terral J-F. New insights into Mediterranean Gallo-Roman farming: a closer look at archaeological wells in Southern France. Archaeological and Anthropological Sciences. 2015; 7: 201–233. https://doi.org/10.1007/s12520-014-0181-8

45. Figueiral I, Bouby L, Buffat L, Petitot H, Terral JF. Archaeobotany, vine growing and wine producing in Roman Southern France: the site of Gasquinoy (Béziers, Hérault). Journal of Archaeological Science. 2010; 37: 139–149. https://doi.org/10.1016/j.jas.2009.09.024

46. R Development Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. Vienna, Austria; 2020.

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

48. Wickham H, Averick M, Bryan J, Chang W, McGowan L, François R, et al. Welcome to the Tidyverse. Journal of Open Source Software. 2019; 4: 1686. https://doi.org/10.21105/joss.01686

49. Rasband W. ImageJ. 2008.

50. Evin A, Cucchi T, Cardini A, Strand Vidarsdottir U, Larson G, Dobney K. The long and winding road: identifying pig domestication through molar size and shape. Journal of Archaeological Science. 2013; 40: 735–743. https://doi.org/10.1016/j.jas.2012.08.005