Tracing the History of Goat Pastoralism: New Clues from Mitochondrial and Y Chromosome DNA in North Africa

Anthropogenic events and responses to environmental stress are shaping the genomes of Ethiopian indigenous goats

Anthropological and biophysical processes have shaped livestock genomes over Millenia and can explain their current geographic distribution and genetic divergence. We analyzed 57 Ethiopian indigenous domestic goat genomes alongside 67 equivalents of east, west, and north-west African, European, South Asian, Middle East, and wild Bezoar goats. Cluster, ADMIXTURE (K = 4) and phylogenetic analysis revealed four genetic groups comprising African, European, South Asian, and wild Bezoar goats. The Middle Eastern goats had an admixed genome of these four genetic groups. At K = 5, the West African Dwarf and Moroccan goats were separated from East African goats demonstrating a likely historical legacy of goat arrival and dispersal into Africa via the coastal Mediterranean Sea and the Horn of Africa. FST, XP-EHH, and Hp analysis revealed signatures of selection in Ethiopian goats overlaying genes for thermo-sensitivity, oxidative stress response, high-altitude hypoxic adaptation, reproductive fitness, pathogen defence, immunity, pigmentation, DNA repair, modulation of renal function and integrated fluid and electrolyte homeostasis. Notable examples include TRPV1 (a nociception gene); PTPMT1 (a critical hypoxia survival gene); RETREG (a regulator of reticulophagy during starvation), and WNK4 (a molecular switch for osmoregulation). These results suggest that human-mediated translocations and adaptation to contrasting environments are shaping indigenous African goat genomes.

Characterization of heterozygosity-rich regions in Italian and worldwide goat breeds

Heterozygosity-rich regions (HRR) are genomic regions of high heterozygosity, which may harbor loci related to key functional traits such as immune response, survival rate, fertility, and other fitness traits. This study considered 30 Italian and 19 worldwide goat breeds genotyped with the Illumina GoatSNP50k BeadChip. The aim of the work was to study inter-breed relationships and HRR patterns using Sliding Window (SW) and Consecutive Runs (CR) detection methods. Genetic relationships highlighted a clear separation between non-European and European breeds, as well as the north-south geographic cline within the latter. The Pearson correlation coefficients between the descriptive HRR parameters obtained with the SW and CR methods were higher than 0.9. A total of 166 HRR islands were detected. CHI1, CHI11, CHI12 and CHI18 were the chromosomes harboring the highest number of HRR islands. The genes annotated in the islands were linked to various factors such as productive, reproductive, immune, and environmental adaptation mechanisms. Notably, the Montecristo feral goat showed the highest number of HRR islands despite the high level of inbreeding, underlining potential balancing selection events characterizing its evolutionary history. Identifying a species-specific HRR pattern could provide a clearer view of the mechanisms regulating the genome modelling following anthropogenic selection combined with environmental interaction.

Whole-genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats.

Multiple genomic solutions for local adaptation in two closely related species (sheep and goats) facing the same climatic constraints

The question of how local adaptation takes place remains a fundamental question in evolutionary biology. The variation of allele frequencies in genes under selection over environmental gradients remains mainly theoretical and its empirical assessment would help understanding how adaptation happens over environmental clines. To bring new insights to this issue we set up a broad framework which aimed to compare the adaptive trajectories over environmental clines in two domesticated mammal species co-distributed in diversified landscapes. We sequenced the genomes of 160 sheep and 161 goats extensively managed along environmental gradients, including temperature, rainfall, seasonality and altitude, to identify genes and biological processes shaping local adaptation. Allele frequencies at putatively adaptive loci were rarely found to vary gradually along environmental gradients, but rather displayed a discontinuous shift at the extremities of environmental clines. Of the 430 candidate adaptive genes identified, only 6 were orthologous between sheep and goats and those responded differently to environmental pressures, suggesting different putative mechanisms involved in local adaptation in these two closely related species. Interestingly, the genomes of the 2 species were impacted differently by the environment, genes related to signatures of selection were most related to altitude, slope and rainfall seasonality for sheep, and summer temperature and spring rainfall for goats. The diversity of candidate adaptive pathways may result from a high number of biological functions involved in the adaptations to multiple eco-climatic gradients, and a differential role of climatic drivers on the two species, despite their co-distribution along the same environmental gradients. This study describes empirical examples of clinal variation in putatively adaptive alleles with different patterns in allele frequency distributions over continuous environmental gradients, thus showing the diversity of genetic responses in adaptive landscapes and opening new horizons for understanding genomics of adaptation in mammalian species and beyond.

Refining the Global Phylogeny of Mitochondrial N1a, X, and HV2 Haplogroups Based on Rare Mitogenomes from Croatian Isolates

Mitochondrial DNA (mtDNA) has been used for decades as a predominant tool in population genetics and as a valuable addition to forensic genetic research, owing to its unique maternal inheritance pattern that enables the tracing of individuals along the maternal lineage across numerous generations. The dynamic interplay between evolutionary forces, primarily genetic drift, bottlenecks, and the founder effect, can exert significant influence on genetic profiles. Consequently, the Adriatic islands have accumulated a subset of lineages that exhibits remarkable absence or rarity within other European populations. This distinctive genetic composition underscores the islands' potential as a significant resource in phylogenetic research, with implications reaching beyond regional boundaries to contribute to a global understanding. In the initial attempt to expand the mitochondrial forensic database of the Croatian population with haplotypes from small isolated communities, we sequenced mitogenomes of rare haplogroups from different Croatian island and mainland populations using next-generation sequencing (NGS). In the next step and based on the obtained results, we refined the global phylogeny of haplogroup N1a, HV2, and X by analyzing rare haplotypes, which are absent from the current phylogenetic tree. The trees were based on 16 novel and 52 previously published samples, revealing completely novel branches in the X and HV2 haplogroups and a new European cluster in the ancestral N1a variant, previously believed to be an exclusively African-Asian haplogroup. The research emphasizes the importance of investigating geographically isolated populations and their unique characteristics within a global context.

