The multifaceted origin of taurine cattle reflected by the mitochondrial genome

Evolution and legacy of East Asian aurochs

Aurochs (Bos primigenius), once widely distributed in Afro-Eurasia, became extinct in the early 1600 s. However, their phylogeography and relative contributions to domestic cattle remain unknown. In this study, we analyzed 16 genomes of ancient aurochs and three mitogenomes of ancient bison (Bison priscus) excavated in East Asia, dating from 43,000 to 3,590 years ago. These newly generated data with previously published genomic information on aurochs as well as ancient/extant domestic cattle worldwide through genome analysis. Our findings revealed significant genetic divergence between East Asian aurochs and their European, Near Eastern, and African counterparts on the basis of both mitochondrial and nuclear genomic data. Furthermore, we identified evidence of gene flow from East Asian aurochs into ancient and present-day taurine cattle, suggesting their potential role in facilitating the environmental adaptation of domestic cattle.

The evolution of contemporary livestock species: Insights from mitochondrial genome

The domestication of animals marks a pivotal moment in human history, profoundly influencing our demographic and cultural progress. This process has led to significant genetic, behavioral, and physical changes in livestock species compared to their wild ancestors. Understanding the evolutionary history and genetic diversity of livestock species is crucial, and mitochondrial DNA (mtDNA) has emerged as a robust marker for investigating molecular diversity in animals. Its highly conserved gene content across animal species, minimal duplications, absence of introns, and short intergenic regions make mtDNA analysis ideal for such studies. Mitochondrial DNA analysis has uncovered distinct cattle domestication events dating back to 8000 years BC in Southwestern Asia. The sequencing of water buffalo mtDNA in 2004 provided important insights into their domestication history. Caprine mtDNA analysis identified three haplogroups, indicating varied maternal origins. Sheep, domesticated 12,000 years ago, exhibit diverse mtDNA lineages, suggesting multiple domestication events. Ovine mtDNA studies revealed clades A, B, C, and a fourth lineage, group D. The origins of domestic pigs were traced to separate European and Asian events followed by interbreeding. In camels, mtDNA elucidated the phylogeographic structure and genetic differentiation between wild and domesticated species. Horses, domesticated around 3500 BC, show significant mtDNA variability, highlighting their diverse origins. Yaks exhibit unique adaptations for high-altitude environments, with mtDNA analysis providing insights into their adaptation. Chicken mtDNA studies supported a monophyletic origin from Southeast Asia's red jungle fowl, with evidence of multiple origins. This review explores livestock evolution and diversity through mtDNA studies, focusing on cattle, water buffalo, goat, sheep, pig, camel, horse, yak and chicken. It highlights mtDNA's significance in unraveling maternal lineages, genetic diversity, and domestication histories, concluding with insights into its potential application in improving livestock production and reproduction dynamics.

The genomic natural history of the aurochs

Now extinct, the aurochs (Bos primigenius) was a keystone species in prehistoric Eurasian and North African ecosystems, and the progenitor of cattle (Bos taurus), domesticates that have provided people with food and labour for millennia1. Here we analysed 38 ancient genomes and found 4 distinct population ancestries in the aurochs-European, Southwest Asian, North Asian and South Asian-each of which has dynamic trajectories that have responded to changes in climate and human influence. Similarly to Homo heidelbergensis, aurochsen first entered Europe around 650 thousand years ago2, but early populations left only trace ancestry, with both North Asian and European B. primigenius genomes coalescing during the most recent glaciation. North Asian and European populations then appear separated until mixing after the climate amelioration of the early Holocene. European aurochsen endured the more severe bottleneck during the Last Glacial Maximum, retreating to southern refugia before recolonizing from Iberia. Domestication involved the capture of a small number of individuals from the Southwest Asian aurochs population, followed by early and pervasive male-mediated admixture involving each ancestral strain of aurochs after domestic stocks dispersed beyond their cradle of origin.

Tracing early pastoralism in Central Europe using sedimentary ancient DNA

Central European forests have been shaped by complex human interactions throughout the Holocene, with significant changes following the introduction of domesticated animals in the Neolithic (∼7.5-6.0 ka before present [BP]). However, understanding early pastoral practices and their impact on forests is limited by methods for detecting animal movement across past landscapes. Here, we examine ancient sedimentary DNA (sedaDNA) preserved at the Velký Mamuťák rock shelter in northern Bohemia (Czech Republic), which has been a forested enclave since the early Holocene. We find that domesticated animals, their associated microbiomes, and plants potentially gathered for fodder have clear representation by the Late Neolithic, around 6.0 ka BP, and persist throughout the Bronze Age into recent times. We identify a change in dominant grazing species from sheep to pigs in the Bronze Age (∼4.1-3.0 ka BP) and interpret the impact this had in the mid-Holocene retrogressions that still define the structure of Central European forests today. This study highlights the ability of ancient metagenomics to bridge archaeological and paleoecological methods and provide an enhanced perspective on the roots of the "Anthropocene."

Late History of Cattle Breeds in Central Europe in Light of Genetic and Archaeogenetic Sources-Overview, Thoughts, and Perspectives

Although Europe was not a primary centre of cattle domestication, its expansion from the Middle East and subsequent development created a complex pattern of cattle breed diversity. Many isolated populations of local historical breeds still carry the message about the physical and genetic traits of ancient populations. Since the way of life of human communities starting from the eleventh millennium BP was strongly determined by livestock husbandry, the knowledge of cattle diversity through the ages is helpful in the interpretation of many archaeological findings. Historical cattle diversity is currently at the intersection of two leading directions of genetic research. Firstly, it is archaeogenetics attempting to recover and interpret the preserved genetic information directly from archaeological finds. The advanced archaeogenetic approaches meet with the population genomics of extant cattle populations. The immense amount of genetic information collected from living cattle, due to its key economic role, allows for reconstructing the genetic profiles of the ancient populations backwards. The present paper aims to place selected archaeogenetic, genetic, and genomic findings in the picture of cattle history in Central Europe, as suggested by archaeozoological and historical records. Perspectives of the methodical connection between the genetic approaches and the approaches of traditional archaeozoology, such as osteomorphology and osteometry, are discussed. The importance, actuality, and effectiveness of combining different approaches to each archaeological find, such as morphological characterization, interpretation of the historical context, and molecular data, are stressed.

Comparison between indicine and taurine cattle DNA methylation reveals epigenetic variation associated to differences in morphological adaptive traits

Indicine and taurine subspecies present distinct morphological traits as a consequence of environmental adaptation and artificial selection. Although the two subspecies have been characterized and compared at genome-wide level and at specific loci, their epigenetic diversity has not yet been explored. In this work, Reduced Representation Bisulphite Sequencing (RRBS) profiling of the taurine Angus (A) and indicine Nellore (N) cattle breeds was applied to identify methylation differences between the two subspecies. Genotyping by sequencing (GBS) of the same animals was performed to detect single nucleotide polymorphisms (SNPs) at cytosines in CpG dinucleotides and remove them from the differential methylation analysis. A total of 660,845 methylated cytosines were identified within the CpG context (CpGs) across the 10 animals sequenced (5 N and 5 A). A total of 25,765 of these were differentially methylated (DMCs). Most DMCs clustered in CpG stretches nearby genes involved in cellular and anatomical structure morphogenesis. Also, sequences flanking DMC were enriched in SNPs compared to all other CpGs, either methylated or unmethylated in the two subspecies. Our data suggest a contribution of epigenetics to the regulation and divergence of anatomical morphogenesis in the two subspecies relevant for cattle evolution and sub-species differentiation and adaptation.

