Ancient genomes reveal tropical bovid species in the Tibetan Plateau contributed to the prevalence of hunting game until the late Neolithic.

Local wild bovids have been determined to be important prey on the northeastern Tibetan Plateau (NETP), where hunting game was a major subsistence strategy until the late Neolithic, when farming lifestyles dominated in the neighboring Loess Plateau. However, the species affiliation and population ecology of these prehistoric wild bovids in the prehistoric NETP remain unknown. Ancient DNA (aDNA) analysis is highly informative in decoding this puzzle. Here, we applied aDNA analysis to fragmented bovid and rhinoceros specimens dating ∼5,200 y B.P. from the Neolithic site of Shannashuzha located in the marginal area of the NETP. Utilizing both whole genomes and mitochondrial DNA, our results demonstrate that the range of the present-day tropical gaur (Bos gaurus) extended as far north as the margins of the NETP during the late Neolithic from ∼29°N to ∼34°N. Furthermore, comparative analysis with zooarchaeological and paleoclimatic evidence indicated that a high summer temperature in the late Neolithic might have facilitated the northward expansion of tropical animals (at least gaur and Sumatran-like rhinoceros) to the NETP. This enriched the diversity of wildlife, thus providing abundant hunting resources for humans and facilitating the exploration of the Tibetan Plateau as one of the last habitats for hunting game in East Asia.

Development of a multivariable prediction model to identify dairy calves too young to be transported to auction markets in Canada using simple physical examination and body weight.

Calves born on Eastern Canadian dairy farms that are not kept in the herds are traditionally sold through auction markets and are raised for meat purposes such as veal calves. Since February 2020, a new Canadian federal regulation has forbidden calves <9 d old to be sold through auction markets. However, in the absence of a real-time birth registry consultation system, it would be of interest to look for predictors that could be associated with age to allow identification of calves too young to be transported. In the current retrospective cross-sectional study, 1,178 calves with a declared birth date (411 calves aged <9 d old; 34.9%) were assessed in 2 large Québec auction sites. Easy-to-record covariates [body weight (BW), breed phenotype, and presence of an umbilical cord remnant] as well as other clinical signs (umbilical swelling, enlargement, umbilical pain, wet umbilicus, skin tent, sunken eyes, ocular and nasal secretion, and hide cleanliness) were assessed. Two logistic regression models using age as a dichotomous dependent variable (<9 d old vs ≥9 d old) were built. The first model (model 1) considered all covariates, which were selected after univariable analyses and a backward stepwise selection process, whereas a more pragmatic model (model 2) only included the 3 easy-to-record variables (i.e., BW, breed, umbilical cord). Both models had similar accuracy to detect calves <9 d old (sensitivity of 38.4 and 37.5%, and specificity of 85.7 and 84.6% for model 1 and 2, respectively). Model 2 was subsequently more specifically studied as it employs a faster and easier assessment. Decision thresholds were tested for their robustness based on misclassification cost term (MCT) analysis with various prevalence of calves <9 d old and various costs of false-negative:false-positive ratio. Despite statistical significance, model accuracy, even if refined with MCT analysis, was limited at the individual level, showing the limits of using physical signs and BW or their combination as a reliable proxy of age. The sensitivity of these models to find calves <9 d old was not to be used for monitoring compliance with the Canadian federal regulation. The relatively high model specificity may help to use this model as a rule-in test (i.e., targeting positive calves for further investigation) rather than a rule-out test (due to its low sensitivity).

Molecular prevalence and risk factors associated with tick-borne pathogens in cattle in western Kenya.

