The Origins of Domesticated Cattle

Genome-wide mining of diversity and evolutionary signatures revealed selective hotspots in Indian Sahiwal cattle

The Sahiwal cattle breed is the best indigenous dairy cattle breed, and it plays a pivotal role in the Indian dairy industry. This is due to its exceptional milk-producing potential, adaptability to local tropical conditions, and its resilience to ticks and diseases. The study aimed to identify selective sweeps and estimate intrapopulation genetic diversity parameters in Sahiwal cattle using ddRAD sequencing-based genotyping data from 82 individuals. After applying filtering criteria, 78,193 high-quality SNPs remained for further analysis. The population exhibited an average minor allele frequency of 0.221 ± 0.119. Genetic diversity metrics, including observed (0.597 ± 0.196) and expected heterozygosity (0.433 ± 0.096), nucleotide diversity (0.327 ± 0.114), the proportion of polymorphic SNPs (0.726), and allelic richness (1.323 ± 0.134), indicated ample genomic diversity within the breed. Furthermore, an effective population size of 74 was observed in the most recent generation. The overall mean linkage disequilibrium (r2) for pairwise SNPs was 0.269 ± 0.057. Moreover, a greater proportion of short Runs of Homozygosity (ROH) segments were observed suggesting that there may be low levels of recent inbreeding in this population. The genomic inbreeding coefficients, computed using different inbreeding estimates (FHOM, FUNI, FROH, and FGROM), ranged from -0.0289 to 0.0725. Subsequently, we found 146 regions undergoing selective sweeps using five distinct statistical tests: Tajima's D, CLR, |iHS|, |iHH12|, and ROH. These regions, located in non-overlapping 500 kb windows, were mapped and revealed various protein-coding genes associated with enhanced immune systems and disease resistance (IFNL3, IRF8, BLK), as well as production traits (NRXN1, PLCE1, GHR). Notably, we identified interleukin 2 (IL2) on Chr17: 35217075-35223276 as a gene linked to tick resistance and uncovered a cluster of genes (HSPA8, UBASH3B, ADAMTS18, CRTAM) associated with heat stress. These findings indicate the evolutionary impact of natural and artificial selection on the environmental adaptation of the Sahiwal cattle population.

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.

A phylogenomic analysis of Limosilactobacillus reuteri reveals ancient and stable evolutionary relationships with rodents and birds and zoonotic transmission to humans

BACKGROUND

Gut microbes play crucial roles in the development and health of their animal hosts. However, the evolutionary relationships of gut microbes with vertebrate hosts, and the consequences that arise for the ecology and lifestyle of the microbes are still insufficiently understood. Specifically, the mechanisms by which strain-level diversity evolved, the degree by which lineages remain stably associated with hosts, and how their evolutionary history influences their ecological performance remain a critical gap in our understanding of vertebrate-microbe symbiosis.

RESULTS

This study presents the characterization of an extended collection of strains of Limosilactobacillus reuteri and closely related species from a wide variety of hosts by phylogenomic and comparative genomic analyses combined with colonization experiments in mice to gain insight into the long-term evolutionary relationship of a bacterial symbiont with vertebrates. The phylogenetic analysis of L. reuteri revealed early-branching lineages that primarily consist of isolates from rodents (four lineages) and birds (one lineage), while lineages dominated by strains from herbivores, humans, pigs, and primates arose more recently and were less host specific. Strains from rodent lineages, despite their phylogenetic divergence, showed tight clustering in gene-content-based analyses. These L. reuteri strains but not those ones from non-rodent lineages efficiently colonize the forestomach epithelium of germ-free mice. The findings support a long-term evolutionary relationships of L. reuteri lineages with rodents and a stable host switch to birds. Associations of L. reuteri with other host species are likely more dynamic and transient. Interestingly, human isolates of L. reuteri cluster phylogenetically closely with strains from domesticated animals, such as chickens and herbivores, suggesting zoonotic transmissions.

