Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

A multidisciplinary approach for investigating dietary and medicinal habits of the Medieval population of Santa Severa (7th-15th centuries, Rome, Italy)

  • Angelo Gismondi ,

    Contributed equally to this work with: Angelo Gismondi, Marica Baldoni

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Laboratory of Botany, Department of Biology, University of Rome “Tor Vergata”, Rome, Italy

  • Marica Baldoni ,

    Contributed equally to this work with: Angelo Gismondi, Marica Baldoni

    Roles Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Centro di Antropologia Molecolare per lo Studio del DNA Antico, Department of Biology, University of Rome “Tor Vergata”, Rome, Italy

  • Micaela Gnes,

    Roles Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Centro di Antropologia Molecolare per lo Studio del DNA Antico, Department of Biology, University of Rome “Tor Vergata”, Rome, Italy

  • Gabriele Scorrano,

    Roles Data curation, Formal analysis, Investigation, Methodology, Writing – review & editing

    Current address: Section for Evolutionary Genomics, GLOBE Institute, University of Copenhagen, Copenhagen K, Denmark.

    Affiliation Centro di Antropologia Molecolare per lo Studio del DNA Antico, Department of Biology, University of Rome “Tor Vergata”, Rome, Italy

  • Alessia D’Agostino,

    Roles Data curation, Formal analysis, Investigation, Methodology, Writing – review & editing

    Affiliation Laboratory of Botany, Department of Biology, University of Rome “Tor Vergata”, Rome, Italy

  • Gabriele Di Marco,

    Roles Formal analysis, Methodology, Writing – review & editing

    Affiliation Laboratory of Botany, Department of Biology, University of Rome “Tor Vergata”, Rome, Italy

  • Giulietta Calabria,

    Roles Formal analysis, Writing – review & editing

    Affiliation Centro di Antropologia Molecolare per lo Studio del DNA Antico, Department of Biology, University of Rome “Tor Vergata”, Rome, Italy

  • Michela Petrucci,

    Roles Formal analysis, Writing – review & editing

    Affiliation Centro di Antropologia Molecolare per lo Studio del DNA Antico, Department of Biology, University of Rome “Tor Vergata”, Rome, Italy

  • Gundula Müldner,

    Roles Data curation, Formal analysis, Methodology, Writing – review & editing

    Affiliation Department of Archaeology, University of Reading, Reading, England, United Kingdom

  • Matthew Von Tersch,

    Roles Formal analysis, Writing – review & editing

    Affiliation Department of Archaeology, BioArCh, University of York, York, England, United Kingdom

  • Alessandra Nardi,

    Roles Data curation, Formal analysis, Writing – review & editing

    Affiliation Department of Mathematics, University of Rome “Tor Vergata”, Rome, Italy

  • Flavio Enei,

    Roles Data curation, Funding acquisition, Writing – review & editing

    Affiliation Museo Civico di Santa Marinella “Museo del Mare e della Navigazione Antica”, Castello di Santa Severa (Roma–Italia)

  • Antonella Canini,

    Roles Conceptualization, Resources, Writing – review & editing

    Affiliation Laboratory of Botany, Department of Biology, University of Rome “Tor Vergata”, Rome, Italy

  • Olga Rickards,

    Roles Conceptualization, Resources, Writing – review & editing

    Affiliation Centro di Antropologia Molecolare per lo Studio del DNA Antico, Department of Biology, University of Rome “Tor Vergata”, Rome, Italy

  • Michelle Alexander,

    Roles Data curation, Methodology, Writing – review & editing

    Affiliation Department of Archaeology, BioArCh, University of York, York, England, United Kingdom

  •  [ ... ],
  • Cristina Martínez-Labarga

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    cristina.martinez@uniroma2.it

    Affiliation Centro di Antropologia Molecolare per lo Studio del DNA Antico, Department of Biology, University of Rome “Tor Vergata”, Rome, Italy

  • [ view all ]
  • [ view less ]

Abstract

A multidisciplinary approach, combining stable isotope analysis from bone proteins and investigations on dental calculus using DNA analysis, light microscopy, and gas chromatography coupled with mass spectrometry, was applied to reconstruct dietary and medicinal habits of the individuals recovered in the cemetery of the Castle of Santa Severa (7th-15th centuries CE; Rome, Italy). Stable isotope analysis was performed on 120 humans, 41 faunal specimens and 8 charred seeds. Dental calculus analyses were carried out on 94 samples. Overall, isotope data indicated an omnivorous diet based on C3-terrestrial protein, although some individuals possessed carbon values indicative of C4 plant consumption. In terms of animal protein, the diet was probably based on cattle, sheep, pig and chicken products, as witnessed by the archaeozoological findings. Evidence from calculus suggested the consumption of C3 cereals, Fabaceae, Fagaceae, milk and dairy products. Secondary metabolites of herbs and wine were also detected. The detection of marine fish ancient DNA, as well as of ω3 fatty acids in calculus, hypothesized the consumption of marine foodstuffs for this coastal population, despite the lack of a clear marine isotopic signal and the presence of polyunsaturated fatty acids in plant tissues. Moreover, the knowledge of ethnopharmacological tradition and the application of medicinal plants (e.g. Punica granatum L., Ephedra sp. L.) were also identified. The detection of artemisinin, known to have antimalarial properties, led to hypothesize the presence of malaria in the area. Altogether, the combined application of microscopy and biomolecular techniques provided an innovative reconstruction of Medieval lifeways in Central Italy.

Introduction

Changes in dietary habits have played a paramount role in most human evolutionary milestones [12]. In particular, the significant socio-economic transformations that impacted on both Italy and Europe throughout the Medieval period, extensively affected food demand and dietary practice [38]. For example, the terrestrial-based subsistence typical of the earlier Middle Ages underwent important changes in the later Medieval period, during which an higher consumption of fish was generally favoured by the Christianity-related abstinence from meat in established periods of the year e.g. Fridays and during Lent [23, 718]. However, Medieval Italian food habits were known to be highly variable and tended to relate both to the social context and local environment of a particular population [7, 1922].

The present investigation focuses on human remains found in the Medieval cemetery located near the Castle of Santa Severa, on the Latium coast, 50 Km North of Rome (42°01′31″N 11°56′54″E; Central Italy). The Medieval castle was built on the settlement of Pyrgi, an Etruscan harbour of the city of Caere, today known as Cerveteri. Archaeological excavation unearthed two cemetery areas, known as “Casa del Nostromo” and “Piazza della Rocca”, in use from the 7th to the 15th century CE, as confirmed by radiocarbon analysis on four skeletal specimens performed at the Centre for Diagnostic and Dating (CEDAD) of the University of Salento [23]. The skeletal collection has undergone osteological analysis and consists of 455 individuals, comprising 188 non-adults and 267 adults, which, overall, showed a high degree of biomechanical stress related to daily activities [24]. Degenerative and infectious diseases were also documented in some individuals. Regarding dental health status, low frequency of abscesses, caries and ante-mortem tooth loss were documented. Based on the results of musculoskeletal stress markers obtained by the osteological analysis, the economy of this community was probably based on both farming and agricultural activities. The individuals were mainly buried wrapped in shrouds into simple earthen graves or in sarcophagi made up of re-used materials from pre-existing buildings. Grave goods were recovered only in some non-adult burials. The only exception is represented by the individual NS SU 321 which was buried in a sarcophagus where a cross and a stone cushion could be detected (see [24] for further details).

This research aims to provide an in-depth analysis of diet, lifeways and medicinal habits of the Medieval population of Santa Severa, through the integration of isotopic analysis of bone proteins and cut-edge archaeobotanical and molecular technologies on dental calculus, both highly informative in palaeodiet reconstruction [e.g. 2528]. Carbon isotopes are useful to distinguish the consumption of C3 versus C4 plants and the contribution of marine foods in the diet [2932]. The interpretation of the human data, is highly dependent on the local environmental context and in that perspective, the analysis of coeval faunal and plant remains, represents a fundamental requisite as the isotopic composition of foodstuffs varies both geographically and temporally [2, 1718, 33].

Dental calculus, or mineralized oral plaque, is a dense mineral matrix that adheres to tooth surfaces [27, 3435]. This deposit develops from inorganic salts, deriving from saliva, which trap organic molecules belonging to oral microbiota, foods, and non-dietary microremains (e.g. pollen grains), accidentally or intentionally inhalated/ingested during everyday activities (providing information on the surrounding environment) [27, 36]. Dental calculus tends to be well preserved in archaeological contexts [34] and is increasingly recognised as a valuable resource to investigate past human diet and phytotherapeutic practices [22, 3743]. The analysis of dental calculus was mainly performed by light microscopy (henceforth LM); only rarely gas-chromatography-mass spectrometry (henceforth GC-MS) and occasionally genetic techniques have been combined [20, 22]. This research used this combination of techniques to evaluate the role of plant species in phytotherapic practices and the extent of aquatic resources consumed by the coastal population of Santa Severa.

Materials and methods

Ethics statement

Specimen numbers of the skeletal remains analysed in the present research are reported in Tables 14 and in S1 Table. The archaeological collection is housed at Museo Civico di Santa Marinella “Museo del Mare e della Navigazione Antica”, Castello di Santa Severa (Roma–Italia; responsible Dr. Flavio Enei). No permits were required for the described study, which complied with all relevant regulations. The responsible of the remains (Dr. Flavio Enei) is co- author of the present research.

thumbnail
Table 1. Results of human and animal samples subjected to stable isotope analysis and quality indicators of protein extraction procedure were reported.

For each human sample, the examined bone and the biological information were provided. Faunal samples were classified at genus level (Sus scrofa: species; Pisces: not identified fish; Sparidae: family level). The specimens excluded from the present study were highlighted in red. For sex assessment, M indicates males, F indicates females, ND was used for not determined individuals (sex determination was impossible for the lack of diagnostic elements due to their poor state of preservation), while IND was applied to indeterminate individuals (non-adults whose sex was impossible to determine because of their sexual immaturity). GAs (generic adults) were the individuals whose age at death was impossible to determine for the lack of diagnostic elements due to their poor preservation state.

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

thumbnail
Table 2. Stable isotope data of seed samples and relative weights and quality indicators.

https://doi.org/10.1371/journal.pone.0227433.t002

thumbnail
Table 3. Results of the DNA barcode analysis on dental calculus; the symbol “+” indicated positive PCR amplification.

https://doi.org/10.1371/journal.pone.0227433.t003

thumbnail
Table 4. Plant microremains found in dental calculus samples by LM.