Ancient Mongolian aurochs genomes reveal sustained introgression and management in East Asia

Societies in East Asia have utilized domesticated cattle for over 5000 years, but the genetic history of cattle in East Asia remains understudied. Genome-wide analyses of 23 ancient Mongolian cattle reveal that East Asian aurochs and ancient East Asian taurine cattle are closely related, but neither are closely related to any modern East Asian breeds. We observe binary variation in aurochs diet throughout the early Neolithic, and genomic evidence shows millennia of sustained male-dominated introgression. We identify a unique connection between ancient Mongolian aurochs and the European Hereford breed. These results point to the likelihood of human management of aurochs in Northeast Asia prior to and during the initial adoption of taurine cattle pastoralism.

Ancient DNA confirms diverse origins of early post-Columbian cattle in the Americas

Before the arrival of Europeans, domestic cattle (Bos taurus) did not exist in the Americas, and most of our knowledge about how domestic bovines first arrived in the Western Hemisphere is based on historical documents. Sixteenth-century colonial accounts suggest that the first cattle were brought in small numbers from the southern Iberian Peninsula via the Canary archipelago to the Caribbean islands where they were bred locally and imported to other circum-Caribbean regions. Modern American heritage cattle genetics and limited ancient mtDNA data from archaeological colonial cattle suggest a more complex story of mixed ancestries from Europe and Africa. So far little information exists to understand the nature and timing of the arrival of these mixed-ancestry populations. In this study we combine ancient mitochondrial and nuclear DNA from a robust sample of some of the earliest archaeological specimens from Caribbean and Mesoamerican sites to clarify the origins and the dynamics of bovine introduction into the Americas. Our analyses support first arrival of cattle from diverse locales and potentially confirm the early arrival of African-sourced cattle in the Americas, followed by waves of later introductions from various sources over several centuries.

The Complete Mitochondrial Genome and Phylogenetic Analyses of To Chicken in Vietnam

Indigenous chicken breeds have both cultural significance and economic value since they possess unique genetic characteristics that enable them to adapt to the local environment and contribute to biodiversity, food security, and sustainable agriculture in Vietnam. To (Tò in Vietnamese) chicken, a Vietnamese indigenous chicken breed, is popularly raised in Thai Binh province; however, little known is about the genetic diversity of this breed. In this study, we sequenced the complete mitochondrial genome of To chicken for a better understanding of the diversity and origin of the breed. The results of sequencing showed that the mitochondrial genome of To chicken spans a total length of 16,784 base pairs and comprises one non-coding control region (known as the displacement-loop (D-loop) region), two ribosomal RNA genes, 13 protein-coding genes, and 22 transfer RNA genes. The phylogenetic tree analyses and estimated genetic distances based on 31 complete mitochondrial genome sequences indicated that To chicken has a close genetic distance with the Laotian native chicken breed, Lv'erwu breed in China, and Nicobari black and Kadaknath breeds in India. The result of the current study might be important for conservation, breeding, and further genetic studies of To chicken.

Mitochondrial Inheritance Following Nuclear Transfer: From Cloned Animals to Patients with Mitochondrial Disease

Mitochondria are indispensable power plants of eukaryotic cells that also act as a major biochemical hub. As such, mitochondrial dysfunction, which can originate from mutations in the mitochondrial genome (mtDNA), may impair organism fitness and lead to severe diseases in humans. MtDNA is a multi-copy, highly polymorphic genome that is uniparentally transmitted through the maternal line. Several mechanisms act in the germline to counteract heteroplasmy (i.e., coexistence of two or more mtDNA variants) and prevent expansion of mtDNA mutations. However, reproductive biotechnologies such as cloning by nuclear transfer can disrupt mtDNA inheritance, resulting in new genetic combinations that may be unstable and have physiological consequences. Here, we review the current understanding of mitochondrial inheritance, with emphasis on its pattern in animals and human embryos generated by nuclear transfer.

Genetic and phenotypic associations of mitochondrial DNA copy number, SNP, and haplogroups with growth and carcass traits in beef cattle

Mitochondrial DNA copy number (mtDNA CN) is heritable and easily obtained from low-pass sequencing (LPS). This study investigated the genetic correlation of mtDNA CN with growth and carcass traits in a multi-breed and crossbred beef cattle population. Blood, leucocyte, and semen samples were obtained from 2,371 animals and subjected to LPS that resulted in nuclear DNA (nuDNA) and mtDNA sequence reads. Mitochondrial DNA CN was estimated as the ratio of mtDNA to nuDNA coverages. Variant calling was performed from mtDNA, and 11 single nucleotide polymorphisms (SNP) were identified in the population. Samples were classified in taurine haplogroups. Haplogroup and mtDNA type were further classified based on the 11 segregating SNP. Growth and carcass traits were available for between 7,249 and 60,989 individuals. Associations of mtDNA CN, mtDNA haplogroups, mtDNA types, and mtDNA SNP with growth and carcass traits were estimated with univariate animal models, and genetic correlations were estimated with a bivariate animal model based on pedigree. Mitochondrial DNA CN tended (P-value ≤0.08) to be associated with birth weight and weaning weight. There was no association (P-value >0.10) between mtDNA SNP, haplogroups, or types with growth and carcass traits. Genetic correlation estimates of mtDNA CN were -0.30 ± 0.16 with birth weight, -0.31 ± 0.16 with weaning weight, -0.15 ± 0.14 with post-weaning gain, -0.11 ± 0.19 with average daily dry-matter intake, -0.04 ± 0.22 with average daily gain, -0.29 ± 0.13 with mature cow weight, -0.11 ± 0.13 with slaughter weight, -0.14 ± 0.13 with carcass weight, -0.07 ± 0.14 with carcass backfat, 0.14 ± 0.14 with carcass marbling, and -0.06 ± 0.14 with ribeye area. In conclusion, mtDNA CN was negatively correlated with most traits investigated, and the genetic correlation was stronger with growth traits than with carcass traits.

Legacies of domestication, Neolithic diffusion and trade between Indian subcontinent and Island Southeast Asia shape maternal genetic diversity of Andaman cattle

Andaman cattle is a precious indigenous livestock species endemic to Andaman and Nicobar Islands, India. Till date, origin and genetic makeup of the breed which is warranted for breed conservation is not known. Moreover, the spread of zebu cattle from Indus valley to different parts of Island Southeast Asia (ISEA) is not properly understood. Here, we report the genetic diversity, population structure of Andaman cattle and their evolution in the context of epicentre of zebu domestication and ISEA. High genetic diversity in complete mitochondrial D-loop sequences indicated the ability of the breed to withstand impending climate change. Total 81 haplotypes were detected and all of them except three belonged to Bos indicus. The presence of taurine haplotypes in Andaman cattle indicate introgression by European-derived cattle. A poor phylogenetic signal of Andaman cattle with genetic affinities with cattle of Indian subcontinent and ISEA was observed. The poor phylogenetic structure may be due to multidirectional gene flow from Indian subcontinent and ISEA, with which Andaman shares a close cultural and trade relationship from Neolithic age. We hypothesize that Andaman cattle is the outcome of Neolithic diffusion from centre of zebu domestication along with multidirectional commercial exchange between Indian subcontinent and ISEA.

