Search

Zhoushan cattle are an endangered cattle breed in the Zhoushan islands in China. Since Zhoushan cattle have been bred in isolation, they show unique characteristics, such as dark black coat colour. However, no studies have been conducted on the genome of Zhoushan cattle. Here, we performed whole-genome sequencing of seven individuals of Zhoushan cattle and nine cattle in Wenling, geographically close to the Zhoushan islands. By integrating our data and publicly-available data, we found that Zhoushan cattle are genetically highly similar to Bos indicus cattle in south-eastern China. Furthermore, by identifying the genomic regions shared between Zhoushan cattle and Angus cattle, a Bos taurus breed, we found that the p.F195L mutation in melanocyte-stimulating hormone receptor (MC1R) could be associated with their dark black coat colour. Taken together, our results provide a valuable resource for characterising the uniqueness of Zhoushan cattle.

Background While the adoption of genomic evaluations in livestock has increased genetic gain rates, its effects on genetic diversity and accumulation of inbreeding have raised concerns in cattle populations. Increased inbreeding may affect fitness and decrease the mean performance for economically important traits, such as fertility and growth in beef cattle, with the age of inbreeding having a possible effect on the magnitude of inbreeding depression. The purpose of this study was to determine changes in genetic diversity as a result of the implementation of genomic selection in Angus cattle and quantify potential inbreeding depression effects of total pedigree and genomic inbreeding, and also to investigate the impact of recent and ancient inbreeding. Results We found that the yearly rate of inbreeding accumulation remained similar in sires and decreased significantly in dams since the implementation of genomic selection. Other measures such as effective population size and the effective number of chromosome segments show little evidence of a detrimental effect of using genomic selection strategies on the genetic diversity of beef cattle. We also quantified pedigree and genomic inbreeding depression for fertility and growth. While inbreeding did not affect fertility, an increase in pedigree or genomic inbreeding was associated with decreased birth weight, weaning weight, and post-weaning gain in both sexes. We also measured the impact of the age of inbreeding and found that recent inbreeding had a larger depressive effect on growth than ancient inbreeding. Conclusions In this study, we sought to quantify and understand the possible consequences of genomic selection on the genetic diversity of American Angus cattle. In both sires and dams, we found that, generally, genomic selection resulted in decreased rates of pedigree and genomic inbreeding accumulation and increased or sustained effective population sizes and number of independently segregating chromosome segments. We also found significant depressive effects of inbreeding accumulation on economically important growth traits, particularly with genomic and recent inbreeding.

An amendment to this paper has been published and can be accessed via the original article.

This paper seeks to evaluate the relationship between the trace element concentrations and the proximate composition of different cuts of young beef to provide helpful information to consumers in the selection of meat that could fit different nutritional requirements. Ten commercially cuts of young beef (rib boneless entrecote, tenderloin, eye round, thick flank, tail of rump, chuck tender, shin, upper chuck, flank and brisket) together with two internal muscles (diaphragm and cardiac muscle) from ten male Galician blonde calves aged approximately 9 months at slaughter and a carcass weight of 242 ± 2 kg were analyzed. A strong negative association between the main essential trace elements (Co, Cu, Fe, Mn and Se) and the protein concentration of the muscles was found, which could be related to the variable predominance of slow- or fast-twitch fibers in the different muscles. Since trace mineral concentrations in muscle are partly genetically determined and related to palatability traits, understanding the relationships between the trace element concentrations and the proximate composition could be a valuable tool for selective breeding of beef to improve the nutritional value of meat. Xunta de Galicia The Barrié Foundation

Background If unmanaged, high rates of inbreeding in livestock populations adversely impact their reproductive fitness. In beef cattle, historical selection strategies have increased the frequency of several segregating fatal autosomal recessive polymorphisms. Selective breeding has also decreased the extent of haplotypic diversity genome-wide. By identifying haplotypes for which homozygotes are not observed but would be expected based on their frequency, candidates for developmentally lethal recessive loci can be localized. This analysis comes without the need for observation of the loss-associated phenotype (e.g., failure to implant, first trimester abortion, deformity at birth). In this study, haplotypes were estimated for 3961 registered Angus individuals using 52,545 SNP loci using findhap v2, which exploited the complex pedigree among the individuals in this population. Results Seven loci were detected to possess haplotypes that were not observed in homozygous form despite a sufficiently high frequency and pedigree-based expectation of homozygote occurrence. These haplotypes were identified as candidates for harboring autosomal recessive lethal alleles. Of the genotyped individuals, 109 were resequenced to an average 27X depth of coverage to identify putative loss-of-function alleles genome-wide and had variants called using a custom in-house developed pipeline. For the candidate lethal-harboring haplotypes present in these bulls, sequence-called genotypes were used to identify concordant variants. In addition, whole-genome sequence imputation of variants was performed into the set of 3961 genotyped animals using the 109 resequenced animals to identify candidate lethal recessive variants at the seven loci. Following the imputation, no variants were identified that were fully concordant with the marker-based diplotypes. Conclusions Selective breeding programs could utilize the predicted lethal haplotypes associated with SNP genotypes. Sequencing and other methods for identifying the causal variants underlying these haplotypes can allow for more efficient methods of management such as gene editing. These two methods in total will reduce the negative impacts of inbreeding on fertility and maximize overall genetic gains.

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.