BACKGROUND: Tick-borne pathogens (TBPs) are of global importance, especially in sub-Saharan Africa where they represent a major constraint to livestock production. Their association with human disease is also increasingly recognized, signalling their zoonotic importance. It is therefore crucial to investigate TBPs prevalence in livestock populations and the factors associated with their presence. We set out to identify TBPs present in cattle and to determine associated risk factors in western Kenya, where smallholder livestock production is important for subsistence and market-driven income. RESULTS: Tick-borne pathogen infections in blood samples collected from cattle at livestock markets and slaughterhouses between May 2017 and January 2019 were identified by high-resolution melting analysis and sequencing of PCR products of genus-specific primers. Of the 422 cattle sampled, 30.1% (127/422) were infected with at least one TBP, while 8.8% (37/422) had dual infections. Anaplasma spp. (19.7%) were the most prevalent, followed by Theileria (12.3%), Ehrlichia (6.6%), and Babesia (0.2%) spp. Sequence analysis of the TBPs revealed them to be Anaplasma platys-like organisms (13.5%), Theileria velifera (7.4%), Anaplasma marginale (4.9%), Theileria mutans (3.1%), Theileria parva (1.6%), and Babesia bigemina (0.2%). Ehrlichia ruminantium, Rickettsia spp., and arboviruses were not detected. Exotic breeds of cattle were more likely to be infected with A. marginale compared to local breeds (OR: 7.99, 95% CI: 3.04-22.02, p <  0.001). Presence of ticks was a significant predictor for Anaplasma spp. (OR: 2.18, 95% CI: 1.32-3.69, p = 0.003) and Ehrlichia spp. (OR: 2.79, 95% CI: 1.22-7.23, p = 0.022) infection. Cattle sampled at slaughterhouses were more likely to be positive for Anaplasma spp. (OR: 1.64, 95% CI: 1.01-2.70, p = 0.048) and A. marginale (OR: 3.84, 95% CI: 1.43-12.21, p = 0.012), compared to those sampled at livestock markets. CONCLUSION: This study reports TBP prevalence and associated risk factors in western Kenya, factors which are key to informing surveillance and control measures.

Host genetics and the rumen microbiome jointly associate with methane emissions in dairy cows.

Cattle and other ruminants produce large quantities of methane (~110 million metric tonnes per annum), which is a potent greenhouse gas affecting global climate change. Methane (CH4) is a natural by-product of gastro-enteric microbial fermentation of feedstuffs in the rumen and contributes to 6% of total CH4 emissions from anthropogenic-related sources. The extent to which the host genome and rumen microbiome influence CH4 emission is not yet well known. This study confirms individual variation in CH4 production was influenced by individual host (cow) genotype, as well as the host's rumen microbiome composition. Abundance of a small proportion of bacteria and archaea taxa were influenced to a limited extent by the host's genotype and certain taxa were associated with CH4 emissions. However, the cumulative effect of all bacteria and archaea on CH4 production was 13%, the host genetics (heritability) was 21% and the two are largely independent. This study demonstrates variation in CH4 emission is likely not modulated through cow genetic effects on the rumen microbiome. Therefore, the rumen microbiome and cow genome could be targeted independently, by breeding low methane-emitting cows and in parallel, by investigating possible strategies that target changes in the rumen microbiome to reduce CH4 emissions in the cattle industry.

Comparative analysis of five different methods to design a breed-specific SNP panel for cattle.

Single nucleotide polymorphisms (SNPs) have now replaced microsatellite markers in several species for various genetic investigations like parentage assignment, genetic breed composition, assessment for individuality and, most popularly, as a useful tool in genomic selection. However, such a resource, which can offer to assist breed identification in a cost-effective manner is still not explored in cattle breeding programs. In our study, we have tried to describe methods for reducing the number of SNPs to develop a breed-specific panel. We have used SNP data from Dryad open public access repository. Starting from a global dataset of 178 animals belonging to 10 different breeds, we selected five panels each comprising of similar number of SNPs using different methods i.e., Delta, Pairwise Wright's FST, informativeness for assignment, frequent item feature selection (FIFS) and minor allele frequency-linkage disequilibrium (MAF-LD) based method. MAF-LD based method has been recently developed by us for construction of breed-specific SNP panels. The STRUCTURE software analysis of MAF-LD based method showed appropriate clustering in comparison to other panels. Later, the panel of 591 breed-specific SNPs was called to their respective breeds using Venny 2.1.0 and UGent web tools software. Breed-specific SNPs were later annotated by using various Bioinformatics softwares.

Single nucleotide polymorphisms of leptin gene in five Ethiopian indigenous cattle breeds and the Korean Hanwoo breed.

