Copy number variation of bovine DYNC1I2 gene is associated with body conformation traits in chinese beef cattle

Previous, studies have shown that the dynein transporter compound has a role in diseases such as intellectual disability and cerebral malformations. However, the study of CNV in DYNC1I2 gene has not been reported. Q-PCR and data association analysis were used for DYNC1I2 gene copy in this study. In this study, blood samples were collected from five breeds of Chinese cattle (Qingchuan cattle, Xianan cattle, Yunling cattle, Pinan cattle and Guyuan cattle) for DYNC1I2 gene CNV type detection. SPSS 20.0 software and method of ANOVA were used to analyzed the association between types of CNV and growth traits. Results reveal that the distribution of different copy number types in different cattle breeds is different. Association analysis indicate that CNV of DYNC1I2 gene showed a positive effect in cattle growth: in XN cattle, individuals with deletion types showed better performance on height at hip cross (P < 0.05); individuals with duplication types have better performance on body length (P < 0.05) in PN cattle; individuals with deletion types was significantly correlated with chest width and Hucklebone width (P < 0.05) in QC cattle; individuals with duplication types in Yunling cattle were better than the normal types, and there was a significant correlation between copy number variant and chest depth (P < 0.05). The results showed that CNV markers closely related to cattle production traits were detected at DNA level, which could be used as an important candidate molecular marker for marker-assisted selection of growth traits in Chinese cattle, and provided a new research basis for genetics and breeding of Chinese beef cattle.

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.

Grouping crossbred Holstein x Gyr heifers according to different feed efficiency indexes and its effects on energy and nitrogen partitioning, blood metabolic variables and gas exchanges