Maternal genetic diversity of ancient goats in Finland and Estonia and comparison with extant northern European goat breeds

Countries in the northern Baltic Sea region have been historically linked for thousands of years, and animal husbandry is one of the many information sources that enables the study of historical connections. Although goat husbandry in this part of Europe has been evidenced by scant archaeological materials, zooarchaeological and historical evidence has revealed its continuity, at least since the Late Iron Age. To explore the historical relationship between goat lineages and investigate affinities between the past and present-day populations in the Baltic Sea region, we analysed a 476-bp fragment of the mitochondrial DNA control region in 14 ancient goats from Finland and Estonia and 10 extant goats from Finland. The results revealed high mitochondrial diversity among the ancient goats. Two maternal lineages were shared between the Late Iron Age and medieval individuals from Finland and Estonia. Moreover, ancient Finnish and Estonian goats showed maternal affinity to extant Finngoat and Swedish Landrace breeds. Overall, the analysis of maternal goat lineages confirmed tight historical connections in the region.

Whole-genome sequencing and genomic analysis of Norduz goat (Capra hircus)

Artificial and natural selective breeding of goats has resulted in many different goat breeds all around the world. Norduz goat is one of these breeds, and it is a local goat breed of Turkey. The goats are favorable due to pre-weaning viability and reproduction values compared to the regional breeds. Development in sequencing technologies has let to understand huge genomic structures and complex phenotypes. Until now, such a comprehensive study has not been carried out to understand the genomic structure of the Norduz goats, yet. In the study, the next-generation sequencing was carried out to understand the genomic structure of Norduz goat. Real-time PCR was used to evaluate prominent CNVs in the Norduz goat individuals. Whole genome of the goat was constructed with an average of 33.1X coverage level. In the stringent filtering condition, 9,757,980 SNPs, 1,536,715 InDels, and 290 CNVs were detected in the Norduz goat genome. Functional analysis of high-impact SNP variations showed that the classical complement activation biological process was affected significantly in the goat. CNVs in the goat genome were found in genes related to defense against viruses, immune response, and cell membrane transporters. It was shown that GBP2, GBP5, and mammalian ortholog GBP1, which are INF-stimulated GTPases, were found to be high copy numbers in the goats. To conclude, genetic variations mainly in immunological response processes suggest that Norduz goat is an immunologically improved goat breed and natural selection could take an important role in the genetical improvements of the goats.

Peeping into Mitochondrial Diversity of Andaman Goats: Unveils Possibility of Maritime Transport with Diversified Geographic Signaling

Andaman and Nicobar Islands, a part of South-East Asia, is enriched with the presence of native breeds of livestock (cattle, pig, goat) and poultry. There are two native goat breeds, viz., Andaman local goat and Teressa goat in Andaman and Nicobar Islands. However, to date, the origin and genetic makeup of these two breeds have not been detailed. Therefore, the present study describes the genetic makeup of Andaman goats through analysis of mitochondrial D-loop sequence for sequence polymorphism, phylogeographical signaling and population expansion events. The genetic diversity of the Teressa goat was less compared to the Andaman local goat due to its sole presence on Teressa Island. Out of 38 well-defined haplotypes of Andaman goats, the majority of haplotypes belonged to haplogroup A followed by haplogroup B and haplogroup D. The result of mismatch distribution and neutrality tests indicated no population expansion event of haplogroup A and B. Finally, based on poor geographical signaling, we hypothesize that Andaman goats have been imported to these Islands either through multidirectional diffusion or unidirectional diffusion. We justify our hypothesis of multidirectional diffusion on the basis of observation of the haplotype and nucleotide diversity of Andaman goats. Simultaneously, the probability of unidirectional diffusion of goats in these islands from the Indian subcontinent in different spells of domestication events through maritime routes cannot be ignored.

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.

Transcriptomics and Selection Pressure Analysis Reveals the Influence Mechanism of PLIN1 Protein on the Development of Small Size in Min Pigs

Body size is an important biological phenotypic trait that has attracted substantial attention. Small domestic pigs can serve as excellent animal models for biomedicine and also help meet sacrificial culture needs in human societies. Although the mechanisms underlying vertebral development regulating body size variation in domestic pigs during the embryonic period have been well described, few studies have examined the genetic basis of body size variation in post embryonic developmental stages. In this study, seven candidate genes-PLIN1, LIPE, PNPLA1, SCD, FABP5, KRT10 and IVL-significantly associated with body size were identified in Min pigs, on the basis of weighted gene co-expression network analysis (WGCNA), and most of their functions were found to be associated with lipid deposition. Six candidate genes except for IVL were found to have been subjected to purifying selection. PLIN1 had the lowest ω value (0.139) and showed heterogeneous selective pressure among domestic pig lineages with different body sizes (p < 0.05). These results suggested that PLIN1 is an important genetic factor regulating lipid deposition and consequently affecting body size variation in pigs. The culture of whole pig sacrifice in Manchu during the Qing Dynasty in China might have contributed to the strong artificial domestication and selection of Hebao pigs.

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.

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY, DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes Y1A, Y1B, Y2A and Y2B with a marked geographic partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and 7 wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during the migrations into northern Europe, eastern and southern Asia and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19^th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographic range.

How have sheep breeds differentiated from each other in Morocco? Genetic structure and geographical distribution patterns

Background

Based on the relatively homogeneous origin of the sheep breeds in Morocco that originate mainly from Iberia, it is highly relevant to address the question of how these very diverse sheep populations differentiated from each other. The Mountains of the High Atlas and Middle Atlas are expected to constitute North–South and West–East geographical barriers, respectively, which could have shaped the history of the differentiation of sheep breeds. The aim of this study was to test this hypothesis by considering the genetic structure and the spatial distribution of five major breeds (Sardi, Timahdite, Beni Guil, Boujaad and D’man) and one minor breed (Blanche de Montagne), by analysing the mtDNA control region, using 30 individuals per breed.

