Sheep Mitochondrial DNA Variation in European, Caucasian, and Central Asian Areas

The evolution of contemporary livestock species: Insights from mitochondrial genome

The domestication of animals marks a pivotal moment in human history, profoundly influencing our demographic and cultural progress. This process has led to significant genetic, behavioral, and physical changes in livestock species compared to their wild ancestors. Understanding the evolutionary history and genetic diversity of livestock species is crucial, and mitochondrial DNA (mtDNA) has emerged as a robust marker for investigating molecular diversity in animals. Its highly conserved gene content across animal species, minimal duplications, absence of introns, and short intergenic regions make mtDNA analysis ideal for such studies. Mitochondrial DNA analysis has uncovered distinct cattle domestication events dating back to 8000 years BC in Southwestern Asia. The sequencing of water buffalo mtDNA in 2004 provided important insights into their domestication history. Caprine mtDNA analysis identified three haplogroups, indicating varied maternal origins. Sheep, domesticated 12,000 years ago, exhibit diverse mtDNA lineages, suggesting multiple domestication events. Ovine mtDNA studies revealed clades A, B, C, and a fourth lineage, group D. The origins of domestic pigs were traced to separate European and Asian events followed by interbreeding. In camels, mtDNA elucidated the phylogeographic structure and genetic differentiation between wild and domesticated species. Horses, domesticated around 3500 BC, show significant mtDNA variability, highlighting their diverse origins. Yaks exhibit unique adaptations for high-altitude environments, with mtDNA analysis providing insights into their adaptation. Chicken mtDNA studies supported a monophyletic origin from Southeast Asia's red jungle fowl, with evidence of multiple origins. This review explores livestock evolution and diversity through mtDNA studies, focusing on cattle, water buffalo, goat, sheep, pig, camel, horse, yak and chicken. It highlights mtDNA's significance in unraveling maternal lineages, genetic diversity, and domestication histories, concluding with insights into its potential application in improving livestock production and reproduction dynamics.

Analysing the genetic diversity of three sheep breeds in Turkey and nearby countries using 50 K SNPs data

This study analysed the genetic diversity and population structure of eight sheep breeds in Turkey and nearby countries. Moderate genetic diversity was observed, with the Sakiz (SKZ) exhibiting the highest diversity based on heterozygosity and allelic richness (AR) values. Genetic distances revealed differentiation between the populations, with the most significant divergence between the Cyprus Fat Tail (CFT) and SKZ breeds. PCA demonstrated SKZ and Chios (CHI) clustering together, indicating genetic similarity. Karakas (KRS), Norduz (NDZ), Afshari (AFS), Moghani (MOG) and others showed overlap, reflecting genetic relationships. Ancestry analysis found that KRS was predominantly inherited from the second ancestral population, while SKZ and NDZ were primarily derived from the first and second ancestral lineages. This illustrated the populations' diverse origins. Most genetic variation (96.84%) was within, not between, populations. The phi-statistic (PhiPT) indicated moderate differentiation overall. Phylogenetic analysis further demonstrated the genetic distinctiveness of the SKZ breed. ROH and FROH analyses showed that SKZ exhibited the highest homozygosity and inbreeding, while KRS displayed the lowest. This study elucidates these breeds' genetic diversity, structure and relationships. Key findings include moderate diversity, evidence of differentiation between breeds, diverse ancestral origins and distinct ROH patterns. This provides insights into the population's genetic characteristics and conservation requirements.

Whole-genome resequencing of Hu sheep identifies candidate genes associated with agronomic traits

The phenotypic diversity resulting from artificial or natural selection of sheep has made a significant contribution to human civilization. Hu sheep are a local sheep breed unique to China with high reproductive rates and rapid growth. Genomic selection signatures have been widely used to investigate the genetic mechanisms underlying phenotypic variation in livestock. Here, we conduct whole-genome sequencing of 207 Hu sheep and compare them with the wild ancestors of domestic sheep (Asiatic mouflon) to investigate the genetic characteristics and selection signatures of Hu sheep. Based on six signatures of selection approaches, we detect genomic regions containing genes related to reproduction (BMPR1B, BMP2, PGFS, CYP19, CAMK4, GGT5, and GNAQ), vision (ALDH1A2, SAG, and PDE6B), nervous system (NAV1), and immune response (GPR35, SH2B2, PIK3R3, and HRAS). Association analysis with a population of 1299 Hu sheep reveals that those missense mutations in the GPR35 (GPR35 g.952651 A>G; GPR35 g.952496 C>T) and NAV1 (NAV1 g.84216190 C>T; NAV1 g.84227412 G>A) genes are significantly associated (P < 0.05) with immune and growth traits in Hu sheep, respectively. This research offers unique insights into the selection characteristics of Hu sheep and facilitates further genetic improvement and molecular investigations.

Review of sheep breeding and genetic research in Türkiye

The livestock industry in Türkiye is vital to the country's agricultural sector and economy. In particular, sheep products are an important source of income and livelihood for many Turkish smallholder farmers in semi-arid and highland areas. Türkiye is one of the largest sheep producers in the world and its sheep production system is heavily dependent on indigenous breeds. Given the importance of the sheep industry in Türkiye, a systematic literature review on sheep breeding and genetic improvement in the country is needed for the development and optimization of sheep breeding programs using modern approaches, such as genomic selection. Therefore, we conducted a comprehensive literature review on the current characteristics of sheep populations and farms based on the most up-to-date census data and breeding and genetic studies obtained from scientific articles. The number of sheep has increased in recent years, mainly due to the state's policy of supporting livestock farming and the increase in consumer demand for sheep dairy products with high nutritional and health benefits. Most of the genetic studies on indigenous Turkish sheep have been limited to specific traits and breeds. The use of genomics was found to be incipient, with genomic analysis applied to only two major breeds for heritability or genome-wide association studies. The scope of heritability and genome-wide association studies should be expanded to include traits and breeds that have received little or no attention. It is also worth revisiting genetic diversity studies using genome-wide single nucleotide polymorphism markers. Although there was no report of genomic selection in Turkish sheep to date, genomics could contribute to overcoming the difficulties of implementing traditional pedigree-based breeding programs that require accurate pedigree recording. As indigenous sheep breeds are better adapted to the local environmental conditions, the proper use of breeding strategies will contribute to increased income, food security, and reduced environmental footprint in a sustainable manner.

Phenotype transition from wild mouflon to domestic sheep

