Estimates of autozygosity derived from runs of homozygosity: empirical evidence from selected cattle populations

Using genome‐wide SNP data, we calculated genomic inbreeding coefficients (F_ROH > _1 Mb, F_{ROH > 2 Mb}, F_{ROH > 8 Mb} and F_{ROH > 16 Mb}) derived from runs of homozygosity (ROH) of different lengths (>1, >2, >8 and > 16 Mb) as well as from levels of homozygosity (F_HOM). We compared these values of inbreeding coefficients with those calculated from pedigrees (F_PED) of 1422 bulls comprising Brown Swiss (304), Fleckvieh (502), Norwegian Red (499) and Tyrol Grey (117) cattle breeds. For all four breeds, population inbreeding levels estimated by the genomic inbreeding coefficients F_{ROH > 8 Mb} and F_{ROH > 16 Mb} were similar to the levels estimated from pedigrees. The lowest values were obtained for Fleckvieh (F_PED = 0.014, F_{ROH > 8 Mb} = 0.019 and F_{ROH > 16 Mb} = 0.008); the highest, for Brown Swiss (F_PED = 0.048, F_{ROH > 8 Mb} = 0.074 and F_{ROH > 16 Mb} = 0.037). In contrast, inbreeding estimates based on the genomic coefficients F_{ROH > 1 Mb} and F_{ROH > 2 Mb} were considerably higher than pedigree‐derived estimates. Standard deviations of genomic inbreeding coefficients were, on average, 1.3–1.7‐fold higher than those obtained from pedigrees. Pearson correlations between genomic and pedigree inbreeding coefficients ranged from 0.50 to 0.62 in Norwegian Red (lowest correlations) and from 0.64 to 0.72 in Tyrol Grey (highest correlations). We conclude that the proportion of the genome present in ROH provides a good indication of inbreeding levels and that analysis based on ROH length can indicate the relative amounts of autozygosity due to recent and remote ancestors.     

Relationship between inbreeding and milk production traits in two Italian dairy sheep breeds

The effects of inbreeding in livestock species breeds have been well documented and they have a negative impact on profitability. The objective of this study was to evaluate the levels of inbreeding in Sarda (SAR, n = 785) and Valle del Belice (VdB, n = 473) dairy sheep breeds and their impact on milk production traits. Two inbreeding coefficients (F) were estimated: using pedigree (F_PED), or runs of homozygosity (ROH; F_ROH) at different minimum ROH lengths and different ROH classes. After the quality control, 38,779 single nucleotide polymorphisms remained for further analyses. A mixed-linear model was used to evaluate the impact of inbreeding coefficients on production traits within each breed. VdB showed higher inbreeding coefficients compared to SAR, with both breeds showing lower estimates as the minimum ROH length increased. Significant inbreeding depression was found only for milk yield, with a loss of around 7 g/day (for SAR) and 9 g/day (VdB) for a 1% increase of F_ROH. The present study confirms how the use of genomic information can be used to manage intra-breed diversity and to calculate the effects of inbreeding on phenotypic traits.     

Estimation of inbreeding depression on female fertility in the Finnish Ayrshire population

Single nucleotide polymorphism (SNP) data enable the estimation of inbreeding at the genome level. In this study, we estimated inbreeding levels for 19,075 Finnish Ayrshire cows genotyped with a low‐density SNP panel (8K). The genotypes were imputed to 50K density, and after quality control, 39,144 SNPs remained for the analysis. Inbreeding coefficients were estimated for each animal based on the percentage of homozygous SNPs (F_PH), runs of homozygosity (F_ROH) and pedigree (F_PED). Phenotypic records were available for 13,712 animals including non‐return rate (NRR), number of inseminations (AIS) and interval from first to last insemination (IFL) for heifers and up to three parities for cows, as well as interval from calving to first insemination (ICF) for cows. Average F_PED was 0.02, F_ROH 0.06 and F_PH 0.63. A correlation of 0.71 was found between F_PED and F_ROH, 0.66 between F_PED and F_PH and 0.94 between F_ROH and F_PH. Pedigree‐based inbreeding coefficients did not show inbreeding depression in any of the traits. However, when F_ROH or F_PH was used as a covariate, significant inbreeding depression was observed; a 10% increase in F_ROH was associated with 5 days longer IFL0 and IFL1, 2 weeks longer IFL3 and 3 days longer ICF2 compared to non‐inbred cows.     

