Targeted high‐throughput growth hormone 1 gene sequencing reveals high within‐breed genetic diversity in South African goats

This study assessed the genetic diversity in the growth hormone 1 gene (GH1) within and between South African goat breeds. Polymerase chain reaction‐targeted gene amplification together with Illumina MiSeq next‐generation sequencing (NGS) was used to generate the full length (2.54 kb) of the growth hormone 1 gene and screen for SNPs in the South African Boer (SAB) (n = 17), Tankwa (n = 15) and South African village (n = 35) goat populations. A range of 27–58 SNPs per population were observed. Mutations resulting in amino acid changes were observed at exons 2 and 5. Higher within‐breed diversity of 97.37% was observed within the population category consisting of SA village ecotypes and the Tankwa goats. Highest pairwise F_ST values ranging from 0.148 to 0.356 were observed between the SAB and both the South African village and Tankwa feral goat populations. Phylogenetic analysis indicated nine genetic clusters, which reflected close relationships between the South African populations and the other international breeds with the exception of the Italian Sarda breeds. Results imply greater potential for within‐population selection programs, particularly with SA village goats.     

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

Biodiversity studies are more efficient when large numbers of breeds belonging to several countries are involved, as they allow for an in‐depth analysis of the within‐ and between‐breed components of genetic diversity. A set of 21 microsatellites was used to investigate the genetic composition of 24 Creole goat breeds (910 animals) from 10 countries to estimate levels of genetic variability, infer population structure and understand genetic relationships among populations across the American continent. Three commercial transboundary breeds were included in the analyses to investigate admixture with Creole goats. Overall, the genetic diversity of Creole populations (mean number of alleles = 5.82 ± 1.14, observed heterozygosity = 0.585 ± 0.074) was moderate and slightly lower than what was detected in other studies with breeds from other regions. The Bayesian clustering analysis without prior information on source populations identified 22 breed clusters. Three groups comprised more than one population, namely from Brazil (Azul and Graúna; Moxotó and Repartida) and Argentina (Long and shorthair Chilluda, Pampeana Colorada and Angora‐type goat). Substructure was found in Criolla Paraguaya. When prior information on sample origin was considered, 92% of the individuals were assigned to the source population (threshold q ≥ 0.700). Creole breeds are well‐differentiated entities (mean coefficient of genetic differentiation = 0.111 ± 0.048, with the exception of isolated island populations). Dilution from admixture with commercial transboundary breeds appears to be negligible. Significant levels of inbreeding were detected (inbreeding coefficient > 0 in most Creole goat populations, P < 0.05). Our results provide a broad perspective on the extant genetic diversity of Creole goats, however further studies are needed to understand whether the observed geographical patterns of population structure may reflect the mode of goat colonization in the Americas.     

Genetic diversity and population structure in multiple Chinese goat populations using a SNP panel

Information about genetic diversity and population structure among goat breeds is essential for genetic improvement, understanding of environmental adaptation as well as utilization and conservation of goat breeds. Here, we measured genetic diversity and population structure in multiple Chinese goat populations, namely, Nanjiang, Qinggeda, Arbas Cashmere, Jining Grey, Luoping Yellow and Guangfeng goats. A total of 193 individuals were genotyped for about 47 401 autosomal single nucleotide polymorphisms (SNPs). We found a high proportion of informative SNPs, ranging from 69.5% in the Luoping Yellow to 93.9% in the Jining Grey goat breeds with an average mean of 84.7%. Diversity, as measured by expected heterozygosity, ranged from 0.371 in Luoping Yellow to 0.405 in Jining Grey goat populations. The average estimated pair‐wise genetic differentiation (F_ST) among the populations was 8.6%, ranging from 0.2% to 16% and indicating low to moderate genetic differentiation. Principal component analysis, genetic structure and phylogenetic tree analysis revealed a clustering of six Chinese goat populations according to geographic distribution. The results from this study can contribute valuable genetic information and can properly assist with within‐breed diversity, which provides a good opportunity for sustainable utilization of and maintenance of genetic resource improvements in the Chinese goat populations.     

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.     

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

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

Two missense mutations in exon 9 of caprine PRLR gene were associated with litter size

Guanzhong (n = 321) and Boer (n = 191) goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the coding regions of the prolactin receptor (PRLR) gene by DNA sequencing and PCR‐RFLP. Two SNPs (c.1457G>A and c.1645G>A) were identified that caused amino acid variations p.Ser485Asn and p.Val548Met respectively. Statistical results indicated that the c.1457G>A and c.1645G>A SNPs were significantly associated with litter size in Boer and Guanzhong goat breeds. Further analysis revealed that combined genotype C4 (GGGG) and haplotype G‐G were better than the others for litter size in both goat breeds. These results might contribute to goat genetic resources and breeding.     

