Genetic diversity and patterns of population structure in Creole goats from the Americas

Assessing the population structure and genetic variability of Kenyan native goats under extensive production system

Indigenous goats are important to many livelihoods. Despite this, they are subjected to indiscriminate crossbreeding. This affects their genetic variability which is needed to survive in current regime of climate change. The study assessed population structure and genetic diversity of Galla and Small East African goats (SEA) using pedigree information. A total of 7384 animals, 5222 of the Galla and 2162 of the SEA breeds, born between the years 1983 and 2022, were utilized. Individuals with known parents were defined as reference population. From the results, the maximum generation traced for Galla and SEA populations was 14.6 and 14.5, respectively. However, only 6 and 5 generations for Galla and SEA were complete. Pedigree completeness increased with the increasing number of generations. The average generation interval (GI) for Galla and SEA was 3.84 ± 0.04 and 4.4 ± 0.13 years. The average increase in the rate of inbreeding per generation for Galla and SEA was 0.04 and 0.05, with the effective ancestors and founders (fa/fe) ratio being same (1.00) for both breeds. Fifty percent (50%) of the genetic variability in the populations was contributed by 3 and 1 ancestor for Galla SEA, respectively. The effective population size (Ne) was 5.19 and 4.77 for Galla and SEA. Therefore, the current breeding programs should be changed to avoid future genetic bottlenecks in this population. These findings offer an opportunity to enhance the current genetic status and management of Kenyan native goats and other regions with similar production systems.

Milk Composition of Creole Goats Raised at Different Altitudes in an Extensive Production System in Northeast Mexico

Goat milk composition is affected by feeding, and in semiarid rangeland, information on Creole goat milk physicochemical composition is lacking. For the fulfillment of this objective, three agroecological regions (AR) considering altitude (lowland 87, highland 779, and mountain 1309 m above sea level) with different botanical compositions were chosen. Every AR analyzed accounted for 30 goat herds, with a total of 90 herds. The results demonstrated that altitude had an influence mainly on density and protein. Milk density increases as altitude increases; conversely, milk protein increases as altitude decreases. On the other hand, in the mountain and lowland ARs, the salts and solids not fat (SNF) percentages were higher compared to that of the highland AR (p < 0.05). The freezing point (FP) was higher at highland altitudes compared to that of mountain and lowland ARs (p < 0.01). In the milk fatty acids (FA) profile, only the C14:1 value was affected by altitude, whereas goat milk at lowland and mountain altitudes had higher values compared to that at highland altitudes (p < 0.05). Additionally, late lactation stage fat, FP, and pH values were higher compared to early lactation values. The opposite effect was observed for salts and SNF. In the FA profile, late lactation values were higher for C10:0 and C8:0 compared to early lactation values. The opposite trend was observed for C18:2n6t. The thrombogenic index was significantly higher at lowland altitudes compared to highland altitudes, and similar to the mountain AR. These goat milk characteristics could be explained as a consequence of animal nutrition, as well as the goat's meat-type phenotype.

Trajectory of livestock genomics in South Asia: A comprehensive review

Livestock plays a central role in sustaining human livelihood in South Asia. There are numerous and distinct livestock species in South Asian countries. Several of them have experienced genetic development in recent years due to the application of genomic technologies and effective breeding programs. This review discusses genomic studies on cattle, buffalo, sheep, goat, pig, horse, camel, yak, mithun, and poultry. The frontiers covered in this review are genetic diversity, admixture studies, selection signature research, QTL discovery, genome-wide association studies (GWAS), and genomic selection. The review concludes with recommendations for South Asian livestock systems to increasingly leverage genomic technologies, based on the lessons learned from the numerous case studies. This paper aims to present a comprehensive analysis of the dichotomy in the South Asian livestock sector and argues that a realistic approach to genomics in livestock can ensure long-term genetic advancements.

Status quo of genetic improvement in local goats: a review

This review aims to summarize and synthesize the fragmented information available on the genetic improvement of local goats (criollo, indigenous, native) on the American and other continents, where populations with these goats have an important role in food security and the economy of rural communities, as well as in conservation of biodiversity and productivity improvement. Topics such as the current state of goat production globally, conservation programs, resistance to parasites and diseases, use of phenotypical characteristics and genomic information, and molecular markers for genetic improvement are addressed. The main challenges, opportunities, and limitations described in recent literature concerning local goats in the immediate future are discussed.