Results

Phylogenetic and network analyses did not indicate any clear separation among the studied breeds and discriminant component principal analysis showed some overlap between them, which indicates a common genetic background. The calculated pairwise F_ST values and Nei’s genetic distances revealed that most breeds showed a moderate genetic differentiation. The lowest and highest degrees of differentiation were retrieved in the Beni Guil and Boujaad breeds, respectively. Analysis of molecular variance (AMOVA) indicated that more than 95% of the genetic diversity occurs within individuals, while between- and within-population variabilities represent only 1.332% and 2.881%, respectively. Isolation-by-distance, spatial Principal Component Analysis (sPCA), and spatial AMOVA analyses evidenced clear examples of geographical structuration among the breeds, both between and within breeds. However, several enigmatic relationships remain, which suggest the occurrence of complex events leading to breed differentiation.

Conclusions

The approaches used here resulted in a convergent view on the hypothetic events that could have led to the progressive differentiation between the Moroccan breeds. The major split seems to be linked to the West–East barrier of the Middle Atlas, whereas the influence of the High Atlas is less obvious and incompletely resolved. The study of additional breeds that have settled near the High Atlas should clarify the relationships between the breeds of the West part of the country, in spite of their small population size.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00679-2.

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.

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.

Characterizing the Mitochondrial Diversity of Arbi Goats from Tunisia

Arbi is one of the main local goat breeds in Tunisia, representing an important economic resource in arid and hot areas where cattle and sheep cannot thrive successfully. In the current work, we have characterized the mitochondrial diversity of 26 Arbi goats by partially sequencing the mitochondrial D-loop region. These sequences plus 10 retrieved from GenBank were analyzed with the DnaSP v.5.10.1, evidencing the existence of 12 different haplotypes. Nucleotide and haplotype diversities were 0.02 and 0.96. Moreover, median-joining network analysis showed that all D-loop sequences from Arbi goats correspond to haplogroup A and that in general they do not cluster with sequences from other goat breeds. The high diversity that has been observed in North African goats is compatible with the maritime diffusion of the Neolithic package 10,000-7000 YBP. Moreover, there are evidences that local Tunisian breeds have been extensively crossed with highly productive transboundary breeds in order to improve meat and milk yields. These uncontrolled crossing practices may lead to the loss of alleles that play key roles in the adaptation of Tunisian local breeds to a harsh environment.

Genetic Variation in Y-Chromosome Genes of Sika Deer (Cervus nippon) in Japan

Sika deer (Cervus nippon) in Japan are classified into southern and northern groups. However, previous studies primarily relied on maternally inherited mitochondrial DNA (mtDNA). The paternally inherited Y-chromosome is useful for analyzing the contribution of males to the population genetic history of sika deer. In total, approximately 16 kb of partial sequences of four Y-chromosomal genes, Y-linked, sex-determining region Y, DEAD-box helicase 3 Y-linked, and Zinc finger protein Y-linked, were sequenced to investigate intraspecific variation. As a result, we identified nine intronic single nucleotide polymorphisms (SNPs) in 478 sika deer samples collected over the entire Japanese archipelago from Hokkaido to Kyushu. SNP genotyping revealed 10 distinct haplotypes (SYH1-SYH10). The most common haplotype (SYH1) was present in all populations and was the most abundant haplotype, identified in 80.3% of the sampled individuals. The remaining haplotypes were unique to a single locality. SYH1 was also central to all other haplotypes that diverged by a SNP, resulting in this haplotype being the core of a star-like cluster topography. We found that contrary to mtDNA patterns, there was no clear differentiation of Y-chromosome markers between the southern and the northern populations. Due to the female philopatry of sika deer, mtDNA may provide a highly structured differentiation of populations. On the other hand, the male-biased gene flow may provide a reduced differentiation of populations. Our findings revealed that the genetic structure of the Japanese sika deer is more complex than previously thought based on mtDNA-based phylogeographic studies.

Diversity Analysis and Genetic Relationships among Local Brazilian Goat Breeds Using SSR Markers

Simple Summary

This study aimed to evaluate the genetic diversity of six groups of native Brazilian goats using a panel of single sequence repeats (SSRs). Results indicated a definite genetic differentiation among the Brazilian goat herd, which indicates the existence of at least four breeds according to the international concepts (Moxotó and Repartida; the Grauna and Serrana Azul; Canindé and Marota breeds).

Abstract

The genetic diversity of six Brazilian native goats was reported using molecular markers. Hair samples of 332 animals were collected from different goat breeds (Moxotó, Canindé, Serrana Azul, Marota, Repartida, and Graúna) from five states of Northeast Brazil (Paraíba, Pernambuco, Rio Grande do Norte, Bahia, and Piauí). A panel of 27 microsatellites or single sequence repeats (SSRs) were selected and amplified using a polymerase chain reaction (PCR) technique. All populations showed an average allele number of over six. The mean observed heterozygosity for Brazilian breeds was superior to 0.50. These results demonstrated the high genetic diversity in the studied populations with values ranging from 0.53 (Serrana Azul) to 0.62 (Repartida). The expected average heterozygosity followed the same trend ranging from 0.58 (Serrana Azul) to 0.65 (Repartida), and the values obtained are very similar for all six breeds. The fixation index (Fis) had values under 10% except for the Moxotó breed (13%). The mean expected heterozygosity of all Brazilian populations was over 0.50. Results indicated a within-breed genetic variability in the Brazilian breeds based on the average number of alleles and the average observed heterozygosity. The interbreed genetic diversity values showed proper genetic differentiation among local Brazilian goat breeds.

