Mitochondrial diversity and the origins of African and European cattle

The evolution of contemporary livestock species: Insights from mitochondrial genome

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

Analysing the genetic diversity and population structure of five native Turkish cattle breeds using SNP data

The conservation and sustainable utilization of cattle genetic resources necessitate a comprehensive understanding of their genetic diversity and population structure. This study provides an analysis of five native Turkish cattle breeds: Anatolian Black (ANB), Turkish Grey (TUR), Anatolian Southern Yellow (ASY), East Anatolian Red (EAR), and South Anatolian Red (SAN) using 50 K SNP data. These breeds were compared with three European breeds, Simmental (SIM), Holstein (HOL), and Jersey (JER), and three Asian Zebu breeds: Arabic Zebu (ZAR), Nelore (NEL), and Red Sindhi (RSI). Genetic diversity indices demonstrated moderate heterogeneity among the breeds, with TUR exhibiting the highest observed heterozygosity (Ho = 0.35). Wright's Fst values indicated significant genetic differentiation, particularly between Turkish breeds and both European (Fst = 0.035-0.071) and Asian breeds (Fst = 0.025-0.150). Principal component analysis distinguished the unique genetic profiles of each breed cluster. Admixture analysis revealed degrees of shared genetic ancestry, suggesting historical gene flow between Turkish, European, and Asian breeds. Analysis of molecular variance (AMOVA) attributed approximately 58% of the variation to population differences. Nei's genetic distances highlighted the closer genetic relatedness within Turkish breeds (distance ranges between 0.032 and 0.046) and suggested a more relative affinity of TUR with European breeds. The study's phylogenetic assessments elucidate the nuanced genetic relationships among these breeds, with runs of homozygosity (ROH) analysis indicating patterns of ancestral relatedness and moderate levels of inbreeding, particularly evident in Turkish breeds. Our findings provide valuable insights into the genetic landscape of Turkish cattle, offering a crucial foundation for informed conservation and breeding strategies aimed at preserving these breeds' genetic integrity and heritage.

Uncovering the architecture of selection in two Bos taurus cattle breeds

Directional selection alters the genome via hard sweeps, soft sweeps, and polygenic selection. However, mapping polygenic selection is difficult because it does not leave clear signatures on the genome like a selective sweep. In populations with temporally stratified genotypes, the Generation Proxy Selection Mapping (GPSM) method identifies variants associated with generation number (or appropriate proxy) and thus variants undergoing directional allele frequency changes. Here, we use GPSM on two large datasets of beef cattle to detect associations between an animal's generation and 11 million imputed SNPs. Using these datasets with high power and dense mapping resolution, GPSM detected a total of 294 unique loci actively under selection in two cattle breeds. We observed that GPSM has a high power to detect selection in the very recent past (<10 years), even when allele frequency changes are small. Variants identified by GPSM reside in genomic regions associated with known breed-specific selection objectives, such as fertility and maternal ability in Red Angus, and carcass merit and coat color in Simmental. Over 60% of the selected loci reside in or near (<50 kb) annotated genes. Using haplotype-based and composite selective sweep statistics, we identify hundreds of putative selective sweeps that likely occurred earlier in the evolution of these breeds; however, these sweeps have little overlap with recent polygenic selection. This makes GPSM a complementary approach to sweep detection methods when temporal genotype data are available. The selected loci that we identify across methods demonstrate the complex architecture of selection in domesticated cattle.

Mapping restricted introgression across the genomes of admixed indigenous African cattle breeds

BACKGROUND

The genomes of indigenous African cattle are composed of components with Middle Eastern (taurine) and South Asian (indicine) origins, providing a valuable model to study hybridization and to identify genetic barriers to gene flow. In this study, we analysed indigenous African cattle breeds as models of hybrid zones, considering taurine and indicine samples as ancestors. In a genomic cline analysis of whole-genome sequence data, we considered over 8 million variants from 144 animals, which allows for fine-mapping of potential genomic incompatibilities at high resolution across the genome.

RESULTS

We identified several thousand variants that had significantly steep clines ('SCV') across the whole genome, indicating restricted introgression. Some of the SCV were clustered into extended regions, with the longest on chromosome 7, spanning 725 kb and including 27 genes. We found that variants with a high phenotypic impact (e.g. indels, intra-genic and missense variants) likely represent greater genetic barriers to gene flow. Furthermore, our findings provide evidence that a large proportion of breed differentiation in African cattle could be linked to genomic incompatibilities and reproductive isolation. Functional evaluation of genes with SCV suggest that mitonuclear incompatibilities and genes associated with fitness (e.g. resistance to paratuberculosis) could account for restricted gene flow in indigenous African cattle.

CONCLUSIONS

To our knowledge, this is the first time genomic cline analysis has been applied to identify restricted introgression in the genomes of indigenous African cattle and the results provide extended insights into mechanisms (e.g. genomic incompatibilities) contributing to hybrid differentiation. These results have important implications for our understanding of genetic incompatibilities and reproductive isolation and provide important insights into the impact of cross-breeding cattle with the aim of producing offspring that are both hardy and productive.

Inference of Admixture Origins in Indigenous African Cattle

Present-day African cattle retain a unique genetic profile composed of a mixture of the Bos taurus and Bos indicus populations introduced into the continent at different time periods. However, details of the admixture history and the exact origins of the source populations remain obscure. Here, we infer the source of admixture in the earliest domestic cattle in Africa, African taurine. We detect a significant contribution (up to ∼20%) from a basal taurine lineage, which might represent the now-extinct African aurochs. In addition, we show that the indicine ancestry of African cattle, although most closely related to so-far sampled North Indian indicine breeds, has a small amount of additional genetic affinity to Southeast Asian indicine breeds. Our findings support the hypothesis of aurochs introgression into African taurine and generate a novel hypothesis that the origin of indicine ancestry in Africa might be different indicine populations than the ones found in North India today.

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.

Concurrently mapping quantitative trait loci associations from multiple subspecies within hybrid populations

Many of the world's agriculturally important plant and animal populations consist of hybrids of subspecies. Cattle in tropical and sub-tropical regions for example, originate from two subspecies, Bos taurus indicus (Bos indicus) and Bos taurus taurus (Bos taurus). Methods to derive the underlying genetic architecture for these two subspecies are essential to develop accurate genomic predictions in these hybrid populations. We propose a novel method to achieve this. First, we use haplotypes to assign SNP alleles to ancestral subspecies of origin in a multi-breed and multi-subspecies population. Then we use a BayesR framework to allow SNP alleles originating from the different subspecies differing effects. Applying this method in a composite population of B. indicus and B. taurus hybrids, our results show that there are underlying genomic differences between the two subspecies, and these effects are not identified in multi-breed genomic evaluations that do not account for subspecies of origin effects. The method slightly improved the accuracy of genomic prediction. More significantly, by allocating SNP alleles to ancestral subspecies of origin, we were able to identify four SNP with high posterior probabilities of inclusion that have not been previously associated with cattle fertility and were close to genes associated with fertility in other species. These results show that haplotypes can be used to trace subspecies of origin through the genome of this hybrid population and, in conjunction with our novel Bayesian analysis, subspecies SNP allele allocation can be used to increase the accuracy of QTL association mapping in genetically diverse populations.