An Overview of the Use of Genotyping Techniques for Assessing Genetic Diversity in Local Farm Animal Breeds

Simple Summary

The number of local farm animal breeds is declining worldwide. However, these breeds have different degrees of genetic diversity. Measuring genetic diversity is important for the development of conservation strategies and, therefore, various genomic analysis techniques are available. The aim of the present work was to shed light on the use of these techniques in diversity studies of local breeds. In summary, a total of 133 worldwide studies that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that over time, almost all available genomic techniques were used and various diversity parameters were calculated. Therefore, the present results provide a comprehensive overview of the application of these techniques in the field of local breeds. This can provide helpful insights into the advancement of the conservation of breeds with high genetic diversity.

Abstract

Globally, many local farm animal breeds are threatened with extinction. However, these breeds contribute to the high amount of genetic diversity required to combat unforeseen future challenges of livestock production systems. To assess genetic diversity, various genotyping techniques have been developed. Based on the respective genomic information, different parameters, e.g., heterozygosity, allele frequencies and inbreeding coefficient, can be measured in order to reveal genetic diversity between and within breeds. The aim of the present work was to shed light on the use of genotyping techniques in the field of local farm animal breeds. Therefore, a total of 133 studies across the world that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that diversity of cattle was most often investigated with microsatellite use as the main technique. Furthermore, a large variety of diversity parameters that were calculated with different programs were identified. For 15% of the included studies, the used genotypes are publicly available, and, in 6%, phenotypes were recorded. In conclusion, the present results provide a comprehensive overview of the application of genotyping techniques in the field of local breeds. This can provide helpful insights to advance the conservation of breeds.

Assessing Genetic Diversity and Estimating the Inbreeding Effect on Economic Traits of Inner Mongolia White Cashmere Goats Through Pedigree Analysis

Objective: The purpose of this study was to discover the population structure and genetic diversity of Inner Mongolia White Cashmere goats (IMCGs) and demonstrate the effect of inbreeding on the live body weight (LBW), cashmere yield (CY), fiber length (FL), and fiber diameter (FD) of IMCGs.

Materials and Methods: All data were collected from pedigree information and production performance records of IMCGs from 1983 to 2019. The population structure and genetic diversity were analyzed by Endog 4.8 software. Inbreeding coefficients were obtained by the pedigree package in R. Then, a linear regression model was used to analyze how inbreeding influences economic traits in IMCGs. Four levels of inbreeding coefficients (F_i) were classified in this study, including F_i = 0, 0< F_i ≤ 6.25, 6.25< F_i ≤ 12.5 and F_i≥12.5. Variance analysis was performed to determine whether inbreeding levels had a significant effect on economic traits in IMCGs.

Results: The proportions of rams and dams in IMCGs for breeding were relatively small, with values of 0.8 and 20.5%, respectively. The proportion of inbred animals in the entire population was high, with values up to 68.6%; however, the average inbreeding coefficient and relatedness coefficient were 4.50 and 8.48%, respectively. To date, the population has experienced 12 generations. The average generation interval obtained in the present study was 4.11 ± 0.01 years. The ram-to-son pathway was lowest (3.97 years), and the ewe-to-daughter pathway was highest (4.24 years). It was discovered that the LBW, CY, and FL increased by 3.88 kg, 208.7 g, and 1.151 cm, respectively, with every 1% increase in the inbreeding coefficient, and the FD decreased by 0.819 μm with every 1% increase in the inbreeding coefficient. Additionally, multiple comparison analysis indicated that when the inbreeding coefficient was higher than 6.25%, the LBW showed an obvious decreasing trend. The threshold value of inbreeding depression in the CY is 12.5%. However, inbreeding depression has not been observed in the FL and FD.

Conclusion: Pedigree completeness needs to be further strengthened. The degree of inbreeding in this flock should be properly controlled when designing breeding programs.

The Portuguese Serrana goat breed: a review

Goats were among the first animals to be domesticated over 10,000 years ago and are part of human societies since the beginning of agriculture. Goats play a major role both in commercial farming systems and in subsistence agriculture systems, particularly in tropical, subtropical and Mediterranean regions where they are crucial for the supply of meat, milk, fibre and dung. This review concerns the Serrana breed, the most important and numerous indigenous goat breed from Portugal that was furthermore exported to other regions of the world, notably South America during the Portuguese colonization. Herein, we describe the origin and history of the breed as well as the productive performance and most common production systems. Finally, we address the local and traditional PDO (protected denomination of origin) and PGI (protected geographical indication) that are produced from these animals.

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.