The objectives of this study were: i) to classify animals into groups of high and low feed efficiency (FE) using three FE indexes (Residual feed intake (RFI), Residual weight gain (RG) and Feed conversion efficiency (FCE)), and ii) to evaluate whether crossbreed Holstein x Gyr heifers divergent for FE indexes exhibit differences in nutrient intake and digestibility, energy partitioning, heat production, methane emissions, nitrogen partitioning and blood parameters. Thirty-five heifers were housed in a tie-stall, received ad libitum TMR (75:25, corn silage: concentrate) and were ranked and classified into high (HE) or low efficiency (LE) for RFI, RG and FCE. The number of animals for each HE group were 13 (< 0.5 standard deviation (SD) for RFI, 11 for RG and 11 for FCE (> 0.5 SD) and for the LE were 10 (> 0.5 SD) for RFI, 11 for RG and 12 for FCE (< 0.5 SD). Gas exchanges (O₂ consumption, CO₂ and CH₄ production) in open-circuit respiratory chambers and whole tract digestibility trial was performed. A completely randomized experimental design was used and the data were analyzed by ANOVA and correlation study. High efficiency animals for RFI produced less CO_2, consumed less O₂ and had lower heat production (HP). Methane production was positively correlated with RFI. High efficiency RG had higher O₂ consumption and CO₂ production in relation to LE-RG. High efficiency FCE had greater NFC digestibility, higher positive energy balance (EB) and excreted (11.4 g/d) less nitrogen in urine. High efficiency RG and FCE groups emitted less CH₄ per kg of weight gain than LE animals. Animals HE for RFI and FCE had lower β-hydroxybutyrate and higher glucose concentrations, respectively. The differences in intake, digestibility, energy and nitrogen partition, CH₄ emission, blood metabolic variables and heat production between the HE and LE groups varied according to the efficiency indexes adopted. The HP (kcal/d/BW^0.75) was lower for HE animals for RFI and FCE indexes.

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.

Uncovering Sub-Structure and Genomic Profiles in Across-Countries Subpopulations of Angus Cattle

Highlighting genomic profiles for geographically distinct subpopulations of the same breed may provide insights into adaptation mechanisms to different environments, reveal genomic regions divergently selected, and offer initial guidance to joint genomic analysis. Here, we characterized similarities and differences between the genomic patterns of Angus subpopulations, born and raised in Canada (N = 382) and Brazil (N = 566). Furthermore, we systematically scanned for selection signatures based on the detection of autozygosity islands common between the two subpopulations, and signals of divergent selection, via FST and varLD tests. The principal component analysis revealed a sub-structure with a close connection between the two subpopulations. The averages of genomic relationships, inbreeding coefficients, and linkage disequilibrium at varying genomic distances were rather similar across them, suggesting non-accentuated differences in overall genomic diversity. Autozygosity islands revealed selection signatures common to both subpopulations at chromosomes 13 (63.77-65.25 Mb) and 14 (22.81-23.57 Mb), which are notably known regions affecting growth traits. Nevertheless, further autozygosity islands along with FST and varLD tests unravel particular sites with accentuated population subdivision at BTAs 7 and 18 overlapping with known QTL and candidate genes of reproductive performance, thermoregulation, and resistance to infectious diseases. Our findings indicate overall genomic similarity between Angus subpopulations, with noticeable signals of divergent selection in genomic regions associated with the adaptation in different environments.

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.

Genomic divergence reveals unique populations among Indian Yaks

The present study focused upon identification of genome-wide SNPs through the reduced representation approach and to study the genomic divergence of the Indian yak populations. A total of 80 samples belonging to Arunachali yak (N = 20), Himachali yak (N = 20), Ladakhi yak (N = 20) and Sikkimi yak (N = 20) of India were used in the study. The results of the study revealed a total of 579575 high quality SNPs along with 50319 INDELs in the Indian yaks. The observed heterozygosity was found to be high in Himachali yak, followed by Arunachali yak, Ladakhi yak and Sikkimi yaks. The Sikkimi yaks was found to be genetically distant, followed by Ladakhi yaks which was observed to have some few individuals from Arunachali and Himachali yaks. Arunachali and Himachali yaks are found to get clustered together and are genetically similar. The study provides evidence about the genomic diversity in the Indian yak populations and information generated in the present study may help to formulate a suitable breeding plan for endangered Indian yaks. Moreover, the unique yak populations identified in the study will further help to focus attention for future characterization and prioritization of the animals for conservation purposes through the ddRAD approach.

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.

Grazing behaviour of dairy cows on biodiverse mountain pastures is more influenced by slope than cow breed

The aim of this study was to determine how cows with different genetic merit behave and perform when grazing biodiverse and heterogeneous mountain pastures with different slopes. Three groups of 12 cows in late lactation, each composed of four Holstein, four Montbéliarde and four Valdostana Red Pied cows, breeds of increasing presumed robustness and decreasing milk yield (MY) potential. Cows grazed without concentrate either on a low-diversity flat pasture or on two species-rich mountainous pastures having slopes of either 7° or 22°. Milk yield, BW and grazing behaviour were monitored two times in the first and once in the second grazing cycle. Cows of different breeds had similar behaviour on all pastures. The Montbéliarde cows performed close to their production potential; Holstein and Valdostana cows produced less milk than anticipated. No breed difference in terms of BW loss was found. The Valdostana cows exhibited the least selective behaviour with respect to plant species and plant growth stage. Still, all cows searched for the most palatable vegetation regardless of pasture diversity. On the steep pasture, cows optimised the trade-off between ingesting and saving energy to obtain feed. They remained longer at the lowest zone and selected forbs, whereas cows on the flatter pasture went to the upper zone to select grasses. The present study gave no evidence for a superior short-term adaptation to harsh grazing conditions through an optimised feeding behaviour of the Valdostana breed compared to Montbéliarde and Holstein cows.

Health, behaviour and growth performance of Charolais and Limousin bulls fattened on different types of flooring

Intensive fattening of late-maturing breeds on concrete or rubberized slatted floors is the prevalent beef production system in mainland Europe. The rationale behind this study is that specific beef breeds with different slaughter weights might have a diverse response to different flooring systems. The study aimed at assessing whether growth performance, health, behaviour and claw condition of two beef breeds, Charolais (CH) and Limousine (LIM), would be affected by their housing on concrete (CS) or rubber-covered (RCS) fully slatted floor. A total of 228 CH (116 on CS; 112 on RCS) and 115 LIM (57 on CS; 58 on RCS) were housed in four and two commercial farms, respectively, in groups of 9.0 ± 2.1 animals/pen with an average space allowance of 3.1 ± 0.2 m2. Draining gaps of CS and RCS pens were 16.9 ± 1.7% and 11.6 ± 1.2% of the total surface, respectively. Bulls of both breeds had similar initial body weight (429.4 ± 31.5 kg for CH; 369.6 ± 31.7 kg for LIM), and they were slaughtered when they reached suitable finishing. Charolais had a higher final body weight (BW) than LIM (750.8 ± 8.6 v. 613.7 ± 10.9 kg; P < 0.01), and bulls of both breeds on RCS had higher average daily gain than on CS (1.47 ± 0.02 v. 1.39 ± 0.02 kg/day; P < 0.05). The percentage of bulls early culled or treated for locomotor disorders were reduced by RCS only for LIM, while RCS tended to prevent the occurrence of bursitis for both breeds. During two 8-h behavioural observations, bulls on RCS performed more head butt/displacements and chases than on CS, and they reduced the frequency of abnormal lying down events. The use of RCS increased mounts' frequency only in LIM, while its reduced drainage capacity impaired only the cleanliness of CH. Postmortem hoof inspection showed longer claw dorsal wall and diagonal lengths, and sharper toe angles for CH on RCS than LIM on both floors. Results of this study point out that fully slatted floors, regardless of being rubberized or not, are not suitable for bulls finished at a final BW above 700 kg due to their detrimental effects on health and welfare. The use of RCS could be recommended as an alternative to CS only if bulls are slaughtered at a lower final BW (around 600 kg), like in the case of LIM breed.

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.

Brangus cows have ovarian reserve parameters more like Brahman than Angus cows

Bos indicus females have more surface antral follicles than Bos taurus females; however, histological studies demonstrated no difference in total number of primordial follicles between these two biological types of cattle. Primordial follicle density in the ovary was less in Nelore ovaries compared to Angus ovaries, but no studies have examined the primordial follicle density in Bos indicus cross-bred females. It, therefore, was hypothesized that primordial follicle density in the ovary would decrease as percentage Bos indicus increased. Ovaries were collected from cross-bred Angus (n = 32, no Bos indicus influence), Brangus (n = 15), or Brahman (n = 9) cows and prepared for histological evaluation. There was no difference in total number of primordial follicles per ovary between breeds (P > 0.10). When numbers of primordial follicles were expressed on a per gram of ovarian tissue basis, there were fewer primordial follicles per gram of ovarian tissue in Brangus and Brahman cows than in Angus cows (P < 0.05). Brangus cows did not differ from Brahman cows in primordial follicle density (P > 0.10). Differences in primordial follicle density could indicate differences in capacity of ovarian stroma to produce factors necessary for oogonial proliferation and primordial follicle formation among breeds. Identifying these factors could improve the aprroach for culturing pre-antral follicles of cattle. Furthermore, these results explain why ultrasonographic antral follicle counts may need to be adjusted to a greater threshold to predict size of the ovarian reserve and determine ovarian reserve related reproductive traits in Bos indicus females.

Combined effects of type and depth of counting chamber, and rate of image frame capture, on bull sperm motility and kinematics

Semen quality assessment requires accurate, reliable and objective methods for examination of sperm variables including sperm motility. For preparation of semen samples for artificial insemination, as a genetic resource, samples that are used for insemination need to have the capacity to result in a highly acceptable fertility rate. Several methods have been developed for evaluation of bull sperm in laboratory conditions and for preparation of doses for artificial insemination. Computer-assisted semen analyses can provide objective information on various sperm variables. Nevertheless, this equipment requires fine-calibrations considering differences among species, breeds and conditions for sample evaluation and data analyses. In the present study, there was examination of the interaction between factors such as image frame rate and type and depth of counting chamber in which sperm were evaluated, together with differences between bulls of four breeds. The use of the Spermtrack® reusable 10 μm-depth chamber provided more reliable results than results obtained using disposable chambers (10 and 20 μm depth). A capture rate of at least 90 fps is required for assessment of sperm motility percentage, whereas a rate of 250 fps is needed for obtaining consistent kinematic data. Differences among breeds in the present study indicate conditions for sperm analyses should include specific equipment calibrations for each breed. These results contribute to development of more precise conditions for assessments of bull sperm quality taking into account breed differences and the requirement each breed has for the adequate evaluation and preparation of samples for artificial insemination.

Economic performance evaluation of Brangus and Nellore cattle breed naturally infested with Rhipicephalus microplus in an extensive production system in Central-West Brazil

The objective of the present study was to evaluate the performance of Brangus and Nellore cattle in an extensive production system. Sixty growing bulls, 30 of the Brangus and 30 of the Nellore breeds that were naturally infested with Rhipicephalus microplus were subdivided into four groups: control Nellore (15 animals), treated Nellore (15), control Brangus (15), and treated Brangus (15). The animals in the treated groups underwent acaricidal treatments for 1 year. Tick counts, acaricidal treatments and animal weighing were performed every 18 days, and the costs of acaricidal, anti-myiasis, and preventive treatments for tick fever were recorded for cost evaluation. The treated Brangus and Nellore groups did not show a significant difference in weight gain, whereas Nellore weight gain was superior in the control groups. The cost of acaricidal treatment throughout the experimental period was $494 US. The costs of the preventive treatment for tick fever and myiasis were $98 US and $15 US, respectively. The highest rates of tick infestation were found in the control group of the Brangus animals, which served as the basis for the suggested implementation of a strategic control program for animals in the growth phase. Nellore animals showed low rates of infestation. Under the conditions of this study, the Nellore animals were more efficient than the Brangus animals because they achieved satisfactory weight gain, similar to the Brangus, and a low tick count even in the control group. This finding demonstrates that expenses related to the acaricidal, anti-myiasis and tick fever treatments are unnecessary in Nellore cattle, making these animals more profitable.

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.

Genomic variants identified from whole-genome resequencing of indicine cattle breeds from Pakistan

The primary goal of cattle genomics is the identification of genome-wide polymorphism associated with economically important traits. The bovine genome sequencing project was completed in 2009. Since then, using massively parallel sequencing technologies, a large number of Bos taurus cattle breeds have been resequenced and scanned for genome-wide polymorphisms. As a result, a substantial number of single nucleotide polymorphisms (SNPs) have been discovered across European Bos taurus genomes, whereas extremely less number of SNPs are cataloged for Bos indicus breeds. In this study, we performed whole-genome resequencing, reference-based mapping, functional annotation and gene enrichment analysis of 20 sires representing eleven important Bos indicus (indicine) breeds of Pakistan. The breeds sequenced here include: Sahiwal, Red Sindhi, Tharparkar and Cholistani (tropically adapted dairy and dual purpose breeds), Achai, Bhagnari, Dajal and Lohani (high altitude adapted dual and drought purpose breeds); Dhanni, Hisar Haryana and Gabrali (dairy and light drought purpose breeds). A total of 17.4 billion QC passed reads were produced using BGISEQ-500 next generation sequencing platform to generate 9 to 27-fold genome coverage (average ~16×) for each of the 20 sequenced sires. A total of 67,303,469 SNPs were identified, of which 3,850,365 were found novel and 1,083,842 insertions-deletions (InDels) were detected across the whole sequenced genomes (491,247 novel). Comparative analysis using coding region SNPs revealed a close relationship between the best milking indicine breeds; Red Sindhi and Sahiwal. On the other hand, Bhagnari and Tharparkar being popular for their adaptation to dry and extremely hot climates were found to share the highest number of SNPs. Functional annotation identified a total of 3,194 high-impact (disruptive) SNPs and 745 disruptive InDels (in 275 genes) that may possibly affect economically important dairy and beef traits. Functional enrichment analysis was performed and revealed that high or moderate impact variants in wingless-related integration site (Wnt) and vascular smooth muscle contraction (VSMC) signaling pathways were significantly over-represented in tropically adapted heat tolerant Pakistani-indicine breeds. On the other hand, vascular endothelial growth factor (VEGF) and hypoxia-inducible factor 1 (HIF-1) signaling pathways were found over-represented in highland adapted Pakistani-indicine breeds. Similarly, the ECM-receptor interaction and Jak-STAT signaling pathway were significantly enriched in dairy and beef purpose Pakistani-indicine cattle breeds. The Toll-like receptor signaling pathway was significantly enriched in most of the Pakistani-indicine cattle. Therefore, this study provides baseline data for further research to investigate the molecular mechanisms of major traits and to develop potential genomic markers associated with economically important breeding traits, particularly in indicine cattle.

A Bayesian approach for analysis of ordered categorical responses subject to misclassification

Ordinal categorical responses are frequently collected in survey studies, human medicine, and animal and plant improvement programs, just to mention a few. Errors in this type of data are neither rare nor easy to detect. These errors tend to bias the inference, reduce the statistical power and ultimately the efficiency of the decision-making process. Contrarily to the binary situation where misclassification occurs between two response classes, noise in ordinal categorical data is more complex due to the increased number of categories, diversity and asymmetry of errors. Although several approaches have been presented for dealing with misclassification in binary data, only limited practical methods have been proposed to analyze noisy categorical responses. A latent variable model implemented within a Bayesian framework was proposed to analyze ordinal categorical data subject to misclassification using simulated and real datasets. The simulated scenario consisted of a discrete response with three categories and a symmetric error rate of 5% between any two classes. The real data consisted of calving ease records of beef cows. Using real and simulated data, ignoring misclassification resulted in substantial bias in the estimation of genetic parameters and reduction of the accuracy of predicted breeding values. Using our proposed approach, a significant reduction in bias and increase in accuracy ranging from 11% to 17% was observed. Furthermore, most of the misclassified observations (in the simulated data) were identified with a substantially higher probability. Similar results were observed for a scenario with asymmetric misclassification. While the extension to traits with more categories between adjacent classes is straightforward, it could be computationally costly. For traits with high heritability, the performance of the methodology would be expected to improve.

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.

A copy number variant scan in the autochthonous Valdostana Red Pied cattle breed and comparison with specialized dairy populations

Copy number variants (CNVs) are an important source of genomic structural variation, recognized to influence phenotypic variation in many species. Many studies have focused on identifying CNVs within and between human and livestock populations alike, but only few have explored population-genetic properties in cattle based on CNVs derived from a high-density SNP array. We report a high-resolution CNV scan using Illumina’s 777k BovineHD Beadchip for Valdostana Red Pied (VRP), an autochthonous Italian dual-purpose cattle population reared in the Alps that did not undergo strong selection for production traits. After stringent quality control and filtering, CNVs were called across 108 bulls using the PennCNV software. A total of 6,784 CNVs were identified, summarized to 1,723 CNV regions (CNVRs) on 29 autosomes covering a total of ~59 Mb of the UMD3.1 assembly. Among the mapped CNVRs, there were 812 losses, 832 gains and 79 complexes. We subsequently performed a comparison of CNVs detected in the VRP and those available from published studies in the Italian Brown Swiss (IBS) and Mexican Holstein (HOL). A total of 171 CNVRs were common to all three breeds. Between VRP and IBS, 474 regions overlapped, while only 313 overlapped between VRP and HOL, indicating a more similar genetic background among populations with common origins, i.e. the Alps. The principal component, clustering and admixture analyses showed a clear separation of the three breeds into three distinct clusters. In order to describe the distribution of CNVs within and among breeds we used the pair V_ST statistic, considering only the CNVRs shared to more than 5 individuals (within breed). We identified unique and highly differentiated CNVs (n = 33), some of which could be due to specific breed selection and adaptation. Genes and QTL within these regions were characterized.

Metabolomics Profiling of Serum and Urine in Three Beef Cattle Breeds Revealed Different Levels of Tolerance to Heat Stress

This study was to determine differences in the global metabolic profiles of serum and urine of Xuanhan yellow cattle, Simmental crossbred cattle (Simmental × Xuanhan yellow cattle), and cattle-yaks (Jersey × Maiwa yak) under heat stress (temperature-humidity index remained above 80 for 1 week). A total of 55 different metabolites associated with the three breeds were identified in the serum and urine samples by gas chromatography-mass spectrometry. The metabolic adaptations to heat stress are heterogeneous. Cattle-yaks mobilize a greater amount of body protein to release glucogenic amino acids to supply energy, whereas the tricarboxylic acid cycle is inhibited. Simmental crossbred cattle mobilize a greater amount of body fat to use free fatty acids as an energy source. In comparison with Simmental crossbred cattle and cattle-yaks, Xuanhan yellow cattle have higher glycolytic activity and possess a stronger antioxidant defense system and are, in conclusion, more adapted to hot and humid environments.

A cerebellar substrate for cognition evolved multiple times independently in mammals

Given that complex behavior evolved multiple times independently in different lineages, a crucial question is whether these independent evolutionary events coincided with modifications to common neural systems. To test this question in mammals, we investigate the lateral cerebellum, a neurobiological system that is novel to mammals, and is associated with higher cognitive functions. We map the evolutionary diversification of the mammalian cerebellum and find that relative volumetric changes of the lateral cerebellar hemispheres (independent of cerebellar size) are correlated with measures of domain-general cognition in primates, and are characterized by a combination of parallel and convergent shifts towards similar levels of expansion in distantly related mammalian lineages. Results suggest that multiple independent evolutionary occurrences of increased behavioral complexity in mammals may at least partly be explained by selection on a common neural system, the cerebellum, which may have been subject to multiple independent neurodevelopmental remodeling events during mammalian evolution.

Genotype by environment interaction effects in genetic evaluation of preweaning gain for Line 1 Hereford cattle from Miles City, Montana

It has long been recognized that genotype × environment interaction potentially influences genetic evaluation of beef cattle. However, this recognition has largely been ignored in systems for national cattle evaluation. The objective of this investigation was to determine if direct and maternal genetic effects on preweaning gain would be reranked depending on an environmental gradient as determined by year effects. Data used were from the 76-yr selection experiment with the Line 1 Hereford cattle raised at Miles City, MT. The data comprised recorded phenotypes from 7,566 animals and an additional 1,862 ancestral records included in the pedigree. The presence of genotype × environment interaction was examined using reaction norms wherein year effects on preweaning gain were hypothesized to linearly influence the EBV. Estimates of heritability for direct and maternal effects, given the average environment, were 10 ± 2 and 26 ± 3%, respectively. In an environment that is characterized by the 5th (95th) percentile of the distribution of year effects, the corresponding estimates of heritability were 18 ± 3 (22 ± 3%) and 30 ± 3% (30 ± 3%), respectively. Rank correlations of direct and maternal EBV appropriate to the 5th and 95th percentiles of the year effects were 0.67 and 0.92, respectively. In the average environment, the genetic trends were 255 ± 1 g/yr for direct effects and 557 ± 3 g/yr for maternal effects. In the fifth percentile environment, the corresponding estimates of genetic trend were 271 ± 1 and 540 ± 3 g/yr, respectively, and in the 95th percentile environment, they were 236 ± 1 and 578 ± 3 g/yr, respectively. Linear genetic trends in environmental sensitivity were observed for both the direct (-8.06 × 10 ± 0.49 × 10) and maternal (8.72 × 10 ± 0.43 × 10) effects. Therefore, changing systems of national cattle evaluation to more fully account for potential genotype × environment interaction would improve the assessment of breeding stock, particularly for direct effects. Estimates of environmental sensitivity parameters could also facilitate identification of genetic limitations to production.

Genetic parameters of carcass and meat quality traits in different muscles (longissimus dorsi and semimembranosus) of Hanwoo (Korean cattle)

We estimated heritability () and genetic and phenotypic correlations for carcass and meat quality traits of longissimus dorsi (LD) and semimembranosus (SM) muscles in 30-mo-old Hanwoo steers. Variance and covariance components were estimated using REML procedures under univariate and bivariate models. The mean carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT), and marbling score (MS) were 428.20 ± 46.30 kg, 87.38 ± 8.54 cm2, 13.00 ± 5.14 mm, and 5.21 ± 1.56, respectively. The mean CIE reflectance of meat lightness (L*), redness (a*), and yellowness (b*) were 40.01 ± 2.73, 22.37 ± 2.18, and 10.35 ± 1.46, respectively, in LD muscles and 36.33 ± 2.44, 22.91 ± 2.43, and 10.25 ± 1.65, respectively, in SM muscles. The mean Warner-Bratzler shear force (WBSF), intramuscular fat content (IMF), water-holding capacity (WHC), and protein and ash content in LD and SM muscles were 3.84 ± 0.96 and 6.52 ± 1.21 kg, 15.91 ± 4.39 and 5.10 ± 1.94%, 62.07 ± 3.38 and 71.61 ± 2.06%, 20.01 ± 1.39 and 21.34 ± 0.89%, and 0.80 ± 0.10 and 0.93 ± 0.07, respectively. The estimates of CWT, EMA, BFT, and MS were 0.51 ± 0.13, 0.45 ± 0.13, 0.29 ± 0.09, and 0.22 ± 0.08, respectively. The estimates were moderate for meat quality traits and were 0.37 ± 0.12, 0.40 ± 0.12, 0.33 ± 0.10, 0.33 ± 0.10, 0.30 ± 0.11, and 0.24 ± 0.09 for L*, WBSF, IMF, WHC, and protein and ash content, respectively, in LD muscle; estimates from SM muscle were comparatively low (0.08 ± 0.06 to 0.25 ± 0.09). Estimates of for a* and b* were also low (0.08 ± 0.06 to 0.13 ± 0.07). Carcass weight had a moderate, positive genetic correlation with EMA (0.60 ± 0.13) and a weak correlation with MS and BFT. The genetic correlations among the 3 colorimeter variants were strong and positive within and between muscles. Intramuscular fat content had moderate to strong and negative genetic correlations with WBSF (-0.49 ± 0.18), WHC (-0.99 ± 0.01), and protein (-0.93 ± 0.04) and ash content (-0.98 ± 0.06) in LD muscle, whereas the associations were less pronounced in SM muscle. In general, CWT and EMA had low genetic and phenotypic correlations with meat quality traits, which suggests that the traits are independent and have distinct genetic contributions in each muscle. Conversely, with few exceptions, meat quality traits had genetic and phenotypic correlations with MS and BFT. In conclusion, the estimated genetic parameters for carcass and meat quality traits could be used for genetic evaluation and breeding programs in Korean Hanwoo cattle populations.

