SNP-based association mapping of the polled gene in divergent cattle breeds

Naturally, hornless cattle are called polled. Although the POLL locus could be assigned to a c. 1.36‐Mb interval in the centromeric region of BTA1, the underlying genetic basis for the polled trait is still unknown. Here, an association mapping design was set up to refine the candidate region of the polled trait for subsequent high‐throughput sequencing. The case group comprised 101 homozygous polled animals from nine divergent cattle breeds, the majority represented by Galloway, Angus, Fleckvieh and Holstein Friesian. Additionally, this group included some polled individuals of Blonde d’Aquitaine, Charolais, Hereford, Jersey and Limousin breeds. The control group comprised horned Belgian Blue, Fleckvieh, Holstein Friesian and Illyrian Bu?a cattle. A genome‐wide scan using 49 163 SNPs was performed, which revealed one shared homozygous haplotype block consisting of nine neighbouring SNPs in all polled animals. This segment defines a 381‐kb interval on BTA1 that we consider to be the most likely location of the POLL mutation. Our results further demonstrate that the polled‐associated haplotype is also frequent in horned animals included in this study, and thus the haplotype as such cannot be used for population‐wide genetic testing. The actual trait‐associated haplotype may be revealed by using higher‐density SNP arrays. For the final identification of the causal mutation, we suggest high‐throughput sequencing of the entire candidate region, because the identification of functional candidate genes is difficult owing to the lack of a comparable model.     

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.     

Mapping the four‐horned locus and testing the polled locus in three Chinese sheep breeds

Four‐horned sheep are an ideal animal model for illuminating the genetic basis of horn development. The objective of this study was to locate the genetic region responsible for the four‐horned phenotype and to verify a previously reported polled locus in three Chinese breeds. A genome‐wide association study (GWAS) was performed using 34 two‐horned and 32 four‐horned sheep from three Chinese indigenous breeds: Altay, Mongolian and Sishui Fur sheep. The top two significant single nucleotide polymorphisms (SNPs) associated with the four‐horned phenotype were both located in a region spanning positions 132.6 to 132.7 Mb on sheep chromosome 2. Similar locations for the four‐horned trait were previously identified in Jacob, Navajo‐Churro, Damara and Sishui Fur sheep, suggesting a common genetic component underlying the four‐horned phenotype. The two identified SNPs were both downstream of the metaxin 2 (MTX2) gene and the HOXD gene cluster. For the top SNP—OAR2:g.132619300G>A—the strong associations of the AA and AG genotypes with the four‐horned phenotype and the GG genotype with the two‐horned phenotype indicated the dominant inheritance of the four‐horned trait. No significant SNPs for the polled phenotype were identified in the GWAS analysis, and a PCR analysis for the detection of the 1.8‐kb insertion associated with polled sheep in other breeds failed to verify the association with polledness in the three Chinese breeds. This study supports the hypothesis that two different loci are responsible for horn existence and number. This study contributes to the understanding of the molecular regulation of horn development and enriches the knowledge of qualitative traits in domestic animals.     

The brown coat colour of Coppernecked goats is associated with a non‐synonymous variant at the TYRP1 locus on chromosome 8

The recent development of a goat SNP genotyping microarray enables genome‐wide association studies in this important livestock species. We investigated the genetic basis of the black and brown coat colour in Valais Blacknecked and Coppernecked goats. A genome‐wide association analysis using goat SNP50 BeadChip genotypes of 22 cases and 23 controls allowed us to map the locus for the brown coat colour to goat chromosome 8. The TYRP1 gene is located within the associated chromosomal region, and TYRP1 variants cause similar coat colour phenotypes in different species. We thus considered TYRP1 as a strong positional and functional candidate. We resequenced the caprine TYRP1 gene by Sanger and Illumina sequencing and identified two non‐synonymous variants, p.Ile478Thr and p.Gly496Asp, that might have a functional impact on the TYRP1 protein. However, based on the obtained pedigree and genotype data, the brown coat colour in these goats is not due to a single recessive loss‐of‐function allele. Surprisingly, the genotype distribution and the pedigree data suggest that the ⁴⁹⁶Asp allele might possibly act in a dominant manner. The ⁴⁹⁶Asp allele was present in 77 of 81 investigated Coppernecked goats and did not occur in black goats. This strongly suggests heterogeneity underlying the brown coat colour in Coppernecked goats. Functional experiments or targeted matings will be required to verify the unexpected preliminary findings.     

