# Microsatellite genotyping of medieval cattle from central Italy suggests an old origin of Chianina and Romagnola cattle

Maria Gargani, Lorraine Pariset, Johannes A. Lenstra, Elisabetta de Menicis, Katayoun Moazami-Goudarzi, Denis Laloë, Alessio Valentini

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

► To cite this version:

Maria Gargani, Lorraine Pariset, Johannes A. Lenstra, Elisabetta de Menicis, Katayoun Moazami-Goudarzi, et al.. Microsatellite genotyping of medieval cattle from central Italy suggests an old origin of Chianina and Romagnola cattle. Frontiers in Genetics, 2015, 6 (68), pp.1-6. hal-01606647

HAL Id: hal-01606647

https://hal.science/hal-01606647

Submitted on 27 May 2020

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

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

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

# Microsatellite genotyping of medieval cattle from central Italy suggests an old origin of Chianina and Romagnola cattle

Maria Gargani<fnref n="1" />*, Lorraine Pariset<fnref n="1" />, Johannes A. Lenstra<fnref n="2" />, Elisabetta De Minicis<fnref n="3" />, European Cattle Genetic Diversity Consortium† and Alessio Valentini<fnref n="1" />

**Edited by:**

Stéphane Joost, École
Polytechnique Fédérale de
Lausanne, Switzerland

**Reviewed by:**

Pablo Orozco ter Wengel, Cardiff
University, UK
Letizia Nicoloso, Università degli
Studi di Milano, Italy
Florian Alberto, CNRS, France

**\*Correspondence:**

Maria Gargani, Department for
Innovation in Biological, Agro-food
and Forest systems, University of
Tuscia, Via S. Camillo de Lellis,
01100 Viterbo, Italy
e-mail: maria.gargani@unitus.it

Analysis of DNA from archeological remains is a valuable tool to interpret the history of ancient animal populations. So far most studies of ancient DNA target mitochondrial DNA (mtDNA), which reveals maternal lineages, but only partially the relationships of current breeds and ancient populations. In this study we explore the feasibility of nuclear DNA analysis. DNA was extracted from 1000-years old cattle bone collected from Ferento, an archeological site in central Italy. Amplification of 15 microsatellite FAO-recommended markers with PCR products yielded genotypes for four markers. Expected heterozygosity was comparable with values of modern breeds, but observed heterozygosity was underestimated due to allelic loss. Genetic distances suggested a position intermediate between (1) Anatolian, Balkan, Sicilian and South-Italian cattle and (2) the Iberian, North-European and Central-European cattle, but also a clear relationship with two central-Italian breeds, Chianina and Romagnola. This suggests that these breeds are derived from medieval cattle living in the same area. Our results illustrate the potential of ancient DNA for reconstructing the history of local cattle husbandry.

**Keywords:** microsatellite, ancient DNA, cattle, Ferento, NeighborNet

## INTRODUCTION

The study of ancient DNA (aDNA) has developed since about 28 years from the analysis of short segments of mitochondrial DNA (mtDNA) to the spectacular whole-genome analysis of close relatives of the *Homo sapiens* (Higuchi et al., 1984; Pääbo, 1985; Reich et al., 2010). It complements archeological and historical investigations by providing ancestral links between ancient samples and present organisms (Vernesi et al., 2004; Caramelli et al., 2007). The majority of aDNA studies still target mtDNA, which is facilitated by its high copy number per cell. This allowed to

determine the sequences of mtDNA control regions extracted from cattle remains excavated at several sites (Lenstra et al., 2014), but also the entire mitochondrial genome of two *Bos primige-nius* specimens (Edwards et al., 2010; Lari et al., 2011). However, the genealogy inferred from the mtDNA not necessarily reflect that inferred with nuclear markers. Analysis of nuclear DNA from ancient remains is more challenging, but has been reported in sev-eral studies (Greenwood et al., 1999, 2001; Noonan et al., 2005; Pariset et al., 2007; Green et al., 2010; Reich et al., 2010). Ancient genotypes of nuclear microsatellites markers are informative for genetic diversity, subdivision and geographical origin (Edwards et al., 2003; Allentoft et al., 2011; Ishida et al., 2012; Nyström et al., 2012). Moreover, there are several studies carried out with the so called FAO list of microsatellites in almost every continent on a huge amount of individuals of several breeds (Ajmone-Marsan and The GLOBALDIV Consortium, 2010). This valuable database could be related to ancient remains in order i. to dis-cover the closeness of ancient individuals to extant populations to understand the type of husbandry and production of the time and ii. to infer the origin, migration and admixture of present breeds.