Estimation of inbreeding and effective population size of full‐blood wagyu cattle registered with the American Wagyu Cattle Association

The objective of this research was to examine the population structure of full‐blood (100%) Wagyu cattle registered in the United States with the American Wagyu Association, with the aim of estimating and comparing the levels of inbreeding from both pedigree and genotypic data. A total of 4132 full‐blood Wagyu cattle pedigrees were assessed and used to compute the inbreeding coefficients (F_IT and F_ST) and the effective population size (N_e) from pedigree data for the period 1994 to 2011. In addition to pedigree analysis, 47 full‐blood Wagyu cattle representing eight prominent sire lines in the American Wagyu cattle population were genotyped using the Illumina BovineSNP50 BeadChip. Genotypic data were then used to estimate genomic inbreeding coefficients (F_ROH) by calculating runs of homozygosity. The mean inbreeding coefficient based on the pedigree data was estimated at 4.80%. The effective population size averaged 17 between the years 1994 and 2011 with an increase of 42.9 in 2000 and a drop of 1.8 in 2011. Examination of the runs of homozygosity revealed that the 47 Wagyu cattle from the eight prominent sire lines had a mean genomic inbreeding coefficient (F_ROH) estimated at 9.08% compared to a mean inbreeding coefficient based on pedigree data of 4.8%. These data suggest that the mean genotype inbreeding coefficient of full‐blood Wagyu cattle exceeds the inbreeding coefficient identified by pedigree. Inbreeding has increased slowly at a rate of 0.03% per year over the past 17 years. Wagyu breeders should continue to utilize many sires from divergent lines and consider outcrossing to other breeds to enhance genetic diversity and minimize the adverse effects of inbreeding in Wagyu.     

Genome-wide identification of runs of homozygosity islands in the Gyr breed (Bos indicus)

Runs of homozygosity (ROH) are contiguous homozygous regions of the genome. These regions can be used to identify genes associated with traits of economic interest, as well as inbreeding levels. The aim of the present study was to analyse the length and distribution of ROH islands in Gyr cattle and to identify genes within these regions. A population of 173 animals selected for beef production and a population of 291 animals selected for dairy production were used. Differences in the number of short ROH (ROH_{1-2 Mb}) were observed between the two populations, while the number of long ROH (ROH_{>16 Mb}) was similar. ROH islands with the highest incidences (>0.50) overlapped in several segments of the genome in the two populations. The genes identified were associated with milk production, growth, reproduction, immune response and resistance traits. Our results contribute to the understanding of how selection can shape the distribution of ROH and ROH islands within the same breed when animals are selected for different purposes such as dairy or beef production.     

Genome‐wide analysis of homozygosity regions in european simmental bulls

The study of Runs of Homozygosity (ROH) is a useful approach for the characterization of the genome of livestock populations. Due to their high relationship with autozygosity, ROH allow to make inference about population genetic history, to estimate the level of inbreeding, to assess within breed heterogeneity and to detect the footprints of selection on livestock genomes. Aim of this study was to investigate the distribution of runs of homozygosity in bulls belonging to five European Simmental populations and to assess the relationship between three production traits (milk yield, fat and protein contents) and autozygosity. ROH count, distribution and ROH‐based coefficient of inbreeding (F_ROH) were calculated for 3,845 Simmental bulls of five different European countries: Austria (AT), Switzerland (CH), Czech Republic (CZ), Germany (DE) and Italy (IT). Average values of ROH number per animal, and total genome length covered by ROH were 77.8 ± 20.7 and 205 ± 74.4 Mb, respectively. Bulls from AT, DE and IT exhibited similar ROH characteristics. Swiss animals showed the highest (12.6%), while CZ the lowest (4.6%) F_ROH coefficient. The relationship between ROH occurrence and milk production traits was investigated through a genome‐wide ROH‐traits association analysis (GWRA). A total of 34 regions previously associated with milk traits (yield and/or composition) were identified by GWRA. Results of the present research highlight a mixed genetic background in the 5 European Simmental populations, with the possible presence of three subgroups. Moreover, a strong relationship between autozygosity and production traits has been detected.     

Detection of runs of homozygosity in conserved and commercial pig breeds in Poland