Iron age genomic data from Althiburos - Tunisia renew the debate on the origins of African taurine cattle

The Maghreb is a key region for understanding the dynamics of cattle dispersal and admixture with local aurochs following their earliest domestication in the Fertile Crescent more than 10,000 years ago. Here, we present data on autosomal genomes and mitogenomes obtained for four archaeological specimens of Iron Age (∼2,800 cal BP-2,000 cal BP) domestic cattle from the Eastern Maghreb, i.e. Althiburos (El Kef, Tunisia). D-loop sequences were obtained for an additional eight cattle specimens from this site. Maternal lineages were assigned to the elusive R and ubiquitous African-T1 haplogroups found in two and ten Althiburos specimens, respectively. Our results can be explained by post-domestication hybridization of Althiburos cattle with local aurochs. However, we cannot rule out an independent domestication in North Africa considering the shared ancestry of Althiburos cattle with the pre-domestic Moroccan aurochs and present-day African taurine cattle.

ddRAD sequencing based genotyping of six indigenous dairy cattle breeds of India to infer existing genetic diversity and population structure

The present investigation aimed to identify genome wide SNPs and to carry out diversity and population structure study using ddRAD-seq based genotyping of 58 individuals of six indigenous milch cattle breeds (Bos indicus) such as Sahiwal, Gir, Rathi, Tharparkar, Red Sindhi and Kankrej of India. A high percentage of reads (94.53%) were mapped to the Bos taurus (ARS-UCD1.2) reference genome assembly. Following filtration criteria, a total of 84,027 high quality SNPs were identified across the genome of 6 cattle breeds with the highest number of SNPs observed in Gir (34,743), followed by Red Sindhi (13,092), Kankrej (12,812), Sahiwal (8956), Tharparkar (7356) and Rathi (7068). Most of these SNPs were distributed in the intronic regions (53.87%) followed by intergenic regions (34.94%) while only 1.23% were located in the exonic regions. Together with analysis of nucleotide diversity (π = 0.373), Tajima's D (D value ranging from - 0.295 to 0.214), observed heterozygosity (HO ranging from 0.464 to 0.551), inbreeding coefficient (FIS ranging from - 0.253 to 0.0513) suggested for the presence of sufficient within breed diversity in the 6 major milch breeds of India. The phylogenetic based structuring, principal component and admixture analysis revealed genetic distinctness as well as purity of almost all of the 6 cattle breeds. Overall, our strategy has successfully identified thousands of high-quality genome wide SNPs that will further enrich the Bos indicus representation basic information about genetic diversity and structure of 6 major Indian milch cattle breeds which should have implications for better management and conservation of valuable indicine cattle diversity.

Profiling the immune epigenome across global cattle breeds

BACKGROUND

Understanding the variation between well and poorly adapted cattle breeds to local environments and pathogens is essential for breeding cattle with improved climate and disease-resistant phenotypes. Although considerable progress has been made towards identifying genetic differences between breeds, variation at the epigenetic and chromatin levels remains poorly characterized. Here, we generate, sequence and analyse over 150 libraries at base-pair resolution to explore the dynamics of DNA methylation and chromatin accessibility of the bovine immune system across three distinct cattle lineages.

RESULTS

We find extensive epigenetic divergence between the taurine and indicine cattle breeds across immune cell types, which is linked to the levels of local DNA sequence divergence between the two cattle sub-species. The unique cell type profiles enable the deconvolution of complex cellular mixtures using digital cytometry approaches. Finally, we show distinct sub-categories of CpG islands based on their chromatin and methylation profiles that discriminate between classes of distal and gene proximal islands linked to discrete transcriptional states.

CONCLUSIONS

Our study provides a comprehensive resource of DNA methylation, chromatin accessibility and RNA expression profiles of three diverse cattle populations. The findings have important implications, from understanding how genetic editing across breeds, and consequently regulatory backgrounds, may have distinct impacts to designing effective cattle epigenome-wide association studies in non-European breeds.

Revealing the diversity of internal body temperature and panting response for feedlot cattle under environmental thermal stress

Core body temperature (CBT) regulation is crucial for mammalian wellbeing and survival. Cattle pant to dissipate excess heat to regulate CBT when ambient conditions exceed thermoneutral zones. However, to date, neither the variability in cattle heat response, the lagged response of CBT to thermal indices, nor the diurnal patterns of thermal indices, CBT and panting have been reported in the literature. We decomposed thermal indices, CBT and panting time-series data for 99 feedlot heifers across three discrete heat events into diurnal, trend and residual components. Both raw and decomposed data were analysed to explore the lagged CBT and panting responses and the association between series. We show ambient thermal conditions impact CBT with a 1-h lag despite a lag of between 1.5 to 3 h from raw data. Average individual panting scores were used to identify heat-susceptible and heat-tolerant cattle. Heat-susceptible cattle showed greater CBT (P < 0.01) between 8:00 and 23:00 and greater panting duration (P < 0.05) between 10:00 and 18:00 than heat-tolerant cattle under the same thermal conditions and these variations followed a similar pattern despite differences in cattle breed. This new information enables targeted amelioration and selection of individuals against heat susceptibility.

Peeping into Mitochondrial Diversity of Andaman Goats: Unveils Possibility of Maritime Transport with Diversified Geographic Signaling

Andaman and Nicobar Islands, a part of South-East Asia, is enriched with the presence of native breeds of livestock (cattle, pig, goat) and poultry. There are two native goat breeds, viz., Andaman local goat and Teressa goat in Andaman and Nicobar Islands. However, to date, the origin and genetic makeup of these two breeds have not been detailed. Therefore, the present study describes the genetic makeup of Andaman goats through analysis of mitochondrial D-loop sequence for sequence polymorphism, phylogeographical signaling and population expansion events. The genetic diversity of the Teressa goat was less compared to the Andaman local goat due to its sole presence on Teressa Island. Out of 38 well-defined haplotypes of Andaman goats, the majority of haplotypes belonged to haplogroup A followed by haplogroup B and haplogroup D. The result of mismatch distribution and neutrality tests indicated no population expansion event of haplogroup A and B. Finally, based on poor geographical signaling, we hypothesize that Andaman goats have been imported to these Islands either through multidirectional diffusion or unidirectional diffusion. We justify our hypothesis of multidirectional diffusion on the basis of observation of the haplotype and nucleotide diversity of Andaman goats. Simultaneously, the probability of unidirectional diffusion of goats in these islands from the Indian subcontinent in different spells of domestication events through maritime routes cannot be ignored.

Detection of genomic regions that differentiate Bos indicus from Bos taurus ancestral breeds for milk yield in Indian crossbred cows