The domestication of animals started around 12,000 years ago in the Near East region. This "endless process" is characterized by the gradual accumulation of changes that progressively marked the genetic, phenotypic and physiological differences between wild and domesticated species. The main distinctive phenotypic characteristics are not all directly attributable to the human-mediated selection of more productive traits. In the last decades, two main hypotheses have been proposed to clarify the emergence of such a set of phenotypic traits across a variety of domestic species. The first hypothesis relates the phenotype of the domesticated species to an altered thyroid hormone-based signaling, whereas the second one relates it to changes in the neural crest cells induced by selection of animals for tameness. These two hypotheses are not necessarily mutually exclusive since they may have contributed differently to the process over time and space. The adaptation model induced by domestication can be adopted to clarify some aspects (that are still controversial and debated) of the long-term evolutionary process leading from the wild Neolithic mouflon to the current domestic sheep. Indeed, sheep are among the earliest animals to have been domesticated by humans, around 12,000 years ago, and since then, they have represented a crucial resource in human history. The aim of this review is to shed light on the molecular mechanisms and the specific genomic variants that underlie the phenotypic variability between sheep and mouflon. In this regard, we carried out a critical review of the most recent studies on the molecular mechanisms that are most accredited to be responsible for coat color and phenotype, tail size and presence of horns. We also highlight that, in such a complicate context, sheep/mouflon hybrids represent a powerful and innovative model for studying the mechanism by which the phenotypic traits related to the phenotypic responses to domestication are inherited. Knowledge of these mechanisms could have a significant impact on the selection of more productive breeds. In fact, as in a journey back in time of animal domestication, the genetic traits of today's domestic species are being progressively and deliberately shaped according to human needs, in a direction opposite to that followed during domestication.

Ancient mitochondrial genome depicts sheep maternal dispersal and migration in Eastern Asia

Sheep have been one of the most important groups of animals since ancient times. However, the knowledge of their migration routes and genetic relationships is still poorly understood. To investigate sheep maternal migration histories alongside Eurasian communications routes, in this study, we obtain mitochondrial genomes (mitogenomes) from 17 sheep remains in 6 Chinese sites and 1 Uzbekistan site dated 4429-3100 years before present (BP). By obtaining the mitogenomes from the sheep (4429-3556 BP) found in the Tongtian Cave site in Xinjiang, Altai region of northwest China, our results support the emergence of haplogroup C sheep in Xinjiang as early as 4429-3556 BP. The combined phylogenetic analyses with extant ancient and modern sheep mitogenomes suggest that the Uzbekistan-Altai region may have been a migration hub for early sheep in eastern Asia. At least two migration events have taken place for sheep crossing Eurasia to China, one passing by Uzbekistan and Northwest China to the middle and lower reaches of the Yellow River at approximately 4000 BP and another following the Altai region to middle Inner Mongolia from 4429 BP to 2500 BP. Overall, this study provides further evidence for early sheep utilization and migration patterns in Eastern Asia.

Comparison of genetic diversity and phylogenetic structure of BRCA1 gene of some domestic and wild sheep breeds in different countries

BRCA1 gene plays an important role in DNA damage repair, cell cycle, and transcription process regulation; hence it's called gate keeper. The current research aims to perform bioinformatics analyzes of the BRCA1 gene of different breeds of domestic and wild sheep from 49 breeds in 14 countries using the NCBI genome database. The desired sequences were aligned using MEGA11 software and a phylogenetic tree was drawn by Neighbor-Joining method. The number of mutations, nucleotide diversity, and haploid diversity were also analyzed using Dnaspv5 software. The analyses showed 296 polymorphisms, which led to the creation of 45 different haplotypes with a haplotype diversity of 0.035. Nucleotide diversity and average nucleotide differences among breeds were estimated as 0.259 and 0.052, respectively. The average genetic distance within the population of countries was calculated as 0.052. The amount of sequence conservation in this research was 0.313 on average, which indicates the high polymorphism of this gene and the emergence of new proteins. Tajima's D value in Tajima's neutrality test was -2.421, which was significant (p < 0.05). One of the reasons for the high genetic diversity in Iran's wild sheep population is the existence of forests and open environments, which prevent genetic drift and reduce inbreeding.

Assessment of genetic diversity of the fat-tailed Dumba sheep of India by mitochondrial and microsatellite markers

India has a centuries-old tradition of sheep production and breeding that accomplish economic, agricultural, and religious roles. In addition to the 44 registered sheep breeds, there is a fat-tailed sheep population referred to as Dumba. This study evaluated genetic variation in Dumba sheep and its differentiation from other Indian sheep breeds using mitochondrial DNA and genomic microsatellite loci. Haplotype and nucleotide diversity based on mitochondrial DNA analysis revealed substantially high maternal genetic diversity in Dumba sheep. Major ovine haplogroups A and B observed in sheep populations across the globe registered their presence in the Dumba sheep. The molecular genetic analysis using microsatellite markers also showed high measures of allele (10.125 ± 0.762) and gene diversity (0.749 ± 0.029). Results correspond to the non-bottleneck population that is near mutation-drift equilibrium despite some deficiency in the number of heterozygotes (FIS = 0.043 ± 0.059). Phylogenetic clustering confirmed Dumba to be a distinct population. Results of this study endow authorities with critical information imperative for sustainable utilization and conservation of Indian fat-tailed sheep, which is considered to be an untapped genetic resource contributing to the food security, livelihood, and economic sustainability of rural households in marginal areas of the country.

Phylogenetic Analysis of Russian Native Sheep Breeds Based on mtDNA Sequences

Eurasia is represented by all climatic zones and various environments. A unique breed variety of farm animals has been developed in Russia, whose territory covers a large area of the continent. A total of 69 local breeds and types of dairy, wool, and meat sheep (Ovis aries) are maintained here. However, the genetic diversity and maternal origin of these local breeds have not been comprehensively investigated. In this study, we describe the diversity and phylogeny of Russian sheep breeds inhabiting different geographical regions based on the analysis of complete sequences of mitochondrial genomes (mtDNA). Complete mtDNA sequences of the studied sheep were obtained using next-generation sequencing technology (NGS). All investigated geographical groups of sheep were characterized by high haplotype (Hd = 0.9992) and nucleotide diversity (π = 0.00378). Analysis of the AMOVA results showed that genetic diversity was majorly determined by within-population differences (77.87%). We identified 128 haplotypes in all studied sheep. Haplotypes belonged to the following haplogroups: B (64.8%), A (28.9%), C (5.5%), and D (0.8%). Haplogroup B was predominant in the western part of Russia. A high level of mtDNA polymorphism in the studied groups of local sheep indicates the presence of a significant reserve of unique genotypes in Russia, which is to be explored.

Harvesting cereals at Tappeh Sang-e Chakhmaq and the introduction of farming in Northeastern Iran during the Neolithic