Runs of homozygosity (ROH) are continuous segments of the genome that arose as a result of inbreeding, resulting in the inheritance of identical haplotypes from both parents who shared a common ancestor. In the present study, we performed a detailed characterization and comparison of ROH in four pig breeds, including intensively selected Polish Landrace as well as native unselected animals of Puławska and two Złotnicka breeds (White and Spotted). We used a medium-density PorcineSNP60 BeadChip assay (Illumina) and cgaTOH software to detect ROH covering a minimum of 30 adjacent SNPs and maintaining a size over 1 Mb. By analysing ROH distribution and frequency across the genome, we also identified genomic regions with high ROH frequency (so-called “ROH hotspots”). The obtained results showed that the analysed conserved breeds were characterized by a higher ROH span and higher ROH-based inbreeding coefficients (F_ROH), which likely result from past population bottlenecks, increasing the overall inbreeding level within these populations. The analysis of ROH distribution across the genomes revealed the presence of both shared and breed-specific ROH hotspots. These hotspots, presumably representing genome regions under selection, overlapped with a variety of genes associated with processes connected with immune system functioning, reproduction, glucose homeostasis and metabolism. The genome regions with ROH hotspots overlapping in all analysed populations, located on SSC4 (51.9–55.9 Mb) and 13 (92.6–97.8 Mb), covered thirty-one different genes, including MMP16, SLC7A13, ATP6V0D2, CNGB3, WWiP1, RiMDN1 and CPNE3. These genes are primarily associated with biological regulation and metabolism, processes that could be responsible for the variety of the selected production and functional features.     

The assessment of genetic diversity within and among the eight subpopulations of Japanese Black cattle using 52 microsatellite markers

Japanese Black cattle are at risk for genetic homogeneity due to intensive use of a few sires. Therefore, assessment of the actual genetic diversity of this breed is important for future breeding plans. In the present study, we investigated the genetic diversity within and among eight subpopulations of Japanese Black cattle using 52 microsatellite markers. The parameters for genetic diversity of Japanese Black cattle were comparable to those of other cattle breeds, suggesting that the relatively high genetic diversity of the breed. However, upon comparison among the eight subpopulations, the Hyogo subpopulation showed markedly low genetic diversity. The results of the pairwise F_ST values, phylogenetic network and structure analysis indicated that the Hyogo population has remarkably high level of genetic differentiation from other populations, while Yamagata, Niigata, Hiroshima and Kagawa populations have low levels of genetic differentiation. Furthermore, multidimensional scaling plots indicated that individuals in some subpopulations were separated from individuals in the other subpopulations. We conclude that while the overall genetic diversity of Japanese Black cattle is still maintained at a relatively high level, that of a particular subpopulation is significantly reduced, and therefore the effective population size of the breed needs to be controlled by correct mating strategies.     

Genome-wide scan for runs of homozygosity in the composite Montana Tropical® beef cattle

The aim of this study was to assess the distribution of runs of homozygosity (ROH) and autozygosity islands in the composite Montana Tropical^® beef cattle to explore hotspot regions which could better characterize the different biological types within the composite breed. Montana animals (n = 1,436) were genotyped with the GGP-LD BeadChip (~30,000 markers). ROH was identified in every individual using the plink v1.90 software. Medium and long ROH prevailed in the genome, which accounted for approximately 74% of all ROH detected. On an average, 2.0% of the genome was within ROH, agreeing with the pedigree-based inbreeding coefficient. The Montana cattle with a higher proportion of productive breed types showed the highest number of autozygosity islands (n = 17), followed by those with a higher proportion of breeds adapted to tropical environments (n = 15). Enriched terms (p < .05) associated with the immune and inflammatory response, homeostasis, reproduction, mineral absorption, and lipid metabolism were described within the autozygosity islands. In this regard, over-represented GO terms and KEGG pathways described in this population may play a key role in providing information to explore the genetic and biological mechanisms together with the genomic regions underlying each biological type that favoured their optimal performance ability in tropical and subtropical regions.     

Evaluation of inbreeding and genetic diversity in Japanese Shorthorn cattle by pedigree analysis

The Japanese Shorthorn is a Japanese Wagyu breed maintained at a small population size. We assessed the degree of inbreeding and genetic diversity among Japanese Shorthorn cattle using pedigree analysis. We analyzed the pedigree records of registered Japanese Shorthorn born between 1980 and 2018, after evaluating the pedigree completeness. The average of the actual inbreeding coefficients increased at the same rates annually from approximately 1.5% in 1980 to 4.2% in 2018 and was higher than the expected inbreeding coefficients over time. The effective population size based on the individual coancestry rate largely decreased from 127.8 in 1980 to 82.6 in 1999, and then remained almost constant at approximately 90. Three effective numbers of ancestors decreased over time until 1995, then remained almost constant. In particular, the effective number of founder genomes (N_ge) decreased from 43.8 in 1980 to 11.9 in 2018. The index of genetic diversity based on N_ge decreased from 0.99 in 1980 to 0.96 in 2018 due to genetic drift in non-founder generations. Changes in inbreeding and genetic diversity parameters were similar between Japanese Shorthorn and other Japanese Wagyu breeds, but the magnitude of the changes was lower in the Japanese Shorthorn.     