Toll-like receptor 4 (TLR4) is a receptor protein that binds pathogen ligands, which are mainly associated with Gram-negative bacteria. The objective of this study was to investigate the association of nucleotide polymorphisms in TLR4 with infectious bovine keratoconjunctivitis (IBK), or pinkeye, incidence in American Angus cattle. Animals with previously calculated breeding values for IBK susceptibility were used to identify two SNPs in TLR4; Int1 (A/G) in intron1 (-26 Ex2 position) and Ex3 (C/T) in exon3 (1,678 position). To investigate the possible role of these SNPs in IBK susceptibility, the disease incidence information was collected on 370 calves raised in Iowa at two time points-June or August (disease season) and October (at weaning) and genotyped using PCR-RFLP protocols. In statistical models including year, pasture management group, and SNP, the Int1 SNP had a significant effect on IBK infection rates both in-season (P 0.05) and at weaning (P 0.01), whereas the Ex3 SNP was not significant (P 0.79) at either time point. Furthermore, the Int1 SNP alone could account for 2.1% of phenotypic variation in IBK infection during the disease season and 3.0% of phenotypic variation in IBK infection at the time of weaning. These data indicate that there is a relationship between Int1 genotype and the rate of IBK infection in American Angus cattle.

Data on South African Angus cattle consisting of 45 259 records on weaning weight (WWT), 4 360 records on average daily gain from on-farm test (ADG-D) and 1 118 from centralized test (ADG-C) were analyzed to evaluate the effect of pre-weaning selection on estimates of genetic parameters and subsequent estimated breeding values (EBV) for post-weaning average daily gain. (Co)variance components and genetic parameters for weaning weight (WWT), ADG-C and ADG-D were estimated by REML procedures fitting three different animal models. Model 1 was a univariate model of WWT, ADG-C or ADG-D and did not account for the effect of pre-weaning selection on post-weaning ADG. Model 2 was a two-trait model of WWT and either ADG-C or ADG-D. Model 3 was a multi-trait animal model including WWT, ADG-C and ADG-D. Estimates of heritability for ADG-C were 0.39 +/- 0.080, 0.42 +/- 0.060 and 0.44 +/- 0.010 from Model 1, 2 and 3 respectively. Corresponding estimates for ADG-D were 0.18 +/- 0.020, 0.19 +/- 0.020 and 0.21 +/- 0.020 respectively. Rank correlations based on EBVs for ADG-C for all bulls were 0.92, 0.83 and 0.94 for Model 1 vs. Model 2, Model 1 vs. Model 3, and Model 2 vs. Model 3, respectively and they indicated a possible re-ranking of bulls when including or excluding a correlated pre-weaning trait. Rank correlations for ADG-D also followed a similar trend. Inclusion of pre-weaning information in genetic analysis for post-weaning average daily gain is necessary to account for selection at weaning.

The aim of the present study was twofold: first, to design a panel of 96 sires that reflects the breadth of genetic diversity in U.S. beef cattle, and second, to use this panel to discover nucleotide sequence diversity and haplotype structures of interleukin (IL)-8 in commercial populations. The latter is a requisite for epidemiological studies designed to test whether IL8 alleles are risk factors for acquiring or maintaining bacterial infections in production environments. IL-8 encodes a proinflammatory cytokine that plays a central role in cell-mediated immunity by attracting and activating neutrophils in the early stages of host defense against bacterial invasion. Seven single-nucleotide polymorphism (SNP) markers were identified by sequencing two IL8 DNA segments amplified from the panel of 17 popular cattle breeds (MARC beef cattle diversity panel, version 2.1). Assays for automated genotype scoring by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) were developed to independently verify the seven SNP alleles in the 96 bulls and 313 cattle from the MARC reference population. Five haplotype structures, spanning the two IL8 DNA segments, were unambiguously defined for the set of seven IL8 SNPs. Based on the breadth of germplasm in bovine diversity panel, the five haplotype structures for IL8 are estimated to represent >98% of those present in these DNA segments in commercial populations of U.S. beef cattle. The frequencies of the five respective haplotypes in the eight Angus sires of the diversity panel (0.75, 0.25, 0.00, 0.00, 0.00) were similar to those scored in 150 purebred Angus cattle from six herds in four Midwestern states (0.82, 0.18, 0.01, 0.00, 0.00), suggesting that the diversity panel may also be useful for estimating allele frequencies in commercial populations.

