Temporal variation in coat colour (genotypes) supports major changes in the Nordic cattle population after Iron Age

Variation in coat colour genotypes of archaeological cattle samples from Finland was studied by sequencing 69 base pairs of the extension locus (melanocortin 1‐receptor, MC1R) targeting both a transition and a deletion defining the three main alleles, such as dominant black (E^D), wild type (E⁺) and recessive red (e). The 69‐bp MC1R sequence was successfully analysed from 23 ancient (1000–1800 AD) samples. All three main alleles and genotype combinations were detected with allele frequencies of 0.26, 0.17 and 0.57 for E^D, E⁺ and e respectively. Recessive red and dominant black alleles were detected in both sexes. According to the best of our knowledge, this is the first ancient DNA study defining all three main MC1R alleles. Observed MC1R alleles are in agreement with calculated phenotype frequencies from historical sources. The division of ancient Finnish cattle population into modern Finnish breeds with settled colours was dated to the 20th century. From the existing genotyped populations in Europe (43 breeds, n = 2360), the closest match to ancient MC1R genotype frequencies was with the Norwegian native multicoloured breeds. In combined published genotype data of ancient (n = 147) and genotypes and phenotypes of modern Nordic cattle (n = 738), MC1R allele frequencies showed temporal changes similar to neutral mitochondrial DNA and Y‐chromosomal haplotypes analysed earlier. All three markers indicate major change in genotypes in Nordic cattle from the Late Iron Age to the Medieval period followed by slower change through the historical periods until the present.     

Molecular diversity and landscape genomics of the crop wild relative Triticum urartu across the Fertile Crescent

Modern plant breeding can benefit from the allelic variation that exists in natural populations of crop wild relatives that evolved under natural selection in varying pedoclimatic conditions. In this study, next‐generation sequencing was used to generate 1.3 million genome‐wide single nucleotide polymorphisms (SNPs) on ex situ collections of Triticum urartu L., the wild donor of the A^u subgenome of modern wheat. A set of 75 511 high‐quality SNPs were retained to describe 298 T. urartu accessions collected throughout the Fertile Crescent. Triticum urartu showed a complex pattern of genetic diversity, with two main genetic groups distributed sequentially from west to east. The incorporation of geographical information on sampling points showed that genetic diversity was correlated to the geographical distance (R² = 0.19) separating samples from Jordan and Lebanon, from Syria and southern Turkey, and from eastern Turkey, Iran and Iraq. The wild emmer genome was used to derive the physical positions of SNPs on the seven chromosomes of the A^u subgenome, allowing us to describe a relatively slow decay of linkage disequilibrium in the collection. Outlier loci were described on the basis of the geographic distribution of the T. urartu accessions, identifying a hotspot of directional selection on chromosome 4A. Bioclimatic variation was derived from grid data and related to allelic variation using a genome‐wide association approach, identifying several marker–environment associations (MEAs). Fifty‐seven MEAs were associated with altitude and temperature measures while 358 were associated with rainfall measures. The most significant MEAs and outlier loci were used to identify genomic loci with adaptive potential (some already reported in wheat), including dormancy and frost resistance loci. We advocate the application of genomics and landscape genomics on ex situ collections of crop wild relatives to efficiently identify promising alleles and genetic materials for incorporation into modern crop breeding.     

Detecting population structure and recent demographic history in endangered livestock breeds: the case of the Italian autochthonous donkeys