Characterization of the <Emphasis Type="Italic">GHR</Emphasis> gene genetic variation in Chinese indigenous goat breeds

The aim of the present work was to investigate single nucleotide polymorphism (SNP) of growth hormone receptor (GHR) gene exon 10, characterize the genetic variation in three Chinese indigenous goat breeds, and search for its potential association with cashmere traits. In this study, a polymerase chain reaction-single strand conformation polymorphism (PCR–SSCP) protocol has been developed for rapid genotyping of the GHR gene in goats. One hundred seventy-eight goats from Liaoning Cashmere (96), Inner Mongolia White Cashmere (40), and Chengdu Grey (42) breeds in China were genotyped at GHR locus using the protocol developed. In all goat breeds investigated, a SNP in exon 10 of GHR gene has been identified by analyzing genomic DNA. The polymorphism consists of a single nucleotide substitution A → G, resulting in two alleles named, respectively, A and G based on the nucleotide at the position. The allele A was found to be more common in the animals investigated, and seems to be more consistent with cattle and zebu at this polymorphic site found in goats. The Hardy–Weinberg equilibrium of genotype distributions of GHR locus was verified in Liaoning Cashmere, and Inner Mongolia White Cashmere breeds. According to the classification of polymorphism information content (PIC), Chengdu Grey was less polymorphic than Liaoning Cashmere and Inner Mongolia White Cashmere breeds at this locus. The phylogenetic tree of different species based on the nucleotide sequences of GHR gene exon 10 is generally in agreement with the known species relationship. No significant association was found between the polymorphism revealed and the cashmere traits analyzed in present work.     

Polymorphisms of caprine <Emphasis Type="Italic">GDF</Emphasis>9 gene and their association with litter size in Jining Grey goats

The exons 1, 2 and flanking region of growth differentiation factor 9 (GDF9) gene in five randomly selected does of Jining Grey, Boer and Liaoning Cashmere goats were amplified and analyzed. Thirteen nucleotide differences were identified in GDF9 gene between sheep (AF078545) and goats. Four SNPs (G3288A in intron 1, G423A, A959C [Gln320Pro] and G1189A [Val397Ile] in exon 2) were detected in four goat breeds with different prolificacy, in which G3288A was a new SNP in goats. The results showed that loci 3288, 423 and 1189 in Boer goats, loci 3288 and 423 in Guizhou White goats, loci 423 and 1189 in Liaoning Cashmere goats were all in complete linkage disequilibrium (D′ = 1, r ² = 1), respectively. In moderate (Boer goat) and low prolificacy (Liaoning Cashmere goat) breeds, linkage analysis indicated that there were more fervent linkage disequilibrium among loci 3288, 423 and 1189 than high prolificacy (Jining Grey and Guizhou White goats) breeds. For the 959 locus, the genotype distribution showed obvious difference between high prolificacy breeds and moderate or low prolificacy breeds (P < 0.05 or P < 0.01). The Jining Grey goat does with genotype CC or AC had 0.81 (P < 0.01) or 0.63 (P < 0.01) kids more than those with genotype AA, respectively. The present study preliminarily showed an association between allele C at 959 locus of GDF9 gene and high litter size in Jining Grey goats. These results provide further evidence that the GDF9 gene may be significantly correlated with high prolificacy in goats.     

Wattles in goats are associated with the FMN1/GREM1 region on chromosome 10