Introduction: In India, crossbred cows incorporate the high production of B. taurus dairy breeds and the environmental adaptation of local B. indicus cattle. Adaptation to different environments and selection in milk production have shaped the genetic differences between B. indicus and B. taurus cattle. The aim of this paper was to detect, for milk yield of crossbred cows, quantitative trait loci (QTL) that differentiate B. indicus from B. taurus ancestry, as well as QTL that are segregating within the ancestral breeds. Methods: A total of 123,042 test-day milk records for 4,968 crossbred cows, genotyped with real and imputed 770 K SNP, were used. Breed origins were assigned to haplotypes of crossbred cows, and from that, were assigned to SNP alleles. Results: At a false discovery rate (FDR) of 30%, a large number of genomic regions showed significant effects of B. indicus versus B. taurus origin on milk yield, with positive effects coming from both ancestors. No significant regions were detected for Holstein Friesian (HF) versus Jersey effects on milk yield. Additionally, no regions for SNP alleles segregating within indigenous, within HF, and within Jersey were detected. The most significant effects, at FDR 5%, were found in a region on BTA5 (43.98-49.44 Mbp) that differentiates B. indicus from B. taurus, with an estimated difference between homozygotes of approximately 10% of average yield, in favour of B. indicus origin. Discussion: Our results indicate that evolutionary differences between B. indicus and B. taurus cattle for milk yield, as expressed in crossbred cows, occur at many causative loci across the genome. Although subject to the usual first estimation bias, some of the loci appear to have large effects that might make them useful for genomic selection in crossbreds, if confirmed in subsequent studies.

Genomics, Origin and Selection Signals of Loudi Cattle in Central Hunan

Due to the geographical, cultural and environmental variability in Xiangxi, China, distinctive indigenous cattle populations have formed. Among them, Loudi cattle and Xiangxi cattle are the local cattle in Hunan, and the environment in Loudi is relatively more enclosed and humid than that in Xiangxi. To study the genome and origin of Loudi cattle in hot and humid environments, 29 individuals were collected and sequenced by whole-genome resequencing. In addition, genomic data were obtained from public databases for 96 individuals representing different cattle breeds worldwide, including 23 Xiangxi cattle from western Hunan. Genetic analysis indicated that the genetic diversity of Loudi cattle was close to that of Chinese cattle and higher than that of other breeds. Population structure and ancestral origin analysis indicated the relationship between Loudi cattle and other breeds. Loudi has four distinctive seasons, with a stereoscopic climate and extremely rich water resources. Selective sweep analysis revealed candidate genes and pathways associated with environmental adaptation and homeostasis. Our findings provide a valuable source of information on the genetic diversity of Loudi cattle and ideas for population conservation and genome-associated breeding of local cattle in today's extreme climate environment.

Genome-wide local ancestry and evidence for mitonuclear coadaptation in African hybrid cattle populations (Bos taurus/indicus)

The phenotypic diversity of African cattle reflects adaptation to a wide range of agroecological conditions, human-mediated selection preferences, and complex patterns of admixture between the humpless Bos taurus (taurine) and humped Bos indicus (zebu) subspecies, which diverged 150-500 thousand years ago. Despite extensive admixture, all African cattle possess taurine mitochondrial haplotypes, even populations with significant zebu biparental and male uniparental nuclear ancestry. This has been interpreted as the result of human-mediated dispersal ultimately stemming from zebu bulls imported from South Asia during the last three millennia. Here, we assess whether ancestry at mitochondrially targeted nuclear genes in African admixed cattle is impacted by mitonuclear functional interactions. Using high-density SNP data, we find evidence for mitonuclear coevolution across hybrid African cattle populations with a significant increase of taurine ancestry at mitochondrially targeted nuclear genes. Our results, therefore, support the hypothesis of incompatibility between the taurine mitochondrial genome and the zebu nuclear genome.

Unraveling genetic admixture in the Indian crossbred cattle by different approaches using Bovine 50K BeadChip

With the upsurge of crossbreeding in India, the admixture levels are highly unpredictable in the composite breeds. Hence, in the present study, 72 Vrindavani animals were assessed for the level of admixture from their known ancestors that are Holstein-Friesian, Jersey, Brown Swiss, and Hariana, through three different software, namely, STRUCTURE, ADMIXTURE, and frappe. The genotype data for ancestral breeds were obtained from a public repository, i.e., DRYAD. The Frieswal crossbred cattle along with ancestral breeds like Holstein-Friesian and Sahiwal were also investigated for the level of admixture with the help of the above-mentioned software. The Frieswal population was found to comprise an average of 62.49, 61.12, and 61.21% of Holstein-Friesian and 37.50, 38.88, and 38.80% of Sahiwal estimated through STRUCTURE, ADMIXTURE, and frappe, respectively. The Vrindavani population was found to consist of on average 39.5, 42.4, and 42.3% of Holstein-Friesian; 22.9, 22.3, and 21.7% of Jersey; 10.7, 10.6, and 11.9% of Brown Swiss; and 26.9, 24.7, and 24.1% of Hariana blood estimated through STRUCTURE, ADMIXTURE, and frappe, respectively. A greater degree of variation was noted in the results from STRUCTURE vs. frappe, STRUCTURE vs. ADMIXTURE than in ADMIXTURE vs. frappe. From this study, we conclude that the admixture analysis based on a single software should be validated through the use of many different approaches for better prediction of admixture levels.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Mutant alleles differentially shape fitness and other complex traits in cattle

Mutant alleles (MAs) that have been classically recognised have large effects on phenotype and tend to be deleterious to traits and fitness. Is this the case for mutations with small effects? We infer MAs for 8 million sequence variants in 113k cattle and quantify the effects of MA on 37 complex traits. Heterozygosity for variants at genomic sites conserved across 100 vertebrate species increase fertility, stature, and milk production, positively associating these traits with fitness. MAs decrease stature and fat and protein concentration in milk, but increase gestation length and somatic cell count in milk (the latter indicative of mastitis). However, the frequency of MAs decreasing stature and fat and protein concentration, increasing gestation length and somatic cell count were lower than the frequency of MAs with the opposite effect. These results suggest bias in the mutations direction of effect (e.g. towards reduced protein in milk), but selection operating to reduce the frequency of these MAs. Taken together, our results imply two classes of genomic sites subject to long-term selection: sites conserved across vertebrates show hybrid vigour while sites subject to less long-term selection show a bias in mutation towards undesirable alleles.

Demographic expansion and high level of matrilineal diversity in two populations of East African Baggara cattle

Western Baggara cattle breed (WBCB) is an East African zebu inhabiting Sudan, well-known as beef-producing cattle. We investigated herein two phenotypically and geographically distinct populations of this breed, namely Nyalawi and Daeinawi, which are renowned for their unique meat production capabilities and adaptation attributes, with the aim to contribute to our understanding of their maternal genetic diversity and demography dynamics. Genetic polymorphism analysis of the full-length D-loop mtDNA region revealed 44 and 35 polymorphic sites defining 28 and 24 distinct haplotypes in the Nyalawi and the Daeinawi, respectively. Observed genetic diversity is high within the population with a low level of genetic differentiation between populations. Approximate Bayesian computation via the calculation of Bayesian skyline plots and neutrality tests support past expansion with a higher maternal effective population size (N_e ) in Nyalawi compared with the Daeinawi population and a population expansion beginning around 4,500 YBP and 3,500 YBP, respectively, before the arrival of zebu into the continent.

BoLA-DRB3 gene haplotypes show divergence in native Sudanese cattle from taurine and indicine breeds

