On the origin of cattle: How aurochs became cattle and colonized the world

Genome-Wide Identification and Evolutionary and Mutational Analysis of the Bos taurus Pax Gene Family

Bos taurus is known for its tolerance of coarse grains, adaptability, high temperature, humidity, and disease resistance. Primarily, cattle are raised for their meat and milk, and pinpointing genes associated with traits relevant to meat production can enhance their overall productivity. The aim of this study was to identify the genome, analyze the evolution, and explore the function of the Pax gene family in B. taurus to provide a new molecular target for breeding in meat-quality-trait cattle. In this study, 44 Pax genes were identified from the genome database of five species using bioinformatics technology, indicating that the genetic relationships of bovids were similar. The Pax3 and Pax7 protein sequences of the five animals were highly consistent. In general, the Pax gene of the buffalo corresponds to the domestic cattle. In summary, there are differences in affinity between the Pax family genes of buffalo and domestic cattle in the Pax1/9, Pax2/5/8, Pax3/7, and Pax4/6 subfamilies. We believe that Pax1/9 has an effect on the growth traits of buffalo and domestic cattle. The Pax3/7 gene is conserved in the evolution of buffalo and domestic animals and may be a key gene regulating the growth of B. taurus. The Pax2/5/8 subfamily affects coat color, reproductive performance, and milk production performance in cattle. The Pax4/6 subfamily had an effect on the milk fat percentage of B. taurus. The results provide a theoretical basis for understanding the evolutionary, structural, and functional characteristics of the Pax family members of B. taurus and for molecular genetics and the breeding of meat-production B. taurus species.

Evolutionary adaptations generally reverse phenotypic plasticity to restore ancestral phenotypes during new environment adaptation in cattle

Phenotype plasticity and evolution adaptations are the two main ways in which allow populations to deal with environmental changes, but the potential relationship between them remains controversial. Using a reciprocal transplant approach with cattle adapted to the Tibetan Plateau and adjacent lowlands, we aim to investigate the relative contributions of evolutionary processes and phenotypic plasticity in driving both phenotypic and transcriptomic changes under natural conditions. We observed that while numerous genetic transcriptomic changes were evident during the forward adaptation to highland environments, plastic changes predominantly facilitate the transformation of transcriptomes into a preferred state when Tibetan cattle are reintroduced to lowland habitats. Genes with ancestral plasticity are generally reversed by evolutionary adaptations and show a closer expression level to the ancestral stage in evolved Tibetan cattle. A similar trend was also observed at the phenotypes level, with a majority of biochemical and hemorheology phenotypes showing a tendency to revert to their ancestral patterns, suggesting the restoration of ancestral expression levels is a widespread evolutionary trend during adaptation. The findings of our study contribute to the debate regarding the relative contributions of plasticity and genetic changes in mammal environment adaptation. Furthermore, we highlight that the restoration of ancestral phenotypes represents a general pattern in cattle new environment adaptation.

Skin transcriptome analysis in Brangus cattle under heat stress

Heat stress is a major factor that negatively affects animal welfare and production systems. Livestock should adapt to tropical and subtropical areas and to meet this, composite breeds have been developed. This work aimed to evaluate gene expression profiles in the skin of Brangus cattle under heat stress using a case-control design, and to correlate this with skin histological characteristics. Two groups of bulls were set using rectal temperature as a criterion to define stress conditions: stressed (N = 5) and non-stressed (N = 5) groups. Skin transcriptomics was performed and correlations between breed composition, phenotypic and skin histological traits were evaluated. Results showed 4309 differentially expressed genes (P < 0.01), 2113 downregulated and 2196 upregulated. Enrichment and ontology analyses revealed 132 GO terms and 67 pathways (P < 0.01), including thermogenesis, glycolysis, gluconeogenesis, mitochondrial activity, antioxidant and immune response, and apoptosis. The identity of the terms and pathways indicated the diversity of mechanisms directed to relieve the animals' suffering, acting from simple passive mechanisms (conduction, convection and radiation) to more complex active ones (behavioural changes, evaporation, vasodilation and wheezing). Furthermore, significant differences between phenotypic and skin histological traits and correlations between pairs of traits suggested a direction towards heat dissipation processes. In this sense, number of vessels was positively correlated with number of sweat glands (P < 0.001) and both were positively correlated with zebuine genetic content (P < 0.05 and P < 0.01, respectively), gland size was positively correlated with epidermal thickness and negatively with hair length (P < 0.05), and epidermal thickness was negatively correlated with gland-epidermis distance (P < 0.0005). These results support the notion that response to heat stress is physiologically complex, producing significant changes in the expression of genes involved in several biological pathways, while the animal's ability to face it depends greatly on their skin features.

MicroRNA breed and parent-of-origin effects provide insights into biological pathways differentiating cattle subspecies in fetal liver

Introduction: MicroRNAs (miRNAs) play a crucial role in regulating gene expression during key developmental processes, including fetal development. Brahman (Bos taurus indicus) and Angus (Bos taurus taurus) cattle breeds represent two major cattle subspecies with strikingly different phenotypes. Methods: We analyzed miRNA expression in liver samples of purebred and reciprocal crosses of Angus and Brahman to investigate breed and parent-of-origin effects at the onset of accelerated fetal growth. Results: We identified eight novel miRNAs in fetal liver samples and 14 differentially expressed miRNAs (DEMs) between purebred samples. Correlation of gene expression modules and miRNAs by breed and parent-of-origin effects revealed an enrichment of genes associated with breed-specific differences in traits such as heat tolerance (Brahman) and fat deposition (Angus). We demonstrate that genes predicted to be targets of DEMs were more likely to be differentially expressed than non-targets (p-value < 0.05). We identified several miRNAs (bta-miR-187, bta-miR-216b, bta-miR-2284c, bta-miR-2285c, bta-miR-2285cp, bta-miR-2419-3p, bta-miR-2419-5p, and bta-miR-11984) that showed similar correlation patterns as bta-miR-2355-3p, which has been associated with the glutamatergic synapse pathway, a key facilitator of heat tolerance. Furthermore, we report Angus-breed-specific miRNAs (bta-miR-2313-5p, btamiR-490, bta-miR-2316, and bta-miR-11990) that may be involved in fat deposition. Finally, we showed that the DEMs identified in fetal liver are involved in Rap1, MAPK, and Ras signalling pathways, which are important for fetal development, muscle development and metabolic traits such as fat metabolism. Conclusion: Our work sheds light on the miRNA expression patterns that contribute to gene expression differences driving phenotypic differences in indicine and taurine cattle.

Comprehensive selection signature analyses in dairy cattle exploiting purebred and crossbred genomic data