The presence of congenital appendages (wattles) on the throat of goats is supposed to be under genetic control with a dominant mode of inheritance. Wattles contain a cartilaginous core covered with normal skin resembling early stages of extremities. To map the dominant caprine wattles (W) locus, we collected samples of 174 goats with wattles and 167 goats without wattles from nine different Swiss goat breeds. The samples were genotyped with the 53k goat SNP chip for a subsequent genome‐wide association study. We obtained a single strong association signal on chromosome 10 in a region containing functional candidate genes for limb development and outgrowth. We sequenced the whole genomes of an informative family trio containing an offspring without wattles and its heterozygous parents with wattles. In the associated goat chromosome 10 region, a total of 1055 SNPs and short indels perfectly co‐segregate with the W allele. None of the variants were perfectly associated with the phenotype after analyzing the genome sequences of eight additional goats. We speculate that the causative mutation is located in one of the numerous gaps in the current version of the goat reference sequence and/or represents a larger structural variant which influences the expression of the FMN1 and/or GREM1 genes. Also, we cannot rule out possible genetic or allelic heterogeneity. Our genetic findings support earlier assumptions that wattles are rudimentary developed extremities.     

A single nucleotide polymorphism on chromosome 10 is highly predictive for the polled phenotype in Australian Merino sheep

The aim of this study was to fine map the genomic location of the Horns locus in the Australian Merino sheep population and to identify markers that can be used to predict the horn phenotype. A linkage disequilibrium analysis of horn data from Australian Merino sheep mapped the Horns locus to a small region on chromosome 10. A single nucleotide polymorphism in the region was found to be highly predictive for the polled phenotype in an experimental population of Merino sheep. This was owing to a dominance effect of one of the alleles when inherited maternally. It was suggested that a genetic test would provide a good predictor of the polled phenotype. Finally, an evaluation of industry data showed that the SNP is at very different frequencies in Poll Merino sheep that have been bred for polledness (based on phenotype alone) compared with the Merino sheep breed.     

A Single Nucleotide Polymorphism within the Interferon Gamma Receptor 2 Gene Perfectly Coincides with Polledness in Holstein Cattle

Polledness is a high impact trait in modern milk and beef production to meet the demands of animal welfare and work safety. Previous studies have mapped the polled-locus to the proximal region of the bovine chromosome 1 (BTA1) and narrowed it down to approximately 1 Mb. Sequencing of the positional candidate genes within the 1 Mb polled region and whole genome sequencing of Holsteins revealed a single nucleotide polymorphism (SNP) AC000158: g.1390292G>A within intron 3 of the interferon gamma receptor 2 gene (IFNGR2) in perfect co-segregation with polledness in Holsteins. This complete association was validated in 443 animals of the same breed. This SNP allows reliable genotyping of horned, heterozygous and homozygous polled Holsteins, even in animals that could not be resolved using the previously published haplotype for Holstein.     

群体遗传不平衡条件下的结构基因遗传共适应特性

本研究以柴达木山羊、柴达木绒山羊和辽宁绒山羊三个群体共147只山羊为材料,利用聚丙烯酰胺凝胶电泳（PAGE）技术检测了5种血液蛋白质（酶）基因座的遗传多态性,并进行了结构基因遗传共适应的研究,结果发现：45个基因座组合中有10个基因座组合处于遗传不平衡状态,并且这些遗传不平衡皆单纯由遗传共适应差异造成;除辽宁绒山羊Tf-PA-3组合的遗传不平衡包含非等位基因间的遗传共适应差异外,其他基因座组合的遗传不平衡皆由等位基因间的共适应差异,即单基因座的遗传不平衡造成;LAP-EsD组合的共适应差异在群体间有遗传传递现象。 Abstract:With the technology of PAGE,the genetic polymorphism of blood protein and enzyme was investigated,and genetic co-adaptability among structural genes was studied in three goat populations(147 goats) including Chaidamu goat(CS),Chaidamu Cashmere goat(CRS) and Liaoning Cashmere goat(LRS) in Qinghai Province,China.The results were showed that the genetic disequilibrium of 10 locus combinations was found among 45 locus combinations in the three goat populations,and these genetic disequilibria were caused only by the difference of genetic co-adaptability among genes,because there didn′t exist the linkage disequilibrium among non-allelic genes.The genetic disequilibrium including the difference of genetic co-adaptability between non-allelic genes was only found at Tf-PA-3 locus combinations in LRS population,the other ones were all caused by the genetic disequilibrium at a single locus.The difference of genetic co-adaptability of LAP-EsD locus combinations could be messaged among different populations.     