Analysis of the earliest complete mtDNA genome of a Caribbean colonial horse (Equus caballus) from 16th-century Haiti

Unlike other European domesticates introduced in the Americas after the European invasion, equids (Equidae) were previously in the Western Hemisphere but were extinct by the late Holocene era. The return of equids to the Americas through the introduction of the domestic horse (Equus caballus) is documented in the historical literature but is not explored fully either archaeologically or genetically. Historical documents suggest that the first domestic horses were brought from the Iberian Peninsula to the Caribbean in the late 15th century CE, but archaeological remains of these early introductions are rare. This paper presents the mitochondrial genome of a late 16th century horse from the Spanish colonial site of Puerto Real (northern Haiti). It represents the earliest complete mitogenome of a post-Columbian domestic horse in the Western Hemisphere offering a unique opportunity to clarify the phylogeographic history of this species in the Americas. Our data supports the hypothesis of an Iberian origin for this early translocated individual and clarifies its phylogenetic relationship with modern breeds in the Americas.

Influence of environmental factors and genetic variation on mitochondrial DNA copy number

Mitochondrial DNA copy number (mtDNA CN) has been shown to be highly heritable and associated with traits of interest in humans. However, studies are lacking in the literature for livestock species such as beef cattle. In this study, 2,371 individuals from a crossbred beef population comprising the Germplasm Evaluation program from the U.S. Meat Animal Research Center had samples of blood, leucocyte, or semen collected for low-pass sequencing (LPS) that resulted in both nuclear DNA (nuDNA) and mitochondrial DNA (mtDNA) sequence reads. Mitochondrial DNA CN was estimated based on the ratio of mtDNA to nuDNA coverages. Genetic parameters for mtDNA CN were estimated from an animal model based on a genomic relationship matrix (~87K SNP from the nuDNA). Different models were used to test the effects of tissue, sex, age at sample collection, heterosis, and breed composition. Maternal effects, assessed by fitting a maternal additive component and by fitting eleven SNP on the mtDNA, were also obtained. As previously reported, mtDNA haplotypes were used to classify individuals into Taurine haplogroups (T1, T2, T3/T4, and T5). Estimates of heritability when fitting fixed effects in addition to the intercept were moderate, ranging from 0.11 to 0.31 depending on the model. From a model ignoring contemporary group, semen samples had the lowest mtDNA CN, as expected, followed by blood and leucocyte samples (P ≤ 0.001). The effect of sex and the linear and quadratic effects of age were significant (P ≤ 0.02) depending on the model. When significant, females had greater mtDNA CN than males. The effects of heterosis and maternal heterosis were not significant (P ≥ 0.47). The estimates of maternal and mtDNA heritability were near zero (≤0.03). Most of the samples (98%) were classified as haplogroup T3. Variation was observed in the mtDNA within Taurine haplogroups, which enabled the identification of 24 haplotypes. These results suggest that mtDNA CN is under nuclear genetic control and would respond favorably to selection.

Mitochondrial DNA diversity of the Sardinian local cattle stock

The aim of this research was to characterize the genetic diversity of the Sarda (Sa, n = 131), Sardo Bruna (SB, n = 44) and Sardo Modicana (SM, n = 26) cattle breeds, reared in the island of Sardinia (Italy). A portion of the mitochondrial DNA hypervariable region was sequenced, in order to identify a potential signature of African introgression. The F_ST coefficients among populations ranged between 0.056 for Sa vs SB and 0.167 for SB vs SM. AMOVA analysis indicated there was a significant differentiation of the three breeds, although most of diversity was gathered at the within-breed level. The Median Joining Network of the Sardinian sequences showed a potential founder effect signature. A MJ network including Sardinian cattle plus African, Italian, Iberian and Asian sequences, revealed the presence of haplogroup T3, already detected in Sa cattle, and the presence of Hg T1 and Hg T1′2′3, in Sa and SB. The presence of a private haplotype belonging to haplogroup T1, which is characteristic of African taurine breeds, may be due to the introgression of Sardinian breeds with African cattle, either directly (most probable source: North African cattle) or indirectly (through a Mediterranean intermediary already introgressed with African blood).

Evaluation of Maternal Genetic Background of Two Hungarian Autochthonous Sheep Breeds Coming from Different Geographical Directions

The aim of our research was the evaluation of the maternal genetic background of two Hungarian autochthonous sheep breeds of different geographical origin. A major argument for the preservation of endangered animal breeds is their documented past and historical importance. These also include the registration of pedigree data. This is the first study to evaluate and compare Tsigai and Cikta sheep in Hungary. Our investigation is based on two complete sequences of mitochondrial DNA (cytochrome b gene and control region). Our research was performed on these two sheep breeds with markedly different breed histories and breed characteristics to determine a possible common maternal genetic background, as ultimately the origin of both breeds can be traced back to Asia Minor. Between 2015 and 2017, a total of 203 biological samples were taken using a newly introduced founder sampling method. We found that the prevailing haplogroup B accounted for over 80% of both breeds, strengthening the common ancestral root. However, the pairwise genetic differentiation estimates (KST) calculated using the sequence-based statistics for cytochrome b gene and control region were 0.034 and 0.021, respectively (both at level p < 0.05); thus, revealing genetic differentiation in both sequences between the Tsigai and Cikta. We note that the known different history of the breeds is clearly justified by the currently studied deviations in their maternal genetic background.

Mitonuclear incompatibility as a hidden driver behind the genome ancestry of African admixed cattle

Background

Africa is an important watershed in the genetic history of domestic cattle, as two lineages of modern cattle, Bos taurus and B. indicus, form distinct admixed cattle populations. Despite the predominant B. indicus nuclear ancestry of African admixed cattle, B. indicus mitochondria have not been found on the continent. This discrepancy between the mitochondrial and nuclear genomes has been previously hypothesized to be driven by male-biased introgression of Asian B. indicus into ancestral African B. taurus. Given that this hypothesis mandates extreme demographic assumptions relying on random genetic drift, we propose a novel hypothesis of selection induced by mitonuclear incompatibility and assess these hypotheses with regard to the current genomic status of African admixed cattle.

Results

By analyzing 494 mitochondrial and 235 nuclear genome sequences, we first confirmed the genotype discrepancy between mitochondrial and nuclear genome in African admixed cattle: the absence of B. indicus mitochondria and the predominant B. indicus autosomal ancestry. We applied approximate Bayesian computation (ABC) to assess the posterior probabilities of two selection hypotheses given this observation. The results of ABC indicated that the model assuming both male-biased B. indicus introgression and selection induced by mitonuclear incompatibility explains the current genomic discrepancy most accurately. Subsequently, we identified selection signatures at autosomal loci interacting with mitochondria that are responsible for integrity of the cellular respiration system. By contrast with B. indicus-enriched genome ancestry of African admixed cattle, local ancestries at these selection signatures were enriched with B. taurus alleles, concurring with the key expectation of selection induced by mitonuclear incompatibility.

Conclusions