The main objective of the current research was to locate, annotate, and highlight specific areas of the bovine genome that are undergoing intense positive selection. Here, we are analyzing selection signatures in crossbred (Bos taurus X Bos indicus), taurine (Bos taurus), and indicine (Bos indicus) cattle breeds. Indicine cattle breeds found throughout India are known for their higher heat tolerance and disease resilience. More breeds and more methods can provide a better understanding of the selection signature. So, we have worked on nine distinct cattle breeds utilizing seven different summary statistics, which is a fairly extensive approach. In this study, we carried out a thorough genome-wide investigation of selection signatures using bovine 50K SNP data. We have included the genotyped data of two taurine, two crossbreds, and five indicine cattle breeds, for a total of 320 animals. During the 1950s, these indicine (cebuine) cattle breeds were exported with the aim of enhancing the resilience of taurine breeds in Western countries. For this study, we employed seven summary statistics, including intra-population, i.e., Tajima's D, CLR, iHS, and ROH and inter-population statistics, i.e., FST, XP-EHH, and Rsb. The NCBI database, PANTHER 17.0, and CattleQTL database were used for annotation after finding the important areas under selection. Some genes, including EPHA6, CTNNA2, NPFFR2, HS6ST3, NPR3, KCNIP4, LIPK, SDCBP, CYP7A1, NSMAF, UBXN2B, UGDH, UBE2K, and DAB1, were shown to be shared by three or more different approaches. Therefore, it gives evidence of the most intense selection in these areas. These genes are mostly linked to milk production and adaptability traits. This study also reveals selection regions that contain genes which are crucial to numerous biological functions, including those associated with milk production, coat color, glucose metabolism, oxidative stress response, immunity and circadian rhythms.

The genome-wide relationships of the critically endangered Quadricorna sheep in the Mediterranean region

Livestock European diffusion followed different human migration waves from the Fertile Crescent. In sheep, at least two diffusion waves have shaped the current breeds' biodiversity generating a complex genetic pattern composed by either primitive or fine-wool selected breeds. Nowadays most of the sheep European breeds derive from the second wave which is supposed to have largely replaced oldest genetic signatures, with the exception of several primitive breeds confined on the very edge of Northern Europe. Despite this, some populations also in the Mediterranean region are characterised by the presence of phenotypic traits considered ancestral such as the policeraty, large horns in the ram, short tail, and a moulting fleece. Italy is home of a large number of local breeds, albeit some are already extinct, others are listed as critically endangered, and among these there is the Quadricorna breed which is a four-horned sheep characterised by several traits considered as ancestral. In this context we genotyped 47 individuals belonging to the Quadricorna sheep breed, a relict and endangered breed, from Central and Southern Italy. In doing so we used the Illumina OvineSNP50K array in order to explore its genetic diversity and to compare it with other 41 breeds from the Mediterranean region and Middle-East, with the specific aim to reconstruct its origin. After retaining 32,862 SNPs following data filtering, the overall genomic architecture has been explored by using genetic diversity indices, Principal Component Analysis (PCA) and admixture analysis, while the genetic relationships and migration events have been inferred using a neighbor-joining tree based on Reynolds' distances and by the maximum likelihood tree as implemented in treemix. The Quadricorna breed exhibit genetic diversity indices comparable with those of most of the other analysed breeds, however, the two populations showed opposing patterns of genetic diversity suggesting different levels of genomic inbreeding and drift (FIS and FROH). In general, all the performed genome-wide analyses returned complementary results, indicating a westward longitudinal cline compatible with human migrations from the Middle-East and several additional genetic footprints which might mirror more recent historical events. Interestingly, among the Italian breeds, the original Quadricorna (QUAD_SA) first separated showing its own ancestral component. In addition, the admixture analysis does not suggest any signal of recent gene exchange with other Italian local breeds, highlighting a rather ancestral purity of this population. On the other hand, both the neighbor-joining tree and the treemix analysis seem to suggest a proximity of the Quadricorna populations to breeds of South-Eastern Mediterranean origin. Although our results do not support a robust link between the genetics of the first wave and the presence of primitive traits, the observed genetic uniqueness together with the inferred phylogeograpic reconstruction would suggest an ancient presence of the Quadricorna breed in the Italian Peninsula. Because of this singularity, urgent conservation actions are needed in order to keep the breed and all related cultural products alive.

Genome-wide diversity and admixture of five indigenous cattle populations from the Tigray region of northern Ethiopia

The Tigray region, where we found around eight per cent of the indigenous cattle population of Ethiopia, is considered as the historic centre of the country, with the ancient pre-Aksumite and Aksumite civilisations in contact with the civilisations of the Fertile Crescent and the Indian subcontinent. Here, we used whole genome sequencing data to characterise the genomic diversity, relatedness, and admixture of five cattle populations (Abergelle, Arado, Begait, Erob, and Raya) indigenous to the Tigray region of Ethiopia. We detected 28 to 29 million SNPs and 2.7 to 2.9 million indels in each population, of which 7% of SNPs and 34% of indels were novel. Functional annotation of the variants showed around 0.01% SNPs and 0.22%-0.27% indels in coding regions. Enrichment analysis of genes overlapping missense private SNPs revealed 20 significant GO terms and KEGG pathways that were shared by or specific to breeds. They included important genes associated with morphology (SCN4A, TAS1R2 and KCNG4), milk yield (GABRG1), meat quality (MMRN2, VWC2), feed efficiency (PCDH8 and SLC26A3), immune response (LAMC1, PCDH18, CELSR1, TLR6 and ITGA5), heat resistance (NPFFR1 and HTR7) and genes belonging to the olfactory gene family, which may be related to adaptation to harsh environments. Tigray indigenous cattle are very diverse. Their genome-wide average nucleotide diversity ranged from 0.0035 to 0.0036. The number of heterozygous SNPs was about 0.6-0.7 times higher than homozygous ones. The within-breed average number of ROHs ranged from 777.82 to 1000.45, with the average sum of the length of ROHs ranging from 122.01 Mbp to 163.88 Mbp. The genomic inbreeding coefficients differed among animals and breeds, reaching up to 10% in some Begait and Raya animals. Tigray indigenous cattle shared a common ancestry with Asian indicine (85.6%-88.7%) and African taurine (11.3%-14.1%) cattle, with very small, if any, European taurine introgression. This study identified high within-breed genetic diversity representing an opportunity for breeding improvement programs and, also, significant novel variants that could increase the number of known cattle variants, an important contribution to the knowledge of domestic cattle genetic diversity.

A Postmortem Case Study-An Analysis of microRNA Patterns in a Korean Native Male Calf (Bos taurus coreanae) That Died of Fat Necrosis

Korean native cattle are highly valued for their rich marbling and flavor. Nonetheless, endeavors to enhance marbling levels can result in obesity, a prevalent contributor to fat necrosis. Fat necrosis is characterized by the formation of necrotic fat masses in the abdominal cavity, which physically puts pressure on affected organs, causing physical torsion or obstruction, resulting in death and consequent economic loss. Pancreatic injuries or diabetes mellitus were reported as factors of fat necrosis in humans; however, the pathogenesis in animals has not been established. In this study, we identified fat necrosis in a 6-month-old Korean native cow and investigated its potential underlying causes. Serum samples were utilized for a microarray analysis of bovine miRNA. Comparative examination of miRNA expression levels between cattle afflicted with fat necrosis and healthy cattle unveiled notable variances in 24 miRNAs, such as bta-miR-26a, bta-miR-29a, bta-miR-30a-5p and bta-miR-181a. Upon conducting miRNA-mediated KEGG pathway analysis, several pathways including the prolactin signal pathway, insulin resistance, autophagy, the insulin-signaling pathway and the FoxO-signaling pathway were found to be significantly enriched in the calf affected by fat necrosis. As a result, this study potentially indicates a potential connection between fat necrosis and diabetes in Korean native cattle.

Ancient Mitochondrial Genomes Provide New Clues to the Origin of Domestic Cattle in China