Design and Characterization of a 52K SNP Chip for Goats

The success of Genome Wide Association Studies in the discovery of sequence variation linked to complex traits in humans has increased interest in high throughput SNP genotyping assays in livestock species. Primary goals are QTL detection and genomic selection. The purpose here was design of a 50–60,000 SNP chip for goats. The success of a moderate density SNP assay depends on reliable bioinformatic SNP detection procedures, the technological success rate of the SNP design, even spacing of SNPs on the genome and selection of Minor Allele Frequencies (MAF) suitable to use in diverse breeds. Through the federation of three SNP discovery projects consolidated as the International Goat Genome Consortium, we have identified approximately twelve million high quality SNP variants in the goat genome stored in a database together with their biological and technical characteristics. These SNPs were identified within and between six breeds (meat, milk and mixed): Alpine, Boer, Creole, Katjang, Saanen and Savanna, comprising a total of 97 animals. Whole genome and Reduced Representation Library sequences were aligned on >10 kb scaffolds of the de novo goat genome assembly. The 60,000 selected SNPs, evenly spaced on the goat genome, were submitted for oligo manufacturing (Illumina, Inc) and published in dbSNP along with flanking sequences and map position on goat assemblies (i.e. scaffolds and pseudo-chromosomes), sheep genome V2 and cattle UMD3.1 assembly. Ten breeds were then used to validate the SNP content and 52,295 loci could be successfully genotyped and used to generate a final cluster file. The combined strategy of using mainly whole genome Next Generation Sequencing and mapping on a contig genome assembly, complemented with Illumina design tools proved to be efficient in producing this GoatSNP50 chip. Advances in use of molecular markers are expected to accelerate goat genomic studies in coming years.     

Contribution of domestic animals to the identification of new genes involved in sex determination.

Among farm animals, two species present an intersex condition at a relatively high frequency: pig and goat. Both are known to contain XX sex-reversed individuals which are genetically female but with a true hermaphrodite or male phenotype. It has been clearly demonstrated that the SRY gene is not involved in these phenotypes. Consequently, autosomal or X-linked mutations in the sex-determining pathway may explain these sex-reversed phenotypes. A mutation referred to as "polled" has been characterized in goats by the suppression of horn formation and abnormal sexual differentiation. The Polled Intersex Syndrome locus (PIS) was initially located in the distal region of goat chromosome 1. The homologous human region has been precisely identified as an HSA 3q23 DNA segment containing the Blepharophimosis Ptosis Epicanthus locus (BPES), a syndrome combining Premature Ovarian Failure (POF) and an excess of epidermis of the eyelids. In order to isolate genes involved in pig intersexuality, a similar genetic approach was attempted in pigs using genome scanning of resource families. Genetic analyses suggest that pig intersexuality is controlled multigenically. Parallel to this work, gonads of fetal intersex animals have been studied during development by light and electron microscopy. The development of testicular tissue and reduction of germ cell number by apoptosis, which simultaneously occurs as soon as 50 days post co?tum, also suggests that several separate genes could be involved in pig intersexuality.     

FSHβ 基因PCR-SSCP多态性及其与济宁青山羊高繁殖力关系的研究