Signatures of positive selection in African Butana and Kenana dairy zebu cattle

Butana and Kenana are two types of zebu cattle found in Sudan. They are unique amongst African indigenous zebu cattle because of their high milk production. Aiming to understand their genome structure, we genotyped 25 individuals from each breed using the Illumina BovineHD Genotyping BeadChip. Genetic structure analysis shows that both breeds have an admixed genome composed of an even proportion of indicine (0.75 ± 0.03 in Butana, 0.76 ± 0.006 in Kenana) and taurine (0.23 ± 0.009 in Butana, 0.24 ± 0.006 in Kenana) ancestries. We also observe a proportion of 0.02 to 0.12 of European taurine ancestry in ten individuals of Butana that were sampled from cattle herds in Tamboul area suggesting local crossbreeding with exotic breeds. Signatures of selection analyses (iHS and Rsb) reveal 87 and 61 candidate positive selection regions in Butana and Kenana, respectively. These regions span genes and quantitative trait loci (QTL) associated with biological pathways that are important for adaptation to marginal environments (e.g., immunity, reproduction and heat tolerance). Trypanotolerance QTL are intersecting candidate regions in Kenana cattle indicating selection pressure acting on them, which might be associated with an unexplored level of trypanotolerance in this cattle breed. Several dairy traits QTL are overlapping the identified candidate regions in these two zebu cattle breeds. Our findings underline the potential to improve dairy production in the semi-arid pastoral areas of Africa through breeding improvement strategy of indigenous local breeds.