Our findings support the current genome status of African admixed cattle as a potential outcome of male-biased B. indicus introgression, where mitonuclear incompatibility exerted selection pressure against B. indicus mitochondria. This study provides a novel perspective on African cattle demography and supports the role of mitonuclear incompatibility in the hybridization of mammalian species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01206-x.

[Signatures of selection and candidate genes for adaptation to extreme environmental factors in the genomes of Turano-Mongolian cattle breeds]

Changes in the environment force populations of organisms to adapt to new conditions, either through phenotypic plasticity or through genetic or epigenetic changes. Signatures of selection, such as specific changes in the frequency of alleles and haplotypes, as well as the reduction or increase in genetic diversity, help to identify changes in the cattle genome in response to natural and artificial selection, as well as loci and genetic variants directly affecting adaptive and economically important traits. Advances in genetics and biotechnology enable a rapid transfer of unique genetic variants that have originated in local cattle breeds in the process of adaptation to local environments into the genomes of cosmopolitan high-performance breeds, in order to preserve their outstanding performance in new environments. It is also possible to use genomic selection approach to increase the frequency of already present adaptive alleles in cosmopolitan breeds. The review examines recent work on the origin and evolution of Turano-Mongolian cattle breeds, adaptation of Turano-Mongolian cattle to extreme environments, and summarizes available information on potential candidate genes for climate adaptation of Turano-Mongolian breeds, including cold resistance genes, immune response genes, and high-altitude adaptation genes. The authors conclude that the current literature data do not provide preference to one of the two possible scenarios of Turano-Mongolian breed origins: as a result of the domestication of a wild aurochs at East Asia or as a result of the migration of taurine proto-population from the Middle East. Turano-Mongolian breeds show a high degree of adaptation to extreme climatic conditions (cold, heat, lack of oxygen in the highlands) and parasites (mosquitoes, ticks, bacterial and viral infections). As a result of high-density genotyping and sequencing of genomes and transcriptomes, prospective candidate genes and genetic variants involved in adaptation to environmental factors have recently been identified.

Kouprey (Bos sauveli) genomes unveil polytomic origin of wild Asian Bos

The evolution of the genera Bos and Bison, and the nature of gene flow between wild and domestic species, is poorly understood, with genomic data of wild species being limited. We generated two genomes from the likely extinct kouprey (Bos sauveli) and analyzed them alongside other Bos and Bison genomes. We found that B. sauveli possessed genomic signatures characteristic of an independent species closely related to Bos javanicus and Bos gaurus. We found evidence for extensive incomplete lineage sorting across the three species, consistent with a polytomic diversification of the major ancestry in the group, potentially followed by secondary gene flow. Finally, we detected significant gene flow from an unsampled Asian Bos-like source into East Asian zebu cattle, demonstrating both that the full genomic diversity and evolutionary history of the Bos complex has yet to be elucidated and that museum specimens and ancient DNA are valuable resources to do so.

•

We generated two genomes from the likely extinct kouprey (Bos sauveli)

•

Extensive mt and nuclear-genome-wide incomplete lineage sorting across wild Asian Bos

•

Initial polytomic diversification of the wild Asian Bos—kouprey, banteng, and gaur

Deleterious protein-coding variants in diverse cattle breeds of the world

The domestication of wild animals has resulted in a reduction in effective population sizes, which can affect the deleterious mutation load of domesticated breeds. In addition, artificial selection contributes to the accumulation of deleterious mutations because of an increased rate of inbreeding among domesticated animals. Since founder population sizes and artificial selection differ between cattle breeds, their deleterious mutation load can vary. We investigated this question by using whole-genome data from 432 animals belonging to 54 worldwide cattle breeds. Our analysis revealed a negative correlation between genomic heterozygosity and nonsynonymous-to-silent diversity ratio, which suggests a higher proportion of single nucleotide variants (SNVs) affecting proteins in low-diversity breeds. Our results also showed that low-diversity breeds had a larger number of high-frequency (derived allele frequency (DAF) > 0.51) deleterious SNVs than high-diversity breeds. An opposite trend was observed for the low-frequency (DAF ≤ 0.51) deleterious SNVs. Overall, the number of high-frequency deleterious SNVs was larger in the genomes of taurine cattle breeds than of indicine breeds, whereas the number of low-frequency deleterious SNVs was larger in the genomes of indicine cattle than in those of taurine cattle. Furthermore, we observed significant variation in the counts of deleterious SNVs within taurine breeds. The variations in deleterious mutation load between taurine and indicine breeds could be attributed to the population sizes of the wild progenitors before domestication, whereas the variations observed within taurine breeds could be due to differences in inbreeding level, strength of artificial selection, and/or founding population size. Our findings imply that the incidence of genetic diseases can vary between cattle breeds.

Population differentiated copy number variation of Bos taurus, Bos indicus and their African hybrids

Background

CNV comprises a large proportion in cattle genome and is associated with various traits. However, there were few population-scale comparison studies on cattle CNV.

Results

Here, autosome-wide CNVs were called by read depth of NGS alignment result and copy number variation regions (CNVRs) defined from 102 Eurasian taurine (EAT) of 14 breeds, 28 Asian indicine (ASI) of 6 breeds, 22 African taurine (AFT) of 2 breeds, and 184 African humped cattle (AFH) of 17 breeds. The copy number of every CNVRs were compared between populations and CNVRs with population differentiated copy numbers were sorted out using the pairwise statistics V_ST and Kruskal-Wallis test. Three hundred sixty-two of CNVRs were significantly differentiated in both statistics and 313 genes were located on the population differentiated CNVRs.

Conclusion

For some of these genes, the averages of copy numbers were also different between populations and these may be candidate genes under selection. These include olfactory receptors, pathogen-resistance, parasite-resistance, heat tolerance and productivity related genes. Furthermore, breed- and individual-level comparison was performed using the presence or copy number of the autosomal CNVRs. Our findings were based on identification of CNVs from short Illumina reads of 336 individuals and 39 breeds, which to our knowledge is the largest dataset for this type of analysis and revealed important CNVs that may play a role in cattle adaption to various environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07808-7.

Harmful mutation load in the mitochondrial genomes of cattle breeds

Objective

Domestication of wild animals results in a reduction in the effective population size, and this could affect the deleterious mutation load of domesticated breeds. Furthermore, artificial selection will also contribute to the accumulation of deleterious mutations due to the increased rate of inbreeding among these animals. The process of domestication, founder population size, and artificial selection differ between cattle breeds, which could lead to a variation in their deleterious mutation loads. We investigated this using mitochondrial genome data from 364 animals belonging to 18 cattle breeds of the world.

Results

Our analysis revealed more than a fivefold difference in the deleterious mutation load among cattle breeds. We also observed a negative correlation between the breed age and the proportion of deleterious amino acid-changing polymorphisms. This suggests a proportionally higher deleterious SNPs in young breeds compared to older breeds. Our results highlight the magnitude of difference in the deleterious mutations present in the mitochondrial genomes of various breeds. The results of this study could be useful in predicting the rate of incidence of genetic diseases in different breeds.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-021-05664-y.

On the origin and diversification of Podolian cattle breeds: testing scenarios of European colonization using genome-wide SNP data

Background

