On the History of Cattle Genetic Resources

Identification of miR-29a as a novel biomarker for lumpy skin disease virus exposure in cattle

Micro RNAs (miRNAs) have been implicated in the regulation of maturation, proliferation, differentiation, and activation of immune cells. In this study, we demonstrated that miR-29a antagonizes IFN-γ production at early times post-LSDV infection in cattle. miR-29a was predicted to target upstream IFN-γ regulators, and its inhibition resulted in enhanced IFN-γ production in sensitized peripheral blood mononuclear cells (PBMCs). Further, stimulation of PBMCs with LSDV antigen exhibited lower levels of miR-29a, concomitant with a potent cell-mediated immune response (CMI), characterized by an increase in LSDV-specific CD8+ T cell counts and enhanced levels of IFN-γ, which eventually facilitated virus clearance. In addition, a few immunocompromised cattle (developed secondary LSDV infection at ~ 6 months) that failed to mount a potent cell-mediated immune response, were shown to maintain higher miR-29a levels. Furthermore, as compared to the sensitized crossbred cattle, PBMCs from sensitized Rathi (a native Indian breed) animals exhibited lower levels of miR-29a along with an increase in CD8+ T cell counts and enhanced levels of IFN-γ. Finally, we analysed that a ≥ 60% decrease in miR-29a expression levels in the PBMCs of sensitized cattle correlated with a potent CMI response. In conclusion, miR-29a expression is involved in antagonizing the IFN-γ response in LSDV-infected cattle and may serve as a novel biomarker for the acute phase of LSDV infection, as well as predicting the functionality of T cells in sensitized cattle. In addition, Rathi cattle mount a more potent CMI response against LSDV than crossbred cattle.

Whole-genome resequencing uncovers diversity and selective sweep in Kazakh cattle

The Kazakh cattle in the Xinjiang Uygur Autonomous Region of China are highly adaptable and have multiple uses, including milk and meat production, and use as draft animals. They are an excellent original breed that could be enhanced by breeding and hybrid improvement. However, the genomic diversity and signature of selection underlying the germplasm characteristics require further elucidation. Herein, we evaluated 26 Kazakh cattle genomes in comparison with 103 genomes of seven other cattle breeds from regions around the world to assess the Kazakh cattle genetic variability. We revealed that the relatively low linkage disequilibrium at large SNP distances was strongly correlated with the largest effective population size among Kazakh cattle. Using population structural analysis, we next demonstrated a taurine lineage with restricted Bos indicus introgression among Kazakh cattle. Notably, we identified putative selected genes associated with resistance to disease and body size within Kazakh cattle. Together, our findings shed light on the evolutionary history and breeding profile of Kazakh cattle, as well as offering indispensable resources for germplasm resource conservation and crossbreeding program implementation.

Whole-genome resequencing reveals diversity and selective signals in the Wuxue goat

Animal genetic resources are crucial for ensuring global food security. However, in recent years, a noticeable decline in the genetic diversity of livestock has occurred worldwide. This decline is pronounced in developing countries, where the management of these resources is insufficient. In the current study, we performed whole genome sequencing for 20 Wuxue (WX) and five Guizhou White (GW) goats. Additionally, we utilized the published genomes of 131 samples representing five different goat breeds from various regions in China. We investigated and compared the genetic diversity and selection signatures of WX goats. Whole genome sequencing analysis of the WX and GW populations yielded 120 425 063 SNPs, which resided primarily in intergenic and intron regions. Population genetic structure revealed that WX exhibited genetic resemblance to GW, Chengdu Brown, and Jintang Black and significant differentiation from the other goat breeds. In addition, three methods (nucleotide diversity, linkage disequilibrium decay, and runs of homozygosity) showed moderate genetic diversity in WX goats. We used nucleotide diversity and composite likelihood ratio methods to identify within-breed signatures of positive selection in WX goats. A total of 369 genes were identified using both detection methods, including genes related to reproduction (GRID2, ZNF276, TCF25, and SPIRE2), growth (HMGA2 and GJA3), and immunity (IRF3 and SRSF3). Overall, this study explored the adaptability of WX goats, shedding light on their genetic richness and potential to thrive in challenges posed by climatic changes and diseases. Further investigations are warranted to harness these insights to enhance more efficient and sustainable goat breeding initiatives.

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.

Candidate signatures of positive selection for environmental adaptation in indigenous African cattle: A review

Environmental adaptation traits of indigenous African cattle are increasingly being investigated to respond to the need for sustainable livestock production in the context of unpredictable climatic changes. Several studies have highlighted genomic regions under positive selection probably associated with adaptation to environmental challenges (e.g. heat stress, trypanosomiasis, tick and tick-borne diseases). However, little attention has focused on pinpointing the candidate causative variant(s) controlling the traits. This review compiled information from 22 studies on signatures of positive selection in indigenous African cattle breeds to identify regions under positive selection. We highlight some key candidate genome regions and genes of relevance to the challenges of living in extreme environments (high temperature, high altitude, high infectious disease prevalence). They include candidate genes involved in biological pathways relating to innate and adaptive immunity (e.g. BoLAs, SPAG11, IL1RL2 and GFI1B), heat stress (e.g. HSPs, SOD1 and PRLH) and hypoxia responses (e.g. BDNF and INPP4A). Notably, the highest numbers of candidate regions are found on BTA3, BTA5 and BTA7. They overlap with genes playing roles in several biological functions and pathways. These include but are not limited to growth and feed intake, cell stability, protein stability and sweat gland development. This review may further guide targeted genome studies aiming to assess the importance of candidate causative mutations, within regulatory and protein-coding genome regions, to further understand the biological mechanisms underlying African cattle's unique adaption.

Global genetic diversity, introgression, and evolutionary adaptation of indicine cattle revealed by whole genome sequencing

Indicine cattle, also referred to as zebu (Bos taurus indicus), play a central role in pastoral communities across a wide range of agro-ecosystems, from extremely hot semiarid regions to hot humid tropical regions. However, their adaptive genetic changes following their dispersal into East Asia from the Indian subcontinent have remained poorly documented. Here, we characterize their global genetic diversity using high-quality whole-genome sequencing data from 354 indicine cattle of 57 breeds/populations, including major indicine phylogeographic groups worldwide. We reveal their probable migration into East Asia was along a coastal route rather than inland routes and we detected introgression from other bovine species. Genomic regions carrying morphology-, immune-, and heat-tolerance-related genes underwent divergent selection according to Asian agro-ecologies. We identify distinct sets of loci that contain promising candidate variants for adaptation to hot semi-arid and hot humid tropical ecosystems. Our results indicate that the rapid and successful adaptation of East Asian indicine cattle to hot humid environments was promoted by localized introgression from banteng and/or gaur. Our findings provide insights into the history and environmental adaptation of indicine cattle.

The History of Bovine Genital Campylobacteriosis in the Face of Political Turmoil and Structural Change in Cattle Farming in Germany

Contagious bovine genital campylobacteriosis (BGC), also known as bovine venereal campylobacteriosis, is a disease relevant to international trade listed by the World Organization for Animal Health (WOAH). It is caused by Campylobacter fetus subsp. venerealis (Cfv), one of three subspecies of Campylobacter fetus. Bulls are the reservoir but BGC may also be spread by artificial insemination (AI). BGC is characterized by severe reproductive losses such as infertility, early embryonic death and abortion with considerable economic losses. This significant economic impact has prompted several countries to adopt stringent eradication and surveillance measures to contain the disease. While there are commercial and autologous vaccines available, scientific evidence for the effectiveness of vaccination is still lacking. In Germany, BCG was already found to be endemic in the 1920s, shortly after the agent and the disease had been described for the first time. It can be assumed that BCG had already circulated uncontrolled for a long time in the predecessor states of Germany, influenced only by the political situation and trading networks of the time. After WW II, BCG was eradicated in the German Democratic Republic due to industrialized cattle production based on AI but it was still endemic at low levels in the Federal Republic of Germany with its diverse cattle production. There has been a steady decline in BGC incidence in re-unified Germany over the past 28 years. A single genetic Cfv lineage was identified which probably emerged in the 19th century and diversified over time. Interestingly, no recurrent cross-border introduction became evident. This review gives insight into the history of bovine genital campylobacteriosis considering the structural change in cattle farming in Germany and reflecting on the political background of the time.