Autochthonous Sudanese cattle breeds, namely Baggara for beef and Butana and Kenana for dairy, are characterized by their adaptive characteristics and high performance in hot and dry agro-ecosystems. They are thus used largely by nomadic and semi-nomadic pastoralists. We analyzed the diversity and genetic structure of the BoLA-DRB3 gene, a genetic locus linked to the immune response, for the indigenous cattle of Sudan and in the context of the global cattle repository. Blood samples (n = 225) were taken from three indigenous breeds (Baggara; n = 113, Butana; n = 60 and Kenana; n = 52) distributed across six regions of Sudan. Nucleotide sequences were genotyped using the sequence-based typing method. We describe 53 alleles, including seven novel alleles. Principal component analysis (PCA) of the protein pockets implicated in the antigen-binding function of the MHC complex revealed that pockets 4 and 9 (respectively) differentiate Kenana-Baggara and Kenana-Butana breeds from other breeds. Venn analysis of Sudanese, Southeast Asian, European and American cattle breeds with 115 alleles showed 14 were unique to Sudanese breeds. Gene frequency distributions of Baggara cattle showed an even distribution suggesting balancing selection, while the selection index (ω) revealed the presence of diversifying selection in several amino acid sites along the BoLA-DRB3 exon 2 of these native breeds. The results of several PCA were in agreement with clustering patterns observed on the neighbor joining (NJ) trees. These results provide insight into their high survival rate for different tropical diseases and their reproductive capacity in Sudan's harsh environment.

Meta-Analysis of Mitochondrial DNA Control Region Diversity to Shed Light on Phylogenetic Relationship and Demographic History of African Sheep (Ovis aries) Breeds

To improve sheep breeding and conservation of genetic resources, the mitochondrial DNA control region (mtDNA CR) of 399 sequences of African indigenous sheep breeds from previously published research articles were meta-analyzed to elucidate their phylogenetic relationship, diversity, and demographic history. A total of 272 haplotypes were found, of which 207 were unique and a high level of mtDNA CR variability was observed. Generally, the number of polymorphic sites, nucleotide and haplotype diversity were high (284, 0.254 ± 0.012 and 0.993 ± 0.002, respectively). The median-joining (MJ) network of haplotypes produced three major haplogroups (A, B and C), with haplogroup B being dominant. A mixture of populations suggests a common matrilineal origin and lack of and/or a weak phylogeographic structure. Mismatch analysis showed recent expansion of North African breeds, whereas East African and continental populations exhibited selection pressures for adaptation. A slight historical genetic difference was also observed between the fat tail and thin tail sheep breeds. However, further investigations are required using more samples and long sequence segments to achieve deeper levels of conclusions on the African sheep phylogenetic relationship. The present meta-analysis results contribute to the general understanding of African native sheep populations for improved management of sheep genetic resources.

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.

Genetic diversity and effective population sizes of thirteen Indian cattle breeds

Background

The genetic structure of a diverse set of 15 Indian indigenous breeds and non-descript indigenous cattle sampled from eight states was examined, based on 777 k single nucleotide polymorphism (SNP) genotypes obtained on 699 animals, with sample sizes ranging from 17 to 140 animals per breed. To date, this is the largest and most detailed assessment of the genetic diversity of Indian cattle breeds.

Results

Admixture analyses revealed that 109 of the indigenous animals analyzed had more than 1% Bos taurus admixture of relatively recent origin. Pure indigenous animals were defined as having more than 99% Bos indicus ancestry. Assessment of the genetic diversity within and between breeds using principal component analyses, F statistics, runs of homozygosity, the genomic relationship matrix, and maximum likelihood clustering based on allele frequencies revealed a low level of genetic diversity among the indigenous breeds compared to that of Bos taurus breeds. Correlations of SNP allele frequencies between breeds indicated that the genetic variation among the Bos indicus breeds was remarkably low. In addition, the variance in allele frequencies represented less than 1.5% between the Indian indigenous breeds compared to about 40% between Bos taurus dairy breeds. Effective population sizes (N_e) increased during a period post-domestication, notably for Ongole cattle, and then declined during the last 100 generations. Although we found that most of the identified runs of homozygosity are short in the Indian indigenous breeds, indicating no recent inbreeding, the high F_ROH coefficients and low F_IS values point towards small population sizes. Nonetheless, the N_e of the Indian indigenous breeds is currently still larger than that of Bos taurus dairy breeds.

Conclusions

The changes in the estimates of effective population size are consistent with domestication from a large native population followed by consolidation into breeds with a more limited population size. The surprisingly low genetic diversity among Indian indigenous cattle breeds might be due to their large N_e since their domestication, which started to decline only 100 generations ago, compared to approximately 250 to 500 generations for Bos taurus dairy cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00640-3.

Genetic Diversity of Taenia solium and its Relation to Clinical Presentation of Cysticercosis

In this perspectives paper, we discuss fertilization strategies for Taenia saginata and Taenia saginata asiatica as well as heterogeneity in Taenia solium, the causative agent of human cysticercosis. Two different genotypes of T. solium (Asian and Afro/American) were confirmed by mitochondrial DNA analysis approximately two decades ago. Since then, outcrossings of the two genotypes have been identified in Madagascar where the two genotypes are distributed sympatrically. Outcrossings were confirmed by the presence of discordance between mitochondrial and nuclear DNA. Since multiple tapeworm infections are common in endemic areas, outcrossing events likely occur quite frequently. Therefore, mitochondrial DNA from T. solium specimens collected from humans and pigs in endemic areas should be analyzed. If variations are found between specimens, nuclear DNA analysis should be performed to confirm the presence of discordance between mitochondrial and nuclear genes. Additional outcrossings likely add complexity to understanding the existing genetic diversity. Serological surveys are also recommended since serodiagnostic glycoprotein can also differentiate between the two genotypes. Viable eggs from different genotypes or from hybrids of two different genotypes should be used for experimental infection of pigs or dogs in order to observe any pathological heterogeneity in cysticercosis development. Although genetic diversity of T. solium is expected to result in clinical heterogeneity of cysticercosis in humans and pigs, there is currently no evidence showing that this occurs. There are also no comparative experimental studies on this topic. Therefore, studies evaluating the link between parasite heterogeneity and clinical outcome are warranted.

What is a rare disease in animal paleopathology?

OBJECTIVE

This paper seeks to review the possibilities and difficulties in identifying rare diseases in ancient animals. Rare diseases are important in human medicine but seldom discussed in modern veterinary practice.

MATERIALS AND METHODS

A total of 1348 pathological archaeozoological remains, most published in the literature, are evaluated. With the exception of a few animal burials, the majority were recovered from food refuse: deposits composed of co-mingled, single bone fragments.

RESULTS

Determination of population-based prevalence is not applicable in animal paleopathology, as almost all lesions occur on isolated bones. Moreover, grave hereditary diseases are rarely detected on animal bones because animals with such disorders seldom survived, except when humans promoted their rare inherited traits.

SIGNIFICANCE

Rare diseases form a special category in human pathology posing both therapeutic and ethical challenges. While in wild animals natural selection tends to prevent the inheritance of such conditions, curious cases of animal morbidity have been brought about by domestication. Humans sheltered animals of lesser vitality and sometimes even promoted their negative traits. Understanding these phenomena in animal paleopathology will help fine-tuning the rare disease paradigm.