Anaplasma marginale infection in a Japanese Black cow 13 years after eradication of Rhipicephalus (Boophilus) microplus in Okinawa, Japan

In October 2007, a 15-year-old Japanese Black cow on Ishigaki Island, Okinawa, Japan, was diagnosed with Anaplasma marginale infection based on clinical symptoms, blood examination, smear observation, 16S rRNA and groEL gene sequence analysis, and the result of a CF test. The cow was introduced into the farm from mainland Japan as a calf in 1993, one year before the eradication of Rhipicephalus (Boophilus) microplus, the main vector of A. marginale in Okinawa Prefecture. It is possible that the cow was first infected with A. marginale as a calf in Ishigaki Island and had been persistently infected since then. This is the first reported clinical case of A. marginale infection of cattle since the eradication of R. microplus in Okinawa Prefecture. Additional analysis of major surface protein 1α amino acid sequences revealed that the A. marginale Okinawa strain presented four new repeat forms which were not seen in other strains. This indicates that the Okinawa strain may be a unique geographical variant of A. marginale.     

Genotypes and allele frequencies of buried SNPs in a bovine single‐nucleotide polymorphism array in Japanese Black cattle

Single nucleotide polymorphism (SNP) arrays are widely used for genetic and genomic analyses in cattle breeding; thus, data derived from SNP arrays have accumulated on a large scale nationwide. Commercial SNP arrays contain a considerable number of unassigned SNPs on the chromosome/position on the genome; these SNPs are excluded in subsequent analyses. Notably, the position‐unassigned SNPs, or “buried SNPs” include some of the markers associated with genetic disease. In this study, we identified the position of buried SNPs using the Basic Local Alignment Search Tool against the surrounding sequences and characterized the relationship between SNPs and genetic diseases in Online Mendelian Inheritance in Animals based on the genomic position. We determined the position of 285 buried SNPs on the genome and surveyed the genotype and allele frequencies of these SNPs in 5,955 individual Japanese Black cattle. Eleven SNPs associated with genetic disease, which contained five buried SNPs, were found in the population with the risk allele frequency ranging from 0.00008396 to 0.46. These results indicate that buried SNPs in the bovine SNP array can be utilized to identify associations with genetic disorders from large scale accumulated SNP genotype data in Japanese Black cattle.     

Association of a single nucleotide polymorphism in ribosomal protein L27a gene with marbling in Japanese Black beef cattle

Marbling, defined by the amount and distribution of intramuscular fat, is an economically important trait of beef cattle in Japan. The c2‐11#2 expressed sequence tag (EST) has been previously shown to possess expression difference in musculus longissimus muscle between low‐marbled and high‐marbled steer groups, and to be located within genomic region of a quantitative trait locus for marbling. Thus, the ribosomal protein L27a (RPL27A) gene containing the c2‐11#2 EST sequence was considered as a positional candidate for the gene responsible for marbling. In the present study, a single nucleotide polymorphism (SNP) in the promoter region of the RPL27A, referred to as g.3109537C>T, was detected between the 2 steer groups. The SNP was associated with the predicted breeding value for beef marbling standard number by the analyses using Japanese Black beef cattle population. The effect of genotypes of the SNP on the predicted breeding value for subcutaneous fat thickness was not statistically significant. These findings suggest that the RPL27A SNP may be useful for effective marker‐assisted selection to increase the levels of marbling in Japanese Black beef cattle.     

SNP芯片评估柯尔克孜羊群体遗传多样性和遗传结构