Since its domestication, about 5000 years ago, the donkey (Equus asinus) has been extensively used as a work or draft animal in agricultural activities and for the transportation of people and goods. In the last century, technology improvement and growing mechanization strongly affected agriculture and the management and use of this livestock species in the industrialized countries. Nowadays, the use of donkeys for work or transport has almost disappeared, together with the need for mules or hinny breeding. During the last five decades, Italian autochthonous donkey populations suffered from a severe reduction in population size, which led to the extinction of several breeds. At present, eight breeds remain, all classified by FAO as critically endangered or endangered: Asinara, Pantesco, Grigio Siciliano, Romagnolo, Amiatino, Sardo Grigio, Martina Franca, and Ragusano. To evaluate the extant genetic variability of Italian donkeys, we typed 16 microsatellite loci in 258 individuals from these breeds. The results highlighted moderate levels of inbreeding ( F _IS = 0.127) and a significant partition of genetic variation into breeds, as suggested by fixation index ( F _ST = 0.109) and analysis of molecular variance (10.86% of total variation assigned to the between‐breeds level) analyses. This was confirmed by a Bayesian clustering procedure that also highlighted a further partitioning at lower hierarchical levels corresponding to the farms of origin. This evidence suggests that an effective management strategy for Italian donkey populations should focus on breeds as conservation units. However, this requires a synergic management strategy at the farm level to maintain diversity and avoid inbreeding.     

Genetic structure and domestication footprints of the tusk,coat color,and ear morphology in East Chinese pigs

The domestication and artificial selection of wild boars have led to dramatic morphological and behavioral changes, especially in East Chinese (ECN) pigs. Here, we provide insights into the population structure and current genetic diversity of representative ECN pig breeds. We identify a 500-kb region containing six tooth development-relevant genes with almost completely different haplotypes between ECN pigs and Chinese wild boars or European domestic pigs. Notably, the c.195A>G missense mutation in exon 2 of AMBN may cause alterations in its protein structure associated with tusk degradation in ECN pigs. In addition, ESR1 may play an important role in the reproductive performance of ECN pigs. A major haplotype of the large lop ear-related MSRB3 gene and eight alleles in the deafness-related GRM7 gene may affect ear morphology and hearing in ECN pigs. Interestingly, we find that the two-end black (TEB) coat color in Jinhua pigs is most likely caused by EDNRB with genetic mechanisms different from other Chinese TEB pigs. This study identifies key loci that may be artificially selected in Chinese native pigs related to the tusk, coat color, and ear morphology, thus providing new insights into the genetic mechanisms of domesticated pigs.     

Polymorphism pattern at a miniature inverted‐repeat transposable element locus downstream of the domestication gene Teosinte‐branched1 in wild and domesticated pearl millet

Unravelling the mechanisms involved in adaptation to understand plant morphological evolution is a challenging goal. For crop species, identification of molecular causal polymorphisms involved in domestication traits is central to this issue. Pearl millet, a domesticated grass mostly found in semi‐arid areas of Africa and India, is an interesting model to address this topic: the domesticated form shares common derived phenotypes with some other cereals such as a decreased ability to develop basal and axillary branches in comparison with the wild phenotype. Two recent studies have shown that the orthologue of the maize gene Teosinte‐Branched1 in pearl millet (PgTb1) was probably involved in branching evolution during domestication and that a miniature inverted‐repeat transposable element (MITE) of the Tuareg family was inserted in the 3′ untranslated region of PgTb1. For a set of 35 wild and domesticated populations, we compared the polymorphism patterns at this MITE and at microsatellite loci. The Tuareg insertion was nearly absent in the wild populations, whereas a strong longitudinal frequency cline was observed in the domesticated populations. The geographical pattern revealed by neutral microsatellite loci clearly demonstrated that isolation by distance does not account for the existence of this cline. However, comparison of population differentiation at the microsatellite and the MITE loci and analyses of the nucleotide polymorphism pattern in the downstream region of PgTb1 did not show evidence that the cline at the MITE locus has been shaped by selection, suggesting the implication of a neutral process. Alternative hypotheses are discussed.     

Mutations in ASIP and MC1R: dominant black and recessive black alleles segregate in native Swedish sheep populations