Evaluation of breeding objectives, breeding practices and reproductive performance of indigenous dairy cows in selected districts of Kaffa Zone, South West Ethiopia

BACKGROUND

Breeding objectives are designed to achieve targeted dairy cow production goals, which can be affected by production type, farmer preferences, environmental factors and genetic factors individually or in combination. Breeding practices, such as both controlled and uncontrolled, and artificial insemination (AI) are the tools used to obtain the desired breeding objectives. The lower reproductive performance of indigenous dairy cows affects the total milk production and calf crops that are produced during their lifetime. Designing appropriate breeding objectives and breeding practices can improve the reproductive performance of dairy cows and their overall production performance.

MATERIALS AND METHODS

The current study was conducted with the objective of evaluating the breeding, practices and performance of indigenous dairy cattle in the south western part of Ethiopia. The districts of Gesha and Chena were purposefully chosen. The study design for the 384 household surveys was a cross-sectional survey with a simple random sample approach. Data analysis was carried out by MS-Excel (2010) and the general linear model procedure of SAS of 2008.

RESULTS

The current study revealed that methods of breeding were predominantly natural-controlled mating, followed by natural-uncontrolled mating and AI in descending order. Breeding objectives were input function, output function, sociocultural and economic functions and assets and security functions in decreasing order of rank. Reproduction performance indexes of indigenous dairy cows age at first service (3.72 ± 0.05 years), age at first calving (AFC) (4.71 ± 0.07 years), calving interval (CI) (1.58 ± 0.03 years), days open (DO) (4.26 ± 0.11 months), services per conception in natural mating (1.4 ± 0.08) and AI (2.73 ± 0.14), age of bull at maturity (4.17 ± 0.74 years), interoestrus interval (23.18 ± 0.61 days), calves crop (7.53 ± 0.22) and the life span of indigenous dairy cow (11.94 ± 0.26 years) were significant (p < 0.01) between two districts, whereas the values of age of bull at maturity and number of services per conception in natural mating were significant (p < 0.05) between districts.

CONCLUSIONS

Using AI and major reproduction performances, such as AFC, CI and DO of indigenous dairy cows in the study area, were very low. Therefore, concerned bodies should intervene to improve reproduction performance through the utilization of AI techniques, with the integration of forage development activities and improvements in livestock health care.

Exploring the microbiome of two uterine sites in cows

Bacterial communities in the mammalian reproductive system can be rich and diverse, differing in structure and quantity depending on location. In addition, its microbiome is associated with the state of health of this tract and reproductive success. This study evaluated the microbiome composition of the uterine body (UB) and uterine horn mucosa (UH) samples using 16S rRNA sequencing of samples extracted from cows in the Amazon region. It was observed that four main phyla were shared between the uterine sites: Actinobacteria, Bacteroidetes, Firmicutes, and Proteobacteria. Linear discriminant analysis effect size and heat tree analysis showed that members of Lachnospiraceae (NK3A20 group) and Oscillospiraceae were significantly more abundant in the UB than in UH. In addition, there are more unique genera in the UB than in the UH. A higher bacterial load in UB than in UH is expected because of the exposure to external factors of UB. However, comparing the site's communities through beta diversity did not generate well-defined clustering. Thus, it can be attributed to the closeness of the sites, which would make the niches similar ecologically and microbiologically. Therefore, this research provides knowledge to understand biomarkers in the prior reproduction period.

A Chinese indicine pangenome reveals a wealth of novel structural variants introgressed from other Bos species

Chinese indicine cattle harbor a much higher genetic diversity compared with other domestic cattle, but their genome architecture remains uninvestigated. Using PacBio HiFi sequencing data from 10 Chinese indicine cattle across southern China, we assembled 20 high-quality partially phased genomes and integrated them into a multiassembly graph containing 148.5 Mb (5.6%) of novel sequence. We identified 156,009 high-confidence nonredundant structural variants (SVs) and 206 SV hotspots spanning ∼195 Mb of gene-rich sequence. We detected 34,249 archaic introgressed fragments in Chinese indicine cattle covering 1.93 Gb (73.3%) of the genome. We inferred an average of 3.8%, 3.2%, 1.4%, and 0.5% of introgressed sequence originating, respectively, from banteng-like, kouprey-like, gayal-like, and gaur-like Bos species, as well as 0.6% of unknown origin. Introgression from multiple donors might have contributed to the genetic diversity of Chinese indicine cattle. Altogether, this study highlights the contribution of interspecies introgression to the genomic architecture of an important livestock population and shows how exotic genomic elements can contribute to the genetic variation available for selection.

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.

Effects of local domestication warrant attention in honey bee population genetics

Honey bees, Apis mellifera, have for millennia been managed and exploited by humans and introduced into most suitable regions worldwide. However, given the lack of records for many introduction events, treating A. mellifera populations as native would predictably bias genetic studies regarding origin and evolution. Here, we used the Dongbei bee, a well-documented population, introduced beyond the natural distribution range approximately 100 years ago, to elucidate the effects of local domestication on animal population genetic analyses. Strong domestication pressure was detected in this population, and the genetic divergence between Dongbei bee and its ancestral subspecies was found to have occurred at the lineage level. Results of phylogenetic and time divergence analyses could consequently be misinterpreted. Proposing new subspecies or lineages and performing analyses of origin should thus strive to eliminate anthropogenic effects. We highlight the need for definitions of landrace and breed in honey bee sciences and make preliminary suggestions.

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.

Pedigree and genome-based patterns of homozygosity in the South African Ayrshire, Holstein, and Jersey breeds

The erosion of genetic diversity limits long-term genetic gain and impedes the sustainability of livestock production. In the South African (SA) dairy industry, the major commercial dairy breeds have been applying estimated breeding values (EBVs) and/or have been participating in Multiple Across Country Evaluations (MACE). The transition to genomic estimated breeding values (GEBVs) in selection strategies requires monitoring of the genetic diversity and inbreeding of current genotyped animals, especially considering the comparatively small population sizes of global dairy breeds in SA. This study aimed to perform a homozygosity-based evaluation of the SA Ayrshire (AYR), Holstein (HST), and Jersey (JER) dairy cattle breeds. Three sources of information, namely 1) single nucleotide polymorphism (SNP) genotypes (3,199 animals genotyped for 35,572 SNPs) 2) pedigree records (7,885 AYR; 28,391 HST; 18,755 JER), and 3) identified runs of homozygosity (ROH) segments were used to quantify inbreeding related parameters. The lowest pedigree completeness was for the HST population reducing from a value of 0.990 to 0.186 for generation depths of one to six. Across all breeds, 46.7% of the detected ROH were between 4 megabase pairs (Mb) and 8 Mb in length. Two conserved homozygous haplotypes were identified in more than 70% of the JER population on Bos taurus autosome (BTA) 7. The JER breed displayed the highest level of inbreeding across all inbreeding coefficients. The mean (± standard deviation) pedigree-based inbreeding coefficient (FPED) ranged from 0.051 (±0.020) for AYR to 0.062 (±0.027) for JER, whereas SNP-based inbreeding coefficients (FSNP) ranged from 0.020 (HST) to 0.190 (JER) and ROH-based inbreeding coefficients, considering all ROH segment coverage (FROH), ranged from 0.053 (AYR) to 0.085 (JER). Within-breed Spearman correlations between pedigree-based and genome-based estimates ranged from weak (AYR: 0.132 between FPED and FROH calculated for ROH &lt;4Mb in size) to moderate (HST: 0.584 between FPED and FSNP). Correlations strengthened between FPED and FROH as the ROH length category was considered lengthened, suggesting a dependency on breed-specific pedigree depth. The genomic homozygosity-based parameters studied proved useful in investigating the current inbreeding status of reference populations genotyped to implement genomic selection in the three most prominent South African dairy cattle breeds.