Genome‐wide SNP analysis unveils genetic structure and phylogeographic history of snow sheep (Ovis nivicola) populations inhabiting the Verkhoyansk Mountains and Momsky Ridge (northeastern Siberia)

Insights into the genetic characteristics of a species provide important information for wildlife conservation programs. Here, we used the OvineSNP50 BeadChip developed for domestic sheep to examine population structure and evaluate genetic diversity of snow sheep (Ovis nivicola) inhabiting Verkhoyansk Range and Momsky Ridge. A total of 1,121 polymorphic SNPs were used to test 80 specimens representing five populations, including four populations of the Verkhoyansk Mountain chain: Kharaulakh Ridge–Tiksi Bay (TIK, n = 22), Orulgan Ridge (ORU, n = 22), the central part of Verkhoyansk Range (VER, n = 15), Suntar‐Khayata Ridge (SKH, n = 13), and Momsky Ridge (MOM, n = 8). We showed that the studied populations were genetically structured according to a geographic pattern. Pairwise F_ST values ranged from 0.044 to 0.205. Admixture analysis identified K = 2 as the most likely number of ancestral populations. A Neighbor‐Net tree showed that TIK was an isolated group related to the main network through ORU. TreeMix analysis revealed that TIK and MOM originated from two different ancestral populations and detected gene flow from MOM to ORU. This was supported by the f3 statistic, which showed that ORU is an admixed population with TIK and MOM/SKH heritage. Genetic diversity in the studied groups was increasing southward. Minimum values of observed (Ho) and expected (He) heterozygosity and allelic richness (Ar) were observed in the most northern population—TIK, and maximum values were observed in the most southern population—SKH. Thus, our results revealed clear genetic structure in the studied populations of snow sheep and showed that TIK has a different origin from MOM, SKH, and VER even though they are conventionally considered a single subspecies known as Yakut snow sheep (Ovis nivicola lydekkeri). Most likely, TIK was an isolated group during the Late Pleistocene glaciations of Verkhoyansk Range.     

Searching new signals for production traits through gene‐based association analysis in three Italian cattle breeds

Genome‐wide association studies (GWAS) have been widely applied to disentangle the genetic basis of complex traits. In cattle breeds, classical GWAS approaches with medium‐density marker panels are far from conclusive, especially for complex traits. This is due to the intrinsic limitations of GWAS and the assumptions that are made to step from the association signals to the functional variations. Here, we applied a gene‐based strategy to prioritize genotype–phenotype associations found for milk production and quality traits with classical approaches in three Italian dairy cattle breeds with different sample sizes (Italian Brown n = 745; Italian Holstein n = 2058; Italian Simmental n = 477). Although classical regression on single markers revealed only a single genome‐wide significant genotype–phenotype association, for Italian Holstein, the gene‐based approach identified specific genes in each breed that are associated with milk physiology and mammary gland development. As no standard method has yet been established to step from variation to functional units (i.e., genes), the strategy proposed here may contribute to revealing new genes that play significant roles in complex traits, such as those investigated here, amplifying low association signals using a gene‐centric approach.     

Taking advantage from phenotype variability in a local animal genetic resource: identification of genomic regions associated with the hairless phenotype in Casertana pigs

Casertana is an endangered autochthonous pig breed (raised in south‐central Italy) that is considered to be the descendant of the influential Neapolitan pig population that was used to improve British breeds in the 19th century. Casertana pigs are characterized by a typical, almost complete, hairless phenotype, even though a few Casertana pigs are normal haired. In this work, using Illumina PorcineSNP60 BeadChip data, we carried out a genome‐wide association study and an F_ST analysis with this breed by comparing animals showing the classical hairless phenotype (n = 81) versus pigs classified as haired (n = 15). Combining the results obtained with the two approaches, we identified two significant regions: one on porcine chromosome (SSC) 7 and one on SSC15. The SSC7 region contains the forkhead box N3 (FOXN3) gene, the most plausible candidate gene of this region, considering that mutations in another gene of the same family (forkhead box N1; Foxn1 or FOXN1) are responsible for the nude locus in rodents and alopecia in humans. Another potential candidate gene, rho guanine nucleotide exchange factor 10 (ARHGEF10), is located in the SSC15 region. FOXN3 and ARHGEF10 have been detected as differentially expressed in androgenetic and senescent alopecia respectively. This study on an autochthonous pig breed contributes to shed some light on novel genes potentially involved in hair development and growth and demonstrates that local animal breeds can be valuable genetic resources for disclosing genetic factors affecting unique traits, taking advantage of phenotype variability segregating in small populations.     