Tappeh Sang-e Chakhmaq is the only Neolithic site in Northeastern Iran, characterised by aceramic and ceramic levels corresponding to an occupation of 1500 years from the eighth to the end of the sixth millennium BCE. The Western and Eastern Mounds represent the oldest and longest occupation among the sites identified East of the Zagros, providing a unique context to explore the origin and spread of farming outside the core area of the Eastern Fertile Crescent. We present data about the first harvesting activities in the Northeastern Iranian Central Plateau by applying usewear and microtexture analysis through confocal microscopy on sickle gloss blades. Our results indicate a community of pioneer farmers who settled down in the area carrying with them both domestic cereals as well as advanced techniques of cereal cultivation. We demonstrate that most of the tools were used for harvesting cereals in a fully ripened state collected near the ground, indicating a well-established cereal cultivation strategy. The use of straight shafts with parallel inserts in Tappeh Sang-e Chakhmaq, as known in some sites in the Zagros, suggests the dispersal of farming practices and technologies from the Eastern Fertile Crescent north-eastward across Iran. We observe an evolution in the degree of ripeness of harvested cereals along the first four levels of occupation of the Western Mound, where semi-ripe harvesting is relatively important, suggesting that domestic cereals to be harvested before full maturity were introduced into the village. From the topmost of the Western Mound and along the occupation of the Eastern Mound, ripe harvesting is dominant, showing a well-established cultivation strategy of fully mature cereal. This shift could indicate an in-situ evolution towards a better-established agricultural technology, including harvesting riper crops, that would have resulted in higher yields, as cereals were collected when the grain was fully formed.

Genomic Landscape of Copy Number Variations and Their Associations with Climatic Variables in the World's Sheep

Sheep show characteristics of phenotypic diversity and adaptation to diverse climatic regions. Previous studies indicated associations between copy number variations (CNVs) and climate-driven adaptive evolution in humans and other domestic animals. Here, we constructed a genomic landscape of CNVs (n = 39,145) in 47 old autochthonous populations genotyped at a set of high-density (600 K) SNPs to detect environment-driven signatures of CNVs using a multivariate regression model. We found 136 deletions and 52 duplications that were significantly (Padj. < 0.05) associated with climatic variables. These climate-mediated selective CNVs are involved in functional candidate genes for heat stress and cold climate adaptation (e.g., B3GNTL1, UBE2L3, and TRAF2), coat and wool-related traits (e.g., TMEM9, STRA6, RASGRP2, and PLA2G3), repairing damaged DNA (e.g., HTT), GTPase activity (e.g., COPG), fast metabolism (e.g., LMF2 and LPIN3), fertility and reproduction (e.g., SLC19A1 and CCDC155), growth-related traits (e.g., ADRM1 and IGFALS), and immune response (e.g., BEGAIN and RNF121) in sheep. In particular, we identified significant (Padj. < 0.05) associations between probes in deleted/duplicated CNVs and solar radiation. Enrichment analysis of the gene sets among all the CNVs revealed significant (Padj. < 0.05) enriched gene ontology terms and pathways related to functions such as nucleotide, protein complex, and GTPase activity. Additionally, we observed overlapping between the CNVs and 140 known sheep QTLs. Our findings imply that CNVs can serve as genomic markers for the selection of sheep adapted to specific climatic conditions.

Genetic diversity and relationship among Indonesian local sheep breeds on Java Island based on mitochondrial cytochrome b gene sequences

BACKGROUND

The cytochrome b (Cyt b) gene is one of the most studied mitochondrial DNA (mtDNA) genes to determine sheep's genetic profile. This study aimed to determine the genetic diversity and relationships of several Indonesian local sheep populations on Java Island, Indonesia, based on the mtDNA Cyt b gene sequences. Blood samples were collected from forty-one individual sheep in seven populations of Indonesia local sheep breeds on Java Island (Priangan = 6, Garut = 6, Batur = 7, Wonosobo = 5, Javanese Thin-Tailed/JTT = 7, Javanese Fat-Tailed/JFT = 5, and Sapudi = 5). DNA extraction was performed on blood samples, and the mtDNA Cyt b gene was amplified using specific primers (Alek-CBF: 5'-CAACCCCACCACTTACAA-3' and Alek-CBR: 5'-CCTTGAGTCTTAGGGAGGTT-3'). The polymerase chain reaction (PCR) products were then sequenced, and data were analyzed using the MEGA version 7.0, DNA SP version 6.0, and NTSYS-pc version 2.11 software.

RESULTS

A total of 1140 bp complete mtDNA Cyt b gene sequences in this study obtained 1134 monomorphic sites (I), six polymorphic sites (V), one segregation site (S), and five parsimony informative sites (P) with a nucleotide diversity (Pi), the average number of nucleotide differences (K), and sequence conservation (SC) were 0.00119, 1.35610, and 0.9947, respectively. There were six haplotypes consisting of two unique haplotypes and four shared haplotypes with a haplotype diversity (Hd) of 0.5415. The genetic distance within and between populations ranged from 0.0000 to 0.0016 and 0.0000 to 0.0020, respectively. Wonosobo, JFT, and Sapudi sheep have the closest relationship, and then these three breeds were close to JTT sheep, followed by Batur, Priangan, and Garut sheep. Two haplogroups have been found based on the Ovine haplogroup clustering. All Wonosobo, JTT, JFT, Sapudi sheep, and most Batur sheep were clustered into haplogroup B. In contrast, Garut sheep were mostly clustered into haplogroup A, while Priangan sheep were clustered into both haplogroups with the same percentage.

CONCLUSION

Seven Indonesian local sheep breeds on Java Island have a close relationship clustered into two haplogroups, namely haplogroups A and B. Most Indonesian local sheep breeds on Java Island in this study were clustered into haplogroup B, except for Garut and Priangan sheep.

Genetic variation and demographic history of Sudan desert sheep reveal two diversified lineages

More than 400 million sheep are raised on the African continent, the majority of which are indigenous and are primarily reared for sustenance. They have effectively adapted to various climatic and production environments, surviving and flourishing. The genetic relationships among these sheep populations remain understudied. Herein, we sequenced the entire mitochondrial DNA control region of 120 animals from Hamary and Kabashi and their crossbreed (Hamary x Kabashi) of Sudan desert sheep (SDS) to understand their maternal-inherited genetic variation and demographic history profiles and relate those to the history of sheep pastoralism on the African continent. The results show a diversified and predominant D- loop haplogroup B (n = 102, 85%), with all other sequences belonging to haplogroup A. Most of the maternal genetic variation was partitioned between haplogroup (76.3%) while within haplogroup accounted for 23.7% of the variation. However, little genetic differentiation was observed among the two breeds and their crosses, with our results supporting a Hamari maternal origin for the crossbreed. Bayesian coalescent-based analysis reveals distinct demographic history between the two haplogroups, two breeds and their crosses. Comparison of the two haplogroup showed that haplogroup B experienced an earlier expansion than haplogroup A. Unlike the breed-based comparison, the expansion of the two breeds started roughly at the same time, around 6500 years ago, with Kabashi having a slightly greater effective population size. The maternal ancestors of SDS may have diverged before their introduction to the African continent. This study provides novel insights into the early history of these two main breeds of Sudan desert sheep and their crosses.

Shedding Light on the Origin of Egyptian Sheep Breeds by Evolutionary Comparison of Mitochondrial D-Loop