LIMITATIONS

The definition of rare disease in animal paleopathology can only be assessed based on ancient and modern human correlates, and rare variation could become cultivated traits, ultimately developed into part of "normal" variability as trademark breed characteristics. Taphonomic limitations in recovering osteological evidence of debilitating hereditary diseases in animals are unlikely to improve. Suggestions for further research: Further research will need to focus on the scarce osteological evidence for rare conditions in light of differential diagnoses. The concept of "rare disease" may be worth expanding beyond hereditary conditions in animals considering fundamental differences between past animal and human lifeways.

Production, reproduction and some adaptation characteristics of Boran cattle breed under changing climate: A systematic review and meta-analysis

Introduction

Climate change affects livestock production and productivity, which could threaten livestock-based food security in pastoral and agro-pastoral production systems of the tropics and sub-tropics. Boran cattle breed is one of the hardiest Zebu cattle reared by Borana Oromo pastoralists for milk and meat production. However, there is limited comprensive information on production, reproduction and adaption traits of the Boran cattle in Ethiopia. Thus, this paper aims to compile the main production, reproduction and some adaptation traits of Boran cattle based on systematic review and meta-analysis of peer reviewed published and unpublished literature.

Methodology

A combination of systematic review and meta-analysis based on PRISMA guideline was employed. Accordingly, out of 646 recorded articles identified through database searching, 64 were found to be eligible for production, reproduction and adaptation characteristics of the Boran cattle, 28 articles were included in qualitative systematic review while 36 articles were used for quantitative meta-analysis.

Result

The Boran cattle breed has the ability to survive, produce and reproduce under high ambient temperature, utilize low quality forage resources, and resist water shortage or long watering intervals and tick infestations. The review revealed that the breed employs various adaptation responses (morphological, physiological, biochemical, metabolic, cellular and molecular responses) to cope with harsh environmental conditions including climate change, rangeland degradation, seasonal feed and water shortages and high incidences of tick infestations. The meta-analysis using a random-effects model allowed provision of pooled estimates of heritability and genetic correlations for reproduction and production traits, which could be used to solve genetic prediction equations under a population level in purebred Boran cattle. In addition, heritability and genetic-correlation estimates found in the present study suggest that there is high genetic variability for most traits in Boran cattle, and that genetic progress is possible for all studied traits in this breed.

Conclusion

The Boran cattle breed has the ability to survive, produce and reproduce under high ambient temperature, utilize low quality forage resources, and resist water shortage or long watering intervals and tick infestations. However, currently there are several challenges such as recurrent droughts, pasture deterioration and lack of systematic selection and breeding programs that play to undermine the realization of the potential of the breed. Thus, we recommend systematic selection for enhancing the reproductive and production performances without compromising the adaptation traits of the breed coupled with improved management of rangelands.

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

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

Insight into the introduction of domestic cattle and the process of Neolithization to the Spanish region Galicia by genetic evidence

Domestic cattle were brought to Spain by early settlers and agricultural societies. Due to missing Neolithic sites in the Spanish region of Galicia, very little is known about this process in this region. We sampled 18 cattle subfossils from different ages and different mountain caves in Galicia, of which 11 were subject to sequencing of the mitochondrial genome and phylogenetic analysis, to provide insight into the introduction of cattle to this region. We detected high similarity between samples from different time periods and were able to compare the time frame of the first domesticated cattle in Galicia to data from the connecting region of Cantabria to show a plausible connection between the Neolithization of these two regions. Our data shows a close relationship of the early domesticated cattle of Galicia and modern cow breeds and gives a general insight into cattle phylogeny. We conclude that settlers migrated to this region of Spain from Europe and introduced common European breeds to Galicia.

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

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

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

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

SNP panels for the estimation of dairy breed proportion and parentage assignment in African crossbred dairy cattle

BACKGROUND

Understanding the relationship between genetic admixture and phenotypic performance is crucial for the optimization of crossbreeding programs. The use of small sets of informative ancestry markers can be a cost-effective option for the estimation of breed composition and for parentage assignment in situations where pedigree recording is difficult. The objectives of this study were to develop small single nucleotide polymorphism (SNP) panels that can accurately estimate the total dairy proportion and assign parentage in both West and East African crossbred dairy cows.

METHODS

Medium- and high-density SNP genotype data (Illumina BovineSNP50 and BovineHD Beadchip) for 4231 animals sampled from African crossbreds, African Bos taurus, European Bos taurus, Bos indicus, and African indigenous populations were used. For estimating breed composition, the absolute differences in allele frequency were calculated between pure ancestral breeds to identify SNPs with the highest discriminating power, and different combinations of SNPs weighted by ancestral origin were tested against estimates based on all available SNPs. For parentage assignment, informative SNPs were selected based on the highest minor allele frequency (MAF) in African crossbred populations assuming two Scenarios: (1) parents were selected among all the animals with known genotypes, and (2) parents were selected only among the animals known to be a parent of at least one progeny.

RESULTS

For the medium-density genotype data, SNPs selected for the largest differences in allele frequency between West African indigenous and European Bos taurus breeds performed best for most African crossbred populations and achieved a prediction accuracy (r2) for breed composition of 0.926 to 0.961 with 200 SNPs. For the high-density dataset, a panel with 70% of the SNPs selected on their largest difference in allele frequency between African and European Bos taurus performed best or very near best across all crossbred populations with r2 ranging from 0.978 to 0.984 with 200 SNPs. In all African crossbred populations, unambiguous parentage assignment was possible with ≥ 300 SNPs for the majority of the panels for Scenario 1 and ≥ 200 SNPs for Scenario 2.

CONCLUSIONS

The identified low-cost SNP assays could overcome incomplete or inaccurate pedigree records in African smallholder systems and allow effective breeding decisions to produce progeny of desired breed composition.

Selection and introgression facilitated the adaptation of Chinese native endangered cattle in extreme environments

Although persistent efforts have identified and characterized a few candidate genes and related biological processes with potential functions in the adaptation of many species to extreme environments, few works have been conducted to determine the genomic basis of adaptation in endangered livestock breeds that have been living in extreme conditions for more than thousands of years. To fill this gap, we sequenced the whole genomes of nine individuals from three Chinese native endangered cattle breeds that are living in high-altitude or arid environments. Phylogenetic and evolutionary history analyses of these three and other six breeds showed that the genetic structure of the cattle populations is primarily related to geographic location. Interestingly, we identified pervasive introgression from the yak to Zhangmu cattle (ZMC) that cover several genes (e.g., NOS2, EGLN1 and EPAS1) involved in the hypoxia response and previously identified as positive selection genes in other species, which suggested that the adaptive introgression from yak may have contributed to the adaptation of ZMC to high-altitude environments. In addition, by contrasting the breeds in opposite living conditions, we revealed a set of candidate genes with various functions from hypoxia response, water metabolism, immune response and body shape change to embryo development and skeletal system development, etc., that may be related to high-altitude or arid adaptation. Our research provides new insights into the recovery and adaptation of endangered native cattle and other species in extreme environments and valuable resources for future research on animal husbandry to cope with climate change.

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.

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.

