Mitochondrial DNA variants of Podolian cattle breeds testify for a dual maternal origin

Abstract

Background

Over the past 15 years, 300 out of 6000 breeds of all farm animal species identified by the Food and Agriculture Organization of the United Nations (FAO) have gone extinct. Among cattle, many Podolian breeds are seriously endangered in various European areas. Podolian cattle include a group of very ancient European breeds, phenotypically close to the aurochs ancestors (Bos primigenius). The aim of the present study was to assess the genetic diversity of Podolian breeds and to reconstruct their origin.

Methodology

The mitochondrial DNA (mtDNA) control-regions of 18 Podolian breeds have been phylogenetically assessed. Nine non-Podolian breeds have been also included for comparison.

Conclusion

The overall analysis clearly highlights some peculiarities in the mtDNA gene pool of some Podolian breeds. In particular, a principal component analysis point to a genetic proximity between five breeds (Chianina, Marchigiana, Maremmana, Podolica Italiana and Romagnola) reared in Central Italy and the Turkish Grey. We here propose the suggestive hypothesis of a dual ancestral contribution to the present gene pool of Podolian breeds, one deriving from Eastern European cattle; the other arising from the arrival of Middle Eastern cattle into Central Italy through a different route, perhaps by sea, ferried by Etruscan boats. The historical migration of Podolian cattle from North Eastern Europe towards Italy has not cancelled the mtDNA footprints of this previous ancient migration.

About the Authors:

Piera Di Lorenzo

Contributed equally to this work with: Piera Di Lorenzo, Hovirag Lancioni

Roles Formal analysis, Writing – original draft

Affiliation: Dipartimento di Scienze Agrarie, Alimentari e Ambientali, Università degli Studi di Perugia, Perugia, Italy

Hovirag Lancioni

Contributed equally to this work with: Piera Di Lorenzo, Hovirag Lancioni

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

* E-mail: [email protected] (HL); [email protected] (AA).

Affiliation: Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Perugia, Italy

Simone Ceccobelli

Roles Formal analysis, Validation

Affiliation: Dipartimento di Scienze Agrarie, Alimentari e Ambientali, Università degli Studi di Perugia, Perugia, Italy

Licia Colli

Roles Formal analysis, Methodology, Supervision

Affiliations Institute of Zootechnics, Università Cattolica del S. Cuore, Piacenza, Italy, Biodiversity and Ancient DNA Research Center–BioDNA, Università Cattolica del S. Cuore, Piacenza, Italy

Irene Cardinali

Roles Formal analysis

Affiliation: Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Perugia, Italy

Taki Karsli

Roles Data curation, Formal analysis

Affiliation: Department of Animal Science, Faculty of Agriculture, University of Akdeniz, Antalya, Turkey

Marco Rosario Capodiferro

Roles Data curation, Formal analysis

Affiliation: Dipartimento di Biologia e Biotecnologie “L. Spallanzani”, Università di Pavia, Pavia, Italy

Emine Sahin

Roles Data curation

Affiliation: Korkuteli Vocational School, University of Akdeniz, Antalya, Turkey

Luca Ferretti

Roles Data curation, Supervision, Writing – review & editing

Affiliation: Dipartimento di Biologia e Biotecnologie “L. Spallanzani”, Università di Pavia, Pavia, Italy

Paolo Ajmone Marsan

Roles Data curation, Resources, Writing – review & editing

Affiliations Institute of Zootechnics, Università Cattolica del S. Cuore, Piacenza, Italy, Biodiversity and Ancient DNA Research Center–BioDNA, Università Cattolica del S. Cuore, Piacenza, Italy

Francesca Maria Sarti

Roles Formal analysis, Supervision

Affiliation: Dipartimento di Scienze Agrarie, Alimentari e Ambientali, Università degli Studi di Perugia, Perugia, Italy

Emiliano Lasagna

Roles Conceptualization, Data curation, Writing – review & editing

Affiliation: Dipartimento di Scienze Agrarie, Alimentari e Ambientali, Università degli Studi di Perugia, Perugia, Italy

Francesco Panella

Roles Conceptualization, Supervision, Writing – review & editing