Simple Summary

Egypt is a carrefour between North African Maghreb countries, East Tropical Africa, Arabian Peninsula, and Near East countries including Turkey. The aim of the present work was to describe the genetic relationships between several domestic sheep populations of these regions. All of the present-day Egyptian breeds are coarse woolen and fat-tailed, but archeologists indicate that this tail form was acquired in the Near East later from populations with a thin tail. To test this idea, we used a phylogenetic-derived program to compare the control region of mtDNA of 37 breeds with fat or thin tails. We showed that most breeds seemed to fit with the archeologist hypothesis, whereas one breed indicates a direct migration of a fat-tailed breed from Turkey. Unexpectedly, one of the breeds from South Egypt was strongly linked to the thin-tailed desert breeds of Sudan, raising the question of the events leading to this situation.

Abstract

(1) Background: It has been recognized that the origin of fat-tailed sheep occurred within coarse wool breeds and that this character was introgressed several times into thin-tailed populations. However, no study has investigated this idea for Egyptian breeds using mtDNA analyses. (2) Methods: Using new sequences of the control region, we constructed a database of 467 sequences representing 37 breeds including fat- and thin-tailed ones with 80 Egyptian individuals belonging to six local breeds (Barki, Fallahi, Ossimi, Rahmani, Saidi, Sohagi). The phylogenetic tree obtained with the maximum likelihood method was submitted to the Newick Extra program to count the direct and indirect links between the individuals of each breed. (3) Results: Several Egyptian breeds were strongly connected to “primitive” thin-tailed breeds from Europe, indicating a clear genetic background of the “thin tail” breed type that supports the view of archeologists. In several cases, we suspected Western Asian breeds to be involved in the introgression of the fat tail character. In contrast, the Ossimi breed showed a high affinity to a fat-tailed breed of Western Asia, suggesting a direct migration and no thin tail ancestors. The Saidi is unique as our analyses revealed its strong connection with thin-tailed Sudanese breeds.

Historical Westward Migration Phases of Ovis aries Inferred from the Population Structure and the Phylogeography of Occidental Mediterranean Native Sheep Breeds

In this study, the genetic relationship and the population structure of western Mediterranean basin native sheep breeds are investigated, analyzing Maghrebian, Central Italian, and Venetian sheep with a highly informative microsatellite markers panel. The phylogeographical analysis, between breeds’ differentiation level (Wright’s fixation index), gene flow, ancestral relatedness measured by molecular coancestry, genetic distances, divergence times estimates and structure analyses, were revealed based on the assessment of 975 genotyped animals. The results unveiled the past introduction and migration history of sheep in the occidental Mediterranean basin since the early Neolithic. Our findings provided a scenario of three westward sheep migration phases fitting properly to the westward Neolithic expansion argued by zooarcheological, historical and human genetic studies.

Sheep Post-Domestication Expansion in the Context of Mitochondrial and Y Chromosome Haplogroups and Haplotypes

Mitochondrial DNA and nonrecombinant parts of Y-chromosome DNA are a great tool for looking at a species’ past. They are inherited for generations almost unaffected because they do not participate in recombination; thus, the time of occurrence of each mutation can be estimated based on the average mutation rate. Thanks to this, male and female haplogroups guide confirming events in the distant past (potential centers of domestication, settlement of areas, trade connections) as well as in modern breeding (crossbreeding, confirmation of paternity). This research focuses mainly on the development of domestic sheep and its post-domestication expansion, which has occurred through human trade from one continent to another. So far, five mitochondrial and five Y-chromosome haplogroups and dozens of their haplotypes have been detected in domestic sheep through studies worldwide. Mitochondrial DNA variability is more or less correlated with distance from the domestication center, but variability on the recombinant region of the Y chromosome is not. According to available data, central China shows the highest variability of male haplogroups and haplotypes.

Broad maternal geographic origin of domestic sheep in Anatolia and the Zagros

We investigated the controversial origin of domestic sheep (Ovis aries) using large samples of contemporary and ancient domestic individuals and their closest wild relatives: the Asiatic mouflon (Ovis gmelini), the urial (Ovis vignei) and the argali (Ovis ammon). A phylogeny based on mitochondrial DNA, including 213 new cytochrome-b sequences of wild Ovism confirmed that O. gmelini is the maternal ancestor of sheep and precluded mtDNA contributions from O. vignei (and O. gmelini × O. vignei hybrids) to domestic lineages. We also produced 54 new control region sequences showing shared haplogroups (A, B, C and E) between domestic sheep and wild O. gmelini which localized the domestication center in eastern Anatolia and central Zagros, excluding regions further east where exclusively wild haplogroups were found. This overlaps with the geographic distribution of O. gmelini gmelini, further suggesting that the maternal origin of domestic sheep derives from this subspecies. Additionally, we produced 57 new CR sequences of Neolithic sheep remains from a large area covering Anatolia to Europe, showing the early presence of at least three mitochondrial haplogroups (A, B and D) in Western colonization routes. This confirmed that sheep domestication was a large-scale process that captured diverse maternal lineages (haplogroups).

Maternal Relationships among Ancient and Modern Southern African Sheep: Newly Discovered Mitochondrial Haplogroups

Simple Summary

The genetic diversity of southern African sheep remains under-studied. We present here the complete mitochondrial genomes of archaeological southern African sheep, as well as the genomes from three indigenous southern African breeds—Damara, Namaqua Afrikaner, and Ronderib Afrikaner. We show that southern African sheep exhibit limited genetic diversity which is consistent with our understanding of their migration south from northernmost Africa. Intriguingly, many of the modern sheep show close relationships with the archaeological sheep, implying an ancestor-descendant relationship. Similarly, the sheep that do not exhibit a close relationship with the archaeological sheep nonetheless cluster closely with each other and do not show a close relationship with European and Asian sheep. This suggests that they too are descendants of indigenous sheep and not the product of historic introductions of exotic breeds.

Abstract

We investigated the genetic diversity and historic relationships among southern African sheep as well as the relationships between them and sheep outside the continent by sourcing both archaeological and modern sheep samples. Archaeological sheep samples derived from the site Die Kelders 1, near Cape Town, date to approximately 1500 years ago. The modern samples were taken as ear snips from Damara, Namaqua Afrikaner, and Ronderib Afrikaner sheep on a farm in Prieska in the Northern Cape. Illumina sequencing libraries were constructed for both ancient and modern specimens. Ancient specimens were enriched for the mitochondrial genome using an in-solution hybridization protocol and modern specimens were subjected to shotgun sequencing. Sequences were mapped to the Ovis aries reference genome, assigned to haplogroups and subhaplogroups, and used to calculate a phylogenetic tree using previously published, geographically dispersed mitochondrial genome sheep sequences. Genetic diversity statistics show that southern African sheep have lower diversity than sheep in other regions. Phylogenetic analysis reveals that many modern southern African sheep are likely descended from prehistoric indigenous sheep populations and not from sheep imported from Europe during the historic period.