Cattle are one of the six livestock species that have occupied an important place in Chinese history. Previous ancient DNA studies have indicated that Chinese taurine cattle (Bos taurus taurus) are exotic, but the exact route and diffusion by which they were introduced to China is unknown. In this study, we extracted the mitochondrial genomes of 34 cases of ancient taurine cattle (from the late Neolithic to Qin and Han dynasties) excavated from sites in northern China and the eastern Eurasian steppe, and successfully obtained 14 mitochondrial genome sequences. The results of ancient DNA analysis reveal that with cultural exchange and trade, there was close genetic exchange between domestic taurine cattle in different regions. The haplotypes shared by domestic cattle have genetic continuity, reflecting the strong cultural influence of the large capital city sites such as Taosi, Shimao and Erlitou on the surrounding areas. This study suggests that ancient northern Chinese taurine cattle may have accompanied the westward transmission of agricultural or painted pottery culture and thus had a maternal genetic contribution to modern Tibetan cattle.

Global dispersal and adaptive evolution of domestic cattle: a genomic perspective

Domestic cattle have spread across the globe and inhabit variable and unpredictable environments. They have been exposed to a plethora of selective pressures and have adapted to a variety of local ecological and management conditions, including UV exposure, diseases, and stall-feeding systems. These selective pressures have resulted in unique and important phenotypic and genetic differences among modern cattle breeds/populations. Ongoing efforts to sequence the genomes of local and commercial cattle breeds/populations, along with the growing availability of ancient bovid DNA data, have significantly advanced our understanding of the genomic architecture, recent evolution of complex traits, common diseases, and local adaptation in cattle. Here, we review the origin and spread of domestic cattle and illustrate the environmental adaptations of local cattle breeds/populations.

Genetic Structure Analysis of 155 Transboundary and Local Populations of Cattle (Bos taurus, Bos indicus and Bos grunniens) Based on STR Markers

Every week, 1-2 breeds of farm animals, including local cattle, disappear in the world. As the keepers of rare allelic variants, native breeds potentially expand the range of genetic solutions to possible problems of the future, which means that the study of the genetic structure of these breeds is an urgent task. Providing nomadic herders with valuable resources necessary for life, domestic yaks have also become an important object of study. In order to determine the population genetic characteristics, and clarify the phylogenetic relationships of modern representatives of 155 cattle populations from different regions of the world, we collected a large set of STR data (10,250 individuals), including unique native cattle, 12 yak populations from Russia, Mongolia and Kyrgyzstan, as well as zebu breeds. Estimation of main population genetic parameters, phylogenetic analysis, principal component analysis and Bayesian cluster analysis allowed us to refine genetic structure and provided insights in relationships of native populations, transboundary breeds and populations of domestic yak. Our results can find practical application in conservation programs of endangered breeds, as well as become the basis for future fundamental research.

The demographic history of house mice (Mus musculus domesticus) in eastern North America

The Western European house mouse (Mus musculus domesticus) is a widespread human commensal that has recently been introduced to North America. Its introduction to the Americas is thought to have resulted from the transatlantic movements of Europeans that began in the early 16th century. To study the details of this colonization history, we examine population structure, explore relevant demographic models, and infer the timing of divergence among house mouse populations in the eastern United States using published exome sequences from five North American populations and two European populations. For North American populations of house mice, levels of nucleotide variation were lower, and low-frequency alleles were less common than for European populations. These patterns provide evidence of a mild bottleneck associated with the movement of house mice into North America. Several analyses revealed that one North American population is genetically admixed, which indicates at least two source populations from Europe were independently introduced to eastern North America. Estimated divergence times between North American and German populations ranged between ∼1,000 and 7,000 years ago and overlapped with the estimated divergence time between populations from Germany and France. Demographic models comparing different North American populations revealed that these populations diverged from each other mostly within the last 500 years, consistent with the timing of the arrival of Western European settlers to North America. Together, these results support a recent introduction of Western European house mice to eastern North America, highlighting the effects of human migration and colonization on the spread of an invasive human commensal.

Temporal and geographic dynamics of bovine viral diarrhea virus in American countries

Bovine viral diarrhea virus (BVDV) is a worldwide distributed pathogen of livestock classified into three species, BVDV-1 (Pestivirus A), BVDV-2 (Pestivirus B), and HoBi-like pestivirus (HoBiPeV; Pestivirus H). Despite being considered endemic in several regions of the Americas, the spatiotemporal distribution of BVDV is scarcely known. This study aimed to reconstruct the population dynamics of BVDV in American countries. The analyses performed with the partial 5´UTR gene showed that BVDV-1 and -2 would have started their diversification in the 1670s and 1790s in the United States, whereas HoBiPeV probably emerged in the 1980s in Brazil. No evident geographic clustering was observed in the Bayesian trees, which may indicate that multiple introductions events would have occurred following the first introduction. This study provides new insights into BVDV dynamics, although further analyses including sequences from other American countries and continents will help to expand the knowledge of BVDV evolution and transmission.

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.

Animal board invited review: Genomic-based improvement of cattle in response to climate change

Climate change brings challenges to cattle production, such as the need to adapt to new climates and pressure to reduce greenhouse emissions (GHG). In general, the improvement of traits in current breeding goals is favourably correlated with the reduction of GHG. Current breeding goals and tools for increasing cattle production efficiency have reduced GHG. The same amount of production can be achieved by a much smaller number of animals. Genomic selection (GS) may offer a cost-effective way of using an efficient breeding approach, even in low- and middle-income countries. As climate change increases the intensity of heatwaves, adaptation to heat stress leads to lower efficiency of production and, thus, is unfavourable to the goal of reducing GHG. Furthermore, there is evidence that heat stress during cow pregnancy can have many generation-long lowering effects on milk production. Both adaptation and reduction of GHG are among the difficult-to-measure traits for which GS is more efficient and suitable than the traditional non-genomic breeding evaluation approach. Nevertheless, the commonly used within-breed selection may be insufficient to meet the new challenges; thus, cross-breeding based on selecting highly efficient and highly adaptive breeds may be needed. Genomic introgression offers an efficient approach for cross-breeding that is expected to provide high genetic progress with a low rate of inbreeding. However, well-adapted breeds may have a small number of animals, which is a source of concern from a genetic biodiversity point of view. Furthermore, low animal numbers also limit the efficiency of genomic introgression. Sustainable cattle production in countries that have already intensified production is likely to emphasise better health, reproduction, feed efficiency, heat stress and other adaptation traits instead of higher production. This may require the application of innovative technologies for phenotyping and further use of new big data techniques to extract information for breeding.

The ancestral origin of the critically endangered Quadricorna sheep as revealed by genome-wide analysis

