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

Who were the miners of Allumiere? A multidisciplinary approach to reconstruct the osteobiography of an Italian worker community

  • Marica Baldoni ,

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

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

    Affiliations Laboratorio di Antropologia Forense e Biologia dello Scheletro, Dipartimento di Biologia Università degli Studi di Roma “Tor Vergata”, Roma, Italia, Laboratorio di Medicina Legale, Dipartimento di Biomedicina e Prevenzione, Università degli Studi di Roma “Tor Vergata”, Roma, Italia

  • Gabriele Scorrano ,

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

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

    Current address: Natural History Museum of Denmark; University of Copenhagen, Copenhagen, Denmark.

    Affiliation Centro di Antropologia Molecolare per lo Studio del DNA antico, Dipartimento di Biologia Università degli Studi di Roma “Tor Vergata”, Roma, Italia

  • Angelo Gismondi ,

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

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

    Affiliation Laboratorio di Botanica, Dipartimento di Biologia Università degli Studi di Roma “Tor Vergata”, Roma, Italia

  • Alessia D’Agostino,

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

    Affiliation Laboratorio di Botanica, Dipartimento di Biologia Università degli Studi di Roma “Tor Vergata”, Roma, Italia

  • Michelle Alexander,

    Roles Formal analysis, Methodology, Writing – review & editing

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

  • Luca Gaspari,

    Roles Formal analysis, Methodology, Writing – review & editing

    Affiliation Centro di Antropologia Molecolare per lo Studio del DNA antico, Dipartimento di Biologia Università degli Studi di Roma “Tor Vergata”, Roma, Italia

  • Fabrizio Vallelonga,

    Roles Data curation, Methodology, Resources, Writing – review & editing

    Affiliation Sapienza Università di Roma, Dipartimento di Scienza dell'Antichità, Cattedra di Archeologia Cristiana e Medievale, Roma, Italia

  • Antonella Canini,

    Roles Conceptualization, Resources, Writing – review & editing

    Affiliation Laboratorio di Botanica, Dipartimento di Biologia Università degli Studi di Roma “Tor Vergata”, Roma, Italia

  • Olga Rickards,

    Roles Conceptualization, Resources, Writing – review & editing

    Affiliation Centro di Antropologia Molecolare per lo Studio del DNA antico, Dipartimento di Biologia Università degli Studi di Roma “Tor Vergata”, Roma, Italia

  • Cristina Martínez-Labarga

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

    martine@uniroma2.it

    Affiliations Laboratorio di Antropologia Forense e Biologia dello Scheletro, Dipartimento di Biologia Università degli Studi di Roma “Tor Vergata”, Roma, Italia, Centro di Antropologia Molecolare per lo Studio del DNA antico, Dipartimento di Biologia Università degli Studi di Roma “Tor Vergata”, Roma, Italia

Abstract

This research presents an in-depth study of the skeletal remains collected from the archaeological site of Allumiere (15th-16th centuries CE; Rome, Italy). A multidisciplinary approach was used, combining skeletal biology, molecular anthropology and archaeobotany with the aim of reconstructing the osteobiography of the alum miners buried at the site. Since 1460, the area of the Tolfa Mountains was significant for the exploitation of alum which was used for a wide range of purposes in the Middle Ages, ranging from woven production to medical practice. A total of 70 individuals (63 adults and 7 juveniles) were studied. The sex ratio of the community indicated a higher prevalence of males with respect to females. Morphological examination indicated occupational musculoskeletal stress markers, which might reflect the specific phase of alum production that each individual was occupied in. Dietary reconstruction was primarily performed through carbon and nitrogen stable isotope analysis with integration of the results obtained by microscopic, genetic and GC-MS investigations on dental calculus. The diet was omnivorous, indicating a reliance on C3-terrestrial protein and evidence for limited C4 consumption by some individuals. Herbivores, such as sheep and cattle, appear to have contributed to the diet more than pigs and chickens. Consumption of Fagaceae and Poaceae species was predominant; moreover, indicators of Brassicaceae and milk and its derivatives were abundantly recurrent in the population, followed by plant oils and theophylline. Furthermore, the detection of pharmacological alkaloids indicated the knowledge and application of medicinal plants by the community. The novel use of multiple techniques based on cutting-edge technologies has provided a unique window on the lifestyles of individuals from one of the first Italian settlements of alum workers.

Introduction

The exploitation of alum in Italy started in 1460, when Giovanni di Castro, a commissioner of the Pontifical State, identified the presence of alunite, a mineral from which alum could be extracted, in the territory of the Tolfa Mountains (Rome, Italy) [13]. Alum is a salt made up of ammonium sulfate and potassium associated with twenty-four molecules of crystallization water, whose applications ranged from textile production to medicine. Because of its water solubility, alum is not directly available in nature, but it is obtained through transformation of the less soluble aluminum minerals, such as alunite. The discovery of alunite on the Tolfa Mountains was of vital importance for the papal coffers and influenced the rise of the Western European textile industry. This source became particularly important after the conquest of Constantinople by the Turks in 1453 meaning that alunite deposits located in the Eastern Mediterranean area became difficult to access [2]. The extraction of the alum was entrusted to companies of contractors while the Apostolic Chamber handled its marketing. From the beginning, the Medici family secured the control of commercialization of the product in collaboration with Genoese merchants. In 1499, the banker Agostino Chigi was responsible for the organization of the mining enterprise and settlement of the area and the birth of the village that would later become Allumiere, is likely down to his actions [2].

Extraction and processing of the alum required the involvement of specialized manpower. The production cycle of alum, described in detail by historical sources [1, 3], began with the excavation of the mineral from the rocks through the use of picks. Then, alunite stones were heated at high temperature (i.e. 600–700°C) in special furnaces for 12–14 hours. The "roasted" mineral was finally treated with water, in order to obtain a doughy solution, which was then heated again and concentrated until alum crystals began to separate.

In the area of La Bianca, several excavation campaigns have unearthed a Medieval cemetery and a church named Cappella dei Minatori (Fig 1). It has a single nave, is East-West oriented and is 19 m long and 8 m wide. The structure and dating of the church suggest that the archaeological site was related to one of the first human settlements in the area associated with alum production. A total of 70 burials were found, most of them located outside the church (Fig 1).

thumbnail
Fig 1.

A) Map of the geographical area of La Bianca B) archaeological plan of the church named Cappella dei Minatori with the position of the skeletal remains.

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

Archaeologically, the stratigraphic sequence of the site is challenging to define, as is common in Medieval cemeteries. In the majority, of the taphonomic observations, the burials were made up of simple earthen graves with individuals wrapped in shrouds, with the exception of a wooden coffin identified through the presence of in situ nails (SU 296). A chronological sequence of burials was indicated by earlier burials were being laid in a N-S orientation and subsequent later individuals buried on an E-W axis.

Items of jewelry and coins were recovered with a few individuals. In SU 110, three silver coins were found close to the left hand; SU 135 was buried with a religious medal, portraying the Immaculate Conception and SU 245 had eleven gold and two silver coins close to the left hand. All these artifacts confirmed the dating of the cemetery; in particular, the silver coins of SU 110 were dated to 1476–1483 [4]. Further archaeological analyses revealed various origins for these coins. Some coins possessed a “fleur-de-lis”, symbol of the city of Florence (Italy). Others were also identified as three ducats of Ferdinando I from Naples, one coin from the Pontifical State and three from Venice. Others were of Iberian origin, specifically, a doble excelente de la Granada and an excelente of the kings Fernando de Aragon and Isabel de Castilla and a Portuguese cruzado of King Giovanni II d'Aviz. There was also a duchy from Rhodes. The dating of the coins ranged between 1464–1523. The variety of coins and their origins are an important proof of the wealth of commercial traffic linked to alum.

The present research aims to reconstruct the osteobiography of these alum miners of the Tolfa Mountains through a multi-proxy approach, combining skeletal biology, molecular anthropology and archaeobotany. Osteological analysis allows the reconstruction of the biological features of the Medieval workers, while molecular and archeobotanical analyses are used to shed light on the dietary patterns and medicinal habits of the community of Allumiere.

There is a general lack of information about the dietary habits of the lower social classes in historical texts from the Medieval period. Therefore, to better understand the dietary pattern of the studied community, we applied carbon (δ13C) and nitrogen (δ15N) stable isotope analysis to human and animal bone collagen. Stable isotope analysis of bone collagen represents an established technique to identify the main protein sources in the diet of archaeological populations. This allows us to investigate marine and terrestrial and C3 and C4 sources of protein in addition to the trophic position an individual is feeding at [57]. Analysis of dental calculus was also carried out to further investigate which vegetal and animal species were present in miners’ diet. Indeed, the analysis of dental calculus is a very informative archaeobiological approach to reconstruct past food habits [810]. Dental calculus is a dense mineral matrix made up of inorganic salts and organic molecules derived from ingested foods, crystallized remains of oral microbiota and accidentally inhaled microremains [1112]. Light microscopy, genetic and gas-chromatographic mass-spectrometry approaches were carried out on dental calculus to provide direct evidence of the main categories of foods and plant drugs introduced, at least once, in the lifetime of our ancient exemplars [1323]. The present work provides new information on an important period of the Italian Medieval history, shedding light on the lives of those who lived in one of the first Italian settlement of alum miners.

Materials and methods