Evaluation of Maternal Genetic Background of Two Hungarian Autochthonous Sheep Breeds Coming from Different Geographical Directions

The aim of our research was the evaluation of the maternal genetic background of two Hungarian autochthonous sheep breeds of different geographical origin. A major argument for the preservation of endangered animal breeds is their documented past and historical importance. These also include the registration of pedigree data. This is the first study to evaluate and compare Tsigai and Cikta sheep in Hungary. Our investigation is based on two complete sequences of mitochondrial DNA (cytochrome b gene and control region). Our research was performed on these two sheep breeds with markedly different breed histories and breed characteristics to determine a possible common maternal genetic background, as ultimately the origin of both breeds can be traced back to Asia Minor. Between 2015 and 2017, a total of 203 biological samples were taken using a newly introduced founder sampling method. We found that the prevailing haplogroup B accounted for over 80% of both breeds, strengthening the common ancestral root. However, the pairwise genetic differentiation estimates (KST) calculated using the sequence-based statistics for cytochrome b gene and control region were 0.034 and 0.021, respectively (both at level p < 0.05); thus, revealing genetic differentiation in both sequences between the Tsigai and Cikta. We note that the known different history of the breeds is clearly justified by the currently studied deviations in their maternal genetic background.

An ancient positively selected BMPRIB missense variant increases litter size of Mongolian sheep populations following latitudinal gradient

New gene mutation origination is a driving force for the evolution of organisms. The effect of FecB mutation in BMPRIB gene on the litter size of sheep has been well known for a long time, each copy of the mutant allele increases litter size by 0.4-0.5. However, the origin and adaptive evolution mechanism of FecB mutation are still unclear. Here we carried on the thorough analysis on evolutionary features of BMPRIB gene and found that 150 species as a whole is under purifying selection while sheep lineage shows evidence of positive selection. The results of allele age estimation revealed that the FecB mutation in Mongolian sheep of China originated in Mongolian Plateau at about 5000 years ago. Due the two shape drops in temperature subsequently, Mongolian sheep migrated from north to south following the northern nomadic people. Accordingly, the FecB mutant allele frequency increased, with the lowest in sheep locating at Mongolian plateau (0.01) and the highest in sheep locating at Yangtze River valley (0.96). In conclusion, the FecB mutation in Mongolian sheep of China originated in Mongolian Plateau at about 5000 years ago, and the differentiated litter size of Mongolian sheep might be the result of adaptation to various environments during the migration following latitudinal gradient. This study may well exemplify selection on an ancient variation triggered by drastic ecological shifts, and is also helpful to analyze the adaptive evolution mechanism of economic traits of domestic animals and identify major genes and molecular markers.

Whole-genome resequencing of worldwide wild and domestic sheep elucidates genetic diversity, introgression and agronomically important loci

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3′-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.

Archaeogenetic analysis of Neolithic sheep from Anatolia suggests a complex demographic history since domestication

Sheep were among the first domesticated animals, but their demographic history is little understood. Here we analyzed nuclear polymorphism and mitochondrial data (mtDNA) from ancient central and west Anatolian sheep dating from Epipaleolithic to late Neolithic, comparatively with modern-day breeds and central Asian Neolithic/Bronze Age sheep (OBI). Analyzing ancient nuclear data, we found that Anatolian Neolithic sheep (ANS) are genetically closest to present-day European breeds relative to Asian breeds, a conclusion supported by mtDNA haplogroup frequencies. In contrast, OBI showed higher genetic affinity to present-day Asian breeds. These results suggest that the east-west genetic structure observed in present-day breeds had already emerged by 6000 BCE, hinting at multiple sheep domestication episodes or early wild introgression in southwest Asia. Furthermore, we found that ANS are genetically distinct from all modern breeds. Our results suggest that European and Anatolian domestic sheep gene pools have been strongly remolded since the Neolithic.

Evidence for independent domestication of sheep mtDNA lineage A in India and introduction of lineage B through Arabian sea route

India ranks the second in the world in terms of its sheep population with approximately 74.26 million represented by 44 well-described breeds in addition to several non-descript populations. Genetic diversity and phylogeography of Indian sheep breeds remain poorly understood, particularly for south Indian breeds. To have a comprehensive view of the domestication history of Indian sheep, we sequenced the mitochondrial DNA (mtDNA) control region (D-loop) and cytochrome b gene (CYTB) of 16 Indian domestic sheep breeds, most of them (13) from the south India. We analysed these sequences along with published data of domestic and wild sheep from different countries, including India. The haplotype diversity was relatively high in Indian sheep, which were classified into the three known mtDNA lineages, namely A, B and C. Lineage A was predominant among Indian sheep whereas lineages B and C were observed at low frequencies but C was restricted to the breeds of north and east India. The median joining network showed five major expanding haplogroups of lineage A (A1–A5). Out of which, A2, A4 and A5 were more frequent in Indian sheep in contrast to breeds from other parts of the world. Among the 27 Indian sheep breeds analysed, Mandya and Sonadi breeds were significantly different from other Indian breeds in the MDS analyses. This was explained by a very high contribution of lineage B into these two breeds. The Approximate Bayesian Computation (ABC) provided evidence for the domestication of lineage A sheep in the Indian subcontinent. Contrary to the current knowledge, we also found strong support for the introduction of lineage B into Indian subcontinent through sea route rather than from the Mongolian Plateau. The neighbour-joining tree of domestic and wild sheep revealed the close genetic relationship of Indian domestic sheep with Pakistani wild sheep O. vignei blanfordi. Based on our analyses and archaeological evidences, we suggest the Indian subcontinent as one of the domestication centres of the lineage A sheep, while lineage B sheep might have arrived into India from elsewhere via Arabian sea route. To the best of our knowledge, this is the first comprehensive study on Indian sheep where we have analysed more than 740 animals belonging to 27 sheep breeds raised in various regions of India. Our study provides insight into the understanding of the origin and migratory history of Indian sheep.

Mitochondrial DNA analysis of ancient sheep from Kazakhstan: evidence for early sheep introduction

Kazakhstan covers a vast territory, and it has always been a land of nomadic pastoralism, where domesticated horses and sheep were moved by nomadic people across the steppe. Previous reports suggest that sheep breeds from Kazakhstan have an intermediate genetic composition between Asian and European breeds; however, this data appears to be limited. Therefore, we studied the genetic diversity of ancient domestic sheep from two Late Bronze Age settlements, Toksanbai and Kent, located in the Pre-Caspian region of Kazakhstan and central Kazakhstan, respectively. We have applied ZooMS analysis for taxonomic identification of small ruminant remains to select ancient specimens of domestic sheep (Ovis aries). To assign sheep mitochondrial DNA (mtDNA) haplogroups, the single nucleotide polymorphisms (SNPs) from the control region were analyzed by real-time PCR and direct sequencing. Identical distribution of mtDNA haplogroups A (8/14; 57%), B (5/14; 36%), and C (1/14; 7%) was observed in the specimens from Toksanbai (n = 14) and Kent (n = 14). Ovine haplogroup A was predominant in both settlements. Both archeological sites had similar patterns of haplogroup distribution, indicating early sheep introduction into the region. These results are important to gain a better understanding of sheep migrations in the Eurasian steppe and highlight the importance of genomic analysis of earlier local lineages.