Genetic Characterization of the "Chusca Lojana", a Creole Goat Reared in Ecuador, and Its Relationship with Other Goat Breeds

The largest population of goats (62%) in Ecuador is in the dry forest region in the south of the country. A Creole goat, named "Chusca Lojana", has adapted to the dry forest region where environmental conditions are warm-dry, with sparse vegetation. Knowledge of the genetic information of the Creole goat is important to determine intra-racial diversity, the degree of genetic distance among other breeds of goats, and the possible substructure of the population, which is valuable for the conservation of such a species' genetic resources. A total of 145 samples of the Creole goat was taken from the four biotypes previously identified. Genetic analyses were performed using 38 microsatellites recommended for studies of goat genetic diversity (FAO-ISAG). The results of within-breed genetic diversity showed a mean number of alleles per locus (MNA) of 8, an effective number of alleles (Ae) of 4.3, an expected heterozygosity (He) of 0.71, an observed heterozygosity (Ho) of 0.63, polymorphic information content (PIC) of 0.67, and an FIS value of 0.11. Between-breed genetic diversity among 43 goat populations (native of Spain, American Creole, Europeans, and Africans) showed the following values: FIS = 0.087, FIT = 0.176, and FST = 0.098. Regarding the analysis of the population structure, the results showed that the Creole Chusca Lojana goat population is homogeneous and no genetic separation was observed between the different biotypes (FST = 0.0073). In conclusion, the Chusca Lojana goat has a high genetic diversity, without exhibiting a genetic substructure. Therefore, it should be considered as a distinct population because crossbreeding with other breeds was not detected.

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.

Conservation Strategies for Local Breed Biodiversity

Breeds of domesticated animals are often overlooked as contributing to biodiversity. Their unique role at the junction of natural and human-influenced environments makes them important potential contributors to biodiversity and ecosystem function. Effective conservation of livestock breeds rests upon valid recognition of breeds as repeatable genomic packages with a reasonably high level of predictability for performance. Local or landrace breeds that lack breeder organization are especially difficult to conserve due to lack of formal recognition as breeds. Achieving success with them involves three major steps: Discover, Secure, and Sustain. Early in the process an evaluation of candidate populations for status as genetic resources is essential. This process is aided by a phenotypic matrix which can be used alongside historical investigations and genetic (DNA) studies. The goal is to include all qualifying animals and to exclude all those that do not qualify. Securing some populations depends on careful rescue protocols for maximizing the recovery of genetic variation, and this can then be followed by breeding protocols that provide for maintaining the population’s production potential along with a healthy and viable genetic structure for long-term survival and use. Sustaining breeds for the long term is also enabled by assuring market demand for the breed and its products.

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.

Characterization of Duck (Anas platyrhynchos) Short Tandem Repeat Variation by Population-Scale Genome Resequencing

Short tandem repeats (STRs) are usually associated with genetic diseases and gene regulatory functions, and are also important genetic markers for analysis of evolutionary, genetic diversity and forensic. However, for the majority of STRs in the duck genome, their population genetic properties and functional impacts remain poorly defined. Recent advent of next generation sequencing (NGS) has offered an opportunity for profiling large numbers of polymorphic STRs. Here, we reported a population-scale analysis of STR variation using genome resequencing in mallard and Pekin duck. Our analysis provided the first genome-wide duck STR reference including 198,022 STR loci with motif size of 2–6 base pairs. We observed a relatively uneven distribution of STRs in different genomic regions, which indicates that the occurrence of STRs in duck genome is not random, but undergoes a directional selection pressure. Using genome resequencing data of 23 mallard and 26 Pekin ducks, we successfully identified 89,891 polymorphic STR loci. Intensive analysis of this dataset suggested that shorter repeat motif, longer reference tract length, higher purity, and residing outside of a coding region are all associated with an increase in STR variability. STR genotypes were utilized for population genetic analysis, and the results showed that population structure and divergence patterns among population groups can be efficiently captured. In addition, comparison between Pekin duck and mallard identified 3,122 STRs with extremely divergent allele frequency, which overlapped with a set of genes related to nervous system, energy metabolism and behavior. The evolutionary analysis revealed that the genes containing divergent STRs may play important roles in phenotypic changes during duck domestication. The variation analysis of STRs in population scale provides valuable resource for future study of genetic diversity and genome evolution in duck.