Authentic cattle breeds in central and southern Italy are of the Podolian type (Felius, 1995). However, it is not yet clear if these are of ancient origin or have been imported (Felius et al., 2014). The aim of this study was to gain insight into his-tory of Italian cattle by comparing 1000-year old cattle DNA

<footnote n="1">¹ Department for Innovation in Biological, Agro-food and Forest systems, University of Tuscia, Viterbo, Italy</footnote>

<footnote n="2">² Faculty of Veterinary Medicine, Utrecht University, Utrecht, Netherlands</footnote>

<footnote n="3">³ Department of Sciences of Cultural Heritage (DISBEC), University of Tuscia, Viterbo, Italy</footnote>

<footnote>†The following members of the European Cattle Genetic Diversity Consortium contributed to this study: K. Moazami-Goudarzi, D. Laloë, INRA, Jouy-en-Josas, France; J. L. Williams (presently at the Parco Tecnologica Padano, Lodi, Italy), P. Wiener, Roslin Institute, UK; D. G. Bradley, Trinity College, Dublin, Ireland; G. Erhardt, Justus-Liebig Universität, Giessen, Germany; B. Harlizius, School of Veterinary Medicine, Hannover, Germany; I. Medugorac, Ludwig-Maximilians-Universität, Munich C. Looft, E. Kalm, Christian-Albrechts-Universität, Kiel, Germany; J. Cañón, S. Dunner, Universidad Complutense de Madrid, Spain; C. Rodellar, I. Martín-Burriel, Veterinary Faculty, Zaragoza, Spain; A. Sanchez, Universitat Autonoma de Barcelona; A. M. Martínez, Universidad de Córdoba, Córdoba, Spain; C. Ginja, Instituto Superior de Agronomía, Lisboa, Portugal; P. Ajmone-Marsan, Università Cattolica del S. Cuore, Piacenza, Italy; F. Pilla, A. Bruzzone, Università del Molise, Campobasso, Italy; D. Marletta, Università degli Studi di Catania, Catania, Italy; L. E. Holm, Danish Institute of Agricultural Sciences, Tjele, Denmark; G. Dolf, University of Berne, Switzerland.</footnote>
of Italian origin with DNA from modern cattle. DNA was extracted from bones collected in Ferento, an archeological site near Viterbo in central Italy inhabited since Bronze Age, but developed mainly during Roman and Medieval ages. It rises on a triangular plateau made of tuff stone [Pianicara, IGM F.137 II NE] and delimited by two ditches: Guzzarella (or Vezzarella) and Acqua Rossa, subaffluents of the Tiber. The archeological investigations in Ferento revealed five areas, named assay I, assay II, assay III, assay IV, and assay V. The assay I comprises the area northward the Decumano in the medieval quarter; the assay II is localized on west area at the limit of the plateau, where an important necropolis was found; the assay III was occupied with a residential building on the southern side of the Decumano and with an important system of tanks, leant against the Theater; the assay IV, placed in the northwestern area of the plateau, and the assay V, near the assay I close to the fortification.

We compared the microsatellite genotypes obtained for four microsatellites with those of modern cattle breeds in order to uncover the origins of the Central-Italian medieval cattle.

# MATERIALS AND METHODS

## ANCIENT AND MODERN SAMPLES

Thirty bone samples of medieval cattle (~1000 BP) were collected at Ferento, an archeological site near Viterbo in central Italy (Table 1). The analyzed remains were recovered from five different excavation areas. The dating was made on the basis of pottery found in the same layer. Moreover, one samples were sent to CEDAD (CEntre for DAting and Diagnostics) belonging to the University of Salento (Department of Innovation Engineering) and were carbon-dated using high resolution mass spectrometry. The results confirmed the age estimated by pottery layers. Species identification, for the selection of the samples to be analyzed, was based on morphology and dimensions of the specimens.

In order to avoid contamination and degradation by sudden change in environmental conditions, samples were picked up using latex gloves, put immediately in vacuum-pack, stored at -20°C (Pariset et al., 2007) and transferred to a lab dedicated to ancient DNA analysis. In a dedicated laboratory room, the outer surface was removed with sand paper and bone powder was collected by drilling. DNA was extracted according to Yang et al. (1998) and quantified using a DTX Multimode Detector