Affiliation: Dipartimento di Scienze Agrarie, Alimentari e Ambientali, Università degli Studi di Perugia, Perugia, Italy

Alessandro Achilli

Roles Conceptualization, Funding acquisition, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

* E-mail: [email protected] (HL); [email protected] (AA).

Affiliation: Dipartimento di Biologia e Biotecnologie “L. Spallanzani”, Università di Pavia, Pavia, Italy

Introduction

Over the past 15 years, 300 out of 6000 livestock breeds identified by Food and Agriculture Organization of the United Nations (FAO) have gone extinct [1]. Risk factors for farm animal breeds are mainly: i) a reduction of genetic variability due to strict selection processes; ii) a strong economic pressure focused on specific traits, such as milk production, which leads to the replacement of local less productive breeds with highly productive industrial breeds; iii) an unrestricted and indiscriminate cross-breeding, especially in developing countries [2]. Bos taurus is one of the most economically important livestock species [3]. Both in historic and current societies it has fulfilled agricultural, economic, cultural, and even religious key roles, often paralleling human evolution [4]. Among cattle, many Podolian breeds are seriously endangered in various European countries [5]; [6]; [7]. Podolian cattle include a group of very ancient European breeds, with a grey coat colour and long horns, phenotypically close to the aurochs (Bos primigenius). According to many traditional notes the name Podolian refers to a common ancestral origin in Podolia (the modern western Ukraine). However place of origin and timing of spread out of the source area are both debated. Alternative hypotheses have been proposed: Podolian cattle might have spread from the eastern steppe southward into Anatolia and westward into the Balkans and Italy in historical times (3rd-5th century AD) along with East-European Barbarian people [8]; other authors suggest a more ancient migration (~3 kya BP) from the Near East to Central Italy through the Mediterranean Sea [9], together with a possible contribution from local wild aurochs through a secondary local domestication/introgression events[10]; [11].

Nowadays, some phenotypic distinctions stand out among Podolian cattle [12]. The noble aurochs-shaped ancient breeds with long horns (such as Hungarian Grey, Katerini, Podolsko, Slavonian Syrmian and Maremmana) are considered as the only true Podolian breeds by some scholars. However, some local breeds (i.e. Podolica Italiana, Ukrainian Grey, Turkish Grey and other Balkan breeds) do not necessarily show the long horns, but maintained some distinctive Podolian traits such as a red coat in calves and light grey in adults [13]. In general, a commercial trait shared by all Podolian cattle is that they are more suitable for beef production rather than for dairy. Because of that some improved beef breeds (Chianina, Marchigiana, Romagnola and Piemontese) are also considered within the Podolian group, although the inclusion of Chianina and Piemontese is still debated [14].

During the last decades, mitochondrial DNA (mtDNA) has been widely used as a molecular tool to investigate genetic origin, history and diversity of livestock species [15] [16] [17] [18] [19] [20]. Following this trend, the aim of the present study is to re-assess the mtDNA diversity of the major Podolian cattle breeds (ten of which classified as endangered or critical by FAO) to obtain additional information on their ancestral origin and ancient dispersal routes.

Results

We have analysed the mitochondrial DNA of 18 Podolian cattle breeds (Fig 1, Table 1, S1–S3 Tables).

[Figure omitted. See PDF.]

Fig 1. Prevalent locations and frequency distributions of mitochondrial haplogroups in 18 Podolian (P, underlined) and 9 non-Podolian (non-P) breeds analyzed in this study.

Breed codes as in Table 1 (see also Table 2 and S3 Table for details). Note that RP and IB are also widespread in Italy. Maps (www.histgeo.ac-aix-marseille.fr/ancien_site/carto/) were used for illustrative purposes only and largely modified by the authors.

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

[Figure omitted. See PDF.]

Table 1. Estimates of genetic diversitya on the 27 breeds analyzed in this work.

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