Genomic analyses of tropical beef cattle fertility based on genotyping pools of Brahman cows with unknown pedigree

We introduce an innovative approach to lowering the overall cost of obtaining genomic EBV (GEBV) and encourage their use in commercial extensive herds of Brahman beef cattle. In our approach, the DNA genotyping of cow herds from 2 independent properties was performed using a high-density bovine SNP chip on DNA from pooled blood samples, grouped according to the result of a pregnancy test following their first and second joining opportunities. For the DNA pooling strategy, 15 to 28 blood samples from the same phenotype and contemporary group were allocated to pools. Across the 2 properties, a total of 183 pools were created representing 4,164 cows. In addition, blood samples from 309 bulls from the same properties were also taken. After genotyping and quality control, 74,584 remaining SNP were used for analyses. Pools and individual DNA samples were related by means of a "hybrid" genomic relationship matrix. The pooled genotyping analysis of 2 large and independent commercial populations of tropical beef cattle was able to recover significant and plausible associations between SNP and pregnancy test outcome. We discuss 24 SNP with significant association ( < 1.0 × 10) and mapped within 40 kb of an annotated gene. We have established a method to estimate the GEBV in young herd bulls for a trait that is currently unable to be predicted at all. In summary, our novel approach allowed us to conduct genomic analyses of fertility in 2 large commercial Brahman herds managed under extensive pastoral conditions.

Single nucleotide variants and InDels identified from whole-genome re-sequencing of Guzerat, Gyr, Girolando and Holstein cattle breeds