CONCLUSIONS

Overall, this study demonstrates that the evolutionary relationship of a vertebrate gut symbiont can be stable in particular hosts over time scales that allow major adaptations and specialization, but also emphasizes the diversity of symbiont lifestyles even within a single bacterial species. For L. reuteri, symbiont lifestyles ranged from autochthonous, likely based on vertical transmission and stably aligned to rodents and birds over evolutionary time, to allochthonous possibly reliant on zoonotic transmission in humans. Such information contributes to our ability to use these microbes in microbial-based therapeutics.

Ancient Faunal History Revealed by Interdisciplinary Biomolecular Approaches

Starting four decades ago, studies have examined the ecology and evolutionary dynamics of populations and species using short mitochondrial DNA fragments and stable isotopes. Through technological and analytical advances, the methods and biomolecules at our disposal have increased significantly to now include lipids, whole genomes, proteomes, and even epigenomes. At an unprecedented resolution, the study of ancient biomolecules has made it possible for us to disentangle the complex processes that shaped the ancient faunal diversity across millennia, with the potential to aid in implicating probable causes of species extinction and how humans impacted the genetics and ecology of wild and domestic species. However, even now, few studies explore interdisciplinary biomolecular approaches to reveal ancient faunal diversity dynamics in relation to environmental and anthropogenic impact. This review will approach how biomolecules have been implemented in a broad variety of topics and species, from the extinct Pleistocene megafauna to ancient wild and domestic stocks, as well as how their future use has the potential to offer an enhanced understanding of drivers of past faunal diversity on Earth.

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

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

Genomic scans for selection signatures revealed candidate genes for adaptation and production traits in a variety of cattle breeds

Domestication and selection are the major driving forces responsible for the determinative genetic variability in livestock. These selection patterns create unique genetic signatures within the genome. BovineSNP50 chip data from 236 animals (seven indicine and five taurine cattle breeds) were analyzed in the present study. We implemented three complementary approaches viz. iHS (Integrated haplotype score), ROH (Runs of homozygosity), and F_ST, to detect selection signatures. A total of 179, 56, and 231 regions revealed 518, 277, and 267 candidate genes identified by iHS, ROH, and F_ST methods, respectively. We found several candidate genes (e.g., NCR3, ARID5A, HIST1H2BN, DEFB4, DEFB7, HSPA1L, HSPA1B, and DNAJB4) related to production traits and the adaptation of indigenous breeds to local environmental constraints such as heat stress and disease susceptibility. However, further studies are warranted to refine the findings using a larger sample size, whole-genome sequencing, and/or high density genotyping.

Evolution and Biogeography of Haemonchus contortus: Linking Faunal Dynamics in Space and Time