The molecular analysis of 221 base pairs of the control region (from np 16042 to np 16262) on the entire dataset of 1,957 samples revealed a total of 247 distinct haplotypes (from four to 70 haplotypes per breed) and 91 polymorphic sites (S), all represented by single nucleotide polymorphisms (SNPs). The average nucleotide diversity (π) was comparable between Podolian and non-Podolian breeds (~0.010–0.011; Table 1), while haplotype diversity was significantly lower (P-value < 0.01) in Podolian (Hd = 0.837 ± 0.010) than in non-Podolian breeds (Hd = 0.879 ± 0.013). Both indices varied largely across breeds as already seen in previous mtDNA studies [11]; [21]; [22]; [23]. Among all Podolian, the highest Hd values (≥0.90) were identified in Chianina, Marchigiana and Turkish Grey, while the lowest values of Hd (<0.70), as well as of nucleotide diversity (≤0.008), were scored in Bianca di Val Padana, Calvana, and Slavonian Syrmian Podolian. As for 1,617 sequences, we were able to extend the analysis to a longer control-region fragment encompassing 731 bps (Table 1). The results largely confirmed the same trend, with the only notable exception of Piemontese, Romagnola and Maremmana showing a higher Hd (>0.960) on this extended fragment. It is also interesting that the highest Hd values were identified in Chianina, and Maremmana, which showed values (>0.970) comparable to the Turkish Grey (0.971).

All control-region haplotypes have been classified in haplogroups and sub-haplogroups through an accurate analysis of mutational motifs (Table 2 and S3 Table), according to previously published classification criteria [11]; [24]; [25]; [26]. Haplogroup T3 was the most common (83%) in all breeds, with the highest value in MP and PO (both 100%), followed by PI (96%) (acronyms are listed in Table 1). The second and third most common haplogroups (both 7%) were T1 and T2, which were missing in MP, PO and SS. The frequency of T2 is lower in non-Podolian (3.85%) than in Podolian breeds (8.35%) with extraordinary high peaks in three breeds, KA (42%), RG (24%) and BG (22%), followed by UK (16%) and MA and TK (both 15%). T1 haplogroup was predominantly found among breeds from Central and Southern Italy, both in Podolian (10%) and non-Podolian (19%) groups (Fig 1). Haplogroup T5 was found exclusively in non-Podolian breeds, and was restricted to IB, RP and VA except for one sequence found in SS and one in PI. Finally, haplogroups Q and R showed very low incidences restricted to Italian non-Podolian (Q = 1.15%, and R = 0.58%) and Podolian (Q = 0.77%, and R = 0.49%) breeds. Overall, the haplogroup distribution differed significantly between the Podolian and non-Podolian groups of breeds included in the current analysis (Table 2; chi-square P-value < 0.001) with the highest contribution given by the T2 haplogroup. This result was also verified by considering haplogroup frequencies based on different haplotypes (chi-square P-value < 0.001) in order to mitigate the effect of inbreeding and recent founder effects.

[Figure omitted. See PDF.]

Table 2. Sources and haplogroup affiliation for the Podolian and non-Podolian mtDNA sequences.

Haplogroup frequencies (%) are in parentheses.

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

Thus, we performed a principal component analysis (PCA) to graphically display the different haplogroup distributions among the Podolian breeds. In order to consider as many populations as possible, the dataset based on the short fragment was included. After variables reduction to PCs, the coordinates of the observations for the 18 breeds were displayed in a two-dimensional plot representing the European Podolian genetic landscape (Fig 2). PC1 clearly separated the Turkish Grey and the five most important Central and Southern Italian beef cattle breeds (CH, RO, MR, MA and IP; S1 Fig) from all the remaining populations, while PC2 contributed to separate the Hungarian and the Turkish Grey.

[Figure omitted. See PDF.]

Fig 2. Principal Component Analysis (PCA) of all Podolian breeds.

Below is the plot of the contribution of each haplogroup to the first and second PC (projections of the axes of the original variables).

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

Because of the peculiar position of some Italian breeds, we used an analysis of molecular variance (AMOVA) to investigate fixation indices in three (artificially created) population groups, one including the Italian Podolian breeds, the other encompassing the European Podolian breeds and the Turkish Grey, and the third group covering the Italian non-Podolian breeds. Most of the variance (about 98%) observed in the Italian Podolian populations explained differences among samples within breeds, while less than 2% represented differences between breeds (Table 3), a value three times lower than in the other Podolian breeds.

[Figure omitted. See PDF.]