Considering the escalating number of scientific reports on the association between the leptin gene and diverse physiological traits and performance of cattle populations, this study was directed towards identifying SNPs in the leptin gene among five indigenous cattle breeds of Ethiopia. DNA samples were extracted from the nasal swabs of the Ethiopian indigenous cattle breeds: Arsi (n = 18), Horro (n = 20), Begait (n = 21), Boran (n = 19), and Fogera (n = 17) and the Korean Hanwoo (a representative taurine breed) (n = 20), followed by PCR amplification of exon 2 and exon 3 regions of the leptin gene and sequence analysis of the PCR products. Five SNPs, two (generating missense mutations) on exon 2 and three (generating silent mutations) on exon 3 regions, were explicated in this study. Allele frequency and genotype frequency distribution pertaining to the SNPs were recorded for the studied cattle breeds besides the minor allele frequency and deviation from the Hardy-Weinberg equilibrium. Positive FIS index values were recorded for all the markers except SNP2, illustrative of heterozygote deficiency. MEGA X software-based evolutionary divergence analysis of the phylogenetic tree based on the SNP data revealed that the large-sized breeds, Hanwoo, Begait, Boran, and Fogera, were more closely clustered compared to the small-sized Arsi breed. Among the seven haplotypes documented from the various breeds, sequence analysis was suggestive of haplotypes 1 and 2 to be ancestral haplotypes for the leptin gene. This study is envisaged to accelerate molecular breeding programs for the genetic improvement of the Ethiopian cattle breeds.

Comparative analyses of copy number variations between Bos taurus and Bos indicus.

BACKGROUND: Bos taurus and Bos indicus are two main sub-species of cattle. However, the differential copy number variations (CNVs) between them are not yet well studied. RESULTS: Based on the new high-quality cattle reference genome ARS-UCD1.2, we identified 13,234 non-redundant CNV regions (CNVRs) from 73 animals of 10 cattle breeds (4 Bos taurus and 6 Bos indicus), by integrating three detection strategies. While 6990 CNVRs (52.82%) were shared by Bos taurus and Bos indicus, large CNV differences were discovered between them and these differences could be used to successfully separate animals into two subspecies. We found that 2212 and 538 genes uniquely overlapped with either indicine-specific CNVRs and or taurine-specific CNVRs, respectively. Based on FST, we detected 16 candidate lineage-differential CNV segments (top 0.1%) under selection, which overlapped with eight genes (CTNNA1, ENSBTAG00000004415, PKN2, BMPER, PDE1C, DNAJC18, MUSK, and PLCXD3). Moreover, we obtained 1.74 Mbp indicine-specific sequences, which could only be mapped on the Bos indicus reference genome UOA_Brahman_1. We found these sequences and their associated genes were related to heat resistance, lipid and ATP metabolic process, and muscle development under selection. We further analyzed and validated the top significant lineage-differential CNV. This CNV overlapped genes related to muscle cell differentiation, which might be generated from a retropseudogene of CTH but was deleted along Bos indicus lineage. CONCLUSIONS: This study presents a genome wide CNV comparison between Bos taurus and Bos indicus. It supplied essential genome diversity information for understanding of adaptation and phenotype differences between the Bos taurus and Bos indicus populations.

The genetic heritage of Alpine local cattle breeds using genomic SNP data.

BACKGROUND: Assessment of genetic diversity and population structure provides important control metrics to avoid genetic erosion, inbreeding depression and crossbreeding between exotic and locally-adapted cattle breeds since these events can have deleterious consequences and eventually lead to extinction. Historically, the Alpine Arc represents an important pocket of cattle biodiversity with a large number of autochthonous breeds that provide a fundamental source of income for the entire regional economy. By using genotype data from medium-density single nucleotide polymorphism (SNP) arrays, we performed a genome-wide comparative study of 23 cattle populations from the Alpine Arc and three cosmopolitan breeds. RESULTS: After filtering, we obtained a final genotyping dataset consisting of 30,176 SNPs for 711 individuals. The local breeds showed high or intermediate values of genetic diversity compared to the highly selected cosmopolitan breeds. Patterns of genetic differentiation, multidimensional scaling, admixture analysis and the constructed phylogenetic tree showed convergence, which indicates the presence of gene flow among the breeds according to both geographic origin and historical background. Among the most differentiated breeds, we identified the modern Brown cattle. In spite of admixture events, several local breeds have preserved distinctive characteristics, which is probably due to differences in genetic origin and geographic location. CONCLUSIONS: This study represents one of the most comprehensive genome-wide analysis of the Alpine cattle breeds to date. Using such a large dataset that includes the majority of the local breeds found in this region, allowed us to expand knowledge on the evaluation and status of Alpine cattle biodiversity. Our results indicate that although many of the analyzed local breeds are listed as endangered, they still harbor a large amount of genetic diversity, even when compared to some cosmopolitan breeds. This finding, together with the reconstruction of the phylogeny and the relationships between these Alpine Arc cattle breeds, provide crucial insights not only into the improvement of genetic stocks but also into the implementation of future conservation strategies.