Paternal Origins and Migratory Episodes of Domestic Sheep

The domestication and subsequent global dispersal of livestock are crucial events in human history, but the migratory episodes during the history of livestock remain poorly documented [1-3]. Here, we first developed a set of 493 novel ovine SNPs of the male-specific region of Y chromosome (MSY) by genome mapping. We then conducted a comprehensive genomic analysis of Y chromosome, mitochondrial DNA, and whole-genome sequence variations in a large number of 595 rams representing 118 domestic populations across the world. We detected four different paternal lineages of domestic sheep and resolved, at the global level, their paternal origins and differentiation. In Northern European breeds, several of which have retained primitive traits (e.g., a small body size and short or thin tails), and fat-tailed sheep, we found an overrepresentation of MSY lineages y-HC and y-HB, respectively. Using an approximate Bayesian computation approach, we reconstruct the demographic expansions associated with the segregation of primitive and fat-tailed phenotypes. These results together with archaeological evidence and historical data suggested the first expansion of early domestic hair sheep and the later expansion of fat-tailed sheep occurred ∼11,800-9,000 years BP and ∼5,300-1,700 years BP, respectively. These findings provide important insights into the history of migration and pastoralism of sheep across the Old World, which was associated with different breeding goals during the Neolithic agricultural revolution.

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.

Organization and Management of Conservation Programs and Research in Domestic Animal Genetic Resources

Creating national committees for domestic animal genetic resources within genetic resource national commissions is recommended to organize in situ and ex situ conservation initiatives. In situ conservation is a high priority because it retains traditional zootechnical contexts and locations to ensure the long-term survival of breeds. In situ actions can be based on subsidies, technical support, structure creation, or trademark definition. Provisional or permanent relocation of breeds may prevent immediate extinction when catastrophes, epizootics, or social conflicts compromise in situ conservation. Ex situ in vivo (animal preservation in rescue or quarantine centers) and in vitro methods (germplasm, tissues/cells, DNA/genes storage) are also potential options. Alert systems must detect emergencies and summon the national committee to implement appropriate procedures. Ex situ coordinated centers must be prepared to permanently or provisionally receive extremely endangered collections. National germplasm banks must maintain sufficient samples of national breeds (duplicated) in their collections to restore extinct populations at levels that guarantee the survival of biodiversity. A conservation management survey, describing national and international governmental and non-governmental structures, was developed. Conservation research initiatives for international domestic animal genetic resources from consortia centralize the efforts of studies on molecular, genomic or geo-evolutionary breed characterization, breed distinction, and functional gene identification. Several consortia also consider ex situ conservation relying on socioeconomic or cultural aspects. The CONBIAND network (Conservation for the Biodiversity of Local Domestic Animals for Sustainable Rural Development) exemplifies conservation efficiency maximization in a low-funding setting, integrating several Latin American consortia with international cooperation where limited human, material, and economic resources are available.

Morphostructural Characterization of the Black Creole Goat Raised in Central Mexico, a Currently Threatened Zoogenetic Resource

In order to evaluate the morphostructural variability of the Black Creole goat (BCG), the present study was carried out in a population of 226 animals from eight localities and 14 morphometric variables were taken. Descriptive statistics for the variables were obtained and 10 of these presented variation coefficients of less than 10%. The degree of harmony in the morphology of the population was determined by the number of positive correlations with significant differences (p < 0.05), including a correlation test using Spearman's method. In order to reduce the matrix of variables, a principal components analysis was performed, and it was evaluated based on Kaiser's criteria (eigenvalue > 1). Finally, a hierarchical analysis of conglomerates using Ward's method was performed using the Euclidean distance to evaluate the distances among localities. Morphometric variables were also included to visualize the relationship among the localities and their average per variable. The results showed that the animals evaluated presented a certain degree of homogeneity and maintained a highly harmonic model. The BCG population showed a high aptitude for milk production, which confirmed the zootechnical purpose of the breed. The BCG populations evaluated maintain similar morphostructural profiles specific to them that can distinguish this population from other animal breeds.

Genome-wide diversity and demographic dynamics of Cameroon goats and their divergence from east African, north African, and Asian conspecifics

Indigenous goats make significant contributions to Cameroon’s national and local economy, but little effort has been devoted to identifying the populations. Here, we assessed the genetic diversity and demographic dynamics of Cameroon goat populations using mitochondrial DNA (two populations) and autosomal markers (four populations) generated with the Caprine 50K SNP chip. To infer genetic relationships at continental and global level, genotype data on six goat populations from Ethiopia and one population each from Egypt, Morocco, Iran, and China were included in the analysis. The mtDNA analysis revealed 83 haplotypes, all belonging to haplogroup A, in Cameroon goats. Four haplotypes were shared between goats found in Cameroon, Mozambique, Namibia, Zimbabwe, Kenya, and Ethiopia. Analysis of autosomal SNPs in Cameroon goats revealed the lowest H_O (0.335±0.13) and H_E (0.352±0.15) in the North-west Highland and Central Highland populations, respectively. Overall, the highest H_O (0.401±0.12) and H_E (0.422±0.12) were found for Barki and Iranian goats, respectively. Barki goats had the highest average MAF, while Central Highland Cameroon goats had the lowest. Overall, Cameroon goats demonstrated high F_IS. AMOVA revealed that 13.29% of the variation was explained by genetic differences between the six population groups. Low average F_ST (0.01) suggests intermixing among Cameroon goats. All measures indicated that Cameroon goats are closer to Moroccan goats than to other goat populations. PCA and STRUCTURE analyses poorly differentiated the Cameroon goats, as did genetic distance, Neighbor-Net network, and neighbor-joining tree analyses. The haplotype analysis of mtDNA showed the initial dispersion of goats to Cameroon and central Africa from north-east Africa following the Nile Delta. Whereas, the approximate Bayesian computation indicated Cameroon goats were separated from Moroccan goats after 506 generations in later times (~1518 YA), as supported by the phylogenetic net-work and admixture outputs. Overall, indigenous goats in Cameroon show weak phylogenetic structure, suggesting either extensive intermixing.