旨在了解柯尔克孜羊群体的遗传多样性和遗传结构,有效地保护和利用其遗传资源。本研究利用绵羊SNP 50K v3芯片检测61只柯尔克孜种羊(31只公羊、30只母羊)个体的单核苷酸多态性(single nucleotide polymorphism, SNP);Plink(V1.90)软件对数据进行质控,计算群体有效含量、多态标记的比例、观测杂合度、期望杂合度、多态信息含量、有效等位基因数、最小等位基因频率,分析群体的遗传多样性;Plink计算连续性纯合片段(runs of homozygosity, ROH)和近交系数F_ROH;构建状态同源距离矩阵(identical by state, IBS),并采用Gmatrix软件构建G矩阵,解析柯尔克孜羊群的遗传距离和亲缘关系;使用Mega X软件构建种公羊进化树,分析群体家系结构。结果显示,61只柯尔克孜羊共得到64 734个SNPs标记,通过质检的SNPs为56 763个;平均多态信息含量为0.273±0.112,平均观察杂合度和平均期望杂合度分别为0.368±0.140和0.368±0.130,平均最小等位基因频率为0...     

A genome‐wide association study for fat‐related traits computed by image analysis in Japanese Black cattle

The objective of this study was to identify genomic regions associated with fat‐related traits using a Japanese Black cattle population in Hyogo. From 1836 animals, those with high or low values were selected on the basis of corrected phenotype and then pooled into high and low groups (n = 100 each), respectively. DNA pool‐based genome‐wide association study (GWAS) was performed using Illumina BovineSNP50 BeadChip v2 with three replicate assays for each pooled sample. GWAS detected that two single nucleotide polymorphisms (SNPs) on BTA7 (ARS‐BFGL‐NGS‐35463 and Hapmap23838‐BTA‐163815) and one SNP on BTA12 (ARS‐BFGL‐NGS‐2915) significantly affected fat percentage (FAR). The significance of ARS‐BFGL‐NGS‐35463 on BTA7 was confirmed by individual genotyping in all pooled samples. Moreover, association analysis between SNP and FAR in 803 Japanese Black cattle revealed a significant effect of SNP on FAR. Thus, further investigation of these regions is required to identify FAR‐associated genes and mutations, which can lead to the development of DNA markers for marker‐assisted selection for the genetic improvement of beef quality.     

Monitoring of genetic diversity in the Japanese Black cattle population by the use of pedigree information

The gene pool of the Japanese Black cattle has been completely closed to foreign breeds during the last 100 years. Genetic diversity of the Japanese Black cattle from 1960 to 2000 was monitored with three estimates of effective number of ancestors. Founder genome equivalent (N_ge) accounts for all the causes of reduction of diversity. Effective number of founders (N_ef) and non‐founders (N_enf) explain reduced diversity because of unequal genetic contributions of founders and random genetic drift in non‐founders, respectively. Further examination using gene dropping simulation was conducted to obtain information on survival of founder alleles. Unique founder alleles were dropped down along the actual pedigree with Monte Carlo procedure following Mendelian segregation rules, and generated genotypes of all the current live animals (612 959 heads). Pedigree records consisted of 2 075 188 animals was used for these analysis. The estimates of three effective numbers (N_ef, N_ge, and N_enf) decreased from 418.6 to 50.3, 86.6 to 7.3, and 109.2 to 8.5, respectively, during the period 1960–2000. The increasing differences between two kinds of genetic diversity indices derived from N_ge and N_ef showed that large part of the reduced diversity from 1980 was attributed to genetic drift caused by the intensive use of particular limited number of sires. In gene dropping analysis, probabilities of extinction of founder alleles were derived from their distributions of frequency in the current animals. Several founders showed low probabilities of allele extinction, irrespective of their relatively low genetic contributions. This suggests that these founders have lineages through which their alleles are surely transmitted to the current breed. The use of these founders as a strategy for recovering the genetic diversity was discussed.     

Estimation of dominance genetic variances for reproductive traits and growth traits of calves in Japanese Black cattle

The dominance genetic effects for reproductive and calf growth abilities in the practical Japanese Black populations were examined using average information (AI) algorithm restricted maximum likelihood (REML) under animal models. The reproductive traits were observed in Japanese Black cattle maintained at Tottori and Okinawa prefectures, and growth traits of calves were observed in cattle at Okinawa. The average of dominance relationships in Tottori ranged from 0.2 to 0.4%, while the level in Okinawa was lower and sparse compared with Tottori. The proportions of the dominance variances to sum of additive and dominance variances () were all 0.02 for reproductive traits in Tottori. In contrast, the proportion was 0.02–0.64 in Okinawa regardless of the level of dominance relationships. These proportions suggested that the dominance might affect the expression of calving interval, days open and gestation length in Okinawa, where breeding units were spread over many islands. Although the dominance variances could not estimate birthweight, w as 0.34 for calf market weight and 0.27 for average daily gain from birth to calf market in Okinawa. These values also suggested that the dominance might affect the early growth of calves. In the near future, genetic relationships will become more complicated with continuation of the current selection and mating systems. Therefore, genetic evaluation accounting for dominance effects would be necessary for particular traits and populations.     