<table><caption>Table 1 | Medieval samples.</caption><thead><tr><th>Sample</th><th>Element</th><th>Age</th><th>Excavation areas</th><th>GPS</th><th>Amplification</th></tr></thead><tbody><tr><td>Ferento 1</td><td>Astragalus</td><td>~1000 BP</td><td>5</td><td>42°29/21.2"N, 12°07/55.8"E</td><td>No</td></tr><tr><td>Ferento 2</td><td>Phalanx II</td><td>~1000 BP</td><td>3</td><td>42°29/18.4"N, 12°07/53.9"E</td><td>Yes</td></tr><tr><td>Ferento 3</td><td>Phalanx I</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>No</td></tr><tr><td>Ferento 4</td><td>Heel</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>Yes</td></tr><tr><td>Ferento 5</td><td>Molar 3</td><td>~1000 BP</td><td>3</td><td>42°29/18.4"N, 12°07/53.9"E</td><td>Yes</td></tr><tr><td>Ferento 6</td><td>Phalanx I</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>Yes</td></tr><tr><td>Ferento 7</td><td>Astragalus</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>Yes</td></tr><tr><td>Ferento 8</td><td>Astragalus sx</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>Yes</td></tr><tr><td>Ferento 9</td><td>Jawbone</td><td>~1000 BP</td><td>3</td><td>42°29/18.4"N, 12°07/53.9"E</td><td>No</td></tr><tr><td>Ferento 10</td><td>Phalanx II</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>Yes</td></tr><tr><td>Ferento 11</td><td>Jawbone</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>Yes</td></tr><tr><td>Ferento 12</td><td>Heel</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>No</td></tr><tr><td>Ferento 13</td><td>Jawbone</td><td>~1000 BP</td><td>2</td><td>42°29/19.4"N, 12°07/47.2"E</td><td>Yes</td></tr><tr><td>Ferento 14</td><td>Metatarsus</td><td>~1000 BP</td><td>3</td><td>42°29/18.4"N, 12°07/53.9"E</td><td>No</td></tr><tr><td>Ferento 15</td><td>Heel</td><td>~1000 BP</td><td>5</td><td>42°29/21.2"N, 12°07/55.8"E</td><td>Yes</td></tr><tr><td>Ferento 16</td><td>Phalanx I</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>No</td></tr><tr><td>Ferento 17</td><td>Astragalus</td><td>~1000 BP</td><td>1</td><td>42°29/20.6"N, 12°07/56.7"E</td><td>No</td></tr><tr><td>Ferento 18</td><td>Phalanx I</td><td>~1000 BP</td><td>1</td><td>42°29/20.6"N, 12°07/56.7"E</td><td>Yes</td></tr><tr><td>Ferento 19</td><td>Phalanx I</td><td>~1000 BP</td><td>1</td><td>42°29/20.6"N, 12°07/56.7"E</td><td>No</td></tr><tr><td>Ferento 20</td><td>Metatarsus</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>Yes</td></tr><tr><td>Ferento 21</td><td>Phalanx I</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>Yes</td></tr><tr><td>Ferento 22</td><td>Metatarsus</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>Yes</td></tr><tr><td>Ferento 23</td><td>Phalanx II</td><td>~1000 BP</td><td>3</td><td>42°29/18.4"N, 12°07/53.9"E</td><td>Yes</td></tr><tr><td>Ferento 24</td><td>Heel</td><td>~1000 BP</td><td>5</td><td>42°29/21.2"N, 12°07/55.8"E</td><td>Yes</td></tr><tr><td>Ferento 25</td><td>Metatarsus</td><td>~1000 BP</td><td>2</td><td>42°29/19.4"N, 12°07/47.2"E</td><td>No</td></tr><tr><td>Ferento 26</td><td>Phalanx I</td><td>~1000 BP</td><td>3</td><td>42°29/18.4"N, 12°07/53.9"E</td><td>No</td></tr><tr><td>Ferento 27</td><td>Humerus</td><td>~1000 BP</td><td>2</td><td>42°29/19.4"N, 12°07/47.2"E</td><td>No</td></tr><tr><td>Ferento 28</td><td>Phalanx II</td><td>~1000 BP</td><td>1</td><td>42°29/20.6"N, 12°07/56.7"E</td><td>No</td></tr><tr><td>Ferento 29</td><td>Phalanx I</td><td>~1000 BP</td><td>2</td><td>42°29/19.4"N, 12°07/47.2"E</td><td>No</td></tr><tr><td>Ferento 30</td><td>Metatarsus</td><td>~1000 BP</td><td>4</td><td>42°29/19.2"N, 12°07/52.4"E</td><td>No</td></tr></tbody></table>
880 (Beckman) with Picogreen method (Quant-iT PicoGreen, Invitrogen). At least two independent DNA extractions were performed. One single sample per day was extracted with a negative control. Comprehensive datasets of modern samples of genotypes of cattle breeds from Europe, Africa and Asia (Table S1) have been described previously (European Cattle Genetic Diversity Consortium, 2006; Medugorac et al., 2009; Laloë et al., 2010; Felius et al., 2011; Delgado et al., 2012).