New world goat populations are a genetically diverse reservoir for future use

Western hemisphere goats have European, African and Central Asian origins, and some local or rare breeds are reported to be adapted to their environments and economically important. By-in-large these genetic resources have not been quantified. Using 50 K SNP genotypes of 244 animals from 12 goat populations in United States, Costa Rica, Brazil and Argentina, we evaluated the genetic diversity, population structure and selective sweeps documenting goat migration to the "New World". Our findings suggest the concept of breed, particularly among "locally adapted" breeds, is not a meaningful way to characterize goat populations. The USA Spanish goats were found to be an important genetic reservoir, sharing genomic composition with the wild ancestor and with specialized breeds (e.g. Angora, Lamancha and Saanen). Results suggest goats in the Americas have substantial genetic diversity to use in selection and promote environmental adaptation or product driven specialization. These findings highlight the importance of maintaining goat conservation programs and suggest an awaiting reservoir of genetic diversity for breeding and research while simultaneously discarding concerns about breed designations.

The Eurasian Steppe is an important goat propagation route: A phylogeographic analysis using mitochondrial DNA and Y-chromosome sequences of Kazakhstani goats

Goats (Capra hircus) were domesticated in the Fertile Crescent and propagated all over the world. The Silk Road through the Eurasian Steppe belt is a possible propagation route for domestic goats to Central Asia. Kazakhstan is in close geographical proximity to domestication centers and covers the majority of the Eurasian Steppe belt. In this study, we examined the genetic diversity and phylogeographic structure of Kazakhstani goats. The mtDNA sequences of 141 Kazakhstani goats were categorized into haplogroups A, C, and D, of which haplogroup A was predominant (97%), whereas haplogroups C and D were detected at low frequencies (1.4% each). The Kazakhstani haplotypes C were thzen categorized into Asian mtDNA type. Sequence analysis of the SRY gene on the Y-chromosome in 67 male Kazakhstani goats revealed two haplotypes: Y1A (64%) and Y2A (36%). Analysis of the distribution of mtDNA haplogroups and SRY haplotypes from Eurasia and Africa demonstrated genetic similarity among animals from Kazakhstan, Mongolia, and Northwest China located on the Eurasian Steppe belt. These phylogeographic results suggested that the Eurasian Steppe belt was an important propagation route for goats to Central Asia.

Genome-wide patterns of homozygosity provide clues about the population history and adaptation of goats

Background

Patterns of homozygosity can be influenced by several factors, such as demography, recombination, and selection. Using the goat SNP50 BeadChip, we genotyped 3171 goats belonging to 117 populations with a worldwide distribution. Our objectives were to characterize the number and length of runs of homozygosity (ROH) and to detect ROH hotspots in order to gain new insights into the consequences of neutral and selection processes on the genome-wide homozygosity patterns of goats.

Results

The proportion of the goat genome covered by ROH is, in general, less than 15% with an inverse relationship between ROH length and frequency i.e. short ROH (< 3 Mb) are the most frequent ones. Our data also indicate that ~ 60% of the breeds display low F_ROH coefficients (< 0.10), while ~ 30 and ~ 10% of the goat populations show moderate (0.10 < F_ROH < 0.20) or high (> 0.20) F_ROH values. For populations from Asia, the average number of ROH is smaller and their coverage is lower in goats from the Near East than in goats from Central Asia, which is consistent with the role of the Fertile Crescent as the primary centre of goat domestication. We also observed that local breeds with small population sizes tend to have a larger fraction of the genome covered by ROH compared to breeds with tens or hundreds of thousands of individuals. Five regions on three goat chromosomes i.e. 11, 12 and 18, contain ROH hotspots that overlap with signatures of selection.

Conclusions

Patterns of homozygosity (average number of ROH of 77 and genome coverage of 248 Mb; F_ROH < 0.15) are similar in goats from different geographic areas. The increased homozygosity in local breeds is the consequence of their small population size and geographic isolation as well as of founder effects and recent inbreeding. The existence of three ROH hotspots that co-localize with signatures of selection demonstrates that selection has also played an important role in increasing the homozygosity of specific regions in the goat genome. Finally, most of the goat breeds analysed in this work display low levels of homozygosity, which is favourable for their genetic management and viability.

Electronic supplementary material

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

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 homogeneity of North-African goats

North Africa represents a rich and early reservoir of goat genetic diversity, from which the main African breeds have been derived. In this study, the genetic diversity of four indigenous Algerian goat breeds (i.e., Arabia, Makatia, M’Zabite and Kabyle, with n = 12 for each breed) has been investigated for the first time by genome-wide SNP genotyping; moreover in a broader context, genetic structuration of Algerian and Moroccan goats was explored (via F_ST, MDS, STRUCTURE, FineSTRUCTURE, BAPS, sPCA and DAPC analyses). At national level, the study revealed high level of genetic diversity and a significant phenomenon of admixture affecting all the Algerian breeds. At broader scale, clear global genetic homogeneity appeared considering both Algerian and Moroccan stocks. Indeed, genetic structuration was almost nonexistent among Arabia (from Algeria), Draa, Black and Nord (from Morocco), while the ancestral Kabyle and M’Zabite breeds, reared by Berber peoples, showed genetic distinctness. The study highlighted the threat to the Maghrebin stock, probably induced by unsupervised cross-breeding practices which have intensified in recent centuries. Moreover, it underlined the necessity to deepen our understanding of the genetic resources represented by the resilient North-African goat stock.

Genetic homogeneity of goat malaria parasites in Asia and Africa suggests their expansion with domestic goat host