Livestock European diffusion followed different human migration waves from the Fertile Crescent. In sheep, at least two diffusion waves have shaped the current breeds' biodiversity generating a complex genetic pattern composed by either primitive or fine-wool selected breeds. Among primitive breeds, aside from sharing common ancestral genomic components, they also show several traits such as the policeraty, large horns in the ram, short tail, and a moulting fleece, considered as ancestral. Although most of the primitive breeds characterized by these traits are confined on the very edge of Northern Europe, several residual populations are also scattered in the Mediterranean region. In fact, although in Italy a large number of local breeds are already extinct, others are listed as critically endangered, and among these there is the Quadricorna breed which is a four-horned sheep characterized by several ancestral traits. In this context we genotyped 47 individuals belonging to the Quadricorna sheep breed, a relict and endangered breed, from Central and Southern Italy. In doing so we used the Illumina OvineSNP50K array in order to explore its genetic diversity and to compare it with other 33 primitive traits-related, Mediterranean and Middle-East breeds, with the specific aim to reconstruct its origin. After retaining 35,680 SNPs following data filtering, the overall genomic architecture has been explored by using genetic diversity indices, Principal Component Analysis (PCA) and admixture analysis, while the genetic relationships and migration events have been inferred using a neighbor-joining tree based on Reynolds' distances and by the maximum likelihood tree as implemented in treemix. Multiple convergent evidence from all our population genetics analyses, indicated that the two Quadricorna populations differ from all the other Italian breeds, while they resulted to be very close to the Middle Eastern and primitive European breeds. In addition, the genetic diversity indices highlighted values comparable with those of most of the other analyzed breeds, despite the two populations exhibit slightly different genetic indices suggesting different levels of genomic inbreeding and drift (FIS and FROH). The admixture analysis does not suggest any signal of recent gene exchange with other Italian local breeds, highlighting a rather ancestral purity of the two populations, while on the other hand the treemix analysis seems to suggest an ancient admixture with other primitive European breeds. Finally, all these evidences seem to trace back the residual Quadricorna sheep to an early Neolithic spread, probably following a Mediterranean route and that urgent conservation actions are needed in order to keep the breed and all related cultural products alive.

Exploring AI-2-mediated interspecies communications within rumen microbial communities

BACKGROUND

The rumen is an ecosystem with a complex microbial microflora in which microbes initiate biofilm formation by attaching to plant surfaces for plant degradation and are capable of converting feed to nutrients and energy via microbial processes. Quorum sensing (QS) is a cell-to-cell communication mechanism that allows microbes to synchronize the expression of multiple genes in the group to perform social behaviors such as chemotaxis and biofilm formation using self-synthesized QS signaling molecules. Whereas QS has been extensively studied in model microorganisms under pure culture conditions, QS mechanisms are poorly understood in complex bacterial communities, such as the rumen microflora, in which cell-to-cell communication may be common.

RESULTS

Here, we analyzed 981 rumens bacterial and archaeal genomes from the Joint Genome Institute (JGI) and GenBank databases and identified 15 types of known QS signaling molecule-related genes. The analysis of the prevalence and abundance of genes involved in QS showed that 767 microbial genomes appeared to possess QS-related genes, including 680 bacterial genomes containing autoinducer-2 (AI-2) synthase- or receptor-encoding genes. Prevotella, Butyivibrio, Ruminococcus, Oribacterium, Selenomonas, and Treponema, known abundant bacterial genera in the rumen, possessed the greatest numbers of AI-2-related genes; these genes were highly expressed within the metatranscriptome dataset, suggesting that intra- and interspecies communication mediated by AI-2 among rumen microbes was universal in the rumen. The QS processes mediated by the dCache_1-containing AI-2 receptors (CahRs) with various functional modules may be essential for degrading plants, digesting food, and providing energy and nutrients to the host. Additionally, a universal natural network based on QS revealed how rumen microbes coordinate social behaviors via the AI-2-mediated QS system, most of which may potentially function via AI-2 binding to the extracellular sensor dCache_1 domain to activate corresponding receptors involved in different signal transduction pathways, such as methyl-accepting chemotaxis proteins, histidine kinases, serine phosphatases, c-di-GMP synthases and phosphodiesterases, and serine/threonine kinases in the rumen.

CONCLUSIONS

The exploration of AI-2-related genes, especially CahR-type AI-2 receptors, greatly increased our insight into AI-2 as a potentially "universal" signal mediating social behaviors and will help us better understand microbial communication networks and the function of QS in plant-microbe interactions in complex microecosystems. Video Abstract.

Using animal history to inform current debates in gene editing farm animals: A systematic review

There is growing interest in gene editing farm animals. Some alterations could benefit animal welfare (e.g., improved heat tolerance in cattle with the “slick” gene), the environment (e.g., reducing methane emissions from cattle with induced pluripotent stem cells), and productivity (e.g., higher weight gains in cattle with the “double muscling” gene). Existing scholarship on the acceptability of such modifications has used myriad approaches to identify societal factors that shape the ethics and governance of this technology. We argue that integrating historical approaches—particularly from the relatively new and burgeoning field of animal history—offers a form of “anticipatory knowledge” that can help guide discussions on this topic. We conducted a systematic review of the animal history literature in English, German, and Spanish to identify the influence of political, scientific, economic, social, and cultural factors on the development and acceptance of such technologies. We identified analogous structures and fault lines in past debates about farm animals that provide insights for contemporary discussions about gene editing. Those analogous structures include the market power of meatpackers or the racialized precepts in livestock breeding, and fault lines, like the disconnect between states and citizens over the direction of food systems. Highlighting these similarities demonstrates how external forces have shaped—and will continue to shape—the acceptance or rejection of emerging biotechnologies as applied to farm animals.

Phylogenomics and Spatiotemporal Dynamics of Bovine Leukemia Virus Focusing on Asian Native Cattle: Insights Into the Early Origin and Global Dissemination

Bovine leukemia virus (BLV), the causative agent of enzootic bovine leukosis, is currently one of the most important pathogens affecting the cattle industry worldwide. Determining where and in which host it originated, and how it dispersed across continents will provide valuable insights into its historical emergence as the cattle pathogen. Various species in the Bos genus were domesticated in Asia, where they also diversified. As native cattle (taurine cattle, zebu cattle, yak, and water buffalo) are indigenous and adapted to local environments, we hypothesized that Asian native cattle could have harbored BLV and, therefore, that they were important for virus emergence, maintenance, and spread. In this study, phylogeographic and ancestral trait analyses—including sequences obtained from Asian native cattle—were used to reconstruct the evolutionary history of BLV. It was shown that, since its probable emergence in Asia, the virus spread to South America and Europe via international trade of live cattle. It was inferred that zebu cattle were the hosts for the early origin of BLV, while taurine cattle played the significant role in the transmission worldwide. In addition, the results of positive selection analysis indicate that yak had a substantially minor role in the transmission of this virus. In this study, endogenous deltaretrovirus sequences in bats, collected in Asian countries, were also analyzed on whether these sequences were present in the bat genome. Endogenous deltaretrovirus sequences were detected from bat species endemic to specific regions and geographically isolated for a long time. Endogenous deltaretrovirus sequences from these geographically isolated species represent ancient exogenous deltaretroviruses distributions. The phylogenetic analysis revealed that these newly obtained endogenous deltaretrovirus sequences were closely related to those of BLV from Asian native cattle, indicating that BLV-related ancient deltaretroviruses circulated in Asia long before the emergence of BLV. Together, our analyses provide evidence for origin and spatiotemporal dynamics of BLV.

Comparing the maternal-fetal microbiome of humans and cattle: a translational assessment of the reproductive, placental, and fetal gut microbiomes

An analysis of sites within the maternal reproductive microbiome that potentially contribute to fetal gut microbial colonization, with a special focus on the comparison between humans and cattle.

Recovery of mitogenomes from whole genome sequences to infer maternal diversity in 1883 modern taurine and indicine cattle