## MICROSATELLITE GENOTYPING

Fifteen microsatellites were selected from the FAO list of recommended microsatellites (Table S2) (FAO, 2011). Only for three of these markers the published primers generated PCR products shorter than 160 bp. In order to optimize the amplification with aDNA templates, we redesigned primers for the other 12 microsatellites (Table S2) using Primer3 software. In a separate laboratory room, a first-round of PCR was carried out using a 20 µl reaction volume containing PCR buffer with 2.5 mM MgCl₂, 200 µM of each dNTPs, 0.05 µM of both the forward (fluorescent labeled dyes) and the reverse primer (Proligo), 0.2 units of Bio-x-act short Taq Polymerase (Bioline) and 10 ng of DNA. A 5 min denaturation step at 94°C was followed by 30 cycles of 30 s denaturation at 94°C, 1 min at the respective annealing temperatures (Table S2) and 1 min of extension at 68°C, followed by a final 5 min extension step at 68°C. Products were purified with exosap (Promega). A second round of PCR was performed with as template the product (1 µl) of the first PCR reaction. During each PCR experiment both extraction and the PCR negative controls were checked for the absence of amplification products. At least four amplifications per sample were performed in order to monitor for allelic drop out. PCR products were separated by electrophoresis on a CEQ 8800 sequencer (Beckman Coulter) and sized with the standard CEQ™ DNA Size Standard Kit—400 (Beckman Coulter). The alleles were scored using the proprietary CEQ fragment analysis software.

Microsatellite dataset from modern cattle breeds were obtained from the EU Resgen project. In order to compare the allele lengths of the ancient cattle with those of the modern cattle database, we have included in all plates subject to PCR three modern samples (ITROM10, ITROM11, ITROM20) as a reference. The ancient cattle alleles were recalculated for comparison with the reference individuals. Genotypes were scored as heterozygous when the electropherogram showed a clear biallelic profile or when two different mono-allelic profiles were observed in different PCRs (Allentoft et al., 2011). The allelic profiles have been deposited in the Dryad database (http://dx.doi.org/10.5061/dryad.d4500). Allele frequencies and observed and expected heterozygosity were calculated for each population using MSToolkit. Genetic distances were calculated by using the program MSAT2 (http://genetics.stanford.edu/hpgl/projects/microsat/) and visualized in NeighborNet graphs with SplitsTree4 program (http://www.splittree.org/). The NeighborNet graphs were simplified by combining several breeds in regional pools as approximations of ancestral gene populations. Thus, we made metapopulations of the dairy breeds from the Northwest-European Lowland, Central Brown cattle (including Swiss-Brown, Rendena, Bruna, and Cabannina), a West-Central cluster (Simmental and

its derivatives); Iberian cattle, Sicilian Modicana and Cinesara, Balkan Busha, Anatolian cattle and Indopakistan zebu breeds, respectively. Principal component analysis of a matrix of FST genetic distances was performed using the program Genalex (Peakall and Smouse, 2006). Unsupervised model-based clustering was done with the program Structure (Pritchard et al., 2000).

## RESULTS