Ancient Faunal History Revealed by Interdisciplinary Biomolecular Approaches

Starting four decades ago, studies have examined the ecology and evolutionary dynamics of populations and species using short mitochondrial DNA fragments and stable isotopes. Through technological and analytical advances, the methods and biomolecules at our disposal have increased significantly to now include lipids, whole genomes, proteomes, and even epigenomes. At an unprecedented resolution, the study of ancient biomolecules has made it possible for us to disentangle the complex processes that shaped the ancient faunal diversity across millennia, with the potential to aid in implicating probable causes of species extinction and how humans impacted the genetics and ecology of wild and domestic species. However, even now, few studies explore interdisciplinary biomolecular approaches to reveal ancient faunal diversity dynamics in relation to environmental and anthropogenic impact. This review will approach how biomolecules have been implemented in a broad variety of topics and species, from the extinct Pleistocene megafauna to ancient wild and domestic stocks, as well as how their future use has the potential to offer an enhanced understanding of drivers of past faunal diversity on Earth.

Meta-Analysis of Mitochondrial DNA Control Region Diversity to Shed Light on Phylogenetic Relationship and Demographic History of African Sheep (Ovis aries) Breeds

To improve sheep breeding and conservation of genetic resources, the mitochondrial DNA control region (mtDNA CR) of 399 sequences of African indigenous sheep breeds from previously published research articles were meta-analyzed to elucidate their phylogenetic relationship, diversity, and demographic history. A total of 272 haplotypes were found, of which 207 were unique and a high level of mtDNA CR variability was observed. Generally, the number of polymorphic sites, nucleotide and haplotype diversity were high (284, 0.254 ± 0.012 and 0.993 ± 0.002, respectively). The median-joining (MJ) network of haplotypes produced three major haplogroups (A, B and C), with haplogroup B being dominant. A mixture of populations suggests a common matrilineal origin and lack of and/or a weak phylogeographic structure. Mismatch analysis showed recent expansion of North African breeds, whereas East African and continental populations exhibited selection pressures for adaptation. A slight historical genetic difference was also observed between the fat tail and thin tail sheep breeds. However, further investigations are required using more samples and long sequence segments to achieve deeper levels of conclusions on the African sheep phylogenetic relationship. The present meta-analysis results contribute to the general understanding of African native sheep populations for improved management of sheep genetic resources.

Comparison of mtDNA control region among descendant breeds of the extinct Zaupel sheep revealed haplogroup C and D in Central Europe

The consideration of the descendance is indispensable in the preservation of endangered animal breeds. The authors compared mitochondrial DNA (mtDNA) control region sequence in three descendant breeds of the extinct Zaupel sheep, firstly. Their investigation was carried out in order to prove the common origin of Waldschaf (Austria), Bovec sheep (Slovenia) and Cikta (Hungary). A total of 118 biological samples were taken from non-related representatives of the three breeds between 2015 and 2017. A newly designed primer pair was also used to amplify the segment (1180 bp) to be tested. The total number of haplotypes in the whole study population was 49. The majority of which fell into haplogroup B. The significant negative value of the Fu's Fs statistic (Fs statistic = -3.296, p = 0.013) based on haplotype frequencies demonstrated a moderate foreign gene flow. As a novel observation haplogroups C and D appeared in Cikta and Bovec sheep, respectively. The Tajima D-test value in the entire study population was -0.914 (p > 0.10), meaning that the separation of the three descendant breeds did not cause genetic drift, these are collectively in genetic equilibrium. The genetic information confirmed the common origin of the breeds known from the breed history.

The genetic profiles and maternal origin of local sheep breeds on Java Island (Indonesia) based on complete mitochondrial DNA D-loop sequences

Background and Aim:

Java Island is one of the islands in Indonesia which has local sheep breeds with specific characteristics and native development geography in certain regions. This study aimed to determine the genetic profiles and maternal origin of six local sheep breeds on Java Island.

Materials and Methods:

This study was conducted by identifying the profiles of complete mitochondrial DNA (mtDNA) displacement loop (D-loop) region sequences on a total of 22 individual in six local sheep breeds on Java Island, including Javanese thin-tailed (JTT), Javanese Fat-Tailed (JFT), Batur (BTR), Wonosobo (WSB), Garut (GRT), and Priangan (PRG) sheep. The D-loop region was amplified using specific primers, and the polymerase chain reaction (PCR) was performed. The PCR products were purified and sequenced.

Results:

The mtDNA D-loop analysis identified 21 haplotypes in the analyzed 22 animals with 123 polymorphic sites (V) consisting of 60 singleton variable sites (S) and 63 parsimony informative sites (P). Within all breeds tested, the haplotype diversity, the average number of pairwise differences (K), and nucleotide diversity (Pi) were 0.99567, 25.36364, and 0.02153, respectively. The genetic distance (D) within groups and between groups was 0.001-0.006 and 0.004-0.036, respectively. The phylogeny resulted in the presence of two haplogroups (Hap), which are 5 Hap A and 16 Hap B. All JTT, JFT, BTR, and WSB breeds were in the same cluster in Hap B, whereas GRT and PRG breeds were in clusters in both Hap A and Hap B.

Conclusion:

The high genetic diversity in six local sheep breeds on Java Island suggests that they originated from different genetic sources. JTT sheep have closer genetic relationships to JFT, BTR, and WSB sheep, and they are close to European sheep, whereas GRT sheep have closer genetic relationships to PRG sheep. Both are closer to Asian sheep than to European sheep.

Microsatellite and mitochondrial DNA analyses unveil the genetic structure of native sheep breeds from three major agro-ecological regions of India

Sheep farming has been fundamental to many civilizations in the world and is practiced in India since antiquity. Several thousand years of adaptation to local environmental conditions and selective breeding have evolved 44 sheep breeds in India. They are paramount in terms of economic, scientific, and cultural heritage. Genetic characterization information is imperative for sustainable utilization and conservation of ovine heritage. In this study, the genetic diversity, differentiation, and structure of 11 indigenous sheep breeds from three different agro-ecological zones of India were explored with genomic microsatellite loci and mitochondrial DNA (D loop). The estimated diversity parameters indicated that populations retained high levels of genetic diversity (Na = 8.27 ± 0.17; Ho = 0.65 ± 0.01), which provides an optimistic viewpoint for their survival. However, significant inbreeding was also observed in the nine populations. Moderate genetic differentiation existed among the groups (F_ST = 0.129 ± 0.012), and most likely clusters existing in the dataset are seven. Phylogenetic clustering was in line with the geographical locations of sheep populations. Mitochondrial sequences revealed high haplotype diversity with the existence of maternal haplogroups A, B, and C, and signals of population expansion. Decreased genetic diversity and unique maternal lineage (C) in endangered Tibetan and Bonpala sheep breed, warrant their immediate scientific management. Overall, the quantitative data reported here on the extant variability, and genetic relationships among the Indian sheep breeds, provide critically important inputs that will be valuable for the decision-making process on their management, both for the conservation of endangered breeds, and formulation of breeding programs to check genetic erosion.