Whole-genome re-sequencing, alignment and annotation analyses were undertaken for 12 sires representing four important cattle breeds in Brazil: Guzerat (multi-purpose), Gyr, Girolando and Holstein (dairy production). A total of approximately 4.3 billion reads from an Illumina HiSeq 2000 sequencer generated for each animal 10.7 to 16.4-fold genome coverage. A total of 27,441,279 single nucleotide variations (SNVs) and 3,828,041 insertions/deletions (InDels) were detected in the samples, of which 2,557,670 SNVs and 883,219 InDels were novel. The submission of these genetic variants to the dbSNP database significantly increased the number of known variants, particularly for the indicine genome. The concordance rate between genotypes obtained using the Bovine HD BeadChip array and the same variants identified by sequencing was about 99.05%. The annotation of variants identified numerous non-synonymous SNVs and frameshift InDels which could affect phenotypic variation. Functional enrichment analysis was performed and revealed that variants in the olfactory transduction pathway was over represented in all four cattle breeds, while the ECM-receptor interaction pathway was over represented in Girolando and Guzerat breeds, the ABC transporters pathway was over represented only in Holstein breed, and the metabolic pathways was over represented only in Gyr breed. The genetic variants discovered here provide a rich resource to help identify potential genomic markers and their associated molecular mechanisms that impact economically important traits for Gyr, Girolando, Guzerat and Holstein breeding programs.

Bovine Milk Oligosaccharide Contents Show Remarkable Seasonal Variation and Intercow Variation

Human milk oligosaccharides (OS) play an important role in protecting the neonate. In addition to fructo-oligosaccharides and galacto-oligosaccharides, bovine milk OS have great potential to be used in pediatric food products to mimic the functions of human milk OS. Currently, little is known about the accumulation of OS in bovine milk in relation to genetic and environmental factors. A systematic survey on seasonal variation of 14 major OS was thus conducted with 19 cows over the entire milking season using a liquid chromatography-mass spectrometry technique. This study revealed a number of significant correlations between structurally related and structurally nonrelated OS and a substantial individual animal difference for all 14 OS. Most of the 14 OS displayed a remarkable seasonal variation in abundance (up to 10-fold change), with the highest abundance observed in April and May (i.e., autumn) for the majority of the 19 cows.

Development of a Bioinformatics Framework for the Detection of Gene Conversion and the Analysis of Combinatorial Diversity in Immunoglobulin Heavy Chains in Four Cattle Breeds

We have developed a new bioinformatics framework for the analysis of rearranged bovine heavy chain immunoglobulin (Ig) variable regions by combining and refining widely used alignment algorithms. This bioinformatics framework allowed us to investigate alignments of heavy chain framework regions (FRHs) and the separate alignments of FRHs and heavy chain complementarity determining regions (CDRHs) to determine their germline origin in the four cattle breeds Aubrac, German Black Pied, German Simmental, and Holstein Friesian. Now it is also possible to specifically analyze Ig heavy chains possessing exceptionally long CDR3Hs. In order to gain more insight into breed specific differences in Ig combinatorial diversity, somatic hypermutations and putative gene conversions of IgG, we compared the dominantly transcribed variable (IGHV), diversity (IGHD), and joining (IGHJ) segments and their recombination in the four cattle breeds. The analysis revealed the use of 15 different IGHV segments, 21 IGHD segments, and two IGHJ segments with significant different transcription levels within the breeds. Furthermore, there are preferred rearrangements within the three groups of CDR3H lengths. In the sequences of group 2 (CDR3H lengths (L) of 11–47 amino acid residues (aa)) a higher number of recombination was observed than in sequences of group 1 (L≤10 aa) and 3 (L≥48 aa). The combinatorial diversity of germline IGHV, IGHD, and IGHJ-segments revealed 162 rearrangements that were significantly different. The few preferably rearranged gene segments within group 3 CDR3H regions may indicate specialized antibodies because this length is unique in cattle. The most important finding of this study, which was enabled by using the bioinformatics framework, is the discovery of strong evidence for gene conversion as a rare event using pseudogenes fulfilling all definitions for this particular diversification mechanism.

[Investigation of metabolic parame- ters in high yielding dairy cows in pasture based production systems]

In the present study differences in metabolism between New Zealand Holstein-Friesian (NZHF) and Brown Swiss (CH-BV) or Swiss Holstein-Friesian (CH-HF) were investigated in a grassland based milk production system in Switzerland. Therefore 14 pairs of CH-BV/NZHF and 11 pairs of CH-HF/NZHF were available. The parameters glucose, insulin, non-esterified fatty acids (NEFA), β-hydroxybutyrate (β-HB), urea and cholesterol were analysed at the times 5-3 weeks before the calculated partus and 2, 3, 5, 7, 10 and 18-22 weeks post partum. Only β-HB showed significantly higher concentrations (P = 0.0059) for both Swiss breeds compared to the NZ-HF. Regarding all other physiological parameters during early lactation New Zealand Holstein-Friesians were not different from Swiss breeds.

Genome-Wide Survey and Analysis of Microsatellite Sequences in Bovid Species

Microsatellites or simple sequence repeats (SSRs) have become the most popular source of genetic markers, which are ubiquitously distributed in many eukaryotic and prokaryotic genomes. This is the first study examining and comparing SSRs in completely sequenced genomes of the Bovidae. We analyzed and compared the number of SSRs, relative abundance, relative density, guanine-cytosine (GC) content and proportion of SSRs in six taxonomically different bovid species: Bos taurus, Bubalus bubalis, Bos mutus, Ovis aries, Capra hircus, and Pantholops hodgsonii. Our analysis revealed that, based on our search criteria, the total number of perfect SSRs found ranged from 663,079 to 806,907 and covered from 0.44% to 0.48% of the bovid genomes. Relative abundance and density of SSRs in these Bovinae genomes were non-significantly correlated with genome size (Pearson, r < 0.420, p > 0.05). Perfect mononucleotide SSRs were the most abundant, followed by the pattern: perfect di- > tri- > penta- > tetra- > hexanucleotide SSRs. Generally, the number of SSRs, relative abundance, and relative density of SSRs decreased as the motif repeat length increased in each species of Bovidae. The most GC-content was in trinucleotide SSRs and the least was in the mononucleotide SSRs in the six bovid genomes. The GC-contents of tri- and pentanucleotide SSRs showed a great deal of similarity among different chromosomes of B. taurus, O. aries, and C. hircus. SSR number of all chromosomes in the B. taurus, O.aries, and C. hircus is closely positively correlated with chromosome sequence size (Pearson, r > 0.980, p < 0.01) and significantly negatively correlated with GC-content (Pearson, r < -0.638, p < 0.01). Relative abundance and density of SSRs in all chromosomes of the three species were significantly negatively correlated with GC-content (Pearson, r < -0.333, P < 0.05) but not significantly correlated with chromosome sequence size (Pearson, r < -0.185, P > 0.05). Relative abundances of the same nucleotide SSR type showed great similarity among different chromosomes of B. taurus, O. aries, and C. hircus.

Whole-Genome Resequencing Analysis of Hanwoo and Yanbian Cattle to Identify Genome-Wide SNPs and Signatures of Selection