Extraction of genomic DNA yielded for all 30 ancient samples quantifiable DNA and for 16 of these also amplification products. From the 15 microsatellites, four (INRA5, HEL1, MM12, CSRM60) could be amplified successfully in at least one sample. A total of 760 amplifications were performed on ancient samples. The success rate PCR per sample and per locus is shown in Table S3. Markers INRA5, HEL1, MM12, and CSRM60 yielded genotypes in 6, 8, 10, and 10 samples, respectively. All extraction and PCR controls for different DNA isolation, purification and amplification steps did not yield PCR products. In 35 PCR reactions we observed allelic dropout of 50% of heterozygous genotypes for HEL1, 60% for CSRM60, 75% for MM12, 73% for INRA5 and 64% averaged of the four loci, but never more than two alleles per marker in a given sample. The expected (He) and observed (Ho) heterozygosity values for a panel of European breeds calculated on the basis of 4 microsatellites ranged from 0.48 (Lagunaire breed) to 0.771 (Istrian breed) and from 0.43 (Lagunaire breed) to 0.76 (Istrian breed) respectively. The He and Ho values for pools of breeds ranged from 0.68 (Zebu breed pool) to 0.786 (Sicilian breed pool) and from 0.62 (Zebu breed pool) to 0.79 (Anatolian breed pool) respectively. These values were in the same range of those based on 30 microsatellites: He and Ho for panel of European breeds ranged from 0.544 (Lagunaire breed) to 0.741 (Podolica breed) and from 0.51 (Lagunaire breed) to 0.73 (Bohemian Red breed); He and Ho for pools of breed ranged from 0.67 (Zebu breed pool) and 0.780 (Anatolian breed pool) and from 0.59 (Zebu breed pool) and 0.72 (Busha breed pool) respectively (Table S1). The expected heterozygosity (0.698) in Ferento samples is within the range of the values found for modern breeds. We observed a clear heterozygote deficit (observed heterozygosity 0.588) related to allelic dropout for ancient samples.

Principal component analysis (Figure S1) as well as model-based clustering of Ferento cattle (not shown) combined with a panel of European breeds typed with the same markers differentiated three clusters: (1) African cattle, (2) Northwestern-European Lowland dairy breeds, and (3) Central-European and Iberian cattle as well as the ancient samples from Ferento. However, PCA did not establish relationships between Ferento cattle and individual breeds, since relative positions of breeds could not be reproduced with other panels of four microsatellites (not shown). We are aware that four microsatellites may have a low power in resolving diversity in populations that diverged since very few generations, however, we could not increase the number of markers that successfully and repeatedly amplified ancient DNA within the FAO list that permits to relate ancient to extant breeds.

In order to get a more informative comparison of ancient Ferento genotypes with those of other cattle, we constructed NeighborNet graphs. **Figure 1** shows NeighborNet graphs of the clusters with several modern Podolian breeds and ancient
4 microsatellites

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

30 microsatellites

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

FIGURE 1 | NeighborNet graphs of Reynold's distances of Ferento cattle with modern Podolian breeds, and European, Anatolian and Indian breed clusters based on 4 or 30 microsatellites. Ferento cattle is only shown with 4 microsatellites.

Ferento cattle. Although distances on the basis of only four microsatellites are not expected to be accurate, these distances reproduce the clustering that are generated by the 30 microsatel- lites. Interestingly, Ferento cattle is clearly linked with the Tuscanian Chianina and the related Romagnola cattle. The ancient and modern central Italian cattle is intermediate between (1) the Anatolian, Sicilian and Balkan Podolian breeds, and (2) the Iberian, Central European and Lowland dairy breeds. We observed that none of the individuals from modern breeds were as close to Ferento as the Chianina and Romagnola breeds.

DISCUSSION

In this study we demonstrate the feasibility of genotyping microsatellites for the study of cattle of 1000 years ago. The rates of amplification we have obtained (between 11 and 19%) appear to be lower than those found in in previous studies (Edwards et al., 2003; Nyström et al., 2012), which is proba- bly due to relatively high temperature in central Italy. Although this limits a more elaborate analysis of the ancient genotypes, our results reconstruct the history of Italian cattle. The origin of the Podolian cattle in Italy is controversial (Felius et al., 2014). The first cattle arrived in the Italian peninsula at the start of the Neolithic (ca. 8000 BP). These long-horned cattle were largely replaced by short-horned cattle around 4500–5000 BP. Neolithic cattle gradually decreased in size and, as in the rest of Europe, Bronze and Iron age cattle in Italy had wither heights of typi- cally 110 cm. However, cattle with wither heights of 115–135 cm