Plasmodium was first identified in a goat in Angola in 1923, and only recently characterized by DNA isolation from a goat blood sample in Zambia. Goats were first domesticated in the Fertile Crescent approximately 10,000 years ago, and are now globally distributed. It is not known if the Plasmodium identified in African goats originated from parasites circulating in the local ungulates, or if it co-evolved in the goat before its domestication. To address this question, we performed PCR-based surveillance using a total of 1,299 goat blood samples collected from Sudan and Kenya in Africa, Iran in west Asia, and Myanmar and Thailand in southeast Asia. Plasmodium DNA was detected from all locations, suggesting that the parasite is not limited to Africa, but widely distributed. Whole mitochondrial DNA sequences revealed that there was only one nucleotide substitution between Zambian/Kenyan samples and others, supporting the existence of a goat-specific Plasmodium species, presumably Plasmodium caprae, rather than infection of goats by local ungulate malaria parasites. We also present the first photographic images of P. caprae, from one Kenyan goat sample.

Convergent genomic signatures of domestication in sheep and goats

The evolutionary basis of domestication has been a longstanding question and its genetic architecture is becoming more tractable as more domestic species become genome-enabled. Before becoming established worldwide, sheep and goats were domesticated in the fertile crescent 10,500 years before present (YBP) where their wild relatives remain. Here we sequence the genomes of wild Asiatic mouflon and Bezoar ibex in the sheep and goat domestication center and compare their genomes with that of domestics from local, traditional, and improved breeds. Among the genomic regions carrying selective sweeps differentiating domestic breeds from wild populations, which are associated among others to genes involved in nervous system, immunity and productivity traits, 20 are common to Capra and Ovis. The patterns of selection vary between species, suggesting that while common targets of selection related to domestication and improvement exist, different solutions have arisen to achieve similar phenotypic end-points within these closely related livestock species.

The sheep and goat were domesticated ~10,500 years ago in the same region of the Middle-East. Here, Alberto et al compare the genomes of wild Asiatic mouflon and Bezoar ibex with that of domestics from local, traditional and improved breeds and find common targets of selection related to domestication and improvement in sheep and goats.

Genome-wide population structure and admixture analysis reveals weak differentiation among Ugandan goat breeds

Uganda has a large population of goats, predominantly from indigenous breeds reared in diverse production systems, whose existence is threatened by crossbreeding with exotic Boer goats. Knowledge about the genetic characteristics and relationships among these Ugandan goat breeds and the potential admixture with Boer goats is still limited. Using a medium‐density single nucleotide polymorphism (SNP) panel, we assessed the genetic diversity, population structure and admixture in six goat breeds in Uganda: Boer, Karamojong, Kigezi, Mubende, Small East African and Sebei. All the animals had genotypes for about 46 105 SNPs after quality control. We found high proportions of polymorphic SNPs ranging from 0.885 (Kigezi) to 0.928 (Sebei). The overall mean observed (H_O) and expected (H_E) heterozygosity across breeds was 0.355 ± 0.147 and 0.384 ± 0.143 respectively. Principal components, genetic distances and admixture analyses revealed weak population sub‐structuring among the breeds. Principal components separated Kigezi and weakly Small East African from other indigenous goats. Sebei and Karamojong were tightly entangled together, whereas Mubende occupied a more central position with high admixture from all other local breeds. The Boer breed showed a unique cluster from the Ugandan indigenous goat breeds. The results reflect common ancestry but also some level of geographical differentiation. admixture and f₄ statistics revealed gene flow from Boer and varying levels of genetic admixture among the breeds. Generally, moderate to high levels of genetic variability were observed. Our findings provide useful insights into maintaining genetic diversity and designing appropriate breeding programs to exploit within‐breed diversity and heterozygote advantage in crossbreeding schemes.

Dissection of ancestral genetic contributions to Creole goat populations

Goats have played a key role as source of nourishment for humans in their expansion all over the world in long land and sea trips. This has guaranteed a place for this species in the important and rapid episode of livestock expansion triggered by Columbus' arrival in the Americas in the late 1400s. The aims of this study are to provide a comprehensive perspective on genetic diversity in American goat populations and to assess their origins and evolutionary trajectories. This was achieved by combining data from autosomal neutral genetic markers obtained in more than two thousand samples that encompass a wide range of Iberian, African and Creole goat breeds. In general, even though Creole populations differ clearly from each other, they lack a strong geographical pattern of differentiation, such that populations of different admixed ancestry share relatively close locations throughout the large geographical range included in this study. Important Iberian signatures were detected in most Creole populations studied, and many of them, particularly the Cuban Creole, also revealed an important contribution of African breeds. On the other hand, the Brazilian breeds showed a particular genetic structure and were clearly separated from the other Creole populations, with some influence from Cape Verde goats. These results provide a comprehensive characterisation of the present structure of goat genetic diversity, and a dissection of the Iberian and African influences that gave origin to different Creole caprine breeds, disentangling an important part of their evolutionary history. Creole breeds constitute an important reservoir of genetic diversity that justifies the development of appropriate management systems aimed at improving performance without loss of genomic diversity.

Differential distribution of Y-chromosome haplotypes in Swiss and Southern European goat breeds

The analysis of Y-chromosome variation has provided valuable clues about the paternal history of domestic animal populations. The main goal of the current work was to characterize Y-chromosome diversity in 31 goat populations from Central Eastern (Switzerland and Romania) and Southern Europe (Spain and Italy) as well as in reference populations from Africa and the Near East. Towards this end, we have genotyped seven single nucleotide polymorphisms (SNPs), mapping to the SRY, ZFY, AMELY and DDX3Y Y-linked loci, in 275 bucks from 31 populations. We have observed a low level of variability in the goat Y-chromosome, with just five haplotypes segregating in the whole set of populations. We have also found that Swiss bucks carry exclusively Y1 haplotypes (Y1A: 24%, Y1B1: 15%, Y1B2: 43% and Y1C: 18%), while in Italian and Spanish bucks Y2A is the most abundant haplotype (77%). Interestingly, in Carpathian goats from Romania the Y2A haplotype is also frequent (42%). The high Y-chromosome differentiation between Swiss and Italian/Spanish breeds might be due to the post-domestication spread of two different Near Eastern genetic stocks through the Danubian and Mediterranean corridors. Historical gene flow between Southern European and Northern African goats might have also contributed to generate such pattern of genetic differentiation.

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.