Table 3. Analysis of molecular variance (AMOVA).

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

Discussion

To date, only a limited number of studies have addressed the genetic composition of Podolian breeds. These investigations were generally limited to few breeds and focused on the nuclear genome [8]; [27]; [28]; [29]; [30]. Mitochondrial DNA data have been previously reported by Ivankovic et al. [6] and Ilie et al. [7] who analysed the control-region sequences of some Croatian and Romanian cattle breeds, respectively. The present study extended the analysis of the mitochondrial genetic variation to 18 Podolian breeds by evaluating their haplogroup distributions, which were eventually compared among them and to nine non-Podolian breeds. Genetic diversity considered in terms of number of haplotypes and nucleotide diversity revealed some peculiarities of several Podolian breeds. The low mtDNA diversities of Bianca di Val Padana, Calvana and Slavonian Syrmian Podolian could be due to a combination of factors, such as i) a sampling bias depending on the low consistence of current herds (S1 Table), ii) a bottleneck effect caused by the strong reduction in population size experienced by these breeds during the last decades, iii) genetic drift acting on small populations. On the contrary, the large diffusion of Piemontese, Marchigiana and Chianina cattle probably favoured the accumulation and maintenance of a high level of mtDNA variation, which is evident also in the Maremmana breed in spite of its lower consistency (S1 Table). In general, the most common European haplogroup T3 is predominant (83%) in our dataset. High frequencies of T1 in central and southern Italy might be due to the intensive migrations across the Mediterranean Sea, eased by the proximity to northern Africa, where T1 is prevalent [20]; [25]; [26]; [31] and to the Near East where T1 is also present [32, 33]. It is worth noting that the central-Italian Podolian breeds show higher frequencies of T2, while the presence of Q and R within the Podolian group is limited to three Italian breeds: Romagnola, Chianina and Marchigiana, the latter derived from crossbreeding between the first two in the early 20th century.

The significant differences between the haplogroup distributions of Podolian and non-Podolian and the low genetic differentiation among the Italian Podolian breeds (three times less than in other Podolian breeds or in the non-Podolian group) points to a common (and perhaps peculiar) origin. As a matter of fact, according to the first component of the PCA, five Italian beef cattle breeds formed a clearly separated group and were also closer to the Turkish Grey than to any other Podolian breeds. At first this peculiarity could be explained as the effect of a stronger beef-oriented selection carried out on these breeds compared to the other Podolian populations. However, another important Italian beef cattle, the Piemontese, is placed within the other Podolian cluster, which shares the common feature of a strong grey coat. Thus, an alternative explanation might assume a different ancestral origin for the two groups of Podolian breeds, as summarized in Fig 1: a first group, mostly consisting of breeds from East Europe and northern Italy that share a similar mitochondrial gene pool, may have originated from ancestors migrated through an inland route from Podolia across eastern Europe all the way into northern Italy, in accordance with the great wave of cattle migration occurred during the Barbarian invasions. Whereas, a second group, including the white Podolian cattle closely related to the Turkish Grey, may descend from ancestral bovines brought to Italy through a different and likely maritime route crossing the Mediterranean Sea. A previous study [34] suggested also a genetic link between the Turkish Grey and Bulgarian and Hungarian breeds, but our results do not support such hypothesis, highlight instead a stronger maternal relationship between the Turkish Grey and five central-southern Italian Podolian breeds. Those cattle are bred since the medieval time [29] in an area that largely overlaps with the ancient territory of Etruria. This finding further supports and extends another hypothesis, according to which at least part of the maternal ancestry of those breeds could be related to the Etruscan migration from Lydia, a region on the south-western coast of ancient Anatolia [9]; [35]; [36], [37]; [38]. It is worth noting that the five Podolian breeds are also the main Italian beef cattle together with the Piemontese, and that previous studies suggested a possible common genetic origin [39]. Our findings suggest that, in spite of a stronger beef-oriented selection, their mitochondrial gene pool still preserves genetic traces of a different maternal origin, confirming that the selection practices were mostly male-mediated and enforcing the importance of the mtDNA screening to reconstruct the ancestry and history of current breeds.

Material and methods

Ethics statement