Application of highly differentiated SNPs between Japanese Black and Holstein to a breed assignment test between Japanese Black and F1 (Japanese Black x Holstein) and Holstein

Two taurine breeds, Japanese Black and Holstein, established from geographically distant origins and selected for different uses, beef and dairy, were extensively genotyped using a genome‐wide single nucleotide polymorphism (SNP) chip with more than 1000 animals of each breed. The genetic structure was examined by principal component analysis, in which the first principal component clearly separated the two breeds and explained more than 15% of the variance. Highly differentiated SNPs were detected throughout the genome, some of which were clustered within small regions on BTA4 (79.2–79.7 Mb, Btau4.0) and BTA26 (22.2–23.6 Mb). A breed assignment test was developed using 18 highly differentiated SNPs to distinguish Japanese Black from F₁ (Japanese Black × Holstein) and Holstein. The error rate that an F₁ or Holstein animal is misjudged as Japanese Black was expected to be < 0.8%, while the error rate that a Japanese Black animal is misjudged as F₁ or Holstein was expected to be < 0.001%. This test provides a reliable and powerful method to detect breed label falsification in retail beef.     

Genome-wide detection of signatures of selection in three Valdostana cattle populations

The Valdostana is a local dual purpose cattle breed developed in Italy. Three populations are recognized within this breed, based on coat colour, production level, morphology and temperament: Valdostana Red Pied (VPR), Valdostana Black Pied (VPN) and Valdostana Chestnut (VCA). Here, we investigated putative genomic regions under selection among these three populations using the Bovine 50K SNP array by combining three different statistical methods based either on allele frequencies (F_ST) or extended haplotype homozygosity (iHS and Rsb). In total, 8, 5 and 8 chromosomes harbouring 13, 13 and 16 genomic regions potentially under selection were identified by at least two approaches in VPR, VPN and VCA, respectively. Most of these candidate regions were population-specific but we found one common genomic region spanning 2.38 Mb on BTA06 which either overlaps or is located close to runs of homozygosity islands detected in the three populations. This region included inter alia two well-known genes: KDR, a well-established coat colour gene, and CLOCK, which plays a central role in positive regulation of inflammatory response and in the regulation of the mammalian circadian rhythm. The other candidate regions identified harboured genes associated mainly with milk and meat traits as well as genes involved in immune response/inflammation or associated with behavioural traits. This last category of genes was mainly identified in VCA, which is selected for fighting ability. Overall, our results provide, for the first time, a glimpse into regions of the genome targeted by selection in Valdostana cattle. Finally, this study illustrates the relevance of using multiple complementary approaches to identify genomic regions putatively under selection in livestock.     

Genetic structure and relationships of 16 Asian and European cattle populations using DigiTag2 assay

In this study, we genotyped 117 autosomal single nucleotide polymorphisms using a DigiTag2 assay to assess the genetic diversity, structure and relationships of 16 Eurasian cattle populations, including nine cattle breeds and seven native cattle. Phylogenetic and principal component analyses showed that Bos taurus and Bos indicus populations were clearly distinguished, whereas Japanese Shorthorn and Japanese Polled clustered with European populations. Furthermore, STRUCTURE analysis demonstrated the distinct separation between Bos taurus and Bos indicus (K=2), and between European and Asian populations (K=3). In addition, Japanese Holstein exhibited an admixture pattern with Asian and European cattle (K=3‐5). Mongolian (K=13‐16) and Japanese Black (K=14‐16) populations exhibited admixture patterns with different ancestries. Bos indicus populations exhibited a uniform genetic structure at K=2‐11, thereby suggesting that there are close genetic relationships among Bos indicus populations. However, the Bhutan and Bangladesh populations formed a cluster distinct from the other Bos indicus populations at K=12‐16. In conclusion, our study could sufficiently explain the genetic construction of Asian cattle populations, including: (i) the close genetic relationships among Bos indicus populations; (ii) the genetic influences of European breeds on Japanese breeds; (iii) the genetic admixture in Japanese Holstein, Mongolian and Japanese Black cattle; and (iv) the genetic subpopulations in Southeast Asia.