By studying genes associated with coat colour, we can understand the role of these genes in pigmentation but also gain insight into selection history. North European short‐tailed sheep, including Swedish breeds, have variation in their coat colour, making them good models to expand current knowledge of mutations associated with coat colour in sheep. We studied ASIP and MC1R, two genes with known roles in pigmentation, and their association with black coat colour. We did this by sequencing the coding regions of ASIP in 149 animals and MC1R in 129 animals from seven native Swedish sheep breeds in individuals with black, white or grey fleece. Previously known mutations in ASIP [recessive black allele: g.100_105del (D₅) and/or g.5172T>A] were associated with black coat colour in Klövsjö and Roslag sheep breeds and mutations in both ASIP and MC1R (dominant black allele: c.218T>A and/or c.361G>A) were associated with black coat colour in Swedish Finewool. In Gotland, Gute, Värmland and Helsinge sheep breeds, coat colour inheritance was more complex: only 11 of 16 individuals with black fleece had genotypes that could explain their black colour. These breeds have grey individuals in their populations, and grey is believed to be a result of mutations and allelic copy number variation within the ASIP duplication, which could be a possible explanation for the lack of a clear inheritance pattern in these breeds. Finally, we found a novel missense mutation in MC1R (c.452G>A) in Gotland, Gute and Värmland sheep and evidence of a duplication of MC1R in Gotland sheep.     

Not another type of potato: MC1R and the russet coloration of Burmese cats

The Burmese is a breed of domestic cat that originated in Southeast Asia and was further developed in the United States. Variants in melanocortin 1 receptor (MC1R) causes common coat colour phenotypes in a variety of mammalian species but only limited colour variation in the domestic cat. Known as the extension (E) locus, melanocortin 1 receptor (MC1R) interacts with the agouti locus to produce the eumelanin and pheomelanin pigments. Recently, a novel reddish coloration, which is termed russet, was identified in the Burmese cat breed. Because this russet Burmese coloration changes with aging, MC1R was suggested as candidate gene. The similar colouration in specific lineages of Norwegian Forest cat known as amber (e) (c.250G>A; p.Asp84Asn) was excluded for this Burmese phenotype. The complete 954‐bp coding region of MC1R was directly sequenced in russet Burmese and suspected carriers. A 3‐bp deletion (c.439_441del) associated with the deletion of a phenyalanine (p.Phe146del) in the protein sequence was identified. All russet coloured cats were homozygous for the variant, and all obligate carriers were heterozygous, confirming that the deletion segregated concordantly with colouring in Burmese cats from the New Zealand foundation lineage. The variant was not identified in 442 cats from 26 different breeds and random‐bred cats. Twenty‐six Burmese from the USA did not have the variant. This MC1R variant defines a unique coloration and the second breed‐specific MC1R variant in cats. The interactions of the two recessive feline MC1R alleles (E  >  e, e^r) is unknown.     

High domestic pig contribution to the local gene pool of free-living European wild boar: a case study in Poland

Rates of hybridization between wild and domesticated animals appear to be increasing worldwide. Recent results suggest that genetic introgression from domestic swine into European wild boar is much more common in local populations than expected, based on pan-European studies. Thus, we screened the genetic purity of 265 free-living wild boars from two hunting areas in Poland by genotyping the melanocortin receptor 1 gene (MC1R) for polymorphism. Unexpectedly, high numbers of individuals with domestic genes (24%) were identified. This suggests that mixed ancestry may be common in Polish wild boar. Among admixed individuals, backcrosses with domestic pig and/or introgressed wild boars were detected (2%). Multiple commercial domestic pig breeds are possibly involved in the introgression observed in the study populations. In addition, the absence of significant differences in the frequency of wild-type allele among two hunting areas suggests high dispersal of individuals and gene flow among populations. We conclude that further study is needed to better understand the mechanisms and sources of introgression in wild boars in Poland.     

Largely unlinked gene sets targeted by selection for domestication syndrome phenotypes in maize and sorghum

