Whole mitochondrial genomes unveil the impact of domestication on goat matrilineal variability

Mitochondrial Genomes of Korean Native Black Goats Reveal Shared Phylogeographic Patterns and Demographic History

This study explores the phylogeny of Korean native black goats through analysis of their complete mitochondrial DNA. The National Institute of Animal Science has gathered genetic material on purebred goats from isolated regions such as Tongyeong, Dangjin, and Jangsu, and is actively breeding them on a national level. These populations, however, are small and exhibit high inbreeding rates, highlighting the urgent need to preserve genetic diversity. The haplotype diversity within this native group is 0.659, with 39 haplotypes identified. By contrast, including international breeds in the analysis increases the overall haplotype diversity to 0.925 with 203 haplotypes identified, highlighting the limited genetic diversity among native black goats. For phylogenetic assessment, a neighbor-joining tree and median-joining network were constructed using identified haplogroups (A, B, C, D, G, and F) from prior studies. The results pinpoint the native black goats as closely related to, but distinct from, Haplogroup A with a bootstrap value of 98, establishing them as a separate clade (A'). This supports the notion of a shared ancestry with various global populations. This research provides essential data on the origins and evolutionary history of Korean native black goats, supporting conservation and breeding efforts aimed at enhancing genetic diversity.

Extensive Population Structure Highlights an Apparent Paradox of Stasis in the Impala (Aepyceros melampus)

Impalas are unusual among bovids because they have remained morphologically similar over millions of years-a phenomenon referred to as evolutionary stasis. Here, we sequenced 119 whole genomes from the two extant subspecies of impala, the common (Aepyceros melampus melampus) and black-faced (A. m. petersi) impala. We investigated the evolutionary forces working within the species to explore how they might be associated with its evolutionary stasis as a taxon. Despite being one of the most abundant bovid species, we found low genetic diversity overall, and a phylogeographic signal of spatial expansion from southern to eastern Africa. Contrary to expectations under a scenario of evolutionary stasis, we found pronounced genetic structure between and within the two subspecies with indications of ancient, but not recent, gene flow. Black-faced impala and eastern African common impala populations had more runs of homozygosity than common impala in southern Africa, and, using a proxy for genetic load, we found that natural selection is working less efficiently in these populations compared to the southern African populations. Together with the fossil record, our results are consistent with a fixed-optimum model of evolutionary stasis, in which impalas in the southern African core of the range are able to stay near their evolutionary fitness optimum as a generalist ecotone species, whereas eastern African impalas may struggle to do so due to the effects of genetic drift and reduced adaptation to the local habitat, leading to recurrent local extinction in eastern Africa and re-colonisation from the South.

Mitochondrial DNA control-region and coding-region data highlight geographically structured diversity and post-domestication population dynamics in worldwide donkeys

Donkeys (Equus asinus) have been used extensively in agriculture and transportations since their domestication, ca. 5000-7000 years ago, but the increased mechanization of the last century has largely spoiled their role as burden animals, particularly in developed countries. Consequently, donkey breeds and population sizes have been declining for decades, and the diversity contributed by autochthonous gene pools has been eroded. Here, we examined coding-region data extracted from 164 complete mitogenomes and 1392 donkey mitochondrial DNA (mtDNA) control-region sequences to (i) assess worldwide diversity, (ii) evaluate geographical patterns of variation, and (iii) provide a new nomenclature of mtDNA haplogroups. The topology of the Maximum Parsimony tree confirmed the two previously identified major clades, i.e. Clades 1 and 2, but also highlighted the occurrence of a deep-diverging lineage within Clade 2 that left a marginal trace in modern donkeys. Thanks to the identification of stable and highly diagnostic coding-region mutational motifs, the two lineages were renamed as haplogroup A and haplogroup B, respectively, to harmonize clade nomenclature with the standard currently adopted for other livestock species. Control-region diversity and population expansion metrics varied considerably between geographical areas but confirmed North-eastern Africa as the likely domestication center. The patterns of geographical distribution of variation analyzed through phylogenetic networks and AMOVA confirmed the co-occurrence of both haplogroups in all sampled populations, while differences at the regional level point to the joint effects of demography, past human migrations and trade following the spread of donkeys out of the domestication center. Despite the strong decline that donkey populations have undergone for decades in many areas of the world, the sizeable mtDNA variability we scored, and the possible identification of a new early radiating lineage further stress the need for an extensive and large-scale characterization of donkey nuclear genome diversity to identify hotspots of variation and aid the conservation of local breeds worldwide.

Ancient mitogenomes from Pre-Pottery Neolithic Central Anatolia and the effects of a Late Neolithic bottleneck in sheep (Ovis aries)

Occupied between ~10,300 and 9300 years ago, the Pre-Pottery Neolithic site of Aşıklı Höyük in Central Anatolia went through early phases of sheep domestication. Analysis of 629 mitochondrial genomes from this and numerous sites in Anatolia, southwest Asia, Europe, and Africa produced a phylogenetic tree with excessive coalescences (nodes) around the Neolithic, a potential signature of a domestication bottleneck. This is consistent with archeological evidence of sheep management at Aşıklı Höyük which transitioned from residential stabling to open pasturing over a millennium of site occupation. However, unexpectedly, we detected high genetic diversity throughout Aşıklı Höyük's occupation rather than a bottleneck. Instead, we detected a tenfold demographic bottleneck later in the Neolithic, which caused the fixation of mitochondrial haplogroup B in southwestern Anatolia. The mitochondrial genetic makeup that emerged was carried from the core region of early Neolithic sheep management into Europe and dominates the matrilineal diversity of both its ancient and the billion-strong modern sheep populations.

Mitochondrial DNA D-Loop Polymorphisms among the Galla Goats Reveals Multiple Maternal Origins with Implication on the Functional Diversity of the HSP70 Gene

Despite much attention given to the history of goat evolution in Kenya, information on the origin, demographic history, dispersal route, and genetic diversity of Galla goats remains unclear. Here, we examined the genetic background, diversity, demographic history, and population genetic variation of Galla goats using mtDNA D-loop and HSP70 single-nucleotide polymorphism markers. The results revealed 90 segregating sites and 68 haplotypes in a 600-bp mtDNA D-loop sequence. The overall mean mitochondrial haplotype diversity was 0.993. The haplotype diversities ranged between 0.8939 ± 0.0777 and 1.0000 ± 0.0221 in all populations supporting high genetic diversity. Mitochondrial phylogenetic analysis revealed three Galla goat haplogroups (A, G, and D), supporting multiple maternal ancestries, of which haplogroup A was the most predominant. Analysis of molecular variance (AMOVA) showed considerable variation within populations at 94.39%, evidence of high genetic diversity. Bimodal mismatch distribution patterns were observed while most populations recorded negative results for Tajima and Fu's Fs neutrality tests supporting population expansion. Genetic variation among populations was also confirmed using HSP70 gene fragment sequences, where six polymorphic sites which defined 21 haplotypes were discovered. Analysis of molecular variance revealed a significant FST index value of 0.134 and a high FIS index value of 0.746, an indication of inbreeding. This information will pave the way for conservation strategies and informed breeding to improve Galla or other goat breeds for climate-smart agriculture.

The Mitogenome Relationships and Phylogeography of Barn Swallows (Hirundo rustica)

The barn swallow (Hirundo rustica) poses a number of fascinating scientific questions, including the taxonomic status of postulated subspecies. Here we obtained and assessed the sequence variation of 411 complete mitogenomes, mainly from the European H. r. rustica, but other subspecies as well. In almost every case, we observed subspecies-specific haplogroups, which we employed together with estimated radiation times to postulate a model for the geographical and temporal worldwide spread of the species. The female barn swallow carrying the Hirundo rustica ancestral mitogenome left Africa (or its vicinity) around 280 thousand years ago (kya), and her descendants expanded first into Eurasia and then, at least 51 kya, into the Americas, from where a relatively recent (< 20 kya) back migration to Asia took place. The exception to the haplogroup subspecies specificity is represented by the sedentary Levantine H. r. transitiva that extensively shares haplogroup A with the migratory European H. r. rustica and, to a lesser extent, haplogroup B with the Egyptian H. r. savignii. Our data indicate that rustica and transitiva most likely derive from a sedentary Levantine population source that split at the end of the Younger Dryas (11.7 kya). Since then, however, transitiva received genetic inputs from and admixed with both the closely related rustica and the adjacent savignii. Demographic analyses confirm this species' strong link with climate fluctuations and human activities making it an excellent indicator for monitoring and assessing the impact of current global changes on wildlife.

Analysis of World-Scale Mitochondrial DNA Reveals the Origin and Migration Route of East Asia Goats

Despite much attention on the history of goat evolution, information on origin, demographic history, and expansion route remains controversial. To address these questions, we collected 4,189 published goat DNA sequences including 1,228 sequences from 57 breeds in China and 2,961 sequences including 193 goat breeds from 71 other countries and carried out an integrated analysis. We found goat breeds from South China had the highest genetic diversity of lineage B, and subclades B2 only were found in Southwest China, suggesting that lineage B (particularly, subclade B2) probably originated from Southwest China and its surrounding areas. In addition, in this study, we found that lineage A from South China also presented higher genetic diversity and earlier expansion time (10, 606 years ago), even earlier than breeds from the Middle East. Hence, we speculated that South China and surrounding areas were the origin of lineage B and also the transportation hub for lineage A spreading to North China and Southwest Asia. Furthermore, according to the analysis of correlation between genetic differentiation value λ1 and λ2 and geographical distance, we further confirmed two phases of migration in goat breeds of North China. These results will contribute to a better understanding of the origin and migration history of domestic goat.

Mitochondrial DNA D-loop sequence analysis reveals high variation and multiple maternal origins of indigenous Tanzanian goat populations

The Small East African (SEA) goat are widely distributed in different agro-ecological zones of Tanzania. We report the genetic diversity, maternal origin, and phylogenetic relationship among the 12 Tanzanian indigenous goat populations, namely Fipa, Songwe, Tanga, Pwani, Iringa, Newala, Lindi, Gogo, Pare, Maasai, Sukuma, and Ujiji, based on the mitochondrial DNA (mtDNA) D-loop. High haplotype (H d = 0.9619-0.9945) and nucleotide (π = 0.0120-0.0162) diversities were observed from a total of 389 haplotypes. The majority of the haplotypes (n = 334) belonged to Haplogroup A which was consistent with the global scenario on the genetic pattern of maternal origin of all goat breeds in the world. Haplogroup G comprised of 45 haplotypes drawn from all populations except the Ujiji goat population while Haplogroup B with 10 haplotypes was dominated by Ujiji goats (41%). Tanzanian goats shared four haplotypes with the Kenyan goats and two with goats from South Africa, Namibia, and Mozambique. There was no sharing of haplotypes observed between individuals from Tanzanian goat populations with individuals from North or West Africa. The indigenous goats in Tanzania have high genetic diversity defined by 389 haplotypes and multiple maternal origins of haplogroup A, B, and G. There is a lot of intermixing and high genetic variation within populations which represent an abundant resource for selective breeding in the different agro-ecological regions of the country.

SNP-Based Genotyping Provides Insight Into the West Asian Origin of Russian Local Goats

Specific local environmental and sociocultural conditions have led to the creation of various goat populations in Russia. National goat diversity includes breeds that have been selected for down and mohair production traits as well as versatile local breeds for which pastoralism is the main management system. Effective preservation and breeding programs for local goat breeds are missing due to the lack of DNA-based data. In this work, we analyzed the genetic diversity and population structure of Russian local goats, including Altai Mountain, Altai White Downy, Dagestan Downy, Dagestan Local, Karachaev, Orenburg, and Soviet Mohair goats, which were genotyped with the Illumina Goat SNP50 BeadChip. In addition, we addressed genetic relationships between local and global goat populations obtained from the AdaptMap project. Russian goats showed a high level of genetic diversity. Although a decrease in historical effective population sizes was revealed, the recent effective population sizes estimated for three generations ago were larger than 100 in all studied populations. The mean runs of homozygosity (ROH) lengths ranged from 79.42 to 183.94 Mb, and the average ROH number varied from 18 to 41. Short ROH segments (<2 Mb) were predominant in all breeds, while the longest ROH class (>16 Mb) was the least frequent. Principal component analysis, Neighbor-Net graph, and Admixture clustering revealed several patterns in Russian local goats. First, a separation of the Karachaev breed from other populations was observed. Moreover, genetic connections between the Orenburg and Altai Mountain breeds were suggested and the Dagestan breeds were found to be admixed with the Soviet Mohair breed. Neighbor-Net analysis and clustering of local and global breeds demonstrated the close genetic relations between Russian local and Turkish breeds that probably resulted from past admixture events through postdomestication routes. Our findings contribute to the understanding of the genetic relationships of goats originating in West Asia and Eurasia and may be used to design breeding programs for local goats to ensure their effective conservation and proper management.

Herded and hunted goat genomes from the dawn of domestication in the Zagros Mountains

The Aceramic Neolithic (∼9600 to 7000 cal BC) period in the Zagros Mountains, western Iran, provides some of the earliest archaeological evidence of goat (Capra hircus) management and husbandry by circa 8200 cal BC, with detectable morphological change appearing ∼1,000 y later. To examine the genomic imprint of initial management and its implications for the goat domestication process, we analyzed 14 novel nuclear genomes (mean coverage 1.13X) and 32 mitochondrial (mtDNA) genomes (mean coverage 143X) from two such sites, Ganj Dareh and Tepe Abdul Hosein. These genomes show two distinct clusters: those with domestic affinity and a minority group with stronger wild affinity, indicating that managed goats were genetically distinct from wild goats at this early horizon. This genetic duality, the presence of long runs of homozygosity, shared ancestry with later Neolithic populations, a sex bias in archaeozoological remains, and demographic profiles from across all layers of Ganj Dareh support management of genetically domestic goat by circa 8200 cal BC, and represent the oldest to-this-date reported livestock genomes. In these sites a combination of high autosomal and mtDNA diversity, contrasting limited Y chromosomal lineage diversity, an absence of reported selection signatures for pigmentation, and the wild morphology of bone remains illustrates domestication as an extended process lacking a strong initial bottleneck, beginning with spatial control, demographic manipulation via biased male culling, captive breeding, and subsequently phenotypic and genomic selection.

The assembly of caprine Y chromosome sequence reveals a unique paternal phylogenetic pattern and improves our understanding of the origin of domestic goat

The mammalian Y chromosome offers a unique perspective on the male reproduction and paternal evolutionary histories. However, further understanding of the Y chromosome biology for most mammals is hindered by the lack of a Y chromosome assembly. This study presents an integrated in silico strategy for identifying and assembling the goat Y-linked scaffolds using existing data. A total of 11.5 Mb Y-linked sequences were clustered into 33 scaffolds, and 187 protein-coding genes were annotated. We also identified high abundance of repetitive elements. A 5.84 Mb subset was further ordered into an assembly with the evidence from the goat radiation hybrid map (RH map). The existing whole-genome resequencing data of 96 goats (worldwide distribution) were utilized to exploit the paternal relationships among bezoars and domestic goats. Goat paternal lineages were clearly divided into two clades (Y1 and Y2), predating the goat domestication. Demographic history analyses indicated that maternal lineages experienced a bottleneck effect around 2,000 YBP (years before present), after which goats belonging to the A haplogroup spread worldwide from the Near East. As opposed to this, paternal lineages experienced a population decline around the 10,000 YBP. The evidence from the Y chromosome suggests that male goats were not affected by the A haplogroup worldwide transmission, which implies sexually unbalanced contribution to the goat trade and population expansion in post-Neolithic period.

Insights into the Mitochondrial and Nuclear Genome Diversity of Two High Yielding Strains of Laying Hens

Simple Summary

Mitochondria are commonly known as “the powerhouse of the cell”, influencing the fitness, lifespan and metabolism of eukaryotic organisms. In our study we examined mitochondrial and nuclear genomic diversity in two high yielding strains of laying hens. We tested if the mitochondrial genome affects functional traits such as body weight and phosphorus utilization. We discovered a surprisingly low mitochondrial genetic diversity and an unequal distribution of the haplotypes among both strains, leading to limitations of robust links to phenotypic traits. In contrast, we found similar levels of nuclear genome diversity in both strains. Our study explores the potential influence of the mitochondrial genome on phenotypic traits and thus contributes to a better understanding of the function of this organelle in laying hens. Further, we focus on its usefulness as a genetic marker, which is often underestimated in breeding approaches, given the different inheritance mechanism compared to the nuclear genome.

Abstract

Mitochondria are essential components of eukaryotes as they are involved in several organismic key processes such as energy production, apoptosis and cell growth. Despite their importance for the metabolism and physiology of all eukaryotic organisms, the impact of mitochondrial haplotype variation has only been studied for very few species. In this study we sequenced the mitochondrial genome of 180 individuals from two different strains of laying hens. The resulting haplotypes were combined with performance data such as body weight, feed intake and phosphorus utilization to assess their influence on the hens in five different life stages. After detecting a surprisingly low level of genetic diversity, we investigated the nuclear genetic background to estimate whether the low mitochondrial diversity is representative for the whole genetic background of the strains. Our results highlight the need for more in-depth investigation of the genetic compositions and mito-nuclear interaction in individuals to elucidate the basis of phenotypic performance differences. In addition, we raise the question of how the lack of mitochondrial variation developed, since the mitochondrial genome represents genetic information usually not considered in breeding approaches.

Identifying the true number of specimens of the extinct blue antelope (Hippotragus leucophaeus)

Native to southern Africa, the blue antelope (Hippotragus leucophaeus) is the only large African mammal species known to have become extinct in historical times. However, it was poorly documented prior to its extinction ~ 1800 AD, and many of the small number of museum specimens attributed to it are taxonomically contentious. This places limitations on our understanding of its morphology, ecology, and the mechanisms responsible for its demise. We retrieved genetic information from ten of the sixteen putative blue antelope museum specimens using both shotgun sequencing and mitochondrial genome target capture in an attempt to resolve the uncertainty surrounding the identification of these specimens. We found that only four of the ten investigated specimens, and not a single skull, represent the blue antelope. This indicates that the true number of historical museum specimens of the blue antelope is even smaller than previously thought, and therefore hardly any reference material is available for morphometric, comparative and genetic studies. Our study highlights how genetics can be used to identify rare species in natural history collections where other methods may fail or when records are scarce. Additionally, we present an improved mitochondrial reference genome for the blue antelope as well as one complete and two partial mitochondrial genomes. A first analysis of these mitochondrial genomes indicates low levels of maternal genetic diversity in the ‘museum population’, possibly confirming previous results that blue antelope population size was already low at the time of the European colonization of South Africa.

Indonesian native goats (Capra hircus) reveal highest genetic frequency of mitochondrial DNA haplogroup B in the world

The objective of this study was to determine mtDNA sequences of the 481 bp HV1 region from two Indonesian native goat breeds, Kacang and Marica, to confirm the phylogeographic distribution of caprine haplogroup B in Southeast Asia. Based on these sequences, 12 haplotypes were observed and categorized into the predominant haplogroup B and minority haplogroup A, indicating that Indonesian native goats present the highest frequency (0.950) of the haplogroup B in the world. These results strongly emphasize previous observations of the haplogroup B frequencies tending to increase southeastward in Southeast Asia. Additionally, this suggests that goats primary bred into Southeast Asia might predominantly carry haplogroup B. This could essentially contribute to the understanding of the origin, propagation route and/or introgression history of Southeast Asian goats.

Assessing the Diversity and Population Substructure of Sarda Breed Bucks by Using Mtdna and Y-Chromosome Markers

A sample of 146 Sarda bucks from eight subregions of Sardinia, Italy (Nuorese, Barbagia, Baronia, Ogliastra, Sarrabus, Guspinese, Iglesiente, Sulcis) were characterized for Y-chromosome and mtDNA markers to assess the levels of population substructure. Five polymorphic loci (SRY, AMELY, ZFY, and DDX3Y) on the Y-chromosome were genotyped. The control region of mtDNA was sequenced as a source of complementary information. Analysis of Y-chromosome data revealed the segregation of 5 haplotypes: Y1A (66.43%), Y2 (28.57%), Y1C (3.57%), Y1B1 (0.71%), and Y1B2 (0.71%). High levels of Y-chromosome diversity were observed in populations from Southwest Sardinia. The FST values based on Y-chromosome and mtDNA data were low, although a paternal genetic differentiation was observed when comparing the Nuorese and Barbagia populations (Central Sardinia) with the Sulcis, Iglesiente, and Sarrabus populations (Southern Sardinia). AMOVA analysis supported the lack of population substructure. These results suggest the occurrence of a historical and extensive gene flow between Sarda goat populations from different locations of Sardinia, despite the fact that this island is covered by several large mountain ranges. Introgression with foreign caprine breeds in order to improve milk production might have also contributed to avoiding the genetic differentiation amongst Sarda populations.

Complete mitochondrial genomes of Karchaev goat (Capra hircus)

Karachaev goat (Capra hircus) is a local breed from North-Caucasus region, Russia. Here we present complete mitochondrial genome of Karachaev goat from the republic of Karachaevo-Cherkessia, Russia. The length of the studied sequence was 16,624 bp in size. It was shown that the studied specimen belonged to haplogroup A.

Maternal Origins and Haplotype Diversity of Seven Russian Goat Populations Based on the D-loop Sequence Variability

Simple Summary

Russia has diverse specifically selected and multipurpose goat resources. However, the origin of the local goats is still enigmatic. In this study, we sequenced and analyzed mitochondrial DNA (mtDNA) fragments of seven Russian local goat populations to provide the first insight into their maternal lineage.

Abstract

The territory of modern Russia lies on the crossroads of East and West and covers various geographical environments where diverse groups of local goats originated. In this work, we present the first study on the maternal origin of Russian local goats, including Altai Mountain (n = 9), Dagestan Downy (n = 18), Dagestan Local (n = 12), Dagestan Milk (n = 15), Karachaev (n = 21), Orenburg (n = 10), and Soviet Mohair (n = 7) breeds, based on 715 bp D-loop mitochondrial DNA (mtDNA) sequences. Saanen goats (n = 5) were used for comparison. Our findings reveal a high haplotype (HD = 0.843–1.000) and nucleotide diversity (π = 0.0112–0.0261). A total of 59 haplotypes were determined in the Russian goat breeds, in which all differed from the haplotypes of the Saanen goats. The haplotypes identified in Altai Mountain, Orenburg, Soviet Mohair, and Saanen goats were breed specific. Most haplotypes (56 of 59) were clustered together with samples belonging to haplogroup A, which was in accordance with the global genetic pattern of maternal origin seen in most goats worldwide. The haplotypes that were grouped together with rare haplogroups D and G were found in the Altai Mountain breed and haplogroup C was detected in the Soviet Mohair breed. Thus, our findings revealed that local goats might have been brought to Russia via various migration routes. In addition, haplotype sharing was found in aboriginal goat populations from overlapping regions, which might be useful information for their official recognition status.

Comprehensive analysis of mitochondrial DNA based genetic diversity in Indian goats

Thirty four distinct breeds and many non-descript populations represent the caprine diversity of India. Genetic characterization of breeds is an essential element in designing breeding strategies and preserving genetic diversity. Considering the popularity of mitochondrial DNA for phylogeographical studies, this study involved an extensive analysis of population structure and genetic diversity of 28 defined breeds and 5 lesser known populations representing all four major agro-climatic zones of India using mitochondrial DNA markers. Analysis of hypervariable region 1 of mtDNA control region in 443 goats together with 22 reference sequences, delineated 341 distinct haplotypes belonging to four maternal haplogroups; A, B, C and D, with haplogroup A representing 90% of the individuals. The haplotype and nucleotide diversity indices of Indian goats were 0.998 ± 0.001 and 0.028 ± 0.001, respectively indicating abundant genetic variability. Estimates of population demographic parameters from mismatch analysis suggested a relatively good fit to the model of either spatial or demographic expansion of Indian goats. AMOVA analysis and topology of MJ network suggested lack of phylogeographic structure in domestic goats, which can be attributed to unstructured animal breeding, dwindling pastures and nomadic pastoralism. Genetic differentiation between goats from different agro-ecological regions was in accordance with their geographical propinquity.

Illuminating Genetic Mysteries of the Dead Sea Scrolls

The discovery of the 2,000-year-old Dead Sea Scrolls had an incomparable impact on the historical understanding of Judaism and Christianity. "Piecing together" scroll fragments is like solving jigsaw puzzles with an unknown number of missing parts. We used the fact that most scrolls are made from animal skins to "fingerprint" pieces based on DNA sequences. Genetic sorting of the scrolls illuminates their textual relationship and historical significance. Disambiguating the contested relationship between Jeremiah fragments supplies evidence that some scrolls were brought to the Qumran caves from elsewhere; significantly, they demonstrate that divergent versions of Jeremiah circulated in parallel throughout Israel (ancient Judea). Similarly, patterns discovered in non-biblical scrolls, particularly the Songs of the Sabbath Sacrifice, suggest that the Qumran scrolls represent the broader cultural milieu of the period. Finally, genetic analysis divorces debated fragments from the Qumran scrolls. Our study demonstrates that interdisciplinary approaches enrich the scholar's toolkit.

High mitochondrial diversity of domesticated goats persisted among Bronze and Iron Age pastoralists in the Inner Asian Mountain Corridor

Goats were initially managed in the Near East approximately 10,000 years ago and spread across Eurasia as economically productive and environmentally resilient herd animals. While the geographic origins of domesticated goats (Capra hircus) in the Near East have been long-established in the zooarchaeological record and, more recently, further revealed in ancient genomes, the precise pathways by which goats spread across Asia during the early Bronze Age (ca. 3000 to 2500 cal BC) and later remain unclear. We analyzed sequences of hypervariable region 1 and cytochrome b gene in the mitochondrial genome (mtDNA) of goats from archaeological sites along two proposed transmission pathways as well as geographically intermediary sites. Unexpectedly high genetic diversity was present in the Inner Asian Mountain Corridor (IAMC), indicated by mtDNA haplotypes representing common A lineages and rarer C and D lineages. High mtDNA diversity was also present in central Kazakhstan, while only mtDNA haplotypes of lineage A were observed from sites in the Northern Eurasian Steppe (NES). These findings suggest that herding communities living in montane ecosystems were drawing from genetically diverse goat populations, likely sourced from communities in the Iranian Plateau, that were sustained by repeated interaction and exchange. Notably, the mitochondrial genetic diversity associated with goats of the IAMC also extended into the semi-arid region of central Kazakhstan, while NES communities had goats reflecting an isolated founder population, possibly sourced via eastern Europe or the Caucasus region.

Georgian cattle, sheep, goats: are they of Near-Eastern origins?

The main aim of this research was to establish the nucleotide sequence of the highly variable region of the D loop of the mitochondrial DNA of some Georgian domestic animal species (cattle, goat, sheep) as well as their phylogenetic position among the worldwide set of domestic animals. In this study, a total of 5 haplogroups (T – 5; T3 – 7; T1 – 1; T2 – 2; T5 – 2) in 17 Georgian Mountain cattle (GMC), 4 haplogroups (A – 15; A2a1 – 3; A1a – 1; A6 – 3) in 22 Georgian goats and 3 haplogroups (A – 10; B – 16; C -15) in 41 Georgian sheeps (15 Imeretian and 26 Tushetian) were detected. This study represents the first attempt of Genetic study of native Georgian livestocks. The GMC, Georgian (Megrelian) goat, Georgian (Imeretian and Tushetian) sheep mitogenomes were grouped phylogenetically in the haplogroups indicating the closeness to the Near Eastern animals.

Mitochondrial DNA variation and phylogeography of native Mongolian goats

OBJECTIVE

Mongolia is one of a few countries that supports over 25 million goats, but genetic diversity, demographic history, and the origin of goat populations in Mongolia have not been well studied. This study was conducted to assess the genetic diversity, phylogenetic status and population structure of Mongolian native goats, as well as to discuss their origin together with other foreign breeds from different countries using hypervariable region 1 (HV1) in mtDNA.

METHODS

In this study, we examined the genetic diversity and phylogenetic status of Mongolian native goat populations using a 452 base-pair long fragment of HVI of mitochondrial DNA from 174 individuals representing 12 populations. In addition, 329 previously published reference sequences from different regions were included in our phylogenetic analyses.

RESULTS

Investigated native Mongolian goats displayed relatively high genetic diversities. After sequencing, we found a total of 109 polymorphic sites that defined 137 haplotypes among investigated populations. Of these, haplotype and nucleotide diversities of Mongolian goats were calculated as 0.997±0.001 and 0.0283±0.002, respectively. These haplotypes clearly clustered into four haplogroups (A, B, C, and D), with the predominance of haplogroup A (90.8%). Estimates of pairwise differences (Fst) and the analysis of molecular variance values among goat populations in Mongolia showed low genetic differentiation and weak geographical structure. In addition, Kazakh, Chinese (from Huanghuai and Leizhou), and Arabian (Turkish and Baladi breeds) goats had smaller genetic differentiation compared to Mongolian goats.

CONCLUSION

In summary, we report novel information regarding genetic diversity, population structure, and origin of Mongolian goats. The findings obtained from this study reveal that abundant haplogroups (A to D) occur in goat populations in Mongolia, with high levels of haplotype and nucleotide diversity.

Abundant Genetic Diversity of Yunling Cattle Based on Mitochondrial Genome

Simple Summary

Mitochondrial DNA (mtDNA) analysis is a critical tool in assessing the maternal origin, phylogeny and population structure of domestic animals. Yunling cattle are a composite breed that was created by local Yunnan cattle breeds (as a maternal line) 30 years ago. It can be said that Yunling cattle represent important reservoirs of genetic diversity of local Yunnan cattle. Yunling cattle are characterized as a tropical/subtropical breed with fast growth, heat tolerance and parasite resistance, and can survive in a harsh environment and on low-quality roughage. Assessing the genetic characteristics of Yunling cattle is of particular importance for reasonable breeding strategies for Yunling cattle and the design of a local Yunnan cattle conservation program.

Abstract

Yunling cattle are a composite beef cattle breed, combining Brahman (1/2), Murray Grey (1/4) and Local Yunnan cattle (1/4), that was developed in Yunnan, China in the 1980s. Understanding the genetic information of Yunling cattle is of great significance to the development of reasonable breeding strategies for this breed. In the present study, we assessed the current genetic status of Yunling cattle in Yunnan Province (China) by analyzing the variability of the whole mitochondrial genome of 129 individuals. Altogether, 129 sequences displayed 47 different haplotypes. The haplotype diversity and the average number of nucleotide differences were 0.964 and 128.074, respectively. Phylogenetic analyses classified Yunling cattle into seven haplogroups: T1, T2, T3, T4, T6, I1 and I2. Haplogroup I1 was found to be predominant (41.86%), followed by T3 (28.68%). Furthermore, we also identified a novel haplogroup, T6, and defined the sub-haplogroup I1a in Yunling cattle. According to the formation process of Yunling cattle (local Yunnan cattle as the maternal line), the high genetic diversities in the mitochondria of Yunling cattle could be due to the complex maternal origin of local Yunnan cattle. Further studies about local Yunnan breeds are necessary to determine the exact source of haplogroup T6 in Yunling cattle. Our results will be useful for the evaluation and effective management of Yunling cattle.

Maternal genetic and phylogenetic characteristics of domesticated cattle in northwestern China

Northwestern China, an important part of the Silk Road, was the birthplace of the Zhouzu farming culture. The domestication of cattle as an important aspect of farming culture has had a long history in northwestern China. In this study, we assessed the maternal structure and phylogeny of cattle by analyzing the mitochondrial DNA hypervariable segment I (HVS-I) in 698 native cattle from eight areas of northwestern China. The phylogenetic analyses revealed two highly divergent mtDNA clades: clade T, which had four sub-clades (Ta—Td), and clade I. The cattle domesticated from Bos taurus showed a clear dominant distribution pattern in northwestern China. The nucleotide diversity of the Bos indicus clade was lower than that of clades from Bos taurus. In summary, our results suggest that the native cattle of northwestern China were domesticated from two different maternal ancestors, Bos taurus and Bos indicus, which migrated to the central plains of China from the north and south, respectively, with Bos taurus remaining at the edges of the region. The population expansion of the cattle domesticated from Bos taurus occurred in the Longdong region of Gansu Province, and these cattle formed four relatively independent evolutionary branches. Subsequent to this expansion event, Bos indicus migrated from southern to northern China.

Genetic diversity analysis of French goat populations reveals selective sweeps involved in their differentiation

After domestication 11 000 years ago in Asia Minor, the goat followed human migration to Europe and Asia. It was then introduced in Africa and is now raised all over the world. In this study, we exploited a dataset composed of 54 000 SNPs (Illumina goat DNA chip) to analyze the genetic diversity of 223 individuals belonging to eight French breeds (Alpine, Angora, Corse, Fossés, Poitevine, Provençale, Pyrénées and Saanen). Analyses carried out included individual-based approaches (principal component analysis and population structure) and population-based approaches (phylogenetic tree constructions). The results of the genetic diversity analyses revealed that French breeds are clearly differentiated, in particular, the Angora breed that originates from south west Asia. The Provençale breed shows a very original genetic pattern that could be the result of ancient admixture. Then, selection signatures were detected by identifying regions of outlying genetic differentiation between populations. Five genomic regions were detected under selection on chromosomes 5, 6, 11, 13 and 20, revealing mainly soft selective sweeps and a few hard selective sweeps and highlighting candidate genes that had been selected for during the evolutionary history of these breeds. Among them, two coat coloration genes (ADAMTS20 and ASIP) and one gene related to milk composition (CSN1S1) were involved.

Genome-wide SNP profiling of worldwide goat populations reveals strong partitioning of diversity and highlights post-domestication migration routes

BACKGROUND

Goat populations that are characterized within the AdaptMap project cover a large part of the worldwide distribution of this species and provide the opportunity to assess their diversity at a global scale. We analysed genome-wide 50 K single nucleotide polymorphism (SNP) data from 144 populations to describe the global patterns of molecular variation, compare them to those observed in other livestock species, and identify the drivers that led to the current distribution of goats.

RESULTS

A high degree of genetic variability exists among the goat populations studied. Our results highlight a strong partitioning of molecular diversity between and within continents. Three major gene pools correspond to goats from Europe, Africa and West Asia. Dissection of sub-structures disclosed regional gene pools, which reflect the main post-domestication migration routes. We also identified several exchanges, mainly in African populations, and which often involve admixed and cosmopolitan breeds. Extensive gene flow has taken place within specific areas (e.g., south Europe, Morocco and Mali-Burkina Faso-Nigeria), whereas elsewhere isolation due to geographical barriers (e.g., seas or mountains) or human management has decreased local gene flows.

CONCLUSIONS

After domestication in the Fertile Crescent in the early Neolithic era (ca. 12,000 YBP), domestic goats that already carried differentiated gene pools spread to Europe, Africa and Asia. The spread of these populations determined the major genomic background of the continental populations, which currently have a more marked subdivision than that observed in other ruminant livestock species. Subsequently, further diversification occurred at the regional level due to geographical and reproductive isolation, which was accompanied by additional migrations and/or importations, the traces of which are still detectable today. The effects of breed formation were clearly detected, particularly in Central and North Europe. Overall, our results highlight a remarkable diversity that occurs at the global scale and is locally partitioned and often affected by introgression from cosmopolitan breeds. These findings support the importance of long-term preservation of goat diversity, and provide a useful framework for investigating adaptive introgression, directing genetic improvement and choosing breeding targets.

Genetic diversity and maternal lineages of south Indian goats

In India, goats are considered to be one of the important livestock species that reinforce the rural economy. Even though India has 23 well-recognized goat breeds, the knowledge about their genetic diversity and domestication history is limited. In this study, we have analyzed the genetic diversity of 104 goats representing five different south Indian breeds using mtDNA D-loop region. The haplotype diversity of the breeds ranged from 0.9524 ± 0.0403 (Malabari) to 0.9921 ± 0.0154 (Kanni aadu). Analysis of molecular variance showed only 1.05% variation among breeds. On other hand, the variation within breed was remarkably high (98.95%) which suggested the weak phylogeographic structure of south Indian goats. The phylogenetic analysis revealed three haplogroups representing maternal lineages namely A, B and D. The analysis of 466 Indian goat sequences showed an additional lineage C. As reported in the previous studies, a major fraction of analyzed goats fell into haplogroup A. Our study confirms the presence of three maternal lineages for south Indian domestic goats.

Ancient goat genomes reveal mosaic domestication in the Fertile Crescent

Current genetic data are equivocal as to whether goat domestication occurred multiple times or was a singular process. We generated genomic data from 83 ancient goats (51 with genome-wide coverage) from Paleolithic to Medieval contexts throughout the Near East. Our findings demonstrate that multiple divergent ancient wild goat sources were domesticated in a dispersed process that resulted in genetically and geographically distinct Neolithic goat populations, echoing contemporaneous human divergence across the region. These early goat populations contributed differently to modern goats in Asia, Africa, and Europe. We also detect early selection for pigmentation, stature, reproduction, milking, and response to dietary change, providing 8000-year-old evidence for human agency in molding genome variation within a partner species.

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

BACKGROUND

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

METHODOLOGY

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

CONCLUSION

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

Mitochondrial DNA variation reveals maternal origins and demographic dynamics of Ethiopian indigenous goats

The Horn of Africa forms one of the two main historical entry points of domestics into the continent and Ethiopia is particularly important in this regard. Through the analysis of mitochondrial DNA (mtDNA) d-loop region in 309 individuals from 13 populations, we reveal the maternal genetic variation and demographic dynamics of Ethiopian indigenous goats. A total of 174 variable sites that generated 231 haplotypes were observed. They defined two haplogroups that were present in all the 13 study populations. Reference haplotypes from the six globally defined goat mtDNA haplogroups show the two haplogroups present in Ethiopia to be A and G, the former being the most predominant. Although both haplogroups are characterized by an increase in effective population sizes (Ne) predating domestication, they also have experienced a decline in Ne at different time periods, suggesting different demographic histories. We observed seven haplotypes, six were directly linked to the central haplotypes of the two haplogroups and one was central to haplogroup G. The seven haplotypes were common between Ethiopia, Kenya, Egypt, and Saudi Arabia populations, suggesting common maternal history and the introduction of goats into East Africa via Egypt and the Arabian Peninsula, respectively. While providing new mtDNA data from a historically important region, our results suggest extensive intermixing of goats mediated by human socio-cultural and economic interactions. These have led to the coexistence of the two haplogroups in different geographic regions in Ethiopia resulting in a large caprine genetic diversity that can be exploited for genetic improvement.

Genetic origin of goat populations in Oman revealed by mitochondrial DNA analysis

The Sultanate of Oman has a complex mosaic of livestock species and production systems, but the genetic diversity, demographic history or origins of these Omani animals has not been expensively studied. Goats might constitute one of the most abundant and important domestic livestock species since the Neolithic transition. Here, we examined the genetic diversity, origin, population structure and demographic history of Omani goats. Specifically, we analyzed a 525-bp fragment of the first hypervariable region of the mitochondrial DNA (mtDNA) control region from 69 Omani individuals and compared this fragment with 17 mtDNA sequences from Somalia and Yemen as well as 18 wild goat species and 1,198 previously published goat sequences from neighboring countries. The studied goat breeds show substantial diversity. The haplotype and nucleotide diversities of Omani goats were found equal to 0.983 ± 0.006 and 0.0284 ± 0.014, respectively. The phylogenetic analyses allowed us to classify Omani goats into three mtDNA haplogroups (A, B and G): haplogroup A was found to be predominant and widely distributed and accounted for 80% of all samples, and haplogroups B and G exhibited low frequencies. Phylogenetic comparisons with wild goats revealed that five of the native Omani goat populations originate from Capra aegagrus. Furthermore, most comparisons of pairwise population F_ST values within and between these five Omani goat breeds as well as between Omani goats and nine populations from nearby countries were not significant. These results suggest strong gene flow among goat populations caused by the extensive transport of goats and the frequent movements of human populations in ancient Arabia. The findings improve our understanding of the migration routes of modern goats from their region of domestication into southeastern Arabia and thereby shed light on human migratory and commercial networks during historical times.

Capturing goats: documenting two hundred years of mitochondrial DNA diversity among goat populations from Britain and Ireland

The domestic goat (Capra hircus) plays a key role in global agriculture, being especially prized in regions of marginal pasture. However, the advent of industrialized breeding has seen a dramatic reduction in genetic diversity within commercial populations, while high extinction rates among feral herds have further depleted the reservoir of genetic variation available. Here, we present the first survey of whole mitochondrial genomic variation among the modern and historical goat populations of Britain and Ireland using a combination of mtDNA enrichment and high throughput sequencing. Fifteen historical taxidermy samples, representing the indigenous 'Old Goat' populations of the islands, were sequenced alongside five modern Irish dairy goats and four feral samples from endangered populations in western Ireland. Phylogenetic and network analyses of European mitochondrial variation revealed distinct groupings dominated by historical British and Irish samples, which demonstrate a degree of maternal genetic structure between the goats of insular and continental Europe. Several Irish modern feral samples also fall within these clusters, suggesting continuity between these dwindling populations and the ancestral 'Old Goats' of Ireland and Britain.

Goat domestication and breeding: a jigsaw of historical, biological and molecular data with missing pieces

Domestic goats (Capra hircus) are spread across the five continents with a census of 1 billion individuals. The worldwide population of goats descends from a limited number of bezoars (Capra aegagrus) domesticated 10 000 YBP (years before the present) in the Fertile Crescent. The extraordinary adaptability and hardiness of goats favoured their rapid spread over the Old World, reaching the Iberian Peninsula and Southern Africa 7000 YBP and 2000 YBP respectively. Molecular studies have revealed one major mitochondrial haplogroup A and five less frequent haplogroups B, C, D, F and G. Moreover, the analysis of autosomal and Y-chromosome markers has evidenced an appreciable geographic differentiation. The implementation of new molecular technologies, such as whole-genome sequencing and genome-wide genotyping, allows for the exploration of caprine diversity at an unprecedented scale, thus providing new insights into the evolutionary history of goats. In spite of a number of pitfalls, the characterization of the functional elements of the goat genome is expected to play a key role in understanding the genetic determination of economically relevant traits. Genomic selection and genome editing also hold great potential, particularly for improving traits that cannot be modified easily by traditional selection.

Mitogenome Sequencing in the Genus Camelus Reveals Evidence for Purifying Selection and Long-term Divergence between Wild and Domestic Bactrian Camels

The genus Camelus is an interesting model to study adaptive evolution in the mitochondrial genome, as the three extant Old World camel species inhabit hot and low-altitude as well as cold and high-altitude deserts. We sequenced 24 camel mitogenomes and combined them with three previously published sequences to study the role of natural selection under different environmental pressure, and to advance our understanding of the evolutionary history of the genus Camelus. We confirmed the heterogeneity of divergence across different components of the electron transport system. Lineage-specific analysis of mitochondrial protein evolution revealed a significant effect of purifying selection in the concatenated protein-coding genes in domestic Bactrian camels. The estimated dN/dS < 1 in the concatenated protein-coding genes suggested purifying selection as driving force for shaping mitogenome diversity in camels. Additional analyses of the functional divergence in amino acid changes between species-specific lineages indicated fixed substitutions in various genes, with radical effects on the physicochemical properties of the protein products. The evolutionary time estimates revealed a divergence between domestic and wild Bactrian camels around 1.1 [0.58-1.8] million years ago (mya). This has major implications for the conservation and management of the critically endangered wild species, Camelus ferus.

Whole Mitogenomes Reveal the History of Swamp Buffalo: Initially Shaped by Glacial Periods and Eventually Modelled by Domestication

The newly sequenced mitochondrial genomes of 107 Asian swamp buffalo (Bubalus bubalis carabensis) allowed the reconstruction of the matrilineal divergence since ~900 Kya. Phylogenetic trees and Bayesian skyline plots suggest a role of the glacial periods in the demographic history of swamp buffalo. The ancestral swamp-buffalo mitogenome is dated ~232 ± 35 Kya. Two major macro-lineages diverged during the 2^nd Pleistocene Glacial Period (~200-130 Kya), but most (~99%) of the current matrilines derive from only two ancestors (SA1'2 and SB) that lived around the Last Glacial Maximum (~26-19 Kya). During the late Holocene optimum (11-6 Kya) lineages differentiated further, and at least eight matrilines (SA1, SA2, SB1a, SB1b, SB2a, SB2b, SB3 and SB4) were domesticated around 7-3 Kya. Haplotype distributions support an initial domestication process in Southeast Asia, while subsequent captures of wild females probably introduced some additional rare lineages (SA3, SC, SD and SE). Dispersal of domestic buffaloes created local population bottlenecks and founder events that further differentiated haplogroup distributions. A lack of maternal gene flow between neighboring populations apparently maintained the strong phylogeography of the swamp buffalo matrilines, which is the more remarkable because of an almost complete absence of phenotypic differentiation.

Polymorphisms, differentiation, and phylogeny of 10 Tibetan goat populations inferred from mitochondrial D-loop sequences

China has abundant population of Tibetan goats, but available information regarding genetic structure and phylogenetic status remains limited. Here, 130 mitochondrial D-loop sequences of individuals from 10 Tibetan goat populations located in distinct plateau areas were analyzed. Eighty-six haplotypes were defined, among which 97.7% were group-specific haplotypes. Haplotype and nucleotide diversity indices were 0.990 ± 0.003 and 0.0145 ± 0.0013, respectively. The pairwise Wright's F-statistics ranged from -0.028 to 0.385, and over half of them were greater than 0.05, indicating apparent genetic differentiation among the populations. AMOVA analysis (F_ST = 0.0858) manifested that the genetic structure has become weak. Phylogenetic trees revealed four haplogroups (A, B, C, and D), suggesting that Tibetan goats had four origins. Mismatch distribution analyses and neutrality tests indicated that at least one population expansion event occurred during the demographic history of Tibetan goat. These results will provide a more complete understanding of Tibetan goat genetic resources.

The Valdostana goat: a genome-wide investigation of the distinctiveness of its selective sweep regions

The Valdostana goat is an alpine breed, raised only in the northern Italian region of the Aosta Valley. This breed's main purpose is to produce milk and meat, but is peculiar for its involvement in the "Batailles de Chèvres," a recent tradition of non-cruel fight tournaments. At both the genetic and genomic levels, only a very limited number of studies have been performed with this breed and there are no studies about the genomic signatures left by selection. In this work, 24 unrelated Valdostana animals were screened for runs of homozygosity to identify highly homozygous regions. Then, six different approaches (ROH comparison, Fst single SNPs and windows based, Bayesian, Rsb, and XP-EHH) were applied comparing the Valdostana dataset with 14 other Italian goat breeds to confirm regions that were different among the comparisons. A total of three regions of selection that were also unique among the Valdostana were identified and located on chromosomes 1, 7, and 12 and contained 144 genes. Enrichment analyses detected genes such as cytokines and lymphocyte/leukocyte proliferation genes involved in the regulation of the immune system. A genetic link between an aggressive challenge, cytokines, and immunity has been hypothesized in many studies both in humans and in other species. Possible hypotheses associated with the signals of selection detected could be therefore related to immune-related factors as well as with the peculiar battle competition, or other breed-specific traits, and provided insights for further investigation of these unique regions, for the understanding and safeguard of the Valdostana breed.

Understanding the evolution of holoparasitic plants: the complete plastid genome of the holoparasite Cytinus hypocistis (Cytinaceae)

Background and Aims Plant plastid genomes are highly conserved in size, gene content and structure; however, parasitic plants are a noticeable exception to this evolutionary stability. Although the evolution of parasites could help to better understand plastome evolution in general, complete plastomes of parasites have been sequenced only for some lineages so far. Here we contribute to filling this gap by providing and analysing the complete plastome sequence of Cytinus hypocistis, the first parasite sequenced for Malvales and a species suspected to have an extremely small genome. Methods We sequenced and assembled de novo the plastid genome of Cytinus hypocistis using a shotgun approach on genomic DNA. Phylogenomic analyses based on coding regions were performed on Malvidae. For each coding region present in Cytinus, we tested for relaxation or intensification of selective pressures in the Cytinus lineage compared with autotrophic Malvales. Key Results Cytinus hypocistis has an extremely divergent genome that is among the smallest sequenced to date (19·4 kb), with only 23 genes and no inverted repeat regions. Phylogenomic analysis confirmed the position of Cytinus within Malvales. All coding regions of Cytinus plastome presented very high substitution rates compared with non-parasitic Malvales. Conclusions Some regions were inferred to be under relaxed negative selection in Cytinus, suggesting that further plastome reduction is occurring due to relaxed purifying selection associated with the loss of photosynthetic activity. On the other hand, increased selection intensity and strong positive selection were detected for rpl22 in the Cytinus lineage, which might indicate an evolutionary role in the host-parasite arms race, a point that needs further research.

A whole mitochondria analysis of the Tyrolean Iceman's leather provides insights into the animal sources of Copper Age clothing

The attire of the Tyrolean Iceman, a 5,300-year-old natural mummy from the Ötzal Italian Alps, provides a surviving example of ancient manufacturing technologies. Research into his garments has however, been limited by ambiguity surrounding their source species. Here we present a targeted enrichment and sequencing of full mitochondrial genomes sampled from his clothes and quiver, which elucidates the species of production for nine fragments. Results indicate that the majority of the samples originate from domestic ungulate species (cattle, sheep and goat), whose recovered haplogroups are now at high frequency in today’s domestic populations. Intriguingly, the hat and quiver samples were produced from wild species, brown bear and roe deer respectively. Combined, these results suggest that Copper Age populations made considered choices of clothing material from both the wild and domestic populations available to them. Moreover, these results show the potential for the recovery of complete mitochondrial genomes from degraded prehistoric artefacts.