The present research deals with the skeletal remains recovered in the area of La Bianca in Allumiere (Rome, Italy) (Fig 1A) in a graveyard close to the church named as Cappella dei Minatori (Fig 1B). The archaeological excavation started in 2010 directed by Dr. Fabrizio Vallelonga authorized by the “Comune di Allumiere” (Municipality of Allumiere). The research was carried out at the Department of Biology of the University of Rome “Tor Vergata” and directed by Dr. Cristina Martínez-Labarga who received the authorization for the analysis of the skeletal remains from La Bianca (Allumiere, Rome, Italy) in 2015. The complete list of the specimens is provided in S1 Table.

Osteological examination

Seventy individuals were analysed, both adults and non-adults collected from the cemetery area of La Bianca (Allumiere, Rome, Italy). The preservation index for the individuals was calculated following the method proposed by Walker et al. [24]. The age estimation for adult individuals (from ca. 18 year old) followed methods based on morphological changes in the pubic symphysis [2527], in the auricular surface of the ilium [28] and in the sternal end of the fourth rib [2930]. Secondarily, dental wear [3132] and obliteration of the cranial sutures [33] were also observed. Age estimation for infant and juvenile skeletal remains (until ca. 18 years old) was carried out through the diaphyseal length of bones [3436], and tooth eruption [37]. Moreover, secondary ossification centers were taken into account, following the methods proposed by different authors and summarized in Minozzi and Canci [38].

Sex diagnosis was performed only on adult samples, as proposed by Acsàdi and Nemeskèri [39] and revised by Ferembach et al. [40] and Phenice’s [41]. Metric analysis complemented the morphological data, as the state of preservation of the sample allowed the measurement of sexually dimorphic bones, utilizing univariate and multivariate techniques [4246].

Osteometrics was applied following the methods and standards proposed by Martin and Saller [47] and Borrini [4849]. Cranial and post-cranial indices were calculated following the guidelines proposed by several authors [38, 43, 5051]. Living stature was estimated as the average value obtained using multiple methods [5256]. When skeletal remains were fragmented, Steele’s formulae [57] were applied to estimate long bone length. Moreover, working activities were determined analyzing enthesal changes (EC) as proposed by Mariotti et al. [5859] and by Borgognini Tarli and Reale [60]. The paleopathological survey was performed through morphological observation of the skeletal remains. In one pathological case, X-ray monitoring was carried out at the “Dipartimento di Diagnostica per Immagini, Imaging molecolare, Radiologia, Interventistica e Radioterapia, Azienda Ospedaliera Universitaria Policlinico Tor Vergata”, using GMM Opera Swing, in order to achieve a differential diagnosis.

Particular attention was paid to evidence of osteoarthritis, which was scored according to the standards proposed by Buikstra and Ubelaker [61]. These authors evaluated intensity and distribution of lipping, porosity and eburnation of the osteoarthritis. A scale system from 0 -absence of injury on bone tissue- to 3 -marked modification of the joint surface- was used. In order to better evaluate the position in which the lesion occurred, Prieto’s joint division scheme was adopted [62]. This method distinguishes proximal and distal epiphyses, dividing them in four areas: two superiors (anterior and posterior) and two inferiors (anterior and posterior).

Statistical analyses were carried out using the statistical software R (v. 3.4.1) for the Student’s t-test and chi square test with Yates correction.

Isotope analyses

Collagen extraction was carried out on 68 human and 12 faunal skeletal remains from the archaeological cistern (1 Hystrix cristata, 1 Cervus elaphus, 2 Equus asinus, 2 Felix catus, 1 Canis lupus 2 Sus domesticus, 1 Equus caballus and 2 unidentified carnivores) in order to conduct carbon and nitrogen stable isotope analysis. Human samples comprised of rib bones, whereas animal samples were taken from the various skeletal elements that were available.

Collagen extraction followed a modified Longin method [63]. Initially, to remove potential contaminants a sterile surgical blade was used on the outer surface of the bone samples and c. 500 mg of bone was subsequently pulverized using a mortar and pestle. To demineralize the bone, 8 mL of HCl 0.6 M at 4°C was added to the powder and left at 4°C on a horizontal mixer for two days, changing the acid after 24 hours.

Once the mineral component of the bone removed, samples were rinsed three times with ddH2O, until the pH level became neutral. The resultant pellet was gelatinized at 75°C, for 24–48 hours, with HCl pH 3.0 (0.001 M). The solution was then frozen at -80°C for four hours and then freeze-dried for two days. A simultaneous extraction on modern bovine bone was performed and used as reference control. Approximately 0.8–1.2 mg of collagen was weighed and analysed in duplicate by EA-IRMS on a Sercon GSL analyser coupled to a Sercon 20–22 Mass Spectrometer at the University of York. The analytical error, calculated from repeated measurements of each sample, an internal laboratory control (fish gelatine), and international standards, was <0.2‰ (1σ) for both δ13C and δ15N. International Atomic Energy Agency (IAEA) standards were N-2, and 600 for nitrogen and International Atomic Energy Agency IAEA-600 and Iso-Analytical R006 for carbon. Isotope data are reported as delta (δ) values relative to V-PDB (Vienna Pee Dee Belemnite) for carbon and AIR (Atmospheric air) for nitrogen. Carbon content (%C), nitrogen content (%N), protein yield, and C/N ratios were checked to monitor the diagenesis of bone [64] and determine bone protein quality for paleodietary reconstruction, according to DeNiro [65] and van Klinken [66].

Molecular and archaeobotanical analyses on dental calculus

Sample collection and decontamination.

Dental calculus study was performed on 35 individuals, who were the only ones presenting maxillary and mandibular bones and teeth. Each analysis was carried out on different aliquots of calculus. Sodium hypochlorite 5% and UV exposure were employed to treat tools and working surfaces to restrict contamination during the analysis. In addition, different laboratories were used to carry out the various phases of this work [67]. Deposits of dental calculus were generally slight, according to Brothwell [31], on all dentition. Initially, to eliminate any contamination before sampling, a scraping action was applied by a sterile surgical blade on the outer surface of the dental calculus. All collecting procedures were carried out under a sterile vertical laminar flow hood (Heraeus HERAsafe HS12 Type). In order to remove any environmental contaminants from the calculus surface, samples were UV-treated for 10 minutes on each side and soaked in 5% sodium hypochlorite for 15 minutes. Lastly, the calculus was washed in sterilized bidistilled water and rinsed in 100% ethanol to eliminate the aqueous components before being left to dry [22, 68]. To validate sterilization protocols, five dental calculus samples previously subjected to decontamination procedures, were randomly selected and washed with 200 μL of water. These last washing solutions were subjected to light microscopy, genetic and GC-MS analysis. No microremains, nucleic acid or chromatographic signals were detected, confirming the efficacy of the cleaning methods.

DNA extraction, amplification and sequencing.

All criteria and precautions for the study and analysis of ancient DNA (aDNA) were applied [6971]. For aDNA extraction a modified protocol suggested by Warinner and collaborators [10] was used. For each sample (50 mg of pulverized dental calculus), 600 μL of extraction buffer (100 mM Tris-HCl pH 8, 100 mM NaCl, 10 mM EDTA and 2% SDS) and 50 μL of proteinase K (20 mg/mL) were added. The samples were incubated in a shaking water bath at 56°C, for 6 hours, adding 20 μL of fresh proteinase K (20 mg/mL) every 2 hours. The sample was then incubated overnight at 37°C. Following this, the sample was centrifuged for 5 minutes at 13.000 rpm and the supernatant was transferred to a new 2 mL Eppendorf tube. The supernatant was mixed with 500 μL of phenol/chloroform/isoamyl alcohol (25:24:1) and, after a centrifugation of 5 minutes at the 13.000 rpm, transferred into a new 15 mL falcon tube. DNA was purified by QIAquick PCR purification kit following the manufacturer’s procedure, and eluted into 50 μL of the elution buffer. Ancient DNA extracts were stored at 4°C.

For each sample, different polymorphic regions of mitochondrial genome of several animal species were amplified (35 PCR cycles), using pairs of primers that only amplified for the target species (S2 Table). For ovine and pig, the COX 1 gene was amplified, while 12S rRNA and 16S rRNA were analysed for bovine and chicken, respectively. These primers were proposed by Natonek-Wiśniewska and colleagues [72]. For fish investigations, the DNA Mini-Barcoding System reported in Shokralla and collaborators [73] was applied. Amplifications were also performed for negative and positive controls. Amplifications of modern and ancient DNA were performed separately in two different laboratories. In particular 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) which has all the facilities to minimize potential contamination with extant DNA [71]. In order to detect the correct activity of the different reagents, the positive controls that contained modern faunal DNA were prepared in the laboratory for modern DNA processing in the department of Biology, of the University of Rome “Tor Vergata”.

To verify the amplifications, PCR products were separated by electrophoresis on 1.5% agarose/Tris-acetate EDTA gel and the amplicons were then purified by enzymatic digestion with HT ExoSAP-ITVR (Affymetrix) following the manufacturer’s procedure and sequenced using the same PCR primers by an ABI PRISM 3130 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA), as reported in Gismondi et al. [74]. The BLAST website tool (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastn&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome) was used to confirm the identity of the sequences found in the dental calculus. The presence of fish DNA was detected by polyacrylamide gel, using a molecular weight marker (GelPilot 50bp Ladder Qiagen).

