Serveur d'exploration sur le cobalt au Maghreb

Attention, ce site est en cours de développement !
Attention, site généré par des moyens informatiques à partir de corpus bruts.
Les informations ne sont donc pas validées.

Phylogeny and Patterns of Diversity of Goat mtDNA Haplogroup A Revealed by Resequencing Complete Mitogenomes

Identifieur interne : 000458 ( Pmc/Corpus ); précédent : 000457; suivant : 000459

Phylogeny and Patterns of Diversity of Goat mtDNA Haplogroup A Revealed by Resequencing Complete Mitogenomes

Auteurs : Maria Grazia Doro ; Daniela Piras ; Giovanni Giuseppe Leoni ; Giuseppina Casu ; Simona Vaccargiu ; Debora Parracciani ; Salvatore Naitana ; Mario Pirastu ; Andrea Novelletto

Source :

RBID : PMC:3999278

Abstract

We sequenced to near completion the entire mtDNA of 28 Sardinian goats, selected to represent the widest possible diversity of the most widespread mitochondrial evolutionary lineage, haplogroup (Hg) A. These specimens were reporters of the diversity in the island but also elsewhere, as inferred from their affiliation to each of 11 clades defined by D-loop variation. Two reference sequences completed the dataset. Overall, 206 variations were found in the full set of 30 sequences, of which 23 were protein-coding non-synonymous single nucleotide substitutions. Many polymorphic sites within Hg A were informative for the reconstruction of its internal phylogeny. Bayesian and network clustering revealed a general similarity over the entire molecule of sequences previously assigned to the same D-loop clade, indicating evolutionarily meaningful lineages. Two major sister groupings emerged within Hg A, which parallel distinct geographical distributions of D-loop clades in extant stocks. The pattern of variation in protein-coding genes revealed an overwhelming role of purifying selection, with the quota of surviving variants approaching neutrality. However, a simple model of relaxation of selection for the bulk of variants here reported should be rejected. Non-synonymous diversity of Hg's A, B and C denoted that a proportion of variants not greater than that allowed in the wild was given the opportunity to spread into domesticated stocks. Our results also confirmed that a remarkable proportion of pre-existing Hg A diversity became incorporated into domestic stocks. Our results confirm clade A11 as a well differentiated and ancient lineage peculiar of Sardinia.


Url:
DOI: 10.1371/journal.pone.0095969
PubMed: 24763315
PubMed Central: 3999278

Links to Exploration step

PMC:3999278

Le document en format XML