were imported from Epirus in the centuries before the Roman
era and are probably the ancestors of the large Roman cat-
tle with wither heights of 150 cm or more (Kron, 2002). The
analysis of mitochondrial DNA from modern taurine cattle has
detected five mtDNA haplogroup with a specific geographic dis-
tribution (Lenstra et al., 2014). The relatively high frequencies of
the T and T2 haplogroups of Italian Podolian cattle (Bonfiglio
et al., 2010; Lenstra et al., 2014) may be explained by these
imports of cattle. In an alternative scenario Anatolian cattle were
brought to Tuscany by the founders of the Etruscan civilization
(Pellecchia et al., 2007). After the Roman era, the large cattle
disappeared from the European fossil record and cattle became
even smaller than in the Iron age with a typical wither heights of
95 cm for the medieval Ferento cattle. Since the Renaissance cattle
became larger and the Tuscan Chianina is now the largest cat-
tle of the world. This breed and two other central-Italian breeds,
Romagnola and Marchigiana, are related to the South Italian and
Balkan long-hornd Podolian breeds (Felius, 1995). These cat-
tle are believed to have emerged after the 12th century by local
selection for long horns (Bökönyi, 1974). There have been sev-
eral opportunities for gene flow between the Balkans and Italy,
e.g., the consecutive migrations of Wisigoths, Osthrogoths and
Lombards and the large scale import of Hungarian Gray cattle
as “meat on the hoof” after the renaissance via Venice. The last
imports have been well documented and involved not only oxes,
but also fertile bulls (Appuhn, 2010). The deviating mtDNA hap-
logroup distribution of the Italian Podolian cattle indicates that
any introgression of Balkan cattle was via bulls only (Lenstra et al., 2014). Alternatively, it has also been proposed that Italian Podolian cattle descends from ancient local cattle or even local aurochs (Ciani and Matassino, 2001).

The present data contributes to a historic reconstruction by indicating a link of the modern large, white and short-horned Chianina and Romagnola with local cattle living in the same area 1000 years earlier. This argues against an origin of the Chianina, Romagnola and the close relative Marchigiana from the postmedieval imports. Since the Germanic invasions are not likely to have replaced completely the local cattle, it is plausible that Ferento cattle descends from cattle kept during the Roman era. The large white central-Italian and Podolian cattle do not form a tight cluster in the network, which leaves open the possibility of an eastern origin of the long-horned Italian Podolian cattle. This may include the long-horned and feral Maremmana, for which no genotypes are available for all four microsatellites, but is in other networks close to the Hungarian Gray.

In conclusion our study demonstrates the feasibility and the potentiality of the analysis of nuclear DNA from ancient livestock remains.

## AUTHOR CONTRIBUTIONS

MG designed the study, carried out microsatellite amplifications experiments, and drafted the manuscript, LP contributed to the design of the work, JL contributed substantially to the analysis and interpretation of data and revised the work, ED provided ancient samples bones, ECGDC provided modern dataset, AV participated in developing ideas and revised the manuscript.

## ACKNOWLEDGMENTS

This work has been supported by the European Commission (project ResGEN 98-118). The content of this publication does not represent the views of the Commission or its services.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fgen.2015.00068/abstract

## REFERENCES

* Ajmone-Marsan, P., and The GLOBALDIV Consortium. (2010). A global view of livestock biodiversity and conservation—GLOBALDIV. *Anim. Genet.* 41, 1–5. doi: 10.1111/j.1365-2052.2010.02036.x

* Allenftoß, M. E., Oskam, C., Houston, J., Hale, M. L., Gilbert, M. T. P., Rasmussen, M., et al. (2011). Profiling the dead: generating microsatellite data from fossil bones of extinct megafauna—protocols, problems, and prospects. *PLoS ONE* 6:e16670. doi: 10.1371/journal.pone.0016670

* Appuhn, K. (2010). Ecologies of beef: eighteenth-century epizootics and the environmental history of early modern Europe. *Environ. Hist.* 15, 268–287. doi: 10.1093/envhis/emq049

* Bökönyi, S. (1974). *History of Domestic Mammals in Central and Eastern Europe*. Budapest: Akadémiai Kiadó.