All experimental procedures were reviewed and approved by the Animal Research Ethics Committee of the Universities of Perugia and Pavia in accordance with the European Union Directive 86/609.

Samples

The entire dataset analyzed in this study encompasses 1,957 mtDNA control-region sequences including 1,321 from our previous studies [11]; [24]; [40], 428 retrieved from GenBank, and 208 additional samples specifically collected for this study (Table 1, S1–S3 Tables). Piemontese (also called Piedmontese) (PI, n = 72), Romagnola (RO, n = 225), Marchigiana (MR, n = 146), Chianina (CH, n = 369), Maremmana (MA, n = 75), Podolica Italiana (also known as Italian Podolian) (IP, n = 125), Mucco Pisano (MP, n = 33), Calvana (CA, n = 35), Bianca di Val Padana (BP, n = 45), Hungarian Grey (HG, n = 93), Bulgarian Grey (BG, n = 36), Istrian cattle (IC, n = 17), Katerini (KA, n = 12), Romanian Grey (RG, n = 17), Slavonian Syrmian Podolian (SS, n = 9), Turkish Grey (TK, n = 85), Ukrainian Grey (UK, n = 32), Podolsko (PO, n = 11). Moreover, nine unrelated non-Podolian breeds from Italy, are included as a control group: Valdostana (VA, n = 54), Bruna Italiana (also known as Italian Brown) (IB, n = 34), Grigio Alpina (also called Grey Alpine) (GA, n = 45), Pezzata Rossa Italiana (also known as Italian Red Pied) (RP, n = 136), Modicana (MO, n = 41), Reggiana (RE, n = 38), Agerolese (AG, n = 36), Cinisara (CI, n = 81), Cabannina (CB, n = 55).

DNA extraction, amplification and sequencing

As for the 208 novel mtDNAs, blood samples were collected from the jugular vein of each animal in vacutainer tubes, containing EDTA as anticoagulant. These animals were chosen in different farms in order to avoid closely related individuals and gather a representative sample of the breeds. Whole blood was stored at -20°C until DNA extraction. DNA was isolated using the GenElute Blood Genomic DNA kit (Sigma Aldrich, St. Louis, MO, USA) and stored at -20°C until genotyping. PCR amplification of the control region was performed using forward and reverse primers (5’-CCTAAGACTCAAGGAAGAAACTGC-3’ and 3’-AACCTAGAGGGCATTCTCACTG-5’ respectively) specifically designed on the Bovine Reference Sequence (BRS; GenBank V00654). The 1138 bp PCR fragment encompassed the mtDNA control region from np 15718 to 517. Amplicons were first purified using exonuclease I and alkaline phosphatase (ExoSAP-IT® enzymatic system-USB Corporation, Cleveland, OH, USA), then sequenced with the primer 15757F (5’-CCCCAAAGCTGAAGTTCTAT-3’), as previously described [40]. A dataset of 1,321 sequences was already available in our laboratories. All data were recorded in GenBank with accession numbers MF474376-MF475904 (S3 Table) and compared to those retrieved from the database (S2 Table).

Data analyses

Sequences were aligned to the Bovine Reference Sequence (BRS; V00654) using the software SequencherTM 5.0. For a total of 1,617 samples, we were able to analyse a 731-bps fragment trimmed from np 15823 to np 215, while only a short fragment (221 bps, from np 16042 to np 16262) was considered in order to include the widest possible number of samples (N = 1,957). Haplotypes were classified in haplogroups and sub-haplogroups according to previously identified mutational motifs [24]. Indices of molecular variation were calculated using the DNAsp 5.1 software [41], while an analysis of molecular variance was computed using AMOVA program implemented in the ARLEQUIN 3.01 package [42]. Finally, principal component analyses (PCA) were performed using Excel software implemented by XLSTAT, as described elsewhere [43, 44]. The PCA is a widely used dimension-reduction method that summarizes the variance of multivariate data in a smaller number of variables (the principal components, PCs), which are linear functions of the original variables, here expressed as haplogroup and sub-haplogroup frequencies, The rarest haplogroups were phylogenetically grouped and frequencies were calculated by considering only different haplotypes within the same breed.

Supporting information

[Figure omitted. See PDF.]