Molecular characterization and phylogenetic analyses of Fasciola gigantica of buffaloes and goats in Punjab, Pakistan

Fasciola gigantica is considered to be a major pathogen causing fasciolosis in the Indian subcontinent, resulting in production losses of millions of dollars in the livestock industry. Understading the dispersal origin and the patterns of spread of F. gigantica is important. A total of 53 Fasciola flukes collected from buffaloes and goats in Punjab, Pakistan between 2017 and 2018 were identified as F. gigantica based on the multiplex PCR for the phosphoenolpyruvate carboxykinase (pepck) and the PCR-restriction fragment length polymorphism (RFLP) for DNA polymerase delta (pold). A significant genetic difference between F. gigantica from buffaloes and goats was indicated by the genetic analyses of mitochondrial markers, NADH dehydrogenase subunit 1 (nad1) and cytochrome C oxidase subunit 1 (cox1). Phylogenetic analysis of the seventeen nad1 haplotypes of F. gigantica from Pakistan with those in neighbouring countries of the Indian subcontinent revealed that all the haplotypes identified in Pakistan were clustered in haplogroup A. fasciola gigantica with the eight haplotypes might be expanded in Pakistan from Indian origin, along with the migration of the domestic animals, since they were related to Indian haplotypes. In contrast, the remaining nine haplotypes were not shared with any neighbouring countries, suggesting independent origin, probably from neighbouring Middle East countries. However, cautious interpretation is required due to the very limited samples size of this study. Our study provides a proof of concept for a method that could be used to investigate the epidemiology of F. gigantica.

Effects of nutrient restriction on the metabolic profile of Bos indicus-influenced and B. taurus suckled beef cows

Recent research from our group demonstrated that Bos indicus-influenced suckled beef cows had greater resilience to withstand nutrient restriction and establish pregnancy compared with B. taurus cows exposed to the same conditions. To further understand these findings, differences in metabolic profile between these same B. indicus-influenced and B. taurus females were explored. Suckled beef cows (n = 134) were enrolled in a completely randomized design with a 2 × 2 factorial arrangement of treatments. On day -21, Angus (AN; Bos taurus) and Brangus (BN; B. indicus-influenced) cows were randomly assigned to 1) a diet that met daily energy maintenance requirements (MAINT), or 2) a diet that restricted intake to 70% of the daily energy maintenance requirements (RESTR). Cows were exposed to an estrus synchronization protocol and received an embryo 7 d after ovulation was pharmacologically induced on day 0. Blood samples were collected on days -21 and 19 to determine circulating concentrations of non-esterified fatty acids (NEFA), β-hydroxybutyrate (BHB), insulin, glucose, and IGF-1. Pregnancy status after embryo transfer was determined on day 28. As a consequence of the proposed diets, cows in the RESTR diet had less body condition score (BCS) on day 19 (P = 0.008) across breed types. Moreover, BCS change from day -21 to 19 was included as independent covariate into subsequent analyses, allowing for the comparison of breed types under an equivalent level of body reserve mobilization. A breed × diet interaction was observed for plasma insulin (P = 0.03) and IGF-1 (P = 0.04) on day 19, where AN-RESTR cows had less plasma concentrations on day 19 compared with AN-MAINT cows. Diets did not impact (P > 0.10) plasma insulin and IGF-1 concentrations in BN cows. No diet or breed effects were observed in circulating concentrations of NEFA, BHB, and glucose (P > 0.10). Across breed types and nutritional treatment, there was positive linear effect (P ≤ 0.04) of plasma concentrations of insulin and IGF-1 on the probability of pregnancy to fixed-time embryo transfer. In summary, the negative impacts of nutrient restriction on the somatotropic axis, independently of body tissue mobilization, were heightened in Bos taurus females compared with B. indicus-influenced cohorts, which corroborate with the differences observed in fertility between these subspecies.

Discovery of Genomic Characteristics and Selection Signatures in Southern Chinese Local Cattle

Chinese local cattle with a high level of genetic diversity mainly originate from two subspecies; the cattle in northern China are primarily Bos Taurus, and the cattle in southern China are primarily Bos indicus. Cattle from southern China are characterized by a specific phenotype and adapted to the local environment. This study explored the genetic diversity, degree of admixture, and selection signature in eight local cattle breeds in southern China. The lowest level of heterozygosity was found in Hainan and Nandan cattle from Hainan and Guangxi province, respectively, whereas the highest level of heterozygosity was detected in Zhaotong cattle from Yunnan province. A neighbor-joining phylogenetic tree analysis clearly separated Lufeng cattle from other breeds, whereas Leiqiong and Hainan cattle have some crossover. Based on linkage disequilibrium-filtered single nucleotide polymorphisms (SNPs), the admixture analysis revealed two clusters corresponding to the taurine and indicine cattle lineages, and the local cattle breeds from southern China showed a certain degree of admixture. When K = 4 and 9, we found a slight separation among Leiqiong, Lufeng, and Hainan cattle. Meanwhile, we performed a selection signature analysis in Hainan, Leiqiong, and Lufeng cattle distributed in the extreme south of China, using the integrated haplotype score (iHS), Rsb statistic, and BayeScan software. Using the iHS approach, we identified 251, 270, and 256 candidate regions in Lufeng, Leiqiong, and Hainan cattle, respectively. Moreover, we identified 184, 174, and 146 candidate regions in pairwise comparisons of Leiqiong vs. Lufeng, Leiqiong vs. Hainan, and Hainan vs. Lufeng cattle using the Rsb approach. In addition, we identified 76 loci with a total of 48 genes under selection, based on the F_ST approach. Several candidate genes under selection were found to be related to meat quality, immunity, and adaptation to the local environment in southern China. Our results provide significant information about the genetic differences among the cattle breeds from southern China and the possible cause of difference in breed-specific characteristics. Selection signature analysis identified a few candidate SNPs and genes related to certain important traits of these cattle. In general, our results provide valuable insights into the genetic basis of specific traits under selection in certain local cattle breeds.

Milk Containing A2 β-Casein ONLY, as a Single Meal, Causes Fewer Symptoms of Lactose Intolerance than Milk Containing A1 and A2 β-Caseins in Subjects with Lactose Maldigestion and Intolerance: A Randomized, Double-Blind, Crossover Trial

Acute-feeding and multiple-day studies have demonstrated that milk containing A2 β-casein only causes fewer symptoms of lactose intolerance (LI) than milk containing both A1 and A2 β-caseins. We conducted a single-meal study to evaluate the gastrointestinal (GI) tolerance of milk containing different concentrations of A1 and A2 β-casein proteins. This was a randomized, double-blind, crossover trial in 25 LI subjects with maldigestion and an additional eight lactose maldigesters who did not meet the QLCSS criteria. Subjects received each of four types of milk (milk containing A2 β-casein protein only, Jersey milk, conventional milk, and lactose-free milk) after overnight fasting. Symptoms of GI intolerance and breath hydrogen concentrations were analyzed for 6 h after ingestion of each type of milk. In an analysis of the 25 LI subjects, total symptom score for abdominal pain was lower following consumption of milk containing A2 β-casein only, compared with conventional milk (p = 0.004). Post hoc analysis with lactose maldigesters revealed statistically significantly improved symptom scores (p = 0.04) and lower hydrogen production (p = 0.04) following consumption of milk containing A2 β-casein only compared with conventional milk. Consumption of milk containing A2 β-casein only is associated with fewer GI symptoms than consumption of conventional milk in lactose maldigesters.