Novel Y‐chromosome short tandem repeats in Sus scrofa and their variation in European wild boar and domestic pig populations

Y‐chromosome markers are important tools for studying male‐specific gene flow within and between populations, hybridization patterns and kinship. However, their use in non‐human mammals is often hampered by the lack of Y‐specific polymorphic markers. We identified new male‐specific short tandem repeats (STRs) in Sus scrofa using the available genome sequence. We selected four polymorphic loci (5–10 alleles per locus), falling in one duplicated and two single‐copy regions. A total of 32 haplotypes were found by screening 211 individuals from eight wild boar populations across Europe and five domestic pig populations. European wild boar were characterized by significantly higher levels of haplotype diversity compared to European domestic pigs (H_D = 0.904 ± 0.011 and H_D = 0.491 ± 0.077 respectively). Relationships among STR haplotypes were investigated by combining them with single nucleotide polymorphisms at two linked genes (AMELY and UTY) in a network analysis. A differentiation between wild and domestic populations was observed (F_ST = 0.229), with commercial breeds sharing no Y haplotype with the sampled wild boar. Similarly, a certain degree of geographic differentiation was observed across Europe, with a number of local private haplotypes and high diversity in northern populations. The described Y‐chromosome markers can be useful to track male inheritance and gene flow in wild and domestic populations, promising to provide insights into evolutionary and population genetics in Sus scrofa.     

Genetic characterization and breed assignment in five Italian sheep breeds using microsatellite markers

The knowledge of the genetic relationship and admixture among neighbouring populations is crucial for conservation efforts. The aim of this study was to analyse the genetic diversity of five Italian sheep breeds (Appenninica, Garfagnina Bianca, Massese, Pomarancina and Zerasca) using a panel of 24 microsatellite markers. Blood samples from 226 individuals belonging to the aforementioned populations were obtained and genotyped. All the investigated breeds showed a significant heterozygote deficiency caused by the high level of inbreeding indicated also by the high level of F_IS (0.146). Genetic differentiation between breeds was moderate (F_ST = 0.05) but significant and the individuals could be assigned to their breeds with an high success rate even if the inter-individual distances showed that few animals clustered separately from the other individuals of the same breed, especially for Pomarancina breed. The genetic distances reflect the historical knowledge of these breeds and some patterns of ancestral and recent gene flow between neighbour populations arise. The clustering analysis detects the presence of six clusters. Massese and Zerasca breeds were grouped together as well as Appenninica and Pomarancina with the latter forming two distinct clusters equally represented. The formation of this last breed is occurred with the absorption of individuals of the Appenninica breed and the gene flow probably continued in these recent years allowing the presence of a population substructure for Pomarancina breed. Such substructure supports the high level of heterozygote deficiency found for this breed despite the relatively high population size. The five populations analysed presented some genetic similarities but a clear uniqueness of the populations has been showed for almost all of them. Special attention to monitor genetic variability and to organize mating plans should be given especially for the three endangered breeds.     

Genetic diversity in Iranian indigenous sheep vis-à-vis selected exogenous sheep breeds and wild mouflon

The heterogeneity of climate and different agro-ecological conditions in Iran have resulted in development of 27 indigenous sheep breeds. Wild Asiatic mouflon (Ovis orientalis) is believed to be the ancestor of Iranian sheep. Evaluation of genetic diversity and population structure within and among domestic breeds has important implications for animal breeding programs and genetic resources management. Based on 50K SNP genotype data, we studied the genetic diversity of five indigenous Iranian sheep breeds: Afshari (n = 37), Moghani (n = 34), Qezel (n = 35), Zel (n = 46) and Lori-Bakhtiari (n = 46), and Asiatic mouflon (n = 8) sampled from Iran. Furthermore, genetic diversity and the breed admixture of Iranian sheep were assessed on a larger geographic scale using a reference panel comprising: three indigenous Afghan breeds – Arabi (n = 15), Balouchi (n = 15) and Gadik (n = 15); three indigenous breeds from Turkey and Cyprus – Cyprus Fat Tail (n = 30), Karakas (n = 18) and Norduz (n = 20); and three commercial European breeds – Suffolk (n = 19), Comisana (n = 24) and Engadine Red Sheep (n = 24). The results revealed that the investigated breeds are divided into five genetically distinct clusters according to their geographic origin. Afshari was closest to the local mouflon population and showed signs of mouflon admixture. Qezel was identified as a hybrid sheep breed. Much evidence supported the Afghan breeds being identical. Inbreeding values, which were estimated based on ROHs, were highest for Suffolk (F_ROH = 0.0544) and lowest for Balouchi (F_ROH = 0.0078). In conclusion, analysis of selected breeds from neighboring countries along with Asiatic mouflon gave a deeper insight into the evolutionary history and origin of Iranian sheep with important implications for future breed management.     