History is the foundation that informs about the nuances of faunal assembly that are essential in understanding the dynamic nature of the host-parasite interface. All of our knowledge begins and ends with evolution, ecology and biogeography, as these interacting facets determine the history of biodiverse systems. These components, relating to Haemonchus, can inform about the complex history of geographical distribution, host association and the intricacies of host-parasite associations that are played out in physiological and behavioural processes that influence the potential for disease and our capacity for effective control in a rapidly changing world. Origins and evolutionary diversification among species of the genus Haemonchus and Haemonchus contortus occurred in a complex crucible defined by shifts in environmental structure emerging from cycles of climate change and ecological perturbation during the late Tertiary and through the Quaternary. A history of sequential host colonization associated with waves of dispersal bringing assemblages of ungulates from Eurasia into Africa and processes emerging from ecosystems in collision and faunal turnover defined the arena for radiation among 12 recognized species of Haemonchus. Among congeners, the host range for H. contortus is exceptionally broad, including species among artiodactyls of 40 genera representing 5 families (and within 12 tribes of Bovidae). Broad host range is dramatically reflected in the degree to which translocation, introduction and invasion with host switching, has characterized an expanding distribution over time in North America, South America, southern Eurasia, Australia and New Zealand, coincidental with agriculture, husbandry and global colonization by human populations driven particularly by European exploration after the 1500s. African origins in xeric to mesic habitats of the African savannah suggest that historical constraints linked to ecological adaptations (tolerances and developmental thresholds defined by temperature and humidity for larval stages) will be substantial determinants in the potential outcomes for widespread geographical and host colonization which are predicted to unfold over the coming century. Insights about deeper evolutionary events, ecology and biogeography are critical as understanding history informs us about the possible range of responses in complex systems under new regimes of environmental forcing, especially, in this case, ecological perturbation linked to climate change. A deeper history of perturbation is relevant in understanding contemporary systems that are now strongly structured by events of invasion and colonization. The relaxation of abiotic and biotic controls on the occurrence of H. contortus, coincidental with inception and dissemination of anthelmintic resistance may be synergistic, serving to exacerbate challenges to control parasites or to limit the socioeconomic impacts of infection that can influence food security and availability. Studies of haemonchine nematodes contribute directly to an expanding model about the nature of diversity and the evolutionary trajectories for faunal assembly among complex host-parasite systems across considerable spatial and temporal scales.

Effects of anthelmintic treatment and feed supplementation on parasite infections and morbidity parameters in Cambodian cattle

Helminth infections are the cause of morbidity in Cambodian cattle but other factors such as nutritional deficiencies and concurrent diseases may enhance the effects of parasites. The present study aimed to investigate the impact of anthelmintic treatment, feed supplementation, or both on gastrointestinal strongyle (GIS) and trematode infections as well as on morbidity parameters in Cambodian village cattle. At the beginning of the dry season, cattle populations in six villages were randomly assigned to a group: (A) receiving anthelmintic treatment (ivermectin+clorsulon) at week 0; (P) feed pellet supplementation during week 0-13 or both (AP). On five visits (week 0-29), faecal and blood samples were obtained for parasitological examination and haematocrit determination, respectively. Body condition (BCS), hind quarter fouling (HQFS), diarrhoea (DS), and conjunctiva colour (FAMACHA^©) were scored and heart girth circumference was determined. To investigate the impact of treatment over time (week 0-29), a mixed model was used with treatment, time, and their interaction as fixed effects, and animal and village as random factors. At baseline, the proportion of GIS positive animals was high (67.9%), whereas trematode infections were low (Paramphistomum: 8.8%; Fasciola: 2.6%). Very thin to emaciated cattle (BCS 1-2) were more prevalent (11.4%) and FAMACHA^© scores of ≤3 or below (65.8%) less prevalent than in an earlier study in the region. A Time ⨯ Treatment interaction was present for faecal egg counts (FEC) of GIS, GIS prevalence (both p<0.0001), PCV (p=0.0034), DS (p=0.0086) and HQFS (p=0.0241). For GIS FEC, treatment groups differed at a specific time point, with levels of treatment group P being higher than in A at week 6 (p=0.0054). For Paramphistomum prevalence as well as FAMACHA^© scoring, heart girth and BCS, the interaction between treatment and time was not significant, yet, time in itself had a significant impact on all (p<0.0001). The beneficial effects of protein supplementation were unclear from the current study.

Trace element distribution in selected edible tissues of zebu (Bos indicus) cattle slaughtered at Jimma, SW Ethiopia