Light microscopy analysis (LM).

A protocol based on Hardy et al. [75] method was used to extract starch granules and other microremains from dental calculus. Twenty mg of each sample were resuspended in 500 μL of 1 M HCl and sonicated for 10 minutes (Falc Instruments MOD: LBS1 34) before being left under agitation for 24 h. After centrifugation at maximum speed for 10 minutes, the pellet was subjected to 3 consecutive washings with bidistilled water. The last pellet was resuspended in 100 μL of bidistilled water, sonicated for 5 minutes and observed at OM (Nikon Eclipse E100), under white and polarized light. The whole volume of the sample was examined and microfossils were photographed at 100X magnification using software for capturing images (ProgRes CapturePro 2.9.0.1) and measured (i.e. for starch granule, maximum length through the hilum and maximum width anywhere along the perpendicular axis) by the SuperAmpelo 2.0 program. Taxonomic identification of pollens, phytoliths and starch granules was carried out by direct comparison with a modern experimental collection (hosted in Laboratory of Botany and Botanical Gardens of the Department of Biology at University of Rome “Tor Vergata”) or data from literature, including a pollen Atlas [76].

Gas-chromatography mass-spectrometry (GC-MS) analysis.

This analysis was carried out in qualitative and non-quantitative terms, according to Gismondi et al. [77]. Ten mg of calculus for each individual were dissolved in 1 mL of 6% HCl. Once left in agitation for three days, 1 mL of hexane was added and shaken for two hours. After centrifugation for 10 minutes at 10.000 rpm, the supernatant hexane fraction was recovered and dried by a speed-vac system (Eppendorf AG 22331 Hamburg, Concentration Plus). The dried pellet was then derivatized by resuspension with 50 μL of hexane and 50 μL of the Methyl-8-Reagent (Thermo Scientific), in thermostated bath, at 60°C for 20 minutes. At least 3 analyses for each sample, by injecting 2 μL of extract in GC-MS (QP2010, Shimadzu, Japan; column DB-5 Phenomenex; helium as gas carrier; splitless modality), were carried out. The run was conducted by a temperature gradient: initial oven temperature was set at 60°C for 5 min; then, increasing temperature at a fixed rate of 6°C per minute, the column was heated up at 150°C for 5 minutes, 250°C for 5 minutes and 330°C for 25 minutes. Parameters and conditions relating to mass spectrometry were: ion source temperature 230°C; interface temperature 320°C; solvent cut time 6 minutes; ionization mode: EI; ionization voltage: 0.70 eV. The identification of each molecule (similarity values were considered acceptable only if higher than 85%) was performed comparing their mass spectrum with those registered in the NIST Library 14 loaded on detection software. No significant differences among replicates were detected.

Results

Skeletal biology

The present series of Medieval samples consisted of 70 individuals, 63 adults and 7 non-adults, sampled from primary burials in the cemetery area of La Bianca. The preservation index was calculated for each skeleton (ranging from a minimum of 0 to a maximum of 93%), showing an average preservation value of 59.4%. Individuals were sorted into 9 age classes. The composition by sex and age at death is summarized in Table 1. It is worth pointing out the under-representation of sub-adults with respect to adults and the high prevalence of males with respect to females among the adults. The imbalance between sexes is indicated by the elevated sex ratio (M:F = 8.0) which is significantly different from the range of expected proportions (χ2, p<0.05).

thumbnail
Table 1. Sex and age at death determination in the analyzed series from Allumiere (Rome. Italy).

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

For adult individuals, the highest mortality falls between the young adult and adult classes. Data on the juveniles’ mortality rate data should be considered cautiously due to the small sample size.

Living stature was estimated on 43 individuals (37 males and 6 females). The mean stature for males was 169.22 ± 6.26 cm, while females exhibited an average stature of 156.70 ± 5.96 cm. A Student’s t-test indicated a statistically significant difference between the sexes (p-value = 3.17x10-7). Post-cranial indices revealed high skeletal robusticity and skeletal asymmetry in males. In particular, the right side seemed to have been subjected to a higher biomechanical stress than the left. The results obtained for all the indices are shown in S3 Table. Pathological assessment was carried out on 57 adult individuals (S4 Table). Six individuals were excluded due to poor preservation.

The highest frequencies were recorded for degenerative diseases (86%), Schmörl’s nodules (75%), periostitis (32%), and fractures (27%). Degenerative diseases were considered, jointly with enthesal changes, able to associate each individual to the most plausible phase in alum production. Fractures were relatively common within the community. For example, SU 245 demonstrated an ante-mortem compound fracture which affected both tibia and fibula. This trauma (Fig 2) caused the loss of ca. 5 cm in limb length.

thumbnail
Fig 2. Ante-mortem compound fracture affecting SU 245.

A) compound fracture on both tibia and fibula; B) X-rays analysis (X-ray exposure: 78 kV).

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

The analysis of musculoskeletal stress markers was conducted on male individuals as historical evidence indicates [1,3] that they would have been involved in alum extraction and production. No detection attempt was made on juveniles because the high plasticity of their skeletal remains could provide unreliable results. Only 20 male specimens were considered suitable for all further analyses on enthesal changes. The remainder were excluded either due to their state of preservation, which prevented macroscopic observation of muscular insertion sites, or because, as was the case for SU 245, a trauma or skeletal disease was present that could alter interpretation. Indeed, the musculoskeletal stress markers of SU 245 seemed to be influenced by skeletal lesions. The enthesal changes observed in the upper limbs of this individual (i.e. a greater muscle development on the right humerus than on the left) led us to hypothesize that they would have used a crutch to partially solve the walking difficulty caused by the loss of the original length of the limb as consequence of the ante-mortem fracture.

On the basis of the development of occupational stress markers, the individuals seemed to cluster in four groups, corresponding to the four different phases of alum production (excavation, calcination, moistening and lixiviation) [78]. These data were also corroborated by the results of the post-cranial indices, which revealed different patterns of skeletal strength as well as a functional asymmetry. Musculoskeletal stress data were finally integrated with the analysis of degenerative pattern and joint diseases. Osteoarthritis could be observed on 45 samples (S5 Table).

Osteoarthritis confirmed all previous suppositions about musculoskeletal stresses and alum working activity of the individuals. The excavation phase caused a marked enthesal development of shoulders, elbows, hips and feet; in the calcination the most involved functional groups are shoulder, elbow, hip and knee; the individuals involved in moistening exhibit a high degree of biomechanical stress that affected shoulders, hips, knee and feet, while those employed in lixiviation show a great muscle development in the shoulders, forearms, hip and foot.

All the working phases are characterized by well-developed musculoskeletal markers on the shoulders and the hip however, the combination of other functional groups significantly diverges amongst the other clusters. For example, the skeletons that show high robusticity of m. pectoralis major (horizontal flexion of the upper limb), m. deltoideus (abduction, extension and flexion of the limb) and m. biceps brachii (flexion of the forearm on the arm) could be related to the excavation phase according to the upper limb’s movements involved in the use of the picks. Even if, as mentioned, all phases of alum production cause a high biomechanical stress on muscles, the movement related to each one as well as the functional groups involved are different, and the excavation was observed to be the most strenuous phase involving the strain on both the shoulder and the elbow [1, 3].

Furthermore, a high degree of biomechanical stress on the knee could be observed in both moistening and calcination phases. However, moistening should involve a higher robusticity of the quadriceps tendon, whereas in the calcination m. vastus medialis should appear more distinctly stressed. The combined results of musculoskeletal stress markers and degenerative diseases are provided in Table 2.

thumbnail
Table 2. Combined results of musculoskeletal stress markers and osteoarthritic pattern (nd: not determinable).

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

Unfortunately, a clear relationship between musculoskeletal markers, osteoarthritis and performed working activity could be identified only for twenty-one individuals because the others exhibited musculoskeletal stress markers or osteoarthritis patterns that could made them suitable for more than one working phase.

Stable isotope analysis from bone collagen and molecular analysis from dental calculus

Isotopic data and quality indicators are presented in Table 3.

thumbnail
Table 3. Carbon and nitrogen stable isotope values and collagen quality indicators of animals and humans from La Bianca.

The samples that fell outside the quality range parameters are indicated in red and were excluded from subsequent analysis.

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

Ten humans (SU 213, 330, 198, 325, 320, 307, 179, 385.3, 385.5, 416) and one animal (SU 101 A Equus caballus), were excluded from the analysis due to their C/N ratio being outside the satisfactory range. The remaining samples were acceptable according to criteria proposed by Ambrose and Norr [64] and van Klinken [66]. Fig 3 shows the plot of δ13C versus δ15N isotope values for both faunal and human samples.

thumbnail
Fig 3. Plot of carbon vs nitrogen isotopic values for faunal and human remains from Allumiere.

Human data are plotted individually and distinguished by age (adults and non-adults) and sex. Animal data are plotted as means and error bars represent ±1σ.

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