Multilocus genotypic data reveal high genetic diversity and low population genetic structure of Iranian indigenous sheep

Iranian livestock diversity is still largely unexplored, in spite of the interest in the populations historically reared in this country located near the Fertile Crescent, a major livestock domestication centre. In this investigation, the genetic diversity and differentiation of 10 Iranian indigenous fat‐tailed sheep breeds were investigated using 18 microsatellite markers. Iranian breeds were found to host a high level of diversity. This conclusion is substantiated by the large number of alleles observed across loci (average 13.83, range 7–22) and by the high within‐breed expected heterozygosity (average 0.75, range 0.72–0.76). Iranian sheep have a low level of genetic differentiation, as indicated by the analysis of molecular variance, which allocated a very small proportion (1.67%) of total variation to the between‐population component, and by the small fixation index (F_ST = 0.02). Both Bayesian clustering and principal coordinates analysis revealed the absence of a detectable genetic structure. Also, no isolation by distance was observed through comparison of genetic and geographical distances. In spite of high within‐breed variation, signatures of inbreeding were detected by the F_IS indices, which were positive in all and statistically significant in three breeds. Possible factors explaining the patterns observed, such as considerable gene flow and inbreeding probably due to anthropogenic activities in the light of population management and conservation programmes, are discussed.     

Genetic variation of two species with different life‐history traits in the endangered renosterveld of South Africa – a comparative analysis of Eriocephalus africanus and Hemimeris racemosa

We tested the effects of life‐history traits on genetic variation and conducted a comparative analysis of two plant species with differing life‐history traits co‐occurring in the highly endangered renosterveld of South Africa. We selected eighteen renosterveld remnants with varying degrees of size and isolation where populations of the herbaceous, annual and insect‐pollinated Hemimeris racemosa and the shrubby perennial and both wind‐ and insect‐pollinated Eriocephalus africanus occurred. We postulated a lower genetic variation within populations and increased genetic variation between populations in the annual than in the perennial species. Genetic variation was lower within populations of H. racemosa than within E. africanus, as is typical for annual compared to perennial species. Variation within populations was, however, not correlated with fragment size or distance in either of the two species and genetic variation between populations of the two species was comparable (Φ_ST = 0.10, 0.09).     

Analysis of genetic variability at codon 42 within caprine <Emphasis Type="Italic">prion protein</Emphasis> gene in relation to production traits in Chinese domestic breeds

In this study, genetic variability at codon 42 within prion protein (PRNP) gene and its associations with production traits were investigated in 2002 goats from four Chinese domestic breeds. The frequencies of allele “A” ranged from 0.353 to 0.562 in analyzed goat breeds with Hardy–Weinberg equilibrium (P > 0.05) except Xinong Sannen (XNSN) dairy breed. The establishment of relationships between different genotypes and growth traits was performed in Inner Mongolia white Cashmere (IMWC) breed and revealed an association of the polymorphism with body weight at 7-year-old goats (P = 0.033). The individuals with genotype GG showed heavier body weight than those with genotype AA. Moreover, association analysis detected two significant associations between different genotypes and cashmere yield and fiber length in IMWC breed (P = 0.009, P = 0.048, respectively). In addition, three significant associations of different genotypes with density of milk (a.m. and p.m.), solids-not-fat of milk (P = 0.013, P = 0.009 and P = 0.002), respectively, were found in XNSN breed. Genotype GG had better milk quality than others. These findings suggested that the polymorphism of codon 42 within PRNP was a useful DNA marker for eliminating or selecting excellent individuals in relation to production traits in marker-assist selection breeding of goat.     

Genomic diversity and population structure of Mexican and Spanish bovine Lidia breed

The Lidia bovine breed is distinguished for its low genetic exchangeability given its selection on aggressive behavior, its management uniqueness and its subdivided structure. In this study, we present a comprehensive genome‐wide analysis of genetic diversity, population structure and admixture of 468 animals from Mexican and Spanish Lidia breed populations and 64 samples belonging to 10 Spanish native and American‐creole breeds using 37 148 single nucleotide polymorphisms. We found similar average inbreeding values in the Lidia breed, with different distributions within groups; variability of inbreeding values among Spanish lineages was significant and no differences were found among the Mexican sub‐populations. Together, the high F_IS of the lineages and the behavior of the runs of homozygosity are consequences of the lineage's small effective population sizes, contributing to their inbreeding increase. Population admixture analysis discarded any influence on the genetic structure of the Lidia populations from the Spanish native and American‐creole breeds. In addition, both Lidia populations depicted different genetic origins, with the exception of some Mexican individuals whose origins traced back to recent Spanish importations.     