The domestication of diverse grain crops from wild grasses was a result of artificial selection for a suite of overlapping traits producing changes referred to in aggregate as ‘domestication syndrome’. Parallel phenotypic change can be accomplished by either selection on orthologous genes or selection on non‐orthologous genes with parallel phenotypic effects. To determine how often artificial selection for domestication traits in the grasses targeted orthologous genes, we employed resequencing data from wild and domesticated accessions of Zea (maize) and Sorghum (sorghum). Many ‘classic’ domestication genes identified through quantitative trait locus mapping in populations resulting from wild/domesticated crosses indeed show signatures of parallel selection in both maize and sorghum. However, the overall number of genes showing signatures of parallel selection in both species is not significantly different from that expected by chance. This suggests that while a small number of genes will extremely large phenotypic effects have been targeted repeatedly by artificial selection during domestication, the optimization part of domestication targeted small and largely non‐overlapping subsets of all possible genes which could produce equivalent phenotypic alterations.     

Detection of hybrids between wild boars (Sus scrofa scrofa) and domestic pigs (Sus scrofa f. domestica) in Greece,using the PCR-RFLP method on melanocortin-1 receptor (MC1R) mutations

The melanocortin-1 receptor (MC1R) regulates melanogenesis in mammals within the mammalian melanocyte and the hair follicle. Common variations (polymorphisms) in the MC1R gene are associated with normal differences in skin and hair colour. So far, a unique MC1R allele (E+) has been identified in European wild boar (Sus scrofa scrofa), associated with the wild-type coat colour (variable shades of brown) that is not found in any of the domestic breeds. In addition, a series of alleles found in pigs, some of which observed only in particular breeds, have been proposed as markers in breed traceability systems. The current study is an attempt to detect possible hybrids between wild boars and domestic pig breeds as well as to identify races of pig that are not purebred. For this purpose, wild boars were analysed against Large White pigs, applying the PCR-restriction fragment length polymorphism (RFLP) method. A high percentage (16.7%) of hybrids was detected within a breeding station compared with the percentage of hybrids within the populations of free-ranging wild boar (5.0%). These results should be taken into consideration for future restocking operations to avoid the chance of outbreeding depression, which is more intense when local populations are introgressed by gene pools from domesticated, usually inbred, animals.     

The genetics of brown coat color and white spotting in domestic yaks (Bos grunniens)

Domestic yaks (Bos grunniens) exhibit two major coat color variations: a brown vs. wild‐type black pigmentation and a white spotting vs. wild‐type solid color pattern. The genetic basis for these variations in color and distribution remains largely unknown and may be complicated by a breeding history involving hybridization between yaks and cattle. Here, we investigated 92 domestic yaks from China using a candidate gene approach. Sequence variations in MC1R, PMEL and TYRP1 were surveyed in brown yaks; TYRP1 was unassociated with the coloration and excluded. Recessive mutations from MC1R, or p.Gln34*, p.Met73Leu and possibly p.Arg142Pro, are reported in bovids for the first time and accounted for approximately 40% of the brown yaks in this study. The remaining 60% of brown individuals correlated with a cattle‐derived deletion mutation from PMEL (p.Leu18del) in a dominant manner. Degrees of white spotting found in yaks vary from color sidedness and white face, to completely white. After examining the candidate gene KIT, we suggest that color‐sided and all‐white yaks are caused by the serial translations of KIT (Cs₆ or Cs₂₉) as reported for cattle. The white‐faced phenotype in yaks is associated with the KIT haplotype S^wf. All KIT mutations underlying the serial phenotypes of white spotting in yaks are identical to those in cattle, indicating that cattle are the likely source of white spotting in yaks. Our results reveal the complex genetic origins of domestic yak coat color as either native in yaks through evolution and domestication or as introduced from cattle through interspecific hybridization.     

Evolutionary history of almond tree domestication in the Mediterranean basin