The δ13C values of the three herbivore samples (Cervus elaphus and Equus asinuus) ranged from -21.4‰ to -21.1‰ (mean -21.3 ± 0.1‰), while δ15N values ranged between 3.5‰ and 5.9‰ (mean 4.5 ± 1.3‰). Omnivores (Sus domesticus and Hystrix cristata,) possessed values from 21.8‰ to -20.6‰ (mean -21.0 ± 0.7‰) for δ13C and from 3.6‰ to 4.7‰ (mean 4.1 ± 0.6‰) for δ15N. The carnivore (Felix catus, Canis familiaris and unidentified carnivores) data ranged from -20.8‰ to -19.5‰ (mean -20.1 ± 0.46‰) and from 7.5‰ to 8.5‰ (mean 7.8 ± 0.4‰) for δ15N. Isotope data for all animals are in a range typically associated with the consumption of C3 plants. In terms of nitrogen, herbivores have δ15N values within the range found for primary consumers in the Mediterranean area [79]. Omnivore δ15N values are similar to those of the herbivores, which indicates that at Allumiere they consumed similar dietary resources, i.e. mainly plants with little/no animal protein. The Hystrix cristata (porcupine) show isotope values that are of an essentially herbivorous and it seems that any occasional consumption of insects and small vertebrates [80] did not play a key role in the stable isotope values for this animal. The analysis of the pig diet is useful to understand the variation in husbandry practices [8188] among the different communities. Indeed, pigs raised in a home-based system are expected to have a more controlled diet than free-range animals. For Medieval northwestern Europe communities, pigs had a diet based on terrestrial plant and human refuse [87], whereas our data show a comparable δ15N values between the Sus domesticus and the other herbivores analyzed. During the Middle Ages, in some areas (including Lazio) pigs were free to roam in the fallow land close to the city [82], where they would typically consume acorns (pannaga) [89]. Our results seem to suggest that the Allumiere pigs were predominantly fed with plant products, which might indicate that they wandered about in forested environments. Therefore, the difference between Allumiere and data from published northwestern European populations [87] is probably related to the different social status of the communities: the Allumiere population was made up of workers of low social status, while the other data derives from predominantly coastal and elite/urban communities with different husbandry practices.

The human remains from La Bianca (Allumiere) possessed δ13C values between -20.6‰ and -17.1‰ (mean -19.2 ± 1‰) and δ15N values between 5.9‰ and 12.1‰ (mean 8.6 ± 1.4‰). Both nitrogen and carbon stable isotope values for humans showed a high variability, in particular two individuals (SU 245 and 317) have elevated δ15N values of 12.1‰ and 11.9‰ and one (SU 388) has a particularly low δ15N value of 5.9‰. Despite the variability in isotopic values, no statistically significant difference was observed, between different sexes or age at death. The majority of humans demonstrate an enrichment of 15N (about 3–5‰, with a nitrogen offset equal 3.5‰) in comparison to the animals (Fig 3), reflecting a typical trophic level shift indicative of the consumption of animal protein. Most of the individuals have diets with low input of animal protein: the wide range in δ15N values for the population (6.2‰) indicates that animal protein made a greater or lesser contribution to the diet of particular individuals. On average, the human-animal offset in is 3.5‰, which is at the lower end of the accepted 3–5‰ between trophic levels [90] and spacing of up to 6‰ has been suggested for humans [91]. Therefore, for many individuals, it seems animal protein did not contribute a major part of the diet, which is in keeping with the low status nature of the population. However, for those possessing δ15N values of upwards of 9–10‰, animal protein will have had a major input and for the most elevated values, fish may also have made a contribution. High δ13C values (> 18‰) for some individuals (e.g. SU 362, 182, 192 and 319) indicate the consumption of C4 plants or marine fish consumption. However, the fact that these individuals also tend to have amongst the lowest δ15N values would suggest that the former is most likely.

Comparisons of Allumiere with other Medieval published data [9299] were performed (Fig 4). Significant statistical differences (Kruskal-Wallis test p<0.05), were found between δ13C values, mostly with northeastern Italian sites (Cividale, Mainizza, Romans d’Isonzo, Siena), Cosa from entral Italy and Montella (AV) from Campania in southern Italy. For northeastern and central Italy sites, the consumption of C4 or marine resources were interpreted to be responsible for enriched δ13C values at these sites [9596]. As regards to Montella [96] the sample was from a Franciscan friary so it is possible to hypothesize that they followed a different dietary plan. Statistically significant differences (Krustal-Wallis test p<0.05) were found between δ15N with the sites of Siena and Montella. However, only the site of Siena [93] showed higher mean δ15N values, although this could also be due to the small sample size (N = 19). These differences describe a heterogeneous dietary landscape of the Medieval Italian communities probably related to geographical position and chronology.

thumbnail
Fig 4. Plot of carbon vs nitrogen mean isotopic values and error bars represent ±1σ for human remains from Allumiere and the other Medieval Italian sites.

Triangles represent the north Italian populations and circles the central Italian populations. Trino Vercellese (Piedmont) [92]; Pava Piene Siena (Tuscany) [93]; Piazza Madonna di Loreto (Rome) [94]; Cividale, Romans d'Isonzo and Mainizza (Friuli-Venezia Giulia) [95]; Cosa Grosseto (Tuscany) [96]; Albano Laziale Rome (Lazio) [97]; Montella Avellino (Campania) [98]; Colonna Rome (Lazio) [99].

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

Analysis of aDNA from dental calculus was carried out to attempt to determine which animals were relied upon and to identify the possibility of fish consumption. Successful results were achieved for 19 individuals out of a total of 35 (Table 4).

In accordance with literature data from different Medieval sites [100102], among the La Bianca community the majority of individuals seemed to consume sheep and cattle products with a minimal amount of fish, pig and chicken (Fig 5).

thumbnail
Fig 5. Percentages of consumption of the different animal dietary resources.

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

The use of the DNA analysis from the dental calculus was useful not only to confirm the stable isotope results but also to understand which of the animal species were predominantly consumed.

Light microscopy analysis

The presence of microfossils in dental calculus was detected in 66% of the individuals. In Table 5, the number of starch granules, pollen grains and other microremains (i.e. phytoliths, calcium oxalate crystals) detected in each sample (lines) was reported, after their taxonomic classification (columns).

thumbnail
Table 5. Amount of starch granules (per taxonomic group or total) and other microremains detected in dental calculus.

The numbers reflect the quantity of each microremain counted by direct microscopy observation in the samples.

https://doi.org/10.1371/journal.pone.0205362.t005

The 277 starch granules detected by LM were clustered in 9 different morphotypes, on the basis of morphometric and morphological parameters (Table 6).

thumbnail
Table 6. Starch granule morphotypes recovered from the dental calculus of the studied individuals.

https://doi.org/10.1371/journal.pone.0205362.t006

Moreover, unidentified remains showing no distinctive features were also recorded. The highest percentage of individuals presented starch of morphotype IX (Fagaceae), followed by miners who consumed morphotype I, II, VI and VII (Poaceae). Particular attention was paid to SU 198 individual, whose dental calculus presented an aggregate containing 220 starch granules, essentially ascribable to morphotype VII (Triticum sp., 192 granules) and to morphotype III (Myristica fragrans Houtt., 24 granules).

Beyond starch granules, microscopic analysis also revealed, the presence of pollen grains of Asteraceae and Laurus nobilis L. in the dental calculus of the individuals. In addition, 8 phytoliths belonging to Poaceae family, 6 calcium oxalates crystals and a Cucurbitaceae fruit epicarp were observed. Examples of images of microremains found in calculus samples are indicated in Fig 6.

thumbnail
Fig 6. Examples of plant microremains found in the dental calculus of the samples.

Starch granule of Triticum sp. and relative polarized image (A); starch granule of Fabaceae and relative polarized image (B); Poaceae phytoliths (C); not determined starch granule (D); fragment of Cucurbitaceae fruit epicarp (E); pollen grain of Asteraceae (F); Poaceae phytoliths (G); calcium oxalate crystal (H); starch granules of Myristica fragrans (I) and relative modern experimental reference (i); starch granule of Avena sp. (J) and modern experimental reference of Avena fatua (j); starch granules of Triticum sp. (K) and modern experimental reference of Triticum durum (k). The black bar indicates 15 μm.

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

GC-MS analysis

GC-MS analysis produced results for 23 individuals. The analysis was repeated for each sample at least 3 times and, in all cases, it always showed similar profiles. In S6 Table the chemical compounds identified in the calculus by GC-MS are listed, clustered by biochemical class. For each chromatographic profile, specific foods (i.e. Artemisia dracunculus L.) or food categories (i.e. milk and derivatives, plant oils), consumed in life at least once by the community of La Bianca were extrapolated, associating the detected compounds as reported in literature and scientific food databases [103104]. Generally, markers of Brassicaceae (i.e. desulfoglucobrassicin; isothiocyanatoacetaldehyde dimethyl acetal; isothiocyanic acid, propyl ester; 13-docosenoic acid, methyl ester) and milk and its derivatives (i.e. beta-lactose) were the most recurrent molecules in the dental calculus of the population, followed by plant oils (i.e. oleic acid; 8,11,14-eicosatrienoic acid, methyl ester, (Z,Z,Z); 9,12-octadecadienoic acid, methyl ester, (E,E)) and theophylline, an alkaloid produced in leaves of Camellia sinensis (L.) Kuntze. The detection of other plant chemical compounds (i.e. lactones, terpens) suggested the key role of plants (i.e. spices, herbs, fruits of Rosaceae) in the miners’ diet. Moreover, evidence of pharmacologically active alkaloids (i.e. stachydrine; securinine; conhydrin) supported the knowledge and application of medicinal plants (i.e. Ephedra sp., Stachys officinalis (L.) Trevis., Artemisia dracunculus, Conium maculatum L. and Securinega suffruticosa (Pall.) Rehd).