The amount of starch granules detected by direct microscopy observation in each sample were reported and classified per morphotype. Moreover, the total of the other microdebris identified in ancient dental calculus were indicated as follows: Oleaceae pollen grain (O), Chenopodiaceae pollen grain (C), Urticaceae pollen grain (U), Poaceae phytoliths (P), Calcium oxalate crystals (Z), Asteraceae inflorescence (A) and not determined pollen grains (ND).

https://doi.org/10.1371/journal.pone.0227433.t004

Protein extraction and isotope analysis

Adult individuals (n = 112), aged approximately 18 years and over, were selected for stable isotope analysis on the basis of their good preservation and the possibility to determine age and sex. Eight non-adults were also examined, for comparison between adult and non-adult diets. To better interpret the human isotope data, 41 adult faunal bones (6 Bos sp., 1 Bubalus sp., 1 Equus sp., 6 Ovis sp., 2 Cervus sp., 4 Gallus sp., 5 Sus sp., 1 Sus scrofa, 2 Canis sp., 1 Felis sp. and 12 fish, including 6 unidentified species—Pisces, 1 Thunnus sp., 1 Galeorhinus sp., 2 specimens of the Sparidae family, 1 Labrus sp. and 1 Sparus sp.; [44]) and 8 charred seeds (2 Triticum sp. L., 3 Hordeum vulgare L., 2 Vicia sp. L. and 1 Vicia faba L.; [45]) were also analysed. For each sample, fragments of bone (reported in Table 1) were cleaned by scraping, to remove potential external contaminants, and pulverized using pestle and mortar. Protein extraction was carried out as reported in Longin [46] with modifications [20, 47]. In brief, 0.5 g of powdered bones were demineralised in 0.6 M HCl for 2 days at 4°C on a horizontal mixer, replacing the HCl every 24 hours. Once all minerals were dissolved, samples were rinsed three times with bi-distilled water, until the pH became neutral. Then, sample gelatinization was performed in presence of 1 mM HCl, at 70°C, for 48 hours. The liquid fraction, containing gelatinized proteins, was frozen at -80°C for 4 hours and lyophilized. To evaluate the extraction efficiency, bone proteins of a modern bovine (with a known isotopic composition) were used as reference control for all protein extractions. For the seeds, no pre-treatment was used [48], they were homogenised using a pestle and mortar prior to analysis. In order to determine carbon (δ13C) and nitrogen (δ15N) isotope ratios for each specimen, 0.8–1.2 mg of proteins or 2 mg of homogenised grain were analysed, in duplicate, by continuous flow isotope ratio mass spectrometry (CF-IRMS) at the University of Reading (UK) or by EA/IRMS in a GSL analyser coupled to a 20–22 mass spectrometer (Sercon, Crewe, UK) at the University of York (UK). Samples were run alongside internal and international reference standards at both laboratories (in-house Fish gelatine: δ13C -15.5 ± 0.1‰, δ15N 14.3 ± 0.2‰; Cane sugar IA-R006: δ13C -11.8 ± 0.1‰; Caffeine IAEA 600: δ13C -27.8 ± 0.1‰, δ15N 0.8 ± 0.1‰; Ammonium Sulfate IAEA N2: δ15N 20.4 ± 0.2‰). The accepted analytical error was ≤0.2‰, for both carbon and nitrogen. To test reliability and exclude contamination events by exogenous carbon and nitrogen sources, the samples were assessed against established criteria [4951]. Wilcoxon and Kruskal-Wallis tests were performed to statistically evaluate the obtained isotopic results. Average linkage cluster analysis was used to identify the presence of subgroups. Since variables with large variance tend to have a larger effect on the resulting clusters, variables were standardised before performing the analysis. All analyses were undertaken using SAS version 9.4 (SAS Institute, Cary NC).

Molecular and archaeobotanical analyses on dental calculus

Sample collection and decontamination.

Dental calculus flakes were collected from 94 individuals with an autoclaved dental pick, recovering the maximum amount of sample available per individual. Biological details of the individuals, calculus location and weight, and the analysis performed per each sample are reported in S1 Table. Dental calculus was frequently collected from multiple different teeth per individual, pooled and subsequently divided into sub-samples to be investigated. The samples were placed in sterile microcentrifuge tubes and preserved at the Department of Biology of the University of Rome ‘Tor Vergata’ (Italy). According to Crowther et al. [52], Gismondi et al. [53] and Soto et al. [54], an intensive regime of cleaning, using 5% sodium hypochlorite (NaClO), 5% sodium hydroxide (NaOH) and boiling sterilized water, was applied on surfaces, instruments and floor of all workspaces. Horizontal slide traps placed in several areas of the laboratories were monitored for the evaluation of the contaminants (S2 Table). Decontamination and sterilization protocols were conducted on the mineralised plaque, to eliminate soil particles still adhering to the external surface of the samples by a sterile acupuncture needle under a stereomicroscope (Leica ZOOM 2000, at a magnification of 30X). Each sample was treated with UV light for 10 min and immersed in 1 mL of 2% NaOH for 15 min. The pellet was then washed twice in sterile bidistilled water (40°C), rinsed in 50 μL of 100% ethanol, and left to evaporate, at 37°C, under a sterile vertical laminar flow hood (Heraeus HERAsafe HS12 Type). Before the application of the cleaning procedure, ten random dental calculus samples were selected and washed by sterile water; the latter was examined by optic microscopy. The results of these observations are reported in S3 Table. To confirm the efficacy of the proposed decontamination method, the same samples, after sterilization, were washed again by sterile water and subjected to microscopy analysis. No microremains were detected in this second washing water.

DNA extraction, amplification and sequencing.