<record>
<TEI>
<teiHeader>
<fileDesc>
<titleStmt>
<title xml:lang="en">Phylogeny and Patterns of Diversity of Goat mtDNA Haplogroup A Revealed by Resequencing Complete Mitogenomes</title>
<author>
<name sortKey="Doro, Maria Grazia" sort="Doro, Maria Grazia" uniqKey="Doro M" first="Maria Grazia" last="Doro">Maria Grazia Doro</name>
<affiliation>
<nlm:aff id="aff1">
<addr-line>Institute of Population Genetics, National Research Council (CNR), Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Piras, Daniela" sort="Piras, Daniela" uniqKey="Piras D" first="Daniela" last="Piras">Daniela Piras</name>
<affiliation>
<nlm:aff id="aff1">
<addr-line>Institute of Population Genetics, National Research Council (CNR), Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="aff2">
<addr-line>Experimental Zoo-prophylactic Institute of Sardinia, Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Leoni, Giovanni Giuseppe" sort="Leoni, Giovanni Giuseppe" uniqKey="Leoni G" first="Giovanni Giuseppe" last="Leoni">Giovanni Giuseppe Leoni</name>
<affiliation>
<nlm:aff id="aff3">
<addr-line>Faculty of Veterinary Sciences, University of Sassari, Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Casu, Giuseppina" sort="Casu, Giuseppina" uniqKey="Casu G" first="Giuseppina" last="Casu">Giuseppina Casu</name>
<affiliation>
<nlm:aff id="aff1">
<addr-line>Institute of Population Genetics, National Research Council (CNR), Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Vaccargiu, Simona" sort="Vaccargiu, Simona" uniqKey="Vaccargiu S" first="Simona" last="Vaccargiu">Simona Vaccargiu</name>
<affiliation>
<nlm:aff id="aff1">
<addr-line>Institute of Population Genetics, National Research Council (CNR), Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Parracciani, Debora" sort="Parracciani, Debora" uniqKey="Parracciani D" first="Debora" last="Parracciani">Debora Parracciani</name>
<affiliation>
<nlm:aff id="aff4">
<addr-line>Genetic Park of Ogliastra, Perdasdefogu, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Naitana, Salvatore" sort="Naitana, Salvatore" uniqKey="Naitana S" first="Salvatore" last="Naitana">Salvatore Naitana</name>
<affiliation>
<nlm:aff id="aff3">
<addr-line>Faculty of Veterinary Sciences, University of Sassari, Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Pirastu, Mario" sort="Pirastu, Mario" uniqKey="Pirastu M" first="Mario" last="Pirastu">Mario Pirastu</name>
<affiliation>
<nlm:aff id="aff1">
<addr-line>Institute of Population Genetics, National Research Council (CNR), Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="aff4">
<addr-line>Genetic Park of Ogliastra, Perdasdefogu, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Novelletto, Andrea" sort="Novelletto, Andrea" uniqKey="Novelletto A" first="Andrea" last="Novelletto">Andrea Novelletto</name>
<affiliation>
<nlm:aff id="aff5">
<addr-line>Department of Biology, University “Tor Vergata”, Rome, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt>
<idno type="wicri:source">PMC</idno>
<idno type="pmid">24763315</idno>
<idno type="pmc">3999278</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3999278</idno>
<idno type="RBID">PMC:3999278</idno>
<idno type="doi">10.1371/journal.pone.0095969</idno>
<date when="2014">2014</date>
<idno type="wicri:Area/Pmc/Corpus">000458</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000458</idno>
</publicationStmt>
<sourceDesc>
<biblStruct>
<analytic>
<title xml:lang="en" level="a" type="main">Phylogeny and Patterns of Diversity of Goat mtDNA Haplogroup A Revealed by Resequencing Complete Mitogenomes</title>
<author>
<name sortKey="Doro, Maria Grazia" sort="Doro, Maria Grazia" uniqKey="Doro M" first="Maria Grazia" last="Doro">Maria Grazia Doro</name>
<affiliation>
<nlm:aff id="aff1">
<addr-line>Institute of Population Genetics, National Research Council (CNR), Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Piras, Daniela" sort="Piras, Daniela" uniqKey="Piras D" first="Daniela" last="Piras">Daniela Piras</name>
<affiliation>
<nlm:aff id="aff1">
<addr-line>Institute of Population Genetics, National Research Council (CNR), Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="aff2">
<addr-line>Experimental Zoo-prophylactic Institute of Sardinia, Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Leoni, Giovanni Giuseppe" sort="Leoni, Giovanni Giuseppe" uniqKey="Leoni G" first="Giovanni Giuseppe" last="Leoni">Giovanni Giuseppe Leoni</name>
<affiliation>
<nlm:aff id="aff3">
<addr-line>Faculty of Veterinary Sciences, University of Sassari, Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Casu, Giuseppina" sort="Casu, Giuseppina" uniqKey="Casu G" first="Giuseppina" last="Casu">Giuseppina Casu</name>
<affiliation>
<nlm:aff id="aff1">
<addr-line>Institute of Population Genetics, National Research Council (CNR), Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Vaccargiu, Simona" sort="Vaccargiu, Simona" uniqKey="Vaccargiu S" first="Simona" last="Vaccargiu">Simona Vaccargiu</name>
<affiliation>
<nlm:aff id="aff1">
<addr-line>Institute of Population Genetics, National Research Council (CNR), Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Parracciani, Debora" sort="Parracciani, Debora" uniqKey="Parracciani D" first="Debora" last="Parracciani">Debora Parracciani</name>
<affiliation>
<nlm:aff id="aff4">
<addr-line>Genetic Park of Ogliastra, Perdasdefogu, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Naitana, Salvatore" sort="Naitana, Salvatore" uniqKey="Naitana S" first="Salvatore" last="Naitana">Salvatore Naitana</name>
<affiliation>
<nlm:aff id="aff3">
<addr-line>Faculty of Veterinary Sciences, University of Sassari, Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Pirastu, Mario" sort="Pirastu, Mario" uniqKey="Pirastu M" first="Mario" last="Pirastu">Mario Pirastu</name>
<affiliation>
<nlm:aff id="aff1">
<addr-line>Institute of Population Genetics, National Research Council (CNR), Sassari, Italy</addr-line>
</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="aff4">
<addr-line>Genetic Park of Ogliastra, Perdasdefogu, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Novelletto, Andrea" sort="Novelletto, Andrea" uniqKey="Novelletto A" first="Andrea" last="Novelletto">Andrea Novelletto</name>
<affiliation>
<nlm:aff id="aff5">
<addr-line>Department of Biology, University “Tor Vergata”, Rome, Italy</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series>
<title level="j">PLoS ONE</title>
<idno type="eISSN">1932-6203</idno>
<imprint>
<date when="2014">2014</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc>
<textClass></textClass>
</profileDesc>
</teiHeader>
<front>
<div type="abstract" xml:lang="en">
<p>We sequenced to near completion the entire mtDNA of 28 Sardinian goats, selected to represent the widest possible diversity of the most widespread mitochondrial evolutionary lineage, haplogroup (Hg) A. These specimens were reporters of the diversity in the island but also elsewhere, as inferred from their affiliation to each of 11 clades defined by D-loop variation. Two reference sequences completed the dataset. Overall, 206 variations were found in the full set of 30 sequences, of which 23 were protein-coding non-synonymous single nucleotide substitutions. Many polymorphic sites within Hg A were informative for the reconstruction of its internal phylogeny. Bayesian and network clustering revealed a general similarity over the entire molecule of sequences previously assigned to the same D-loop clade, indicating evolutionarily meaningful lineages. Two major sister groupings emerged within Hg A, which parallel distinct geographical distributions of D-loop clades in extant stocks. The pattern of variation in protein-coding genes revealed an overwhelming role of purifying selection, with the quota of surviving variants approaching neutrality. However, a simple model of relaxation of selection for the bulk of variants here reported should be rejected. Non-synonymous diversity of Hg's A, B and C denoted that a proportion of variants not greater than that allowed in the wild was given the opportunity to spread into domesticated stocks. Our results also confirmed that a remarkable proportion of pre-existing Hg A diversity became incorporated into domestic stocks. Our results confirm clade A11 as a well differentiated and ancient lineage peculiar of Sardinia.</p>
</div>
</front>
<back>
<div1 type="bibliography">
<listBibl>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Diamond, J" uniqKey="Diamond J">J Diamond</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Zeder, Ma" uniqKey="Zeder M">MA Zeder</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Luikart, G" uniqKey="Luikart G">G Luikart</name>
</author>
<author>
<name sortKey="Gielly, L" uniqKey="Gielly L">L Gielly</name>
</author>
<author>
<name sortKey="Excoffier, L" uniqKey="Excoffier L">L Excoffier</name>
</author>
<author>
<name sortKey="Vigne, J D" uniqKey="Vigne J">J-D Vigne</name>
</author>
<author>
<name sortKey="Bouvet, J" uniqKey="Bouvet J">J Bouvet</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Naderi, S" uniqKey="Naderi S">S Naderi</name>
</author>
<author>
<name sortKey="Rezaei, H R" uniqKey="Rezaei H">H-R Rezaei</name>
</author>
<author>
<name sortKey="Taberlet, P" uniqKey="Taberlet P">P Taberlet</name>
</author>
<author>
<name sortKey="Zundel, S" uniqKey="Zundel S">S Zundel</name>
</author>
<author>
<name sortKey="Rafat, S A" uniqKey="Rafat S">S-A Rafat</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Naderi, S" uniqKey="Naderi S">S Naderi</name>
</author>
<author>
<name sortKey="Rezaei, H R" uniqKey="Rezaei H">H-R Rezaei</name>
</author>
<author>
<name sortKey="Pompanon, F" uniqKey="Pompanon F">F Pompanon</name>
</author>
<author>
<name sortKey="Blum, Mgb" uniqKey="Blum M">MGB Blum</name>
</author>
<author>
<name sortKey="Negrini, R" uniqKey="Negrini R">R Negrini</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hughes, S" uniqKey="Hughes S">S Hughes</name>
</author>
<author>
<name sortKey="Fernandez, H" uniqKey="Fernandez H">H Fernàndez</name>
</author>
<author>
<name sortKey="Cucchi, T" uniqKey="Cucchi T">T Cucchi</name>
</author>
<author>
<name sortKey="Duffraisse, M" uniqKey="Duffraisse M">M Duffraisse</name>
</author>
<author>
<name sortKey="Casabianca, F" uniqKey="Casabianca F">F Casabianca</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Taberlet, P" uniqKey="Taberlet P">P Taberlet</name>
</author>
<author>
<name sortKey="Valentini, A" uniqKey="Valentini A">A Valentini</name>
</author>
<author>
<name sortKey="Rezaei, Hr" uniqKey="Rezaei H">HR Rezaei</name>
</author>
<author>
<name sortKey="Naderi, S" uniqKey="Naderi S">S Naderi</name>
</author>
<author>
<name sortKey="Pompanon, F" uniqKey="Pompanon F">F Pompanon</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Taberlet, P" uniqKey="Taberlet P">P Taberlet</name>
</author>
<author>
<name sortKey="Coissac, E" uniqKey="Coissac E">E Coissac</name>
</author>
<author>
<name sortKey="Pansu, J" uniqKey="Pansu J">J Pansu</name>
</author>
<author>
<name sortKey="Pompanon, F" uniqKey="Pompanon F">F Pompanon</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Achilli, A" uniqKey="Achilli A">A Achilli</name>
</author>
<author>
<name sortKey="Bonfiglio, S" uniqKey="Bonfiglio S">S Bonfiglio</name>
</author>
<author>
<name sortKey="Olivieri, A" uniqKey="Olivieri A">A Olivieri</name>
</author>
<author>
<name sortKey="Malusa, A" uniqKey="Malusa A">A Malusà</name>
</author>
<author>
<name sortKey="Pala, M" uniqKey="Pala M">M Pala</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Lancioni, H" uniqKey="Lancioni H">H Lancioni</name>
</author>
<author>
<name sortKey="Di Lorenzo, P" uniqKey="Di Lorenzo P">P Di Lorenzo</name>
</author>
<author>
<name sortKey="Ceccobelli, S" uniqKey="Ceccobelli S">S Ceccobelli</name>
</author>
<author>
<name sortKey="Perego, Ua" uniqKey="Perego U">UA Perego</name>
</author>
<author>
<name sortKey="Miglio, A" uniqKey="Miglio A">A Miglio</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Piras, D" uniqKey="Piras D">D Piras</name>
</author>
<author>
<name sortKey="Doro, Mg" uniqKey="Doro M">MG Doro</name>
</author>
<author>
<name sortKey="Casu, G" uniqKey="Casu G">G Casu</name>
</author>
<author>
<name sortKey="Melis, Mp" uniqKey="Melis M">MP Melis</name>
</author>
<author>
<name sortKey="Vaccargiu, S" uniqKey="Vaccargiu S">S Vaccargiu</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Pereira, F" uniqKey="Pereira F">F Pereira</name>
</author>
<author>
<name sortKey="Queir S, S" uniqKey="Queir S S">S Queiròs</name>
</author>
<author>
<name sortKey="Gusmao, L" uniqKey="Gusmao L">L Gusmào</name>
</author>
<author>
<name sortKey="Nijman, Ij" uniqKey="Nijman I">IJ Nijman</name>
</author>
<author>
<name sortKey="Cuppen, E" uniqKey="Cuppen E">E Cuppen</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hughes, Al" uniqKey="Hughes A">AL Hughes</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bjornerfeldt, S" uniqKey="Bjornerfeldt S">S Björnerfeldt</name>
</author>
<author>
<name sortKey="Webster, Mt" uniqKey="Webster M">MT Webster</name>
</author>
<author>
<name sortKey="Vila, C" uniqKey="Vila C">C Vilà</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Nomura, K" uniqKey="Nomura K">K Nomura</name>
</author>
<author>
<name sortKey="Yonezawa, T" uniqKey="Yonezawa T">T Yonezawa</name>
</author>
<author>
<name sortKey="Mano, S" uniqKey="Mano S">S Mano</name>
</author>
<author>
<name sortKey="Kawakami, S" uniqKey="Kawakami S">S Kawakami</name>
</author>
<author>
<name sortKey="Shedlock, Am" uniqKey="Shedlock A">AM Shedlock</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Vacca, Gm" uniqKey="Vacca G">GM Vacca</name>
</author>
<author>
<name sortKey="Daga, C" uniqKey="Daga C">C Daga</name>
</author>
<author>
<name sortKey="Pazzola, M" uniqKey="Pazzola M">M Pazzola</name>
</author>
<author>
<name sortKey="Carcangiu, V" uniqKey="Carcangiu V">V Carcangiu</name>
</author>
<author>
<name sortKey="Dettori, Ml" uniqKey="Dettori M">ML Dettori</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Mcdonald, Jh" uniqKey="Mcdonald J">JH McDonald</name>
</author>
<author>
<name sortKey="Kreitman, M" uniqKey="Kreitman M">M Kreitman</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Zeder, Ma" uniqKey="Zeder M">MA Zeder</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Sakuta, R" uniqKey="Sakuta R">R Sakuta</name>
</author>
<author>
<name sortKey="Goto, Y" uniqKey="Goto Y">Y Goto</name>
</author>
<author>
<name sortKey="Nonaka, I" uniqKey="Nonaka I">I Nonaka</name>
</author>
<author>
<name sortKey="Horai, S" uniqKey="Horai S">S Horai</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Pereira, L" uniqKey="Pereira L">L Pereira</name>
</author>
<author>
<name sortKey="Soares, P" uniqKey="Soares P">P Soares</name>
</author>
<author>
<name sortKey="Radivojac, P" uniqKey="Radivojac P">P Radivojac</name>
</author>
<author>
<name sortKey="Li, B" uniqKey="Li B">B Li</name>
</author>
<author>
<name sortKey="Samuels, Dc" uniqKey="Samuels D">DC Samuels</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Fernandez, H" uniqKey="Fernandez H">H Fernàndez</name>
</author>
<author>
<name sortKey="Hughes, S" uniqKey="Hughes S">S Hughes</name>
</author>
<author>
<name sortKey="Vigne, J D" uniqKey="Vigne J">J-D Vigne</name>
</author>
<author>
<name sortKey="Helmer, D" uniqKey="Helmer D">D Helmer</name>
</author>
<author>
<name sortKey="Hodgins, G" uniqKey="Hodgins G">G Hodgins</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Amills, M" uniqKey="Amills M">M Amills</name>
</author>
<author>
<name sortKey="Ramirez, O" uniqKey="Ramirez O">O Ramírez</name>
</author>
<author>
<name sortKey="Tomas, A" uniqKey="Tomas A">A Tomàs</name>
</author>
<author>
<name sortKey="Badaoui, B" uniqKey="Badaoui B">B Badaoui</name>
</author>
<author>
<name sortKey="Marmi, J" uniqKey="Marmi J">J Marmi</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Nabholz, B" uniqKey="Nabholz B">B Nabholz</name>
</author>
<author>
<name sortKey="Glemin, S" uniqKey="Glemin S">S Glèmin</name>
</author>
<author>
<name sortKey="Galtier, N" uniqKey="Galtier N">N Galtier</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Soares, P" uniqKey="Soares P">P Soares</name>
</author>
<author>
<name sortKey="Abrantes, D" uniqKey="Abrantes D">D Abrantes</name>
</author>
<author>
<name sortKey="Rito, T" uniqKey="Rito T">T Rito</name>
</author>
<author>
<name sortKey="Thomson, N" uniqKey="Thomson N">N Thomson</name>
</author>
<author>
<name sortKey="Radivojac, P" uniqKey="Radivojac P">P Radivojac</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Woolley, S" uniqKey="Woolley S">S Woolley</name>
</author>
<author>
<name sortKey="Johnson, J" uniqKey="Johnson J">J Johnson</name>
</author>
<author>
<name sortKey="Smith, Mj" uniqKey="Smith M">MJ Smith</name>
</author>
<author>
<name sortKey="Crandall, Ka" uniqKey="Crandall K">KA Crandall</name>
</author>
<author>
<name sortKey="Mcclellan, Da" uniqKey="Mcclellan D">DA McClellan</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Ballard, Jwo" uniqKey="Ballard J">JWO Ballard</name>
</author>
<author>
<name sortKey="Melvin, Rg" uniqKey="Melvin R">RG Melvin</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kang, L" uniqKey="Kang L">L Kang</name>
</author>
<author>
<name sortKey="Zheng, H X" uniqKey="Zheng H">H-X Zheng</name>
</author>
<author>
<name sortKey="Chen, F" uniqKey="Chen F">F Chen</name>
</author>
<author>
<name sortKey="Yan, S" uniqKey="Yan S">S Yan</name>
</author>
<author>
<name sortKey="Liu, K" uniqKey="Liu K">K Liu</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Brandano, P" uniqKey="Brandano P">P Brandano</name>
</author>
<author>
<name sortKey="Piras, B" uniqKey="Piras B">B Piras</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Drummond, Aj" uniqKey="Drummond A">AJ Drummond</name>
</author>
<author>
<name sortKey="Rambaut, A" uniqKey="Rambaut A">A Rambaut</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hassanin, A" uniqKey="Hassanin A">A Hassanin</name>
</author>
<author>
<name sortKey="Bonillo, C" uniqKey="Bonillo C">C Bonillo</name>
</author>
<author>
<name sortKey="Nguyen, Bx" uniqKey="Nguyen B">BX Nguyen</name>
</author>
<author>
<name sortKey="Cruaud, C" uniqKey="Cruaud C">C Cruaud</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bandelt, Hj" uniqKey="Bandelt H">HJ Bandelt</name>
</author>
<author>
<name sortKey="Forster, P" uniqKey="Forster P">P Forster</name>
</author>
<author>
<name sortKey="Rohl, A" uniqKey="Rohl A">A Rohl</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kumar, S" uniqKey="Kumar S">S Kumar</name>
</author>
<author>
<name sortKey="Nei, M" uniqKey="Nei M">M Nei</name>
</author>
<author>
<name sortKey="Dudley, J" uniqKey="Dudley J">J Dudley</name>
</author>
<author>
<name sortKey="Tamura, K" uniqKey="Tamura K">K Tamura</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Thompson, Jd" uniqKey="Thompson J">JD Thompson</name>
</author>
<author>
<name sortKey="Gibson, T" uniqKey="Gibson T">T Gibson</name>
</author>
<author>
<name sortKey="Plewniak, F" uniqKey="Plewniak F">F Plewniak</name>
</author>
<author>
<name sortKey="Jeanmougin, F" uniqKey="Jeanmougin F">F Jeanmougin</name>
</author>
<author>
<name sortKey="Higgins, Dg" uniqKey="Higgins D">DG Higgins</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article">
<pmc-dir>properties open_access</pmc-dir>
<front>
<journal-meta>
<journal-id journal-id-type="nlm-ta">PLoS One</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group>
<journal-title>PLoS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher>
<publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="pmid">24763315</article-id>
<article-id pub-id-type="pmc">3999278</article-id>
<article-id pub-id-type="publisher-id">PONE-D-14-02342</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0095969</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Research Article</subject>
</subj-group>
<subj-group subj-group-type="Discipline-v2">
<subject>Biology and Life Sciences</subject>
<subj-group>
<subject>Agriculture</subject>
<subj-group>
<subject>Animal Management</subject>
</subj-group>
</subj-group>
<subj-group>
<subject>Evolutionary Biology</subject>
<subj-group>
<subject>Evolutionary Processes</subject>
<subj-group>
<subject>Natural Selection</subject>
</subj-group>
</subj-group>
<subj-group>
<subject>Population Genetics</subject>
<subj-group>
<subject>Haplotypes</subject>
</subj-group>
</subj-group>
<subj-group>
<subject>Evolutionary Genetics</subject>
</subj-group>
</subj-group>
<subj-group>
<subject>Genetics</subject>
<subj-group>
<subject>Animal Genetics</subject>
</subj-group>
</subj-group>
<subj-group>
<subject>Population Biology</subject>
</subj-group>
<subj-group>
<subject>Veterinary Science</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Phylogeny and Patterns of Diversity of Goat mtDNA Haplogroup A Revealed by Resequencing Complete Mitogenomes</article-title>
<alt-title alt-title-type="running-head">Goat mtDNA Haplogroup A</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Doro</surname>
<given-names>Maria Grazia</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Piras</surname>
<given-names>Daniela</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Leoni</surname>
<given-names>Giovanni Giuseppe</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Casu</surname>
<given-names>Giuseppina</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Vaccargiu</surname>
<given-names>Simona</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Parracciani</surname>
<given-names>Debora</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Naitana</surname>
<given-names>Salvatore</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Pirastu</surname>
<given-names>Mario</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Novelletto</surname>
<given-names>Andrea</given-names>
</name>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
<xref ref-type="corresp" rid="cor1">
<sup>*</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<label>1</label>
<addr-line>Institute of Population Genetics, National Research Council (CNR), Sassari, Italy</addr-line>
</aff>
<aff id="aff2">
<label>2</label>
<addr-line>Experimental Zoo-prophylactic Institute of Sardinia, Sassari, Italy</addr-line>
</aff>
<aff id="aff3">
<label>3</label>
<addr-line>Faculty of Veterinary Sciences, University of Sassari, Sassari, Italy</addr-line>
</aff>
<aff id="aff4">
<label>4</label>
<addr-line>Genetic Park of Ogliastra, Perdasdefogu, Italy</addr-line>
</aff>
<aff id="aff5">
<label>5</label>
<addr-line>Department of Biology, University “Tor Vergata”, Rome, Italy</addr-line>
</aff>
<contrib-group>
<contrib contrib-type="editor">
<name>
<surname>Yao</surname>
<given-names>Yong-Gang</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1">
<addr-line>Kunming Institute of Zoology, Chinese Academy of Sciences, China</addr-line>
</aff>
<author-notes>
<corresp id="cor1">* E-mail:
<email>novelletto@bio.uniroma2.it</email>
</corresp>
<fn fn-type="conflict">
<p>
<bold>Competing Interests: </bold>
The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con">
<p>Conceived and designed the experiments: MGD D. Piras MP. Performed the experiments: MGD D. Piras GC. Analyzed the data: MGD D. Piras AN. Contributed reagents/materials/analysis tools: GGL SV D. Parracciani SN. Wrote the paper: AN.</p>
</fn>
</author-notes>
<pub-date pub-type="collection">
<year>2014</year>
</pub-date>
<pub-date pub-type="epub">
<day>24</day>
<month>4</month>
<year>2014</year>
</pub-date>
<volume>9</volume>
<issue>4</issue>
<elocation-id>e95969</elocation-id>
<history>
<date date-type="received">
<day>16</day>
<month>1</month>
<year>2014</year>
</date>
<date date-type="accepted">
<day>1</day>
<month>4</month>
<year>2014</year>
</date>
</history>
<permissions>
<copyright-year>2014</copyright-year>
<copyright-holder>Doro et al</copyright-holder>
<license>
<license-p>This is an open-access article distributed under the terms of the
<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>
, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<abstract>
<p>We sequenced to near completion the entire mtDNA of 28 Sardinian goats, selected to represent the widest possible diversity of the most widespread mitochondrial evolutionary lineage, haplogroup (Hg) A. These specimens were reporters of the diversity in the island but also elsewhere, as inferred from their affiliation to each of 11 clades defined by D-loop variation. Two reference sequences completed the dataset. Overall, 206 variations were found in the full set of 30 sequences, of which 23 were protein-coding non-synonymous single nucleotide substitutions. Many polymorphic sites within Hg A were informative for the reconstruction of its internal phylogeny. Bayesian and network clustering revealed a general similarity over the entire molecule of sequences previously assigned to the same D-loop clade, indicating evolutionarily meaningful lineages. Two major sister groupings emerged within Hg A, which parallel distinct geographical distributions of D-loop clades in extant stocks. The pattern of variation in protein-coding genes revealed an overwhelming role of purifying selection, with the quota of surviving variants approaching neutrality. However, a simple model of relaxation of selection for the bulk of variants here reported should be rejected. Non-synonymous diversity of Hg's A, B and C denoted that a proportion of variants not greater than that allowed in the wild was given the opportunity to spread into domesticated stocks. Our results also confirmed that a remarkable proportion of pre-existing Hg A diversity became incorporated into domestic stocks. Our results confirm clade A11 as a well differentiated and ancient lineage peculiar of Sardinia.</p>
</abstract>
<funding-group>
<funding-statement>Work was supported by grants Istituto di Genetica delle Popolazioni - Progetti di ricerca e innovazione al di sotto della soglia de minimis - Programma Operativo Regionale Sardegna 2000/2006 Asse 3 Risorse Umane Misura 3.1.3. – Ricerca e sviluppo tecnologico nelle imprese e territorio – “Sviluppo di un sistema informatico integrato per l'identificazione, la valorizzazione e la conservazione del patrimonio caprino in Sardegna” to M.P.; L. R. 7 agosto 2007, n. 7 - Promozione della ricerca scientifica e dell' innovazione tecnologica in Sardegna (Bando 2008) – “Diversità genetica, conservazione e diffusione del germoplasma di specie di interesse faunistico presenti in Sardegna” to G.G.L.; Programmi di Ricerca Scientifica di Rilevante Interesse Nazionale - Ministero dell' Istruzione dell' Università e della Ricerca 2012JA4BTY_003 to A.N. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
<counts>
<page-count count="11"></page-count>
</counts>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>The domestication of the goat (
<italic>Capra hircus</italic>
) has been one of the major landmarks in the evolution of human modes of subsistence
<xref rid="pone.0095969-CluttonBrock1" ref-type="bibr">[1]</xref>
,
<xref rid="pone.0095969-Diamond1" ref-type="bibr">[2]</xref>
,
<xref rid="pone.0095969-Zeder1" ref-type="bibr">[3]</xref>
. Besides its valuable dairy products, its portability has made the goat the key element to ensure a stable human settlement in otherwise inhospitable habitats.</p>
<p>The genetic traces of domestication in this species have been first addressed
<xref rid="pone.0095969-Luikart1" ref-type="bibr">[4]</xref>
by analyzing the diversity of the mitochondrial D-loop among individuals collected worldwide. These authors showed that the matrilines of the current breeds coalesce to three founders of the so-called haplogroups (Hg's) A, B and C, indicative of a limited number but not a singular drawing of female founders from wild populations. Later, the number of haplogroups grew to 6 (A, B, C, D, F, G
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
), with the newly described Hg F as the most basal one. Each of these found close matches in the alleged ancestor of the domestic goat, the bezoar (
<italic>C. aegagrus</italic>
)
<xref rid="pone.0095969-Naderi2" ref-type="bibr">[6]</xref>
.</p>
<p>While Hg A is the most derived lineage, it is also the most widespread and displays a fast and markedly star-shaped radiation, to be interpreted as the result of the success of animals bred in maternal isolation from wild population(s). Given the presence of this Hg in many breeds currently grown worldwide, its phylogeography could potentially be the essential piece of information to reconstruct the spread of maternal lineages. In particular it would help in distinguishing between ancient and recent displacements of stocks. This seems particularly useful as a further characterization of some locales, in addition to the assignment of a reasonable window of time for their introduction in a given area. This additional information can also be an added value whenever such locales are particularly adapted to a specific habitat, are valued as a recreational resource, or when the labeling of the products as “genuinely local” may enhance their economic relevance, thus calling for specific protection and conservation actions
<xref rid="pone.0095969-Hughes1" ref-type="bibr">[7]</xref>
,
<xref rid="pone.0095969-Taberlet1" ref-type="bibr">[8]</xref>
,
<xref rid="pone.0095969-Taberlet2" ref-type="bibr">[9]</xref>
,
<xref rid="pone.0095969-Achilli1" ref-type="bibr">[10]</xref>
,
<xref rid="pone.0095969-Lancioni1" ref-type="bibr">[11]</xref>
.</p>
<p>In this context, however, the internal phylogeny of haplogroup A is poorly understood, mainly because of the hypermutability of the D-loop of mtDNA, which causes many positions to display the same nucleotide by the repeated occurrence of the same mutational event (identity by state) rather than by descent. The net result is a set of inextricable reticulations when one attempts to reconstruct internal, evolutionarily meaningful lineages. As an intermediate step forward, we
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
took advantage of the demography underlying a large proportion of the population of the Italian island of Sardinia. By sampling a large number of animals we identified, within Hg A, subgroups of D-loop sequences related in a star-like fashion, reminiscent of common descent from a limited number of founders. We then defined 11 “ad interim” clades of Hg A, and listed a reduced number of mutations, among those occurring in the D-loop, that can distinguish them. Their frequencies ranged 1.8–21% among >1,500 Hg A individuals. By comparing our results with those of other screenings
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
,
<xref rid="pone.0095969-Pereira1" ref-type="bibr">[13]</xref>
we also realized that, with the exception of two clades (A-9 and A-11), the same groupings accounted for a large proportion of the stocks grown elsewhere. While this result can be explained by a relatively recent (historical) import of matrilines (possibly conveyed by different breeds) in Sardinia, the presence of peculiar lineages might be due to a much earlier introduction. These lineages might have been rare even in their homeland at the time of introduction in Sardinia and now lost outside the Island, or they might be the result of new mutations that occurred and were selected
<italic>in loco.</italic>
</p>
<p>In summary, Sardinia recapitulates the D-loop phylogeny of a non trivial portion of Hg A and thus offers the opportunity of working out a clearer picture of divergence within the entire haplogroup. We then selected animals carrier of each of the 11 Hg A clades, plus an Hg C outgroup, and sequenced their mtDNA to near completion. In this way we could test our initial hypothesis that D-loop clades are predictors of evolutionarily meaningful lineages. This question is of general interest as far as the D-loop is the DNA region most commonly addressed to initially approach a number of issues of species' biology (see e.g. Chapter 21 in ref.
<xref rid="pone.0095969-Frankham1" ref-type="bibr">[14]</xref>
). Also, finding a large number of variants along the entire mtDNA sequence can identify major evolutionary forces affecting livestock populations after domestication. This includes a better understanding of the selective regime(s) that impacted the mitochondrial genome. We then asked: i) whether our set of variants fitted a model of relaxed selection allowing for otherwise mildly deleterious alleles, as hypothesized for other species
<xref rid="pone.0095969-Hughes2" ref-type="bibr">[15]</xref>
,
<xref rid="pone.0095969-Bjrnerfeldt1" ref-type="bibr">[16]</xref>
, and ii) whether a reduced number of mutational changes or assemblages of mutational changes was subjected to positive selection, possibly due to their potential physiological relevance.</p>
<p>In the context of the resulting phylogenetic reconstruction, the clade A-11 of Hg A emerges as an ancient lineage, currently peculiar of Sardinia, characterized by many non-synonymous substitutions.</p>
</sec>
<sec id="s2">
<title>Results</title>
<p>The sequence coverage attained in the 28 individuals selected for complete mtDNA resequencing (Table S1 in
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
) is summarized in
<xref ref-type="supplementary-material" rid="pone.0095969.s001">Figure S1</xref>
, with colors indicating Hg A clade affiliation. With the exception of a single clade A8 individual (id. 304) sequence coverage was always above 80%. For all clades at least one individual was covered at >89%. For the single Hg C individual only a 336 bp portion of the sequence was missing (98% coverage).</p>
<p>Only two regions remained poorly covered in several individuals, i.e. those at the beginning of amplicon 3 and at the end of amplicon 31. In the first case the presence of a C
<sub>10</sub>
stretch often caused slippage that generated overlapping shifted profiles further downstream (this problem was not observed in carriers of an interrupting T at pos. 1114). In the second case the complete coverage of the large amplicon 31 was attempted by sequencing in both directions, but profiles obtained with the reverse primer turned out to be often short or of poor quality and then discarded.</p>
<p>As one of our aims was to explore the phylogenetic reliability of the relationships between Hg A clades defined by D-loop sequences, our entire analysis took into account only variations in the region spanning reference positions 1–15,430 (i.e. outside the D-loop). In addition to our Hg C individual (id. 985), all our analyses included sequences GU068049 (representative of Hg A from Mongolia) and GU295658 (representative of Hg B from South East Asia), thus reaching a total of 30 sequences hereafter referred to as the “full set”.</p>
<p>Overall, 206 variations were found in the full set of 30 sequences (Table S2 in
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
). Of these, 113 were contributed by the Hg C and B sequences, whereas 93 were found among Hg A sequences. This confirmed Hg's C and B as strongly divergent lineages, but the number of residual variable positions denoted that many polymorphic sites within Hg A could be informative for the reconstruction of its internal phylogeny. In fact, 59 of the 93 variations were singletons in Hg A individuals and the remaining 34 were shared by two or more individuals. Eleven variable positions were private for the reference sequence GU068049, 7 of which residing in regions covered in ref.
<xref rid="pone.0095969-Nomura1" ref-type="bibr">[17]</xref>
and also reported therein.</p>
<p>Among single nucleotide substitutions 193 and 8 were transitions and transversions, respectively. Among transitions T/C changes outnumbered A/G changes (146 vs. 47). As to indels, two consecutive thymine residues reported for the reference sequence GU068049, were not found in any other sequence and were considered as two additional adjacent private single nucleotide indels (positions 177, 178). Two variations consisted in length changes in low complexity A
<sub>8</sub>
(pos. 11560-7) and C
<sub>10</sub>
(pos. 1111-20) stretches. In particular this latter segment was C
<sub>11</sub>
in Hg B but shortened and affected by other substitutions in Hg C. Finally, an A insertion was found in a single individual, between positions 1238 and 1239.</p>
<sec id="s2a">
<title>Reconstructing Hg A phylogeny</title>
<p>We used two independent methods to obtain a reconstruction of affinities between mtDNA sequences. The first method was a Bayesian tree reconstruction. The results are displayed in
<xref ref-type="supplementary-material" rid="pone.0095969.s002">Figure S2</xref>
. The relative order of branching of Hg's C and B replicates that of refs.
<xref rid="pone.0095969-Luikart1" ref-type="bibr">[4]</xref>
,
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
,
<xref rid="pone.0095969-Nomura1" ref-type="bibr">[17]</xref>
with highly supported nodes, and identifies Hg C as the most basal lineage. The next node (support 100%) defines the radiation of all Hg A lineages and separates the reference sequence GU068049 from the rest of Hg A sequences. A series of closely spaced nodes then separates small groupings. Each of these is formed by sequences clustered together with high support and previously assigned to the same D-loop clade. The only exceptions are represented by an A5 sequence (id. 1495) clustering with two A7 sequences, another A5 (id. 331) defining a single branch outside the clade of the remaining A5 sequences, and the two A8 sequences (id's 830 and 304) which did not cluster together. Of the three sequences not previously assigned to any of the D-loop clades, two (id's 1332 and 480) were highly similar in the D-loop and clustered together in the tree, while the third (id. 692) was divergent in the D-loop and clustered with A9 in the tree.</p>
<p>In summary, this clustering scheme denoted a general similarity over the entire molecule of sequences previously assigned to the same D-loop clade. As to the relationships among clades, these remained less strongly supported and partly contradicted our previous reconstruction (see
<xref ref-type="fig" rid="pone-0095969-g002">Fig. 2</xref>
in
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
), most likely because of the star-like and rapid divergence between them, and the presence of nucleotide states generated by recurrent mutational events (see below). The 27 Sardinian Hg A sequences here examined clustered into two sister groupings downstream to a poorly (0.21) supported node. However we noticed two features of this partition that may strengthen its significance. First, it separated the sequences in the same manner as a major hiatus in the D-loop network (clades A1 to A6 plus A9 in
<xref ref-type="supplementary-material" rid="pone.0095969.s002">Fig. S2</xref>
, bottom vs. A7 to A11 in
<xref ref-type="supplementary-material" rid="pone.0095969.s002">Fig. S2</xref>
, top) and second, it had a clear geographic correlate (
<xref ref-type="supplementary-material" rid="pone.0095969.s002">Fig. S2</xref>
). In particular, the sequences of the first group turned out to be associated with D-loop clades represented mainly in Europe when inspecting three datasets
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
,
<xref rid="pone.0095969-Pereira1" ref-type="bibr">[13]</xref>
,
<xref rid="pone.0095969-Vacca1" ref-type="bibr">[18]</xref>
, with only a minority of observations from Asia. Conversely, the sequences of the second group appeared to be associated with three clades (A7, A8 and A10) abundant throughout Asia. Sequences corresponding to D-loop clade A11 were assigned to this latter group (
<xref ref-type="supplementary-material" rid="pone.0095969.s002">Fig. S2</xref>
, top). Indeed a single instance of the haplotypic arrangement of D-loop variants characteristic of clade A11 was represented among the 2488 individuals in the datasets of refs.
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
,
<xref rid="pone.0095969-Pereira1" ref-type="bibr">[13]</xref>
, and it was sampled in Italy. Conversely, its common presence in Sardinia was confirmed in an independent study
<xref rid="pone.0095969-Vacca1" ref-type="bibr">[18]</xref>
.</p>
<p>We also constructed a network representing the arrangement of allele states among our sequences (
<xref ref-type="fig" rid="pone-0095969-g001">Fig. 1</xref>
). Since we excluded some positions, the topology of the network turned out to be simpler than the Bayesian tree. Nevertheless, the branching order of Hg's C, B and A was replicated (
<xref ref-type="fig" rid="pone-0095969-g001">Fig. 1</xref>
, inset), with a large multifurcation of all Hg A sequences. The network also replicated the clustering of sequences according to their D-loop clade affiliation already evident in the Bayesian tree, with a single exception represented by id. 304 (clade A8) which now appeared at the tip of a long branch distal to an A5 sequence (
<xref ref-type="fig" rid="pone-0095969-g001">Fig. 1</xref>
, bottom right). Interestingly, one short branch of the multifurcation led to a node basal to all sequences of clades A8, A7, A10 and A11 (plus A5, id 1495), i.e. one of the major groupings in the Bayesian tree. The network construction identified 7 positions with recurrent mutational events (Table S2 in
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
).</p>
<fig id="pone-0095969-g001" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0095969.g001</object-id>
<label>Figure 1</label>
<caption>
<title>Median joining network of 30 sequences showing the occurrence of 23 non-synonymous substitutions.</title>
<p>Each circle represents a sequence. Branch length is proportional to the number of mutations. Note that two A11 and two uncharacterized sequences are lumped into single (larger) nodes (see Materials and Methods). The haplogroup or D-loop clade (for Hg A) affiliation of clusters of sequences is shown in large empty characters. Aminoacid replacements are shown beside each branch in their polarized form. The Hg A portion of the network is magnified while the original network is in the inset. A large lozenge indicates the likely position of the root and was added in the drawing to partition three variants peculiar of Hg C. Their polarization was obtained by comparison with
<italic>C. ibex</italic>
.</p>
</caption>
<graphic xlink:href="pone.0095969.g001"></graphic>
</fig>
</sec>
<sec id="s2b">
<title>Exploring possible selection regimes</title>
<p>We searched for genetic traces left by selection regimes in the whole set of lineages here considered by analyzing separately protein-coding and RNA-coding genes.</p>
<p>As to protein-coding genes, we argued that if variants in our set were largely neutral, they should be found in all genes and in numbers proportional to the targetable positions in each gene. This expectation was verified, with a linear and highly correlated increase of variants as a function of gene length, both including (R
<sup>2</sup>
 = 0.71; p = 3×10
<sup>−4</sup>
) and excluding (R
<sup>2</sup>
 = 0.52; p = 5×10
<sup>−3</sup>
) the large number of variants contributed by the Hg B and C sequences (
<xref ref-type="fig" rid="pone-0095969-g002">Fig. 2A</xref>
). This observation was in large part attributable to variants in third codon positions, which accumulated linearly (
<xref ref-type="fig" rid="pone-0095969-g002">Fig. 2B</xref>
) with gene length in both the full set of 30 sequences (R
<sup>2</sup>
 = 0.81; p = 3×10
<sup>−5</sup>
) as well as in the 28 Hg A sequences (R
<sup>2</sup>
 = 0.51; p = 6×10
<sup>−3</sup>
). In both cases
<italic>coxIII</italic>
displayed an unusual accumulation of variants.</p>
<fig id="pone-0095969-g002" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0095969.g002</object-id>
<label>Figure 2</label>
<caption>
<title>Scatterplot of number of variants recorded as a function of length of each of the 13 protein-coding mtDNA genes.</title>
<p>Filled dots: variants scored among all 30 sequences considered. Empty dots: variants scored among the 28 Hg A sequences. Panel A: all positions. Panel B: only 3
<sup>rd</sup>
codon positions. Panel C: only non-synonymous variants (note that the three black dots at 0 variants overlap and obscure empty dots).
<italic>coxIII</italic>
is indicated by an arrow.</p>
</caption>
<graphic xlink:href="pone.0095969.g002"></graphic>
</fig>
<p>The linearity between the number of variants and gene length was somewhat disrupted for non-synonymous variants (
<xref ref-type="fig" rid="pone-0095969-g002">Fig. 2C</xref>
), with low correlations for both the full set of 30 sequences (R
<sup>2</sup>
 = 0.06, p = 0.41) and the 28 Hg A sequences (R
<sup>2</sup>
 = 0.08, p = 0.33). The highest absolute numbers of non-synonymous variants were in