The patterns of admixture, divergence, and ancestry of African cattle populations determined from genome-wide SNP data

Background

Humpless Bos taurus cattle are one of the earliest domestic cattle in Africa, followed by the arrival of humped Bos indicus cattle. The diverse indigenous cattle breeds of Africa are derived from these migrations, with most appearing to be hybrids between Bos taurus and Bos indicus. The present study examines the patterns of admixture, diversity, and relationships among African cattle breeds.

Methods

Data for ~ 40 k SNPs was obtained from previous projects for 4089 animals representing 35 African indigenous, 6 European Bos taurus, 4 Bos indicus, and 5 African crossbred cattle populations. Genetic diversity and population structure were assessed using principal component analyses (PCA), admixture analyses, and Wright’s F statistic. The linkage disequilibrium and effective population size (Ne) were estimated for the pure cattle populations.

Results

The first two principal components differentiated Bos indicus from European Bos taurus, and African Bos taurus from other breeds. PCA and admixture analyses showed that, except for recently admixed cattle, all indigenous breeds are either pure African Bos taurus or admixtures of African Bos taurus and Bos indicus. The African zebu breeds had highest proportions of Bos indicus ancestry ranging from 70 to 90% or 60 to 75%, depending on the admixture model. Other indigenous breeds that were not 100% African Bos taurus, ranged from 42 to 70% or 23 to 61% Bos indicus ancestry. The African Bos taurus populations showed substantial genetic diversity, and other indigenous breeds show evidence of having more than one African taurine ancestor. Ne estimates based on r² and r²_adj showed a decline in Ne from a large population at 2000 generations ago, which is surprising for the indigenous breeds given the expected increase in cattle populations over that period and the lack of structured breeding programs.

Conclusion

African indigenous cattle breeds have a large genetic diversity and are either pure African Bos taurus or admixtures of African Bos taurus and Bos indicus. This provides a rich resource of potentially valuable genetic variation, particularly for adaptation traits, and to support conservation programs. It also provides challenges for the development of genomic assays and tools for use in African populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07270-x.

Sequence analysis in Bos taurus reveals pervasiveness of X-Y arms races in mammalian lineages

Studies of Y Chromosome evolution have focused primarily on gene decay, a consequence of suppression of crossing-over with the X Chromosome. Here, we provide evidence that suppression of X-Y crossing-over unleashed a second dynamic: selfish X-Y arms races that reshaped the sex chromosomes in mammals as different as cattle, mice, and men. Using super-resolution sequencing, we explore the Y Chromosome of Bos taurus (bull) and find it to be dominated by massive, lineage-specific amplification of testis-expressed gene families, making it the most gene-dense Y Chromosome sequenced to date. As in mice, an X-linked homolog of a bull Y-amplified gene has become testis-specific and amplified. This evolutionary convergence implies that lineage-specific X-Y coevolution through gene amplification, and the selfish forces underlying this phenomenon, were dominatingly powerful among diverse mammalian lineages. Together with Y gene decay, X-Y arms races molded mammalian sex chromosomes and influenced the course of mammalian evolution.

Coxiella burnetii in Dromedary Camels (Camelus dromedarius): A Possible Threat for Humans and Livestock in North Africa and the Near and Middle East?

The "One Health" concept recognizes that human health is connected to animal health and to the ecosystems. Coxiella burnetii-induced human Q fever is one of the most widespread neglected zoonosis. The main animal reservoirs responsible for C. burnetii transmission to humans are domesticated ruminants, primarily goats, sheep, and cattle. Although studies are still too sparse to draw definitive conclusions, the most recent C. burnetii serosurvey studies conducted in herds and farms in Africa, North Africa, Arabian Peninsula, and Asia highlighted that seroprevalence was strikingly higher in dromedary camels (Camelus dromedarius) than in other ruminants. The C. burnetii seroprevalence in camel herds can reach more than 60% in Egypt, Saudi Arabia, and Sudan, and 70 to 80% in Algeria and Chad, respectively. The highest seroprevalence was in female camels with a previous history of abortion. Moreover, C. burnetii infection was reported in ticks of the Hyalomma dromedarii and Hyalomma impeltatum species collected on camels. Even if dromedary camels represent <3% of the domesticated ruminants in the countries of the Mediterranean basin Southern coast, these animals play a major socioeconomic role for millions of people who live in the arid zones of Africa, Middle East, and Asia. In Chad and Somalia, camels account for about 7 and 21% of domesticated ruminants, respectively. To meet the growing consumers demand of camel meat and milk (>5 million tons/year of both raw and pasteurized milk according to the Food and Agriculture Organization) sustained by a rapid increase of population (growth rate: 2.26-3.76 per year in North Africa), dromedary camel breeding tends to increase from the Maghreb to the Arabic countries. Because of possible long-term persistence of C. burnetii in camel hump adipocytes, this pathogen could represent a threat for herds and breeding farms and ultimately for public health. Because this review highlights a hyperendemia of C. burnetii in dromedary camels, a proper screening of herds and breeding farms for C. burnetii is urgently needed in countries where camel breeding is on the rise. Moreover, the risk of C. burnetii transmission from camel to human should be further evaluated.

The distribution of four missense mutations in SPTBN5 gene across native Chinese breeds