Discussion

Archaeological evidence suggests that the archaeological site of La Bianca (Allumiere) is relative in date to one of the first recorded human settlements in the area aimed to alum production. Musculoskeletal stress markers, combined with degenerative disease patterns provided important results in the reconstruction of miners’ working activity. Morphological changes depend on repeated daily exercises that stimulate bone remodeling at the attachment sites, increasing blood flow as a consequence [105106]. However, although the macroscopic analysis of bone morphological modifications can be related to muscles subjected to a high biomechanical stress, often the absence of archaeological data or documents confirming the validity of these assumptions does not allow determining the exact activity carried out in the individual’s lifetime [107].

The present work represents a unique case study, providing evidence of a relationship between occupational stress markers of bone tissue and working activities [78]. This oddity is partly due to the accurate and detailed description of the four phases of alum production that is available from historical written sources [13]. The excavation procedure represented the most physically taxing phase, where picks were used to extract alunite from the rocks. This activity may have also involved climbing with narrow strings on the chest. Calcination was another complex phase, where alunite stones were roasted for 12–14 hours at 600–700°C. At the end of this process, alunite was cooled through a generous watering, in order to break the “roasted” stones and obtain a doughy solution: this phase, defined moistening, could last up to 60 days. The last step was the lixiviation during which the moistened solution was purified by mixing for 24 hours. Even if all four phases (i.e. excavation, calcination, moistening and lixiviation) markedly involved the shoulder girdle, the different scheme of enthesal changes and osteoarthritis pattern observed on the specimens seemed to differ in relation to the biomechanical stress caused by the different tasks; thus allowing the classification of the individuals in four main clusters, each one corresponding to the working phase they were probably most involved in. It is known that osteoarthritis may be a result of traumas or infections, certainly the biomechanical stress suffered from the joint represents one of its causative factors [108] even if the aetiology is far from being fully understood [109]. These data were strongly supported by osteometrics, which reported the existence of different strength and asymmetry patterns among individuals. A tall and robust skeleton could be attributable to a probable selection of individuals to be miners in relation to their adaptability to such type of work. Working tasks and related biomechanical stress seemed to have influenced also the prevalence of skeletal traumas and pathologies. In particular, a comparison with the Latium Medieval populations of Colonna [99] and Santa Severa [110] showed significant differences in the observed pathological pattern. In Colonna and Santa Severa, periostitis was detected in 4.2% and 4.4% of the individuals, respectively, vs the value of 32% found in Allumiere. Similarly, degenerative diseases affected only 16.7% and 22.3% of the skeletal remains, in the two coeval cases, compared to 86% of La Bianca community. The morphological examination demonstrated peculiar and unique features on the analyzed sample showing that it is also one of the few Medieval examples that allows a clear identification of the social status of its components. The combined analysis of enthesal changes and osteoarthritis led to interesting results however the exact relationship between morphological variations on bone tissue and the activity carried out should be considered as a plausible hypothesis. It is known that age at death could also play an important role in the morphology of enthesal changes and that the different studies in literature underline the difficulties in finding a real correlation between entheses and activity patterns [111113].

Stable isotope data supported a diet based on terrestrial proteins. The δ13C values also suggested a contribution of C4 plants to the diet, as confirmed by archaeobotanical microscopy evidence. Indeed, the use of C4 plants as food, such as millet, has been documented in Italy since the Bronze Age [114], and continued in the Medieval period at least in the north east of Italy [95].

DNA analysis revealed the consumption of animal proteins deriving from sheep and cattle meat, whereas pig and chicken use seemed to be negligible. Historical evidence indicates that the consumption of chicken remained low in the Middle Ages, corroborating our results [115]. Similarly, as we observed in the Medieval period, the use of pork meat never reached the elevated consumption of the first century of the Roman Imperial Age [79].

The dietary evidence presented by GC-MS and starch analysis will not be quantifiable in relation to the amounts consumed, however they do give an indication of the types of foods ingested during life, particularly when considered at the population level [20]. The microscopic analysis showed a consumption of Fagaceae, followed by C3 Poaceae caryopses (i.e. Hordeum sp., Triticum sp.). Acorns, being rich in proteins, unsaturated fat, carbohydrates, minerals and vitamins, represented a food with high energy value [116118]. Indeed, they were commonly added to other cereal flours (after elimination of indigestive tannins by thermal treatment in water) or prepared as anti-diarrheal and astringent decoctions [119]. The high number of starch granules that were unable to be taxonomically identified could be due to modifications induced by cooking processes, ptialin enzymimatic activity or grinding procedures. Among the miners, individual SU 198 in particular presented a higher number of starch granules of Triticum sp. and Myristica fragrans (nutmeg). This Asian spice was also detected in the dental calculus of other two individuals, suggesting the use of this plant species in the studied community. Nutmeg trades have been documented in the Mediterranean area since the 6th century, although archaeological findings are rare. It was usually grounded and mixed with red wine or used, in medicine, for its anti-inflammatory properties [120121]. GC-MS analysis revealed that more than half of the individuals consumed Brassicaceae and dairy products, confirming the fundamental role of these foodstuffs in the Medieval diet. At the population level, the consumption of plant oils, contained for example in Olea europaea fruits or acorns [119, 122123], was observed. In addition, the detection of theophylline, an alkaloid with diuretic power, suggested the use of plant species native of South-East Asia, such as Camellia sinensis (the tea plant). To a lesser extent, the use of foods containing cholesterol (perhaps meat, cheese and eggs), Artemisia dracunculus (tarragon, spice with antiseptic, anti-inflammatory and digestive properties) [124], Rosaceae fruits (i.e. lactones) and various herbs were also detected. The identification of secondary metabolites, essentially alkaloids, with pharmaceutical properties suggested the knowledge of medicinal plant species, such as Stachys officinalis (astringent, digestive and sedative) [125] and Ephedra sp. (bronchodilator) [126]. The detection of a metabolite typical of Cucurbitaceae (cucumber aldehyde) supported the use of species such as cucumber and squash; this evidence was also confirmed by the presence, in the calculus of a specimen, of a fragment of cucumber fruit epicarp [127].

Individual SU 176 indicated the consumption of Conium maculatum (hemlock) known to be poisonous but used in low doses in the Middle Ages, along with extract of S. officinalis and fennel seeds, as remedy against epilepsy, muscle spasms, impotence, tubercular adenitis and bites from dogs infected with rabies [128]. Individual SU 185 had evidence for Sorbus domestica, a small tree of the Rosaceae family used in the past for its fruits, called sorbs, rich in vitamin C and individual SU 110 showed evidence for Securinega suffruticosa, an oriental plant containing secondary metabolites that act by stimulating the central nervous system [129]. The peculiar use of the S. suffruticosa, C. sinensis and M. fragrans, typical of Asian areas, could be explained by the existence of commercial trades with Eastern regions or the presence, among the miners, of individuals with Eastern origin, which probably introduced these medicinal species in the community. This latter hypothesis is in line with archaeological and historical data on alum production, a process widely developed in Turkey and more in general in the Eastern Europe [13]. It is not possible to exclude the idea that, in order to establish the first Italian alum extraction system, the Pontifical State would have called on foreign experts. This hypothesis is corroborated by the presence of four individuals with shovel teeth, a typical feature of Eastern origin individuals [130].

This present combined approach allowed us to obtain a more detailed appreciation of the diet. In particular, two individuals SU 317 and SU 245 showed an enrichment of δ15N values >3.5‰ than the maximum animal measurements which indicated the consumption of protein enriched in 15N, perhaps from aquatic sources [131132], but note that trophic level enrichment for humans has been reported up to 6‰ [133]. Analysis of DNA from dental calculus, performed only on SU 245, suggested marine fish consumption. We hypothesized that this enrichment of δ15N value could be also associated to breastfeeding in SU 317 (which was a juvenile individual aged between 2–3 years old) [134] and possibly to nutritional stress in SU 245, who suffered from a compound fracture on the tibia and fibula (Fig 2). It is widely documented that food deficiencies activate gluconeogenesis, a biochemical process, which uses non-carbohydrate sources to produce glucose [96, 135137], and increases δ15N values in tissues [138]. Furthermore, GC-MS analysis of the dental calculus of SU 245 revealed the presence of markers of A. dracunculus, a plant documented to be an anesthetic and anti-inflammatory [124]. It is reasonable to believe that a compound fracture, as that detected in SU 245, caused pain and induced the application of drugs and medical treatments to treat this critical condition.

In conclusion, applying an innovative and original multidisciplinary approach, we present a detailed osteobiography of the first Italian community of alum miners. In detail, according to morphological features, we hypothesized the working task of each individual and reconstructed their dietary patterns. Furthermore, archeobotanical analysis demonstrated that miners used spices and herbs endemic of Asia, for therapeutic purposes, revealing their knowledge about Eastern traditional medicine.

Supporting information

S1 Table. List of the skeletal remains analyzed in the present research housed at the Department of Biology of the University of Rome “Tor Vergata”.

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

(DOCX)

S2 Table. Primer sequences used to amplify animal different DNA regions, size of the PCR product obtained and annealing temperature.

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