Genomic diversity and relationship analyses of endangered German Black Pied cattle (DSN) to 68 other taurine breeds based on whole-genome sequencing

German Black Pied cattle (Deutsches Schwarzbuntes Niederungsrind, DSN) are an endangered dual-purpose cattle breed originating from the North Sea region. The population comprises about 2,500 cattle and is considered one of the ancestral populations of the modern Holstein breed. The current study aimed at defining the breeds closest related to DSN cattle, characterizing their genomic diversity and inbreeding. In addition, the detection of selection signatures between DSN and Holstein was a goal. Relationship analyses using fixation index (F_ST), phylogenetic, and admixture analyses were performed between DSN and 68 other breeds from the 1000 Bull Genomes Project. Nucleotide diversity, observed heterozygosity, and expected heterozygosity were calculated as metrics for genomic diversity. Inbreeding was measured as excess of homozygosity (F_Hom) and genomic inbreeding (F_RoH) through runs of homozygosity (RoHs). Region-wide F_ST and cross-population-extended haplotype homozygosity (XP-EHH) between DSN and Holstein were used to detect selection signatures between the two breeds, and RoH islands were used to detect selection signatures within DSN and Holstein. DSN showed a close genetic relationship with breeds from the Netherlands, Belgium, Northern Germany, and Scandinavia, such as Dutch Friesian Red, Dutch Improved Red, Belgian Red White Campine, Red White Dual Purpose, Modern Angler, Modern Danish Red, and Holstein. The nucleotide diversity in DSN (0.151%) was higher than in Holstein (0.147%) and other breeds, e.g., Norwegian Red (0.149%), Red White Dual Purpose (0.149%), Swedish Red (0.149%), Hereford (0.145%), Angus (0.143%), and Jersey (0.136%). The F_Hom and F_RoH values in DSN were among the lowest. Regions with high F_ST between DSN and Holstein, significant XP-EHH regions, and RoH islands detected in both breeds harbor candidate genes that were previously reported for milk, meat, fertility, production, and health traits, including one QTL detected in DSN for endoparasite infection resistance. The selection signatures between DSN and Holstein provide evidence of regions responsible for the dual-purpose properties of DSN and the milk type of Holstein. Despite the small population size, DSN has a high level of diversity and low inbreeding. F_ST supports its relatedness to breeds from the same geographic origin and provides information on potential gene pools that could be used to maintain diversity in DSN.

Copy Number Variants in Two Northernmost Cattle Breeds Are Related to Their Adaptive Phenotypes

Copy number variations (CNVs) are genomic structural variants with potential functional and evolutionary effects on phenotypes. In this study, we report the identification and characterization of CNVs from the whole-genome resequencing data of two northernmost cattle breeds from Russia: the Yakut and Kholmogory cattle and their phylogenetically most related breeds, Hanwoo and Holstein, respectively. Comparisons of the CNV regions (CNVRs) among the breeds led to the identification of breed-specific CNVRs shared by cold-adapted Kholmogory and Yakut cattle. An investigation of their overlap with genes, regulatory domains, conserved non-coding elements (CNEs), enhancers, and quantitative trait loci (QTLs) was performed to further explore breed-specific biology and adaptations. We found CNVRs enriched for gene ontology terms related to adaptation to environments in both the Kholmogory and Yakut breeds and related to thermoregulation specifically in Yakut cattle. Interestingly, the latter has also been supported when exploring the enrichment of breed-specific CNVRs in the regulatory domains and enhancers, CNEs, and QTLs implying the potential contribution of CNVR to the Yakut and Kholmogory cattle breeds’ adaptation to a harsh environment.

How Geography and Climate Shaped the Genomic Diversity of Italian Local Cattle and Sheep Breeds

Simple Summary

In this paper, we study the inter-relationships among geography, climate, and genetics in Italian local cattle and sheep breeds. In terms of genetic diversity, geography (latitude and longitude) appears to play a larger role in sheep (26.4%) than that in cattle (13.8%). Once geography is accounted for, 10.1% of cattle genomic diversity and 13.3% of that of sheep are attributable to climatic effects. Stronger geographic effects in sheep can be related to a combination of higher predomestication genetic variability together with biological and productive specializations. The climate alone seems to have had less impact on the current genetic diversity in both species even if climate and geography are greatly confounded. Results confirm that both species are the result of complex evolutionary histories triggered by interactions between human needs and environmental conditions.

Abstract

Understanding the relationships among geography, climate, and genetics is increasingly important for animal farming and breeding. In this study, we examine these inter-relationships in the context of local cattle and sheep breeds distributed along the Italian territory. To this aim, we used redundancy analysis on genomic data from previous projects combined with geographical coordinates and corresponding climatic data. The effect of geographic factors (latitude and longitude) was more important in sheep (26.4%) than that in cattle (13.8%). Once geography had been partialled out of analysis, 10.1% of cattle genomic diversity and 13.3% of that of sheep could be ascribed to climatic effects. Stronger geographic effects in sheep can be related to a combination of higher pre-domestication genetic variability together with biological and productive specificities. Climate alone seems to have had less impact on current genetic diversity in both species, even if climate and geography are greatly confounded. Results confirm that both species are the result of complex evolutionary histories triggered by interactions between human needs and environmental conditions.

Trajectory of livestock genomics in South Asia: A comprehensive review

Livestock plays a central role in sustaining human livelihood in South Asia. There are numerous and distinct livestock species in South Asian countries. Several of them have experienced genetic development in recent years due to the application of genomic technologies and effective breeding programs. This review discusses genomic studies on cattle, buffalo, sheep, goat, pig, horse, camel, yak, mithun, and poultry. The frontiers covered in this review are genetic diversity, admixture studies, selection signature research, QTL discovery, genome-wide association studies (GWAS), and genomic selection. The review concludes with recommendations for South Asian livestock systems to increasingly leverage genomic technologies, based on the lessons learned from the numerous case studies. This paper aims to present a comprehensive analysis of the dichotomy in the South Asian livestock sector and argues that a realistic approach to genomics in livestock can ensure long-term genetic advancements.

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.

Assessment of genetic diversity and relationship of the two Sanga type cattle of Botswana based on microsatellite markers