A Hu sheep genome with the first ovine Y chromosome reveal introgression history after sheep domestication

The Y chromosome plays key roles in male fertility and reflects the evolutionary history of paternal lineages. Here, we present a de novo genome assembly of the Hu sheep with the first draft assembly of ovine Y chromosome (oMSY), using nanopore sequencing and Hi-C technologies. The oMSY that we generated spans 10.6 Mb from which 775 Y-SNPs were identified by applying a large panel of whole genome sequences from worldwide sheep and wild Iranian mouflons. Three major paternal lineages (HY1a, HY1b and HY2) were defined across domestic sheep, of which HY2 was newly detected. Surprisingly, HY2 forms a monophyletic clade with the Iranian mouflons and is highly divergent from both HY1a and HY1b. Demographic analysis of Y chromosomes, mitochondrial and nuclear genomes confirmed that HY2 and the maternal counterpart of lineage C represented a distinct wild mouflon population in Iran that diverge from the direct ancestor of domestic sheep, the wild mouflons in Southeastern Anatolia. Our results suggest that wild Iranian mouflons had introgressed into domestic sheep and thereby introduced this Iranian mouflon specific lineage carrying HY2 to both East Asian and Africa sheep populations.

Mitochondrial DNA Variations in Colombian Creole Sheep Confirm an Iberian Origin and Shed Light on the Dynamics of Introduction Events of African Genotypes

The genetic origins and diversity of Creole sheep from five regions of Colombia were investigated based on mitochondrial DNA (mtDNA) variations across 89 sequences from five breeds: one wool Creole sheep (CL) and four hair Creole sheep, including Ethiopian (OPCE), Sudan (OPCS), Pelibuey (OPCP) and Wayúu (OPCW). A global comparison was done using 62 haplotypes from Iberian, African, Indian, Caribbean, Mexican, Caucasian and European sheep based on sequences retrieved from GenBank. This study aimed to identify the maternal origin of Colombian Creole sheep and their genetic relationships at a global level. The results showed 31 different haplotypes from Colombian Creole sheep, which can be assigned to maternal lineage B, the most common lineage found in European sheep breeds and the only one found in several Iberian breed (e.g., Churra, Spanish Merino) that most likely participated in the Creole formation. Additional analyses showed that wool and hair sheep retained a broad genetic identity despite being geographically separated. The global-level phylogenetic analysis revealed that Colombian Creole sheep belong to a distinct and defined genetic lineage that is likely the result of a founder effect with ecotypes of Iberian descent and the subsequent introduction of foreign breeds. This is consistent with historical reports on the presence of sheep in South America and, particularly, Colombia.

Phylogeography and Phylogenetic Evolution in Tibetan Sheep Based on MT-CYB Sequences

Simple Summary

The molecular and population genetic evidence of the phylogenetic status of the Tibetan sheep (Ovis aries) is not well understood, and little is known about this species’ genetic diversity. The aim of the present research was to explore phylogeography and phylogenetic evolution of Tibetan sheep populations, on the basis of mitochondrial DNA (mtDNA) gene MT-CYB (1140 base pairs). The results of this study will add information to the Tibetan sheep populations and help enlighten upcoming programs related to conservation of Tibetan sheep living in the Qinghai–Tibetan Plateau.

Abstract

To date, molecular genetics and population studies in Tibetan sheep (Ovis aries) have been limited, and little is known about the phylogenetic evolution and phylogeography of Tibetan sheep populations. The aim of the present research was to explore phylogeography and phylogenetic evolution of Tibetan sheep populations, on the basis of mitochondrial DNA (mtDNA) gene MT-CYB (1140 base pairs). Our dataset consisted of 641 MT-CYB sequences from the same amount of animals belonging to 15 populations of Tibetan sheep living in the Qinghai–Tibetan Plateau, China. Haplotype and nucleotide diversities were 0.748 ± 0.010 and 0.003 ± 0.001, respectively. The analysis of phylogeography revealed the presence of two formerly described haplogroups in 15 populations of Tibetan sheep, however only one haplogroup was present in Awang sheep. Moreover, 641 Tibetan sheep were distributed into a minimum of two clusters by clustering analysis. The 15 Tibetan sheep populations and 19 reference populations of 878 individuals were separated into six main groups based on their substitutions per site, from which we constructed a phylogenetic tree. Minor differences in branching order of various taxa between trees acquired from either gene were observed. This study provides insights on the origins and phylogenetic evolution of populations residing in the Qinghai–Tibetan Plateau, which will aid information of future conservation programs aimed at conserving this valuable genetic resource.

The origin of sheep settlement in Western Mediterranean

The arrival of Neolithic culture in North Africa, especially domestic animals has been essentially documented from archaeological records. As the data relative to sheep are scarce, we studied the genetic relationship between Moroccan sheep breeds and Mediterranean ones using the sequencing of 628 bp of the mitochondrial DNA control region in 193 Moroccan individuals, belonging to six breeds, and 652 sequences from other breeds in Europe and Middle East. Through Network analysis and an original phylogenetically derived method, the connection proportions of each Moroccan breed to foreign ones were estimated, highlighting the strong links between Moroccan and Iberian breeds. The first founders of Moroccan sheep population were issued at 79% from Iberia and 21% from a territory between Middle East and Africa. Their calculated expansion times were respectively 7,100 and 8,600 years B.P. This suggests that Neolithization was introduced by a double influence, from Iberia and from another route, maybe Oriental or Sub-Saharan. The consequence of the environmental changes encountered by founders from Iberia was tested using different neutrality tests. There are significant selection signatures at the level of Moroccan and European breeds settled in elevated altitudes, and an erosion of nucleotide diversity in Moroccan breeds living in arid areas.

Evolution of the sheep coat: the impact of domestication on its structure and development