During the Neolithic expansion, cattle accompanied humans and spread from their domestication centres to colonize the ancient world. In addition, European cattle occasionally intermingled with both indicine cattle and local aurochs resulting in an exclusive pattern of genetic diversity. Among the most ancient European cattle are breeds that belong to the so-called Podolian trunk, the history of which is still not well established. Here, we used genome-wide single nucleotide polymorphism (SNP) data on 806 individuals belonging to 36 breeds to reconstruct the origin and diversification of Podolian cattle and to provide a reliable scenario of the European colonization, through an approximate Bayesian computation random forest (ABC-RF) approach.

Results

Our results indicate that European Podolian cattle display higher values of genetic diversity indices than both African taurine and Asian indicine breeds. Clustering analyses show that Podolian breeds share close genomic relationships, which suggests a likely common genetic ancestry. Among the simulated and tested scenarios of the colonization of Europe from taurine cattle, the greatest support was obtained for the model assuming at least two waves of diffusion. Time estimates are in line with an early migration from the domestication centre of non-Podolian taurine breeds followed by a secondary migration of Podolian breeds. The best fitting model also suggests that the Italian Podolian breeds are the result of admixture between different genomic pools.

Conclusions

This comprehensive dataset that includes most of the autochthonous cattle breeds belonging to the so-called Podolian trunk allowed us not only to shed light onto the origin and diversification of this group of cattle, but also to gain new insights into the diffusion of European cattle. The most well-supported scenario of colonization points to two main waves of migrations: with one that occurred alongside with the Neolithic human expansion and gave rise to the non-Podolian taurine breeds, and a more recent one that favoured the diffusion of European Podolian. In this process, we highlight the importance of both the Mediterranean and Danube routes in promoting European cattle colonization. Moreover, we identified admixture as a driver of diversification in Italy, which could represent a melting pot for Podolian cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00639-w.

Mitochondrial genomes from modern and ancient Turano-Mongolian cattle reveal an ancient diversity of taurine maternal lineages in East Asia

Turano-Mongolian cattle are a group of taurine cattle from Northern and Eastern Asia with distinct morphological traits, which are known for their ability to tolerate harsh environments, such as the Asian steppe and the Tibetan plateau. Through the analysis of 170 mitogenomes from ten modern breeds, two sub-lineages within T3 (T3₁₁₉ and T3₀₅₅) were identified as specific of Turano-Mongolian cattle. These two T3 sub-lineages, together with the previously identified T4, were also present in six Neolithic samples, dated to ~3900 years BP, which might represent the earliest domestic taurine stocks from Southwest Asia. The rare haplogroup Q, found in three Tibetan cattle, testifies for the legacy of ancient migrations from Southwest Asia and suggests that the isolated Tibetan Plateau preserved unique prehistoric genetic resources. These findings confirm the geographic substructure of Turano-Mongolian cattle breeds, which have been shaped by ancient migrations and geographic barriers.

Distinguishing African bovids using Zooarchaeology by Mass Spectrometry (ZooMS): New peptide markers and insights into Iron Age economies in Zambia

Assessing past foodways, subsistence strategies, and environments depends on the accurate identification of animals in the archaeological record. The high rates of fragmentation and often poor preservation of animal bones at many archaeological sites across sub-Saharan Africa have rendered archaeofaunal specimens unidentifiable beyond broad categories, such as “large mammal” or “medium bovid”. Identification of archaeofaunal specimens through Zooarchaeology by Mass Spectrometry (ZooMS), or peptide mass fingerprinting of bone collagen, offers an avenue for identification of morphologically ambiguous or unidentifiable bone fragments from such assemblages. However, application of ZooMS analysis has been hindered by a lack of complete reference peptide markers for African taxa, particularly bovids. Here we present the complete set of confirmed ZooMS peptide markers for members of all African bovid tribes. We also identify two novel peptide markers that can be used to further distinguish between bovid groups. We demonstrate that nearly all African bovid subfamilies are distinguishable using ZooMS methods, and some differences exist between tribes or sub-tribes, as is the case for Bovina (cattle) vs. Bubalina (African buffalo) within the subfamily Bovinae. We use ZooMS analysis to identify specimens from extremely fragmented faunal assemblages from six Late Holocene archaeological sites in Zambia. ZooMS-based identifications reveal greater taxonomic richness than analyses based solely on morphology, and these new identifications illuminate Iron Age subsistence economies c. 2200–500 cal BP. While the Iron Age in Zambia is associated with the transition from hunting and foraging to the development of farming and herding, our results demonstrate the continued reliance on wild bovids among Iron Age communities in central and southwestern Zambia Iron Age and herding focused primarily on cattle. We also outline further potential applications of ZooMS in African archaeology.

The genome of the lowland anoa (Bubalus depressicornis) illuminates the origin of river and swamp buffalo

Two types of domestic water buffalo are currently recognized: the river buffalo from the Indian subcontinent and Mediterranean countries and the swamp buffalo from China and Southeast Asia. To test the hypothesis of two separate species of water buffalo, we sequenced the genome of the lowland anoa, Bubalus depressicornis, which is a dwarf wild buffalo endemic to Sulawesi, and two genomes of swamp buffalo, and made comparisons with 12 additional genomes. Three genomic data sets were constructed to infer phylogenetic relationships: the mitochondrial genome (15,468 bp; maternal transmission), two concatenated Y-chromosomal genes, AMELY and DDX3Y (20,036 bp; paternal transmission), and a selection of 30 nuclear genes representing all cattle chromosomes (364,887 bp; biparental transmission). The comparisons between our 30 nuclear gene sequences obtained by read mapping and those directly extracted from Bos taurus and Bubalus bubalis genome assemblies show that the mapping approach revealed higher levels of heterozygosity at both nucleotide sites and indels (insertions and deletions) (0.09-0.15%), as well as several sequence errors (0.07%). Our phylogenetic and molecular dating analyses provide strong evidence that the lowland anoa, river buffalo, and swamp buffalo are three distinct taxa which separated rapidly from each other during the Pleistocene epoch. We therefore conclude that two species of domestic water buffalo should be distinguished: Bubalus bubalis for the river buffalo and Bubalus kerabau for the swamp buffalo. The new classification can have deep implications for understanding the evolution and selection of domesticated forms and for the conservation and management of wild buffalo populations in South and Southeast Asia.

Genetic variation of Nigerian cattle inferred from maternal and paternal genetic markers

The African cattle provide unique genetic resources shaped up by both diverse tropical environmental conditions and human activities, the assessment of their genetic diversity will shade light on the mechanism of their remarkable adaptive capacities. We therefore analyzed the genetic diversity of cattle samples from Nigeria using both maternal and paternal DNA markers. Nigerian cattle can be assigned to 80 haplotypes based on the mitochondrial DNA (mtDNA) D-loop sequences and haplotype diversity was 0.985 + 0.005. The network showed two major matrilineal clustering: the dominant cluster constituting the Nigerian cattle together with other African cattle while the other clustered Eurasian cattle. Paternal analysis indicates only zebu haplogroup in Nigerian cattle with high genetic diversity 1.000 ± 0.016 compared to other cattle. There was no signal of maternal genetic structure in Nigerian cattle population, which may suggest an extensive genetic intermixing within the country. The absence of Bos indicus maternal signal in Nigerian cattle is attributable to vulnerability bottleneck of mtDNA lineages and concordance with the view of male zebu genetic introgression in African cattle. Our study shades light on the current genetic diversity in Nigerian cattle and population history in West Africa.