Modern cattle belong to two subspecies, Bos taurus and Bos indicus. Since divergence, cattle types have accumulated different genetic variations, which have contributed to highly differentiated phenotypes. The mammalian inner ear possesses functional and morphological innovations that contribute to its unique hearing capacities. The spectrin beta, non-erythrocytic 5 (SPTBN5) gene has been shown to play an important function in the inner ear. Four missense mutations: rs522333459 (c.7232G > C:p.Cys2411Ser), rs718838405 (c.6568A > C:p.Met2190Leu), rs516536785 (c.6283C > T:p.Leu2095Phe) and rs480278206 (c.4201T > C:p.Cys1401Arg) were identified in the bovine SPTBN5 gene by whole genome resequencing (http://animal.nwsuaf.edu.cn/code/index.php/BosVar), which might be candidate mutations related with hearing of both taurine and indicine cattle. In our study, PCR and DNA sequencing were used to explore the allele frequencies of four mutations of 971 individuals belonging to 38 native Chinese cattle breeds. We find that four mutant alleles showing strong geographic distribution, consisting with the ancestry distribution of taurine and indicine in China. In addition, we identified four mutations of SPTBN5 were diverged in taurine and indicine cattle showing signatures of adaptive evolution in two subspecies, which might participate in bovine inner ear development.

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.

The mosaic genome of indigenous African cattle as a unique genetic resource for African pastoralism

Cattle pastoralism plays a central role in human livelihood in Africa. However, the genetic history of its success remains unknown. Here, through whole-genome sequence analysis of 172 indigenous African cattle from 16 breeds representative of the main cattle groups, we identify a major taurine × indicine cattle admixture event dated to circa 750-1,050 yr ago, which has shaped the genome of today's cattle in the Horn of Africa. We identify 16 loci linked to African environmental adaptations across crossbred animals showing an excess of taurine or indicine ancestry. These include immune-, heat-tolerance- and reproduction-related genes. Moreover, we identify one highly divergent locus in African taurine cattle, which is putatively linked to trypanotolerance and present in crossbred cattle living in trypanosomosis-infested areas. Our findings indicate that a combination of past taurine and recent indicine admixture-derived genetic resources is at the root of the present success of African pastoralism.

Whole genome characterization of autochthonous Bos taurus brachyceros and introduced Bos indicus indicus cattle breeds in Cameroon regarding their adaptive phenotypic traits and pathogen resistance

Background

African indigenous taurine cattle display unique adaptive traits shaped by husbandry management, regional climate and exposure to endemic pathogens. They are less productive with respect to milk and meat production which has been associated with amongst others, small size, traditional beliefs, husbandry practices, limited feed resources, disease burden and lack of sustained breeding for trait improvement. This resulted in the severe dwindling of their population size rendering them vulnerable to extinction.

The Namchi taurine cattle breed is referred to as [Namchi (Doayo)] and shows resistance traits against trypanosome infection and exposure to tick infestation. Nonetheless, the historically later introduced Zebu cattle are the main cattle breeds in Africa today, even though they suffer more from locally prevailing pathogens.

By using a whole genome sequencing approach, we sequenced with high depth for the first time the genomes of five cattle breeds from Cameroon in order to provide a valuable genetic resource for future African cattle breeding: the Namchi, an endangered trypano-tolerant taurine breed, the Kapsiki, an indigenous trypano-susceptible taurine breed, and three Zebu (Bos indicus indicus) breeds: Ngaoundere Gudali, White Fulani and Red Fulani.

Results

Approximately 167 Gigabases of raw sequencing data were generated for each breed and mapped to the cattle reference genomes ARS-UCD1.2 and UMD3.1.The coverage was 103 to 140-fold when aligning the reads to ARS-UCD1.2 with an average mapping rate of ~ 99%, and 22 to 30-fold when aligning the reads to UMD3.1 with an average mapping rate of ~ 64%. The single nucleotide polymorphisms (SNPs) obtained from analysis using the genome ARS-UCD1.2 were compared with reference genomes of European Bos taurus Holstein, the Asian Bos indicus Brahman, and the African trypanotolerant N’Dama breeds.

A total of ~ 100 million (M) SNPs were identified and 7.7 M of those were breed-specific. An approximately 11.1 M constituted of small insertions and deletions. By using only breed-specific non-synonymous variants we identified genes as genetic signatures and associated Gene Ontology (GO) terms that could explain certain cattle-breed specific phenotypes such as increased tolerance against trypanosome parasites in the Namchi breed and heat tolerance in the Kapsiki breed. Phylogenetic analysis grouped, except for Namchi, the Bos taurus breeds Kapsiki, N’Dama and Holstein together while the B. indicus breeds White and Red Fulani, Gudali and Brahman clustered separately. The deviating result for Namchi indicates a hybrid status of the selected animal with a recent introgression of Zebu genes into its genome.

Conclusions

The findings provide the first comprehensive set of genome-wide variant data of the most important Cameroonian cattle breeds. The genomic data shall constitute a foundation for breed amelioration whilst exploiting the heritable traits and support conservation efforts for the endangered local cattle breeds.

Caracterização do gene do choque térmico (HSP-70.1) e sua relação com características de produção em bovinos leiteiros criados no semiárido brasileiro

RESUMO Objetivou-se com este trabalho avaliar a diversidade genética do gene HSP-70.1 e associar os polimorfismos encontrados com a performance de vacas leiteiras das raças Holandesa, Girolando (5/8H-G) e Sindi criadas em região do semiárido brasileiro. Os polimorfismos foram identificados e avaliados pela técnica de PCR-RFLP, usando-se a enzima de restrição EcoRII. Avaliou-se a variabilidade genética por meio do índice de diversidade padrão e da análise de variância molecular (AMOVA). Os polimorfismos identificados foram avaliados sobre as características de produção de leite. Foram identificados sete alelos, os quais demonstraram que houve polimorfismo para a região gênica analisada, e alguns alelos foram compartilhados entre os rebanhos. As raças bovinas Holandesa e Sindi foram similares geneticamente para o gene analisado. A AMOVA demonstrou que há variação genética entre os rebanhos e dentro deles, com a maior parte da variação ocorrendo dentro dos rebanhos para todos os grupos avaliados. Houve efeito dos alelos identificados sobre a produção de leite dos rebanhos das raças Holandesa (P<0,0001) e Girolando (P<0,0117). O gene HSP-70.1 foi polimórfico na população de bovinos leiteiros estudada, sendo, portanto, um marcador molecular promissor para avaliar a produção de leite de raças criadas em região semiárida.

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

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

Population genetic structure and maternal lineage of South African crossbred Nguni cattle using the cytochrome b gene in mtDNA

The Nguni cattle breed predominates South Africa and is endowed with traits favourable against environmental stressors such as heat stress and resistance to diseases. Interventions to improve production have led to the erosion of the genetic integrity of local breeds and the introduction of exotic breeds has proved ineffective as they fail to perform well due to different climatic conditions and production systems. In this study, the genetic structure and genetic lineage of Nguni crossbreds from 6 populations were assessed using the mitochondrial cytochrome b gene. Twelve polymorphic sites were detected resulting in 11 haplotypes with haplotype and nucleotide diversities of 0.550 ± 0.135 and 0.0019 ± 0.0011, respectively. Only 2 of the 6 populations displayed recent population expansion events, whereas the majority adhered to neutral evolution. The basal haplotype contained approximately 60% of the studied populations and there were four unique haplotypes that were revealed. A possible Nguni descript haplotype was uncovered, and this haplotype was found in all populations but was however devoid of individuals from around the world. The genetic structure of the populations was rather low (average pairwise F_ST = 0.066 and Slatkins F_ST = 0.094), and approximately 96% of the total genetic variation was accounted for by differences within populations. Phylogenetic analyses supported the clustering of all the samples within the Bos taurus clade and no Bos indicus haplotype was detected. Furthermore, no intermediate haplotype of taurine and indicine was detected. Overall, the maternal lineage of the crossbreds points to a taurine origin and the low genetic diversity depicts the retention of the Nguni genetic pool and possibly its superior adaptive traits.

Unlocking the origins and biology of domestic animals using ancient DNA and paleogenomics

Animal domestication has fascinated biologists since Charles Darwin first drew the parallel between evolution via natural selection and human-mediated breeding of livestock and companion animals. In this review we show how studies of ancient DNA from domestic animals and their wild progenitors and congeners have shed new light on the genetic origins of domesticates, and on the process of domestication itself. High-resolution paleogenomic data sets now provide unprecedented opportunities to explore the development of animal agriculture across the world. In addition, functional population genomics studies of domestic and wild animals can deliver comparative information useful for understanding recent human evolution.