* Bonfiglio, S., Achilli, A., Olivieri, A., Negrini, R., Colli, L., Liotta, L., et al. (2010). The enigmatic origin of bovine mtDNA haplogroup R: sporadic interbreeding or an independent event of Bos primigenius domestication in italy? *PLoS ONE* 5:e15760. doi: 10.1371/journal.pone.0015760

* Caramelli, D., Vernesi, C., Sanna, S., Sampietro, L., Lari, M., Castri, L., et al. (2007). Genetic variation in prehistoric Sardinia. *Hum. Genet.* 122, 327–336. doi: 10.1007/s00439-007-0403-6

* Ciani, F., and Matassino, D. (2001). Il bovino grigio allevato in Italia; origine. Nota 1: il bovino Macrocero. *Taurus Speciale* 13, 89–99.

* Delgado, J. V., Martínez, A. M., Acosta, A., Alvarez, L. A., Armstrong, E., Camacho, E., et al. (2012). Genetic characterization of Latin-American Creole cattle using microsatellite markers. *Anim. Genet.* 43, 2–10. doi: 10.1111/j.1365-2052.2011.02207.x

* Edwards, C. J., Connellan, J., Wallace, P. F., Park, S. D. E., McCormick, F. M., Olsaker, I., et al. (2003). Feasibility and utility of microsatellite markers in archaeological cattle remains from a Viking Age settlement in Dublin. *Anim. Genet.* 34, 410–416. doi: 10.1046/j.0268-9146.2003.01043.x

* Edwards, C. J., Magee, D. A., Park, S. D. E., McGettigan, P. A., Lohan, A. J., Murphy, A., et al. (2010). A complete mitochondrial genome sequence from a mesolithic wild aurochs (Bos primigenius). *PLoS ONE* 5:9255. doi: 10.1371/journal.pone.0009255

* European Cattle Genetic Diversity Consortium. (2006). Marker-assisted conservation of European cattle breeds: an evaluation. *Anim. Genet.* 37, 475–481. doi: 10.1111/j.1365-2052.2006.01511.x

* FAO. (2011). *Molecular Genetic Characterization of Animal Genetic Resources*. Rome: FAO Animal Production and Health Guidelines, No. 9.

* Felius, M. (1995). *Cattle Breeds-An Encyclopedia*. Doetinchem: Misset.

* Felius, M., Beerling, M.-L., Buchanan, D. S., Theunissen, B., Koolmees, P. A., and Lenstra, J. A. (2014). On the history of cattle genetic resources. *Diversity* 6, 705–750. doi: 10.3390/d6040705

* Felius, M., Koolmees, P. A., Theunissen, B., European Cattle Genetic Diversity Consortium, and Lenstra, J. A. (2011). On the breeds of cattle—historic and current classifications. *Diversity* 3, 660–692. doi: 10.3390/d3040660

* Green, R. E., Krause, J., Briggs, A. W., Maricic, T., Stenzel, U., Kircher, M., et al. (2010). A draft sequence of the Neandertal genome. *Science* 328, 710–722. doi: 10.1126/science.1188021

* Greenwood, A. D., Capelli, C., Possnert, G., and Pääbo, S. (1999). Nuclear DNA sequences from late Pleistocene megafauna. *Mol. Biol. Evol.* 16, 1466–1473. doi: 10.1093/oxfordjournals.molbev.a026058

* Greenwood, A. D., Lee, F., Capelli, C., DeSalle, R., and Tikhonov, A. (2001). Evolution of endogenous retro-virus-like elements of the woolly mammoth (Mammuthus primigenius) and its relatives. *Mol. Biol. Evol.* 18, 840–847. doi: 10.1093/oxfordjournals.molbev.a003865

* Higuchi, R., Bowman, B., Freiberger, M., Ryder, O. A., and Wilson, A. C. (1984). DNA sequences from the quagga, an extinct member of the horse family. *Nature* 312, 282–284. doi: 10.1038/312282a0

* Ishida, Y., Roca, A. L., Fratpietro, S., and Greenwood, A. D. (2012). Successful Genotyping of Microsatellites in the Woolly Mammoth. *J. Hered.* 103, 459–464. doi: 10.1093/jhered/esr139

* Kron, G. (2002). Archaeozoological evidence for the productivity of Roman livestock farming. *MBAH* 21, 53–73.