Genetic characterization of indigenous goat breeds in Romania and Hungary with a special focus on genetic resistance to mastitis and gastrointestinal parasitism based on 40 SNPs

Goat breeding has become an important sector in Eastern Europe, with Romania and Hungary being among the major producer countries. Given the limited number of research done up-to-date concerning genetic studies of indigenous goat breeds reared in Romania and Hungary, the current preliminary study aimed to analyze the variability of genes related to mastitis and gastrointestinal parasitism by using Kompetitive Allele Specific PCR (KASP™). We studied 52 single nucleotide polymorphisms (SNPs) belonging to 19 genes in indigenous breeds from both countries, namely Banat's White (n = 36), Carpatina (n = 35) from Romania and Hungarian Milking (n = 79) and identified 16 polymorphic SNPs among 10 genes (PTX3, IL6, CLEC4E, IL8, IL1RN, IL15RA, TNFSF13, SOCS3, TNF and TLR3) in 150 animals. Furthermore, the diversity of the studied breeds was investigated. The PIC values ranged from 0.042 to 0.691. The mean values of observed and expected heterozygosity were 0.235 and 0.246 respectively. The highest observed heterozygosity was obtained for IL15RA g.10343904C>T in Banat's White (0.464), IL15RA g.10354813C>T in Carpatina (0.577) and SOCS3 g.52626440T>G in Hungarian Milking (0.588). Pairwise FST values between the Romanian breeds and Romanian and Hungarian breeds were small (0.009 and 0.015), indicating the close relationship among the studied goat populations. From all the polymorphic SNPs identified, the Hungarian Milking breed showed the highest proportion of polymorphisms (100%), whereas the Carpatina breed had the lowest percentage (87.5%). The highest value of MAF was obtained for SOCS3 g.52626440T>G (0.46), IL15RA g.10343904C>T (0.47), IL15RA g.10344025C>T (0.45), and IL15RA g.10354813C>T (0.42). The 16 polymorphic SNPs identified in a panel of 150 unrelated individuals belonging to three Romanian and Hungarian indigenous goat breeds could be used in future genomic based breeding schemes as markers for genetic resistance to mastitis and gastrointestinal parasitism in goat breeds found in Eastern and Central Europe.

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.

Genetic diversity and structure related to expansion history and habitat isolation: stone marten populating rural-urban habitats

Background

Population genetic diversity and structure are determined by past and current evolutionary processes, among which spatially limited dispersal, genetic drift, and shifts in species distribution boundaries have major effects. In most wildlife species, environmental modifications by humans often lead to contraction of species’ ranges and/or limit their dispersal by acting as environmental barriers. However, in species well adapted to anthropogenic habitat or open landscapes, human induced environmental changes may facilitate dispersal and range expansions. In this study, we analysed whether isolation by distance and deforestation, among other environmental features, promotes or restricts dispersal and expansion in stone marten (Martes foina) populations.

Results

We genotyped 298 martens from eight sites at twenty-two microsatellite loci to characterize the genetic variability, population structure and demographic history of stone martens in Poland. At the landscape scale, limited genetic differentiation between sites in a mosaic of urban, rural and forest habitats was mostly influenced by isolation by distance. Statistical clustering and multivariate analyses showed weak genetic structuring with two to four clusters and a high rate of gene flow between them. Stronger genetic differentiation was detected for one stone marten population (NE1) located inside a large forest complex. Genetic differentiation between this site and all others was 20% higher than between other sites separated by similar distances. The genetic uniqueness index of NE1 was also twofold higher than in other sites. Past demographic history analyses showed recent expansion of this species in north-eastern Poland. A decrease in genetic diversity from south to north, and MIGRAINE analyses indicated the direction of expansion of stone marten.

Conclusions

Our results showed that two processes, changes in species distribution boundaries and limited dispersal associated with landscape barriers, affect genetic diversity and structure in stone marten. Analysis of local barriers that reduced dispersal and large scale analyses of genetic structure and demographic history highlight the importance of isolation by distance and forest cover for the past colonization of central Europe by stone marten. This confirmed the hypothesis that human-landscape changes (deforestation) accelerated stone marten expansion, to which climate warming probably has also been contributing over the last few decades.

Electronic supplementary material

The online version of this article (10.1186/s12898-017-0156-6) contains supplementary material, which is available to authorized users.

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.