Maternal diversity based on a sub-region of mitochondrial genome or variants were commonly used to understand past demographic events in livestock. Additionally, there is growing evidence of direct association of mitochondrial genetic variants with a range of phenotypes. Therefore, this study used complete bovine mitogenomes from a large sequence database to explore the full spectrum of maternal diversity. Mitogenome diversity was evaluated among 1883 animals representing 156 globally important cattle breeds. Overall, the mitogenomes were diverse: presenting 11 major haplogroups, expanding to 1309 unique haplotypes, with nucleotide diversity 0.011 and haplotype diversity 0.999. A small proportion of African taurine (3.5%) and indicine (1.3%) haplogroups were found among the European taurine breeds and composites. The haplogrouping was largely consistent with the population structure derived from alternate clustering methods (e.g. PCA and hierarchical clustering). Further, we present evidence confirming a new indicine subgroup (I1a, 64 animals) mainly consisting of breeds originating from China and characterised by two private mutations within the I1 haplogroup. The total genetic variation was attributed mainly to within-breed variance (96.9%). The accuracy of the imputation of missing genotypes was high (99.8%) except for the relatively rare heteroplasmic genotypes, suggesting the potential for trait association studies within a breed.

Identification of important genomic footprints using eight different selection signature statistics in domestic cattle breeds

In the present study, the population genomic data of different cattle breeds were explored to decipher the genomic regions affected due to selective events and reflected in the productive, reproductive, thermo-tolerance, and health-related traits. To find out these genomic deviations due to selective sweeps, we used eight different statistical tools (Tajima's D, Fu & Li's D*, CLR, ROH, iHS, F_ST, FLK, and hapFLK) on seven indigenous and five exotic cattle breeds. We further performed composite analysis by comparing their covariance matrix. Several candidate genes were found to be related to milk production (ADARB, WDR70, and CA8), reproductive (PARN, FAM134B2, and ZBTB20), and health-related traits (SP110, CXCL2, CLXCL3, CXCL5, IRF8, and MYOM1). The outcome of this investigation provides a basis for detecting selective sweeps that explain the genetic variation of traits. They may possess functional importance for multiple cattle breeds in different subcontinents. However, further studies are required to improve the findings using high-density arrays or whole-genome sequencing with higher resolution and greater sample sizes.

Temporal analysis of bovine pestivirus diversity in Brazil

In this study, phylogenetic and evolutionary analyses of cattle pestiviruses (BVDV-1, 2 and HoBiPeV) originating in Brazil were used to investigate the temporal diversification of subgenotypes in the country. Inferred dated phylogeny and time of the most recent common ancestor (tMRCA) demonstrated that some BVDV subgenotypes (1a, 1b, 1d, 1e, and 2b) and HoBi-like sequences clustered according to the region in which they were collected and that the diversification of subgenotypes appears to have occurred around the introduction of first Bos taurus and then Bos indicus, followed by expansion to form the adapted Brazilian breeds. The present results help to elucidate the temporal facts that led to diversification of ruminant pestiviruses in cattle in Brazil.

A cattle graph genome incorporating global breed diversity

Despite only 8% of cattle being found in Europe, European breeds dominate current genetic resources. This adversely impacts cattle research in other important global cattle breeds, especially those from Africa for which genomic resources are particularly limited, despite their disproportionate importance to the continent's economies. To mitigate this issue, we have generated assemblies of African breeds, which have been integrated with genomic data for 294 diverse cattle into a graph genome that incorporates global cattle diversity. We illustrate how this more representative reference assembly contains an extra 116.1 Mb (4.2%) of sequence absent from the current Hereford sequence and consequently inaccessible to current studies. We further demonstrate how using this graph genome increases read mapping rates, reduces allelic biases and improves the agreement of structural variant calling with independent optical mapping data. Consequently, we present an improved, more representative, reference assembly that will improve global cattle research.

Bos taurus haplotypes segregating in Nellore (Bos indicus) cattle

Brazil is the largest exporter of beef in the world, and most of that beef derives from Nellore cattle. Although considered a zebu breed (Bos indicus), the history of Nellore cattle in Brazil is marked by the importation of bulls from India, the use of a Creole taurine (Bos taurus) maternal lineage to quickly expand the herds and backcrossing to Nellore bulls to recover zebu ancestry. As a consequence, the current Brazilian Nellore population carries an average taurine ancestry of approximately 1%. Although that percentage seems small, some taurine variants deviate substantially from that average, with the better-known cases being the PLAG1-Q haplotype involved with body size variation and the Guarani (P_G ) polled variant producing hornless animals. Here, we report taurine haplotypes in 9074 Nellore animals genotyped for 539 657 imputed SNP markers. Apart from PLAG1-Q and P_G , our analysis further revealed common taurine haplotypes (>3%) spanning genes related to immunity, growth, reproduction and hair and skin phenotypes. Using data from 22 economically important traits, we showed that many of the major QTL previously reported in the breed are at least partially driven by taurine haplotypes. As B. taurus and B. indicus haplotypes are highly divergent, presenting widely different sets of functional variants, our results provide promising targets for future scrutiny in Nellore cattle.

Ancient Mitogenomes Reveal the Domestication and Distribution of Cattle During the Longshan Culture Period in North China

Cattle, as an important tool for agricultural production in ancient China, have a complex history of domestication and distribution in China. Although it is generally accepted that ancient Chinese taurine cattle originated from the Near East, the explanation regarding their spread through China and whether or not this spread was associated with native aurochs during ancient times are still unclear. In this study, we obtained three nearly complete mitochondrial genomes (mitogenomes) from bovine remains dating back ca. 4,000 years at the Taosi and Guchengzhai sites in North China. For the first time at the mitogenome level, phylogenetic analyses confirmed the approximately 4,000-year-old bovines from North China as taurine cattle. All ancient cattle from both sites belonged to the T3 haplogroup, suggesting their origin from the Near East. The high affinity between ancient samples and southern Chinese taurine cattle indicated that ancient Chinese cattle had a genetic contribution to the taurine cattle of South China. A rapid decrease in the female effective population size ca. 4.65 thousand years ago (kya) and a steep increase ca. 1.99 kya occurred in Chinese taurine cattle. Overall, these results provide increasing evidence of the origin of cattle in the middle Yellow River region of China.

Effects of Temperament on the Reproduction of Beef Cattle

Simple Summary

In beef cattle and other livestock species, temperament can be considered as their behavioral responses to human interaction. Temperament evaluation allows classifying cattle according to the level of responsiveness, from calm to excitable. Consistently across studies, beef females classified as more excitable have signs of increased stress response and are less fertile compared to cohorts with calmer temperament. Cattle temperament can be improved via genetic selection and acclimating young animals to human interaction and handling procedures.

Abstract

Temperament is often defined as the behavioral expression of animals in response to human interaction. Cattle temperament can be evaluated using an association of chute score and exit velocity, with cattle then classified as adequate or excitable temperament. To assess the impacts of temperament on various beef systems, these evaluation criteria were associated with productive and reproductive parameters of Bos taurus and B. indicus-influenced cattle. Consistently across studies, excitable cattle had greater plasma cortisol compared to animals with adequate temperament. Studies also reported that excitable beef females have poorer reproductive performance compared to calmer cohorts, including reduced annual pregnancy rates, decreased calving rate, weaning rate, and kg of calf weaned/cow exposed to breeding. Acclimating B. indicus × B. taurus or B. taurus heifers to human handling improved behavioral expression of temperament and hastened puberty attainment. However, similar benefits were not observed when mature cows were acclimated to human handling. Collectively, temperament of beef females measured via behavioral responses upon human handling impacts their reproductive and productive responses independent of breed type, and should be considered for optimal beef cattle production.

Humans, megafauna and landscape structure – Rock engravings from Namibia encourage a comparative approach to central Europe and southern Africa