Genetic diversity of contemporary domesticated species is shaped by both natural and human‐driven processes. However, until now, little is known about how domestication has imprinted the variation of fruit tree species. In this study, we reconstruct the recent evolutionary history of the domesticated almond tree, Prunus dulcis, around the Mediterranean basin, using a combination of nuclear and chloroplast microsatellites [i.e. simple sequence repeat (SSRs)] to investigate patterns of genetic diversity. Whereas conservative chloroplast SSRs show a widespread haplotype and rare locally distributed variants, nuclear SSRs show a pattern of isolation by distance with clines of diversity from the East to the West of the Mediterranean basin, while Bayesian genetic clustering reveals a substantial longitudinal genetic structure. Both kinds of markers thus support a single domestication event, in the eastern side of the Mediterranean basin. In addition, model‐based estimation of the timing of genetic divergence among those clusters is estimated sometime during the Holocene, a result that is compatible with human‐mediated dispersal of almond tree out of its centre of origin. Still, the detection of region‐specific alleles suggests that gene flow from relictual wild preglacial populations (in North Africa) or from wild counterparts (in the Near East) could account for a fraction of the diversity observed.     

Next generation semiconductor based sequencing of bitter taste receptor genes in different pig populations and association analysis using a selective DNA pool‐seq approach

Taste perception in animals affects feed intake and may influence production traits. In particular, bitter is sensed by receptors encoded by the family of TAS2R genes. In this research, using a DNA pool‐seq approach coupled with next generation semiconductor based target resequencing, we analysed nine porcine TAS2R genes (TAS2R1, TAS2R3, TAS2R4, TAS2R7, TAS2R9, TAS2R10, TAS2R16, TAS2R38 and TAS2R39) to identify variability and, at the same time, estimate single nucleotide polymorphism (SNP) allele frequencies in several populations and testing differences in an association analysis. Equimolar DNA pools were prepared for five pig breeds (Italian Duroc, Italian Landrace, Pietrain, Meishan and Casertana) and wild boars (5–10 individuals each) and for two groups of Italian Large White pigs with extreme and divergent back fat thickness (50 + 50 pigs). About 1.8 million reads were obtained by sequencing amplicons generated from these pools. A total of 125 SNPs were identified, of which 37 were missense mutations. Three of them (p.Ile53Phe and p.Trp85Leu in TAS2R4; p.Leu37Ser in TAS2R39) could have important effects on the function of these bitter taste receptors, based on in silico predictions. Variability in wild boars seems lower than that in domestic breeds potentially as a result of selective pressure in the wild towards defensive bitter taste perception. Three SNPs in TAS2R38 and TAS2R39 were significantly associated with back fat thickness. These results may be important to understand the complexity of taste perception and their associated effects that could be useful to develop nutrigenetic approaches in pig breeding and nutrition.     

Evolution of MHC class II SLA‐DRB1 locus in the Croatian wild boar (Sus scrofa) implies duplication and weak signals of positive selection

The wild boar is an ancestor of the domestic pig and an important game species with the widest geographical range of all ungulates. Although a large amount of data are available on major histocompatibility complex (MHC) variability in domestic pigs, only a few studies have been performed on wild boars. Due to their crucial role in appropriate immune responses and extreme polymorphism, MHC genes represent some of the best candidates for studying the processes of adaptive evolution. Here, we present the results on the variability and evolution of the entire MHC class II SLA‐DRB1 locus exon 2 in 133 wild boars from Croatia. Using direct sequencing and cloning methods, we identified 20 SLA‐DRB1 alleles, including eight new variants, with notable divergence. In some individuals, we documented functional locus duplication, and SLA‐DRB1*04:10 was identified as the allele involved in the duplication. The expression of a duplicated locus was confirmed by cloning and sequencing cDNA‐derived amplicons. Based on individual genotypes, we were able to assume that alleles SLA‐DRB1*04:10 and SLA‐DRB1*06:07 are linked as an allelic combination that co‐evolves as a two‐locus haplotype. Our investigation of evolutionary processes at the SLA‐DRB1 locus confirmed the role of intralocus recombination in generating allelic variability, whereas tests of positive selection based on the dN/dS (non‐synonymous/synonymous substitution rate ratio) test revealed atypically weak and ambiguous signals.     

Mitochondrial lineages reveal intense gene flow between Iberian wild boars and South Iberian pig breeds