采用PCR-SSCP技术检测促卵泡素b亚基(follicle-stimulating hormone β, FSHβ)基因5′调控区、外显子1和外显子2在高繁殖力山羊品种(济宁青山羊)和低繁殖力山羊品种(辽宁绒山羊、波尔山羊、安哥拉山羊)中的单核苷酸多态性, 同时研究该基因对济宁青山羊高繁殖力的影响。结果表明: 山羊与绵羊的FSHβ 基因该段核苷酸序列同源性为98%。9对引物中, 只有P9的扩增片段存在多态性。P9的扩增片段在济宁青山羊和辽宁绒山羊中检测到AA、AB和AC 3种基因型; 在波尔山羊中检测到AA、CC和AC 3种基因型; 在安哥拉山羊中检测到AA、BB、CC、AB、AC和BC共6种基因型。测序分析发现BB型与AA型相比在外显子2的第94 bp处有G→A突变, 并引起氨基酸改变(丙氨酸→苏氨酸); CC型与AA型相比在外显子2的第174 bp有一处C→T沉默突变。济宁青山羊AA、AB和AC基因型频率分别为0.686、0.137和0.177。AA基因型济宁青山羊产羔数最小二乘均值比AB基因型的多0.78只(P<0.05), 比AC基因型的多0.64只(P<0.05)。     

催乳素受体基因外显子10多态性及其与济宁青山羊高繁殖力关系的研究

设计5对引物, 采用PCR-SSCP技术检测催乳素受体(prolactin receptor, PRLR)基因外显子10及部分3′非翻译区在高繁殖力山羊品种(济宁青山羊)和低繁殖力山羊品种(辽宁绒山羊、波尔山羊和安哥拉山羊)中的单核苷酸多态性, 同时研究该基因对济宁青山羊高繁殖力的影响。结果表明: 首次拼接出的山羊PRLR基因外显子10及部分3′非翻译区的核苷酸序列长度为1,118 bp, 与已公布的绵羊、牛、人PRLR基因mRNA相应序列的同源性分别为98.33%、93.92%、74.52%, 与已公布的羊驼PRLR基因部分序列的同源性为78.29%。引物P1、P2与P4扩增片段具有多态性, 其余2对引物的扩增片段不存在多态性。对于P1扩增片段, 在济宁青山羊和辽宁绒山羊中检测到AA型和AB型, 在安哥拉山羊中检测到AA型和AC型, 在波尔山羊中只检测到AA型; 克隆测序表明AB型与AA型相比有两处突变(186G→A和220T→C), 分别导致氨基酸由天冬氨酸变为天冬酰胺、亮氨酸变为脯氨酸; AC型与AA型相比有1处突变(140A→G), 该突变没有导致氨基酸变化; 济宁青山羊AA和AB基因型之间产羔数的最小二乘均值差异不显著(P>0.05)。对于P2扩增片段, 在济宁青山羊、辽宁绒山羊和波尔山羊中都检测到DD型和DE型, 而在安哥拉山羊中只检测到DD型; 克隆测序表明DE型和DD型相比有两处突变(52G→A和122G→A), 其中122 bp处的突变导致氨基酸由精氨酸变为甘氨酸; 济宁青山羊DD和DE基因型之间产羔数的最小二乘均值差异不显著(P>0.05)。对于P4扩增片段, 在济宁青山羊中检测到FF型和FG型, 在辽宁绒山羊中检测到FF型和GG型, 在波尔山羊中只检测到FF型, 在安哥拉山羊中检测到FF型、FG型和GG型; 克隆测序表明GG型和FF型相比在扩增片段的143 bp处发生1处碱基突变(A→G), 并导致氨基酸由蛋氨酸变为缬氨酸; FG基因型济宁青山羊产羔数最小二乘均值比FF基因型的多0.76只(P<0.05)。研究结果初步表明: PRLR基因可能是控制济宁青山羊多胎性能的一个主效基因或是与之存在紧密遗传连锁的一个标记。     

Polymorphism of bone morphogenetic protein 4 gene and its relationship with litter size of Jining Grey goats