Over the last 30 years, Hanwoo has been selectively bred to improve economically important traits. Hanwoo is currently the representative Korean native beef cattle breed, and it is believed that it shared an ancestor with a Chinese breed, Yanbian cattle, until the last century. However, these two breeds have experienced different selection pressures during recent decades. Here, we whole-genome sequenced 10 animals each of Hanwoo and Yanbian cattle (20 total) using the Illumina HiSeq 2000 sequencer. A total of approximately 3.12 and 3.07 billion sequence reads were mapped to the bovine reference sequence assembly (UMD 3.1) at an average of approximately 10.71- and 10.53-fold coverage for Hanwoo and Yanbian cattle, respectively. A total of 17,936,399 single nucleotide polymorphisms (SNPs) were yielded, of which 22.3% were found to be novel. By annotating the SNPs, we further retrieved numerous nonsynonymous SNPs that may be associated with traits of interest in cattle. Furthermore, we performed whole-genome screening to detect signatures of selection throughout the genome. We located several promising selective sweeps that are potentially responsible for economically important traits in cattle; the PPP1R12A gene is an example of a gene that potentially affects intramuscular fat content. These discoveries provide valuable genomic information regarding potential genomic markers that could predict traits of interest for breeding programs of these cattle breeds.

Exploratory analysis of differences in sperm morphology in Nelore and Gir (Bos indicus) bulls

The aim of this study was to evaluate the effects of breed and season on semen quality parameters of zebu bulls. Data (1,632 registers) of semen production from Gir (n = 4) and Nelore (n = 15) bulls were collected between October 2005 and November 2009. The ejaculates were collected twice a week during various seasons (summer, fall, winter, and spring) and evaluated for the following semen parameters: ejaculate volume, sperm concentration, sperm motility, forward progressive motility (FPM), and sperm morphology. Factor analysis was used to determine the relationship among variables. The effect of breed (Gir and Nelore) and season and their cross effect on each parameter and extracted factor were tested using ANOVA. A negative correlation (P < 0.05) was observed between FPM and proximal droplet, as well as with abnormal loose head, abnormal small head, pouch formation, abnormal mid-piece, and strongly folded tail. Gir bull sperm showed more major defects, detached acrosome, and minor FPM (P < 0.01), whereas Nelore bulls showed a higher number of sperm with normally loose head.

Copy number variations of MICAL-L2 shaping gene expression contribute to different phenotypes of cattle

Copy number variations (CNVs) recently have been recognized as an important source of genetic variability. Compelling evidence has indicated that CNVs are responsible for phenotypic traits by altering the copy numbers of functional genes. The molecule interacting with CasL-like protein 2 (MICAL-L2) gene plays a critical role in muscle fiber development and has been identified in the CNV region by comparative genomic hybridization array. In the present study, we detected the different distributions of MICAL-L2 gene copy numbers in four Chinese cattle breeds (QC, NY, LX, and CY) and investigated the functional effects of MICAL-L2 CNVs on the gene's expression level and the phenotypic traits in QC and NY cattle. The results showed that the copy number loss (relative to Angus cattle) was more frequent in CY than in the other breeds. The MICAL-L2 gene copy number presented a moderate negative correlation with the transcriptional expression in fetal skeletal muscles (P < 0.05). Statistical analysis revealed that the MICAL-L2 CNVs were significantly associated with body weight, body height, and body length of NY cattle in the early stages (6 and 12 months old), and the copy number loss showed better traits than the gain and/or median groups (P < 0.05). No significance was found at the late stages in QC (24 months old) and NY cattle (18 and 24 months old). These observations provided further insight into the associations between cattle CNVs and economic traits, suggesting that the CNVs may be considered promising markers for the molecular breeding of Chinese beef cattle.

Ancient mtDNA analysis of early 16(th) century Caribbean cattle provides insight into founding populations of New World creole cattle breeds

The Columbian Exchange resulted in a widespread movement of humans, plants and animals between the Old and New Worlds. The late 15^th to early 16^th century transfer of cattle from the Iberian Peninsula and Canary Islands to the Caribbean laid the foundation for the development of American creole cattle (Bos taurus) breeds. Genetic analyses of modern cattle from the Americas reveal a mixed ancestry of European, African and Indian origins. Recent debate in the genetic literature centers on the ‘African’ haplogroup T1 and its subhaplogroups, alternatively tying their origins to the initial Spanish herds, and/or from subsequent movements of taurine cattle through the African slave trade. We examine this problem through ancient DNA analysis of early 16^th century cattle bone from Sevilla la Nueva, the first Spanish colony in Jamaica. In spite of poor DNA preservation, both T3 and T1 haplogroups were identified in the cattle remains, confirming the presence of T1 in the earliest Spanish herds. The absence, however, of “African-derived American” haplotypes (AA/T1c1a1) in the Sevilla la Nueva sample, leaves open the origins of this sub-haplogroup in contemporary Caribbean cattle.

Possible secondary population-level effects of selective harvest of adult male muskoxen

Selective harvest regimes are often focused on males resulting in skewed sex-ratios, and for many ungulate species this strategy is sustainable. However, muskoxen (Ovibos moschatus) are very social and mature bulls (≥4 years old), particularly prime-age bulls (6–10 years old), play important roles in predator defense and recruitment. A year-round social structure incorporating large males into mixed-sex groups could make this species more susceptible to the effects of selective harvest if population composition and sex-ratios influence overall survival and reproductive success. Using detailed data collected on the muskox population occupying the Seward Peninsula, Alaska during 2002–2012, we formulated the hypothesis that the selective harvest of mature bulls may be related to documented changes in population composition and growth rates in this species. In addition, we reviewed existing published information from two other populations in Alaska, the Cape Thompson and Northeastern populations, to compare population growth rates among the three areas under differential harvest rates relative to our hypothesis. We found that on the Seward Peninsula, mature bull:adult cow ratios declined 4–12%/year and short-yearling:adult cow ratios (i.e., recruitment) declined 8–9%/year in the most heavily harvested areas. Growth rates in all 3 populations decreased disproportionately after increases in the number of bulls harvested, and calf:cow ratios declined in the Northeastern population as harvest increased. While lack of appropriate data prevented us from excluding other potential causes such as density dependent effects and changes in predator densities, our results did align with our hypothesis, suggesting that in the interest of conservation, harvest of mature males should be restricted until causal factors can be more definitively identified. If confirmed by additional research, our findings would have important implications for harvest management and conservation of muskoxen and other ungulate species with similar life-histories.

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.

Integrated allele-specific polymerase chain reaction-capillary electrophoresis microdevice for single nucleotide polymorphism genotyping

An integrated allele-specific (AS) polymerase chain reaction (PCR) and capillary electrophoresis (CE) microdevice has been developed for multiplex single nucleotide polymorphism (SNP) genotyping on a portable instrumentation, which was applied for on-site identification of HANWOO (Korean indigenous beef cattle). Twelve sets of primers were designed for targeting beef cattle's eleven SNP loci for HANWOO verification and one primer set for a positive PCR control, and the success rate for identification of HANWOO was demonstrated statistically. The AS PCR and CE separation for multiplex SNP typing was carried out on a glass-based microchip consisting of four layers: a microchannel plate for microfluidic control, a Pt-electrode plate for a resistance temperature detector (RTD), a poly(dimethylsiloxane) (PDMS) membrane and a manifold glass for microvalve function. The operation of the sample loading, AS PCR, microvalve, and CE on a chip was automated with a portable genetic analyzer, and the laser-induced fluorescence detection was performed on a miniaturized fluorescence detector. The blind samples were correctly identified as a HANWOO by showing one or two amplicon peaks in the electropherogram, while the imported beef cattle revealed more than five peaks. Our genetic analysis platform provides rapid, accurate, and on-site multiplex SNP typing.

[The effect of genetic polymorphism of the exon 2 of the beta-lactoglobulin gene on the milk composition in Chinese Holstein]

OBJECTIVE