The study was performed to evaluate genetic variation on two Sanga type cattle found in Botswana: Tswana and Tuli using twelve microsatellite markers. All amplified loci were polymorphic with 75 and 77 alleles genotyped in the Tswana and Tuli populations, respectively. The total number of alleles per locus ranged from 2 (BM1818) to 10 (TGLA227) with total mean of 6.25 for Tswana and 6.43 for Tuli population. Almost all the markers showed high polymorphic information content (PIC) apart from BM1818 (0.375) and INRA23 (0.393) which were moderately informative in Tswana population. Most of the markers were in Hardy-Weinberg equilibrium except for CSSRM60 and CSSM66 loci in Tswana population and ETH10, ETH225 and CSSM66 loci in Tuli population. A total of 103 unique alleles were genotyped across the two breeds with 49-shared, and 26 and 28 were unique to Tswana and Tuli populations, respectively. The expected heterozygosity (He) values were higher than the observed heterozygosity (Ho) in both populations: Tswana (He = 0.7895 ± 0.033 vs Ho = 0.631 ± 0.091) and Tuli (He = 0.8123 ± 0.033 vs Ho = 0.556 ± 0.021). The inbreeding coefficient was 0.200 ± 0.002 and 0.332 ± 0.001 in Tswana and Tuli populations, respectively. Analysis of molecular variance revealed 6.8% of the total genetic variation corresponding to differences between the two breeds and 93.2% within populations. The genetic identity between the two breeds was 56% and there were similar levels of multilocus heterozygosity and allelic diversity in the two breeds. The use of Tswana and Tuli breeds in a crossbreeding programme is likely to result in minimal heterosis and therefore not recommended.

The complete mitochondrial genome of an ancient cattle (Bos taurus) from Taosi site, China, and its phylogenetic assessment

The complete mitochondrial genome (mitogenome) of one 4000-years-old cattle from Taosi site was determined by high throughput sequencing. The mitogenome was 16,336 bp in length and contained 13 protein-coding genes, two rRNA genes, and 22 tRNA genes. The protein-coding genes had two types of start codons (ATG and ATA) and three types of stop codons (TAA, TAG, and AGA). The overall base composition of the genome was 33%-A, 27%-T, 26%-C, 14%-G. The matrilineal genealogical analysis based on mitogenome revealed that the 4000-years-old cattle from Taosi site was domestic taurine cattle. In this study, we not only reported a complete mitogenome for a 4000-years-old bovine remain from the middle Yellow River region but also provided the mitogenomic evidence for the close phylogenetic relationship between the early taurine cattle in Northern China and modern domestic cattle.

Why and How European Farmers Are Dedicated to Breeding the Dwarf Dahomey Cattle

Simple Summary

The introduction of high productive cattle breeds in Africa is well known, but the contribution of African breeds to livestock biodiversity in Europe is generally overlooked. This study reports, for the first time, European farmers’ interests in keeping the Dahomey cattle, and characterizes their management practices. The Dahomey cattle from Benin (West Africa) are the smallest cattle in the world, and they were introduced to Europe in the early 1900s. The findings revealed that European farmers are increasingly interested in keeping Dahomey cattle, because of their suitability for grassland maintenance and meat production, as well as their low management requirements (with regard to feeding, preventive and curative health care and reproduction management). Overall, the study displays the agricultural importance and ecological utilization of Dahomey cattle in European countries. It shows how small-sized cattle can support the promotion of sustainable livestock production and the management of ecosystems, including faunistic and floristic diversity.

Abstract

This study investigates the motivations and breeding practices of farmers keeping Dahomey cattle in European countries. Data were collected using a web-based open-closed questionnaire survey targeting 55 farmers from Germany, Switzerland and Austria. Descriptive analyses revealed that the earliest European Dahomey herds were established in 2005. Moreover, interest in the breed recently increased as 63.7% of the investigated farmers established their herds between 2016 and 2020. The average herd size comprises seven Dahomey cattle, kept for managing grassland (59.3%), for production of meat or as breeding stock (32.1%) and for a hobby (8.6%). The animals are mostly kept in grazing systems throughout the year, partly fattened with supplement feeds. The low disease incidence and no need for extra health care in the herds indicate the robustness of the breed. Furthermore, meat quality, calving ease, small size, calm character and low feed requirements of Dahomey cattle were valued by the farmers. For the preservation of these features, farmers confirmed their enthusiasm to support any breeding and conservation program of this smallholder breed in Europe and Benin. This study highlights the importance of small-sized cattle for sustainable breeding systems and with regard to ecosystem management practices.

Single Nucleotide Polymorphisms in the Bovine TLR2 Extracellular Domain Contribute to Breed and Species-Specific Innate Immune Functionality

Recent evidence suggests that several cattle breeds may be more resistant to infection with the zoonotic pathogen Mycobacterium bovis. Our data presented here suggests that the response to mycobacterial antigens varies in macrophages generated from Brown Swiss (BS) and Holstein Friesian (HF) cattle, two breeds belonging to the Bos taurus family. Whole genome sequencing of the Brown Swiss genome identified several potential candidate genes, in particular Toll-like Receptor-2 (TLR2), a pattern recognition receptor (PRR) that has previously been described to be involved in mycobacterial recognition. Further investigation revealed single nucleotide polymorphisms (SNP) in TLR2 that were identified between DNA isolated from cells of BS and HF cows. Interestingly, one specific SNP, H326Q, showed a different genotype frequency in two cattle subspecies, Bos (B.) taurus and Bos indicus. Cloning of the TLR2 gene and subsequent gene-reporter and chemokine assays revealed that this SNP, present in BS and Bos indicus breeds, resulted in a significantly higher response to mycobacterial antigens as well as tri-acylated lipopeptide ligands in general. Comparing wild-type and H326Q containing TLR2 responses, wild-type bovine TLR2 response showed clear, diminished mycobacterial antigen responses compared to human TLR2, however bovine TLR2 responses containing H326Q were found to be partially recovered compared to human TLR2. The creation of human:bovine TLR2 chimeras increased the response to mycobacterial antigens compared to the full-length bovine TLR2, but significantly reduced the response compared to the full-length human TLR2. Thus, our data, not only present evidence that TLR2 is a major PRR in the mammalian species-specific response to mycobacterial antigens, but furthermore, that there are clear differences between the response seen in different cattle breeds, which may contribute to their enhanced or reduced susceptibility to mycobacterial infection.

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.

Signatures of selection are present in the genome of two close autochthonous cattle breeds raised in the North of Italy and mainly distinguished for their coat colours