Two pairs of primers (P1 and P2) were designed to detect single nucleotide polymorphisms of exon 2 and intron 2 of bone morphogenetic protein 4 (BMP4) gene in both high fecundity breed (Jining Grey goat) and low fecundity breeds (Boer, Angora and Inner Mongolia Cashmere goats) by single strand conformation polymorphism. Results showed that no polymorphism was detected for exon 2 (primer P1) of BMP4 gene in four goat breeds. For intron 2 (primer P2), three genotypes (AA, AB and BB) were detected in Jining Grey and Inner Mongolia Cashmere goats, two genotypes (AB and BB) in Angora goats, and only one genotype (AA) in Boer goats. Sequencing revealed one mutation (2203G>A) of BMP4 gene in the genotype BB in comparison to the genotype AA. The differences of litter size between AA, AB and BB genotypes were not significant (P > 0.05) in Jining Grey goats. A pair of primer (P3) was designed to detect polymorphism in the 3' flanking region of BMP4 gene that contained dinucleotide repeated sequence (CA) in the four goat breeds by microsatellite analysis. For primer P3, three genotypes (CC, CD and DD) were detected in four goat breeds. Sequencing revealed one more CA dinucleotide in genotype DD than in genotype CC. The Jining Grey does with genotype CC had 0.55 (P < 0.05) or 0.72 (P < 0.05) kids more than those with genotype CD or DD. These results preliminarily indicated that allele C of BMP4 gene is a potential DNA marker for improving litter size in goats.     

Polledness in Argentinean Creole cattle,five centuries surviving

Polledness has been shown to have autosomal Mendelian inheritance, with the polled locus being dominant to the horned locus. This trait was mapped to the BTA1 centromeric end in several breeds. One of the distinctive attributes of Creole cattle, such as the Argentinean Creole, is the presence of long, lyre‐shaped horns. However, polled native animals were reported before the introduction of modern selected European breeds. Here, we studied the origin of the polled mutation, either independent or introgressed, in a Creole line from the Creole cattle founder group at the IIACS‐INTA Leales Experimental Station (northwest Argentina). The study sample (65 animals: 26 horned and 39 polled) was genotyped using high‐density SNP microarrays and three previously reported genetic markers (P_202ID, P_80kbID and P_G). A genome‐wide association study, selection signatures, linkage disequilibrium analysis and copy number variations were used to detect the responsible region and the segregating haplotypes/alleles. The interval mapped in the Leales herd (1.23–2.13 Mb) overlapped with the region previously reported in several European cattle breeds, suggesting that the same locus could be segregating in this population. The previously reported variants P_F and P_G were not detected, thus dismissing the Holstein‐Friesian and Nellore origins of the polled phenotype in this native breed. Conversely, the presence of the Celtic variant P_C suggests an almost complete co‐segregation. The cluster analysis rejected the hypothesis of recent introgression, which is compatible with the historical record of polled Creole cattle in northwest Argentina.     

Genome‐wide association study of body morphological traits in Sudanese goats

Long‐term selection of goats for a certain production system and/or different environmental conditions will be reflected in the body morphology of the animals under selection. To investigate the variation contributing to different morphological traits and to identify genomic regions that are associated with body morphological traits in Sudanese goats, we genotyped 96 females belonging to four Sudanese goat breeds with the SNP52 BeadChip. After quality control of the data, the genome‐wide association study was performed using 95 goats and 24 027 informative single nucleotide polymorphisms (SNPs). Bicoastal diameter was significantly associated (LOD = 6.32) with snp10185‐scaffold1365‐620922 on chromosome 2. The minor allele has an additive effect, increasing the bicoastal diameter by 2.6 cm. A second significant association was found between body length and snp56482‐scaffold89‐467312 on chromosome 3 (LOD = 5.65). The minor allele is associated with increased body length. Additionally, five regions were suggestive for cannon bone, head width, rump length and withers height (LOD > 5). Only one gene (CNTNAP5) is located within the 1‐Mb region surrounding the significant SNP for bicoastal diameter on chromosome 2. The body length QTL on chromosome 3 harbors 49 genes. Further research is required to validate the observed associations and to prioritize candidate genes.     

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.