This paper deals with reflections that arose after observing prehistoric rock engravings at different locations in Namibia. These observations stimulated comparative considerations with focus on southern Africa and central Europe. Similar to the Aurignacian rock art of European origin, the most common motifs in the Namibian rock engravings are large animals. While in Europe, the species that served as a blueprint for the illustration of Aurignacian rock art have mostly disappeared, the megafauna illustrated on the rock engravings in Namibia can still be found in the immediate vicinity of the rock art. Against this background, we discuss and further develop a comparative regional approach. We reconstruct and evaluate the suitability of African savannas and still-existing megafauna communities as an appropriate reference-frame for natural European grassland systems and extinct associated warm-adapted megafauna (Eemian Interglacial megafauna). Special attention is laid on the unique situation in Africa in the light of a global extinction wave of megafauna following increasing human activity in the Late Quaternary. This leads us to discuss the use of domesticated ungulates as surrogate taxa to fulfill ecosystem functions in Europe as part of concepts termed “rewilding” or “naturalistic grazing”. After critically examining these concepts, we conclude that using domesticated forms as representatives of extinct or locally disappeared species in Europe has its justification to some extent. If, however, the naturally occurring megaherbivore community still exists (Africa), these naturally occurring species should be given priority due to their organismic abilities and limitations adapted to the harsh conditions in their specific environment. Finally, we discuss the application of (transboundary) protected areas as effective instruments to mitigate human-wildlife conflicts. A holistic approach, including nature conservation and preservation of cultural achievements (domesticated forms, grazing systems), appears promising for the effective protection of the natural African savanna ecosystems with their unique fauna elements, as illustrated in rock engravings that inspired us to write this paper.

Genomic diversity and phylodynamic of bovine viral diarrhea virus in Argentina

Bovine viral diarrhea virus (BVDV) is an important pathogen of ruminants worldwide and is characterized by high genetic diversity and a wide range of clinical presentations. In Argentina, several studies have evaluated the genetic diversity of BVDV but no phylodynamic study has been published yet. In this study, a comprehensive compilation and update of Argentinean BVDV sequences were performed, and the evolutionary history of BVDV was characterized by phylodynamic analyses based on the 5´UTR. Although BVDV-1b and BVDV-1a were the most frequent subtypes, novel subtypes for Argentina, 1e and 1i, were identified. The phylodynamic analysis suggested that BVDV started its diversification in the mid-1650s with an exponential increase in viral diversity since the late 1990s, possibly related to the livestock expansion and intensification in the country. Evolutionary rate in the 5´UTR was faster for BVDV-1a than for BVDV-1b, and both subtypes presented an endemic nature according to the demographic reconstructions. The current study contributes to clarify the evolutionary history of BVDV in the main cattle region of the country and provides useful information about the epidemiology and future development of diagnostic and control tools in Argentina.

A Review of Australian Tick Vaccine Research

Tick vaccine research in Australia has demonstrated leadership worldwide through the development of the first anti-tick vaccine in the 1990s. Australia’s Commonwealth Scientific and Industrial Research Organisation’s (CSIRO) research led to the development of vaccines and/or precursors of vaccines (such as crude extracts) for both the cattle tick and the paralysis tick. CSIRO commercialised the Bm86 vaccine in the early 1990s for Rhipicephalus australis; however, issues with dosing and lack of global conservation led to the market closure of Tick-GARD in Australia. New research programs arose both locally and globally. The Australian paralysis tick Ixodes holocyclus has perplexed research veterinarians since the 1920s; however, not until the 2000s did biotechnology exist to elucidate the neurotoxin—holocyclotoxin family of toxins leading to a proof of concept vaccine cocktail. This review revisits these discoveries and describes tributes to deceased tick vaccine protagonists in Australia, including Sir Clunies Ross, Dr Bernard Stone and Dr David Kemp.

Genome-Wide Analyses Reveal the Genetic Architecture and Candidate Genes of Indicine, Taurine, Synthetic Crossbreds, and Locally Adapted Cattle in Brazil

Cattle population history, breeding systems, and geographic subdivision may be reflected in runs of homozygosity (ROH), effective population size (N_e), and linkage disequilibrium (LD) patterns. Thus, the assessment of this information has become essential to the implementation of genomic selection on purebred and crossbred cattle breeding programs. In this way, we assessed the genotype of 19 cattle breeds raised in Brazil belonging to taurine, indicine, synthetic crossbreds, and Iberian-derived locally adapted ancestries to evaluate the overall LD decay patterns, N_e, ROH, and breed composition. We were able to obtain a general overview of the genomic architecture of cattle breeds currently raised in Brazil and other tropical countries. We found that, among the evaluated breeds, different marker densities should be used to improve the genomic prediction accuracy and power of genome-wide association studies. Breeds showing low N_e values indicate a recent inbreeding, also reflected by the occurrence of longer ROH, which demand special attention in the matting schemes to avoid extensive inbreeding. Candidate genes (e.g., ABCA7, PENK, SPP1, IFNAR1, IFNAR2, SPEF2, PRLR, LRRTM1, and LRRTM4) located in the identified ROH islands were evaluated, highlighting biological processes involved with milk production, behavior, rusticity, and fertility. Furthermore, we were successful in obtaining the breed composition regarding the taurine and indicine composition using single-nucleotide polymorphism (SNP) data. Our results were able to observe in detail the genomic backgrounds that are present in each breed and allowed to better understand the various contributions of ancestor breeds to the modern breed composition to the Brazilian cattle.

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.

Genome-Wide Scanning for Signatures of Selection Revealed the Putative Genomic Regions and Candidate Genes Controlling Milk Composition and Coat Color Traits in Sahiwal Cattle

Background

In the evolutionary time scale, selection shapes the genetic variation and alters the architecture of genome in the organisms. Selection leaves detectable signatures at the genomic coordinates that provide clues about the protein-coding regions. Sahiwal is a valuable indicine cattle adapted to tropical environments with desirable milk attributes. Insights into the genomic regions under putative selection may reveal the molecular mechanisms affecting the quantitative and other important traits. To understand this, the present investigation was undertaken to explore signatures of selection in the genome of Sahiwal cattle using a medium-density genotyping INDUS chip.

Result

De-correlated composite of multiple selection signals (DCMS), which combines five different univariate statistics, was computed in the dataset to detect the signatures of selection in the Sahiwal genome. Gene annotations, Quantitative Trait Loci (QTL) enrichment, and functional analyses were carried out for the identification of significant genomic regions. A total of 117 genes were identified, which affect a number of important economic traits. The QTL enrichment analysis highlighted 14 significant [False Discovery Rate (FDR)-corrected p-value ≤ 0.05] regions on chromosomes BTA 1, 3, 6, 11, 20, and 21. The top three enriched QTLs were found on BTA 6, 20, and 23, which are associated with exterior, health, milk production, and reproduction traits. The present study on selection signatures revealed some key genes related with coat color (PDGFRA, KIT, and KDR), facial pigmentation (LEF), milk fat percent (MAP3K1, HADH, CYP2U1, and SGMS2), sperm membrane integrity (OSTC), lactation persistency (MRPS30, NNT, CCL28, HMGCS1, NIM1K, ZNF131, and CCDC152), milk yield (GHR and ZNF469), reproduction (NKX2-1 and DENND1A), and bovine tuberculosis susceptibility (RNF144B and PAPSS1). Further analysis of candidate gene prioritization identified four hub genes, viz., KIT, KDR, MAP3K1, and LEF, which play a role in coat color, facial pigmentation, and milk fat percentage in cattle. Gene enrichment analysis revealed significant Gene ontology (GO) terms related to breed-specific coat color and milk fat percent.