Background Infectious Bovine Keratoconjunctivitis (IBK) in beef cattle, commonly known as pinkeye, is a bacterial disease caused by M o r a x e l l a bovis. IBK is characterized by excessive tearing and ulceration of the cornea. Perforation of the cornea may also occur in severe cases. IBK is considered the most important ocular disease in cattle production, due to the decreased growth performance of infected individuals and its subsequent economic effects. IBK is an economically important, lowly heritable categorical disease trait. Mass selection of unaffected animals has not been successful at reducing disease incidence. Genome-wide studies can determine chromosomal regions associated with IBK susceptibility. The objective of the study was to detect single-nucleotide polymorphism (SNP) markers in linkage disequilibrium (LD) with genetic variants associated with IBK in American Angus cattle. Results The proportion of phenotypic variance explained by markers was 0.06 in the whole genome analysis of IBK incidence classified as two, three or nine categories. Whole-genome analysis using any categorisation of (two, three or nine) IBK scores showed that locations on chromosomes 2, 12, 13 and 21 were associated with IBK disease. The genomic locations on chromosomes 13 and 21 overlap with QTLs associated with Bovine spongiform encephalopathy, clinical mastitis or somatic cell count. Conclusions Results of these genome-wide analyses indicated that if the underlying genetic factors confer not only IBK susceptibility but also IBK severity, treating IBK phenotypes as a two-categorical trait can cause information loss in the genome-wide analysis. These results help our overall understanding of the genetics of IBK and have the potential to provide information for future use in breeding schemes.

miRNAs are a class of small, single-stranded, non-coding RNAs that perform post-transcriptional repression of target genes by binding to 3’ untranslated regions. Research has found that miRNAs involved in the regulation of many metabolic processes. Here we uncovered that the beef quality of Angus cattle sharply diversified after acute stress. By performing miRNA microarray analysis, 13 miRNAs were significantly differentially expressed in stressed group compared to control group. Using a bioinformatics method, 135 protein-coding genes were predicted as the targets of significant differentially expressed miRNAs. Gene Ontology (GO) term and Ingenuity Pathway Analysis (IPA) mined that these target genes involved in some important pathways, which may have impact on meat quality and beef tenderness.

Background Molecular estimates of breeding value are expected to increase selection response due to improvements in the accuracy of selection and a reduction in generation interval, particularly for traits that are difficult or expensive to record or are measured late in life. Several statistical methods for incorporating molecular data into breeding value estimation have been proposed, however, most studies have utilized simulated data in which the generated linkage disequilibrium may not represent the targeted livestock population. A genomic relationship matrix was developed for 698 Angus steers and 1,707 Angus sires using 41,028 single nucleotide polymorphisms and breeding values were estimated using feed efficiency phenotypes (average daily feed intake, residual feed intake, and average daily gain) recorded on the steers. The number of SNPs needed to accurately estimate a genomic relationship matrix was evaluated in this population. Results Results were compared to estimates produced from pedigree-based mixed model analysis of 862 Angus steers with 34,864 identified paternal relatives but no female ancestors. Estimates of additive genetic variance and breeding value accuracies were similar for AFI and RFI using the numerator and genomic relationship matrices despite fewer animals in the genomic analysis. Bootstrap analyses indicated that 2,500-10,000 markers are required for robust estimation of genomic relationship matrices in cattle. Conclusions This research shows that breeding values and their accuracies may be estimated for commercially important sires for traits recorded in experimental populations without the need for pedigree data to establish identity by descent between members of the commercial and experimental populations when at least 2,500 SNPs are available for the generation of a genomic relationship matrix.

Background Linkage disequilibrium (LD) maps can provide a wealth of information on specific marker-phenotype relationships, especially in areas of the genome where positional candidate genes with similar functions are located. A recently published high resolution radiation hybrid map of bovine chromosome 14 (BTA14) together with the bovine physical map have enabled the creation of more accurate LD maps for BTA14 in both dairy and beef cattle. Results Over 500 Single Nucleotide Polymorphism (SNP) markers from both Angus and Holstein animals had their phased haplotypes estimated using GENOPROB and their pairwise r2 values compared. For both breeds, results showed that average LD extends at moderate levels up to 100 kilo base pairs (kbp) and falls to background levels after 500 kbp. Haplotype block structure analysis using HAPLOVIEW under the four gamete rule identified 122 haplotype blocks for both Angus and Holstein. In addition, SNP tagging analysis identified 410 SNPs and 420 SNPs in Holstein and Angus, respectively, for future whole genome association studies on BTA14. Correlation analysis for marker pairs common to these two breeds confirmed that there are no substantial correlations between r-values at distances over 10 kbp. Comparison of extended haplotype homozygosity (EHH), which calculates the LD decay away from a core haplotype, shows that in Holstein there is long range LD decay away from the DGAT1 region consistent with the selection for milk fat % in this population. Comparison of EHH values for Angus in the same region shows very little long range LD. Conclusion Overall, the results presented here can be applied in future single or haplotype association analysis for both populations, aiding in confirming or excluding potential polymorphisms as causative mutations, especially around Quantitative Trait Loci regions. In addition, knowledge of specific LD information among markers will aid the research community in selecting appropriate markers for whole genome association studies.