The amount of trace elements present in edible bovine tissues is of importance for both animal health and human nutrition. This study presents data on trace element concentrations in semitendinosus and cardiac muscles, livers and kidneys of 60 zebu (Bos indicus) bulls, sampled at Jimma, Ethiopia. From 28 of these bulls, blood samples were also obtained. Deficient levels of copper were found in plasma, livers, kidneys and semitendinosus muscles. Suboptimal selenium concentrations were found in plasma and semitendinosus muscles. Semitendinosus muscles contained high iron concentrations. Trace elements were mainly stored in the liver, except for iron and selenium. Cardiac muscles generally contained higher concentrations of trace elements than semitendinous muscles except for zinc. A strong association was found between liver and kidney concentrations of copper, iron, cobalt and molybdenum. Liver storage was well correlated with storage in semitendinosus muscle for selenium and with cardiac muscle for cobalt and selenium. Plasma concentrations of copper, selenium, cobalt were well related with their respective liver concentrations and for cobalt and selenium, also with cardiac muscle concentrations. The data suggest multiple trace element deficiencies in zebu cattle in South-West Ethiopia, with lowered tissue concentrations as a consequence. Based on the comparison of our data with other literature, trace element concentrations in selected edible tissues of Bos indicus seem quite similar to those in Bos taurus. However, tissue threshold values for deficiency in Bos taurus cattle need to be refined and their applicability for Bos indicus cattle needs to be evaluated.

The genetic diversity of the Nguni breed of African Cattle (Bos spp.): complete mitochondrial genomes of haplogroup T1

Domesticated cattle were commonplace in northern Africa by about 7,000 years ago. Archaeological evidence, however, suggests they were not established in southern Africa until much later, no earlier than 2,000 years ago. Genetic reconstructions have started to shed light on the movement of African cattle, but efforts have been frustrated by a lack of data south of Ethiopia and the nature of the mitochondrial haplogroup T1 which is almost fixed across the continent. We sequenced 35 complete mitochondrial genomes from a South African herd of Nguni cattle, a breed historically associated with Bantu speaking farmers who were among the first to bring cattle to southern Africa. As expected, all individuals in the study were found to be members of haplogroup T1. Only half of the sub-haplogroups of T1 (T1a-T1f) are represented in our sample and the overwhelming majority (94%) in this study belong to subhaplogroup T1b. A previous study of African cattle found frequencies of T1b of 27% in Egypt and 69% in Ethiopia. These results are consistent with serial multiple founder effects significantly shaping the gene pool as cattle were moved from north to south across the continent. Interestingly, these mitochondrial data give no indication that the impacts of the founder effects were ameliorated by gene flow from recently introduced Indian cattle breeds.

Origin and spread of Bos taurus: new clues from mitochondrial genomes belonging to haplogroup T1

Background

Most genetic studies on modern cattle have established a common origin for all taurine breeds in the Near East, during the Neolithic transition about 10 thousand years (ka) ago. Yet, the possibility of independent and/or secondary domestication events is still debated and is fostered by the finding of rare mitochondrial DNA (mtDNA) haplogroups like P, Q and R. Haplogroup T1, because of its geographic distribution, has been the subject of several investigations pointing to a possible independent domestication event in Africa and suggesting a genetic contribution of African cattle to the formation of Iberian and Creole cattle. Whole mitochondrial genome sequence analysis, with its proven effectiveness in improving the resolution of phylogeographic studies, is the most appropriate tool to investigate the origin and structure of haplogroup T1.

Methodology

A survey of >2200 bovine mtDNA control regions representing 28 breeds (15 European, 10 African, 3 American) identified 281 subjects belonging to haplogroup T1. Fifty-four were selected for whole mtDNA genome sequencing, and combined with ten T1 complete sequences from previous studies into the most detailed T1 phylogenetic tree available to date.

Conclusions

Phylogenetic analysis of the 64 T1 mitochondrial complete genomes revealed six distinct sub-haplogroups (T1a–T1f). Our data support the overall scenario of a Near Eastern origin of the T1 sub-haplogroups from as much as eight founding T1 haplotypes. However, the possibility that one sub-haplogroup (T1d) arose in North Africa, in domesticated stocks, shortly after their arrival from the Near East, can not be ruled out. Finally, the previously identified “African-derived American" (AA) haplotype turned out to be a sub-clade of T1c (T1c1a1). This haplotype was found here for the first time in Africa (Egypt), indicating that it probably originated in North Africa, reached the Iberian Peninsula and sailed to America, with the first European settlers.