Conclusion

The key candidate genes and putative genomic regions associated with economic traits were identified in Sahiwal using single nucleotide polymorphism data and the DCMS method. It revealed selection for milk production, coat color, and adaptability to tropical climate. The knowledge about signatures of selection and candidate genes affecting phenotypes have provided a background information that can be further utilized to understand the underlying mechanism involved in these traits in Sahiwal cattle.

High-Level Gene Flow Restricts Genetic Differentiation in Dairy Cattle Populations in Thailand: Insights from Large-Scale Mt D-Loop Sequencing

Domestication and artificial selection lead to the development of genetically divergent cattle breeds or hybrids that exhibit specific patterns of genetic diversity and population structure. Recently developed mitochondrial markers have allowed investigation of cattle diversity worldwide; however, an extensive study on the population-level genetic diversity and demography of dairy cattle in Thailand is still needed. Mitochondrial D-loop sequences were obtained from 179 individuals (hybrids of Bos taurus and B. indicus) sampled from nine different provinces. Fifty-one haplotypes, of which most were classified in haplogroup "I", were found across all nine populations. All sampled populations showed severely reduced degrees of genetic differentiation, and low nucleotide diversity was observed in populations from central Thailand. Populations that originated from adjacent geographical areas tended to show high gene flow, as revealed by patterns of weak network structuring. Mismatch distribution analysis was suggestive of a stable population, with the recent occurrence of a slight expansion event. The results provide insights into the origins and the genetic relationships among local Thai cattle breeds and will be useful for guiding management of cattle breeding in Thailand.

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.

Herd clustering strategies and corresponding genetic evaluations based on social-ecological characteristics for a local endangered cattle breed

The accuracy of breeding values strongly depends on the population and herd structure, i.e., the number of animals considered in genetic evaluations and the size of contemporary groups (CGs). Local breeds are usually kept in small-sized family farms under alternative husbandry conditions. For such herd structure, consideration of classical herd or herd-test-day effects in CG modeling approaches implies only a few records per effect level. In consequence, the present study aimed on methodological evaluations of different herd clustering strategies, considering social–ecological and herd characteristics. In this regard, we considered 19 herds keeping cows from the small local population of German Black Pied cattle (Deutsches Schwarzbuntes Niederungsrind; DSN), 10 herds keeping Holstein Friesian (HF) cows and one mixed herd with HF and DSN cows. Herds were characterized for 106 variables, reflecting farm conditions, husbandry practices, feeding regime, herd management, herd fertility status, herd health status and breeding strategies as well as social–ecological descriptors. The variables were input data for different clustering approaches including agglomerative hierarchical clustering (AHC), partition around medoids (PAM), fuzzy clustering (FZC) and a clustering of variables combined with agglomerative hierarchical clustering (CoVAHC). The evaluation criterion was the average silhouette width (ASW), suggesting a CoVAHC application and consideration of four herd clusters (HCs) for herd allocation (ASW of 0.510). HC1 comprised the larger, half organic and half conventional DSN family farms, which generate their main income from milk production. HC2 consisted of small organic DSN family farms where cows are kept in tie stables. HC3 included the DSN sub-population from former East Germany, reflecting the large-scale farm types. The specialized HF herds were well separated and allocated to HC4. Generalized linear mixed models with appropriate link functions were applied to compare test-day and female fertility traits of 5538 cows (2341 DSN and 3197 HF) from the first three lactations among the four HCs. Least squares means for milk, fat and protein yield (Mkg, Fkg and Pkg) significantly differed between HC. The significant differences among the four HCs clearly indicate the influence of varying herd conditions on cow traits. The similarities of herds within HC suggested the application of HCs in statistical models for genetic evaluations for DSN. In this regard, we found an increase of accuracies of estimated breeding values of cows and sires and of heritabilities for milk yield when applying models with herd-cluster-test-day or herd-cluster-test-month effects compared to classical herd-test-day models. The identified increase for the number of cows and cow records in CG due to HC effects may be the major explanation for the identified superiority.

Variants at the ASIP locus contribute to coat color darkening in Nellore cattle

BACKGROUND

Nellore cattle (Bos indicus) are well-known for their adaptation to warm and humid environments. Hair length and coat color may impact heat tolerance. The Nellore breed has been strongly selected for white coat, but bulls generally exhibit darker hair ranging from light grey to black on the head, neck, hump, and knees. Given the potential contribution of coat color variation to the adaptation of cattle populations to tropical and sub-tropical environments, our aim was to map positional and functional candidate genetic variants associated with darkness of hair coat (DHC) in Nellore bulls.

RESULTS

We performed a genome-wide association study (GWAS) for DHC using data from 432 Nellore bulls that were genotyped for more than 777 k single nucleotide polymorphism (SNP) markers. A single major association signal was detected in the vicinity of the agouti signaling protein gene (ASIP). The analysis of whole-genome sequence (WGS) data from 21 bulls revealed functional variants that are associated with DHC, including a structural rearrangement involving ASIP (ASIP-SV1). We further characterized this structural variant using Oxford Nanopore sequencing data from 13 Australian Brahman heifers, which share ancestry with Nellore cattle; we found that this variant originates from a 1155-bp deletion followed by an insertion of a transposable element of more than 150 bp that may impact the recruitment of ASIP non-coding exons.

CONCLUSIONS

Our results indicate that the variant ASIP sequence causes darker coat pigmentation on specific parts of the body, most likely through a decreased expression of ASIP and consequently an increased production of eumelanin.

Human adaptation, demography and cattle domestication: an overview of the complexity of lactase persistence in Africa

Lactase persistence (LP) is a genetically-determined trait that is prevalent in African, European and Arab populations with a tradition of animal herding and milk consumption. To date, genetic analyses have identified several common variants that are associated with LP. Furthermore, data have indicated that these functional alleles likely have been maintained in pastoralist populations due to the action of recent selection, exemplifying the ongoing evolution of anatomically modern humans. Additionally, demographic history has also played a role in the geographic distribution of LP and associated alleles in Africa. In particular, the migration of ancestral herders and their subsequent admixture with local populations were integral to the spread of LP alleles and the culture of pastoralism across the continent. The timing of these demographic events was often correlated with known major environmental changes and/or the ability of domesticated cattle to resist/avoid infectious diseases. This review summarizes recent advances in our understanding of the genetic basis and evolutionary history of LP, as well as the factors that influenced the origin and spread of pastoralism in Africa.

Researching Human-Cattle Interaction on Rangelands: Challenges and Potential Solutions

Simple Summary

Research investigating the influence of human-animal interactions on beef cattle production and welfare often occurs during stress-inducing contexts, such as restraint and transport. Given that beef cattle which graze on expansive rangelands do not frequently interact with humans, assessment of animal response to humans may be inaccurate if they are only recorded during such stressful context. Future research aimed at identifying the causes and impacts of individual human–cattle relationships may consider conducting experiments while the cattle are in a low stress context, such as while ruminating, resting, or, grazing. By adapting methods used to research wildlife response to humans, researchers may be able to better understand how animal personality and previous experience can influence cattle production and welfare.

Abstract