(DOCX)

S3 Table. Results of the post-cranial indices. (nr: not recordable).

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

(DOCX)

S4 Table. Presence of pathology/stress markers for each individual.

Neoplasia (NP), Infectious diseases (periostitis, PO; osteomyelitis, OM; infections, IN), trauma and stress markers (fractures, FR; grasping GR), degenerative pathologies (degenerative diseases, DG; and axial degereative diseases (Schmörl’s nodes, SN)) and congenital disorders (CO), and inflammation (IF).

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

(DOCX)

S5 Table. Results of the analysis on degenerative disease pattern (nr: not recordable).

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

(DOCX)

S6 Table. The molecules identified in dental calculus by GC-MS analysis were listed and clustered in biochemical classes for each sample.

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

(DOCX)

Acknowledgments

The authors sincerely thank Gruppo Archeologico Romano (GAR) for their valuable support to the excavation of the analyzed skeletal remains, Marianna D’Amico for the cooperation in the osteological study and Domitilla Tibaldi for her technical collaboration in the molecular paleodiet analyses. The authors would also like to thank Dr. Riccardo Polini for the photographs, Dr. Simone Greco and Dr. Guglielmo Manenti for performing X-rays at Azienda Ospedaliera Universitaria Policlinico Tor Vergata, Matt von Tersch (BioArCh, University of York) for isotope mass spectrometry, and Abigail Sequeira and Sophie Gart for their technical assistance. Gabriele Scorrano is supported by the Marie Skłodowska-Curie Individual Fellowship “PALAEO-ENEO”, a project funded by the European Union’s EU Framework Programme for Research and Innovation Horizon 2020 (Grant Agreement number 751349). The authors thank the editor and the anonymous reviewer for their insightful and constructive comments, which contributed improvements to the present paper.