Genomic regions associated with muscularity in beef cattle differ in five contrasting cattle breeds.

BACKGROUND: Linear type traits, which reflect the muscular characteristics of an animal, could provide insight into how, in some cases, morphologically very different animals can yield the same carcass weight. Such variability may contribute to differences in the overall value of the carcass since primal cuts vary greatly in price; such variability may also hinder successful genome-based association studies. Therefore, the objective of our study was to identify genomic regions that are associated with five muscularity linear type traits and to determine if these significant regions are common across five different breeds. Analyses were carried out using linear mixed models on imputed whole-genome sequence data in each of the five breeds, separately. Then, the results of the within-breed analyses were used to conduct an across-breed meta-analysis per trait. RESULTS: We identified many quantitative trait loci (QTL) that are located across the whole genome and associated with each trait in each breed. The only commonality among the breeds and traits was a large-effect pleiotropic QTL on BTA2 that contained the MSTN gene, which was associated with all traits in the Charolais and Limousin breeds. Other plausible candidate genes were identified for muscularity traits including PDE1A, PPP1R1C and multiple collagen and HOXD genes. In addition, associated (gene ontology) GO terms and KEGG pathways tended to differ between breeds and between traits especially in the numerically smaller populations of Angus, Hereford, and Simmental breeds. Most of the SNPs that were associated with any of the traits were intergenic or intronic SNPs located within regulatory regions of the genome. CONCLUSIONS: The commonality between the Charolais and Limousin breeds indicates that the genetic architecture of the muscularity traits may be similar in these breeds due to their similar origins. Conversely, there were vast differences in the QTL associated with muscularity in Angus, Hereford, and Simmental. Knowledge of these differences in genetic architecture between breeds is useful to develop accurate genomic prediction equations that can operate effectively across breeds. Overall, the associated QTL differed according to trait, which suggests that breeding for a morphologically different (e.g. longer and wider versus shorter and smaller) more efficient animal may become possible in the future.

Evaluation of two bovine SNP genotyping arrays for breed clustering and stratification analysis in well-known taurine and indicine breeds.

The present study aimed to evaluate the efficiency of two Bovine SNP genotyping arrays (i.e., 50 K and HD) for breed clustering and stratification related studies in taurine and indicine breeds. The whole-genome SNP data at two densities were assembled into three datasets (A, B and C). Dataset A (N = 213) included 50 K genotypic data for five taurine (Holstein-Friesian, Guernsey, Brown Swiss, Angus and Jersey) and two indicine (Gir and Nellore) breeds. Dataset B (N = 241) included the same breeds with HD density data. Dataset C (N = 299) included 50 K SNP genotypic data for six taurine (Holstein-Friesian, Jersey, Guernsey, Brown Swiss, Angus and Hereford) and six indicine (Hariana, Kankrej, Brahman, Nellore, Sahiwal and Gir) breeds. The analysis was done using ADMIXTURE program (bioinformatics-based) and cross-validation errors and Principal Component Analysis (statistical analysis). The proportion of polymorphic markers and minor allele frequencies were assessed for each breed. The proportion of markers polymorphic was consistently higher in taurine breeds when compared with breeds from indicine group. Minor allele frequency estimates and ADMIXTURE results showed differential patterns for both the lineages. However, no significant increase in the accuracy of genomic clustering was found on moving from 50 K to HD density data.

Breed performance in organic dairy farming in Northern Spain.