Human-animal interaction (HAI) research spans across many scientific fields and animal taxa. For livestock species, HAI research tends to focus on animals that are managed in close proximity with humans such as poultry, dairy cattle, and swine. Given the nature of rangeland cattle production, HAI research with beef cattle often occurs in and around the processing environment. This high arousal context may skew behavioral and physiological responses by the animals due to the potentially negative interaction. The aim of this review is to describe cattle production on rangelands, examine the considerations and limitations of current HAI research used to evaluate interaction quality or traits of rangeland cattle, identify contexts in which rangeland cattle interact with humans, and provide recommendations for improving future HAI research with rangeland cattle. Current research delineating individual differences in response to humans by beef cattle occur during routine husbandry and management on rangelands (pragmatic) and in a research context (experimental). Human-cattle interactions can be distinguished based on the quality and goal of the interaction into four broad categories: human presence, human approach, human contact, and restraint. Limitations of HAI research with rangeland cattle are identified and reconciled by recommendations for HAI research that can take place outside of the processing environment (i.e., while cattle are ruminating, resting or grazing on rangelands).

Fit for the Job: Proportion and the Portrayal of Cattle in Egyptian Old and Middle Kingdom Elite Tomb Imagery

Depictions of the natural world are an intrinsic feature of Egyptian visual culture, with the vast array of imagery documenting animals a testimony to the fundamental role they played. Despite the significance of animals in Egypt, an anthropocentric bias still exists in research on the methods used by practitioners during initial scene composition. To help bridge the divide, the author herein undertook an investigation to determine if proportional guides were in place when rendering animal figures in ancient Egyptian elite tomb imagery of the Old and Middle Kingdoms. A notable outcome of the proportional analysis was the identification of two distinct body-types for domestic cattle (Bos taurus taurus). The aim of the current paper is to further examine these proportional differences to explore if variations in physique (namely the distance between the chest floor and withers) were rendered by Egyptian practitioners to reflect the conditions in which they appeared by considering two overarching factors: (1) biological factors and (2) contextual factors. As such, the study will employ proportional analysis to challenge the prevailing perspective of a deregulated approach when illustrating fauna in elite tomb imagery, highlighting the significance of animals within ancient Egypt.

Genetic Diversity and Signatures of Selection for Thermal Stress in Cattle and Other Two Bos Species Adapted to Divergent Climatic Conditions

Understanding the biological mechanisms of climatic adaptation is of paramount importance for the optimization of breeding programs and conservation of genetic resources. The aim of this study was to investigate genetic diversity and unravel genomic regions potentially under selection for heat and/or cold tolerance in thirty-two worldwide cattle breeds, with a focus on Chinese local cattle breeds adapted to divergent climatic conditions, Datong yak (Bos grunniens; YAK), and Bali (Bos javanicus) based on dense SNP data. In general, moderate genetic diversity levels were observed in most cattle populations. The proportion of polymorphic SNP ranged from 0.197 (YAK) to 0.992 (Mongolian cattle). Observed and expected heterozygosity ranged from 0.023 (YAK) to 0.366 (Sanhe cattle; SH), and from 0.021 (YAK) to 0.358 (SH), respectively. The overall average inbreeding (±SD) was: 0.118 ± 0.028, 0.228 ± 0.059, 0.194 ± 0.041, and 0.021 ± 0.004 based on the observed versus expected number of homozygous genotypes, excess of homozygosity, correlation between uniting gametes, and runs of homozygosity (ROH), respectively. Signatures of selection based on multiple scenarios and methods (F ST, HapFLK, and ROH) revealed important genomic regions and candidate genes. The candidate genes identified are related to various biological processes and pathways such as heat-shock proteins, oxygen transport, anatomical traits, mitochondrial DNA maintenance, metabolic activity, feed intake, carcass conformation, fertility, and reproduction. This highlights the large number of biological processes involved in thermal tolerance and thus, the polygenic nature of climatic resilience. A comprehensive description of genetic diversity measures in Chinese cattle and YAK was carried out and compared to 24 worldwide cattle breeds to avoid potential biases. Numerous genomic regions under positive selection were detected using three signature of selection methods and candidate genes potentially under positive selection were identified. Enriched function analyses pinpointed important biological pathways, molecular function and cellular components, which contribute to a better understanding of the biological mechanisms underlying thermal tolerance in cattle. Based on the large number of genomic regions identified, thermal tolerance has a complex polygenic inheritance nature, which was expected considering the various mechanisms involved in thermal stress response.

Breed Ancestry, Divergence, Admixture, and Selection Patterns of the Simbra Crossbreed

In this study, we evaluated an admixed South African Simbra crossbred population, as well as the Brahman (Indicine) and Simmental (Taurine) ancestor populations to understand their genetic architecture and detect genomic regions showing signatures of selection. Animals were genotyped using the Illumina BovineLD v2 BeadChip (7K). Genomic structure analysis confirmed that the South African Simbra cattle have an admixed genome, composed of 5/8 Taurine and 3/8 Indicine, ensuring that the Simbra genome maintains favorable traits from both breeds. Genomic regions that have been targeted by selection were detected using the linkage disequilibrium-based methods iHS and Rsb. These analyses identified 10 candidate regions that are potentially under strong positive selection, containing genes implicated in cattle health and production (e.g., TRIM63, KCNA10, NCAM1, SMIM5, MIER3, and SLC24A4). These adaptive alleles likely contribute to the biological and cellular functions determining phenotype in the Simbra hybrid cattle breed. Our data suggested that these alleles were introgressed from the breed's original indicine and taurine ancestors. The Simbra breed thus possesses derived parental alleles that combine the superior traits of the founder Brahman and Simmental breeds. These regions and genes might represent good targets for ad-hoc physiological studies, selection of breeding material and eventually even gene editing, for improved traits in modern cattle breeds. This study represents an important step toward developing and improving strategies for selection and population breeding to ultimately contribute meaningfully to the beef production industry.

Whole-genome sequencing reveals insights into the adaptation of French Charolais cattle to Cuban tropical conditions

Background

In the early 20th century, Cuban farmers imported Charolais cattle (CHFR) directly from France. These animals are now known as Chacuba (CHCU) and have become adapted to the rough environmental tropical conditions in Cuba. These conditions include long periods of drought and food shortage with extreme temperatures that European taurine cattle have difficulty coping with.

Results

In this study, we used whole-genome sequence data from 12 CHCU individuals together with 60 whole-genome sequences from six additional taurine, indicus and crossed breeds to estimate the genetic diversity, structure and accurate ancestral origin of the CHCU animals. Although CHCU animals are assumed to form a closed population, the results of our admixture analysis indicate a limited introgression of Bos indicus. We used the extended haplotype homozygosity (EHH) approach to identify regions in the genome that may have had an important role in the adaptation of CHCU to tropical conditions. Putative selection events occurred in genomic regions with a high proportion of Bos indicus, but they were not sufficient to explain adaptation of CHCU to tropical conditions by Bos indicus introgression only. EHH suggested signals of potential adaptation in genomic windows that include genes of taurine origin involved in thermogenesis (ATP9A, GABBR1, PGR, PTPN1 and UCP1) and hair development (CCHCR1 and CDSN). Within these genes, we identified single nucleotide polymorphisms (SNPs) that may have a functional impact and contribute to some of the observed phenotypic differences between CHCU and CHFR animals.

Conclusions

Whole-genome data confirm that CHCU cattle are closely related to Charolais from France (CHFR) and Canada, but also reveal a limited introgression of Bos indicus genes in CHCU. We observed possible signals of recent adaptation to tropical conditions between CHCU and CHFR founder populations, which were largely independent of the Bos indicus introgression. Finally, we report candidate genes and variants that may have a functional impact and explain some of the phenotypic differences observed between CHCU and CHFR cattle.

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.