Autochthonous cattle breeds are genetic resources that, in many cases, have been fixed for inheritable exterior phenotypes useful to understand the genetic mechanisms affecting these breed-specific traits. Reggiana and Modenese are two closely related autochthonous cattle breeds mainly raised in the production area of the well-known Protected Designation of Origin Parmigiano-Reggiano cheese, in the North of Italy. These breeds can be mainly distinguished for their standard coat colour: solid red in Reggiana and solid white with pale shades of grey in Modenese. In this study we genotyped with the GeneSeek GGP Bovine 150k single nucleotide polymorphism (SNP) chip almost half of the extant cattle populations of Reggiana (n = 1109 and Modenese (n = 326) and used genome-wide information in comparative F_ST analyses to detect signatures of selection that diverge between these two autochthonous breeds. The two breeds could be clearly distinguished using multidimensional scaling plots and admixture analysis. Considering the top 0.0005% F_ST values, a total of 64 markers were detected in the single-marker analysis. The top F_ST value was detected for the melanocortin 1 receptor (MC1R) gene mutation, which determines the red coat colour of the Reggiana breed. Another coat colour gene, agouti signalling protein (ASIP), emerged amongst this list of top SNPs. These results were also confirmed with the window-based analyses, which included 0.5-Mb or 1-Mb genome regions. As variability affecting ASIP has been associated with white coat colour in sheep and goats, these results highlighted this gene as a strong candidate affecting coat colour in Modenese breed. This study demonstrates how population genomic approaches designed to take advantage from the diversity between local genetic resources could provide interesting hints to explain exterior traits not yet completely investigated in cattle.

Zoometric Characterization of Creole Cows from the Southern Amazon Region of Peru

Bovine breeding began in Peru more than five centuries ago; since then, adaptation has started. Nowadays, Creole cattle are part of economic activities of Peruvian farmers, as they supply food and profits, among other things, for families in Andean region. Creole cattle have many strengths such as resistance to diseases, low nutrient requirements and easy adaptation to different levels of altitude and environmental conditions. However, even with all these attributes, they are not valued within production systems. Moreover, these valuable genetic resources are disappearing gradually in Perú. For this reason, this study aims to characterize the zoometry of Creole cows from the Southern Amazonas region of Peru. Biotype classification was performed by hierarchical cluster and multivariate factors analysis. Qualitative traits were analyzed with descriptive statistics and Duncan’s test (α = 0.05) was used to compare means among groups. Associations between qualitative traits were determined with Chi-square. We distinguished three biotypes of Creole cows with characteristics for meat and milk production. This information is helpful for future conservation programs for Creole cattle.

Disentangling the Multidimensional Relationship between Livestock Breeds and Ecosystem Services

Simple Summary

Livestock breeds represent the diversity of livestock animals. They participate in the delivery of ecosystem services (ES), i.e., the benefits to humans provided by nature. In recent years, the contribution of livestock breeds to ES has received attention in livestock research. Additionally, there is increasing interest in integrating this knowledge into policies to make agriculture more sustainable. In this work, we elaborate on livestock breed characteristics that are key to the study of livestock breed contributions to ES. Thus, we explore the natural and human factors that have produced livestock breeds as ecologically and culturally mediated entities. In addition, we review the different roles of livestock breeds as biodiversity components. Finally, we examine how livestock breeds participate in livestock system heterogeneity. By integrating these aspects, we might better understand how livestock breeds provide and modulate ES provision and, therefore, how to improve breed conservation and livestock policies toward more sustainable farming.

Abstract

There is an increasing interest in assessing livestock breed contributions to ecosystem services (ES) and including this knowledge in decision making. However, this task has been limited due to the complexity of the multidimensional relationship between livestock diversity and ecosystem services. In this work, we elaborate on the livestock breed characteristics central to developing a comprehensive approach to livestock breed inclusion in the ecosystem services framework. Thus, we explore the multidimensional nature of livestock breeds, i.e., as eco-cultural entities, biodiversity components, and drivers of livestock system heterogeneity and functioning. First, anthropogenic and natural factors have acted jointly to develop breeds as eco-cultural entities. This fact represents an opportunity to move toward farming system sustainability by Nature-Based Solutions and Nature’s Contribution to People paradigms. Second, livestock breeds are components of biodiversity, and as such, can be framed as goods, as final ecosystem services, and as regulators of ecosystem processes. Third, livestock breeds contribute to livestock system heterogeneity and resilience. By integrating these aspects, we might better understand how livestock breeds provide and modulate ecosystem service provision and, therefore, how to improve breed conservation and livestock policies toward farming system sustainability.

Demographic History, Adaptation, and NRAP Convergent Evolution at Amino Acid Residue 100 in the World Northernmost Cattle from Siberia

Native cattle breeds represent an important cultural heritage. They are a reservoir of genetic variation useful for properly responding to agriculture needs in the light of ongoing climate changes. Evolutionary processes that occur in response to extreme environmental conditions could also be better understood using adapted local populations. Herein, different evolutionary histories of the world northernmost native cattle breeds from Russia were investigated. They highlighted Kholmogory as a typical taurine cattle, whereas Yakut cattle separated from European taurines approximately 5,000 years ago and contain numerous ancestral and some novel genetic variants allowing their adaptation to harsh conditions of living above the Polar Circle. Scans for selection signatures pointed to several common gene pathways related to adaptation to harsh climates in both breeds. But genes affected by selection from these pathways were mostly different. A Yakut cattle breed-specific missense mutation in a highly conserved NRAP gene represents a unique example of a young amino acid residue convergent change shared with at least 16 species of hibernating/cold-adapted mammals from six distinct phylogenetic orders. This suggests a convergent evolution event along the mammalian phylogenetic tree and fast fixation in a single isolated cattle population exposed to a harsh climate.

Origin and population structure of Brazilian hair sheep breeds

Brazilian hair sheep constitute a genetic diversity hotspot. These animals are found in the harsh environments of the Brazilian Northwest (semi-arid) region. Genotypes (50K SNP chip) from seven Brazilian sheep breeds (five hair and two coarse wool types) and 87 worldwide breeds were used to test for population structure, admixture and genetic diversity. Moreover, phylogenetic trees evaluating migration events between genetic groups were built. Brazilian Somali, a fat-tailed breed, had a close relationship with East African breeds and clustered distinctly from other Brazilian breeds. Brazilian Blackbelly and Barbados Blackbelly had a close relationship. The Morada Nova breed did not show close relationships with European or African breeds, revealing a single migration event from an Algerian hair breed. Brazilian Fat-tail and Morada Nova share a common ancestor, but the former showed introgressions from Brazilian Somali and Afrikaner breeds, explaining the fat-tail phenotype. The Santa Inês breed received a substantial contribution from Brazilian Bergamasca and showed an admixed origin with recent introgressions from other breeds, mainly from Suffolk. Furthermore, Brazilian Somali and Brazilian Fat-tail are the most endangered sheep genetic resources in Brazil and should be the focus for ex situ conservation programs. In conclusion, Brazilian hair sheep show an African origin and are characterized by diverse genetic composition, reinforcing the need for conservation of these genetic resources, and at the same time, this highly diverse group has variability that can be used in breeding programs.

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.

Farm Animals Are Long Away from Natural Behavior: Open Questions and Operative Consequences on Animal Welfare

Simple Summary

Animal welfare is a very important issue. One of the tasks of researchers is to provide explanations and possible solutions to questions arising from non-experts. This work analyzes part of the extensive literature on relationships between selection and domestic, mainly farm, animals’ behavior and deals with some very important themes, such as the role of regulations, domestication, and selection.

Abstract

The concept of welfare applied to farm animals has undergone a remarkable evolution. The growing awareness of citizens pushes farmers to guarantee the highest possible level of welfare to their animals. New perspectives could be opened for animal welfare reasoning around the concept of domestic, especially farm, animals as partial human artifacts. Therefore, it is important to understand how much a particular behavior of a farm animal is far from the natural one of its ancestors. This paper is a contribution to better understand the role of genetics of the farm animals on their behavior. This means that the naïve approach to animal welfare regarding returning animals to their natural state should be challenged and that welfare assessment should be considered.

A low-density SNP genotyping panel for the accurate prediction of cattle breeds

Genomic tools to better define breed composition in agriculturally important species have sparked scientific and commercial industry interest. Knowledge of breed composition can inform multiple scientifically important decisions of industry application including DNA marker-assisted selection, identification of signatures of selection, and inference of product provenance to improve supply chain integrity. Genomic tools are expensive but can be economized by deploying a relatively small number of highly informative single-nucleotide polymorphisms (SNP) scattered evenly across the genome. Using resources from the 1000 Bull Genomes Project we established calibration (more stringent quality criteria; N = 1,243 cattle) and validation (less stringent; N = 864) data sets representing 17 breeds derived from both taurine and indicine bovine subspecies. Fifteen successively smaller panels (from 500,000 to 50 SNP) were built from those SNP in the calibration data that increasingly satisfied 2 criteria, high differential allele frequencies across the breeds as measured by average Euclidean distance (AED) and high uniformity (even spacing) across the physical genome. Those SNP awarded the highest AED were in or near genes previously identified as important signatures of selection in cattle such as LCORL, NCAPG, KITLG, and PLAG1. For each panel, the genomic breed composition (GBC) of each animal in the validation dataset was estimated using a linear regression model. A systematic exploration of the predictive accuracy of the various sized panels was then undertaken on the validation population using 3 benchmarking approaches: (1) % error (expressed relative to the estimated GBC made from over 1 million SNP), (2) % breed misassignment (expressed relative to each individual's breed recorded), and (3) Shannon's entropy of estimated GBC across the 17 target breeds. Our analyses suggest that a panel of just 250 SNP represents an adequate balance between accuracy and cost-only modest gains in accuracy are made as one increases panel density beyond this point.

Spatial distance between sites of sampling associated with genetic variation among Neospora caninum in aborted bovine foetuses from northern Italy

BACKGROUND

Neospora caninum, a coccidian protozoan, represents an important cause of bovine abortion. Available N. caninum strains show considerable variation in vitro and in vivo, including different virulence in cattle. To which extent sexual recombination, which is possible in the intestines of domestic dogs and closely related carnivores as definitive hosts, contributes to this variation is not clear yet.

METHODS

Aborted bovine foetuses were collected between 2015 and early 2019 from Italian Holstein Friesian dairy herds suffering from reproductive problems. A total of 198 samples were collected from 165 intensive farms located in Lombardy, northern Italy. N. caninum samples were subjected to multilocus-microsatellite genotyping using ten previously established microsatellite markers. In addition to our own data, those from a recent study providing data on five markers from other northern Italian regions were included and analysed.

RESULTS

Of the 55 samples finally subjected to genotyping, 35 were typed at all or 9 out of 10 loci and their individual multilocus-microsatellite genotype (MLMG) determined. Linear regression revealed a statistically significant association between the spatial distance of the sampling sites with the genetic distance of N. caninum MLMGs (P < 0.001). Including data from this and a previous North Italian study into eBURST analysis revealed that several of N. caninum MLMGs from northern Italy separate into four groups; most of the samples from Lombardy clustered in one of these groups. Principle component analysis revealed similar clusters and confirmed MLMG groups identified by eBURST. Variations observed between MLMGs were not equally distributed over all loci, but predominantly observed in MS7, MS6A, or MS10.

CONCLUSIONS

Our findings confirm the concept of local N. caninum subpopulations. The geographic distance of sampling was associated with the genetic distance as determined by microsatellite typing. Results suggest that multi-parental recombination in N. caninum is a rare event, but does not exclude uniparental mating. More comprehensive studies on microsatellites in N. caninum and related species like Toxoplasma gondii should be undertaken, not only to improve genotyping capabilities, but also to understand possible functions of these regions in the genomes of these parasites.

Recent land use and management changes decouple the adaptation of livestock diversity to the environment

Native livestock breeds, i.e. those autochthonous to a specific region, are locally adapted domesticated animals that conserve genetic resources, guaranty food security and provide agroecosystem services. Native breeds are largely threatened worldwide by agricultural intensification and rural areas abandonment processes related to recent changes in production schemes and planning. Yet, our gap of knowledge regarding livestock breed-environment relationships may prevent the design of successful conservation measures. In this work, we analyse the links between livestock diversity -i.e. richness of native breeds- and a selection of environmental factors that express at broad scales, with a temporal perspective. We compare native breeds distributional patterns before and after the agricultural intensification, in the context of land-use change in mainland Spain. Our results confirm the existence of strong associations between the distribution of native livestock breeds and environmental factors. These links, however, weaken for contemporary distributions. In fact, changes in breed distribution reflect a shift towards more productive environments. Finally, we found that the areas having higher breed richness are undergoing land abandonment processes. Succeeding in the conservation of threatened native breeds will require going beyond merely genetic and production-oriented views. Ecological and sociocultural perspectives should also be accounted for as global change processes are determinant for livestock agrobiodiversity.

The Identification of Runs of Homozygosity Gives a Focus on the Genetic Diversity and Adaptation of the "Charolais de Cuba" Cattle

Simple Summary

The Charolais de Cuba cattle is a tropical adapted breed founded in Cuba around 120 years ago from Charolais French specimens. Nowadays, it is still a closed breed and remains as a small population. In this work, we analyzed the inbreeding and diversity patterns, as well as the population size, of this recent adapted breed via a run of homozygosity (ROH) analysis. We found that the genomic inbreeding levels are higher in the Charolais de Cuba breed compared to French and British Charolais populations. Nevertheless, we detected that the effective population size experienced a very similar decline during the last century in the three Charolais populations studied. Finally, a number of regions with exceptional patterns of long homozygosity were identified in this breed, and these could be related to processes of adaptation to tropical conditions.

Abstract

Inbreeding and effective population size (Ne) are fundamental indicators for the management and conservation of genetic diversity in populations. Genomic inbreeding gives accurate estimates of inbreeding, and the Ne determines the rate of the loss of genetic variation. The objective of this work was to study the distribution of runs of homozygosity (ROHs) in order to estimate genomic inbreeding (F_ROH) and an effective population size using 38,789 Single Nucleotide Polymorphisms (SNPs) from the Illumina Bovine 50K BeadChip in 86 samples from populations of Charolais de Cuba (n = 40) cattle and to compare this information with French (n = 20) and British Charolais (n = 26) populations. In the Cuban, French, and British Charolais populations, the average estimated genomic inbreeding values using the F_ROH statistics were 5.7%, 3.4%, and 4%, respectively. The dispersion measured by variation coefficient was high at 43.9%, 37.0%, and 54.2%, respectively. The effective population size experienced a very similar decline during the last century in Charolais de Cuba (from 139 to 23 individuals), in French Charolais (from 142 to 12), and in British Charolais (from 145 to 14) for the ~20 last generations. However, the high variability found in the ROH indicators and F_ROH reveals an opportunity for maintaining the genetic diversity of this breed with an adequate mating strategy, which can be favored with the use of molecular markers. Moreover, the detected ROH were compared to previous results obtained on the detection of signatures of selection in the same breed. Some of the observed signatures were confirmed by the ROHs, emphasizing the process of adaptation to tropical climate experienced by the Charolais de Cuba population.

An insight into the evolutionary history of Indonesian cattle assessed by whole genome data analysis

The domestication of Indonesian cattle was investigated through a study of their genetic diversity, up to the genome level. Little documentation exists regarding the history of domestication of Indonesian cattle and questions remain despite a growing body of molecular evidence. In this study, we genotyped seven Indonesian cattle breeds using an Illumina BovineSNP50 Bead Chip to provide insight into their domestication and demographic history in a worldwide population context. Our analyses indicated the presence of hybrid cattle, with Bos javanicus and Bos indicus ancestries being most prevalent, as well as purebred cattle. We revealed that all the breeds were interconnected through several migration events. However, their demographic status varied widely. Although almost all the Indonesian cattle had an effective population size higher than the minimum level required to ensure breed fitness, efforts are still needed to maintain their genetic variability and purity.

Polymorphism of the Melanocortin 1 Receptor (MC1R) Gene and its Role in Determining the Coat Colour of Central European Cattle Breeds

There are many genes responsible for the appearance of different coat colours, among which the melanocortin 1 receptor gene (MC1R) plays an important role. The aim of the study was to characterize genetic variation in Central European cattle breeds based on polymorphism of the MC1R gene and factors determining their coat colour. The study was conducted on 290 individuals of the following breeds: Polish White-Backed (PW), Lithuanian White-Backed (LW), Polish Red (PR), Lithuanian Red (LR), Carpathian Brown (CB), Ukrainian Grey (UG), and Slovak Pinzgau (SP). Polymorphism at the MC1R gene locus was analysed by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) using two restriction enzymes: Cfr10I and SsiI. The proportions of alleles and genotypes in the MC1R locus indicates a strong relationship between polymorphism and the coat colour of cattle: The ED allele proved to be characteristic for the breeds with a white-backed coat (PW and LW), while the dominant allele in the red breeds (PR and LR) was E+. It is noteworthy that coat colour in the SP population was determined only by the recessive e allele, which resulted in the formation of a separate clade in the phylogenetic tree.

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.

Diversity and Resilience to Socio-Ecological Changes of Smallholder Lagune Cattle Farming Systems of Benin

The indigenous smallholder Lagune breed and the production systems in which it is embedded in Southern Benin have received very little research or policy attention. Consequently, very little information exists on the diversity of these production systems or on their capacity to adapt to ongoing socio-economic and environmental changes. This study aimed to explore and characterize the diversity of Lagune cattle production systems along with farmers’ local knowledge and resilience strategies. A questionnaire was administered to 417 Lagune cattle farmers across two agro-ecological zones: Ouémé Valley (OVZ) and Pobe (PZ). It included, inter alia, questions related to households’ socio-economic conditions, their cattle herd characteristics, and management practices. Categorical principal component analysis and the two-step clustering method were used to classify the production systems which were then compared using the chi-square and ANOVA procedures. Four distinct farm types were identified. This study revealed the important role of agroecology in the diversity of farmers’ breeding practices. Controlled mating was more common in tethering systems whereas uncontrolled mating, widespread in free-roaming systems, has favored Lagune breed admixture with zebus. Opportunities for conserving the genetic diversity within the Lagune breed might be greater in PZ where breed admixture was almost inexistent.

Refining the genetic structure and relationships of European cattle breeds through meta-analysis of worldwide genomic SNP data, focusing on Italian cattle

The availability of genotyping assays has allowed the detailed evaluation of cattle genetic diversity worldwide. However, these comprehensive studies did not include some local European populations, including autochthonous Italian cattle. In this study, we assembled a large-scale, genome-wide dataset of single nucleotide polymorphisms scored in 3,283 individuals from 205 cattle populations worldwide to assess genome-wide autozygosity and understand better the genetic relationships among these populations. We prioritized European cattle, with a special focus on Italian breeds. Moderate differences in estimates of molecular inbreeding calculated from runs of homozygosity (FROH) were observed among domesticated bovid populations from different geographic areas, except for Bali cattle. Our findings indicated that some Italian breeds show the highest estimates of levels of molecular inbreeding among the cattle populations assessed in this study. Patterns of genetic differentiation, shared ancestry, and phylogenetic analysis all suggested the occurrence of gene flow, particularly among populations originating from the same geographical area. For European cattle, we observed a distribution along three main directions, reflecting the known history and formation of the analyzed breeds. The Italian breeds are split into two main groups, based on their historical origin and degree of conservation of ancestral genomic components. The results pinpointed that also Sicilian breeds, much alike Podolian derived-breeds, in the past experienced a similar non-European influence, with African and indicine introgression.

The Pastoral Origin of Semiotically Functional Tonal Organization of Music

This paper presents a new line of inquiry into when and how music as a semiotic system was born. Eleven principal expressive aspects of music each contains specific structural patterns whose configuration signifies a certain affective state. This distinguishes the tonal organization of music from the phonetic and prosodic organization of natural languages and animal communication. The question of music’s origin can therefore be answered by establishing the point in human history at which all eleven expressive aspects might have been abstracted from the instinct-driven primate calls and used to express human psycho-emotional states. Etic analysis of acoustic parameters is the prime means of cross-examination of the typical patterns of expression of the basic emotions in human music versus animal vocal communication. A new method of such analysis is proposed here. Formation of such expressive aspects as meter, tempo, melodic intervals, and articulation can be explained by the influence of bipedal locomotion, breathing cycle, and heartbeat, long before Homo sapiens. However, two aspects, rhythm and melodic contour, most crucial for music as we know it, lack proxies in the Paleolithic lifestyle. The available ethnographic and developmental data leads one to believe that rhythmic and directional patterns of melody became involved in conveying emotion-related information in the process of frequent switching from one call-type to another within the limited repertory of calls. Such calls are usually adopted for the ongoing caretaking of human youngsters and domestic animals. The efficacy of rhythm and pitch contour in affective communication must have been spontaneously discovered in new important cultural activities. The most likely scenario for music to have become fully semiotically functional and to have spread wide enough to avoid extinctions is the formation of cross-specific communication between humans and domesticated animals during the Neolithic demographic explosion and the subsequent cultural revolution. Changes in distance during such communication must have promoted the integration between different expressive aspects and generated the basic musical grammar. The model of such communication can be found in the surviving tradition of Scandinavian pastoral music - kulning. This article discusses the most likely ways in which such music evolved.

Risk Factors for Outbreaks of Lumpy Skin Disease and the Economic Impact in Cattle Farms of Nakuru County, Kenya

Lumpy Skin Disease (LSD) is an emerging disease of cattle that causes substantial economic loss to affected regions. However, factors favouring transmission under field conditions and farm-level impacts are poorly quantified. This was a retrospective case-control study of cattle farms in Nakuru, Kenya to determine risk factors associated with lumpy skin disease and the farm-level economic impacts of an outbreak. Data were collected using questionnaires administered through personal interview. Collected data included herd sizes, age, and sex structures, breeds, sources of replacement stock, grazing systems, and costs (direct and indirect) incurred when LSD outbreaks occurred. Farm-level risk factors were examined through univariable and multivariable logistic regression and a final model built using backward stepwise regression and likelihood ratio tests. The factors associated with LSD outbreaks on univariable analysis included breed (exotic vs. indigenous, OR = 15.01, P = 0.007), source of replacement stock (outside the herd vs. within the herd, OR = 8.38, P < 0.001) and herd size (large [>10 cattle] vs. small [1–3 cattle], OR = 3.51, P = 0.029). In the multivariable logistic regression model, only breed (exotic vs. indigenous, OR = 14.87, 95% CI 1.94–113.97, P = 0.009) and source of replacement stock (outside the herd vs. within the herd OR = 8.7, 95% CI 2.80–27.0, P < 0.001) were associated with outbreaks. The economic impact was compared between farms keeping purely indigenous (n = 10) or exotic (n = 29) breeds of cattle which indicated mean farm-level losses of 12,431 KSH/123 USD and 76,297 KSH/755 USD, respectively. The mean farm-level losses from reduction in milk yield and mortality were estimated at 4,725 KSH/97 USD and 3,103 KSH/31USD for farms keeping indigenous breeds whilst for farms keeping exotic breeds the equivalent losses were 26,886 KSH/266 USD and 43,557 KSH/431 USD, respectively. The indirect losses from treatments and vaccinations were proportionately much higher on farms with indigenous breeds at 4,603 KSH/46 USD making up ~37% of the total costs compared to ~8% (5,855 KSH/58 USD per farm) of the total costs for farms with exotic breeds. These findings indicate that LSD caused significant economic losses at the farm level in Nakuru County. This justifies implementation of disease control measures including quarantine of cattle post-purchase and the need for effective vaccinations of susceptible cattle herds.

Novel haplotypes responsible for prenatal death in Nordic Red and Danish Jersey cattle

Haplotypes that are common in a population but not observed as homotypes in living animals may harbor lethal alleles that compromise embryo survival. In this study, we searched for homozygous-deficient haplotypes in the genomes of 19,309 Nordic Red Dairy (RDC) and 4,291 Danish Jersey (JER) cattle genotyped using the Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA). For statistically significant deficient haplotypes, we evaluated the effect on nonreturn rate in at-risk matings (mating between carrier bull and daughter of carrier sire) versus not-at-risk matings (mating between noncarrier bull and daughter of noncarrier sire). Next, we analyzed whole-genome sequence variants from the 1000 Bull Genomes Project to identify putative causal variants underlying these haplotypes. In RDC, we identified 3 homozygous-deficient regions (HDR) that overlapped with known recessive lethal mutations: a 662-kb deletion on chromosome 12 in RDC [Online Mendelian Inheritance in Animals (OMIA) 001901-9913), a missense mutation in TUBD1, g.11063520T>C, in Braunvieh cattle (OMIA 001939-9913), and a 525-kb deletion on chromosome 23 in RDC (OMIA 001991-9913)]. In addition, we identified 15 novel HDR and their tag haplotypes for the underlying causative variants. The tag haplotype located between 39.2 and 40.3 Mbp on chromosome 18 had a negative effect on nonreturn rate in at-risk mating, confirming embryonic lethality. In Danish Jersey, we identified 12 novel HDR and their tag haplotypes for underlying causative variants. For 3 of these 12 tag haplotypes, insemination records of at-risk mating showed a negative effect on nonreturn rate, confirming the association with early embryonic mortality. Cattle that had both genotype and whole-genome sequence data were analyzed to detect the causative variants underlying each tag haplotype. However, none of the functional variants or deletions showed concordance with carrier status of the novel tag haplotypes. Carrier status of these detected haplotypes can be used to select bulls to reduce the frequencies of lethal alleles in the population and to avoid at-risk matings.