Characterization and phylogenetic evolution of mitochondrial genome in Tibetan chicken

The objective of this study was to characterize mitochondrial genome and investigate phylogenetic evolution in Tibetan chicken. In this study, four haplotypes were identified based on D-loop sequencing in Tibetan chicken (n = 40), and each representative of four haplotypes was selected for total mitochondrial genome sequencing and analyzed together with published mitochondrial genome data of red jungle fowl. Four haplotypes belonged to three previously published clades, i.e., Clade A, clade B and clade E. Based on D-loop sequencing data, the average haplotype diversity and nucleotide diversity were 0.658 ± 0.065 and 0.00442 ± 0.00094, respectively. The mitochondrial genome of Tibetan chicken is 16,785 bp in size, consisting of 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, 13 protein-coding genes and one non-coding control region (CR). Compared with the mitochondrial genome, a phylogenetic tree based on the D-loop sequence had a messy distribution, and no breed cluster pattern was observed in Tibetan chicken. The results indicate that Tibetan chicken populations in our study have relatively low nucleotide and haplotype diversity and likely share multiple maternal lineages. The D-loop sequence has limited power for the resolution of phylogenetic relationships in comparison with the complete mitochondrial genome.

Cattle mitogenome variation reveals a post-glacial expansion of haplogroup P and an early incorporation into northeast Asian domestic herds

Surveys of mitochondrial DNA (mtDNA) variation have shown that worldwide domestic cattle are characterized by just a few major haplogroups. Two, T and I, are common and characterize Bos taurus and Bos indicus, respectively, while the other three, P, Q and R, are rare and are found only in taurine breeds. Haplogroup P is typical of extinct European aurochs, while intriguingly modern P mtDNAs have only been found in northeast Asian cattle. These Asian P mtDNAs are extremely rare with the exception of the Japanese Shorthorn breed, where they reach a frequency of 45.9%. To shed light on the origin of this haplogroup in northeast Asian cattle, we completely sequenced 14 Japanese Shorthorn mitogenomes belonging to haplogroup P. Phylogenetic and Bayesian analyses revealed: (1) a post-glacial expansion of aurochs carrying haplogroup P from Europe to Asia; (2) that all Asian P mtDNAs belong to a single sub-haplogroup (P1a), so far never detected in either European or Asian aurochs remains, which was incorporated into domestic cattle of continental northeastern Asia possibly ~ 3700 years ago; and (3) that haplogroup P1a mtDNAs found in the Japanese Shorthorn breed probably reached Japan about 650 years ago from Mongolia/Russia, in agreement with historical evidence.

Analyses of mitogenomic markers shed light on the divergence, population dynamics, and demographic history of Pakistani chickens

Pakistan is one of a few sites, associated with the earliest known independent domestication event in the evolutionary history of chicken, which is socio-economically and historically the most important poultry bird in the country. However, the divergence, past population dynamics, and demographic history of Pakistani chickens have not been addressed so far. Therefore, we herein investigated the indigenous Pakistani chickens using mitogenomic markers. We first prepared individual DNA samples from the chicken feathers, and generated nucleotide sequence data, which was then subjected to various population genetics analyses. In molecular phylogenetic analysis, the Pakistani chickens were clustered under nine different clades. Among the wild fowls, the Indian red jungle fowl (IRJF) shared very close affinities to Pakistani chickens. The Bayesian skyline plot showed an increase in the effective population size of Pakistani chickens during the last 50 years. Finally, a time-calibrated phylogeny inferred molecular divergence of the Pakistani chickens. A molecular rate of 3.6 × 10^-6 mutations/site/year (95% HPD interval: 2.28 × 10^-8 to 9.32 × 10^-6) was estimated for the data set. In a rooted tree with root-age of 12058 years (95% HPD interval: 1161-38411), the Pakistani chicken haplotypes showed divergence from IRJF haplotypes around 6987 years (95% HPD interval: 1132-20746) ago, and they shared their most recent common ancestor with Gallus gallus spadiceus, and G. g. jabouillei at the root of the tree. Overall, these results suggest that Pakistani chicken haplotypes share their ancestral gene pool with the IRJF as compared to other red jungle fowl subspecies.

Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species - taurine cattle, zebu and yak - and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.

Illuminating Genetic Mysteries of the Dead Sea Scrolls

The discovery of the 2,000-year-old Dead Sea Scrolls had an incomparable impact on the historical understanding of Judaism and Christianity. "Piecing together" scroll fragments is like solving jigsaw puzzles with an unknown number of missing parts. We used the fact that most scrolls are made from animal skins to "fingerprint" pieces based on DNA sequences. Genetic sorting of the scrolls illuminates their textual relationship and historical significance. Disambiguating the contested relationship between Jeremiah fragments supplies evidence that some scrolls were brought to the Qumran caves from elsewhere; significantly, they demonstrate that divergent versions of Jeremiah circulated in parallel throughout Israel (ancient Judea). Similarly, patterns discovered in non-biblical scrolls, particularly the Songs of the Sabbath Sacrifice, suggest that the Qumran scrolls represent the broader cultural milieu of the period. Finally, genetic analysis divorces debated fragments from the Qumran scrolls. Our study demonstrates that interdisciplinary approaches enrich the scholar's toolkit.

Ancient cattle genomics, origins, and rapid turnover in the Fertile Crescent

Genome-wide analysis of 67 ancient Near Eastern cattle, Bos taurus, remains reveals regional variation that has since been obscured by admixture in modern populations. Comparisons of genomes of early domestic cattle to their aurochs progenitors identify diverse origins with separate introgressions of wild stock. A later region-wide Bronze Age shift indicates rapid and widespread introgression of zebu, Bos indicus, from the Indus Valley. This process was likely stimulated at the onset of the current geological age, ~4.2 thousand years ago, by a widespread multicentury drought. In contrast to genome-wide admixture, mitochondrial DNA stasis supports that this introgression was male-driven, suggesting that selection of arid-adapted zebu bulls enhanced herd survival. This human-mediated migration of zebu-derived genetics has continued through millennia, altering tropical herding on each continent.

Mitogenome Diversity and Maternal Origins of Guangxi Cattle Breeds

Guangxi Province is located in the southernmost region of China, adjacent to the hotspot that is considered to be the putative migratory corridor or domestication area for Chinese indicine cattle. Here, we investigated the evolutionary status of Guangxi native breeds, Longlin (n = 21), Nandan (n = 18), and Weizhou cattle (n = 17) using mitogenome sequencing. Our results show that Bos indicus sub-haplogroup I1a predominates in Guangxi cattle breeds. Population structure by multidimensional-scaling analysis significantly differentiates Weizhou from the other two breeds (Longlin and Nandan). Moreover, the mtDNA haplotype composition and FST values indicate that the formation of Longlin and Nandan breeds may have been affected by Indian indicine, whereas, Weizhou island might have preserved pure Chinese indicine cattle due to its geographical isolation. We speculate that following the initial entry of zebu into southern China, the subsequent introgression of Indian indicine may have influenced the matrilineal origin of local breeds in southwestern China.