Organic farm management combines best environmental practices, a high level of biodiversity, preservation of natural resources and high animal welfare standards. To meet these criteria, farmers must have livestock well adapted to local organic conditions and information about how different breeds and crosses perform under different conditions. The objective of this study was to compare the performance of different pure breeds and cross-breeds of cattle in organic dairy systems in Northern Spain. The data analysed were obtained from monthly records kept between 2010 and 2016 on organic farms registered in the regional milk recording system. Analysis of various traits indicated that the Holstein-Friesian breed suits the organic production system in the study region. Although the reproductive performance of Holstein-Friesian cows was poorer (in terms of number of services per conception) than that of cross-breed and Brown Swiss cows, the Holstein-Friesian produced more milk and lived longer. In addition, there was no difference in calving type or calving ease between the different breed groups. The better milk fat and protein yields produced by the crosses may be useful traits for farmers interested in milk transformation. The advantage of continuing to use Holstein-Friesian cattle is that the breed is predominant worldwide, and the genealogy is well documented. If Holstein-Friesian cattle continue to be used, the main priority will be to search for well-adapted bulls (particularly for pasture-based conditions) and to elaborate a genetic merit index for organic and pasture-based systems with the aim of predicting and minimizing genotype × environment interactions.

The effect of breed, sex, and drug concentration on the pharmacokinetic profile of ivermectin in cattle.

Ivermectin (IVM) is one of the most widely used antiparasitic drugs worldwide and has become the drug of choice for anthelmintic and tick treatment in beef cattle production. It is known that pharmacokinetic parameters are fundamental to the rational use of a drug and food safety and these parameters are influenced by different factors. The aim of this study was to evaluate the pharmacokinetic profile of IVM in Bos indicus, Bos taurus, and crossbreed cattle (B. indicus × B. taurus) kept under same field conditions and the possible impacts of sex and IVM formulation (1% and 3.15%). It was observed that IVM concentration was significantly affected by breed. The plasma concentrations of IVM, AUC, Cmax , and t1/2ß were significantly higher in B. indicus compared to B. taurus. Crossbreed animals showed an intermediate profile between European and Indian cattle. No alteration in pharmacokinetics parameters was detected when comparing different gender. Concerning the pharmacokinetic data of IVM formulation, it was verified that Tmax , AUC, and t1/2ß were higher in 3.15% IVM animals than those from 1% IVM formulation. The results clearly indicated that the IVM plasma concentrations in B. indicus were higher than that in B. taurus.

Contribution of connective tissue components, muscle fibres and marbling to beef tenderness variability in longissimus thoracis, rectus abdominis, semimembranosus and semitendinosus muscles.

BACKGROUND: The present study aimed to identify relationships between components of intramuscular connective tissue, proportions of the different fiber types, intramuscular fat and sensory tenderness of beef cooked at 55 °C. Accordingly, four muscles differing in their metabolic and contractile properties, as well as in their collagen content and butcher value, were obtained from dairy and beef cattle of several ages and sexes and were then used to create variability. RESULTS: Correlation analyses and/or stepwise regressions were applied on Z-scores to identify the existing and robust associations. Tenderness scores were further categorized into tender, medium and tough classes using unsupervised learning methods. The findings revealed a muscle-dependant role with respect to tenderness of total and insoluble collagen, cross-links, and type IIB + X and IIA muscle fibers. The longissimus thoracis and semitendinosus muscles that, in the present study, were found to be extreme in their tenderness potential were also very different from each other and from the rectus abdominis (RA) and semimembranosus (SM). RA and SM muscles were very similar regarding their relationship for muscle components and tenderness. A relationship between marbling and tenderness was only present when the results were analysed irrespective of all factors of variation of the experimental model relating to muscle and animal type. CONCLUSION: The statistical approaches applied in the present study using Z-scores allowed identification of the robust associations between muscle components and sensory beef tenderness and also identified discriminatory variables of beef tenderness classes. © 2020 Society of Chemical Industry.

Genetic diversity and relationships among three Southern African Nguni cattle populations.

