The mtDNA haplogroup P of modern Asian cattle: A genetic legacy of Asian aurochs?

BACKGROUND: Aurochs (Bos primigenius) were distributed throughout large parts of Eurasia and Northern Africa during the late Pleistocene and the early Holocene, and all modern cattle are derived from the aurochs. Although the mtDNA haplogroups of most modern cattle belong to haplogroups T and I, several additional haplogroups (P, Q, R, C and E) have been identified in modern cattle and aurochs. Haplogroup P was the most common haplogroup in European aurochs, but so far, it has been identified in only three of >3,000 submitted haplotypes of modern Asian cattle. METHODOLOGY: We sequenced the complete mtDNA D-loop region of 181 Japanese Shorthorn cattle and analyzed these together with representative bovine mtDNA sequences. The haplotype P of Japanese Shorthorn cattle was analyzed along with that of 36 previously published European aurochs and three modern Asian cattle sequences using the hypervariable 410 bp of the D-loop region. CONCLUSIONS: We detected the mtDNA haplogroup P in Japanese Shorthorn cattle with an extremely high frequency (83/181). Phylogenetic networks revealed two main clusters, designated as Pa for haplogroup P in European aurochs and Pc in modern Asian cattle. We also report the genetic diversity of haplogroup P compared with the sequences of extinct aurochs. No shared haplotypes are observed between the European aurochs and the modern Asian cattle. This finding suggests the possibility of local and secondary introgression events of haplogroup P in northeast Asian cattle, and will contribute to a better understanding of its origin and genetic diversity.

Low mitochondrial DNA diversity of Japanese Polled and Kuchinoshima feral cattle.

This study aims to estimate the mitochondrial genetic diversity and structure of Japanese Polled and Kuchinoshima feral cattle, which are maintained in small populations. We determined the mitochondrial DMA (mtDNA) displacement loop (D-loop) sequences for both cattle populations and analyzed these in conjunction with previously published data from Northeast Asian cattle populations. Our findings showed that Japanese native cattle have a predominant, Asian-specific mtDNA haplogroup T4 with high frequencies (0.43-0.81). This excluded Kuchinoshima cattle (32 animals), which had only one mtDNA haplotype belonging to the haplogroup T3. Japanese Polled showed relatively lower mtDNA diversity in the average sequence divergence (0.0020) than other Wagyu breeds (0.0036-0.0047). Japanese Polled have been maintained in a limited area of Yamaguchi, and the population size is now less than 200. Therefore, low mtDNA diversity in the Japanese Polled could be explained by the decreasing population size in the last three decades. We found low mtDNA diversity in both Japanese Polled and Kuchinoshima cattle. The genetic information obtained in this study will be useful for maintaining these populations and for understanding the origin of Japanese native cattle.

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.