Detecting the existence of gene flow between Spanish and North African goats through a coalescent approach

Human-driven migrations are one of the main processes shaping the genetic diversity and population structure of domestic species. However, their magnitude and direction have been rarely analysed in a statistical framework. We aimed to estimate the impact of migration on the population structure of Spanish and African goats. To achieve this goal, we analysed a dataset of 1,472 individuals typed with 23 microsatellites. Population structure of African and Spanish goats was moderate (mean FST = 0.07), with the exception of the Canarian and South African breeds that displayed a significant differentiation when compared to goats from North Africa and Nigeria. Measurement of gene flow with Migrate-n and IMa coalescent genealogy samplers supported the existence of a bidirectional gene flow between African and Spanish goats. Moreover, IMa estimates of the effective number of migrants were remarkably lower than those calculated with Migrate-n and classical approaches. Such discrepancies suggest that recent divergence, rather than extensive gene flow, is the main cause of the weak population structure observed in caprine breeds.

A genome-wide perspective about the diversity and demographic history of seven Spanish goat breeds

Background

The main goal of the current work was to infer the demographic history of seven Spanish goat breeds (Malagueña, Murciano-Granadina, Florida, Palmera, Mallorquina, Bermeya and Blanca de Rasquera) based on genome-wide diversity data generated with the Illumina Goat SNP50 BeadChip (population size, N = 176). Five additional populations from Europe (Saanen and Carpathian) and Africa (Tunisian, Djallonké and Sahel) were also included in this analysis (N = 80) for comparative purposes.

Results

Our results show that the genetic background of Spanish goats traces back mainly to European breeds although signs of North African admixture were detected in two Andalusian breeds (Malagueña and Murciano-Granadina). In general, observed and expected heterozygosities were quite similar across the seven Spanish goat breeds under analysis irrespective of their population size and conservation status. For the Mallorquina and Blanca de Rasquera breeds, which have suffered strong population declines during the past decades, we observed increased frequencies of large-sized (ROH), a finding that is consistent with recent inbreeding. In contrast, a substantial part of the genome of the Palmera goat breed comprised short ROH, which suggests a strong and ancient founder effect.

Conclusions

Admixture with African goats, genetic drift and inbreeding have had different effects across the seven Spanish goat breeds analysed in the current work. This has generated distinct patterns of genome-wide diversity that provide new clues about the demographic history of these populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0229-6) contains supplementary material, which is available to authorized users.

Genetic Diversity and Population Structure in South African, French and Argentinian Angora Goats from Genome-Wide SNP Data

The Angora goat populations in Argentina (AR), France (FR) and South Africa (SA) have been kept geographically and genetically distinct. Due to country-specific selection and breeding strategies, there is a need to characterize the populations on a genetic level. In this study we analysed genetic variability of Angora goats from three distinct geographical regions using the standardized 50k Goat SNP Chip. A total of 104 goats (AR: 30; FR: 26; SA: 48) were genotyped. Heterozygosity values as well as inbreeding coefficients across all autosomes per population were calculated. Diversity, as measured by expected heterozygosity (H_E) ranged from 0.371 in the SA population to 0.397 in the AR population. The SA goats were the only population with a positive average inbreeding coefficient value of 0.009. After merging the three datasets, standard QC and LD-pruning, 15 105 SNPs remained for further analyses. Principal component and clustering analyses were used to visualize individual relationships within and between populations. All SA Angora goats were separated from the others and formed a well-defined, unique cluster, while outliers were identified in the FR and AR breeds. Apparent admixture between the AR and FR populations was observed, while both these populations showed signs of having some common ancestry with the SA goats. LD averaged over adjacent loci within the three populations per chromosome were calculated. The highest LD values estimated across populations were observed in the shorter intervals across populations. The N_e for the Angora breed was estimated to be 149 animals ten generations ago indicating a declining trend. Results confirmed that geographic isolation and different selection strategies caused genetic distinctiveness between the populations.

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

Background

The current extensive use of the domestic goat (Capra hircus) is the result of its medium size and high adaptability as multiple breeds. The extent to which its genetic variability was influenced by early domestication practices is largely unknown. A common standard by which to analyze maternally-inherited variability of livestock species is through complete sequencing of the entire mitogenome (mitochondrial DNA, mtDNA).

Results

We present the first extensive survey of goat mitogenomic variability based on 84 complete sequences selected from an initial collection of 758 samples that represent 60 different breeds of C. hircus, as well as its wild sister species, bezoar (Capra aegagrus) from Iran. Our phylogenetic analyses dated the most recent common ancestor of C. hircus to ~460,000 years (ka) ago and identified five distinctive domestic haplogroups (A, B1, C1a, D1 and G). More than 90 % of goats examined were in haplogroup A. These domestic lineages are predominantly nested within C. aegagrus branches, diverged concomitantly at the interface between the Epipaleolithic and early Neolithic periods, and underwent a dramatic expansion starting from ~12–10 ka ago.

Conclusions