Organization and Management of Conservation Programs and Research in Domestic Animal Genetic Resources

Creating national committees for domestic animal genetic resources within genetic resource national commissions is recommended to organize in situ and ex situ conservation initiatives. In situ conservation is a high priority because it retains traditional zootechnical contexts and locations to ensure the long-term survival of breeds. In situ actions can be based on subsidies, technical support, structure creation, or trademark definition. Provisional or permanent relocation of breeds may prevent immediate extinction when catastrophes, epizootics, or social conflicts compromise in situ conservation. Ex situ in vivo (animal preservation in rescue or quarantine centers) and in vitro methods (germplasm, tissues/cells, DNA/genes storage) are also potential options. Alert systems must detect emergencies and summon the national committee to implement appropriate procedures. Ex situ coordinated centers must be prepared to permanently or provisionally receive extremely endangered collections. National germplasm banks must maintain sufficient samples of national breeds (duplicated) in their collections to restore extinct populations at levels that guarantee the survival of biodiversity. A conservation management survey, describing national and international governmental and non-governmental structures, was developed. Conservation research initiatives for international domestic animal genetic resources from consortia centralize the efforts of studies on molecular, genomic or geo-evolutionary breed characterization, breed distinction, and functional gene identification. Several consortia also consider ex situ conservation relying on socioeconomic or cultural aspects. The CONBIAND network (Conservation for the Biodiversity of Local Domestic Animals for Sustainable Rural Development) exemplifies conservation efficiency maximization in a low-funding setting, integrating several Latin American consortia with international cooperation where limited human, material, and economic resources are available.