Genome sequence and assembly of Bos indicus

Cattle are divided into 2 groups referred to as taurine and indicine, both of which have been under strong artificial selection due to their importance for human nutrition. A side effect of this domestication includes a loss of genetic diversity within each specialized breed. Recently, the first taurine genome was sequenced and assembled, allowing for a better understanding of this ruminant species. However, genetic information from indicine breeds has been limited. Here, we present the first genome sequence of an indicine breed (Nellore) generated with 52X coverage by SOLiD sequencing platform. As expected, both genomes share high similarity at the nucleotide level for all autosomes and the X chromosome. Regarding the Y chromosome, the homology was considerably lower, most likely due to uncompleted assembly of the taurine Y chromosome. We were also able to cover 97% of the annotated taurine protein-coding genes.

Genomic characteristics of cattle copy number variations

Background

Copy number variation (CNV) represents another important source of genetic variation complementary to single nucleotide polymorphism (SNP). High-density SNP array data have been routinely used to detect human CNVs, many of which have significant functional effects on gene expression and human diseases. In the dairy industry, a large quantity of SNP genotyping results are becoming available and can be used for CNV discovery to understand and accelerate genetic improvement for complex traits.

Results

We performed a systematic analysis of CNV using the Bovine HapMap SNP genotyping data, including 539 animals of 21 modern cattle breeds and 6 outgroups. After correcting genomic waves and considering the pedigree information, we identified 682 candidate CNV regions, which represent 139.8 megabases (~4.60%) of the genome. Selected CNVs were further experimentally validated and we found that copy number "gain" CNVs were predominantly clustered in tandem rather than existing as interspersed duplications. Many CNV regions (~56%) overlap with cattle genes (1,263), which are significantly enriched for immunity, lactation, reproduction and rumination. The overlap of this new dataset and other published CNV studies was less than 40%; however, our discovery of large, high frequency (> 5% of animals surveyed) CNV regions showed 90% agreement with other studies. These results highlight the differences and commonalities between technical platforms.

Conclusions

We present a comprehensive genomic analysis of cattle CNVs derived from SNP data which will be a valuable genomic variation resource. Combined with SNP detection assays, gene-containing CNV regions may help identify genes undergoing artificial selection in domesticated animals.

Analysis of copy number variations among diverse cattle breeds

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. Here, we describe the first systematic and genome-wide analysis of copy number variations (CNVs) in modern domesticated cattle using array comparative genomic hybridization (array CGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH). The array CGH panel included 90 animals from 11 Bos taurus, three Bos indicus, and three composite breeds for beef, dairy, or dual purpose. We identified over 200 candidate CNV regions (CNVRs) in total and 177 within known chromosomes, which harbor or are adjacent to gains or losses. These 177 high-confidence CNVRs cover 28.1 megabases or approximately 1.07% of the genome. Over 50% of the CNVRs (89/177) were found in multiple animals or breeds and analysis revealed breed-specific frequency differences and reflected aspects of the known ancestry of these cattle breeds. Selected CNVs were further validated by independent methods using qPCR and FISH. Approximately 67% of the CNVRs (119/177) completely or partially span cattle genes and 61% of the CNVRs (108/177) directly overlap with segmental duplications. The CNVRs span about 400 annotated cattle genes that are significantly enriched for specific biological functions, such as immunity, lactation, reproduction, and rumination. Multiple gene families, including ULBP, have gone through ruminant lineage-specific gene amplification. We detected and confirmed marked differences in their CNV frequencies across diverse breeds, indicating that some cattle CNVs are likely to arise independently in breeds and contribute to breed differences. Our results provide a valuable resource beyond microsatellites and single nucleotide polymorphisms to explore the full dimension of genetic variability for future cattle genomic research.