Wild sheep and many primitive domesticated breeds have two coats: coarse hairs covering shorter, finer fibres. Both are shed annually. Exploitation of wool for apparel in the Bronze Age encouraged breeding for denser fleeces and continuously growing white fibres. The Merino is regarded as the culmination of this process. Archaeological discoveries, ancient images and parchment records portray this as an evolutionary progression, spanning millennia. However, examination of the fleeces from feral, two-coated and woolled sheep has revealed a ready facility of the follicle population to change from shedding to continuous growth and to revert from domesticated to primitive states. Modifications to coat structure, colour and composition have occurred in timeframes and to sheep population sizes that exclude the likelihood of variations arising from mutations and natural selection. The features are characteristic of the domestication phenotype: an assemblage of developmental, physiological, skeletal and hormonal modifications common to a wide variety of species under human control. The phenotypic similarities appeared to result from an accumulation of cryptic genetic changes early during vertebrate evolution. Because they did not affect fitness in the wild, the mutations were protected from adverse selection, becoming apparent only after exposure to a domestic environment. The neural crest, a transient embryonic cell population unique to vertebrates, has been implicated in the manifestations of the domesticated phenotype. This hypothesis is discussed with reference to the development of the wool follicle population and the particular roles of Notch pathway genes, culminating in the specific cell interactions that typify follicle initiation.

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.

Identification of an ancestral haplotype in the mitochondrial phylogeny of the ovine haplogroup B

BACKGROUND

European mouflon (Ovis orientalis musimon) has been reintroduced in mainland Europe since the 18th-century sourcing from the Sardinian and Corsican autochthonous mouflon populations. The European mouflon is currently considered the feral descendent of the Asian mouflon (O. orientalis), and the result of first wave of sheep domestication occurred 11,000 years ago in the Fertile Crescent, and brought to Corsica and Sardinia ca. 6,000 years ago, where they still live as autochthonous populations. However, this phylogeny is based on mitogenome sequences of European mouflon individuals exclusively.

METHODS

We sequenced the first complete mtDNA of the long-time isolated Sardinian mouflon and compared it with several ovine homologous sequences, including mouflon from mainland Europe and samples representative of the five known mitochondrial domestic sheep haplogroups. We applied Bayesian inference, Maximum Likelihood and Integer Neighbour-Joining network methods and provided a robust, fully-resolved phylogeny with strong statistical support for all nodes.

RESULTS

We identified an early split (110,000 years ago) of the Sardinian mouflon haplotype from both sheep and mainland European mouflon belonging to haplogroup B, the latter two sharing a more recent common ancestor (80,000 years ago). Further, the Sardinian mouflon sequence we generated had the largest genetic distance from domestic sheep haplogroups (0.0136 ± 0.004) among mouflon species. Our results suggest the Sardinian mouflon haplotype as the most ancestral in the HPG-B lineage, hence partially redrawing the known phylogeny of the genus Ovis.

The global diversity of Haemonchus contortus is shaped by human intervention and climate

Haemonchus contortus is a haematophagous parasitic nematode of veterinary interest. We have performed a survey of its genome-wide diversity using single-worm whole genome sequencing of 223 individuals sampled from 19 isolates spanning five continents. We find an African origin for the species, together with evidence for parasites spreading during the transatlantic slave trade and colonisation of Australia. Strong selective sweeps surrounding the β-tubulin locus, a target of benzimidazole anthelmintic drug, are identified in independent populations. These sweeps are further supported by signals of diversifying selection enriched in genes involved in response to drugs and other anthelmintic-associated biological functions. We also identify some candidate genes that may play a role in ivermectin resistance. Finally, genetic signatures of climate-driven adaptation are described, revealing a gene acting as an epigenetic regulator and components of the dauer pathway. These results begin to define genetic adaptation to climate in a parasitic nematode.

Based on single worm whole genome sequencing, the authors here characterise the global evolution of the gastrointestinal parasite Haemonchus contortus and identify genes that play a role in drug resistance as well as climate-driven adaptations involving an epigenetic regulator.

Genetic analysis identifies the missing parchment of New Zealand's founding document, the Treaty of Waitangi

Genetic analyses provide a powerful tool with which to identify the biological components of historical objects. Te Tiriti o Waitangi | The Treaty of Waitangi is New Zealand’s founding document, intended to be a partnership between the indigenous Māori and the British Crown. Here we focus on an archived piece of blank parchment that has been proposed to be the missing portion of the lower parchment of the Waitangi Sheet of the Treaty. However, its physical dimensions and characteristics are not consistent with this hypothesis. We perform genetic analyses on the parchment membranes of the Treaty, plus the blank piece of parchment. We find that all three parchments were made from ewes and that the blank parchment is highly likely to be a portion cut from the lower membrane of the Waitangi Sheet because they share identical whole mitochondrial genomes, including an unusual heteroplasmic site. We suggest that the differences in size and characteristics between the two pieces of parchment may have resulted from the Treaty’s exposure to water in the early 20^th century and the subsequent repair work, light exposure during exhibition or the later conservation treatments in the 1970s and 80s. The blank piece of parchment will be valuable for comparison tests to study the effects of earlier treatments and to monitor the effects of long-term display on the Treaty.

Population structure and genetic diversity of 25 Russian sheep breeds based on whole-genome genotyping

Background

Russia has a diverse variety of native and locally developed sheep breeds with coarse, fine, and semi-fine wool, which inhabit different climate zones and landscapes that range from hot deserts to harsh northern areas. To date, no genome-wide information has been used to investigate the history and genetic characteristics of the extant local Russian sheep populations. To infer the population structure and genome-wide diversity of Russian sheep, 25 local breeds were genotyped with the OvineSNP50 BeadChip. Furthermore, to evaluate admixture contributions from foreign breeds in Russian sheep, a set of 58 worldwide breeds from publicly available genotypes was added to our data.

Results

We recorded similar observed heterozygosity (0.354–0.395) and allelic richness (1.890–1.955) levels across the analyzed breeds and they are comparable with those observed in the worldwide breeds. Recent effective population sizes estimated from linkage disequilibrium five generations ago ranged from 65 to 543. Multi-dimensional scaling, admixture, and neighbor-net analyses consistently identified a two-step subdivision of the Russian local sheep breeds. A first split clustered the Russian sheep populations according to their wool type (fine wool, semi-fine wool and coarse wool). The Dagestan Mountain and Baikal fine-fleeced breeds differ from the other Merino-derived local breeds. The semi-fine wool cluster combined a breed of Romanian origin, Tsigai, with its derivative Altai Mountain, the two Romney-introgressed breeds Kuibyshev and North Caucasian, and the Lincoln-introgressed Russian longhaired breed. The coarse-wool group comprised the Nordic short-tailed Romanov, the long-fat-tailed outlier Kuchugur and two clusters of fat-tailed sheep: the Caucasian Mountain breeds and the Buubei, Karakul, Edilbai, Kalmyk and Tuva breeds. The Russian fat-tailed breeds shared co-ancestry with sheep from China and Southwestern Asia (Iran).

Conclusions

In this study, we derived the genetic characteristics of the major Russian local sheep breeds, which are moderately diverse and have a strong population structure. Pooling our data with a worldwide genotyping set gave deeper insight into the history and origin of the Russian sheep populations.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0399-5) contains supplementary material, which is available to authorized users.