Dental calculus was powdered using a pestle and mortar. To minimize contamination, the aDNA analysis was carried out in the aDNA Laboratory in the Departmental Center of Molecular Anthropology for Ancient DNA Studies, University of Rome, “Tor Vergata”, in Villa Mondragone, Monte Porzio Catone, Rome (http://www.bio.uniroma2.it/biologia/laboratori/lab-antropologia/DNAantico/DNA_antico/Facilities.htm) [55]. At least two independent DNA extraction were performed on 52 samples, following the protocol modified from Warinner and collaborators [56], as a minimum of 50 mg of tartar was necessary for this analysis. All criteria and precautions for ancient DNA (aDNA) study were applied [55, 5759]. Negative control extraction and amplifications were performed. To each sample, 600 μL of extraction buffer (100 mM Tris/HCl pH 8, 100 mM NaCl, 10 mM EDTA pH 8 and 2% SDS) and 50 μl of proteinase K (20 mg/mL) were added. Samples were kept at a temperature of 56°C for 6 hours and 10 μL of proteinase K (20 mg/ml) were added every two hours. After that, they were kept at 37°C overnight. After centrifugation at 13000 rpm for 5 minutes, 500 μL of phenol/chloroform/isoamyl alcohol (25:24:1) were added to the liquid phase, followed by another centrifugation. The supernatant was then purified using QIAquick PCR purification kit following the manufacturer’s procedure. Finally, DNA was eluted in 50 μL of EB Buffer. For each investigated region (Ovis aries, Gallus gallus, Bos taurus and Sus scrofa domesticus) at least two independent amplifications by PCR (Polymerase Chain Reaction) were carried out using species-specific primers and barcode primers for marine fish [60]. For animals, amplified sequences represented portions of COXI or 12S/16S rRNA genes of mitochondrial DNA (mtDNA) [61]. Before amplifying aDNA, a “Spike PCR” was performed to verify the absence of polymerase inhibitors. PCR mix was prepared in the following way: 3 μL of template DNA; 2.5 μL of PCR Buffer (10X); 2.5 μL of MgCl2 (25 mM); 2.5 μL of dNTPs (10 mM); 1 μL of forward Primer (1 μM); 1 μL of reverse Primer (1 μM); 0.25 μL of BSA (10 mg/mL); 0.2 μL of Taq polymerase (5 U/μL); 12.05 μL of ddH2O; for a total volume of 25 μL. The primers used for the “Spike PCR” and aDNA amplification are specified in S4 Table. PCR products were evaluated on 1.5% agarose gel stained with GelStar and, then, positive amplifications were purified using the ExoSAP method. After labelling with fluorescent nucleotides (Big Dye Terminator chemistry) and purification by a standard ethanol precipitation technique, samples were finally sequenced. The obtained sequences were compared with reference genes registered in GenBank, using BLAST software (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastn&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome).

In order to detect the correct activity of the different reagents and to check their contamination amplification were also performed on positive and negative controls.

Light microscopy analysis (LM).

In the laboratories of Botany, specifically reserved to the analysis of ancient biomolecules, a decalcification procedure was employed to extrapolate micro-remains from dental calculus. Twenty mg per each of the 94 specimen were degraded in 0.5 mL of 0.2 M hydrochloric acid (HCl), for 8 hours, and, after 3 washings with bidistilled water, the pellet was resuspended in 100 μL of bidistilled water and glycerol (1:1), under a sterile hood, and placed on glass slides to be analysed at OM (Axio Observer 7, ZEISS). The whole sample volume was analysed and each micro-debris was measured and photographed using ZEN imaging software. All microfossils were recognised on the basis of their morphology and distinctive features by direct comparison with a modern experimental collection of starches [62] and literature data [6364].

Gas-Chromatography Mass-Spectrometry analysis (GC-MS).

GC-MS analysis was carried out for the individuals presenting a minimum of 10 mg of tartar (68 individuals). This qualitative approach had already been successfully applied by the authors, both on ancient dental calculus and Roman frescoes [22, 65]. Briefly, once dissolved in 0.5 mL of 3% HCl, samples were incubated with 0.5 mL di hexane and left in agitation for two hours. Then, the supernatant fraction was recovered and dried out by a speed-vac system (Eppendorf AG 22331 Hamburg, Concentration Plus). After resuspension in 60 μL of hexane and derivatization by 40 μL of Methyl-8-Reagent (Thermo Scientific), 2 μL of extract was injected in a GC-MS QP2010 system (Shimadzu, Japan), equipped with a DB-5 capillary column (Phenomenex; length 30 m × diameter 0.25 mm × thickness 0.25μm), at the temperature of 280°C (splitless modality), in triplicate per each dental calculus. The carrier gas was helium (constant flow of 1 mL/min). Column temperature was initially set at 60°C for 5 min (initial oven temperature) and, then, increased at a rate 6°C/min up to 150°C for 5 min, then up to 250°C for 5 min, and to the final temperature of 330°C for 25 min. An electron impact of 70 eV (scanning from 100 to 700 m/z) was used for the ionization (ion source temperature 230°C, interface temperature 320°C, solvent cut time 6 min). Each chemical compound was identified by comparison of its mass spectrum with those registered in the NIST Library 14 loaded in the detection software (similarity values higher than 85%). Specific plants and food categories, ingested at least once in lifetime, were supposed combining the detected analytes with data from the literature and scientific food databases[6667]. No relevant differences among replicates were detected.

Results

Stable isotope analysis

Proteins were successfully extracted from all 120 human and 41 animal bones. All stable isotope results, including the relative quality control indicators, are reported in Table 1.

The protein samples extracted from human and faunal bones exhibited C:N atomic ratios between 3.1 and 3.3 and between 3.1 and 3.6, respectively, and sufficient yield values (for humans, ranging from 2 to 43%; for animals, ranging from 2 to 38%). Twelve humans and seven animals were discarded for this study (specimens highlighted in red in Table 1), because of the lack of consistency between replicate samples (standard deviation of over 0.2‰). The isotopic results for the eight plant seeds from the site are indicated in Table 2.

The isotopic ratios obtained for the faunal samples exhibited high variability even for specimens attributable to the same species (Fig 1).

thumbnail
Fig 1. Human, faunal and seed isotope analysis.

Bivariate plot of the Santa Severa human, faunal and seed stable isotope data; where multiple individuals of a single faunal or plant species were measured, the median and 25th and 75th percentile levels are shown.

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

The specimens placed at the extremes of the plot, deer (Cervus sp.) and unidentified fishes (Pisces), showing the lowest and the highest δ13C and δ15N values among all samples, were considered as examples of herbivorous terrestrial and marine animals, respectively. For this reason, they were used as reference values for the interpretation of the remaining animal and human samples. Overall, herbivores demonstrated δ13C and δ15N mean values of -20.2‰ (S.D. 1.3‰) and 6.3‰ (S.D. 1.3‰), respectively. These values are characteristic of animals feeding on C3 plants. The cattle (Bos sp.), sheep (Ovis sp.), horse (Equus sp.), and buffalo (Bubalus sp.) all possess δ13C and δ15N values enriched with respect to those observed for the deer (Cervus sp.). A single cat represents the only obligate carnivore: its isotopic values (δ13C = -19.2‰, S.D. 0.1‰, δ15N = 9.5‰, S.D. 0.2‰) are enriched by about +1‰ in δ13C and +3‰ in δ15N compared to herbivores, and are coherent with its relative trophic position (Fig 1). Sus sp. samples, expected to possessed isotope values supporting an omnivorous diet, showed measured ratios (mean δ13C = -20.4‰, S.D. 0.5‰, δ15N = 6.4‰, S.D. 1.1‰) very close to those observed for sheep, cattle and horse, suggesting they had more likely a herbivorous diet, with the exception of the wild boar that showed values compatible with those observed for omnivores (Fig 1). As anticipated, four chickens (Gallus sp.), and two dogs (Canis sp.) showed isotope values typical of omnivores, and the marine fishes (3 unidentified samples, 2 specimens of the Sparidae family, 1 Labrus sp. and 1 Sparus sp.) were more enriched in 13C compared to terrestrial species, as to be expected (δ13C = -9.3‰, S.D. 1‰; δ15N = 10.7‰, S.D. 0.7‰). Galeorhinus sp., a shark genus widely distributed in subtropical areas, represented the consumer at the upper trophic level (δ13C = -12.0‰; δ15N = 13.9‰).

Finally, the charred seeds revealed isotope values compatible with C3 photosynthesis and with a basal trophic level. The δ15N value of Vicia (mean δ15N = 2.4‰, S.D. 0.4 ‰), a legume, was low, as expected, since these plants are not dependent on soil for nitrogen uptake [6870].

For humans (n = 108), δ13C values ranged between -19.8‰ to -16.8‰ (mean -19.0, S.D. 0.5‰), while δ15N values ranged from 6.7‰ to 11.4‰ (mean 9.3, S.D. 1.1‰), as reported in Fig 1. In order to evaluate possible dietary variations between sexes and across 5 selected age groups (0–17 years seven individuals; 18–30 years 33 individuals; 31–40 years 22 individuals; 41–50 years 12 individuals; >50 years 4 individuals) isotope data were statistically compared by Wilcoxon and the Kruskal-Wallis tests. No significant differences were found between the isotopic values for males (n = 62) and females (n = 37), for δ13C (Wilcoxon test: p-value = 0.6522) or δ15N (Wilcoxon test: p-value = 0.8203) ratios. Similar results were obtained on the different age groups (for δ13C Kruskal-Wallis test: p-value = 0.4746, for δ15N Kruskal-Wallis test: p-value = 0.2056). To differentiate between individuals with a similar diet, average linkage cluster analysis was performed on the basis of isotope data. A total of five different groups of samples were identified (Fig 2).

thumbnail
Fig 2. Human isotope analysis.

Bivariate plot of average linkage cluster analysis for the human stable isotope data from Santa Severa with five groups identified.

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

The first group (cluster 1) included 71 individuals; in this context, δ13C values ranged from -18.7‰ to -19.8‰ (mean = -19.1‰, S.D. 0.2‰) and δ15N values from 8.3‰ to 10.8‰ (mean = 9.7‰ S.D. 0.7‰). Nitrogen isotope values measured for human specimens were enriched by at least 2.0‰ compared to herbivores, showing a mean human-herbivore difference (Δ15N) of 3.4‰. The second group (cluster 2) was made up of 21 individuals; for these samples, δ13C values ranged between -18.5‰ and -19.7‰ (mean = -19.1‰, S.D. 0.4‰) and δ15N values were between 6.7‰ and 8.9‰ (mean = 7.6‰, S.D. 0.6‰). The mean human-herbivore difference (Δ15N) was 1.3‰. The third set (cluster 3) comprised 3 individuals with high δ13C values; δ13C ranged between -17.9‰ to -17.7‰ (mean = -17.8‰, S.D. 0.1‰) and δ15N from 7.8‰ to 8.1‰ (mean = 8.0‰, S.D. 0.2‰). The mean human-herbivore difference (Δ15N) was 1.7‰. The forth cluster (cluster 4) counts 10 individuals, δ13C ranged from -18.6‰ to -18.0‰ (mean = -18.4‰, S.D. 0.2‰) and δ15N from 9.0‰ to 11.4‰ (mean = 10.2‰ S.D. 0.8‰) with a mean human-herbivore difference (Δ15N) of 3.9‰. The last set (cluster 5) includes 2 individuals with the highest δ13C values, where δ13C ranged between -16.8‰ to -17.5‰ (mean -17.2‰, S.D. 0.5‰) and δ15N from 9.9‰ to 10.4‰ (mean = 10.2‰ S.D. 0.4‰) with a mean human-herbivore difference (Δ15N) of 3.9‰.

To explore potential differences in diet linked with burial practice, individuals buried in earthen graves were statistically compared to those interred in sarcophagi and no differences were found for either δ13C or δ15N (Wilcoxon, δ13C: p-value = 0.5312, δ15N: p-value = 0.7830).

Dental calculus analysis

Dental calculus was collected from 94 individuals. Ancient DNA was extracted from the dental calculus of 52 individuals and specific barcode genes were amplified and sequenced, in order to investigate the consumption of specific animal food sources: sheep, chicken, cattle, pig and fish. The spike PCR control showed that no sample contained PCR inhibitors. For each individual, positive amplicons were reported in Table 3.

Overall, animal aDNA was detected in a large number of samples, sheep was present in 55.8% of the analysed individuals, cattle in 34.6% and pig in 25%. Chicken genes were, however, detected only in 11.5% of the analysed sample. Marine fish aDNA was detected in 88.5% of the specimens. These results did not show differences related to age, sex or social status of the individuals inferred from the type of burial.

LM analysis of 94 individuals indicated the presence of plant micro-remains in 70 samples (74.5%). This investigation revealed the presence of several types of micro-debris, as reported in Table 4.

In total, 22 Poaceae phytoliths, 3 calcium oxalate crystals, 2 fragments of Asteraceae inflorescence, 8 pollen grains (3 Oleaceae, 1 Chenopodiaceae, 1 Urticaceae, 3 not determined) and 390 starch granules were observed in the whole population. Only 327 starches could be taxonomically identified and assigned to eight morphotypes, on the basis of morphological and morphometric criteria using both our experimental collection [62] and the International Code for Starch Nomenclature [71]. Each morphotype is described below.

Morphotype I (tot. 24 micro-remains) showed the typical features of Fabaceae granules. They were reniform, oval to elongated in shape and ranged in size between 7–35 μm in length and 3–30 μm in width. Hilum was not detectable (obscured). Clear concentric lamellae were present and, in some granules, a longitudinal crack in the amorphous central area could be observed. Some grains appeared ascribable to starches of Pisum L. and Vicia L. genera.

Morphotype II (tot. 29) consisted of starch granules elongated, drop or oblong-shaped (dimensional range: 5–28 μm in length, 4–20 μm in width). They presented an invisible and eccentric hilum, faintly visible lamellae and longitudinal fractures, which usually occur in several species of Fagaceae, such as Quercus ilex L.

Morphotype III (tot. 102) was characterised by a bimodal distribution of disc-shaped grains, typical of grasses of the Triticeae tribe. Among large diagnostic granules (size range: 6–35 μm in length and 4–33 μm in width), 18 were reminiscent of Triticum sp. L. (e.g. T. durum Desf., T. aestivum L., T. dicoccon Schrank ex Schübl.), 13 of Hordeum sp. L. (e.g. H. vulgare L.) and 3 of Secale sp. L.

Morphotype IV (tot. 114) consisted of polyhedral-shape units from one side and dome shaped on the other. Each unit showed a size ranging from 3–11 μm both in length and in width, a centric-indistinct hilum and indistinct lamellae. Sometimes, they were found in aggregates. These features were attributable to Avena sp. L.

Morphotype V (tot. 3) revealed a morphology which completely fitted with those of Daucus carota L. grains. In detail, spherical shape, size range from 5–10 μm in length to 4–8 μm in width, peculiar multiple radial fractures, slightly visible concentric lamellae and centric hilum were detected in these granules.

Morphotype VI (tot. 40) was characterized by polyhedral-shape units with blunted edges (size range: 2.5–5 μm both in length and in width). This distinctive trait, associated to a centric-distinct hilum and indistinct lamellae, suggested that these starches belonged to Panicum sp. L.

Morphotype VII (tot. 2) showed polyhedral units (size range from 1 to 2.6 μm both in length and in width) with concave faces, acute edges and typical bright boundaries. Hilum was centric, while lamellae appeared indistinct. All these features indicated Piper sp. L. starches.

Morphotype VIII (tot. 13) consisted of ovoidal granules with some peculiar flattened surfaces. Their size ranged from 5–18 μm in length to 4–14 μm in width and radial fissures starting from the centric hilum. These characteristics were consistent with Sorghum sp. Moench starches.

Morphotypes III and IV (Poaceae starches) were identified in most of the individuals, followed by morphotypes I and II (Fabaceae and Fagaceae starches). Particular attention was paid to individual NS SU 321 whose tartar contained 21 starch granules mainly ascribable to morphotype I and II. No difference was evidenced among individuals of different age at death and sex. Representative images of microfossils found in tartar samples are reported in Fig 3.

thumbnail
Fig 3. Plant microremains and starch granules at LM.

Representative images of microdebris found in dental calculus samples: Oleaceae pollen grain (A); Chenopodiaceae pollen grain (B); Poaceae phytolith (C); fragment of Asteraceae inflorescence perforation plate (D); aggregate of Avena sp. starch granules and relative polarized image (E); Fabaceae starch granule and relative polarized image (F); Sorghum sp. starch granules and relative polarized image (G); Triticeae starch granule and relative polarized image (H); Triticeae starch granules (I and J); Fagaceae starch granule (K); Fabaceae starch granule (L). The scale bar indicates 30 μm.

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

GC-MS analysis was performed only on the samples who presented sufficient dental calculus (68 individuals). In S5 Table, the molecules identified by this approach are listed and clustered in biochemical classes per individual. Food categories or specific plants were deduced by associating the recognized compounds. All chromatographic profiles revealed a significant presence of n-alkanes and n-alkenes (C6-C35). Generally, monounsaturated (e.g. docosenoic and octadecenoic acids) and polyunsaturated (e.g. octadecadienoic and octadecatrienoic acids) fatty acids, ω3-fats (e.g. eicosapentaenoic and docosahexaenoic acids) and lactose were the most recurrent molecules detected in the population, followed by cholesterol and secondary metabolites indicating aromatic herbs (e.g. anisole, alpha-cubebene, estragole, santolina triene, beta-copaene, dehydroelsholtzia ketone). The identification of phytosterols (e.g. stigmasterol, sitosterol, campesterol), vitamins (e.g ascorbic acid and tocopherol) and other classes of plant markers (e.g. lactones, glucosinolates, phenolics) suggested the key role of vegetables and fruits in the diet of the population. Moreover, the identification of tartaric and gallic acids, pyrocathecol, and pyrogallol hypothesized a possible consumption of wine. Evidence of alkaloids (e.g., pseudopelletierine, ephedrine) and terpenic compounds (e.g. bisabolene, scoparone, alpha-bisabolol, dihydroartemisinin, all ascribable to Asteraceae) supported the use of aromatic and medicinal plants. Finally, detection of ergosine and bovinocidin in two samples suggested a fungal contamination of stored foods. No difference was evidenced among individuals showing different age at death, social class and sex.

Discussion

Here we have applied a suite of scientific techniques to the skeletal remains of the Medieval community of Santa Severa (7th-15th century CE) in order to highlight their culture and food habits. This research integrates the available morphological data [24] with the results obtained by the present original and innovative multidisciplinary approach, with the aim to reconstruct the osteobiography of one of the largest Italian Medieval burial populations.

Faunal specimens showed high variability in their isotopic values. The differences observed between the deer (Cervus sp.) specimen and the other herbivores (6 cattle Bos sp., 5 sheep Ovis sp., 1 horse Equus sp., and 1 buffalo Bubalus sp.) could be attributable to the different environments in which these animals lived [31]. Deer tend to live in wooded environments, while domestic animals would have lived in closer proximity to humans and would have been managed by them. Therefore, the observed 15N enrichment in sheep and cattle could be ascribable to both human activities and potentially proximity to the sea [31]. Some crop management practices (e.g. manuring) could increase δ15N values in soil and plants, while plants might also have absorbed nitrogen of marine origin, due to the “sea spray effect” [7273]. Sus sp. samples had dietary signals more similar to those observed for herbivores than omnivores, as in earlier studies [74] although the wild boar (Sus scrofa) had significantly different values and presented an omnivorous diet. In Medieval North-Western Europe, the diet of pigs was mainly based on terrestrial plants and human refuse [17, 18, 75], although in some areas, including Latium, pigs were free to roam in the uncultivated land surrounding the city [17] where they could consume acorns [15, 17].

The stable isotope data for humans indicated an omnivorous diet with a contribution of both animal products and plants. The high variability observed in nitrogen isotopic values (ranging from 6.7‰ to 11.4‰) seemed to witness a probable differential access to food sources within the population.

The δ13C values for humans were compatible with a diet mainly based on consumption of C3 plants; however, the highest δ13C values (δ13C > 18‰) observed in some specimens (cluster 3) could also suggest a contribution of C4 plants to the diet. This is not surprising as during the Middle Ages, C4 crops such as sorghum and millets, were grown and consumed due to their ease of cultivation and relatively high yields [15, 7677]. These results are in line with the archaeobotanical results from dental calculus described below.

Given the coastal location and late Medieval date, it is surprising that the isotope data do not reflect a significant contribution of marine protein to the human diet. The identification of a marine isotopic signal with typical enrichment in both 13C and 15N however, is problematic in the Mediterranean context where the contribution of C4 terrestrial protein exists, as Mediterranean fish tend to have lower δ13C and δ15N values than those of the Atlantic [7880]. Fish consumption can not be totally excluded as marine fish aDNA and ω3-fatty acids were detected in dental calculus, although the latter are also widely abundant in plants. Aside from the issue of C4/marine diets, exhibiting similar values [7880], it may also be possible that the consumption of marine protein sources was not sufficient to induce a significant shift of isotopic values [81].

Only three individuals (NS SU 124 Aa, PR SU 262 Aa and NS SU 321) possessed isotopic values potentially compatible with fish consumption (δ13C -16.8‰, -17.5‰, -18.1‰ and δ15N 10.4‰, 9.9‰, 11.3‰ respectively). However, due to the presence of non-specific stress markers (e.g. cribra), as well as infections (e.g periostitis and osteomyelitis) in two of these specimens (NS SU 124 Aa, NS SU 321) [24], the hypothesis of higher nitrogen values for these individuals being affected by potential nutritional stress cannot be totally excluded [82]. It has been documented that starvation may stimulate gluconeogenesis, and therefore the production of glucose from non-carbohydrate sources [47, 8387], causing 15N enrichment in body tissues [83, 88]. The isotope values for NS SU 321 in particular, however, may be reflective of a high status diet. This individual, atypically buried in a sarcophagus with a cross and a stone cushion [23], presented the markers for all the tested food sources (bovine, swine, ovine, chicken and fish: Table 3), according to genetic analysis of dental calculus. The consumption of high trophic level protein (meat and to some extent fish) was affiliated with a high-status identity in Medieval society [15]. Although no statistical difference was detected between sex and age at death groupings, a probable internal division of the society was suggested by the average linkage cluster analysis (Fig 2). The clusters differ on both carbon and nitrogen isotopic signatures. In particular, as regard the former, three of them (cluster 1, cluster 2, and cluster 4) include individuals possessing isotopic values compatible with a prevalent consumption of C3 plants whereas carbon isotopic values of the other two clusters (cluster 3 and cluster 5) are compatible with the intake of C4 species. Similarities in δ13C values, however do not square with nitrogen isotopic ratios. The lowest nitrogen detected in clusters 2 and 3, both also therefore characterized by the lowest human-herbivore 15N enrichment (1.3‰ and 1.7‰, respectively) suggesting a more vegetarian diet with a reduced animal protein intake.

Alternatively, a higher contribution of animal (mainly terrestrial) protein sources in diet is indicated by the higher δ15N values observed for the individuals grouped into the other three clusters that also posess higher human-herbivore differences (3.4‰, 3.9‰, and 3.9‰ respectively). None of the clusters are defined by individuals showing clear isotopic signatures of fish consumption. The above mentioned three individuals for which fish consumption is most likely are grouped in clusters 4 (NS SU 321) and 5 (NS SU 124 Aa, PR SU 262 Aa). These isotopic clusters, however, are not correlated with the different typologies of burials identified in the archaeological site (e.g. simple earthen graves and sarcophagi made of Etruscan and Roman re-used tuff stones), as they are equally represented within the clusters. The variability of nitrogen isotopic values detected in the population may instead be the result of the wide time span that these burials represent (about 900 years, 7th-15th century CE). During the Early Middle Ages (6th-10th century), the diet, in Italy, was typically associated with the consumption of a large amount of animal protein (meat and fish), in association with plant foods [77, 89]. During the Late Medieval period (11th-15th century), the intake of appreciable amounts of animal protein became a privilege of the upper classes [15, 77]. A potential dietary transition between the Early and Late Middle Ages may also be the result of ecclesiastic law requiring abstinence from meat promoting, instead the consumption of fish during the later period [2, 718]. Unfortunately, we are unable to further explore the chronological element, as it has not been possible to identify which burials derived from the early and later Medieval period.

The isotope data for published Medieval populations from Italy is generally very variable (S6 Table). The highest δ13C values can be observed by North-Eastern sites, dating to the early Medieval period (6th-11th centuries), like Romans d’Isonzo, Cividale Gallo, Cividale Santo Stefano and Mainizza [90], as well as Cosa (11th-13th centuries; [86]) from Central Italy. The observed enrichment in 13C of these samples, with respect to Santa Severa, is attributable to a higher consumption of C4 plants (e.g. millet) [86, 90]. Indeed, North-Eastern Italy, experienced an onward decline in bread consumption which was substituted by soups [90]. With regard to nitrogen isotopes, the population from Santa Severa was more enriched in 15N than several sites across Medieval Italy, namely Romans d’Isonzo and Mainizza in the North East [90], Colonna (8th-10th centuries; [19]) in Central Italy, and Montella (13th-15th centuries; [91]) in Southern Italy. The North-Eastern Medieval populations were characterized by a low animal protein intake, (although freshwater fish intake has been suggested for some [90]), whereas at Montella (Southern Italy), the analysed population was associated with a Franciscan friary and so there is the potential that this monastic population followed a distinctive dietary regime [91]. Although Colonna and Santa Severa are present in the same region of Italy, a differential access to nutritional sources related to sex and age at death in this case, was detected at Colonna [19]. In particular, at Colonna, adult male individuals demonstrated a greater intake of animal protein with respect to females and juveniles, who followed a more vegetarian diet [19].

When considering the results from dental calculus, it should be borne in mind that the number of starch granules and molecules present in dental calculus of an individual is not necessarily proportional to the amount of foods consumed in life. However, these data, considered at population level, represent direct evidence of the substances that entered the oral cavity [9293]. LM analysis indicated the presence of C3 cereals (e.g. Avena sp. and Triticeae), Fabaceae (e.g. Vicia sp., Pisum sp.) and Fagaceae (e.g. Quercus sp. L.) seeds, in more than half of the individuals. Moreover, for some specimens, we hypothesized the intake of C4 caryopses (e.g. Sorghum sp. and Panicum sp.). In general, the finding of these types of starches is in keeping with the carpological remains discovered in an oven of the same archaeological site [23], supporting the use of these plant species as the main sources of carbohydrates. The presence of Fagaceae starches is not surprising because acorns have been recognised as a nutritional resource containing sugars, vitamins and proteins since the Roman period. After the removal of indigestive tannins, acorns were powdered and used to ennoble cereal flours or prepare astringent and anti-diarrheal decoctions [22, 9396]. Similarly, legumes represented both an important food that could be preserved for a long time after drying, and a source of nutraceutical compounds applied in dermatology [9798]. Moreover, these species were also excellent feed crops, already used in the Roman three-field rotation strategy (together with cereals and fallow), for their ability to increase soil nitrogen [99100]. The high number of starch grains that could not be identified may be the result of cooking processes, grinding procedures or exposure to Ptialin enzyme activity [101]. The finding of Oleaceae, Chenopodiaceae, Urticaceae pollen grains and fragments of Asteraceae inflorescence testified the existence of these plant families in the studied area. The possible use of these species in dietary and/or ethnopharmacological traditions cannot be excluded. For instance, Urtica dioica L. (nettle) and Parietaria officinalis L. (upright pellitory) were widely employed for medicinal purposes, as diuretic, emollient and expectorant [102].

The GC-MS analysis on dental calculus revealed a great variety of n-alkanes and n-alkenes, probably deriving from degradation of both plant and animal food molecules or/and oral microbiota [103104]. The detection of mono/poly-unsaturated fatty acids, phytosterols and plant vitamins suggested a diet based on plant foods, such as seeds (e.g. cereals, confirmed by LM), vegetables (e.g. Brassicaceae, also supported by the presence of glucosinolates) and fruits (e.g. of Rosaceae family due to the presence of lactones) [3, 105107]. Moreover, ω3-fatty acids could be associated with ingestion of plant and/or aquatic resources (e.g. dried fruits, seaweeds, molluscs, blue fishes) [108109]. A great part of the community also had lactose present in their dental calculus, the main sugar of milk and dairy products, confirming the fundamental role of these animal derivatives in the Medieval diet [110111]. Some individuals presented secondary metabolites of aromatic plants, such as Apiaceae (e.g. dill) and Lamiaceae (e.g. sage), leading us to hypothesize that they employed these herbs as food preservatives and/or taste regulators [112]. Our data confirmed the consumption of wine (Vitis vinifera L.) [113115]. The pseudopelletierine, an alkaloid found in NS SU 124 Aa sample, proved the ingestion of Punica granatum L. bark [116]. Indeed, various portions of pomegranate were used in folk medicine, including Egyptian one; flowers for treating diabetes; fruits for expelling parasites; seeds and fruit peels for managing diarrhoea; bark and roots as coagulants and anti-ulcer remedies [116118]. Pedanius Dioscorides, a Greek physician (40–90 AD), and several Roman documents report the recipes for the preparation of mouthwashes based on pomegranate rind and bark [119120]. Indeed, it has actually been recently scientifically demonstrated that the extracts of this species prevent the dental plaque development [118, 121]. A typical marker of Elsholtzia sp. Willd. (dehydroelsholtzia ketone) and a peculiar chemical compound of Ephedra sp. L. (ephedrine) indicated a possible intake of decoctions based on these plants, probably to treat respiratory and gastro-intestinal disorders [122123].

In three individuals (NS SU 27 Ac, NS SU 115, NS SU 290) we found the artemisinin, a sesquiterpene lactone typically synthesized in Artemisia annua L. (Asteraceae family) and presenting antimalarial properties [124]. As Santa Severa Castle was surrounded by marshy areas during the Medieval period [125], a condition which favour the survival of the Plasmodium falciparum, our results might sustain the existence of malaria in Central Italy and its therapeutic treatment by the use of this plant. Unfortunately, no distinctive signs on bone remains can document and support the previous hypothesis. As extensively investigated by earlier studies [126127] malaria determines a haemolytic anaemia that may result in several modifications on human bones even though it can be considered one of the possible causes of cribra orbitalia [126129]. Indeed, cribra orbitalia, represents one of the most frequently recorded alterations in archeaeological skeletal collections; it pertains to the presence of small foramina in the orbit vaults, due to the expansion of the diploe accompanied by the narrowing of the outer bone cortex as a consequence of bone marrow hypertrophy [127, 130132]. Cribra orbitalia has been traditionally considered a clue of iron deficiency anaemia [127, 130132], it is worth noting, however, that defining the exact cause of anaemic conditions could be challenging as they could have been determined by several factors, such as iron or ascorbic acid depletion, other nutritional deficiencies, infections, or even by a combination of these elements [47, 127, 133136]. Moreover, a similar bone morphology could be also determined by scurvy and chronic infections [47, 127, 130132]. As reported by Gowland and Western [127] malaria can be directly diagnosed only by the detection of Plasmodium species DNA in human skeletal remains [127], as successfully carried out in Egyptian mummies [127, 137138] and in individuals recovered in an infant cemetery in Umbria (Italy, 5th century CE; [139]) [127]. The use of Asteraceae, as food or medicinal sources [140], in these individuals was also confirmed by the detection of bisabolol (in NS SU 27 Ac), bisabolene (in NS SU 27 Ac and NS SU 115) and the LM survey of fragments of capitulum inflorescence (in NS SU 115 and NS SU 290).

Lastly, the presence of two molecules, the bovinocidin (produced by food-contaminating molds such as Aspergillus sp. P. Micheli ex Haller) [141] and the ergosine (an alkaloid synthesized by Claviceps purpurea (Fr.) Tul., the fungal parasite of Poaceae) [142], could indicate the use of inadequate practices in food conservation and the diffusion of Graminaceae phytopathogens in that historical context.

Conclusions

In conclusion, the cutting-edge multidisciplinary approach applied in this work has enabled a detailed reconstruction of the dietary habits of one of the largest Italian Medieval populations analysed to date. The isotope analyses performed on human, faunal and seed samples, combined with archeobotanical and molecular analysis of dental calculus, revealed that the population generally adopted an omnivorous diet, in which plants and fruits played a prominent role. With regard to animal protein intake, cattle, sheep and pig meat were consumed with a minor contribution of chicken. Consumption of C4 plants (millet and/or sorghum) was also evidenced in a few individuals. Significantly, the application of multiple complementary methodologies revealed a potential intake of marine foodstuffs as dietary resource for this Medieval population. The consumption of marine foods was not clearly detectable using stable isotope evidence, highlighting the limitations of using this method alone in the Mediterranean context. Archaeobotanical analysis further revealed the presence of artemisinin that may be an indication of the presence of malaria in the area of Santa Severa. The dietary (e.g. cereals, Brassicaceae) and medicinal (e.g. pomegranate, Ephedra sp.) uses of other plant species was also hypothesised. This innovative archaeobotanical approach combined with the anthropological data may provide interesting information about presence and diffusion of pathologies and diseases in the past, beyond those are not directly detectable on skeletal remains.

Supporting information

S1 Table. Biological profile of the Medieval individuals subjected to dental calculus analyses.

Each sample was associated to a code, reporting NS (for identifying the archaeological area of “Casa del Nostromo”) and a number (relative to the stratigraphic unit, SU), according to [24]. For each specimen, estimated sex (F, female; M, male; ND, sex determination was not determinable for the lack of diagnostic elements due to their bad conservation; IND, sub-adults with not determined sex because of sexual immaturity), estimated age at death (GA, generic adult individual) and district of the masticatory apparatus used for dental calculus sampling were reported. In particular, sampled teeth and relative surfaces were codified according to the Universal Teeth Numbering System [143] (B, buccal; L, lingual; M, mesial; D, distal). Finally, sample weight was indicated in grams and the analytical techniques applied on each sample were marked with an X.

https://doi.org/10.1371/journal.pone.0227433.s001

(DOCX)

S2 Table. Light microscopy results of the laboratory contamination tests.

https://doi.org/10.1371/journal.pone.0227433.s002

(DOCX)

S3 Table. Light microscopy results obtained from the washing water applied on ancient dental calculus before the cleaning procedure.

https://doi.org/10.1371/journal.pone.0227433.s003

(DOCX)

S4 Table. The primers used for Spike PCR (L15996 and H16401) and detection of aDNA relative to bovine (Bov84/90-F/R), pig (Sus85-F/R; Sus98-F/R), ovine (Ovis-F/R), chicken (Gall-F/R) and fishes (Fish_miniA_F/R; Fish_miniC_F/R) were shown.

In addition, length of PCR amplicons and annealing temperatures were also reported.

https://doi.org/10.1371/journal.pone.0227433.s004

(DOCX)

S5 Table. The chemical compounds detected in dental calculus by GC-MS analysis were listed and clustered in biochemical classes for each sample.

https://doi.org/10.1371/journal.pone.0227433.s005

(DOCX)

S6 Table. Mean of δ13C and δ15N of human samples recovered in Santa Severa archaeological site and in other coeval Italian sites from literature.

For each sample. chronology, sample size and relative reference were reported [7, 8, 19, 20, 21, 86, 90, 91, 144145].

https://doi.org/10.1371/journal.pone.0227433.s006

(DOCX)

Acknowledgments

The authors thank Harry Robson (University of York) for his help in identifying the fish remains and we are grateful to Luke Spindler (University of York) for his help in collagen extraction for fish.

The authors sincerely thank the editor and the anonymous reviewer for the valuable comments that contributed to improve the present paper.

References

  1. 1. Pearson MP. Food, Culture and Identity in the Neolithic and Early Bronze Age. Oxford: BAR International Series 1117; 2003.
  2. 2. Reitsema LJ, Crews DE, Polcyn M. Preliminary evidence for medieval Polish diet from carbon and nitrogen stable isotopes. J. Archaeol. Sci. 2010; 37: 1413–1423.
  3. 3. Montanari M. Gusti del Medioevo: i prodotti, la cucina, la tavola. Roma-Bari: Laterza; 2014.
  4. 4. White L. Medieval technology and social change. Oxford: Oxford University Press; 1962.
  5. 5. Adamson MW. Food in medieval times. Westport, CT: Greenwood; 2004.
  6. 6. Vercellotti G, Stout SD, Boano R, Sciulli P. Intrapopulation variation in stature and body proportions: social status and sex differences in an Italian medieval population (Trino Vercellese, VC). Am. J. Phys. Anthropol. 2011; 145: 203–214. pmid:21312185
  7. 7. Reitsema LJ, Vercellotti G. Stable isotope evidence for sex- and status-based variations in diet and life history at medieval TrinoVercellese, Italy. Am J Phys Anthropol. 2012; 148: 589–600. pmid:22553011
  8. 8. Ciaffi R, Lelli R, Müldner G, Stantcheva K, Fischetti AL, Ghini G, et al. Palaeobiology of the medieval population of Albano (Rome, Italy): a combined morphological and biomolecular approach. Int. J. Osteoarchaeol. 2015; 25: 477–488.
  9. 9. Woolgar C. ‘‘Take this penance now, and afterwards the fare will improve”: seafood and late medieval diet. In: Starkey DJ, Reid C, Ashcroft N, editors. England’s sea fisheries: the commercial sea fisheries of England and Wales since 1300. London, Chatham; 2000. pp. 36–44.
  10. 10. Barrett JH, Richards MP. Identity, gender, religion and economy: new isotope and radiocarbon evidence for marine resource intensification in early historic Orkney, Scotland, UK. Eur. J. Archaeol. 2004; 7: 249–271.
  11. 11. Nada Patrone AMN. Il cibo del ricco ed il cibo del povero: contributo alla storia qualitativa dell’alimentazione; l’area pedemontana negli ultimi secoli del Medio Evo. Torino: Centro Studi Piemontesi; 1981.
  12. 12. Lanconelli A. Gli Statuta pescivendulorum Urbis (1405). Note sul commercio di pesce a Roma fra XIV e XV secolo. Archivio della Società Romana Storia Patria 1985; 108: 84–131.
  13. 13. Lanconelli A. I lavori alla peschiera del Marta. Contributo alla storia della pesca nel Lazio bassomedievale, Scritti in memoria di Giuseppe Marchetti Longhi. Anagni: Istituto di Storia e di Arte del Lazio meridionale; 1990.
  14. 14. Zug Tucci H. Il mondo medievale dei pesci tra realtà e immaginazione. L’uomo di fronte al mondo animale nell’Alto Medioevo. Spoleto: Centro Italiano di Studi sull’Alto Medioevo; 1985.
  15. 15. Montanari M. Alimentazione e cultura nel Medioevo. Roma-Bari: Laterza; 1988.
  16. 16. Salamon M, Coppa A, McCormick M, Rubini M, Vargui R, Tuross N. The consilience of historical and isotopic approaches in reconstructing the medieval Mediterranean diet. J Archaeol Sci. 2008; 35: 1667–1672.
  17. 17. Müldner G, Richards MP. Fast or feast: reconstructing diet in later medieval England by stable isotope analysis. J Archaeol Sci. 2005; 32: 39–48.
  18. 18. Müldner G, Richards MP. Stable isotopic evidence for 1500 years of human diet at the city of York, UK. Am J Phys Anthropol. 2007; 133: 682–697. pmid:17295296
  19. 19. Baldoni M, Nardi A, Müldner G, Lelli R, Gnes M, Ferraresi F, et al. Archaeo-biological reconstruction of the Italian medieval population of Colonna (8th-10th centuries CE). J Archaeol Sci Rep.2016; 10: 483–494.
  20. 20. Baldoni M, Scorrano G, Gismondi A, D’Agostino A, Alexander M, Gaspari L, et al. Who were the miners of Allumiere? A multidisciplinary approach to reconstruct the osteobiography of an Italian worker community. PLoS ONE 2018; 13: e0205362. pmid:30308078
  21. 21. Baldoni M, Scorrano G, Alexander M, Stasolla FR, Marsella LT, Rickards O, et al. The medieval population of Leopoli-Cencelle (Viterbo, Latium): Dietary reconstruction through stable isotope analysis from bone proteins. J Archaeol Sci Rep. 2019; 24: 92–101
  22. 22. Gismondi A, D'Agostino A, Canuti L, Di Marco G, Martínez-Labarga C, Angle M, et al. Dental calculus reveals diet habits and medicinal plant use in the Early Medieval Italian population of Colonna. J Archaeol Sci Rep. 2018a; 20: 556–564.
  23. 23. Enei F. Santa Severa—Tra leggenda e realtà storica. Viterbo: Ceccarelli Editrice S.r.l.; 2013.
  24. 24. Gnes M, Baldoni M, Marchetti L, Basoli F, Leonardi D, Canini A, et al. Bioarchaeological approach to the study of the medieval population of Santa Severa (Rome, 7th-15th centuries). J Archaeol Sci Rep. 2018; 18: 11–25.
  25. 25. Henry AG, Piperno DR. Using plant microfossils from dental calculus to recover human diet: a case study from Tell al-Raqā'i, Syria. J Archaeol Sci 2008; 35: 1943–1950.
  26. 26. Craig OE, Biazzo M, Colonese AC, Di Giuseppe Z, Martinez-Labarga C, Lo Vetro D, et al. Stable isotope analysis of Late Upper Palaeolithic human and faunal remains from Grotta del Romito (Cosenza), Italy. J Archaeol Sci. 2010; 37: 2504–2512.
  27. 27. Cristiani E, Radini A, Borić D, Robson HK, Caricola I, Carra M, et al. Dental calculus and isotopes provide direct evidence of fish and plant consumption in Mesolithic Mediterranean. Sci Rep. 2018; 8, 8147. pmid:29802341
  28. 28. Katzenberg MA, Waters‐Rist AL. Stable isotope analysis: a tool for studying past diet, demography, and life history. In: Katzenberg MA, Saunders SR, editors. Biological anthropology of the human skeleton. Hoboken: Wiley-Blackwell; 2018. pp.467–504.
  29. 29. Katzenberg MA. Stable isotope analyses: a tool for studying past diet, demography and life history. In: Katzenberg MA, Saunders SR, editors. The biological anthropology of the human skeleton. New York: Wiley; 2000. pp. 305–327.
  30. 30. Schoeninger MJ. Diet reconstruction and ecology using stable isotope ratios. In: Larsen CS, editor. A comparison to biological anthropology. Chichester, UK: Wiley-Blackwell; 2011. pp. 445–464.
  31. 31. Hedges REM, Reynard L. Nitrogen isotopes and the trophic level of humans in archaeology. J. Archaeol. Sci. 2007; 34: 1240–1251.
  32. 32. Bocherens H, Drucker D. Trophic level isotopic enrichment of carbon and nitrogen in bone collagen: case studies from recent and ancient terrestrial ecosystems. Int. J. Osteoarchaeol. 2003; 13: 46–53.
  33. 33. Tykot RH. Stable isotopes and diet: you are what you eat. In: Martini M, Milazzo M, Piacentini M, editors. Proceedings of the International School of Physics “Enrico Fermi”, Course CLIV. Amsterdam: IOS Press; 2004. pp. 433–444.
  34. 34. Lieverse AR. Diet and the aetiology of dental calculus. Int J Osteoarchaeol. 1999; 9: 219–232.
  35. 35. Warinner C, Speller C, Collins MJ. A new era in palaeomicrobiology: prospects for ancient dental calculus as a long-term record of the human oral microbiome. Phil. Trans. R. Soc. B 2015; 370: 20130376. pmid:25487328
  36. 36. Radini A, Tromp M, Beach A, Tong E, Speller C, McCormick M, et al. Medieval women’s early involvement in manuscript production suggested by lapis lazuli identification in dental calculus. Sci Adv. 2019; 5: eaau7126. pmid:30662947
  37. 37. Hendy J, Warriner C, Bouwman A, Collins MJ, Fiddyment S, Fisher R, et al. Proteomic evidence of dietary sources in ancient dental calculus. Proc. R. Soc. B 2018; 285: 20180977 pmid:30051838
  38. 38. Cristiani E, Radini A, Edinborough M, Borić D. Dental calculus reveals Mesolithic foragers in the Balkans consumed domesticated plant foods. PNAS 2016; 113: 10298–10303. pmid:27573829
  39. 39. Hardy K, Radini A, Buckley S, Sarig R, Copeland L, Gopher A, et al. Dental calculus reveals potential respiratory irritants and ingestion of essential plant-based nutrients at Lower Palaeolithic Qesem Cave Israel. Quat Int. 2016; 398: 129–135.
  40. 40. Weyrich LS, Duchene S, Soubrier J, Arriola L, Llamas B, Breen J, et al. Neanderthal behaviour, diet, and disease inferred from ancient DNA in dental calculus. Nature 2017; 544: 357–361. pmid:28273061
  41. 41. Cummings LS, Yost C, Sołtysiak A. Plant microfossils in human dental calculus from Nemrik 9, a Pre-Pottery Neolithic site in Northern Iraq. Archaeol Anthropol Sci. 2018; 10: 883–891.
  42. 42. Eerkens JW, Tushingham S, Brownstein KJ, Garibay R, Perez K, Murga E, et al. Dental calculus as a source of ancient alkaloids: Detection of nicotine by LC-MS in calculus samples from the Americas. J Archaeol Sci Rep. 2018; 18: 509–515.
  43. 43. Power RC, Salazar-García DC, Rubini M, Darlas A, Havarti K, Walker M, et al. Dental calculus indicates widespread plant use within the stable Neanderthal dietary niche. J Hum Evol. 2018; 119: 27–41. pmid:29685752
  44. 44. Fatucci M, Cerilli E. Primi risultati dallo studio dei resti faunistici rinvenuti nello scavo effettuato al castello di Santa Severa “Piazza della Rocca”. In: Enei F editor. Santa Severa—Tra leggenda e realtà storica Viterbo: Ceccarelli Editrice S.r.l.; 2013. pp. 224–233.
  45. 45. Susanna F. Studio archeobotanico di macroresti provenienti dal sito di Santa Severa. Primi risultati. In: Enei F, editor. Santa Severa—Tra leggenda e realtà storica. Viterbo, Italy: Ceccarelli Editrice S.r.l; 2013. pp. 234–236.
  46. 46. Longin R. New method of collagen extraction for radiocarbon dating. Nature 1971; 230: 241. pmid:4926713
  47. 47. Scorrano G, Baldoni M, Brilli M, Rolfo MF, Fornaciari G, Rickards O, et al. Effect of the Neolithic transition on an Italian community: Mora Cavorso (Jenne, Rome). Archaeological and Anthropological Sciences 2019; 11:1443–1459.
  48. 48. Brinkkemper O, Braadbaart F, van Os B, van Hoesel A, van Brussel AAN, Fernandes R. Effectiveness of different pre‐treatments in recovering pre‐burial isotopic ratios of charred plants. Rapid Commun Mass Spectrom. 2018; 32: 251–261. pmid:29134706
  49. 49. DeNiro MJ. Postmortem preservation and alteration of in vivo bone collagen isotope ratios in relation to palaeodietary reconstruction. Nature 1985; 317: 806–809.
  50. 50. Ambrose SH. Preparation and characterization of bone and tooth collagen for stable carbon and nitrogen isotope analysis. J Archaeol Sci. 1990; 17: 431–451.
  51. 51. van Klinken GJ. Bone collagen quality indicators for palaeodietary and radiocarbon measurements. J Archaeol Sci. 1999; 26: 687–695.
  52. 52. Crowther A, Haslam M, Oakden N, Walde D, Mercader J. Documenting contamination in ancient starch laboratories. J Archaeol Sci. 2014; 49: 90–104.
  53. 53. Gismondi A, Di Marco G, Martini F, Sarti L, Crespan M, Martínez-Labarga C, et al. Grapevine carpological remains revealed the existence of a Neolithic domesticated Vitis vinifera L. specimen containing ancient DNA partially preserved in modern ecotypes. J Archaeol Sci. 2016; 69: 75–84.
  54. 54. Soto M, Inwood J, Clarke S, Crowther A, Covelli D, Favreau J, et al. Structural characterization and decontamination of dental calculus for ancient starch research. Archaeol Anthropol Sci. 2019;
  55. 55. Scorrano G, Valentini F, Martínez-Labarga C, Rolfo MF, Fiammenghi A, Lo Vetro D, et al. Methodological strategies to assess the degree of bone preservation for ancient DNA studies. Ann Hum Biol. 2015; 42: 10–19. pmid:25231926
  56. 56. Warinner C, Rodrigues JF, Vyas R, Trachsel C, Shved N, Grossmann J, et al. Pathogens and host immunity in the ancient human oral cavity. Nature Genet 2014; 46: 336–434. pmid:24562188
  57. 57. Gilbert MTP, Bandelt HJ, Hofreiter M, Barnes I. Assessing ancient DNA studies. Trends Ecol Evol. 2005; 20: 541–544. pmid:16701432
  58. 58. Willerslev E, Cooper A. Ancient DNA. Proc R Soc Lond [Biol] 2005; 272: 3–16. pmid:15875564
  59. 59. De Angelis F, Scorrano G, Rickards O. DNA antico: principi e metodologie. Medicina dei Secoli 2013; 25: 51–83.
  60. 60. Shokralla S, Hellberg RS, Handy SM, King I, Hajibabaei M. A DNA Mini-Barcoding System for Authentication of Processed Fish Products. Sci Rep. 2015; 5: 15894. pmid:26516098
  61. 61. Natonek-Wiśniewska M, Krzyścin P, Piestrzynska-Kajtoch A. The species identification of bovine, porcine, ovine and chicken components in animal meals, feeds and their ingredients, based on COX I analysis and ribosomal DNA sequences. Food Control 2013; 34: 69–78.
  62. 62. Gismondi A, D’Agostino A, Canuti L, Di Marco G, Basoli F, Canini A. Starch granules: a data collection of 40 food species. Plant Biosyst. 2019; 153: 273–279.
  63. 63. Bucher E. Lo spettro pollinico dei mieli dell'Alto Adige. Agenzia provinciale per la protezione dell'ambiente e la tutela del lavoro, Laboratorio biologico. Alto Adige: Karodruck; 2004.
  64. 64. Martin AC, Harvey WJ. The Global Pollen Project: a new tool for pollen identification and the dissemination of physical reference collections. Methods Ecol Evol. 2017; 8: 892–897.
  65. 65. Gismondi A, Canuti L, Rocco G, Pisani M, Ghelli A, Bonanno M, et al. GC–MS detection of plant pigments and metabolites in Roman Julio-Claudian wall paintings. Phytochem Lett. 2018b; 25: 47–51.
  66. 66. FoodDB, 2013. version 1.0. http://fooddb.ca/.
  67. 67. TGSC, 2015. The Good Scents Company. http://www.thegoodscentscompany.com/.
  68. 68. Virginia RA, Delwiche CC. Natural 15N abundance of presumed N2-fixing plants from selected ecosystems. Oecologia 1982; 54: 317–325. pmid:28309955
  69. 69. DeNiro MJ. Stable isotope and archaeology. Am Sci. 1987; 75: 182–191.
  70. 70. Hoefs J. Stable isotope geochemistry. New York: Springer-Verlag; 1987.
  71. 71. ICSN, 2011. The international code for starch nomenclature (http://fossilfarm.org/ICSN/Code.html) (accessed on 22 April 2018).
  72. 72. Ambrose SH. Effects of diet, climate and physiology on nitrogen isotope abundances in terrestrial foodwebs. J Archaeol Sci. 1991; 18: 293–317.
  73. 73. Britton K, Müldner G, Bell M. Stable isotope evidence for salt-marsh grazing in the Bronze Age Severn Estuary, UK: implications for palaeodietary analysis at coastal sites. J Archaeol Sci. 2008; 35: 2111–2118.
  74. 74. Craig OE, Biazzo M, O’Connell TC, Garnsey P, Martínez-Labarga C, Lelli R, et al. Stable Isotopic Evidence for Diet at the Imperial Roman Coastal Site of Velia (1st and 2nd Centuries AD) in Southern Italy. Am. J. Phys. Anthropol. 2009; 139: 572–583. pmid:19280672
  75. 75. Halley DJ, Rosvold J. Stable isotope analysis and variation in medieval domestic pig husbandry practices in northwest Europe: absence of evidence for a purely herbivorous diet. J Archaeol Sci. 2014; 49: 1–5.
  76. 76. Montanari M. L'alimentazione contadina nell'alto Medioevo. Napoli: Liguori Editore; 1979.
  77. 77. Montanari M. Gusti nel Medioevo. Roma-Bari: Laterza; 2012.
  78. 78. Craig OE, Bondioli L, Fattore L, Higham T, Hedges R. Evaluating marine diets through radiocarbon dating and stable isotope analysis of victims of the AD79 eruption of Vesuvius. Am J Phys Anthropol. 2013; 152: 345–352. pmid:24000142
  79. 79. Alexander MM, Gerrard CM, Gutiérrez A, Millard AR. Diet, society, and economy in late medieval Spain: Stable isotope evidence from Muslims and Christians from Gandía, Valencia. Am J Phys Anthropol. 2015; 156: 263–273. pmid:25351146
  80. 80. Alexander MM, Gutiérrez A, Millard AR, Richards MP, Gerrard CM. Economic and socio-cultural consequences of changing political rule on human and faunal diets in medieval Valencia (c. fifth–fifteenth century AD) as evidenced by stable isotopes. Archaeol Anthropol Sci. 2019;
  81. 81. Hedges RE. Isotopes and red herrings: comments on Milner et al. and Lidén et al. Antiquity 2004; 78: 34–37.
  82. 82. Fuller BT, Fuller JL, Sage NE, Harris DA, O'Connell TC, Hedges REM. Nitrogen balance and d15N: why you're not what you eat during nutritional stress. Rapid Commun Mass Spectrom. 2005; 19: 2497–2506. pmid:16106342
  83. 83. Katzenberg MA, Lovell NC. Stable isotope variation in pathological bone. Int J Osteoarchaeol. 1999; 9: 316–324.
  84. 84. Mekota AM, Grupe G, Ufer, S, Cuntz U. Serial analysis of stable nitrogen and carbon isotopes in hair: monitoring starvation and recovery phases of patients suffering from anorexia nervosa. Rapid Commun Mass Spectrom. 2006; 20: 1604–1610. pmid:16628564
  85. 85. Olsen KC, White CD, Longstaffe FJ, von Heyking K, McGlynn G, Grupe G, et al. Intraskeletal isotopic compositions (δ13C, δ15N) of bone collagen: nonpathological and pathological variation. Am J Phys Anthropol. 2014; 153: 598–604. pmid:24374993
  86. 86. Scorrano G, Brilli M, Martínez-Labarga C, Giustini F, Pacciani E, Chilleri F, et al. Palaeodiet reconstruction in a woman with probable celiac disease: a stable isotope analysis of bone remains from the archaeological site of Cosa (Italy). Am J Phys Anthropol 2014; 154: 349–356. pmid:24706415
  87. 87. Scorrano G. The stable isotope method in human paleopathology and nutritional stress analysis. Archaeology & Anthropology Open Access 2018; 1: 1–3.
  88. 88. Fogel ML, Tuross N, Owsley DW. Nitrogen isotope tracers of human lactation in modern and archeological populations. Annual Report of the Director of the Geophysical Laboratory, Carnegie Instn Washington, 1988–1989. Washington, DC: Geophysical laboratory; 1989. pp 111–117.
  89. 89. Montanari M. Storia, alimentazione e storia dell'alimentazione: le fonti scritte altomedievali. Archeol Med. 1981; 8: 25–37.
  90. 90. Iacumin P, Galli E, Cavalli F, Cecere L. C4-consumers in southern Europe: the case of Friuli V.G. (NE-Italy) during early and central Middle Ages. Am J Phys Anthropol. 2014; 154: 561–574. pmid:24889200
  91. 91. Torino M, Boldsen JL, Tarp P, Rasmussen KL, Skytte L, Nielsen L, et al. Convento di San Francesco a Folloni: the function of a Medieval Franciscan Friary seen through the burials. Heri Sci. 2015; 3: 27.
  92. 92. Leonard C, Vashro L, O'Connell JF, Henry AG. Plant microremains in dental calculus as a record of plant consumption: A test with Twe forager-horticulturalists. J Archaeol Sci Rep. 2015; 2: 449–457.
  93. 93. Radini A, Nikita E, Buckley S, Copeland L, Hardy K. Beyond food: The multiple pathways for inclusion of materials into ancient dental calculus. Am J Phys Anthropol. 2017; 162: 71–83. pmid:28105717
  94. 94. Rakić S, Povrenović D, Tešević V, Simić M, Maletić R. Oak acorn, polyphenols and antioxidant activity in functional food. J Food Eng. 2006; 74: 416–23.
  95. 95. Pignone D, Laghetti G. On sweet acorn (Quercus spp.) cake tradition in Italian cultural and ethnic islands. Genet Resour Crop Ev. 2010; 57: 1261–1266.
  96. 96. Al-Rousan WM, Ajo RY, Al-Ismail KM, Attlee A, Shaker RR, Osaili TM. Characterization of acorn fruit oils extracted from selected mediterranean Quercus species. Grasas Aceites 2013; 64: 554–60.
  97. 97. Staub PO, Casu L, Leonti M. Back to the roots: A quantitative survey of herbal drugs in Dioscorides’ De Materia Medica (ex Matthioli, 1568). Phytomedicine 2016; 23: 1043–1052. pmid:27444350
  98. 98. Hughes RE, Weiberg E, Bonnier A, Finné M, Kaplan JO. Quantifying land use in past societies from cultural practice and archaeological data. Land 2018; 7: 9.
  99. 99. Duranti M. Grain legume proteins and nutraceutical properties. Fitoterapia 2006; 77: 67–82. pmid:16406359
  100. 100. Rubiales D, Mikic A. Introduction: legumes in sustainable agriculture. Crit Rev Plant Sci 2015; 34: 2–3.
  101. 101. Copeland L, Hardy K. Archaeological Starch. Agronomy 2018; 8: 4.
  102. 102. Vitalini S, Puricelli C, Mikerezi I, Iriti M. Plants, people and traditions: ethnobotanical survey in the Lombard Stelvio national park and neighbouring areas (Central Alps, Italy). J Ethnopharmacol. 2015; 173: 435–458. pmid:26031473
  103. 103. Eglinton G, Hamilton RJ, Raphael RA, Gonzalez AG. Hydrocarbon constituents of the wax coatings of plant leaves: a taxonomic survey. Phytochemistry 1962; 1: 89–102.
  104. 104. Buckley SA, Stott AW, Evershed RP. Studies of organic residues from ancient Egyptian mummies using high temperature-gas chromatography-mass spectrometry and sequential thermal desorption-gas chromatography-mass spectrometry and pyrolysis-gas chromatography-mass spectrometry. Analyst 1999; 124: 443–452. pmid:10605875
  105. 105. Ishida M, Hara M, Fukino N, Kakizaki T, Morimitsu Y. Glucosinolate metabolism, functionality and breeding for the improvement of Brassicaceae vegetables. Breed Sci. 2014; 64: 48–59. pmid:24987290
  106. 106. Mercuri AM, Allevato E, Arobba D, Mazzanti MB, Bosi G, Caramiello R, et al. Pollen and macroremains from Holocene archaeological sites: A dataset for the understanding of the bio-cultural diversity of the Italian landscape. Rev Palaeobot Palyno. 2015; 218: 250–266.
  107. 107. Bell L, Wagstaff C. Enhancement of Glucosinolate and Isothiocyanate Profiles in Brassicaceae Crops: Addressing Challenges in Breding for Cultivation, Storage, and Consumer-Related Traits. J Agr Food Chem. 2017; 65: 9379–9403.
  108. 108. Kris-Etherton PM, Harris WS, Appel LJ. Fish consumption, fish oil, omega-3 fatty acids, and cardiovascular disease. Arterioscler Thromb Vasc Biol. 2003; 23: 20–30.
  109. 109. Swanson D, Block R, Mousa SA. Omega-3 fatty acids EPA and DHA: health benefits throughout life. Adv Nutr. 2012; 3: 1–7. pmid:22332096
  110. 110. Urashima T, Fukuda K, Messer M. Evolution of milk oligosaccharides and lactose: A hypothesis. Animal 2012; 6: 369–374. pmid:22436215
  111. 111. Borstad CM, Garvie-Lok S, Katsonopoulou D. Diet at ancient Helike, Achaea, Greece based on stable isotope analysis: From the Hellenistic to the Roman and Byzantine periods. J Archaeol Sci Rep. 2018; 18: 1–10.
  112. 112. Franzak F, Pitta D. New product development at eastern spice and flavorings. J Prod Brand Manag. 2005; 14: 462–467.
  113. 113. Piccinni G. La campagna e la città (secoli XII-XV). In Cortonesi A, Pasquali G, Piccinni G, editors. Uomini e campagne nell'Italia medievale. Roma-Bari: Laterza; 2002. pp. 125–189.
  114. 114. Guasch-Jané MR, Ibern-Gómez M, Andrés-Lacueva C, Jáuregui O, Lamuela Raventós RM. Liquid chromatography with mass spectrometry in tandem mode applied for the identification of wine markers in residues from ancient Egyptian vessels. Anal Chem. 2004; 76: 1672–1677. pmid:15018566
  115. 115. Salvini L, Pecci A, Giorgi G. Cooking activities during the Middle Ages: organic residues in ceramic vessels from the Sant'Antimo Church (Piombino–Central Italy). J Mass Spectrom 2008; 43: 108–115. pmid:17724781
  116. 116. Wu S, Tian L. Diverse phytochemicals and bioactivities in the ancient fruit and modern functional food pomegranate (Punica granatum). Molecules 2017; 22: 1606.
  117. 117. Bhowmik D, Gopinath H, Kumar B, Duraivel S, Aravind G, Kumar K. Medicinal uses of Punica granatum and its health benefits. J Pharmacogn Phytochem 2013; 1: 28–35.
  118. 118. Wadhawan R, Reddy Y, Sharma R, Solanki G. Punica granatum; ancient seeds for modern cure in medical and dental field: An overview. Int J Trad System Med 2015; 2: 1–8.
  119. 119. Connelly E, del Genio CI, Harrison F. Datamining a medieval medical text reveals patterns in ingredient choice that reflect biological activity against the causative agents of specified infections. New York: Cornell University Editor; 2018. Available from arXiv preprint arXiv:1807.07127.
  120. 120. Jadhav RS, Lokhande SS, Nikam RJ. Herbal mouth wash: an update review. WJPPS 2018; 7: 436–445.
  121. 121. Thangavelu A, Elavarasu S, Sundaram R, Kumar T, Rajendran D, Prem F. Ancient seed for modern cure–pomegranate review of therapeutic applications in periodontics. J Pharm Bioallied Sci. 2017; 9: S11. pmid:29284927
  122. 122. Liu A, Lee SM, Wang Y, Du G. Elsholtzia: review of traditional uses, chemistry and pharmacology. J Chin Pharm Sci 2007; 16: 73.
  123. 123. Limberger RP, Jacques ALB, Schmitt GC, Arbo MD. Pharmacological effects of ephedrine. Berlin: Nat. Prod. Springer; 2013.
  124. 124. Desrosiers MR, Weathers PJ. Effect of leaf digestion and artemisinin solubility for use in oral consumption of dried Artemisia annua leaves to treat malaria. J Ethnopharmacol. 2016; 190: 313–318. pmid:27340104
  125. 125. Fatucci M, Cerilli E. Allevamento e caccia al castello di Santa Severa (Roma) durante il basso medioevo: mense signorili e popolari dal XIII al XIV secolo. Sezione di Museologia Scientifica e Naturalistica 2016; 12: 235–242.
  126. 126. Gowland RL, Garnsey P. Skeletal evidence for health, nutritional status and malaria in Rome and the empire. J. Rom. Arch. Suppl. Series 2010; 78: 131–156.
  127. 127. Gowland RL, Western AG. Morbidity in the marshes: Using spatial epidemiology to investigate skeletal evidence for malaria in Anglo‐Saxon England (AD 410–1050). Am J Phys Anthropol. 2012; 147: 301–311. pmid:22183814
  128. 128. Buckley HR, Tayles N. Skeletal pathology in a prehistoric Pacific island skeletal sample: issues in lesion recording, quantification and interpretation. Am. J. Phys. Anthropol. 2003; 122: 303–324. pmid:14614753
  129. 129. Buckley HR. ‘The predators within’: investigating the relationship between malaria and health in the prehistoric Pacific Islands. In: Oxenham M, Tayles N, editors. Bioarchaeology of Southeast Asia. Cambridge: Cambridge University Press; 2006. pp. 309–332.
  130. 130. Brothwell DR. Digging up bones. Ithaca: Cornell University Press; 1981.
  131. 131. Ortner DJ. Identification of pathological conditions in human skeletal remains. San Diego: Academic Press; 2003.
  132. 132. Walker PL, Bathurst RR, Richman R, Gjerdrum T, Andrushko VA. The causes of porotic hyperostosis and cribra orbitalia: a reappraisal of the iron-deficiency-anemia hypothesis. Am. J. Phys. Anthropol. 2009; 139: 109–125. pmid:19280675
  133. 133. Ulijaszec SJ, Hillman G, Boldsen JL, Henry CJ. Human dietary change [and discussion]. Phil T R Soc B. 1991; 334: 271–279. https://doi.org/10.1098/rstb.1991.0116
  134. 134. Papathanasiou A. Health status of the Neolithic population of Alepotrypa Cave, Greece. Am J Phys Anthropol. 2005; 126; 377–390. pmid:15386282
  135. 135. Larsen CS. The agricultural revolution as environmental catastrophe: implications for health and lifestyle in the Holocene. Quat Int. 2006; 150: 12–20.
  136. 136. Eshed V, Gopher A, Pinhasi R, Hershkovitz I. Paleopathology and the origin of agriculture in the Levant. Am J Phys Anthropol 2010; 143: 121–133. pmid:20564538
  137. 137. Taylor GM, Rulland P, Molleson T. A sensitive polymerase chain reaction method of the detection of plasmodium species DNA in ancient human remains. Ancient Biomolecules 1997; 1: 193–203.
  138. 138. Nerlich AG, Schrout B, Dittrich S, Thomas J, Zink AR. Plasmodium falciparum in Ancient Egypt. Emerg Infect Dis. 2008; 14: 1317–1319. pmid:18680669
  139. 139. Sallares R, Bouwman A, Anderung C. The spread of malaria to Southern Europe in antiquity: new approaches to old problems. Med Hist. 2004; 48: 311–328. pmid:16021928
  140. 140. Petrovska BB. Historical review of medicinal plants’ usage. Pharmacogn Rev. 2012; 6: 1–5. pmid:22654398
  141. 141. Burdock GA, Carabin IG, Soni MG. Safety assessment of β-nitropropionic acid: a monograph in support of an acceptable daily intake in humans. Food Chem. 2001; 75: 1–27.
  142. 142. Maruo VM, Bracarense AP, Metayer JP, Vilarino M, Oswald IP, Pinton P. Ergot Alkaloids at Doses Close to EU Regulatory Limits Induce Alterations of the Liver and Intestine. Toxins 2018; 10: 183.
  143. 143. American Dental Association. Current Dental Terminology, 2005, CDT-2005, Chicago, IL: American Dental Association, 2004.
  144. 144. Ricci P, Mongelli V, Vitiello A, Campana S, Sirignano C, Rubino M, et al. The priviledged burial of the PavaPieve (Siena, 8th Century AD). Rapid Commun Mass Spectrom. 2012; 26: 2393–2398. pmid:22976205
  145. 145. Pescucci L, Battistini A, De Angelis F, Catalano P. Vivere al centro di Roma nell'VIII secolo d.C. Indicazioni antropologiche. Bollettino di Archeologia On-Line 2013; 4: 2–4.