* Laloë, D., Moazami-Goudarzi, K., Lenstra, J. A., Ajmone-Marsan, P., Azor, P., Rodero, E., et al. (2010). Spatial trends of domestic ruminants in Europa. *Diversity* 2, 932–945. doi: 10.3390/d2060932

* Lari, M., Rizzi, E., Mona, S., Corti, G., Catalano, G., Chen, K., et al. (2011). The Complete Mitochondrial Genome of an 11,450-year-old Aurochsen (Bos primigenius) from Central Italy. *BMC Evol. Biol.* 11:32. doi: 10.1186/1471-2148-11-32

* Lenstra, J. A., Ajmone-Marsan, P., Beja-Pereira, A., Bollongino, R., Bradley, D. G., Colli, L., et al. (2014). Meta-analysis of mitochondrial DNA reveals several population bottlenecks during worldwide migrations of cattle. *Diversity* 6, 178–187. doi: 10.3390/d6010178

* Medugorac, I., Medugorac, A., Russ, I, Veit-Kensch, C. E., Taberlet, P., Luntz, B., et al. (2009). Genetic diversity of European cattle breeds highlights the conservation value of traditional unselected breeds with high effective population size. *Mol. Ecol.* 18, 3394–3410. doi: 10.1111/j.1365-294X.2009.04286.x

* Noonan, J. P., Hofreiter, M., Smith, D., Priest, J. R., and Rohland, N. (2005). Genomic sequencing of Pleistocene cave bears. *Science* 309, 597–599. doi: 10.1126/science.1113485

* Nyström, V., Humphrey, J., Skoglund, P., Mc Keown, N. J., Vartanyan, S., Shaw, P. W., et al. (2012). Microsatellite genotyping reveals end-Pleistocene decline in mammoth autosomal genetic variation. *Mol. Ecol.* 21, 3391–3402. doi: 10.1111/j.1365-294X.2012.05525.x

* Pääbo, S. (1985). Molecular cloning of ancient Egyptian mummy DNA. *Nature* 314, 644–645. doi: 10.1038/314644a0

* Pariset, L., Gabbianelli, F., De Bernardis, D., Fortunato, M. T., Alhaique, F., De Minicis, E., et al. (2007). Ancient DNA: genomic amplification of roman and medieval bovine bones. *Ital. J. Anim. Sci.* 6, 179–181. doi: 10.4081/ijas.2007.1s.179
Peakall, R., and Smouse, P. E. (2006). GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. *Mol. Ecol. Notes* 6, 288–295. doi: 10.1111/j.1471-8286.2005.01155.x

Pellecchia, M., Negrini, R., Colli, L., Patrini, M., Milanesi, E., Achilli, A., et al. (2007). The mystery of Etruscan origins: novel clues from Bos taurus mitochondrial DNA. *Proc. Biol. Sci.* 274, 1175–1179. doi: 10.1098/rspb. 2006.0258

Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. *Genetics* 155, 945–959.

Reich, D., Green, R. E., Kircher, M., Krause, J., Patterson, N., Durand, E. Y., et al. (2010). Genetic history of an archaic hominin group from Denisova Cave in Siberia. *Nature* 468, 1053–1060. doi: 10.1038/ nature09710

Vernesi, C., Caramelli, D., Dupanloup, I., Bertorelle, G., Lari, M., Cappellini, E., et al. (2004). The Etruscans: a population-genetic study. *Am. J. Hum. Genet.* 74, 694–704. doi: 10.1086/383284

Yang, D. Y., Eng, B., Waye, J. S., Dudar, J. C., and Saunders, S. R. (1998). Technical note: improved DNA extraction from ancient bones using silica-based spin columns. *Am. J. Phys. Anthropol.* 105, 539–543.

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Received: 25 November 2014; accepted: 10 February 2015; published online: 04 March 2015.

Citation: Gargani M, Pariset L, Lenstra JA, De Minicis E, European Cattle Genetic Diversity Consortium and Valentini A (2015) Microsatellite genotyping of medieval cattle from central Italy suggests an old origin of Chianina and Romagnola cattle. *Front. Genet.* 6:68. doi: 10.3389/fgene.2015.00068

This article was submitted to Livestock Genomics, a section of the journal Frontiers in Genetics.

Copyright © 2015 Gargani, Pariset, Lenstra, De Minicis, European Cattle Genetic Diversity Consortium and Valentini. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.