The Nguni cattle breed has distinct populations that are adapted to the different ecological zones of Southern Africa. This study was carried out to assess genetic diversity and establish the relationships among South African (SA), Mozambican (Landim), and Swazi Nguni cattle populations, using 25 microsatellite markers. Genotypic data were generated from deoxyribonucleic acid (DNA) samples of 90 unrelated individuals of the three cattle populations, collected from government conservations and stud herds. DNA profiles of five local beef breeds were used as the reference populations. Most of the 25 microsatellite markers were highly polymorphic across the studied populations, with an overall polymorphic information content (PIC) mean of 0.676. Genetic diversity within populations was high with expected heterozygosity varying from 0.705 ± 0.024 (Landim) to 0.748 ± 0.021 (SA Nguni) and mean number of alleles being highest in the SA Nguni (7.52 ± 0.42). Average observed heterozygosity (0.597 ± 0.046) compared to the expected heterozygosity (0.719 ± 0.022) was lowest for the Swazi Nguni, which also had a high number of Hardy-Weinberg Equilibrium (HWE) deviated loci (13), confirming the relatively high level of inbreeding (0.158 ± 0.058) in that population. Analysis of molecular variance revealed only 9.61% of the total variation between the populations and 90.39% within populations. A short genetic distance (0.299) was observed between Landim and Swazi Nguni, with the SA Nguni (> 0.500) being the most genetically distant population. The distant relationship between SA Nguni and the other two Nguni cattle populations was further confirmed by a principal coordinates analysis. The three Nguni populations clustered independently from each other, despite some evidence of admixture. Therefore, it can be concluded that SA Nguni, Landim, and Swazi Nguni populations in Southern Africa exhibit high levels of genetic diversity and are genetically distant; with the two latter populations being less genetically apart. These results present useful information for the development of strategies for regional management of animal genetic resources, through conservation and sustainable utilisation.

About the inconvenience of handling mixed-breed herds; aspects of social behavior as a potential source of stress and economic losses.

The objective of this study was to assess aspects of the social behavior of a mixed-breed herd of beef cows as a potential source for stress and economic losses. Angus (AN; N = 10), Brahman (BR; N = 10), and Senepol (SE; N = 10) cows were assigned to two groups (N = 15 each containing equal breed numbers) on separate pastures. Agonistic interactions (win/loss) during feeding were recorded daily for 45 days. Dominance values were estimated as the proportion of individuals dominated to total herdmates. From this, individuals were placed into social categories based upon linear ranking as follows: dominants (D), intermediate (I), and subordinates (S). Breed influenced (P < 0.01) social category, with SE cows being dominants (P < 0.05) over AN and BR cows. Interactions between AN and BR cows were less (P < 0.0005) than interactions between AN and SE (53 vs 140, respectively). Within breeds, BR (152) and SE (182) cows had more (P < 0.0005) agonistic interactions than AN (107) cows. Although apparently influenced by breed, agonistic interactions occurred more frequently (P < 0.005) between social categories than within social categories (814 vs 310, respectively). Dominant cows were involved in more agonistic interactions with cows from different social categories than were intermediate and subordinate cows (P < 0.0005). However, intermediate (100) and subordinate (157) cows generated more (P < 0.0005) agonistic interactions within their own social category than dominant cows (53). It was concluded that, in mixed-breed herds, breed influences both social organization and agonistic interactions which could be considered as potential sources of stress and economic losses.

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

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

The effect of breed and diet type on the global transcriptome of hepatic tissue in beef cattle divergent for feed efficiency.

BACKGROUND: Feed efficiency is an important economic and environmental trait in beef production, which can be measured in terms of residual feed intake (RFI). Cattle selected for low-RFI (feed efficient) have similar production levels but decreased feed intake, while also emitting less methane. RFI is difficult and expensive to measure and is not widely adopted in beef production systems. However, development of DNA-based biomarkers for RFI may facilitate its adoption in genomic-assisted breeding programmes. Cattle have been shown to re-rank in terms of RFI across diets and age, while also RFI varies by breed. Therefore, we used RNA-Seq technology to investigate the hepatic transcriptome of RFI-divergent Charolais (CH) and Holstein-Friesian (HF) steers across three dietary phases to identify genes and biological pathways associated with RFI regardless of diet or breed. RESULTS: Residual feed intake was measured during a high-concentrate phase, a zero-grazed grass phase and a final high-concentrate phase. In total, 322 and 33 differentially expressed genes (DEGs) were identified across all diets for CH and HF steers, respectively. Three genes, GADD45G, HP and MID1IP1, were differentially expressed in CH when both the high-concentrate zero-grazed grass diet were offered. Two canonical pathways were enriched across all diets for CH steers. These canonical pathways were related to immune function. CONCLUSIONS: The absence of common differentially expressed genes across all dietary phases and breeds in this study supports previous reports of the re-ranking of animals in terms of RFI when offered differing diets over their lifetime. However, we have identified biological processes such as the immune response and lipid metabolism as potentially associated with RFI divergence emphasising the previously reported roles of these biological processes with respect to RFI.

A global analysis of CNVs in diverse yak populations using whole-genome resequencing.