The phylogeography of wild boars (WB) and domestic pigs (Sus scrofa) has contributed important insights into where and when domestication occurred. The geographic distribution of two core haplotypes (E1a and E1c) of the main European phylogenetic clade suggests that Central Europe was an early domestication centre, although the complexity of the pattern does not exclude the possibility that multiple domestication events occurred in different regions. To investigate the relationships among WB and domestic pig breeds in Iberia, a fragment of the mitochondrial DNA control region from a large sample (n = 409) of WB and local pig breeds was co‐analysed with published sequences from other European populations. The Iberian sample revealed a high frequency of a sub‐cluster (E1c) of the European haplogroup E1 in 77% of total Iberian samples, 96% of WB, 90% of Alentejano (Portugal) and 87% of Iberian breed pigs (Spain; Black Hairy, Black Hairless and Red varieties). Low genetic distance (F’_ST = 0.105) was observed between Alentejano (Portugal) and Iberian breed pigs (Spain). Alentejano and Iberian breed pigs showed low genetic distances to both Iberian and Central European WB (average F’_ST = 0.345 and 0.215, respectively). This pattern suggests that early pig husbandry in the Iberian Peninsula did not solely rely on imported Central European stock, but also included the recruitment of local WB.     

Global investigation of the co‐evolution of MIRNA genes and microRNA targets during soybean domestication

Although the selection of coding genes during plant domestication has been well studied, the evolution of MIRNA genes (MIRs) and the interaction between microRNAs (miRNAs) and their targets in this process are poorly understood. Here, we present a genome‐wide survey of the selection of MIRs and miRNA targets during soybean domestication and improvement. Our results suggest that, overall, MIRs have higher evolutionary rates than miRNA targets. Nonetheless, they do demonstrate certain similar evolutionary patterns during soybean domestication: MIRs and miRNA targets with high expression and duplication status, and with greater numbers of partners, exhibit lower nucleotide divergence than their counterparts without these characteristics, suggesting that expression level, duplication status, and miRNA–target interaction are essential for evolution of MIRs and miRNA targets. Further investigation revealed that miRNA–target pairs that are subjected to strong purifying selection have greater similarities than those that exhibited genetic diversity. Moreover, mediated by domestication and improvement, the similarities of a large number of miRNA–target pairs in cultivated soybean populations were increased compared to those in wild soybeans, whereas a small number of miRNA–target pairs exhibited decreased similarity, which may be associated with the adoption of particular domestication traits. Taken together, our results shed light on the co‐evolution of MIRs and miRNA targets during soybean domestication.     

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.     

Genetic diversity of Chinese cattle revealed by Y‐SNP and Y‐STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y‐SNPs (UTY19 and ZFY10) and Y‐STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y‐chromosome haplotypes that represented three haplogroups. Y2‐104‐158 and Y2‐102‐158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3‐88‐156 predominates in southern China. Haplotypes Y2‐108‐158, Y2‐110‐158, Y2‐112‐158 and Y3‐92‐156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.     

Evaluation of genetic diversity and population structure of west‐central Indian cattle breeds

Evaluations of genetic diversity in domestic livestock populations are necessary to implement region‐specific conservation measures. We determined the genetic diversity and evolutionary relationships among eight geographically and phenotypically diverse cattle breeds indigenous to west‐central India by genotyping these animals for 22 microsatellite loci. A total of 326 alleles were detected, and the expected heterozygosity ranged from 0.614 (Kenkatha) to 0.701 (Dangi). The mean number of alleles among the cattle breeds ranged from 7.182 (Khillar) to 9.409 (Gaolao). There were abundant genetic variations displayed within breeds, and the genetic differentiation was also high between the Indian cattle breeds, which displayed 15.9% of the total genetic differentiation among the different breeds. The genetic differentiation (pairwise F_ST) among the eight Indian breeds varied from 0.0126 for the Kankrej–Malvi pair to 0.2667 for Khillar–Kenkatha pair. The phylogeny, principal components analysis, and structure analysis further supported close grouping of Kankrej, Malvi, Nimari and Gir; Gaolao and Kenkatha, whereas Dangi and Khillar remained at distance from other breeds.