Abundant Genetic Diversity of Yunling Cattle Based on Mitochondrial Genome

Simple Summary

Mitochondrial DNA (mtDNA) analysis is a critical tool in assessing the maternal origin, phylogeny and population structure of domestic animals. Yunling cattle are a composite breed that was created by local Yunnan cattle breeds (as a maternal line) 30 years ago. It can be said that Yunling cattle represent important reservoirs of genetic diversity of local Yunnan cattle. Yunling cattle are characterized as a tropical/subtropical breed with fast growth, heat tolerance and parasite resistance, and can survive in a harsh environment and on low-quality roughage. Assessing the genetic characteristics of Yunling cattle is of particular importance for reasonable breeding strategies for Yunling cattle and the design of a local Yunnan cattle conservation program.

Abstract

Yunling cattle are a composite beef cattle breed, combining Brahman (1/2), Murray Grey (1/4) and Local Yunnan cattle (1/4), that was developed in Yunnan, China in the 1980s. Understanding the genetic information of Yunling cattle is of great significance to the development of reasonable breeding strategies for this breed. In the present study, we assessed the current genetic status of Yunling cattle in Yunnan Province (China) by analyzing the variability of the whole mitochondrial genome of 129 individuals. Altogether, 129 sequences displayed 47 different haplotypes. The haplotype diversity and the average number of nucleotide differences were 0.964 and 128.074, respectively. Phylogenetic analyses classified Yunling cattle into seven haplogroups: T1, T2, T3, T4, T6, I1 and I2. Haplogroup I1 was found to be predominant (41.86%), followed by T3 (28.68%). Furthermore, we also identified a novel haplogroup, T6, and defined the sub-haplogroup I1a in Yunling cattle. According to the formation process of Yunling cattle (local Yunnan cattle as the maternal line), the high genetic diversities in the mitochondria of Yunling cattle could be due to the complex maternal origin of local Yunnan cattle. Further studies about local Yunnan breeds are necessary to determine the exact source of haplogroup T6 in Yunling cattle. Our results will be useful for the evaluation and effective management of Yunling cattle.

Maternal genetic and phylogenetic characteristics of domesticated cattle in northwestern China

Northwestern China, an important part of the Silk Road, was the birthplace of the Zhouzu farming culture. The domestication of cattle as an important aspect of farming culture has had a long history in northwestern China. In this study, we assessed the maternal structure and phylogeny of cattle by analyzing the mitochondrial DNA hypervariable segment I (HVS-I) in 698 native cattle from eight areas of northwestern China. The phylogenetic analyses revealed two highly divergent mtDNA clades: clade T, which had four sub-clades (Ta—Td), and clade I. The cattle domesticated from Bos taurus showed a clear dominant distribution pattern in northwestern China. The nucleotide diversity of the Bos indicus clade was lower than that of clades from Bos taurus. In summary, our results suggest that the native cattle of northwestern China were domesticated from two different maternal ancestors, Bos taurus and Bos indicus, which migrated to the central plains of China from the north and south, respectively, with Bos taurus remaining at the edges of the region. The population expansion of the cattle domesticated from Bos taurus occurred in the Longdong region of Gansu Province, and these cattle formed four relatively independent evolutionary branches. Subsequent to this expansion event, Bos indicus migrated from southern to northern China.

Size does matter: Parallel evolution of adaptive thermal tolerance and body size facilitates adaptation to climate change in domestic cattle

The adaptive potential of livestock under a warming climate is increasingly relevant in relation to the growing pressure of global food security. Studies on heat tolerance demonstrate the interplay of adaptation and acclimatization in functional traits, for example, a reduction in body size and enhanced tolerance in response to a warming climate. However, current lack of understanding of functional traits and phylogenetic history among phenotypically distinct populations constrains predictions of climate change impact. Here, we demonstrate evidence of parallel evolution in adaptive tolerance to heat stress in dwarf cattle breeds (DCB, Bos taurus indicus) and compare their thermoregulatory responses with those in standard size cattle breeds (SCB, crossbred, Bos taurus indicus × Bos taurus taurus). We measured vital physiological, hematological, biochemical, and gene expression changes in DCB and SCB and compared the molecular phylogeny using mitochondrial genome (mitogenome) analysis. Our results show that SCB can acclimatize in the short term to higher temperatures but reach their tolerance limit under prevailing tropical conditions, while DCB is adapted to the warmer climate. Increased hemoglobin concentration, reduced cellular size, and smaller body size enhance thermal tolerance. Mitogenome analysis revealed that different lineages of DCB have evolved reduced size independently, as a parallel adaptation to heat stress. The results illustrate mechanistic ways of dwarfing, body size-dependent tolerance, and differential fitness in a large mammal species under harsh field conditions, providing a background for comparing similar populations during global climate change. These demonstrate the value of studies combining functional, physiological, and evolutionary approaches to delineate adaptive potential and plasticity in domestic species. We thus highlight the value of locally adapted breeds as a reservoir of genetic variation contributing to the global domestic genetic resource pool that will become increasingly important for livestock production systems under a warming climate.

Genetic characterization of Vietnamese Yellow cattle using mitochondrial DNA and Y-chromosomal haplotypes and genes associated with economical traits

Vietnamese Yellow cattle are native cattle well adapted to local tropical environment. The aim of this study was to investigate genetic characteristics of the Yellow cattle using molecular markers. We investigated the nucleotide sequences of mitochondrial DNA and SRY gene on Y chromosome, and genotyped SREBP-1, SCD1, EDG1, NCAPG, DGAT1, MC1R, and HSP70 genes in the Yellow cattle population. The sequence analysis of the mitochondrial DNA showed that most of the cattle possesses zebu (Bos indicus) type I1 haplotype, suggesting relatively low genetic diversity in maternal lineage. The sequence analysis of the SRY gene indicates that while most of the males possess zebu type haplotype, taurine (Bos taurus) type haplotype was also observed, suggesting gene-flow from taurine cattle. The results of the genotyping of the functional genes showed that the NCAPG, SCD, MC1R, and HsSP70 genes are polymorphic in the population, whereas the SREBP-1, EDG1, and DGAT1 genes are monomorphic. Particularly, the presence of the desirable and undesirable alleles of the NCAPG and HSP70 genes, respectively, will be important for the selection of animals by potential performances in meat productivity and fertility. The present findings will be informative for future conservation and breeding of the Vietnamese Yellow cattle.

Some maternal lineages of domestic horses may have origins in East Asia revealed with further evidence of mitochondrial genomes and HVR-1 sequences

Objectives

There are large populations of indigenous horse (Equus caballus) in China and some other parts of East Asia. However, their matrilineal genetic diversity and origin remained poorly understood. Using a combination of mitochondrial DNA (mtDNA) and hypervariable region (HVR-1) sequences, we aim to investigate the origin of matrilineal inheritance in these domestic horses.

Methods

To investigate patterns of matrilineal inheritance in domestic horses, we conducted a phylogenetic study using 31 de novo mtDNA genomes together with 317 others from the GenBank. In terms of the updated phylogeny, a total of 5,180 horse mitochondrial HVR-1 sequences were analyzed.

Results