Evaluation of genetic introgression from domesticated pigs into the Ryukyu wild boar population on Iriomote Island in Japan

We evaluated genetic introgression from domesticated pigs into the Ryukyu wild boar (RWB) population on Iriomote Island based on their genetic structure and diversity. We used a combination of mitochondrial DNA D‐loop region (596 bp) polymorphisms and 23 microsatellite markers. RWBs (n = 130) were collected from 18 locations on Iriomote Island and compared with 66 reference samples of European and Asian domestic pigs. We identified six distinct haplotypes, involving 22 single nucleotide polymorphisms (including one insertion) in the RWB population. The phylogenetic tree had two branches: the RWB group and domestic lineage. Fourteen of 130 RWBs (10.8%) belonged to the European domestic lineage, including 11 RWBs from the Panari Islands, northwest of Iriomote Main Island (IMI). The heterozygosity values, total number of alleles, number of effective alleles and polymorphism information content of the RWB groups were lower than those of the European domestic groups. The RWB population on IMI had a lower heterozygous deficiency index (F_IS = 0.059) than did the other populations, which indicates that this population was more inbred. There was a large genetic distance (F_ST = 0.560) between RWBs on IMI and the Meishan populations. Structure analysis using the 23 microsatellite markers revealed that 16 RWBs had an admixture pattern between RWB and domesticated pig breeds. These results suggest that gene flow may have occurred from domestic pigs to RWBs and demonstrate that there was low genetic variation in the IMI population.     

Farm‐by‐farm analysis of microsatellite,mtDNA and SNP genotype data reveals inbreeding and crossbreeding as threats to the survival of a native Spanish pig breed

The Chato Murciano (CM), a pig breed from the Murcia region in the southeastern region of Spain, is a good model for endangered livestock populations. The remaining populations are bred on approximately 15 small farms, and no herdbook exists. To assess the genetic threats to the integrity and survival of the CM breed, and to aid in designing a conservation program, three genetic marker systems – microsatellites, SNPs and mtDNA – were applied across the majority of the total breeding stock. In addition, mtDNA and SNPs were genotyped in breeds that likely contributed genetically to the current CM gene pool. The analyses revealed the levels of genetic diversity within the range of other European local breeds (H_e = 0.53). However, when the eight farms that rear at least 10 CM pigs were independently analyzed, high levels of inbreeding were found in some. Despite the evidence for recent crossbreeding with commercial breeds on a few farms, the entire breeding stock remains readily identifiable as CM, facilitating the design of traceability assays. The genetic management of the breed is consistent with farm size, farm owner and presence of other pig breeds on the farm, demonstrating the highly ad hoc nature of current CM breeding. The results of genetic diversity and substructure of the entire breed, as well as admixture and crossbreeding obtained in the present study, provide a benchmark to develop future conservation strategies. Furthermore, this study demonstrates that identifying farm‐based practices and farm‐based breeding stocks can aid in the design of a sustainable breeding program for minority breeds.     

Genetic structure of Eurasian cattle (Bos taurus) based on microsatellites: clarification for their breed classification1

We pool three previously published data sets and present population genetic analyses of microsatellite variation in 48 Bos taurus cattle breeds from a wide range of geographical origins in Eurasia, mostly its northern territory. Bayesian model‐based clustering reveals six distinct clusters: besides a single‐population cluster of the Yakutian Cattle from Far Eastern Siberia and a cluster of breeds characteristic of an early origin, the other four major clusters largely correspond to previously defined morphological subgroups of Red Lowland, Lowland Black‐Pied, Longhorned Dairy and North European Polled cattle breeds. The results highlighted past expansion events of the productive breeds such as Danish Red, Angeln, Holstein‐Friesian and Ayrshire in northern and Eastern Europe. Based on genetic assignment of the breeds and the availability of breed information, we provide a preliminary classification of the five breeds that were to date undefined. Furthermore, in the analysis of molecular variance, despite some correspondence between geographical proximity and genetic similarity, the breed classification appears to be a better predictor of genetic structure in the cattle populations (the among‐group variance component: breed classification, 2.47%, P < 0.001; geographical division, 0.77%, P < 0.001).