<italic>coxIII</italic>
and
<italic>atp6</italic>
(total length  = 280 and 226 codons, respectively) whereas the relatively long
<italic>coxI</italic>
(514 codons) harboured no variants.</p>
<p>In order to distinguish whether this result was due to stochastic variations in the small numbers of non-synonymous variants or differential selection across genes, we analyzed the data in a McDonald-Kreitman framework
<xref rid="pone.0095969-McDonald1" ref-type="bibr">[19]</xref>
. For each gene we then compared non-synonymous and synonymous variations found among Hg's A, B and C (polymorphic) or only in Hg A, with those recorded in the divergence with
<italic>C. ibex</italic>
. The results are shown in Table S3 in
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
. For none of the genes the McDonald-Kreitman test turned out to be significant, either considering the three major haplogroups or only Hg A.</p>
<p>Thus, in summary, the absolute number of non-synonymous as compared to synonymous variants was in line with a generalized purifying selective regime on the former ones, possibly with variable intensities across genes, as exemplified by the accumulation above average in
<italic>coxIII</italic>
(21.7% of non-synonymous variants). However, when compared with the wild outgroup, the proportions of the two types of variants showed that, quantitatively, the intensity of purification remained unaltered under domestication.</p>
<p>The above results could not rule out that only a small subset of non-synonymous variants were impacted by the breeding practices. In fact, under domestication, selection on these variants may have become less stringent or null (see e.g.
<xref rid="pone.0095969-Hughes2" ref-type="bibr">[15]</xref>
), possibly as an effect of breeding practices that favored the survival of females and drove the demography of matrilines
<xref rid="pone.0095969-Zeder2" ref-type="bibr">[20]</xref>
. A second hypothesis is that diversifying selection operated by breeders maintained a repertoire of valuable phenotypes influenced by some of the non synonymous variants.</p>
<p>We analyzed the occurrence of non-synonymous variants in the phylogeny, to rule out possible cumulative effects of variants with subtle phenotypic effects. Of the 23 non-synonymous substitutions observed in the full set of 30 sequences (
<xref ref-type="table" rid="pone-0095969-t001">Table 1</xref>
), 3 were peculiar of Hg C (
<xref ref-type="fig" rid="pone-0095969-g001">Fig. 1</xref>
, inset). For one of these Hg C carried the derived aminoacid state (
<italic>atp6</italic>
62S), whereas for the remaining two Hg C carried the ancestral state (
<italic>coxIII</italic>
168T and
<italic>nd3</italic>
29A) and the derived alleles were shared by Hg's B and A (
<italic>coxIII</italic>
168M and
<italic>nd3</italic>
29V). The excess of 7 non-synonymous vs. 20 synonymous changes on the branch leading to Hg B (matching the increased ω ratio reported in this Hg reported in
<xref rid="pone.0095969-Nomura1" ref-type="bibr">[17]</xref>
), was not significant (p = 0.06; Fisher exact test) as compared to the rest of variants in protein-coding genes.</p>
<table-wrap id="pone-0095969-t001" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0095969.t001</object-id>
<label>Table 1</label>
<caption>
<title>Comparative analysis of the 23 non-synonymous substitutions in protein-coding genes inferred in the full set of 30 sequences.</title>
</caption>
<alternatives>
<graphic id="pone-0095969-t001-1" xlink:href="pone.0095969.t001"></graphic>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Position</td>
<td align="left" rowspan="1" colspan="1">Gene</td>
<td align="left" rowspan="1" colspan="1">Allelic variants</td>
<td align="left" rowspan="1" colspan="1">aa. replacement</td>
<td align="left" rowspan="1" colspan="1">Homologous aa. in rCRS</td>
<td align="left" rowspan="1" colspan="1">Variants in MITOMAP
<xref ref-type="table-fn" rid="nt101">1</xref>
</td>
<td align="left" rowspan="1" colspan="1">Occurrences in GB</td>
<td align="left" rowspan="1" colspan="1">MutPred pathogenicity score in human
<xref ref-type="table-fn" rid="nt102">2</xref>
</td>
<td align="left" rowspan="1" colspan="1">mtDNA selection score in human
<xref ref-type="table-fn" rid="nt102">2</xref>
</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">2767</td>
<td align="left" rowspan="1" colspan="1">
<italic>nd1</italic>
</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">F9L</td>
<td align="left" rowspan="1" colspan="1">L</td>
<td align="left" rowspan="1" colspan="1">
<italic>L-L</italic>
</td>
<td align="left" rowspan="1" colspan="1">23</td>
<td align="left" rowspan="1" colspan="1">0.625</td>
<td align="left" rowspan="1" colspan="1">0,768</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">4587</td>
<td align="left" rowspan="1" colspan="1">
<italic>nd2</italic>
</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">T227I</td>
<td align="left" rowspan="1" colspan="1">T</td>
<td align="left" rowspan="1" colspan="1">T-A</td>
<td align="left" rowspan="1" colspan="1">1</td>
<td align="left" rowspan="1" colspan="1">0.586</td>
<td align="left" rowspan="1" colspan="1">0.661</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">7053</td>
<td align="left" rowspan="1" colspan="1">
<italic>coxII</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">T13A</td>
<td align="left" rowspan="1" colspan="1">T</td>
<td align="left" rowspan="1" colspan="1">
<italic>T-T</italic>
</td>
<td align="left" rowspan="1" colspan="1">122</td>
<td align="left" rowspan="1" colspan="1">0.378</td>
<td align="left" rowspan="1" colspan="1">0.297</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">7185</td>
<td align="left" rowspan="1" colspan="1">
<italic>coxII</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">D57N</td>
<td align="left" rowspan="1" colspan="1">D</td>
<td align="left" rowspan="1" colspan="1">
<bold>D-N</bold>
, D-A, D-G, D-D</td>
<td align="left" rowspan="1" colspan="1">10, 1, 1, 3</td>
<td align="left" rowspan="1" colspan="1">0.499</td>
<td align="left" rowspan="1" colspan="1">0.473</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">8056</td>
<td align="left" rowspan="1" colspan="1">
<italic>atp6</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">V43I</td>
<td align="left" rowspan="1" colspan="1">I</td>
<td align="left" rowspan="1" colspan="1">
<bold>I-V</bold>
, I-T, I-I</td>
<td align="left" rowspan="1" colspan="1">9, 5, 1127</td>
<td align="left" rowspan="1" colspan="1">0.325</td>
<td align="left" rowspan="1" colspan="1">0.242</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">8114</td>
<td align="left" rowspan="1" colspan="1">
<italic>atp6</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">N62S</td>
<td align="left" rowspan="1" colspan="1">N</td>
<td align="left" rowspan="1" colspan="1">
<bold>N-S</bold>
</td>
<td align="left" rowspan="1" colspan="1">6</td>
<td align="left" rowspan="1" colspan="1">0.526</td>
<td align="left" rowspan="1" colspan="1">0.525</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">8134</td>
<td align="left" rowspan="1" colspan="1">
<italic>atp6</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">T69A</td>
<td align="left" rowspan="1" colspan="1">S</td>
<td align="left" rowspan="1" colspan="1">
<italic>T-T</italic>
</td>
<td align="left" rowspan="1" colspan="1">7</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">8596</td>
<td align="left" rowspan="1" colspan="1">
<italic>atp6</italic>
</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">H223Y</td>
<td align="left" rowspan="1" colspan="1">H</td>
<td align="left" rowspan="1" colspan="1">H-R</td>
<td align="left" rowspan="1" colspan="1">3</td>
<td align="left" rowspan="1" colspan="1">0.656</td>
<td align="left" rowspan="1" colspan="1">0.865</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">8704</td>
<td align="left" rowspan="1" colspan="1">
<italic>coxIII</italic>
</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">I32T</td>
<td align="left" rowspan="1" colspan="1">A</td>
<td align="left" rowspan="1" colspan="1">A-V, A-A</td>
<td align="left" rowspan="1" colspan="1">8, 26</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">8731</td>
<td align="left" rowspan="1" colspan="1">
<italic>coxIII</italic>
</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">A41V</td>
<td align="left" rowspan="1" colspan="1">T</td>
<td align="left" rowspan="1" colspan="1">T-A, T-M, T-T</td>
<td align="left" rowspan="1" colspan="1">7, 1, 33</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">9015</td>
<td align="left" rowspan="1" colspan="1">
<italic>coxIII</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">V136I</td>
<td align="left" rowspan="1" colspan="1">V</td>
<td align="left" rowspan="1" colspan="1">
<italic>V-V</italic>
</td>
<td align="left" rowspan="1" colspan="1">8</td>
<td align="left" rowspan="1" colspan="1">n.a.
<xref ref-type="table-fn" rid="nt103">3</xref>
</td>
<td align="left" rowspan="1" colspan="1">n.a.
<xref ref-type="table-fn" rid="nt103">3</xref>
</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">9112</td>
<td align="left" rowspan="1" colspan="1">
<italic>coxIII</italic>
</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">M168T</td>
<td align="left" rowspan="1" colspan="1">L</td>
<td align="left" rowspan="1" colspan="1">L-L, L-S</td>
<td align="left" rowspan="1" colspan="1">17, 1</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">9369</td>
<td align="left" rowspan="1" colspan="1">
<italic>coxIII</italic>
</td>
<td align="left" rowspan="1" colspan="1">C/G</td>
<td align="left" rowspan="1" colspan="1">V254L</td>
<td align="left" rowspan="1" colspan="1">V</td>
<td align="left" rowspan="1" colspan="1">V-I, V-V</td>
<td align="left" rowspan="1" colspan="1">117, 9</td>
<td align="left" rowspan="1" colspan="1">0.471</td>
<td align="left" rowspan="1" colspan="1">0.425</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">9548</td>
<td align="left" rowspan="1" colspan="1">
<italic>nd3</italic>
</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">V29A</td>
<td align="left" rowspan="1" colspan="1">G</td>
<td align="left" rowspan="1" colspan="1">G-S</td>
<td align="left" rowspan="1" colspan="1">66</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">10433</td>
<td align="left" rowspan="1" colspan="1">
<italic>nd4</italic>
</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">P89S</td>
<td align="left" rowspan="1" colspan="1">L</td>
<td align="left" rowspan="1" colspan="1">L-P, L-L</td>
<td align="left" rowspan="1" colspan="1">64, 4</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">10493</td>
<td align="left" rowspan="1" colspan="1">
<italic>nd4</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">T109A</td>
<td align="left" rowspan="1" colspan="1">T</td>
<td align="left" rowspan="1" colspan="1">T-P,
<bold>T-A</bold>
, T-T</td>
<td align="left" rowspan="1" colspan="1">n.a., 109, 6</td>
<td align="left" rowspan="1" colspan="1">0.482</td>
<td align="left" rowspan="1" colspan="1">0.443</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">11872</td>
<td align="left" rowspan="1" colspan="1">
<italic>nd5</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">C42Y</td>
<td align="left" rowspan="1" colspan="1">S</td>
<td align="left" rowspan="1" colspan="1">
<italic>S-S</italic>
</td>
<td align="left" rowspan="1" colspan="1">3</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">13100</td>
<td align="left" rowspan="1" colspan="1">
<italic>nd5</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/T</td>
<td align="left" rowspan="1" colspan="1">I451M</td>
<td align="left" rowspan="1" colspan="1">L</td>
<td align="left" rowspan="1" colspan="1">
<italic>L-L</italic>
</td>
<td align="left" rowspan="1" colspan="1">5</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">13390</td>
<td align="left" rowspan="1" colspan="1">
<italic>nd5</italic>
</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">S548L</td>
<td align="left" rowspan="1" colspan="1">L</td>
<td align="left" rowspan="1" colspan="1">
<italic>L-L</italic>
</td>
<td align="left" rowspan="1" colspan="1">62</td>
<td align="left" rowspan="1" colspan="1">n.a.
<xref ref-type="table-fn" rid="nt103">3</xref>
</td>
<td align="left" rowspan="1" colspan="1">n.a.
<xref ref-type="table-fn" rid="nt103">3</xref>
</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">13964</td>
<td align="left" rowspan="1" colspan="1">
<italic>nd6</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">V39A</td>
<td align="left" rowspan="1" colspan="1">V38</td>
<td align="left" rowspan="1" colspan="1">V-L, V-I</td>
<td align="left" rowspan="1" colspan="1">2, 6</td>
<td align="left" rowspan="1" colspan="1">0.544</td>
<td align="left" rowspan="1" colspan="1">0.562</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">14479</td>
<td align="left" rowspan="1" colspan="1">
<italic>cytb</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/T</td>
<td align="left" rowspan="1" colspan="1">F109L</td>
<td align="left" rowspan="1" colspan="1">Y</td>
<td align="left" rowspan="1" colspan="1">Y-N, Y-H, Y-C</td>
<td align="left" rowspan="1" colspan="1">n.a., 33, 1</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">14795</td>
<td align="left" rowspan="1" colspan="1">
<italic>cytb</italic>
</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">T215A</td>
<td align="left" rowspan="1" colspan="1">S</td>
<td align="left" rowspan="1" colspan="1">
<italic>S-S</italic>
</td>
<td align="left" rowspan="1" colspan="1">19</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">14865</td>
<td align="left" rowspan="1" colspan="1">
<italic>cytb</italic>
</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">V238A</td>
<td align="left" rowspan="1" colspan="1">S</td>
<td align="left" rowspan="1" colspan="1">S-P, S-F, S-S</td>
<td align="left" rowspan="1" colspan="1">28, 7, 8</td>
<td align="left" rowspan="1" colspan="1">n.a</td>
<td align="left" rowspan="1" colspan="1">n.a.</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="nt101">
<label>1</label>
<p>Non synonymous variants are bolded. For positions in which only synonymous variants are known italics is used;</p>
</fn>
<fn id="nt102">
<label>2</label>
<p>from ref.
<xref rid="pone.0095969-Pereira2" ref-type="bibr">[22]</xref>
;</p>
</fn>
<fn id="nt103">
<label>3</label>
<p>same aminoacid encoded by a different codon in rCRS and goat;</p>
</fn>
<fn id="nt104">
<label></label>
<p>n.a.  =  not applicable.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>One substitution (
<italic>atp6</italic>
43I>V) was shared by all Hg A sequences but not Hg's B and C and was then assigned to the branch basal to Hg A, with the derived state in this Hg in agreement with previous results
<xref rid="pone.0095969-Nomura1" ref-type="bibr">[17]</xref>
.</p>
<p>Whithin Hg A, a non-synonymous substitution was found on each of the branches basal to sequences affiliated with D-loop clades A1, A2, A3, A7, A10 and A11, while two were on the branch leading to the highly similar uncharacterized sequences id. 1332 and 480. This co-occurrence of aminoacid substitutions in
<italic>atp6</italic>
and
<italic>coxIII</italic>
represents a non-significant excess (p = .38, Fisher exact test) as compared to 4 and 5 substitutions in each gene, respectively, over the 37 branches of the network. Additional private non-synonymous substitutions were found for sequences affiliated with clades A8 (id. 304), A3, A5 and A11. Finally, an additional non-synonymous substitution (
<italic>nd5</italic>
586L>Q) was found in one of the 5 supplementary clade A11 individuals sequenced in amplicon 28. In summary, a non-synonymous substitution was peculiar of all Hg A sequences and additional 12 affected 18 of the 28 Hg A sequences here considered, ruling out the selection of functionally highly diversified lineages.</p>
<p>We also reasoned that the two selective regimes outlined above lead to different expectations as to the composition of the repertoire of non-synonymous variants, with an enrichment of mildly deleterious variants under relaxed selection vs. an absence of detrimental variants under diversifying selection. In order to gain insights on the features of our array of non-synonymous variants, we explored the homologous positions in the largest collection of mtDNA mutations recorded so far for a single species, i.e. the MITOMAP database of
<italic>Homo sapiens</italic>
. For each DNA position this database reports the observed variant(s), the corresponding aminoacid change, the number of occurrences of each variant in GenBank (GB, to be used as a rough estimate of its degree of polymorphism) and the phenotypic consequences of each mutation, if available. The results are reported in columns 5–7 of
<xref ref-type="table" rid="pone-0095969-t001">Table 1</xref>
. In nearly half of cases (12/23) one of the aminoacids observed among our variants was conserved in the orthologous position of the human reference sequence (rCRS). Four of our aminoacid substitutions were replicated exactly in humans. Three of these had 10, 9 and 6 records in GB, respectively, i.e. values above the 60th percentile of frequencies among the 7932 variants. Interestingly the fourth (
<italic>nd4</italic>
T109A) was initially associated with the MELAS syndrome. However it was also found at high frequency in normal controls
<xref rid="pone.0095969-Sakuta1" ref-type="bibr">[21]</xref>
. Its causative role in the disease is now questioned (
<ext-link ext-link-type="uri" xlink:href="http://omim.org/entry/516003">http://omim.org/entry/516003</ext-link>
). Among our animals, this variation was found in only one of the two A3 individuals, and its frequency in the Sardinian stock at large remains to be determined. Overall, in the set of 13 protein-coding genes, the ratio of non-synonymous to synonymous variants recorded in MITOMAP approached that observed here only when rare variants were excluded (considering only variations found >40 times in GenBank,
<xref ref-type="supplementary-material" rid="pone.0095969.s003">Fig. S3</xref>
).</p>
<p>In 8 of the 23 human orthologous positions (34.7%) only synonymous variants were recorded, a finding potentially denoting that aminoacid replacements at these positions are not tolerated. Indeed this observation does not deviate from expectations, as far as the proportions of positions displaying only synonymous variations across the human orthologues of the 10 genes listed in
<xref ref-type="table" rid="pone-0095969-t001">Table 1</xref>
range between 21.9% in
<italic>atp6</italic>
and 63.7% in
<italic>nd4</italic>
.</p>
<p>Finally, we considered the indicators of possible pathogenicity computed in ref.
<xref rid="pone.0095969-Pereira2" ref-type="bibr">[22]</xref>
for all possible aminoacid substitutions in the human rCRS. The 10 substitutions that would replicate those observed in goats all display MutPred pathogenicity scores below the mode of the overall distribution and selection scores lower than 1, indicative of effects much milder than mutations recognized to be pathogenic. Even the
<italic>nd4</italic>
T109A, the sole substitution which appears in the list of 75 pathogenic mutations, has scores as low as 0.482 and 0.443 (
<xref ref-type="table" rid="pone-0095969-t001">Table 1</xref>
).</p>
<p>When feasible, we repeated the same type of comparative analyses for variants in RNA-coding genes (
<xref ref-type="table" rid="pone-0095969-t002">Table 2</xref>
). For variants in rRNA genes, 8 high-confidence orthologous positions displayed high frequency variants in human, too. The fact that for 5 positions there are no variants recorded in MITOMAP is within expectations, as far as this database contains only 271 and 438 variable positions for the 945 and 1569 bp of the 12S and 16S rRNA, respectively.</p>
<table-wrap id="pone-0095969-t002" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0095969.t002</object-id>
<label>Table 2</label>
<caption>
<title>Comparative analysis of the substitutions in RNA-coding genes inferred in the full set of 30 sequences.</title>
</caption>
<alternatives>
<graphic id="pone-0095969-t002-2" xlink:href="pone.0095969.t002"></graphic>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">Position</td>
<td align="left" rowspan="1" colspan="1">Gene</td>
<td align="left" rowspan="1" colspan="1">Allelic variants</td>
<td align="left" rowspan="1" colspan="1">Homologous pos. in rCRS</td>
<td align="left" rowspan="1" colspan="1">Variants in MITOMAP</td>
<td align="left" rowspan="1" colspan="1">Occurrences in GB</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">165</td>
<td align="left" rowspan="1" colspan="1">12S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">743C</td>
<td align="left" rowspan="1" colspan="1">not present</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">177</td>
<td align="left" rowspan="1" colspan="1">12S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/del</td>
<td align="left" rowspan="1" colspan="1">amb. align.
<xref ref-type="table-fn" rid="nt105">1</xref>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">178</td>
<td align="left" rowspan="1" colspan="1">12S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/del</td>
<td align="left" rowspan="1" colspan="1">amb. align.
<xref ref-type="table-fn" rid="nt105">1</xref>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">180</td>
<td align="left" rowspan="1" colspan="1">12S-rRNA</td>
<td align="left" rowspan="1" colspan="1">G/T</td>
<td align="left" rowspan="1" colspan="1">amb. align.
<xref ref-type="table-fn" rid="nt105">1</xref>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">181</td>
<td align="left" rowspan="1" colspan="1">12S-rRNA</td>
<td align="left" rowspan="1" colspan="1">A/T</td>
<td align="left" rowspan="1" colspan="1">amb. align.
<xref ref-type="table-fn" rid="nt105">1</xref>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">196</td>
<td align="left" rowspan="1" colspan="1">12S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">773T</td>
<td align="left" rowspan="1" colspan="1">T-C</td>
<td align="left" rowspan="1" colspan="1">1</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">205</td>
<td align="left" rowspan="1" colspan="1">12S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">782A</td>
<td align="left" rowspan="1" colspan="1">not present</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">213</td>
<td align="left" rowspan="1" colspan="1">12S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">790A</td>
<td align="left" rowspan="1" colspan="1">not present</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">605</td>
<td align="left" rowspan="1" colspan="1">12S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">1187T</td>
<td align="left" rowspan="1" colspan="1">T-C</td>
<td align="left" rowspan="1" colspan="1">18</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">944</td>
<td align="left" rowspan="1" colspan="1">12S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">1520T</td>
<td align="left" rowspan="1" colspan="1">T-C</td>
<td align="left" rowspan="1" colspan="1">13</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">1026</td>
<td align="left" rowspan="1" colspan="1">12S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">1601C</td>
<td align="left" rowspan="1" colspan="1">C-T</td>
<td align="left" rowspan="1" colspan="1">2</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">1114</td>
<td align="left" rowspan="1" colspan="1">16S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">amb. align.
<xref ref-type="table-fn" rid="nt105">1</xref>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">1111-20</td>
<td align="left" rowspan="1" colspan="1">16S-rRNA</td>
<td align="left" rowspan="1" colspan="1">C
<sub>10</sub>
/C
<sub>11</sub>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">1192</td>
<td align="left" rowspan="1" colspan="1">16S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">1764C</td>
<td align="left" rowspan="1" colspan="1">C-T</td>
<td align="left" rowspan="1" colspan="1">1</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">1238/39</td>
<td align="left" rowspan="1" colspan="1">16S-rRNA</td>
<td align="left" rowspan="1" colspan="1">A/del</td>
<td align="left" rowspan="1" colspan="1">amb. align.
<xref ref-type="table-fn" rid="nt105">1</xref>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">1369</td>
<td align="left" rowspan="1" colspan="1">16S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">amb. align.
<xref ref-type="table-fn" rid="nt105">1</xref>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">1671</td>
<td align="left" rowspan="1" colspan="1">16S-rRNA</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">2241A</td>
<td align="left" rowspan="1" colspan="1">not present</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">1850</td>
<td align="left" rowspan="1" colspan="1">16S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">2417C</td>
<td align="left" rowspan="1" colspan="1">C-G</td>
<td align="left" rowspan="1" colspan="1">4</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">2120</td>
<td align="left" rowspan="1" colspan="1">16S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">2686G</td>
<td align="left" rowspan="1" colspan="1">not present</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">2190</td>
<td align="left" rowspan="1" colspan="1">16S-rRNA</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">2756C</td>
<td align="left" rowspan="1" colspan="1">C-A, C-T</td>
<td align="left" rowspan="1" colspan="1">1, 5</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">2208</td>
<td align="left" rowspan="1" colspan="1">16S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">2774C</td>
<td align="left" rowspan="1" colspan="1">not present</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">2619</td>
<td align="left" rowspan="1" colspan="1">16S-rRNA</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">3184C</td>
<td align="left" rowspan="1" colspan="1">C-T</td>
<td align="left" rowspan="1" colspan="1">1</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">3771</td>
<td align="left" rowspan="1" colspan="1">tRNA-Gln</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">4335C</td>
<td align="left" rowspan="1" colspan="1">C-T</td>
<td align="left" rowspan="1" colspan="1">9</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">3892</td>
<td align="left" rowspan="1" colspan="1">tRNA-Met</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">4454T</td>
<td align="left" rowspan="1" colspan="1">T-A, T-C</td>
<td align="left" rowspan="1" colspan="1">90, 80</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">6894</td>
<td align="left" rowspan="1" colspan="1">tRNA-Ser (UCN)</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">7468C</td>
<td align="left" rowspan="1" colspan="1">C-T</td>
<td align="left" rowspan="1" colspan="1">5</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">7723</td>
<td align="left" rowspan="1" colspan="1">tRNA-Lys</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">8315A</td>
<td align="left" rowspan="1" colspan="1">A-G</td>
<td align="left" rowspan="1" colspan="1">1</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">9827</td>
<td align="left" rowspan="1" colspan="1">tRNA-Arg</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">10421C</td>
<td align="left" rowspan="1" colspan="1">C-T</td>
<td align="left" rowspan="1" colspan="1">1</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">9860</td>
<td align="left" rowspan="1" colspan="1">tRNA-Arg</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">10454T</td>
<td align="left" rowspan="1" colspan="1">T-C</td>
<td align="left" rowspan="1" colspan="1">59</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">9863</td>
<td align="left" rowspan="1" colspan="1">tRNA-Arg</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">amb. align.
<xref ref-type="table-fn" rid="nt105">1</xref>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">11560-7</td>
<td align="left" rowspan="1" colspan="1">tRNA-His</td>
<td align="left" rowspan="1" colspan="1">A
<sub>8</sub>
/A
<sub>9</sub>
</td>
<td align="left" rowspan="1" colspan="1">amb. align.
<xref ref-type="table-fn" rid="nt105">1</xref>
</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">11569</td>
<td align="left" rowspan="1" colspan="1">tRNA-His</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">12160A</td>
<td align="left" rowspan="1" colspan="1">not present</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">11624</td>
<td align="left" rowspan="1" colspan="1">tRNA-Ser (AGY)</td>
<td align="left" rowspan="1" colspan="1">T/C</td>
<td align="left" rowspan="1" colspan="1">12214C</td>
<td align="left" rowspan="1" colspan="1">not present</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">11691</td>
<td align="left" rowspan="1" colspan="1">tRNA-Leu(CUN)</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">12279A</td>
<td align="left" rowspan="1" colspan="1">A-T</td>
<td align="left" rowspan="1" colspan="1">1</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">11702</td>
<td align="left" rowspan="1" colspan="1">tRNA-Leu(CUN)</td>
<td align="left" rowspan="1" colspan="1">A/G</td>
<td align="left" rowspan="1" colspan="1">12290A</td>
<td align="left" rowspan="1" colspan="1">A-G</td>
<td align="left" rowspan="1" colspan="1">2</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="nt105">
<label>1</label>
<p>Ambiguous alignment.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Also the orthologues of variants residing in tRNA genes appeared to be highly polymorphic in human, with the replication of the alternative alleles in 7 out of 8 cases.</p>
</sec>
</sec>
<sec id="s3">
<title>Discussion</title>
<p>In this work we proceeded a step forward in the description of the phyletic relationships between lineages belonging to a major component of the goat mitochondrial diversity, Hg A. We obtained this by leveraging the diversity represented in a single Italian region, where the goat is highly valued and its diversity preserved.</p>
<p>We restricted our analysis to a portion of the molecule less prone to recurrent mutations than the highly variable D-loop (115 variable positions in 1213 bp across our 30 sequences). Nevertheless we showed that clusters of sequences defined by a selected subset of D-loop mutations
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
do match clusters defined by variants in the rest of the molecule. The topology connecting these clusters needs additional support, which can be obtained by completion of the typing of missing positions in our set. In particular, dubious assignments involve the lineages associated with the D-loop clades A5 and A8, which probably still contain a large amount of internal diversity. As to A9 we
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
already warned that this has to be considered mainly a paragroup, as it is defined only by ancestral alleles in the D-loop.</p>
<p>Also, the inclusion of high quality complete Hg A mitogenomes from outside Sardinia in future analyses is imperative. In fact, only three (7214T/C, 9221T/C and 14518A/G) of the variable positions here reported were also present among the 62 positions variable in Hg A animals sampled mainly in eastern and south-eastern Asia
<xref rid="pone.0095969-Nomura1" ref-type="bibr">[17]</xref>
, and allele states at two of these were probably homoplasic. Interestingly, also in this latter study the peculiarity of the reference sequence GU068049 (from the Inner Mongolian breed) emerged, with four private variants and three variants shared with a single individual (sample ID G1307). These findings qualify the trees of the two studies as complementary subsections of the entire Hg A phylogeny, but also predict that northern Asia may represent a further and poorly explored reservoir of goat mitochondrial diversity. Only when a robust topology will be attained it will be possible to reconsider the diversity of the D-loop and to enumerate the recurrent events within it.</p>
<p>With the above caveats in mind, a major branching among the Hg A sequences here reported separates lineages currently found mostly in western Eurasia from lineages common in the eastern Old World. Phylogeographic considerations in this species are weakened by its high portability by both land and sea since several millennia
<xref rid="pone.0095969-Fernndez1" ref-type="bibr">[23]</xref>
, which may have led to the seeding of widely distant regions of the world in a manner disconnected from phyletic affinities (see e.g.
<xref rid="pone.0095969-Amills1" ref-type="bibr">[24]</xref>
). Despite this, the present study captured a partition probably reflecting an ancient split in the spread of stocks carrying Hg A from its likely homeland in Eastern Anatolia
<xref rid="pone.0095969-Naderi2" ref-type="bibr">[6]</xref>
. This hypothesis may serve to orient the search of markers in archeological material that documents that process. As far as genetics is concerned, D-loop clade assignment was feasible in at least 89% of 29 Hg A bezoars
<xref rid="pone.0095969-Naderi2" ref-type="bibr">[6]</xref>
and identified 17 and 8 carriers of clades A8 and A5, respectively (
<xref ref-type="table" rid="pone-0095969-t001">Table 1</xref>
), corresponding to each of the major Hg A tree branches. While the former animals were largely sampled in Iran (14 vs. 3 in Turkey) the latter ones were from Iran and Turkey in equal proportions, replicating to some extent the same east-west imbalance in the distribution of domestic goat lineages.</p>
<p>Our results also confirm the view that a remarkable proportion of pre-existing Hg A diversity became incorporated into domestic stocks. A date between 201 and 281 thousand years ago was proposed
<xref rid="pone.0095969-Luikart1" ref-type="bibr">[4]</xref>
for the coalescent of Hg's C, B and A. In our Bayesian tree Hg A occupies approximately 15% of the total tree height, corresponding to a date far predating the beginning of domestication. Using the average rate proposed
<xref rid="pone.0095969-Nabholz1" ref-type="bibr">[25]</xref>
for 3
<sup>rd</sup>
codon positions in wild Caprinae results in a similar estimate (
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
). This discrepancy between the age of the molecules and that of the domestic stocks testifies of a process of domestication not accompanied by strong bottlenecks and loss of matrilinear diversity, extending previous results obtained in the goat
<xref rid="pone.0095969-Naderi2" ref-type="bibr">[6]</xref>
,
<xref rid="pone.0095969-Nomura1" ref-type="bibr">[17]</xref>
. Speculations on which mechanism(s) produced this outcome include multiple drawings also for Hg A, from populations less differentiated that those harboring Hg's C and B. This could have occurred during a long period preceding full-blown domestication, in which semi-wild populations were simply managed
<xref rid="pone.0095969-Zeder1" ref-type="bibr">[3]</xref>
,
<xref rid="pone.0095969-Naderi2" ref-type="bibr">[6]</xref>
,
<xref rid="pone.0095969-Tudge1" ref-type="bibr">[26]</xref>
,
<xref rid="pone.0095969-Shipman1" ref-type="bibr">[27]</xref>
, without interrupting gene flow with wild relatives. It is worth noting that this represents one of the extremes of a full spectrum of possibilities. At the opposite end lies the taurine cattle (
<italic>Bos taurus</italic>
) for which an extremely narrow bottleneck was inferred based on mtDNA data
<xref rid="pone.0095969-Achilli1" ref-type="bibr">[10]</xref>
.</p>
<p>Overall, the pattern of variation in protein-coding genes revealed and confirmed an overwhelming role of purifying selection
<xref rid="pone.0095969-Hughes2" ref-type="bibr">[15]</xref>
,
<xref rid="pone.0095969-Nomura1" ref-type="bibr">[17]</xref>
, with the quota of surviving variants approaching neutrality. The features of the set of 23 non-synonymous variants did not provide any evidence for an enrichment in deleterious aminoacid replacements which have approached neutrality due to relaxed selection. Conversely, they replicated a group of variants which have been found to reach polymorphic frequencies in humans
<xref rid="pone.0095969-Pereira2" ref-type="bibr">[22]</xref>
. Interestingly, a large proportion of non-synonymous substitutions were on branches with two or more sequences downstream (internal branches), an additional evidence of their neutral or quasi-neutral behaviour. The overall proportions of synonymous variants here reported for Hg's A, B, C or Hg A only (0.865 and 0.838, respectively; Table S3 in
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
, bottom) closely approach the asymptotic values reported for a number of mammalian taxa in ref.
<xref rid="pone.0095969-Soares1" ref-type="bibr">[28]</xref>
. These estimate the constant corresponding to the observed proportion of synonymous mutations after the effect of purifying selection has become negligible. In our dataset the applicability of formal tests on the drasticity of aminaocid replacements (e.g. as implemented in TreeSAAP,
<xref rid="pone.0095969-Woolley1" ref-type="bibr">[29]</xref>
) was limited by the reduced number of variants in each gene.</p>
<p>A similar argument applies to RNA-coding genes, as well. Here we found a number of variants per kilobase (11.5, 6.9 and 8.0 for 12S rRNA, 16S rRNA and tRNA's, respectively) even lower than in protein-coding genes (15.0), in line with the strong purifying selection inferred in other species
<xref rid="pone.0095969-Hughes2" ref-type="bibr">[15]</xref>
. However, the few observed variants have polymorphic homologous variants in humans.</p>
<p>In summary, a simple model of relaxation of selection for the bulk of variants here reported should be rejected. Rather, our analysis of non-synonymous diversity (which considered intra- vs. inter-species instead of intra- vs. inter-haplogroup
<xref rid="pone.0095969-Nomura1" ref-type="bibr">[17]</xref>
comparisons) revealed that a proportion of variants not greater that that allowed in the wild was given the opportunity to spread into domesticated stocks. This also excludes a regime of generalized diversifying selection. Mitochondrial phenotypes have been claimed to play a role in a number of life-history traits
<xref rid="pone.0095969-Ballard1" ref-type="bibr">[30]</xref>
,
<xref rid="pone.0095969-Kang1" ref-type="bibr">[31]</xref>
. A reduction of metabolic efficiency was hypothesized to be favourable in a carnivore, the dog
<xref rid="pone.0095969-Bjrnerfeldt1" ref-type="bibr">[16]</xref>
, but this seems hardly transferable to a ruminant, whose valued products are milk, fat and lean body mass and a high reproductive rate. Of main interest here are instead temperature and starvation resistance, in light of the new environmental conditions encountered by the goats along with their breeders' migrations. An approach to reveal possible positive selection on individual variants will be to measure diversity indexes in several sequences from the same clade. This has not been feasible here, due to the highly biased criterion for inclusion of individuals in the study.</p>
<p>Our results now confirmed that clade A11, (previously recognized as peculiar of Sardinia
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
,
<xref rid="pone.0095969-Vacca1" ref-type="bibr">[18]</xref>
where it reaches frequencies as high as 8.6%), indeed identifies a lineage that is well differentiated along the entire mtDNA molecule and is not simply a derivative of one of the lineages confined to Europe. Based on the current distribution of mtDNA lineages mostly related to it, the possible geographical range for the origins of clade A11 sequences is much wider than Europe alone. Moreover, based on the number of mutations distinguishing it, its antiquity is comparable to that of much more common lineages currently spread from Europe to East Asia. This suggests that goats carriers of this lineage were relocated into the Island from populations that are now lost or not yet explored for mtDNA diversity.</p>
</sec>
<sec sec-type="materials|methods" id="s4">
<title>Materials and Methods</title>
<sec id="s4a">
<title>The individuals</title>
<p>All activities carried out by the Regional Association of Sardinian Farmers (ARAS) officers, part of which related to this study, followed ethical guidelines for care and use of agricultural animals for research (EC Directive 86/609/EEC). For this study specific approval by a review board was not necessary, as none of the procedures used here met the criteria to define them experiments as defined in Article 2 of the cited Directive. All animals of the present study belong to the series previously studied
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
. They were not re-sampled and the new data were generated on previously collected DNA, without any additional burden to the animal.</p>
<p>Twenty-eight individuals were selected to represent all the clades reported in
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
, based on D-loop diversity. In order to represent also the within-clade diversity, individuals belonging to the same clade were included only if they differed at some positions in the D-loop. Three individuals were included specifically to represent the unclassified D-loop sequences.</p>
<p>Additional criteria for inclusion in the above group of 28 individuals were: a general phenotype fitting that of the “Sarda” breed
<xref rid="pone.0095969-Brandano1" ref-type="bibr">[32]</xref>
, and a geographical origin preferentially in the Ogliastra sub-region of Sardinia (see
<xref ref-type="fig" rid="pone-0095969-g001">Fig. 1</xref>
in
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
).</p>
<p>In order to confirm/dismiss a subset of mutations, additional 38 individuals belonging to D-loop clades A2, A4, A5, A7, A8 and A11 were considered and sequenced for amplicons 10, 22 and 28 (see below). This group, too, was selected to display the “Sarda” phenotype, but the geographical criterion was relaxed and representatives of other sub-regions were included. The essential data are reported in Table S1 in
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
.</p>
</sec>
<sec id="s4b">
<title>Amplification and sequencing</title>
<p>The entire mtDNA was subdivided into 33 partially overlapping amplicons (
<xref ref-type="supplementary-material" rid="pone.0095969.s001">Fig. S1</xref>
and Table S4 in
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
). Amplifications were performed in a total volume of 15 ul in the presence of 1x reaction buffer, 1.5 mM MgCl2, 7.5 pmoles of each primer, 200 uM dNTPs, 1 ng of genomic DNA and 0.5 U Taq polymerase. Thermal profiles consisted of 5 min. at 96°C, 35 cycles of 30 sec. at 95°C, 30 sec. at the annealing temperatures reported in Table S4 in
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
and 1 min. at 72°C, followed by a final extension of 10 min. at 72°C. Five ul of the PCR product were then incubated 15 min. at 37°C with exosap, followed by 15 min. at 80°C. Sequencing reactions were performed with the ABI PRISM BigDye Terminator kit (v1.1, Applied Biosystems) in the presence of 1.7 pmoles of primer F (and R for amplicon 31) for 10 min. at 96°C, 25 cycles of 5 sec. at 55°C and 2 min. at 60°C. Products were separated in a PRISM ABI 3130xL Genetic Analyzer. Low quality sequence profiles were discarded upon visual inspection. High quality profiles were aligned with the Lasergene SeqMan software together with sequences representing haplogroups A to D (accession nos. GU068049, GU229278, GU295658, GU229279, GU229280, GU229281).</p>
<p>Throughout the text position numbering refers to sequence GU068049. All sequences are deposited in GenBank with Acc. n. KJ192209-KJ192236.</p>
<p>Two individuals (id. 1237 and 1300) displayed ambiguous peaks at pos. 4587 (
<xref ref-type="supplementary-material" rid="pone.0095969.s004">Fig. S4</xref>
). As an unambiguous T was observed in 6 additional individuals it was considered a real variant. However the T/C ambiguity (Y) was retained in all analyses involving id. 1237 and 1300.</p>
</sec>
<sec id="s4c">
<title>Data analysis</title>
<p>A Bayesian tree reconstruction was performed with BEAST
<xref rid="pone.0095969-Drummond1" ref-type="bibr">[33]</xref>
. We used all the 28 nearly complete sequences between positions 1 and 15,430, i.e. outside the D-loop, retaining all uncertainties present in the original data. In this analysis we included the reference sequences GU068049 and GU295658, obtained from an Hg A and an Hg B individual, respectively, but excluded the other sequences reported in ref.
<xref rid="pone.0095969-Hassanin1" ref-type="bibr">[34]</xref>
to be unexpectedly divergent or chimeric (final whole set of 30 sequences). A number of preliminary runs were used to set the priors, a suitable starting tree and to explore the possibility of estimating absolute dates for the tree nodes (see
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
). In the final analysis we used a substitution rate of 1 to obtain branch lengths in mutational units. We then used a GTR+Γ+i substitution model with strict clock and constant population size for a 10-million step run, sampled every 10,000 steps. After discarding a 50% burn-in, the results were inspected with Tracer, condensed with TreeAnnotator and visualized with TreeFIG as recommended in the user manual.</p>
<p>In order to obtain an independent interpretation of the phylogentic relationships between sequences and to precisely locate mutations on branches, we constructed a median joining network
<xref rid="pone.0095969-Bandelt1" ref-type="bibr">[35]</xref>
. In order to avoid a large number of unsupported reticulations, in this analysis we removed all positions not typed in 9 or more of the 28 individuals, plus 6 problematic positions (177, 178, 605, 4606, 9600, 11567) in which the inferred number of state changes did not fit that immediately deduced upon visual inspection (due to the additional presence of missing data in some of the sequences). The statistics file produced by the program was used to identify positions undergoing multiple mutational changes along the network. The assignment of mutations to each branch is reported in Table S2 in
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
, with branches numbered as in
<xref ref-type="supplementary-material" rid="pone.0095969.s005">Fig. S5</xref>
.</p>
<p>The assignment of D-loop sequences from external datasets to the same clades here used was based on the presence of haplotypic arrangements of variants described in
<xref ref-type="table" rid="pone-0095969-t002">Table 2</xref>
of
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
. For each of the cited datasets, the full set of sequences was aligned with a non-redundant subset of those in ref.
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
, and allelic states compared after renumbering.</p>
<p>Synonymous and non-synonymous substitutions were identified by translating the DNA sequences in MEGA
<xref rid="pone.0095969-Kumar1" ref-type="bibr">[36]</xref>
. The counts of non-synonymous and synonymous variants used in the McDonald-Kreitman tests
<xref rid="pone.0095969-McDonald1" ref-type="bibr">[19]</xref>
were obtained by aligning the DNA sequences of the 13 protein-coding genes in GU068049 with the orthologues of
<italic>Capra ibex</italic>
isolate G1253 (Acc. n. AB743814.1 to AB743826.1). From the raw counts we excluded all positions for which the difference between C.
<italic>ibex</italic>
and
<italic>C. hircus</italic>
GU068049 was attributable to the presence of a derived polymorphic allele in the latter, thus leaving only variants fixed before the divergence of Hg's C, B and A.</p>
<p>The complete list of mtDNA variants recorded at MITOMAP (
<ext-link ext-link-type="uri" xlink:href="http://www.mitomap.org/MITOMAP">http://www.mitomap.org/MITOMAP</ext-link>
) as of Aug. 30, 2013 was downloaded and analyzed in a spreadsheet. This contains data for 7932 DNA variants obtained from 18363 GenBank (GB) sequences with size greater than 14 kb. To ensure that orthologous positions at the aminoacid level were compared, alignments of goat and human protein sequences were obtained with COBALT (
<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/tools/cobalt/">http://www.ncbi.nlm.nih.gov/tools/cobalt/</ext-link>
), using the protein id's in GU68049 and in the human revised Cambridge Reference Sequence (rCRS, Acc. n. NC_012920). For RNA-coding genes, a multiple alignment of GU068049, GU295658 and the rCRS sequences was generated with Clustal
<xref rid="pone.0095969-Thompson1" ref-type="bibr">[37]</xref>
, with default settings. For each position variant in the goat, a surrounding range of 40 bp was re-aligned and visually inspected to exclude instances with multiple alignment solutions.</p>
<p>The list of all possible aminoacid substitutions in the rCRS
<xref rid="pone.0095969-Pereira2" ref-type="bibr">[22]</xref>
was downloaded and analyzed in a spreadsheet. The use of this list was possible only for positions in which one of the goat variant aminoacids matched the homologue in rCRS. The list includes the MutPred pathogenicity score, which is determined by a set of features reflecting protein structure and its dynamics, the presence of functional residues, biases of amino acid sequence, and evolutionary conservation at the substitution site and in its neighborhood. This score ranges from 0 (null pathogenicity) to 1 (maximum pathogenicity). Also included is the mtDNA selection score, defined as the ratio between the probability of observing a substitution of a given pathogenicity among all substitutions and the probability of observing it in three major clades of the human mtDNA tree. Values lower than 1 denote substitutions with likely weak or null detrimental effects.</p>
<p>Statistical tests were performed with R.</p>
</sec>
</sec>
<sec sec-type="supplementary-material" id="s5">
<title>Supporting Information</title>
<supplementary-material content-type="local-data" id="pone.0095969.s001">
<label>Figure S1</label>
<caption>
<p>
<bold>Linearized map of goat mtDNA (at bottom).</bold>
The positions of the 33 amplicons used to sequence the entire molecule are shown in the lower panel as black segments. Note that Fragment 33 spans mtDNA position 1. A representation of the coverage of the entire mitogenome in 28 individuals is given in the top panel. Lines are colored according to D-loop clade affiliation (reported on the left, after individual id.). Grayed segments indicate portions not covered by sequencing or producing poor results. Individual id's are as in
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
.</p>
<p>(JPG)</p>
</caption>
<media xlink:href="pone.0095969.s001.jpg">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0095969.s002">
<label>Figure S2</label>
<caption>
<p>
<bold>Phylogenetic tree of 28 goat mtDNAs (positions 1-15,430) plus two reference sequences, obtained with BEAST.</bold>
Individual id's and D-loop clade affiliations are reported at the tips of the tree. The posterior for each node is shown next to it, in italics. The scale bar is in mutation units. The table at right shows the assignment of domestic goat sequences to each clade (based solely on D-loop sequences, according to criteria of ref.
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
) and their geographic provenance. Weu  =  western Europe (ref.
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
); Ceu  =  central Europe (ref.
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
); Meu  =  Mediterranean Europe excluding Sardinia (ref.
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
); Sar  =  Sardinia (refs.
<xref rid="pone.0095969-Piras1" ref-type="bibr">[12]</xref>
+
<xref rid="pone.0095969-Vacca1" ref-type="bibr">[18]</xref>
); Naf  =  north Africa (refs.
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
+
<xref rid="pone.0095969-Pereira1" ref-type="bibr">[13]</xref>
); Ssaf  =  sub-saharan Africa (ref.
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
); Cas  =  central Asia (ref.
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
); Eas  =  eastern Asia (ref.
<xref rid="pone.0095969-Naderi1" ref-type="bibr">[5]</xref>
). Only one line is given for tree tips corresponding to clades A5 and A8. The corresponding assignments for the 29 haplogroup A bezoar (
<italic>C. a.)</italic>
sequences (ref.
<xref rid="pone.0095969-Naderi2" ref-type="bibr">[6]</xref>
) are also reported. 1. includes an ambiguous A7/A9 assignment; 2. includes an ambiguous A8/A9 assignment.</p>
<p>(JPG)</p>
</caption>
<media xlink:href="pone.0095969.s002.jpg">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0095969.s003">
<label>Figure S3</label>
<caption>
<p>
<bold>Ratio between the number of non-synonymous and synonymous variants in the 13 mtDNA protein-coding genes among all variants recorded in MITOMAP (black bars), variants with >40 records in GenBank (grey bars) and variants found in this study (white bars).</bold>
</p>
<p>(JPG)</p>
</caption>
<media xlink:href="pone.0095969.s003.jpg">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0095969.s004">
<label>Figure S4</label>
<caption>
<p>
<bold>Portions of the electropherograms of individuals 1237 and 1300 for amplicon 3, showing a double peak at pos. 4587.</bold>
</p>
<p>(PPTX)</p>
</caption>
<media xlink:href="pone.0095969.s004.ppt">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0095969.s005">
<label>Figure S5</label>
<caption>
<p>
<bold>Same network as in </bold>
<xref ref-type="fig" rid="pone-0095969-g001">
<bold>Fig. 1</bold>
</xref>
<bold> with numbered branches.</bold>
The assignment of the mutations listed in Table S2 in
<xref ref-type="supplementary-material" rid="pone.0095969.s006">File S1</xref>
to each branch follows this numbering system.</p>
<p>(JPG)</p>
</caption>
<media xlink:href="pone.0095969.s005.jpg">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0095969.s006">
<label>File S1</label>
<caption>
<p>
<bold>Supporting tables and text.</bold>
</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0095969.s006.docx">
<caption>
<p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back>
<ack>
<p>We gratefully acknowledge all farmers of Ogliastra, the ARAS (Regional Association of Sardinian Farmers), the local authorities and the municipalities of Ogliastra, who took part in the present study, giving access to the animals and their life history data.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="pone.0095969-CluttonBrock1">
<label>1</label>
<mixed-citation publication-type="other">Clutton-Brock J (1999) A natural history of domesticated mammals. 2nd ed. Cambridge (UK) :Cambridge University Press.</mixed-citation>
</ref>
<ref id="pone.0095969-Diamond1">
<label>2</label>
<mixed-citation publication-type="journal">
<name>
<surname>Diamond</surname>
<given-names>J</given-names>
</name>
(
<year>2002</year>
)
<article-title>Evolution, consequences and future of plant and animal domestication</article-title>
.
<source>Nature</source>
<volume>418</volume>
:
<fpage>700</fpage>
<lpage>707</lpage>
<pub-id pub-id-type="pmid">12167878</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Zeder1">
<label>3</label>
<mixed-citation publication-type="journal">
<name>
<surname>Zeder</surname>
<given-names>MA</given-names>
</name>
(
<year>2008</year>
)
<article-title>Domestication and early agriculture in the Mediterranean Basin: Origins, diffusion, and impact</article-title>
.
<source>Proc Natl Acad Sci USA</source>
<volume>105</volume>
:
<fpage>11597</fpage>
<lpage>11604</lpage>
<pub-id pub-id-type="pmid">18697943</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Luikart1">
<label>4</label>
<mixed-citation publication-type="journal">
<name>
<surname>Luikart</surname>
<given-names>G</given-names>
</name>
,
<name>
<surname>Gielly</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Excoffier</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Vigne</surname>
<given-names>J-D</given-names>
</name>
,
<name>
<surname>Bouvet</surname>
<given-names>J</given-names>
</name>
,
<etal>et al</etal>
(
<year>2001</year>
)
<article-title>Multiple maternal origins and weak phylogeographic structure in domestic goats</article-title>
.
<source>Proc Natl Acad Sci USA</source>
<volume>98</volume>
:
<fpage>5927</fpage>
<lpage>5932</lpage>
<pub-id pub-id-type="pmid">11344314</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Naderi1">
<label>5</label>
<mixed-citation publication-type="journal">
<name>
<surname>Naderi</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Rezaei</surname>
<given-names>H-R</given-names>
</name>
,
<name>
<surname>Taberlet</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Zundel</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Rafat</surname>
<given-names>S-A</given-names>
</name>
,
<etal>et al</etal>
(
<year>2007</year>
)
<article-title>Large-scale mitochondrial DNA analysis of the domestic goat reveals six haplogroups with high diversity</article-title>
.
<source>PLoS ONE</source>
<volume>2</volume>
:
<fpage>e1012</fpage>
<pub-id pub-id-type="pmid">17925860</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Naderi2">
<label>6</label>
<mixed-citation publication-type="journal">
<name>
<surname>Naderi</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Rezaei</surname>
<given-names>H-R</given-names>
</name>
,
<name>
<surname>Pompanon</surname>
<given-names>F</given-names>
</name>
,
<name>
<surname>Blum</surname>
<given-names>MGB</given-names>
</name>
,
<name>
<surname>Negrini</surname>
<given-names>R</given-names>
</name>
,
<etal>et al</etal>
(
<year>2008</year>
)
<article-title>The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals</article-title>
.
<source>Proc Natl Acad Sci USA</source>
<volume>105</volume>
:
<fpage>17659</fpage>
<lpage>17664</lpage>
<pub-id pub-id-type="pmid">19004765</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Hughes1">
<label>7</label>
<mixed-citation publication-type="journal">
<name>
<surname>Hughes</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Fernàndez</surname>
<given-names>H</given-names>
</name>
,
<name>
<surname>Cucchi</surname>
<given-names>T</given-names>
</name>
,
<name>
<surname>Duffraisse</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Casabianca</surname>
<given-names>F</given-names>
</name>
,
<etal>et al</etal>
(
<year>2012</year>
)
<article-title>A dig into the past mitochondrial diversity of Corsican goats reveals the influence of secular herding practices</article-title>
.
<source>PLoS ONE</source>
<volume>7</volume>
:
<fpage>e30272</fpage>
<pub-id pub-id-type="pmid">22299033</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Taberlet1">
<label>8</label>
<mixed-citation publication-type="journal">
<name>
<surname>Taberlet</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Valentini</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Rezaei</surname>
<given-names>HR</given-names>
</name>
,
<name>
<surname>Naderi</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Pompanon</surname>
<given-names>F</given-names>
</name>
,
<etal>et al</etal>
(
<year>2008</year>
)
<article-title>Are cattle, sheep, and goats endangered species?</article-title>
<source>Mol Ecol</source>
<volume>17</volume>
:
<fpage>275</fpage>
<lpage>284</lpage>
<pub-id pub-id-type="pmid">17927711</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Taberlet2">
<label>9</label>
<mixed-citation publication-type="journal">
<name>
<surname>Taberlet</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Coissac</surname>
<given-names>E</given-names>
</name>
,
<name>
<surname>Pansu</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Pompanon</surname>
<given-names>F</given-names>
</name>
(
<year>2011</year>
)
<article-title>Conservation genetics of cattle, sheep, and goats</article-title>
.
<source>Comptes Rendus Biologies</source>
<volume>334</volume>
:
<fpage>247</fpage>
<lpage>254</lpage>
<pub-id pub-id-type="pmid">21377620</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Achilli1">
<label>10</label>
<mixed-citation publication-type="journal">
<name>
<surname>Achilli</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Bonfiglio</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Olivieri</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Malusà</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Pala</surname>
<given-names>M</given-names>
</name>
,
<etal>et al</etal>
(
<year>2009</year>
)
<article-title>The multifaceted origin of taurine cattle reflected by the mitochondrial genome</article-title>
.
<source>PLoS ONE</source>
<volume>4</volume>
:
<fpage>e5753</fpage>
<pub-id pub-id-type="pmid">19484124</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Lancioni1">
<label>11</label>
<mixed-citation publication-type="journal">
<name>
<surname>Lancioni</surname>
<given-names>H</given-names>
</name>
,
<name>
<surname>Di Lorenzo</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Ceccobelli</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Perego</surname>
<given-names>UA</given-names>
</name>
,
<name>
<surname>Miglio</surname>
<given-names>A</given-names>
</name>
,
<etal>et al</etal>
(
<year>2013</year>
)
<article-title>Phylogenetic relationships of three Italian Merino-derived sheep breeds evaluated through a complete mitogenome analysis</article-title>
.
<source>PLoS ONE</source>
<volume>8</volume>
:
<fpage>e73712</fpage>
<pub-id pub-id-type="pmid">24040036</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Piras1">
<label>12</label>
<mixed-citation publication-type="journal">
<name>
<surname>Piras</surname>
<given-names>D</given-names>
</name>
,
<name>
<surname>Doro</surname>
<given-names>MG</given-names>
</name>
,
<name>
<surname>Casu</surname>
<given-names>G</given-names>
</name>
,
<name>
<surname>Melis</surname>
<given-names>MP</given-names>
</name>
,
<name>
<surname>Vaccargiu</surname>
<given-names>S</given-names>
</name>
,
<etal>et al</etal>
(
<year>2012</year>
)
<article-title>Haplotype affinities resolve a major component of goat (
<italic>Capra hircus</italic>
) mtDNA D-loop diversity and reveal specific features of the Sardinian stock</article-title>
.
<source>PLoS ONE</source>
<volume>7</volume>
:
<fpage>e30785</fpage>
<pub-id pub-id-type="pmid">22363488</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Pereira1">
<label>13</label>
<mixed-citation publication-type="journal">
<name>
<surname>Pereira</surname>
<given-names>F</given-names>
</name>
,
<name>
<surname>Queiròs</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Gusmào</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Nijman</surname>
<given-names>IJ</given-names>
</name>
,
<name>
<surname>Cuppen</surname>
<given-names>E</given-names>
</name>
,
<etal>et al</etal>
(
<year>2009</year>
)
<article-title>Tracing the history of goat pastoralism: new clues from mitochondrial and Y chromosome DNA in North Africa</article-title>
.
<source>Mol Biol Evol</source>
<volume>26</volume>
:
<fpage>2765</fpage>
<lpage>2773</lpage>
<pub-id pub-id-type="pmid">19729424</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Frankham1">
<label>14</label>
<mixed-citation publication-type="other">Frankham R, Ballou JD, Briscoe DA (2010) Introduction to Conservation Genetics. 2nd ed. Cambridge, U.K.: Cambridge University Press.</mixed-citation>
</ref>
<ref id="pone.0095969-Hughes2">
<label>15</label>
<mixed-citation publication-type="journal">
<name>
<surname>Hughes</surname>
<given-names>AL</given-names>
</name>
(
<year>2013</year>
)
<article-title>Accumulation of slightly deleterious mutations in the mitochondrial genome: A hallmark of animal domestication</article-title>
.
<source>Gene</source>
<volume>515</volume>
:
<fpage>28</fpage>
<lpage>33</lpage>
<pub-id pub-id-type="pmid">23237775</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Bjrnerfeldt1">
<label>16</label>
<mixed-citation publication-type="journal">
<name>
<surname>Björnerfeldt</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Webster</surname>
<given-names>MT</given-names>
</name>
,
<name>
<surname>Vilà</surname>
<given-names>C</given-names>
</name>
(
<year>2006</year>
)
<article-title>Relaxation of selective constraint on dog mitochondrial DNA following domestication</article-title>
.
<source>Genome Res</source>
<volume>16</volume>
:
<fpage>990</fpage>
<lpage>994</lpage>
<pub-id pub-id-type="pmid">16809672</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Nomura1">
<label>17</label>
<mixed-citation publication-type="journal">
<name>
<surname>Nomura</surname>
<given-names>K</given-names>
</name>
,
<name>
<surname>Yonezawa</surname>
<given-names>T</given-names>
</name>
,
<name>
<surname>Mano</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Kawakami</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Shedlock</surname>
<given-names>AM</given-names>
</name>
,
<etal>et al</etal>
(
<year>2013</year>
)
<article-title>Domestication process of the goat revealed by an analysis of the nearly complete mitochondrial protein-encoding genes</article-title>
.
<source>PLoS ONE</source>
<volume>8</volume>
:
<fpage>e67775</fpage>
<pub-id pub-id-type="pmid">23936295</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Vacca1">
<label>18</label>
<mixed-citation publication-type="journal">
<name>
<surname>Vacca</surname>
<given-names>GM</given-names>
</name>
,
<name>
<surname>Daga</surname>
<given-names>C</given-names>
</name>
,
<name>
<surname>Pazzola</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Carcangiu</surname>
<given-names>V</given-names>
</name>
,
<name>
<surname>Dettori</surname>
<given-names>ML</given-names>
</name>
,
<etal>et al</etal>
(
<year>2010</year>
)
<article-title>D-loop sequence mitochondrial DNA variability of Sarda goat and other goat breeds and populations reared in the Mediterranean area</article-title>
.
<source>J An Breed Genet</source>
<volume>127</volume>
:
<fpage>352</fpage>
<lpage>360</lpage>
</mixed-citation>
</ref>
<ref id="pone.0095969-McDonald1">
<label>19</label>
<mixed-citation publication-type="journal">
<name>
<surname>McDonald</surname>
<given-names>JH</given-names>
</name>
,
<name>
<surname>Kreitman</surname>
<given-names>M</given-names>
</name>
(
<year>1991</year>
)
<article-title>Adaptive protein evolution at the Adh locus in Drosophila</article-title>
.
<source>Nature</source>
<volume>351</volume>
:
<fpage>652</fpage>
<lpage>654</lpage>
<pub-id pub-id-type="pmid">1904993</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Zeder2">
<label>20</label>
<mixed-citation publication-type="journal">
<name>
<surname>Zeder</surname>
<given-names>MA</given-names>
</name>
(
<year>2001</year>
)
<article-title>A metrical analysis of a collection of modern goats (
<italic>Capra hircus aegargus</italic>
and
<italic>C. h. hircus</italic>
) from Iran and Iraq: Implications for the study of Caprine domestication</article-title>
.
<source>J Archaeol Sci</source>
<volume>28</volume>
:
<fpage>61</fpage>
<lpage>79</lpage>
</mixed-citation>
</ref>
<ref id="pone.0095969-Sakuta1">
<label>21</label>
<mixed-citation publication-type="journal">
<name>
<surname>Sakuta</surname>
<given-names>R</given-names>
</name>
,
<name>
<surname>Goto</surname>
<given-names>Y</given-names>
</name>
,
<name>
<surname>Nonaka</surname>
<given-names>I</given-names>
</name>
,
<name>
<surname>Horai</surname>
<given-names>S</given-names>
</name>
(
<year>1993</year>
)
<article-title>An A-to-G transition at nucleotide pair 11084 in the ND4 gene may be an mtDNA polymorphism</article-title>
.
<source>Am J Hum Genet</source>
<volume>53</volume>
:
<fpage>964</fpage>
<lpage>965</lpage>
<pub-id pub-id-type="pmid">8213827</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Pereira2">
<label>22</label>
<mixed-citation publication-type="journal">
<name>
<surname>Pereira</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Soares</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Radivojac</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Li</surname>
<given-names>B</given-names>
</name>
,
<name>
<surname>Samuels</surname>
<given-names>DC</given-names>
</name>
(
<year>2011</year>
)
<article-title>Comparing phylogeny and the predicted pathogenicity of protein variations reveals equal purifying selection across the global human mtDNA diversity</article-title>
.
<source>Am J Hum Genet</source>
<volume>88</volume>
:
<fpage>433</fpage>
<lpage>439</lpage>
<pub-id pub-id-type="pmid">21457906</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Fernndez1">
<label>23</label>
<mixed-citation publication-type="journal">
<name>
<surname>Fernàndez</surname>
<given-names>H</given-names>
</name>
,
<name>
<surname>Hughes</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Vigne</surname>
<given-names>J-D</given-names>
</name>
,
<name>
<surname>Helmer</surname>
<given-names>D</given-names>
</name>
,
<name>
<surname>Hodgins</surname>
<given-names>G</given-names>
</name>
,
<etal>et al</etal>
(
<year>2006</year>
)
<article-title>Divergent mtDNA lineages of goats in an Early Neolithic site, far from the initial domestication areas</article-title>
.
<source>Proc Natl Acad Sci USA</source>
<volume>103</volume>
:
<fpage>15375</fpage>
<lpage>15379</lpage>
<pub-id pub-id-type="pmid">17030824</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Amills1">
<label>24</label>
<mixed-citation publication-type="journal">
<name>
<surname>Amills</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Ramírez</surname>
<given-names>O</given-names>
</name>
,
<name>
<surname>Tomàs</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Badaoui</surname>
<given-names>B</given-names>
</name>
,
<name>
<surname>Marmi</surname>
<given-names>J</given-names>
</name>
,
<etal>et al</etal>
(
<year>2009</year>
)
<article-title>Mitochondrial DNA diversity and origins of South and Central American goats</article-title>
.
<source>Animal Genet</source>
<volume>40</volume>
:
<fpage>315</fpage>
<lpage>322</lpage>
<pub-id pub-id-type="pmid">19292711</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Nabholz1">
<label>25</label>
<mixed-citation publication-type="journal">
<name>
<surname>Nabholz</surname>
<given-names>B</given-names>
</name>
,
<name>
<surname>Glèmin</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Galtier</surname>
<given-names>N</given-names>
</name>
(
<year>2008</year>
)
<article-title>Strong variations of mitochondrial mutation rate across mammals-the longevity hypothesis</article-title>
.
<source>Mol Biol Evol</source>
<volume>25</volume>
:
<fpage>120</fpage>
<lpage>130</lpage>
<pub-id pub-id-type="pmid">17998254</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Tudge1">
<label>26</label>
<mixed-citation publication-type="other">Tudge C (1998) Neanderthals, bandits and farmers: how agriculture really began. New Haven, CT: Yale Univ Press.</mixed-citation>
</ref>
<ref id="pone.0095969-Shipman1">
<label>27</label>
<mixed-citation publication-type="other">Shipman P (2011) The animal connection. New York: Norton.</mixed-citation>
</ref>
<ref id="pone.0095969-Soares1">
<label>28</label>
<mixed-citation publication-type="journal">
<name>
<surname>Soares</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Abrantes</surname>
<given-names>D</given-names>
</name>
,
<name>
<surname>Rito</surname>
<given-names>T</given-names>
</name>
,
<name>
<surname>Thomson</surname>
<given-names>N</given-names>
</name>
,
<name>
<surname>Radivojac</surname>
<given-names>P</given-names>
</name>
,
<etal>et al</etal>
(
<year>2013</year>
)
<article-title>Evaluating purifying selection in the mitochondrial DNA of various mammalian species</article-title>
.
<source>PLoS ONE</source>
<volume>8</volume>
:
<fpage>e58993</fpage>
<pub-id pub-id-type="pmid">23533597</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Woolley1">
<label>29</label>
<mixed-citation publication-type="journal">
<name>
<surname>Woolley</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Johnson</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Smith</surname>
<given-names>MJ</given-names>
</name>
,
<name>
<surname>Crandall</surname>
<given-names>KA</given-names>
</name>
,
<name>
<surname>McClellan</surname>
<given-names>DA</given-names>
</name>
(
<year>2003</year>
)
<article-title>TreeSAAP: Selection on Amino Acid Properties using phylogenetic trees</article-title>
.
<source>Bioinformatics</source>
<volume>19</volume>
:
<fpage>671</fpage>
<lpage>672</lpage>
<pub-id pub-id-type="pmid">12651734</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Ballard1">
<label>30</label>
<mixed-citation publication-type="journal">
<name>
<surname>Ballard</surname>
<given-names>JWO</given-names>
</name>
,
<name>
<surname>Melvin</surname>
<given-names>RG</given-names>
</name>
(
<year>2010</year>
)
<article-title>Linking the mitochondrial genotype to the organismal phenotype</article-title>
.
<source>Mol Ecol</source>
<volume>19</volume>
:
<fpage>1523</fpage>
<lpage>1539</lpage>
<pub-id pub-id-type="pmid">20345689</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Kang1">
<label>31</label>
<mixed-citation publication-type="journal">
<name>
<surname>Kang</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Zheng</surname>
<given-names>H-X</given-names>
</name>
,
<name>
<surname>Chen</surname>
<given-names>F</given-names>
</name>
,
<name>
<surname>Yan</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Liu</surname>
<given-names>K</given-names>
</name>
,
<etal>et al</etal>
(
<year>2013</year>
)
<article-title>mtDNA lineage expansions in Sherpa population suggest adaptive evolution in Tibetan highlands</article-title>
.
<source>Mol Biol Evol</source>
<volume>30</volume>
:
<fpage>2579</fpage>
<lpage>2587</lpage>
<pub-id pub-id-type="pmid">24002810</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Brandano1">
<label>32</label>
<mixed-citation publication-type="journal">
<name>
<surname>Brandano</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Piras</surname>
<given-names>B</given-names>
</name>
(
<year>1978</year>
)
<article-title>La capra Sarda I. I caratteri morfologici (in Italian)</article-title>
.
<source>Annali della Facoltà di Agraria di Sassari</source>
<volume>26</volume>
:
<fpage>232</fpage>
<lpage>265</lpage>
</mixed-citation>
</ref>
<ref id="pone.0095969-Drummond1">
<label>33</label>
<mixed-citation publication-type="journal">
<name>
<surname>Drummond</surname>
<given-names>AJ</given-names>
</name>
,
<name>
<surname>Rambaut</surname>
<given-names>A</given-names>
</name>
(
<year>2007</year>
)
<article-title>BEAST: Bayesian evolutionary analysis by sampling trees</article-title>
.
<source>BMC Evol Biol</source>
<volume>7</volume>
:
<fpage>214</fpage>
<pub-id pub-id-type="pmid">17996036</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Hassanin1">
<label>34</label>
<mixed-citation publication-type="journal">
<name>
<surname>Hassanin</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Bonillo</surname>
<given-names>C</given-names>
</name>
,
<name>
<surname>Nguyen</surname>
<given-names>BX</given-names>
</name>
,
<name>
<surname>Cruaud</surname>
<given-names>C</given-names>
</name>
(
<year>2010</year>
)
<article-title>Comparisons between mitochondrial genomes of domestic goat (
<italic>Capra hircus</italic>
) reveal the presence of numts and multiple sequencing errors</article-title>
.
<source>Mitochondrial DNA</source>
<volume>21</volume>
:
<fpage>68</fpage>
<lpage>76</lpage>
<pub-id pub-id-type="pmid">20540682</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Bandelt1">
<label>35</label>
<mixed-citation publication-type="journal">
<name>
<surname>Bandelt</surname>
<given-names>HJ</given-names>
</name>
,
<name>
<surname>Forster</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Rohl</surname>
<given-names>A</given-names>
</name>
(
<year>1999</year>
)
<article-title>Median-joining networks for inferring intraspecific phylogenies</article-title>
.
<source>Mol Biol Evol</source>
<volume>16</volume>
:
<fpage>37</fpage>
<lpage>48</lpage>
<pub-id pub-id-type="pmid">10331250</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Kumar1">
<label>36</label>
<mixed-citation publication-type="journal">
<name>
<surname>Kumar</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Nei</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Dudley</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Tamura</surname>
<given-names>K</given-names>
</name>
(
<year>2008</year>
)
<article-title>MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences</article-title>
.
<source>Brief Bioinform</source>
<volume>9</volume>
:
<fpage>299</fpage>
<lpage>306</lpage>
<pub-id pub-id-type="pmid">18417537</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0095969-Thompson1">
<label>37</label>
<mixed-citation publication-type="journal">
<name>
<surname>Thompson</surname>
<given-names>JD</given-names>
</name>
,
<name>
<surname>Gibson</surname>
<given-names>T</given-names>
</name>
,
<name>
<surname>Plewniak</surname>
<given-names>F</given-names>
</name>
,
<name>
<surname>Jeanmougin</surname>
<given-names>F</given-names>
</name>
,
<name>
<surname>Higgins</surname>
<given-names>DG</given-names>
</name>
(
<year>1997</year>
)
<article-title>The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools</article-title>
.
<source>Nucleic Acids Res</source>
<volume>24</volume>
:
<fpage>4876</fpage>
<lpage>4882</lpage>
<pub-id pub-id-type="pmid">9396791</pub-id>
</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>

Pour manipuler ce document sous Unix (Dilib)

EXPLOR_STEP=$WICRI_ROOT/Wicri/Terre/explor/CobaltMaghrebV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000458 | SxmlIndent | more

Ou

HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000458 | SxmlIndent | more

Pour mettre un lien sur cette page dans le réseau Wicri

{{Explor lien
   |wiki=    Wicri/Terre
   |area=    CobaltMaghrebV1
   |flux=    Pmc
   |étape=   Corpus
   |type=    RBID
   |clé=     PMC:3999278
   |texte=   Phylogeny and Patterns of Diversity of Goat mtDNA Haplogroup A Revealed by Resequencing Complete Mitogenomes
}}

Pour générer des pages wiki

HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i   -Sk "pubmed:24763315" \
       | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd   \
       | NlmPubMed2Wicri -a CobaltMaghrebV1 

Wicri

This area was generated with Dilib version V0.6.32.
Data generation: Tue Nov 14 12:56:51 2017. Site generation: Mon Feb 12 07:59:49 2024