References

  1. 1. Delumeau J. L’allume di Roma XV-XIX secolo. Civitavecchia: La Litografica; 2003.
  2. 2. Ait I. Dal governo signorile al governo del capitale mercantile: i Monti della Tolfa e 'le lumere'del papa. MEFRM 2014, 126:187–200. Italian.
  3. 3. Santacroce S. Tulpharum alumina lavoro, politica e società sui Monti della Tolfa tra XVI e XVII secolo. Tolfa: Comune di Tolfa; 2014.
  4. 4. Vallelonga F. Ricerche archeologiche nel territorio di Allumiere: gli scavi della Farnesiana e della Bianca. In: Contardi A, editor. Notiziario IX. Allumiere: Petruzzi Stampa. 2012; p. 47–71.
  5. 5. Ambrose SH. Isotopic analysis of paleodiets: methodological and interpretative considerations. In: Sandford K. editor. Investigation of ancient human tissue, chemical analyses in anthropology. Langhorne: Gordon and Breach; 1993. p. 59–130.
  6. 6. 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. p. 305–27.
  7. 7. Schoeninger MJ. Diet reconstruction and ecology using stable isotope ratios. In: Larsen CS, editor. A comparison to biological anthropology. Chichester: Wiley- Blackwell; 2011. p. 445–64.
  8. 8. Henry AG, Brooks AS, Piperno DR. Microfossils in calculus demonstrate consumption of plants and cooked foods in Neanderthal diets (Shanidar III, Iraq; Spy I and II, Belgium). PNAS. 2011; 108: 486–91. pmid:21187393
  9. 9. Buckley S, Usai D, Jakob T, Radini A, Hardy K. Dental calculus reveals unique insights into food items, cooking and plant processing in Prehistoric Central Sudan. PLoS One. 2014; 9(7): e100808. pmid:25028938
  10. 10. Warinner C, Rodrigues JF, Vyas R, Trachsel C, Shved N, Grossmann J, et al. Pathogens and host immunity in the ancient human oral cavity. Nat Genet. 2014; 46: 336–44. pmid:24562188
  11. 11. Lindhe J. Parodontologia ed implantologia dentale. Milano: Edi-Ermes; 1984.
  12. 12. Radini A, Nikita E, Buckley S, Copeland L, Hardy K. Beyond food: The multiple pathways for inclusion of materials into ancient dental calculus. American journal of physical anthropology 2017; 162: 71–83. pmid:28105717
  13. 13. Fumière O, Dubois M, Baeten V, von Holst C, Berben G. Effective PCR detection of animal species in highly processed animal byproducts and compound feeds. Anal Bioanal Chem. 2006; 385: 1045–54. pmid:16761123
  14. 14. 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–50.
  15. 15. Piperno DR, Ranere AJ, Holst I, Iriarte J, Dickau R. Starch grain and phytolith evidence for early ninth millennium BP maize from the Central Balsas River Valley, Mexico. PNAS. 2009; 106: 5019–24. pmid:19307570
  16. 16. Cawthraw S, Saunders GC, Martin TC, Sawyer J, Windl O, Reaney SD. Real-Time PCR detection and identification of prohibited mammalian and avian material in animal feeds. J Food Protect. 2009; 72: 1055e1062.
  17. 17. Yancy HF, Washington JD, Callahan L, Mason JA, Deaver CM, Farrell DE, et al. Development, evaluation, and peer verification of a rapid real-time PCR method for the detection of animal material. J Food Prot. 2009; 72: 2368–74. pmid:19903402
  18. 18. Wesolowski V, de Souza SMFM, Reinhard KJ, Ceccantini G. Evaluating microfossil content of dental calculus from Brazilian sambaquis. J Archaeol Sci. 2010; 37: 1326–38.
  19. 19. Hardy K, Buckley S, Collins MJ, Estalrrich A, Brothwell D, Copeland L, et al. Neanderthal medics? Evidence for food, cooking, and medicinal plants entrapped in dental calculus. Naturwissenschaften. 2012; 99: 617–26. pmid:22806252
  20. 20. 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–57.
  21. 21. Cristiani E, Radini A, Edinborough M, Borić D. Dental calculus reveals Mesolithic foragers in the Balkans consumed domesticated plant foods. PNAS USA 2016; 113: 10298–10303. pmid:27573829
  22. 22. 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. 2016; 1–9.
  23. 23. 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. Quatern Int. 2016; 398: 129–35.
  24. 24. Walker PL, Johnson JR, Lambert PM. Age and Sex Biases in the Preservation of Human Skeletal Remains. Am J Phys Anthropol. 1988; 76: 183–88. pmid:3046371
  25. 25. Todd TW. Age changes in the pubic bone: I. The white male pubis. Am J Phys Anthropol. 1920a; 3: 285–334.
  26. 26. Todd TW. Age changes in the pubic bone. Am J Phys Anthropol. 1920b; 4: 1–70.
  27. 27. Brooks ST, Suchey JM. Skeletal age determination based on the os pubis: a comparison of the Acsádi e Nemeskéri and Suchey-Brooks methods. Hum Evol. 1990; 5: 227–38.
  28. 28. Lovejoy CO, Meindl RS, Pryzbeck TR, Mensforth RP. Chronological metamorphosis of the auricular surface of the ilium: a new method for the determination of adult skeletal age at death. Am J Phys Anthropol. 1985; 68: 47–56.
  29. 29. İşcan MY, Loth SR, Wright RK. Metamorphosis at the sternal rib end: A new method to estimate age at death in white males. Am J Phys Anthropol. 1984; 65: 147–56. pmid:6507605
  30. 30. İşcan MY, Loth SR, Wright RK. Age estimation from the rib by phase analysis: White females. J For Sci.1985; 30: 853–863.
  31. 31. Brothwell DR. Digging up bones: the excavation, treatment, and study of human skeletal remains. Ithaca: Cornell University Press; 1981
  32. 32. Lovejoy CO. Dental wear in the Libben population: its functional pattern and role in the determination of adult skeletal age at death. Am J Phys Anthropol. 1985; 68: 47–56. pmid:4061601
  33. 33. Meindl RS, Lovejoy CO. Ectocranial Suture Closure: A Revised Method for the Determination of Skeletal Age at Death Based on the Lateral-Anterior Sutures. Am J Phys Anthropol. 1985; 68: 57–66. pmid:4061602
  34. 34. Fazecas IG, Kósa F. Forensic Fetal Osteology. Budapest: Akadémiai Kiadó; 1978.
  35. 35. Stloukal M, Hanáková H. Die Länge der Längsknochen altslawischer Bevölkerungen—Unter besonderer Berücksichtigung von Wachstumsfragen. Homo. 1978 29: 53–69. German.
  36. 36. Scheuer L, Black S. Developmental Juvenile Osteology. Cambridge: Cambridge University Press; 2000.
  37. 37. Ubelaker DH. Human Skeletal Remains: Excavation, Analysis, Interpretation 2nd ed. Washington: Taraxacum Washington; 1989.
  38. 38. Minozzi S, Canci A. Archeologia dei resti umani. Dallo scavo al laboratorio. Roma: Carocci Editore; 2015.
  39. 39. Acsádi G, Nemeskéri J. History of Human Life Span and Mortality. Budapest: Akadémiai Kiadó; 1970.
  40. 40. Ferembach D, Schwidetzky I, Stloukal M. Recommandations pour déterminer l'âge et le sexe sur le squelette. Bull Mém Soc Anthropol Paris. 1979; 6: 7–45. French.
  41. 41. Phenice TW. A newly developed visual method of sexing the os pubis. Am J Phys Anthropol. 1969; 30: 297–302. pmid:5772048
  42. 42. Di Bennardo R, Taylor JV. Sex assessment of the femur: a test of a new method. Am J Phys Anthropol. 1979; 50: 635–638. pmid:464035
  43. 43. Bass WM. Human Osteology. A laboratory and field manual, 3rd ed. Columbia: Missouri Archaeological Society; 1987
  44. 44. Berrizbeitia EL. Sex Determination with the Head of the Radius. J For Sci. 1989; 34: 1206–13.
  45. 45. Safont S, Malgosa A, Subirá ME. Sex Assessment on the Basis of Long Bone Circumference. Am J Phys Anthropol. 2000; 113; 317–28. pmid:11042535
  46. 46. Cowal LS, Pastor RF. Dimensional Variation in the Proximal Ulna: Evaluation of a Metric Method for Sex Assessment. Am J Phys Anthropol. 2008; 135: 469–78. pmid:18186505
  47. 47. Martin R, Saller K. Lehrbuch der Antropologie. Stuttgart: G. Fischer; 1957.
  48. 48. Borrini M. Archeologia Forense. Bologna: Editrice Lo Scarabeo; 2007.
  49. 49. Borrini M. Antropologia forense: protocollo e linee guida per il recupero e lo studio dei resti umani. PhD Thesis, Università degli Studi di Roma “Tor Vergata”. 2011.
  50. 50. Borgognini Tarli S, Pacciani E. I resti umani nello scavo archeologico. Metodiche di Recupero e Studio. Roma: Bulzoni Editore; 1993.
  51. 51. Dal Poz M, Ricci F, Reale B, Malvone M, Salvadei L, Manzi G. Paleobiologia della popolazione altomedievale di San Lorenzo di Quingentole, Mantova. Cranio e scheletro postcraniale. In: Manicardi A, editor. San Lorenzo di Quingentole. Archeologia, Storia e Antropologia. Mantova: Sap; 2002. p. 151–198.
  52. 52. Pearson K. Mathematical contributions to the theory of evolution. V. On the reconstruction of stature of prehistoric races. Philos Trans R Soc Lond (Biol). 1899; 192: 169–244.
  53. 53. Telkkä A. On the prediction of human stature from the long bones. Acta Anat (Basel). 1950; 9: 103–17.
  54. 54. De Mendonça MC. Estimation of height from the length of long bones in a Portuguese adult population. Am J Phys Anthropol. 2000; 112: 39–48. pmid:10766942
  55. 55. Rodoinova D, Tenekedjiev K, Yornanov Y. Stature estimation from long bone lengths in Bulgarians. Homo. 2002; 52: 221–32. pmid:12018118
  56. 56. Belmonte Expósito MT. Estimación de la estatura a trevés de la tibia en población española contemporanéa. PhD thesis, Universidad de Granada. 2012. Available from: http://www.aeaof.com/web/blog/tesis-maa-trinidad-belmonte-exposito.pdf
  57. 57. Steele DG. Estimation of stature from fragments of long limb bones. In: Stewart TD, editor. Personal Identification in Mass Disasters. Washington DC: National Museum of Natural History; 1970. p. 85–97.
  58. 58. Mariotti V, Facchini F, Belcastro MG. Enthesopathies-Proposal of a standardized scored method and applications. Coll Antropol. 2004; 28: 145–59. pmid:15636072
  59. 59. Mariotti V, Facchini F, Belcastro MG. The Study of Entheses: Proposal of a Standardised Scoring Method for Twenty-Three Entheses of the Postcranial Skeleton. Coll Antropol. 2007; 31: 291–313. pmid:17598416
  60. 60. Borgognini Tarli S, Reale B. Metodo di analisi degli indicatori non metrici di stress funzionale. Riv Antropol. 1997; 75: 1–39. Italian.
  61. 61. Buikstra JE, Ubelaker DH. Standards for data collection from human skeletal remains. Indianapolis: Western Newspaper Company; 1994.
  62. 62. Prieto I. Protocolo de valoración de patología articular en paleopatología. Estudio de la rodilla en la necropolis hispanomusulmana de San Nicolás(Murcia, s. XI-XIII). Diploma de estudios avanzados, Universidad Autónoma de Madrid, Universidad de Alcalà de Henares. 2009. Available from: https://repositorio.uam.es/handle/10486/3334
  63. 63. Longin R. New method of collagen extraction for radiocarbon dating. Nature 1971; 230: 241. pmid:4926713
  64. 64. Ambrose SH, Norr L. Experimental evidence for the relationship of the carbon isotope ratios of whole diet and dietary protein to those of bone collagen and carbonate. In: Lambert JB, Grupe G, editors. Prehistoric Human Bone: Archaeology at the Molecular Level. New York: Springer-Verlag; 1993. p. 1–37.
  65. 65. DeNiro MJ. Post-mortem preservation and alteration of in vivo bone collagen isotope ratios in ralation to paleodietary reconstruction. Nature. 1985; 317: 806–9.
  66. 66. van Klinken GJ. Bone Collagen Quality Indicators for Palaeodietary and Radiocarbon Measurements. J Archaeol Sci. 1999; 26:687–95.
  67. 67. 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.
  68. 68. 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–61. pmid:28273061
  69. 69. Gilbert MTP, Bandelt HJ, Hofreiter M, Barnes I. Assessing ancient DNA studies. Trends in Ecol Evol. 2005; 20: 541–44.
  70. 70. Willerslev E, Cooper A. Ancient DNA. Proc Biol Sci. 2005; 272: 3–16. pmid:15875564
  71. 71. 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–9. pmid:25231926
  72. 72. 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.
  73. 73. 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
  74. 74. Gismondi A, Rolfo MF, Leonardi D, Rickards O. and Canini A. Identification of ancient Olea europaea L. and Cornus mas L. by DNA barcoding. Comptes Rendus Biologies. 2012; 335: 472–479. pmid:22847014
  75. 75. Hardy K, Blakeney T, Copeland L, Kirkham J, Wrangham R, Collins M. Starch granules, dental calculus and new perspectives on ancient diet. J Archaeol Sci. 2009; 36: 248–55.
  76. 76. Gismondi A, D'Agostino A, Canuti L, Di Marco G, Basoli F, Canini A. Starch granules: a data collection of 40 food species. Plant Biosystems. 2018;
  77. 77. Gismondi A, D’Agostino A, Canuti L, Di Marco G, Martínez -Labarga C, Angle M, Rickards O, Canini A. Dental calculus reveals diet habits and medicinal plant use in the Early Medieval Italian population of Colonna. J Archaeol Sci Reports 2018; 20: 556–564.
  78. 78. Baldoni M, Borrini M, Martínez-Labarga C. Musculoskeletal Stress Markers as Evidence of Labor: the Allumiere Cemetery Case. Proceedings of 69th Meeting American Academy of Forensic Sciences; 2017 Feb 13–18; New Orleans, LA. p. 71.
  79. 79. Craig OE, Biazzo M, O'Connell TC, Garnsey P, Martinez-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–83. pmid:19280672
  80. 80. Mori E, Bozzi R, Laurenzi A. Feeding habits of the crested porcupine Hystrix cristata L. 1758 (Mammalia, Rodentia) in a Mediterranean area of Central Italy. The European Zoological Journal. 2017; 84: 261–265.
  81. 81. Minagawa M, Matsui A, Ishiguro N. Patterns of prehistoric boar Sus scrofa domestication, and inter- islands pig trading across the East China Sea, as determined by carbon and nitrogen isotope analysis. Chem Geol. 2005, 218: 91–102.
  82. 82. 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.
  83. 83. Fuller BT, De Cupere B, Marinova E, Van Neer W, Waelkens M, Richards MP. Isotopic reconstruction of human diet and animal husbandry practices during the Classical-Hellenistic, imperial, and Byzantine periods at Sagalassos, Turkey. Am J Phys Anthropol. 2012, 149: 157–171. pmid:22729657
  84. 84. Hamilton J, Thomas R. Pannage, pulses and pigs: isotopic and zooarchaeological evidence for changing pig management practices in later Medieval England. Mediev Archaeol. 2012; 56: 234–259.
  85. 85. Madgwick R, Mulville J, Stevens RE. Diversity in foddering strategy and herd management in late Bronze Age Britain: an isotopic investigation of pigs and other fauna from two midden sites. Environ Archaeol. 2012, 17: 126–140.
  86. 86. Hammond C, O’Connor T. Pig diet in medieval York: carbon and nitrogen stable isotopes. Archaeol Anthrop Sci. 2013; 5: 123–127.
  87. 87. 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.
  88. 88. Navarrete V, Colonese AC, Tornero C, Antolín CF, Von Tersch M, Eulàlia Subirà M, Comes P, Rosell‐Melé A, Saña M. Feeding Management Strategies among the Early Neolithic Pigs in the NE of the Iberian Peninsula. Int J Osteoarchaeol. 2017, 27: 839–852.
  89. 89. Montanari M. Alimentazione e cultura nel medioevo. Bari: Laterza; 1998.
  90. 90. 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.
  91. 91. O'Connell TC, Kneale CJ, Tasevska N, Kuhnle GG. The diet‐body offset in human nitrogen isotopic values: A controlled dietary study. Am J Phys Anthropol. 2012; 149: 426–434 pmid:23042579
  92. 92. Reitsema L, Vercellotti G. Stable isotope evidence for sex- and status-based variations in diet and life history at medieval Trino Vercellese, Italy. Am J Phys Anthropol. 2012; 148: 589–600. pmid:22553011
  93. 93. Ricci P, Mongelli V, Vitiello A, Campana S, Sirignano C, Rubino M, Fornaciari G, Lubritto C. The priviledged burial of the Pava Pieve (Siena, 8th century AD). Rapid Commun Mass Spectrom., 2012; 26: 2393–2398. pmid:22976205
  94. 94. Pescucci L, Battistini A, De Angelis F, Catalano P. Vivere al centro di Roma nell'VIII secolo d.C. indicazioni antropologiche. Boll Archeol On-Line. 2013; IV: 2–4.
  95. 95. 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–564. pmid:24889200
  96. 96. 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–56. pmid:24706415
  97. 97. Ciaffi R, Lelli R, Müldner G, Stantcheva K, Fischetti AL, Ghini G, Craig O, Milano F, Rickards O, Arcudi G, Martínez-Labarga C. Palaeobiology of the medieval population of Albano (Rome, Italy): a combined morphological and biomolecular approach. Int J Osteoarchaeol. 2015; 25: 477–488.
  98. 98. 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. Herit Sci. 2015; 3:27.
  99. 99. Baldoni M, Nardi A, Müldner G, Lelli R, Gnes M, Ferraresi F, et al. Archaeobiological recontruction of the Italian medieval population of Colonna (8th-10th centuries CE). J Archaeol Sci Rep. 2016; 10: 483–94.
  100. 100. Wilkens B. I resti faunistici, in Farnese, Testimonianze archeologiche di vita quotidiana dai butti del centro storico, a cura del gruppo Archeologico Mediovaldarno. Firenze: All’Insegna del Giglio; 1991.
  101. 101. Colonnelli G, De Grossi Mazzorin J. 2000. Nuovi dati sull’alimentazione a Farnese (VT) nei secoli XV e XVI. Atti del 2° Convegno Nazionale di Archeozoologia; 1997 Nov 14–16; Asti, Italy. Forlì: Abaco; 1997. p. 369–376.
  102. 102. Minniti C. Economia e alimentazione nel Lazio medievale: nuovi dati dalle evidenze archeozoologiche. Archeol Med 2009; XXXVI: 273–83. Italian.
  103. 103. FooDB, 2013. FoodDB version 1.0. http://fooddb.ca/.
  104. 104. TGSC, 2015. The Good Scents Company. http://www.thegoodscentscompany.com/
  105. 105. Hawkey D, Merbs C. Activity-induced muscoloskeletal stress markers (MSM) and subsistence strategy changes among ancient Hudson Bay Eskimos. Int J Osteoarchaeol. 1995; 5: 324–38.
  106. 106. Ruff CB, Holt B, Trinkaus E. Who's afraid of the big bad Wolff?: “Wolff's law” and bone functional adaptation. Am J Phys Anthropol. 2006; 129: 484–98. pmid:16425178
  107. 107. Dutour O. Enthesopathies (Lesions of Muscular Insertions) as indicators of the Activities of Neolithic Saharan Populations. Am J Phys Anthropol.1986; 71: 221–24. pmid:3541646
  108. 108. Weiss E, Jurmain R. Osteoarthritis revised: A contemporary review of Aetiology. Int J Osteoarchaeol. 2007; 17: 437–450.
  109. 109. Ortner DJ. Identification of pathological conditions in human skeletal remains. San Diego: Academic Press; 2003.
  110. 110. Gnes M, Baldoni M, Marchetti L, Basoli F, Leonardi D, Canini A, et al. The medieval population of Santa Severa (7th-15th centuries CE): A morphological approach to investigate its bioarchaeology. J Archaeol Sci Rep. 2018; 18: 11–25.
  111. 111. Milella M, Belcastro MG, Zollikofer CP, Mariotti V. The effect of age, sex, and physical activity on entheseal morphology in a contemporary Italian skeletal collection. Am J Phys Anthropol. 2012; 148: 379–388. pmid:22460619
  112. 112. Milella M, Cardoso FA, Assis S, Lopreno GP, Speith N. Exploring the relationship between entheseal changes and physical activity: A multivariate study. Am J Phys Anthropol. 2015; 156: 215–223. pmid:25330328
  113. 113. Henderson CY, Nikita E. Accounting for multiple effects and the problem of small sample sizes in osteology: a case study focussing on entheseal changes. Archaeol Anthropol Sci. 2016; 8: 805–817.
  114. 114. Tafuri MA, Craig OE, Canci A. Stable isotope evidence for the consumption of millet and other plants in Bronze Age Italy. Am J Phys Anthropol. 2009;, 139: 146–153. pmid:19051259
  115. 115. De Grossi Mazzorin J. 2015. I resti archeozoologici come strumento di conoscenza ell’economia alimentare nell’alto medioevo. L’alimentazione nell’Alto Medioevo: pratiche, simboli, ideologie. Atti delle LXIII Settimane di studio della fondazione centro italiano di studi sull’alto medioevo 2015 Apr 9–14; Spoleto, Italy. Spoleto: CISAM; 2015 p. 182–91
  116. 116. Pignone D, Laghetti G. On sweet acorn (Quercus spp.) cake tradition in Italian cultural and ethnic islands. Genet Resour Crop Evol 2010; 57: 1261–66.
  117. 117. Šálková T, Divišová M, Kadochová Š, Beneš J, Delawská K, Kadlcˇková E, et al. Acorns as a food resource. An experiment with acorn preparation and taste. IANSA. 2011; II: 9–14.
  118. 118. 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.
  119. 119. 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.
  120. 120. Asgarpanah J, Kazemivash N. Phytochemistry and pharmacologic properties of Myristica fragrans Hoyutt.: A review. African Journal of Biotechnology 2012; 11:12787–12793.
  121. 121. Preusza M, Kodýdkováa K, Kočárb P, Vaněčekc Z. Exotic spices in flux: archaeobotanical material from Medieval and Early Modern sites of the Czech Lands (Czech Republic). IANSA. 2015; 6: 223–236.
  122. 122. Di Fiore B, Chianese D. Electric and magnetic tomographic approach for the geophysical investigation of an unexplored area in the archaeological site of Pompeii (Southern Italy). J Archaeol Sci. 2008; 35, 14–25.
  123. 123. Kanthilatha N, Boyd W, Dowell A, Mann A, Chang N, Wohlmuth H, et al. Identification of preserved fatty acids in archaeological floor sediments from prehistoric sites at Ban Non Wat and Nong Hua Raet in northeast Thailand using gas chromatography. J Archaeol Sci. 2014; 46: 353–362.
  124. 124. Mohammad Reza S, Hamideh M, Zahra S. The nociceptive and anti-inflammatory effects of Artemisia dracunculus L. aqueous extract on fructose fed male rats. Evid Based Complement Alternat Med. 2015;
  125. 125. Popović Z, Smiljanić M, Kostić M, Nikić P, Janković S. Wild flora and its usage in traditional phytotherapy (Deliblato Sands, Serbia, South East Europe). Indian J Tradit Know. 2014; 13: 9–35.
  126. 126. Limberger RP, Jacques ALB, Schmitt GC, Arbo MD. Pharmacological effects of ephedrine. Nat Prod 2012; 1217–37.
  127. 127. Piperno DR, Holst I, Wessel-Beaver L, Andres TC. Evidence for the control of phytolith formation in Cucurbita fruits by the hard rind (Hr) genetic locus: archaeological and ecological implications. PNAS. 2002; 99: 10923–8. pmid:12149443
  128. 128. Vetter J. Poison hemlock (Conium maculatum L.). Food Chem Toxicol. 2004; 42: 1373–82. pmid:15234067
  129. 129. Raj D, Łuczkiewicz M. Securinega suffruticosa. Fitoterapia. 2008; 79: 419–27. pmid:18534774
  130. 130. Hillson SW. Dental Anthropology. London: Cambridge University Press; 1996
  131. 131. Hedges REM, Reynard LM. Nitrogen isotopes and the trophic level of humans in archaeology. J Archaeol Sci. 2007; 34: 1240–51.
  132. 132. Muldner G, Richards MP. Stable isotope evidence for 1500 years of human diet at the city of York, UK. Am J Phys Anthropol. 2007; 133:682–97. pmid:17295296
  133. 133. O'Connell TC, Kneale CJ, Tasevska N, Kuhnle GG,. The diet‐body offset in human nitrogen isotopic values: A controlled dietary study. Am J Phys Anthropol. 2012; 149:426–434 pmid:23042579
  134. 134. 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. Geophysical Laboratory. Washington, DC, pp 111–17; 1989.
  135. 135. Katzenberg MA, Lovell NC. Stable Isotope Variation in Pathological Bone. Int J Osteoarchaeol. 1999; 9: 316–24.
  136. 136. 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–10. pmid:16628564
  137. 137. Olsen KC, White CD, Longstaffe FJ, von Heyking K, McGlynn G, Grupe G, et al. Intraskeletal isotopic compositions (d13C, d15N) of bone collagen: nonpathological and pathological variation. Am J Phys Anthropol. 2014; 153:598–604. pmid:24374993
  138. 138. 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