Domestic goat mitogenomes descended from a small number of founding haplotypes that underwent domestication after surviving the last glacial maximum in the Near Eastern refuges. All modern haplotypes A probably descended from a single (or at most a few closely related) female C. aegagrus. Zooarchaelogical data indicate that domestication first occurred in Southeastern Anatolia. Goats accompanying the first Neolithic migration waves into the Mediterranean were already characterized by two ancestral A and C variants. The ancient separation of the C branch (~130 ka ago) suggests a genetically distinct population that could have been involved in a second event of domestication. The novel diagnostic mutational motifs defined here, which distinguish wild and domestic haplogroups, could be used to understand phylogenetic relationships among modern breeds and ancient remains and to evaluate whether selection differentially affected mitochondrial genome variants during the development of economically important breeds.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2342-2) contains supplementary material, which is available to authorized users.

Prospects and challenges for the conservation of farm animal genomic resources, 2015-2025

Livestock conservation practice is changing rapidly in light of policy developments, climate change and diversifying market demands. The last decade has seen a step change in technology and analytical approaches available to define, manage and conserve Farm Animal Genomic Resources (FAnGR). However, these rapid changes pose challenges for FAnGR conservation in terms of technological continuity, analytical capacity and integrative methodologies needed to fully exploit new, multidimensional data. The final conference of the ESF Genomic Resources program aimed to address these interdisciplinary problems in an attempt to contribute to the agenda for research and policy development directions during the coming decade. By 2020, according to the Convention on Biodiversity's Aichi Target 13, signatories should ensure that “…the genetic diversity of …farmed and domesticated animals and of wild relatives …is maintained, and strategies have been developed and implemented for minimizing genetic erosion and safeguarding their genetic diversity.” However, the real extent of genetic erosion is very difficult to measure using current data. Therefore, this challenging target demands better coverage, understanding and utilization of genomic and environmental data, the development of optimized ways to integrate these data with social and other sciences and policy analysis to enable more flexible, evidence-based models to underpin FAnGR conservation. At the conference, we attempted to identify the most important problems for effective livestock genomic resource conservation during the next decade. Twenty priority questions were identified that could be broadly categorized into challenges related to methodology, analytical approaches, data management and conservation. It should be acknowledged here that while the focus of our meeting was predominantly around genetics, genomics and animal science, many of the practical challenges facing conservation of genomic resources are societal in origin and are predicated on the value (e.g., socio-economic and cultural) of these resources to farmers, rural communities and society as a whole. The overall conclusion is that despite the fact that the livestock sector has been relatively well-organized in the application of genetic methodologies to date, there is still a large gap between the current state-of-the-art in the use of tools to characterize genomic resources and its application to many non-commercial and local breeds, hampering the consistent utilization of genetic and genomic data as indicators of genetic erosion and diversity. The livestock genomic sector therefore needs to make a concerted effort in the coming decade to enable to the democratization of the powerful tools that are now at its disposal, and to ensure that they are applied in the context of breed conservation as well as development.

Paternal phylogeography and genetic diversity of East Asian goats

This study was a first analysis of paternal genetic diversity for extensive Asian domestic goats using SRY gene sequences. Sequencing comparison of the SRY 3'-untranslated region among 210 Asian goats revealed four haplotypes (Y1A, Y1B, Y2A and Y2B) derived from four variable sites including a novel substitution detected in this study. In Asian goats, the predominant haplotype was Y1A (62%) and second most common was Y2B (30%). Interestingly, the Y2B was a unique East Asian Y chromosomal variant, which differentiates eastern and western Eurasian goats. The SRY geographic distribution in Myanmar and Cambodia indicated predominant the haplotype Y1A in plains areas and a high frequency of Y2B in mountain areas. The results suggest recent genetic infiltration of modern breeds into South-East Asian goats and an ancestral SRY Y2B haplotype in Asian native goats.

Characterizing neutral genomic diversity and selection signatures in indigenous populations of Moroccan goats (Capra hircus) using WGS data

Since the time of their domestication, goats (Capra hircus) have evolved in a large variety of locally adapted populations in response to different human and environmental pressures. In the present era, many indigenous populations are threatened with extinction due to their substitution by cosmopolitan breeds, while they might represent highly valuable genomic resources. It is thus crucial to characterize the neutral and adaptive genetic diversity of indigenous populations. A fine characterization of whole genome variation in farm animals is now possible by using new sequencing technologies. We sequenced the complete genome at 12× coverage of 44 goats geographically representative of the three phenotypically distinct indigenous populations in Morocco. The study of mitochondrial genomes showed a high diversity exclusively restricted to the haplogroup A. The 44 nuclear genomes showed a very high diversity (24 million variants) associated with low linkage disequilibrium. The overall genetic diversity was weakly structured according to geography and phenotypes. When looking for signals of positive selection in each population we identified many candidate genes, several of which gave insights into the metabolic pathways or biological processes involved in the adaptation to local conditions (e.g., panting in warm/desert conditions). This study highlights the interest of WGS data to characterize livestock genomic diversity. It illustrates the valuable genetic richness present in indigenous populations that have to be sustainably managed and may represent valuable genetic resources for the long-term preservation of the species.

Y-chromosomal variation of local goat breeds of Turkey close to the domestication centre

Genetic variations in chromosome Y are enabling researchers to identify paternal lineages, which are informative for introgressions and migrations. In this study, the male-specific region markers, sex-determining region-Y (SRY), amelogenin (AMELY) and zinc finger (ZFY) were analysed in seven Turkish native goat breeds, Angora, Kilis, Hair, Honamlı, Norduz, Gürcü and Abaza. A SNP in the ZFY gene defined a new haplotype Y2C. All domestic haplogroups originate from Capra aegagrus, while the finding of Y1A, Y1B, Y2A and Y2C in 32, 4, 126 and 2 Turkish domestic goats, respectively, appears to indicate a predomestic origin of the major haplotypes. The occurrence of four haplotypes in the Hair goat and, in contrast, a frequency of 96% of Y1A in the Kilis breed illustrate that Y-chromosomal variants have a more breed-dependent distribution than mitochondrial or autosomal DNA. This probably reflects male founder effects, but a role in adaptation cannot be excluded.