Eightteen haplogroups (Aw-Rw) were uncovered from the analysis of the whole mitochondrial genomes. Most of which have a divergence time before the earliest domestication of wild horses (about 5,800 years ago) and during the Upper Paleolithic (35-10 KYA). The distribution of some haplogroups shows geographic patterns. The Lw haplogroup contained a significantly higher proportion of European horses than the horses from other regions, while haplogroups Jw, Rw, and some maternal lineages of Cw, have a higher frequency in the horses from East Asia. The 5,180 sequences of horse mitochondrial HVR-1 form nine major haplogroups (A-I). We revealed a corresponding relationship between the haplotypes of HVR-1 and those of whole mitochondrial DNA sequences. The data of the HVR-1 sequences also suggests that Jw, Rw, and some haplotypes of Cw may have originated in East Asia while Lw probably formed in Europe.

Conclusions

Our study supports the hypothesis of the multiple origins of the maternal lineage of domestic horses and some maternal lineages of domestic horses may have originated from East Asia.

Variations in mitochondrial cytochrome b region among Ethiopian indigenous cattle populations assert Bos taurus maternal origin and historical dynamics

Objective

This study was carried out to assess the haplotype diversity and population dynamics in cattle populations of Ethiopia.

Methods

We sequenced the complete mitochondrial cytochrome b gene of 76 animals from five indigenous and one Holstein Friesian×Barka cross bred cattle populations.

Results

In the sequence analysis, 18 haplotypes were generated from 18 segregating sites and the average haplotype and nucleotide diversities were 0.7540±0.043 and 0.0010±0.000, respectively. The population differentiation analysis shows a weak population structure (4.55%) among the populations studied. Majority of the variation (95.45%) is observed by within populations. The overall average pair-wise distance (F_ST) was 0.049539 with the highest (F_ST = 0.1245) and the lowest (F_ST = 0.011) F_ST distances observed between Boran and Abigar, and Sheko and Abigar from the indigenous cattle, respectively. The phylogenetic network analysis revealed that all the haplotypes detected clustered together with the Bos taurus cattle and converged to a haplogroup. No haplotype in Ethiopian cattle was observed clustered with the reference Bos indicus group. The mismatch distribution analysis indicates a single population expansion event among the cattle populations.

Conclusion

Overall, high haplotype variability was observed among Ethiopian cattle populations and they share a common ancestor with Bos taurus.

Mitochondrial DNA variants of Podolian cattle breeds testify for a dual maternal origin

BACKGROUND

Over the past 15 years, 300 out of 6000 breeds of all farm animal species identified by the Food and Agriculture Organization of the United Nations (FAO) have gone extinct. Among cattle, many Podolian breeds are seriously endangered in various European areas. Podolian cattle include a group of very ancient European breeds, phenotypically close to the aurochs ancestors (Bos primigenius). The aim of the present study was to assess the genetic diversity of Podolian breeds and to reconstruct their origin.

METHODOLOGY

The mitochondrial DNA (mtDNA) control-regions of 18 Podolian breeds have been phylogenetically assessed. Nine non-Podolian breeds have been also included for comparison.

CONCLUSION

The overall analysis clearly highlights some peculiarities in the mtDNA gene pool of some Podolian breeds. In particular, a principal component analysis point to a genetic proximity between five breeds (Chianina, Marchigiana, Maremmana, Podolica Italiana and Romagnola) reared in Central Italy and the Turkish Grey. We here propose the suggestive hypothesis of a dual ancestral contribution to the present gene pool of Podolian breeds, one deriving from Eastern European cattle; the other arising from the arrival of Middle Eastern cattle into Central Italy through a different route, perhaps by sea, ferried by Etruscan boats. The historical migration of Podolian cattle from North Eastern Europe towards Italy has not cancelled the mtDNA footprints of this previous ancient migration.

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.

Complete mitochondrial genome sequences of Korean native horse from Jeju Island: uncovering the spatio-temporal dynamics

The Korean native horse (Jeju horse) is one of the most important animals in Korean historical, cultural, and economical viewpoints. In the early 1980s, the Jeju horse was close to extinction. The aim of this study is to explore the phylogenomics of Korean native horse focusing on spatio-temporal dynamics. We determined complete mitochondrial genome sequences for the first Korean native (n = 6) and additional Mongolian (n = 2) horses. Those sequences were analyzed together with 143 published ones using Bayesian coalescent approach as well as three different phylogenetic analysis methods, Bayesian inference, maximum likelihood, and neighbor-joining methods. The phylogenomic trees revealed that the Korean native horses had multiple origins and clustered together with some horses from four European and one Middle Eastern breeds. Our phylogenomic analyses also supported that there was no apparent association between breed or geographic location and the evolution of global horses. Time of the most recent common ancestor of the Korean native horse was approximately 13,200-63,200 years, which was much younger than 0.696 My of modern horses. Additionally, our results showed that all global horse lineages including Korean native horse existed prior to their domestication events occurred in about 6000-10,000 years ago. This is the first study on phylogenomics of the Korean native horse focusing on spatio-temporal dynamics. Our findings increase our understanding of the domestication history of the Korean native horses, and could provide useful information for horse conservation projects as well as for horse genomics, emergence, and the geographical distribution.

Identification of a new Y chromosome haplogroup in Spanish native cattle

The aim of this work was to perform a thorough analysis of the diversity of Y-haplotypes in Spanish cattle. A total of 207 Bos taurus males were sampled across 25 European breeds, with a special focus on rare, local Spanish populations. Animals were genotyped with five Y-specific microsatellites (INRA189, UMN0103, UMN0307, BM861 and BYM1), two indels (ZFY10 and USP9Y) and one SNP (UTY19). A new haplogroup, distinct from those described by Götherström et al. (2005), was identified and named Y1.2. Samples representing the three B. taurus Y-haplogroups were genotyped for four additional Y chromosome SNPs (rs121919254, rs121919281, rs121919323 and rs137049553). Among these SNPs, only rs121919281 was informative in B. taurus and helped to confirm the new Y1.2 haplogroup. Analysis of a larger dataset of standardized haplotypes for 1507 individuals from 57 populations from Spain, other European countries and Africa showed the new Y1.2 haplogroup to be found exclusively in Spanish breeds. This finding reinforces the importance of local Spanish cattle as reservoirs of genetic diversity as well as the importance of the Iberian Peninsula in the history of cattle.

Balkan brachicerous cattle - the first domesticated cattle in Europe

The present study aimed to compare mitochondrial diversity among Balkan Neolithic/Chalcolithic cattle and present day Shorthorn Rhodopean cattle (Busha) to throw a new insight into European cattle domestication. The results showed that both ancient- and present-day samples belonged to the macrohaplogroup T. From the 28 sequences (8 ancient and 20 modern), the T1 and T2 haplogroup represent about 3.6% (1/28; 1/28). The T3 haplogroup was with the highest frequency - 57% (16/28). Based on the SNPs on 16057A and 16133C, the new T6 haplogroup was proposed. This haplogroup represents 75% from the ancient and 20% from the present day Bulgarian brachicerous cattle population. The survey in GenBank data base did not find a similar motif, except for the recent Serbian Busha cattle. Overall, these results showed that: (i) The newly named T6 haplogroup is Balkan specific; (ii) The T6 haplogroup survives in present day Bulgarian rhodopean cattle; (iii) The Balkan brachicerous cattle is the oldest European cattle breed.