Runs of Homozygosity and NetView analyses provide new insight into the genome-wide diversity and admixture of three German cattle breeds

Angler (RVA) and Red-and-White dual-purpose (RDN) cattle were in the past decades crossed with influential Red Holstein (RH) sires. However, genome-wide diversity studies in these breeds are lacking. The objective of the present study was to elucidate the genome-wide diversity and population structure of the three German cattle breeds. Using 40,851 single nucleotide polymorphism markers scored in 337 individuals, runs of homozygosity (ROH) were analysed in each breed. Clustering and a high-resolution network visualisation analyses were performed on an extended dataset that included 11 additional (outgroup) breeds. Genetic diversity levels were high with observed heterozygosity above 0.35 in all three breeds. Only RVA had a recent past effective population size (Ne) estimate above 100 at 5 generations ago. ROH length distribution followed a similar pattern across breeds and the majority of ROH were found in the length class of >5 to 10 Mb. Estimates of average inbreeding calculated from ROH (FROH) were 0.021 (RVA), 0.045 (RDN) and 0.053 (RH). Moderate to high positive correlations were found between FROH and pedigree inbreeding (FPED) and between FROH and inbreeding derived from the excess of homozygosity (FHOM), while the intercept of the regression of FROH on FPED was above zero. The population structure analysis showed strong evidence of admixture between RVA and RH. Introgression of RDN with RH genes was minimally detected and for the first time, the study uncovered Norwegian Red Cattle ancestry in RVA. Highly heterogeneous genetic background was found for RVA and RH and as expected, the breeds of the extended dataset effectively differentiated mostly based on geographical origin, validating our findings. The results of this study confirm the impact of RH sires on RVA and RDN populations. Furthermore, a close monitoring is suggested to curb further reduction of Ne in the breeds.