S1 Table. List of Podolian breeds analyzed in this study.

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

(XLSX)

S2 Table. List of samples retrieved from GenBank.

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

(XLSX)

S3 Table. List of samples analysed in this study.

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

(XLSX)

S1 Fig. The five most important Italian beef cattle breeds from central and southern Italy discussed in this paper.

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

(PDF)

Citation: Di Lorenzo P, Lancioni H, Ceccobelli S, Colli L, Cardinali I, Karsli T, et al. (2018) Mitochondrial DNA variants of Podolian cattle breeds testify for a dual maternal origin. PLoS ONE 13(2): e0192567. https://doi.org/10.1371/journal.pone.0192567

References

1. Scherf B. World Watch List for Domestic Animal Diversity. 3rd ed. Rome. 2000.

2. Soysal M, Tuna Y, Gurcan E, Ozkan E. Farms in Turkey: sustainable development in the preservation of animal genetic resources in Turkey and in the world. Trakia J Sci. 2004;2(3):47–53.

3. Cunningham E. Selected Issues in Livestock Industry Development. Washington, DC.1992.

4. Bradley DG, MacHugh DE, Cunningham P, Loftus RT. Mitochondrial diversity and the origins of African and European cattle. Proc Natl Acad Sci U S A. 1996;93(10):5131–5. pmid:8643540; PubMed Central PMCID: PMCPMC39419.

5. Pariset L, Joost S, Marsan PA, Valentini A, Econogene Consortium E. Landscape genomics and biased FST approaches reveal single nucleotide polymorphisms under selection in goat breeds of North-East Mediterranean. BMC Genet. 2009;10:7. Epub 2009/02/19. pmid:19228375; PubMed Central PMCID: PMCPMC2663570.

6. Ivankovic A, Paprika S, Ramljak J, Dovc P, M. K. Mitochondrial DNA-based genetic evaluation of autochthonous cattle breeds in Croatia. Czech J Anim Sci 2014;59(11):519–28.

7. Ilie DE, Cean A, Cziszter LT, Gavojdian D, Ivan A, Kusza S. Microsatellite and Mitochondrial DNA Study of Native Eastern European Cattle Populations: The Case of the Romanian Grey. PLoS One. 2015;10(9):e0138736. Epub 2015/09/23. pmid:26398563; PubMed Central PMCID: PMCPMC4580412.

8. Maretto F, Ramljak J, Sbarra F, Penasa M, Mantovani R, Ivankovic A, et al. Genetic relationship among Italian and Croatian Podolian cattle breeds assessed by microsatellite markers. Livest Sci. 2012;Sept 13;150 (1–3):256–64.

9. Pellecchia M, Negrini R, Colli L, Patrini M, Milanesi E, Achilli A, et al. The mystery of Etruscan origins: novel clues from Bos taurus mitochondrial DNA. Proc Biol Sci. 2007;274(1614):1175–9. pmid:17301019; PubMed Central PMCID: PMCPMC2189563.

10. Beja-Pereira A, Caramelli D, Lalueza-Fox C, Vernesi C, Ferrand N, Casoli A, et al. The origin of European cattle: evidence from modern and ancient DNA. Proc Natl Acad Sci U S A. 2006;103(21):8113–8. pmid:16690747; PubMed Central PMCID: PMCPMC1472438.

11. Bonfiglio S, Achilli A, Olivieri A, Negrini R, Colli L, Liotta L, et al. The enigmatic origin of bovine mtDNA haplogroup R: sporadic interbreeding or an independent event of Bos primigenius domestication in Italy? PLoS One. 2010;5(12):e15760. Epub 2011/01/07. pmid:21209945.

12. Bodò I, Gera I, Koppàny G. The Hungarian Grey cattle breed. Budapest: Association of the hungarian grey cattle breeders; 2004. 128 p p.

13. Bonadonna T. Etnologia zootecnica. Torino: UTET.

14. Bartosiewicz L. Hungarian Grey cattle in search of origins. Hungarian Agricultural research 1996;(3):13–20.

15. Groeneveld LF, Lenstra JA, Eding H, Toro MA, Scherf B, Pilling D, et al. Genetic diversity in farm animals—a review. Anim Genet. 2010;41 Suppl 1:6–31. pmid:20500753.