BACKGROUND: Genomic structural variation represents a source for genetic and phenotypic variation, which may be subject to selection during the environmental adaptation and population differentiation. Here, we described a genome-wide analysis of copy number variations (CNVs) in 16 populations of yak based on genome resequencing data and CNV-based cluster analyses of these populations. RESULTS: In total, we identified 51,461 CNV events and defined 3174 copy number variation regions (CNVRs) that covered 163.8 Mb (6.2%) of yak genome with more "loss" events than both "gain" and "both" events, and we confirmed 31 CNVRs in 36 selected yaks using quantitative PCR. Of the total 163.8 Mb CNVR coverage, a 10.8 Mb region of high-confidence CNVRs directly overlapped with the 52.9 Mb of segmental duplications, and we confirmed their uneven distributions across chromosomes. Furthermore, functional annotation indicated that the CNVR-harbored genes have a considerable variety of molecular functions, including immune response, glucose metabolism, and sensory perception. Notably, some of the identified CNVR-harbored genes associated with adaptation to hypoxia (e.g., DCC, MRPS28, GSTCD, MOGAT2, DEXI, CIITA, and SMYD1). Additionally, cluster analysis, based on either individuals or populations, showed that the CNV clustering was divided into two origins, indicating that some yak CNVs are likely to arisen independently in different populations and contribute to population difference. CONCLUSIONS: Collectively, the results of the present study advanced our understanding of CNV as an important type of genomic structural variation in yak, and provide a useful genomic resource to facilitate further research on yak evolution and breeding.

Consequences of breed formation on patterns of genomic diversity and differentiation: the case of highly diverse peripheral Iberian cattle.

BACKGROUND: Iberian primitive breeds exhibit a remarkable phenotypic diversity over a very limited geographical space. While genomic data are accumulating for most commercial cattle, it is still lacking for these primitive breeds. Whole genome data is key to understand the consequences of historic breed formation and the putative role of earlier admixture events in the observed diversity patterns. RESULTS: We sequenced 48 genomes belonging to eight Iberian native breeds and found that the individual breeds are genetically very distinct with FST values ranging from 4 to 16% and have levels of nucleotide diversity similar or larger than those of their European counterparts, namely Jersey and Holstein. All eight breeds display significant gene flow or admixture from African taurine cattle and include mtDNA and Y-chromosome haplotypes from multiple origins. Furthermore, we detected a very low differentiation of chromosome X relative to autosomes within all analyzed taurine breeds, potentially reflecting male-biased gene flow. CONCLUSIONS: Our results show that an overall complex history of admixture resulted in unexpectedly high levels of genomic diversity for breeds with seemingly limited geographic ranges that are distantly located from the main domestication center for taurine cattle in the Near East. This is likely to result from a combination of trading traditions and breeding practices in Mediterranean countries. We also found that the levels of differentiation of autosomes vs sex chromosomes across all studied taurine and indicine breeds are likely to have been affected by widespread breeding practices associated with male-biased gene flow.

Genome-wide association study and scan for signatures of selection point to candidate genes for body temperature maintenance under the cold stress in Siberian cattle populations.

BACKGROUND: Design of new highly productive livestock breeds, well-adapted to local climatic conditions is one of the aims of modern agriculture and breeding. The genetics underlying economically important traits in cattle are widely studied, whereas our knowledge of the genetic mechanisms of adaptation to local environments is still scarce. To address this issue for cold climates we used an integrated approach for detecting genomic intervals related to body temperature maintenance under acute cold stress. Our approach combined genome-wide association studies (GWAS) and scans for signatures of selection applied to a cattle population (Hereford and Kazakh Whiteheaded beef breeds) bred in Siberia. We utilized the GGP HD150K DNA chip containing 139,376 single nucleotide polymorphism markers. RESULTS: We detected a single candidate region on cattle chromosome (BTA)15 overlapping between the GWAS results and the results of scans for selective sweeps. This region contains two genes, MSANTD4 and GRIA4. Both genes are functional candidates to contribute to the cold-stress resistance phenotype, due to their indirect involvement in the cold shock response (MSANTD4) and body thermoregulation (GRIA4). CONCLUSIONS: Our results point to a novel region on BTA15 which is a candidate region associated with the body temperature maintenance phenotype in Siberian cattle. The results of our research and the follow up studies might be used for the development of cattle breeds better adapted to cold climates of the Russian Federation and other Northern countries with similar climates.