A genomic map of climate adaptation in Mediterranean cattle breeds

Domestic species such as cattle (Bos taurus taurus and B. t. indicus) represent attractive biological models to characterize the genetic basis of short-term evolutionary response to climate pressure induced by their post-domestication history. Here, using newly generated dense SNP genotyping data, we assessed the structuring of genetic diversity of 21 autochtonous cattle breeds from the whole Mediterranean basin and performed genome-wide association analyses with covariables discriminating the different Mediterranean climate subtypes. This provided insights into both the demographic and adaptive histories of Mediterranean cattle. In particular, a detailed functional annotation of genes surrounding variants associated with climate variations highlighted several biological functions involved in Mediterranean climate adaptation such as thermotolerance, UV protection, pathogen resistance or metabolism with strong candidate genes identified (e.g., NDUFB3, FBN1, METTL3, LEF1, ANTXR2 and TCF7). Accordingly, our results suggest that main selective pressures affecting cattle in Mediterranean area may have been related to variation in heat and UV exposure, in food resources availability and in exposure to pathogens, such as anthrax bacteria (Bacillus anthracis). Furthermore, the observed contribution of the three main bovine ancestries (indicine, European and African taurine) in these different populations suggested that adaptation to local climate conditions may have either relied on standing genomic variation of taurine origin, or adaptive introgression from indicine origin, depending on the local breed origins. Taken together, our results highlight the genetic uniqueness of local Mediterranean cattle breeds and strongly support conservation of these populations.

Morphological divergence in the West African shorthorn Lagune cattle populations from Benin

This study explored the potential role of agro-ecology in shaping the morphology of Lagune cattle population of Benin. A total of 708 adult Lagune cattle were sampled randomly from nine provenances in two agro-ecological zones (AEZs) and were assessed for eight qualitative and twelve linear body measurements. Data were analyzed using generalized linear model procedures (PROC GLM) followed by the multiple comparison of least square means (LSMEAN) according to the Tukey-Kramer method and multivariate analytical methods, including canonical discrimination analysis (CDA) and hierarchical ascendant classification. Irrespective of AEZ and sex, the body length (102.3 ± 9.31 cm) was greater than the wither height (93.1 ± 7.39 cm) and the body index smaller than 0.85. However, there were significant differences between the two AEZs for most of the measured morphometric and qualitative traits. Moreover, a male-biased sexual size dimorphism was recorded. The CDA based on only four basic body measurements (rump height, body length, heart girth, and ear length) and the calculated Mahalanobis distances suggest that the populations from the two AEZs are distinct and could be further considered ecotypes. Nevertheless, the overall moderate classification rate (70%) of the individual animals into their group of origin indicates interbreeding between the two populations. The pairwise Mahalanobis distances between provenances in the same AEZ were also significant. Together, these results provide supporting evidence for the existence of subdivisions in the Lagune cattle populations from South Benin. The high morphological diversity in the Lagune cattle recorded in the present study could serve as a starting point for the development of efficient selection and sound subpopulation management strategies but also for further phenotypic and genetic characterizations.