16. Achilli A, Olivieri A, Soares P, Lancioni H, Hooshiar Kashani B, Perego UA, et al. Mitochondrial genomes from modern horses reveal the major haplogroups that underwent domestication. Proc Natl Acad Sci U S A. 2012;109(7):2449–54. pmid:22308342; PubMed Central PMCID: PMCPMC3289334.

17. Lenstra JA, Groeneveld LF, Eding H, Kantanen J, Williams JL, Taberlet P, et al. Molecular tools and analytical approaches for the characterization of farm animal genetic diversity. Anim Genet. 2012;43(5):483–502. pmid:22497351.

18. Colli L, Lancioni H, Cardinali I, Olivieri A, Capodiferro MR, Pellecchia M, et al. Whole mitochondrial genomes unveil the impact of domestication on goat matrilineal variability. BMC Genomics. 2015;16(1):1115. pmid:26714643; PubMed Central PMCID: PMCPMC4696231.

19. Scheu A, Powell A, Bollongino R, Vigne JD, Tresset A, Çakırlar C, et al. The genetic prehistory of domesticated cattle from their origin to the spread across Europe. BMC Genet. 2015;16:54. pmid:26018295; PubMed Central PMCID: PMCPMC4445560.

20. Di Lorenzo P, Lancioni H, Ceccobelli S, Curcio L, Panella F, et al. Uniparental genetic systems: a male and a female perspective in the domestic cattle origin and evolution. Electronic Journal of Biotechnology. 2016;Sept;23:69–78.

21. Lai SJ, Liu YP, Liu YX, Li XW, Yao YG. Genetic diversity and origin of Chinese cattle revealed by mtDNA D-loop sequence variation. Mol Phylogenet Evol. 2006;38(1):146–54. Epub 2005/07/28. pmid:16054846.

22. Dadi H, Tibbo M, Takahashi Y, Nomura K, Hanada H, Amano T. Variation in mitochondrial DNA and maternal genetic ancestry of Ethiopian cattle populations. Anim Genet. 2009;40(4):556–9. Epub 2009/04/03. pmid:19397526.

23. Hristov P, Spassov N, Iliev N, Radoslavov G. An independent event of Neolithic cattle domestication on the South-eastern Balkans: evidence from prehistoric aurochs and cattle populations. Mitochondrial DNA A DNA Mapp Seq Anal. 2017;28(3):383–91. Epub 2015/12/29. pmid:26711535.

24. Achilli A, Bonfiglio S, Olivieri A, Malusà A, Pala M, Hooshiar Kashani B, et al. The multifaceted origin of taurine cattle reflected by the mitochondrial genome. PLoS One. 2009;4(6):e5753. pmid:19484124; PubMed Central PMCID: PMCPMC2684589.

25. Bonfiglio S, Ginja C, De Gaetano A, Achilli A, Olivieri A, Colli L, et al. Origin and spread of Bos taurus: new clues from mitochondrial genomes belonging to haplogroup T1. PLoS One. 2012;7(6):e38601. pmid:22685589; PubMed Central PMCID: PMCPMC3369859.

26. Olivieri A, Gandini F, Achilli A, Fichera A, Rizzi E, Bonfiglio S, et al. Mitogenomes from Egyptian Cattle Breeds: New Clues on the Origin of Haplogroup Q and the Early Spread of Bos taurus from the Near East. PLoS One. 2015;10(10):e0141170. pmid:26513361.

27. Moioli B, Napolitano F, Catillo G. Genetic diversity between Piedmontese, Maremmana, and Podolica cattle breeds. J Hered. 2004;95(3):250–6. pmid:15220392.

28. Pariset L, Mariotti M, Nardone A, Soysal MI, Ozkan E, Williams JL, et al. Relationships between Podolic cattle breeds assessed by single nucleotide polymorphisms (SNPs) genotyping. J Anim Breed Genet. 2010;127(6):481–8. Epub 2010/10/28. pmid:21077972.