In order to study the effect of the polymorphism at the exon2 region of the (3-LG allele gene on milk composition and yield.

METHODS

The single-strand conformation polymorphism method (PCR-SSCP) was used to analyze for polymorphism the exon2 region of the 3-LG gene (NCBI accession number: DQ489319) in Chinese Holstein.

RESULTS

Eight SSCP patterns were detected in the fragments: ab, abc, abd, abe, abcd, abce, abde and abcde, and the patterns frequencies as follows: 0.14, 0.10, 0.27, 0.23, 0.05, 0.04, 0.11 and 0.06 (P < 0.05); Six single nucleotide polymorphism (SNP) were detected in this study: sitel C>T, site2 T>C, site3 C>T, site4 C>C, site5 C> A, site6 A>T or C, and the polymorphism infonnation content (PIC) of these SNPs were in median or high polymorphism (PIC > 0.25).

CONCLUSION

These SNPs at the exon2 region of the beta-LG gene were remarkably and affected milk performance traits (milk yield, protein and fat contents) in Chinese Holstein.

Breed and trait preferences of Sheko cattle keepers in southwestern Ethiopia

Like their smallholder subsistence counterparts in developing countries, breed and trait preferences of Sheko cattle keepers have broad perspectives. Our study has documented breed and trait preferences of Sheko cattle keepers in southwestern Ethiopia--the natural breeding tract of Sheko cattle. Our results showed that due to their multifunctionality, cattle are the most preferred livestock species. Overall, farmers showed slightly more preference to local Zebus over Sheko breed. This is due to voracious feeding behavior of Sheko cattle, which make them less preferable in the face of worsening feed shortage, and due to aggressive temperament of Sheko cattle. This is despite Sheko's outperforming potential over local Zebus in their milk production, draft power, and hardiness. At trait level, overall milk production was consistently reported as the most preferred trait followed by fertility and traction. This trait preference rank has matched with the reported frequency count ranks for Sheko cattle use. However, breed preference rank has not matched with reported trait preference ranks because Sheko excels local Zebus in all the three most preferred traits, but it was ranked second. Therefore, to minimize these conflicting interests, breed management plans for Sheko cattle should target on strategies that help to solve feed shortage problem and to improve feeding practices, and on selection of less aggressive Sheko cattle. Therefore, these strategies should be considered in line with Sheko cattle conservation and genetic improvement programs.

Insights into the genetic history of French cattle from dense SNP data on 47 worldwide breeds

Background

Modern cattle originate from populations of the wild extinct aurochs through a few domestication events which occurred about 8,000 years ago. Newly domesticated populations subsequently spread worldwide following breeder migration routes. The resulting complex historical origins associated with both natural and artificial selection have led to the differentiation of numerous different cattle breeds displaying a broad phenotypic variety over a short period of time.

Methodology/Principal Findings

This study gives a detailed assessment of cattle genetic diversity based on 1,121 individuals sampled in 47 populations from different parts of the world (with a special focus on French cattle) genotyped for 44,706 autosomal SNPs. The analyzed data set consisted of new genotypes for 296 individuals representing 14 French cattle breeds which were combined to those available from three previously published studies. After characterizing SNP polymorphism in the different populations, we performed a detailed analysis of genetic structure at both the individual and population levels. We further searched for spatial patterns of genetic diversity among 23 European populations, most of them being of French origin, under the recently developed spatial Principal Component analysis framework.

Conclusions/Significance

Overall, such high throughput genotyping data confirmed a clear partitioning of the cattle genetic diversity into distinct breeds. In addition, patterns of differentiation among the three main groups of populations—the African taurine, the European taurine and zebus—may provide some additional support for three distinct domestication centres. Finally, among the European cattle breeds investigated, spatial patterns of genetic diversity were found in good agreement with the two main migration routes towards France, initially postulated based on archeological evidence.

Genotype effects on body temperature in dairy cows under grazing conditions in a hot climate including evidence for heterosis

We compared diurnal patterns of vaginal temperature in lactating cows under grazing conditions to evaluate genotype effects on body temperature regulation. Genotypes evaluated were Holstein, Jersey, Jersey x Holstein and Swedish Red x Holstein. The comparison of Holstein and Jersey versus Jersey x Holstein provided a test of whether heterosis effects body temperature regulation. Cows were fitted with intravaginal temperature recording devices that measured vaginal temperature every 15 min for 7 days. Vaginal temperature was affected by time of day (P < 0.0001) and genotype x time (P < 0.0001) regardless of whether days in milk and milk yield were used as covariates. Additional analyses indicated that the Swedish Red x Holstein had a different pattern of vaginal temperatures than the other three genotypes (Swedish Red x Holstein vs others x time; P < 0.0001) and that Holstein and Jersey had a different pattern than Jersey x Holstein [(Holstein + Jersey vs Jersey x Holstein) x time, P < 0.0001]. However, Holstein had a similar pattern to Jersey [(Holstein vs Jersey) x time, P > 0.10]. These genotype x time interactions reflect two effects. First, Swedish Red x Holstein had higher vaginal temperatures than the other genotypes in the late morning and afternoon but not after the evening milking. Secondly, Jersey x Holstein had lower vaginal temperatures than other genotypes in the late morning and afternoon and again in the late night and early morning. Results point out that there are effects of specific genotypes and evidence for heterosis on regulation of body temperature of lactating cows maintained under grazing conditions and suggest that genetic improvement for thermotolerance through breed choice or genetic selection is possible.

Karyotype of Malayan Gaur (Bos gaurus hubbacki), Sahiwal-Friesian cattle and Gaur x cattle hybrid backcrosses

Interspecific hybridization has been reported for a wide variety of vertebrate species either spontaneous or by organized crossing of bovine species. The hybrids were often carrying intermediate characters genetically and phenotypically of the parents. Thus, status information of both aspects is valuable in animal production for selection and breeding management. The Gaur-cattle hybrids was reported to be superior in production value compared to their parent cattle but fertility status was still questionable. The project was abandoned due to their fertility issue and the hybrids were kept within the cattle in a dairy farm. Cytogenetic status and breeding record of the remaining herd were unavailable since then. The herd was then translocated to a deer farm (PTH Lenggong) and kept freely in the paddock. Recently, two female calves were born via inter se mating. Peripheral blood cultures of Malayan Gaur, Sahiwal-Friesian cattle and Gaur x cattle hybrid backcrosses were analyzed via Giemsa stained metaphase. The Gaur and cattle were having diploid chromosome number (2n) of 56 and 60, respectively. Interestingly, the backcrosses from the hybrids by cattle bulls were found to have two chromosome arrangements, which are 2n = 58 and 2n = 60.

Establishment of paternity testing system using microsatellite markers in Chinese Holstein

To estimate the efficiency of microsatellite markers in paternity testing among Chinese Holstein, 30 microsatellite loci were used to differentiate 330 Chinese Holstein genotypes, according to the calculation of the allele frequency, number of alleles, effective number of alleles, genetic heterozygosity, polymorphic information content (PIC), and the exclusion probability in this cattle population. The results demonstrated that the exclusion probability ranged from 0.620 in locus BM1818 to 0.265 in locus INRA005 with the average of 0.472 and 11 microsatellite markers exceeding 0.5. The combined exclusion probability of nine microsatellite markers was over 0.99. The result showed that paternity testing of Chinese Holstein was basically resolved using the nine microsatellite markers selected.

Association of polymorphisms of leptin gene with body weight and body sizes indexes in Chinese indigenous cattle

PCR-SSCP was used to analyze the polymorphism of leptin gene in 539 samples of six cattle breeds, namely Nanyang (NY), Qinchuan (QC), Jiaxianred (JXR), Xizhen (XZ), Luxi (LX), and Holstein cow (HOL) breeds. PCR products with a 330 bp were amplified and sequenced. The results showed that the frequencies of alleles A/B of NY, QC, JXR, XZ, LX, and HOL breeds were 0.558/0.442, 0.492/0.508, 0.571/0.429, 0.658/0.342, 0.591/0.409, and 0.615/0.385, respectively. The association of variations of leptin gene with growth traits in NY, QC, JXR breeds was analyzed. Some indexes of the individuals with genotype BB were higher than that with genotype AA and AB in NY breed, such as the indexes of body length, heart length, body weight, hucklebone width, body height, and average day gain. The height at hip cross of the individuals with genotype BB was higher than that of those with genotype AA and AB in QC breed (P < 0.05). So leptin gene may be one of the candidate genes for growth traits with height at hip cross, but not for body weight, heart length, and body length trait. However, the height at hip cross and hucklebone width of the individuals with genotype AB and BB were higher than that of those with genotype AA in JXR breed (P < 0.05), but the difference was not statistically significant in body weight and body sizes (body height, body length, and heart length). And the polymorphisms in leptin gene were caused by G --> T transversion at the 66th bp position, A --> C transversion at the 67th bp position and G --> T transversion at the 299th bp position. These results may be applied to marker-assisted selection of Chinese cattle breeds.

Chemical classification of cattle. 1. Breed groups

From approximately 1000 papers with data on protein polymorphism in some 216 breeds of cattle, 10 polymorphic proteins were compared in means and variances of gene frequencies (arcsin p 1/2) for ten well-recognized breed groups for 196 of the breeds. The polymorphic proteins were alpha-lactalbumin, beta-lactoglobulin, caseins (alpha s1, beta and chi), serum albumin, transferrin, haemoglobin, amylase I and carbonic anhydrase II. The breed groups were North European, Pied Lowland, European Red brachyceros, Channel Island brachyceros, Upland brachyceros, primigenius-brachyceros mixed, primigenius, Indian Zebu, African Humped (with Zebu admixture), and African Humped (Sanga). The coherence within groups and the differences between groups are often impressive. Only carbonic anhydrase II fails to differentiate at least some of the major breed groups. In some cases paradoxical distributions of rare genetic variants can be explained by a more detailed inspection of breed history. The chemical data support the morphological and geographical divisions of cattle into major breed groups. There are three distinct but related brachyceros groups; for some polymorphisms the two Channel Island breeds, the Jersey and the Guernsey, are quite divergent. Although some authorities have considered the Pied Lowland as primigenius, it is a very distinct breed group.

Chemical classification of cattle. 2. Phylogenetic tree and specific status of the Zebu

Phylogenetic trees for the ten major breed groups of cattle were constructed by Farris's (1972) maximum parsimony method, or Fitch & Margoliash's (1967) method, which averages ou the deviation over the entire assemblage. Both techniques yield essentially identical trees. The phylogenetic tree for the ten major cattle breed groups can be superimposed on a map of Europe and western Asia, the root of the tree being close to the 'fertile crescent' in Asia Minor, believed to be a primary centre of bovine domestication. For some but not all protein variants there is a cline of gene frequencies as one proceeds from the British Isles and northwest Europe towards southeast Europe and Asia Minor, with the most extreme gene frequencies in the Zebu breeds of India. It is not clear to what extent the observed clines are primary or secondary, i.e., consequent to the initial migrations of cattle towards the end of the Pleistocene or consequent to the many migrations of man with his domesticated cattle. Such clines as exist are not in themselves sufficient to prove either selection versus genetic drift or to establish taxonomic ranking. Contrary to some suggestions in the literature, the biochemical evidence supports Linnaeus's original conclusions: Bos taurus and Bos indicus are distinct species.

Genetic variation of mitochondrial D-loop region and evolution analysis in some Chinese cattle breeds

The complete mitochondrial D-loop region from 123 individuals in 12 Chinese cattle breeds and two individuals in Germany Yellow cattle breed was sequenced and analyzed. The results were shown as follows: 93 variations and 57 haplotypes were detected, and the average number of nucleotide difference was 22.708, nucleotide diversity (d) was 0.0251 +/- 0.00479, and haplotype diversity (Hd) was 0.888 +/- 0.026, indicating very high genetic diversity in Chinese cattle breeds. In the Neighbor-Joining tree, 13 cattle breeds were divided into two main clades, Bos taurus and Bos indicus; new Clade ? had only one sequence from Apeijiaza cattle breed in Tibet, which was similar to that of yak at a higher level than other cattle breeds, proving the introgression of genes from the yak. The proportions of Bos taurus and Bos indicus were 64.3% and 35.7% respectively in Xigazê Humped cattle breed, and 50.0% and 50.0% respectively in Apeijiaza cattle breed, which revealed that Tibet cattle included Bos indicus haplotypes. The importance of Yunnan cattle in the origin of Chinese cattle was also confirmed based on their abundant haplotypes. Then, a very special haplotype i1 discovered in 27 Chinese cattle breeds, including 11 breeds in this study and 16 breeds in the GenBank, played the role of a nucleus in Chinese zebu and was further discussed. At the same time, the construction of Chinese zebu core group based on haplotype i1 validated the distinct origin of Bos indicus in Tibet, which was different from that of the other cattle breeds with zebu haplotypes in China.

Sequence variation and molecular evolution of hormone-sensitive lipase genes in species of bovidae

The partial sequences of exon I of hormone-sensitive lipase (HSL) genes in yak (Bos grunniens), cattle (Bos taurus), zebu (Bos indicus), and buffalo (Bubalus bubalis) were analyzed. Comparisons of these sequences and the deduced amino acid sequences with the homologous HSL gene and protein sequences in other mammalian species including pig (Sus scrofa), human (Homo sapiens), mouse (Mus musculus), and rat (Rattus sp.) retrieved from the GenBank were carried out and finally a phylogenetic tree was constructed using the partial DNA sequences of the HSL genes in all species. The results showed that the homologies of the partial exon I sequences of the HSL genes between yak and cattle, zebu, buffalo, pig, human, mouse, and rat were as high as 99.8%, 99.6%, 97.4%, 90.6%, 88.4%, 83.5%, and 82.3%, respectively. This was accompanied by highly homologous amino acid sequences of the HSLs: 100%, 100%, 98.2%, 94.0%, 92.2%, 89.8%, and 89.8% identity, respectively. There are more transitions, less transversions, and no insertion or deletion in variable nucleotides of the HSL genes between the yak and other species. The majority of the variable mutations was synonymous and was found most frequently at the third codon, followed by the first and second codons, a finding that was in accordance with the neutralism hypothesis for molecular evolution. In the phylogenetic tree, the cattle and zebu were clustered together first, followed by the yak, buffalo, pig, human, mouse, and rat. This was in agreement with taxonomy suggesting that the partial sequences of exon I of the HSL genes were useful in constructing the phylogenetic tree of mammalian species. Among the four species of Bovidae, genetic differentiation in the HSL genes between yak and buffalo is equivalent to that between buffalo and cattle and between buffalo and zebu. Furthermore, the genetic distances in the HSL genes are much smaller between yak, cattle, and zebu than those between each of the three species and the buffalo. Therefore, it is reasonable to consider yak as an independent species of the genus Bos.