29. Gargani M, Pariset L, Lenstra JA, De Minicis E, Valentini A, Consortium ECGD. Microsatellite genotyping of medieval cattle from central Italy suggests an old origin of Chianina and Romagnola cattle. Front Genet. 2015;6:68. pmid:25788902; PubMed Central PMCID: PMCPMC4349168.

30. Upadhyay MR, Chen W, Lenstra JA, Goderie CR, MacHugh DE, Park SD, et al. Genetic origin, admixture and population history of aurochs (Bos primigenius) and primitive European cattle. Heredity (Edinb). 2017;118(2):169–76. Epub 2016/09/28. pmid:27677498; PubMed Central PMCID: PMCPMC5234481.

31. Lenstra J, Ajmone , Marsan P, Beja Pereira A, Bollongino R, Bradley DG, et al. Meta-Analysis of Mitochondrial DNA Reveals Several Population Bottlenecks during Worldwide Migrations of Cattle. Diversity. 2014;6 178–87.

32. Cerezo M, Achilli A, Olivieri A, Perego UA, Gómez-Carballa A, Brisighelli F, et al. Reconstructing ancient mitochondrial DNA links between Africa and Europe. Genome Res. 2012;22(5):821–6. pmid:22454235; PubMed Central PMCID: PMCPMC3337428.

33. Ascunce S, Kitchen A, Schmidt P, Miyamoto M, Mulligan C. An Unusual Pattern of Ancient Mitochondrial DNA Haplogroups in Northern African Cattle. Zool Stud. 2007;46(1):123–5.

34. Mason I. A world dictionary of livestock breeds, types and varieties. 4th ed. 1996.

35. Gómez-Carballa A, Pardo-Seco J, Amigo J, Martinón-Torres F, Salas A. Mitogenomes from The 1000 Genome Project reveal new Near Eastern features in present-day Tuscans. PLoS One. 2015;10(3):e0119242. Epub 2015/03/18. pmid:25786119; PubMed Central PMCID: PMCPMC4365045.

36. Achilli A, Olivieri A, Pala M, Metspalu E, Fornarino S, Battaglia V, et al. Mitochondrial DNA variation of modern Tuscans supports the near eastern origin of Etruscans. Am J Hum Genet. 2007;80(4):759–68. Epub 2007/02/06. pmid:17357081; PubMed Central PMCID: PMCPMC1852723.

37. Brisighelli F, Capelli C, Alvarez-Iglesias V, Onofri V, Paoli G, Tofanelli S, et al. The Etruscan timeline: a recent Anatolian connection. Eur J Hum Genet. 2009;17(5):693–6. Epub 2008/12/03. pmid:19050723; PubMed Central PMCID: PMCPMC2986270.

38. Pardo-Seco J, Gómez-Carballa A, Amigo J, Martinón-Torres F, Salas A. A genome-wide study of modern-day Tuscans: revisiting Herodotus's theory on the origin of the Etruscans. PLoS One. 2014;9(9):e105920. Epub 2014/09/17. pmid:25230205; PubMed Central PMCID: PMCPMC4167696.

39. Rognoni G, Pagnacco G. Atlante Etnografico delle popolazioni bovine allevate in Italia. Rome: Consiglio Nazionale delle Ricerche; 1983.

40. Achilli A, Olivieri A, Pellecchia M, Uboldi C, Colli L, Al-Zahery N, et al. Mitochondrial genomes of extinct aurochs survive in domestic cattle. Curr Biol. 2008;18(4):R157–8. pmid:18302915.

41. Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2. Epub 2009/04/03. pmid:19346325.

42. Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online. 2007;1:47–50. Epub 2007/02/23. pmid:19325852; PubMed Central PMCID: PMCPMC2658868.

43. Lindgren G, Backström N, Swinburne J, Hellborg L, Einarsson A, Sandberg K, et al. Limited number of patrilines in horse domestication. Nat Genet. 2004;36(4):335–6. pmid:15034578.

44. Cardinali I, Lancioni H, Giontella A, Capodiferro MR, Capomaccio S, Buttazzoni L, et al. An Overview of Ten Italian Horse Breeds through Mitochondrial DNA. PLoS One. 2016;11(4):e0153004. Epub 2016/04/07. pmid:27054850; PubMed Central PMCID: PMCPMC4824442.