Association study between SNP markers located in meat quality candidate genes with intramuscular fat content in an endangered dual-purpose cattle population

The aim of this study was to associate single nucleotide polymorphisms (SNP) of the bovine calcium-activated neutral protease µ-calpain, calpastatin, diacylglycerol-O-acyltransferase, adipose fatty acid binding protein, retinoic acid receptor-related orphan receptor C (RORC), and thyroglobulin (TG) gene with intramuscular fat content (IMF). Therefore, 542 animals of the cattle breed "Rotes Höhenvieh" (RHV) were phenotyped for IMF. Genotyping of the animals was performed using polymerase chain reaction-restriction fragment length polymorphism tests for six SNP from candidate genes for meat quality traits. In addition, we calculated allele substitution and dominance effects on IMF. A subgroup of animals (n = 44, reduced dataset) with extraordinary high IMF was analyzed separately. The mean IMF content was 2.5% (SD: 2.8) but ranged from 0.02% to 23.9%, underlining the breeds' potential for quality meat production. Allele and genotype frequencies for all SNP were similar in the complete and reduced dataset. Association analyses in the complete dataset revealed the strongest effects of RORC on IMF (P = 0.075). The log-transformed least-squares mean for IMF of genotype g.3290GG was 0.45 ± 0.16, 0.26 ± 0.14 for genotype g.3290GT, and 0.32 ± 0.14 for genotype g.3290TT. In the reduced dataset, we found a significant effect (P < 0.05) of the g.422C>T-SNP of TG on IMF, with highest IMF for genotype CT (0.91 ± 0.17), lowest IMF for genotype TT (0.37 ± 0.25), and medium IMF for genotype CC (0.59 ± 0.16; log-transformed values). Compared to the complete dataset, allele substitution effects increased in the reduced dataset for most of the SNP, possibly due to the selective genotyping strategy, with focus on animals with highest IMF implying strong phenotypic IMF contrast. Dominance effects were small in both datasets, related to the high heritability of IMF. Results indicated RHV breed particularities regarding the effects of meat quality genes on IMF. An explanation might be the breeding history of RHV with focus on adaptation and resilience in harsh outdoor systems. Consequently, it is imperative to develop breed-specific selection strategies. Allele substitution and dominance effects were in a similar direction in both datasets, suggesting the same breeding approaches for different RHV strains in different regions. Nevertheless, a selective genotyping approach (reduced dataset), contributed to more pronounced genotype effect differences on IMF and dominance values.

Association of polymorphisms in lipid and energy metabolism-related genes with fattening performance in Simmental cattle

Lipid and energy metabolism are major constituents of mammal growth and thus fattening performance of cattle. This study was designed to evaluate the effects of polymorphisms in lipid and energy metabolism-related genes including oxidized low-density lipoprotein receptor 1 (OLR1), lactoferrin (LTF), stearoyl-CoA desaturase (SCD), beta-lactoglobulin (LGB), thyroglobulin (TG), annexin A9 (ANXA9), myogenic factor 5 (MYF5), protein kinase AMP-activated non-catalytic subunit gamma 3 (PRKAG3), and pituitary-specific transcriptional factor 1 (PIT1), on fattening performance in Simmental cattle. A total of 72 purebred Simmental bulls with a similar initial age and weight were fattened on the same farm for 10 months. Association analysis was performed using linear mixed models. The OLR1 marker was significantly associated with the final weight (FW), hot carcass weight (HCW), chilled carcass weight (CCW), dressing percentage (DP), and total weight gain (TWG). SCD affected the FW, TWG, and average daily live weight gain (ADWG). The present results clearly demonstrated the significant impact of the TG marker on fattening performance. It was highly significantly associated with the FW, HCW, CCW, and TWG. The SCD × TG and the OLR1 × TG interactions had remarkable effects on the traits analyzed. The GACC and CCCC haplotypes of the SCD × TG and OLR1 × TG, respectively, were found to be powerful markers for fattening performance in Simmentals. Novel associations in this study may be useful for further genetic evaluations to improve beef cattle breeding.

Polymorphism of Genes and Their Impact on Beef Quality

The single-nucleotide polymorphism (SNP) form of genes is a valuable source of information regarding their suitability for use as specific markers of desirable traits in beef cattle breeding. For several decades, breeding work focused on improving production efficiency through optimizing the feed conversion ratio and improving daily gains and meat quality. Many research teams previously undertook research work on single-nucleotide polymorphism in myostatin (MSTN), thyroglobulin (TG), calpain (CAPN), and calpastatin (CAST) proteins. The literature review focuses on the most frequently addressed issues concerning these genes in beef cattle production and points to a number of relevant studies on the genes' polymorphic forms. The four genes presented are worth considering during breeding work as a set of genes that can positively influence productivity and production quality.

Genetic parameters and genome-wide association study of digital cushion thickness in Holstein cows

The digital cushion is linked to the development of claw horn lesions (CHL) in dairy cattle. The objectives of this study were to (1) estimate genetic parameters for digital cushion thickness (DCT), (2) estimate the genetic correlation between DCT and CHL, and (3) identify candidate genes associated with DCT. A cohort of 2,352 Holstein dairy cows were prospectively enrolled on 4 farms and assessed at 4 time points: before calving, immediately after calving, in early lactation, and in late lactation. At each time point, CHL was recorded by veterinary surgeons, and ultrasonographic images of the digital cushion were stored and retrospectively measured at 2 anatomical locations. Animals were genotyped and pedigree details extracted from the national database. Genetic parameters were estimated following a single-step approach implemented in AIREMLF90. Four traits were analyzed: the 2 DCT measurements, sole lesions (sole hemorrhage and sole ulcers), and white line lesions. All traits were analyzed with univariate linear mixed models; bivariate models were fit to estimate the genetic correlation between traits within and between time points. Single-marker and window-based genome-wide association analyses of DCT traits were conducted at each time point; candidate genes were mapped near (<0.2 Mb) or within the genomic markers or windows with the largest effects. Heritability estimates of DCT ranged from 0.14 to 0.44 depending on the location of DCT measurement and assessment time point. The genetic correlation between DCT and sole lesions was generally negative, notably between DCT immediately after calving and sole lesions in early or late lactation, and between DCT in early or late lactation and sole lesion severity in early or late lactation. Digital cushion thickness was not genetically correlated with white line lesions. A polygenic background to DCT was found; genes associated with inflammation, fat metabolism, and bone development were mapped near or within the top markers and windows. The moderate heritability of DCT provides an opportunity to use selective breeding to change DCT in a population. The negative genetic correlation between DCT and sole lesions at different stages of production lends support to current hypotheses of sole lesion pathogenesis. Highlighted candidate genes provide information regarding the complex genetic background of DCT in Holstein cows, but further studies are needed to explore and corroborate these findings.

Weighted single-step GWAS for body mass index and scans for recent signatures of selection in Yorkshire pigs

Controlling extra fat deposition is economically favorable in modern swine industry. Understanding the genetic architecture of fat deposition traits such as body mass index (BMI) can help in improving genomic selection for such traits. We utilized a weighted single-step genome-wide association study (WssGWAS) to detect genetic regions and candidate genes associated with BMI in a Yorkshire pig population. Three extended haplotype homozygosity (EHH)-related statistics were also incorporated within a de-correlated composite of multiple signals (DCMS) framework to detect recent selection signatures signals. Overall, the full pedigree consisted of 7,016 pigs, of which 5,561 had BMI records and 598 pigs were genotyped with an 80 K single nucleotide polymorphism (SNP) array. Results showed that the most significant windows (top 15) explained 9.35% of BMI genetic variance. Several genes were detected in regions previously associated with pig fat deposition traits and treated as potential candidate genes for BMI in Yorkshire pigs: FTMT, SRFBP1, KHDRBS3, FOXG1, SOD3, LRRC32, TSKU, ACER3, B3GNT6, CCDC201, ADCY1, RAMP3, TBRG4, CCM2. Signature of selection analysis revealed multiple candidate genes previously associated with various economic traits. However, BMI genetic variance explained by regions under selection pressure was minimal (1.31%). In conclusion, candidate genes associated with Yorkshire pigs’ BMI trait were identified using WssGWAS. Gene enrichment analysis indicated that the identified candidate genes were enriched in the insulin secretion pathway. We anticipate that these results further advance our understanding of the genetic architecture of BMI in Yorkshire pigs and provide information for genomic selection for fat deposition in this breed.

Associations between gene polymorphisms and selected meat traits in cattle - A review

Maintaining a high level of beef consumption requires paying attention not only to quantitative traits but also to the quality and dietary properties of meat. Growing consumer demands do not leave producers many options for how animals are selected for breeding and animal keeping. Meat and carcass fatness quality traits, which are influenced by multiple genes, are economically important in beef cattle breeding programs. The recent availability of genome sequencing methods and many previously identified molecular markers offer new opportunities for animal breeding, including the use of molecular information in selection programs. Many gene polymorphisms have thus far been analyzed and evaluated as potential candidates for molecular markers of meat quality traits. Knowledge of these markers can be further applied to breeding programs through marker-assisted selection. In this literature review, we discuss the most promising and well-described candidates and their associations with selected beef production traits.

Effect of the bovine TG5 gene polymorphism on milk- and meat-producing ability

Aim:

This study aimed to determine the effect of thyroglobulin (TG5) gene polymorphism on milk and meat productivity in the various cattle breeds currently bred in the Republic of Bashkortostan.

Materials and Methods:

The test was performed on dairy cattle of Black-and-White, Bestuzhev, and Simmental breeds, and meat cattle of Hereford and limousine breeds. The purpose of the test was to search for associations between the polymorphic alleles of the thyroglobulin (TG5) gene and economically useful traits.

Results:

All studied breeds showed a frequency predominance of the TG5^C allele (from 0.56 to 0.71). A clear trend of an effect of the genotypes of the TG5 gene on milk-productivity indicators was revealed; cows with the TG5^TT genotype have the highest milk yield and fat content in milk. The milk of cows of Bestuzhev and Simmental breeds that possessed this genotype was also characterized by higher protein content.

Conclusion:

We identified an effect of the polymorphism of the TG5 gene in the Hereford and limousine breeds on fat metabolism intensity indicators, such as fat output and fat content, in the longissimus muscle and in the general sample of ground beef.

Comprehensive assessment of candidate genes associated with fattening performance in Holstein-Friesian bulls

The objective of this study was to determine the association of single nucleotide polymorphisms (SNPs) in selected candidate genes with fattening performance traits in a commercial cattle herd. Fifteen SNPs in 12 candidate genes (LEP, FABP4, DGAT1, TG, IGF1, IGF1R, MYF5, LGB, CAPN1, CAST, GHR, and OLR1) were evaluated in 296 purebred Holstein-Friesian bulls using PCR-RFLP (polymerase chain reaction - restriction fragment length polymorphism). Associations between each segregating SNP and genetic merit for fattening performance were quantified using linear mixed models. Traits included in the study were fattening period, final weight, dry matter intake, feed conversion rate, and average daily weight gain. Apart from the general determination of the above-mentioned traits, each trait was evaluated based on the fattening periods between five selected target body weights (W1  =  100 kg, W2  =  200 kg, W3  =  300 kg, W4  =  400 kg, W5  =  450 kg). All markers with the exception of CAPN1 530, IGF1R, TG, and DGAT1 were associated with at least one of the traits. Furthermore, novel associations were observed for LEP  ×  GHR, IGF1  ×  LEP, FABP4 3691  ×  FABP4 2834, and FAP4 3533  ×  LEP interactions. The results of this study confirm some previously reported associations. Moreover, novel associations have been identified, which may be incorporated into breeding programs to improve fattening performance.

Transformation of animal genomics by next-generation sequencing technologies: a decade of challenges and their impact on genetic architecture

For more than a quarter of a century, sequencing technologies from Sanger's method to next-generation high-throughput techniques have provided fascinating opportunities in the life sciences. The continuing upward trajectory of sequencing technologies will improve livestock research and expedite the development of various new genomic and technological studies with farm animals. The use of high-throughput technologies in livestock research has increased interest in metagenomics, epigenetics, genome-wide association studies, and identification of single nucleotide polymorphisms and copy number variations. Such studies are beginning to provide revolutionary insights into biological and evolutionary processes. Farm animals, such as cattle, swine, and horses, have played a dual role as economically and agriculturally important animals as well as biomedical research models. The first part of this study explores the current state of sequencing methods, many of which are already used in animal genomic studies, and the second part summarizes the state of cattle, swine, horse, and chicken genome sequencing and illustrates its achievements during the last few years. Finally, we describe several high-throughput sequencing approaches for the improved detection of known, unknown, and emerging infectious agents, leading to better diagnosis of infectious diseases. The insights from viral metagenomics and the advancement of next-generation sequencing will strongly support specific and efficient vaccine development and provide strategies for controlling infectious disease transmission among animal populations and/or between animals and humans. However, prospective sequencing technologies will require further research and in-field testing before reaching the marketplace.

Genome-wide association study for feed efficiency and growth traits in U.S. beef cattle

Background

Single nucleotide polymorphism (SNP) arrays for domestic cattle have catalyzed the identification of genetic markers associated with complex traits for inclusion in modern breeding and selection programs. Using actual and imputed Illumina 778K genotypes for 3887 U.S. beef cattle from 3 populations (Angus, Hereford, SimAngus), we performed genome-wide association analyses for feed efficiency and growth traits including average daily gain (ADG), dry matter intake (DMI), mid-test metabolic weight (MMWT), and residual feed intake (RFI), with marker-based heritability estimates produced for all traits and populations.

Results

Moderate and/or large-effect QTL were detected for all traits in all populations, as jointly defined by the estimated proportion of variance explained (PVE) by marker effects (PVE ≥ 1.0%) and a nominal P-value threshold (P ≤ 5e-05). Lead SNPs with PVE ≥ 2.0% were considered putative evidence of large-effect QTL (n = 52), whereas those with PVE ≥ 1.0% but < 2.0% were considered putative evidence for moderate-effect QTL (n = 35). Identical or proximal lead SNPs associated with ADG, DMI, MMWT, and RFI collectively supported the potential for either pleiotropic QTL, or independent but proximal causal mutations for multiple traits within and between the analyzed populations. Marker-based heritability estimates for all investigated traits ranged from 0.18 to 0.60 using 778K genotypes, or from 0.17 to 0.57 using 50K genotypes (reduced from Illumina 778K HD to Illumina Bovine SNP50). An investigation to determine if QTL detected by 778K analysis could also be detected using 50K genotypes produced variable results, suggesting that 50K analyses were generally insufficient for QTL detection in these populations, and that relevant breeding or selection programs should be based on higher density analyses (imputed or directly ascertained).

Conclusions

Fourteen moderate to large-effect QTL regions which ranged from being physically proximal (lead SNPs ≤ 3Mb) to fully overlapping for RFI, DMI, ADG, and MMWT were detected within and between populations, and included evidence for pleiotropy, proximal but independent causal mutations, and multi-breed QTL. Bovine positional candidate genes for these traits were functionally conserved across vertebrate species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3754-y) contains supplementary material, which is available to authorized users.

Estimates of epistatic and pleiotropic effects of () and () genetic markers on beef heifer performance traits enhanced by selection

Genetic marker effects and type of inheritance are estimated with poor precision when minor marker allele frequencies are low. A stable composite population (MARC II) was subjected to marker assisted selection for 2 yr to equalize and genetic marker frequencies to evaluate the epistatic and pleiotropic effects of these markers on BW, reproduction, and first calf performance traits in replacement beef females ( = 171) managed under 2 postweaning development protocols. Traits evaluated on the heifers were birth BW, weaning BW, 11-mo BW, 12-mo BW, 13-mo BW, first breeding season pregnancy evaluation BW, first calving season BW, 11-mo puberty, 12-mo puberty, 13-mo puberty, first breeding season pregnancy, and first calf weaning rate. Additionally, heifer's first calf performance traits of ordinal calving date, first calf birth BW, and first calf weaning BW (with and without age adjustment) were analyzed. Selection to increase minor allele frequencies and balanced sampling across genotype classes enhanced the ability to detect all genetic effects except dominance × dominance epistasis. The × genotype effect was significant ( < 0.05) for 11-mo BW and 12-mo BW and tended to be significant ( = 0.08) for 13-mo BW. Consistently, for all 3 traits, the most significant effect among epistatic × genotype effects was the additive effect, with the G allele decreasing BW. There were no associations between × genotype and fertility related traits ( ≥ 0.46) in this study. Additionally, there were no × genotype associations with first progeny performance traits ( ≥ 0.14). The large effect of the additive × additive interaction on first calf weaning BW was imprecisely estimated, which may warrant further investigation.

Exploring evidence of positive selection reveals genetic basis of meat quality traits in Berkshire pigs through whole genome sequencing

BACKGROUND

Natural and artificial selection following domestication has led to the existence of more than a hundred pig breeds, as well as incredible variation in phenotypic traits. Berkshire pigs are regarded as having superior meat quality compared to other breeds. As the meat production industry seeks selective breeding approaches to improve profitable traits such as meat quality, information about genetic determinants of these traits is in high demand. However, most of the studies have been performed using trained sensory panel analysis without investigating the underlying genetic factors. Here we investigate the relationship between genomic composition and this phenotypic trait by scanning for signatures of positive selection in whole-genome sequencing data.

RESULTS

We generated genomes of 10 Berkshire pigs at a total of 100.6 coverage depth, using the Illumina Hiseq2000 platform. Along with the genomes of 11 Landrace and 13 Yorkshire pigs, we identified genomic variants of 18.9 million SNVs and 3.4 million Indels in the mapped regions. We identified several associated genes related to lipid metabolism, intramuscular fatty acid deposition, and muscle fiber type which attribute to pork quality (TG, FABP1, AKIRIN2, GLP2R, TGFBR3, JPH3, ICAM2, and ERN1) by applying between population statistical tests (XP-EHH and XP-CLR). A statistical enrichment test was also conducted to detect breed specific genetic variation. In addition, de novo short sequence read assembly strategy identified several candidate genes (SLC25A14, IGF1, PI4KA, CACNA1A) as also contributing to lipid metabolism.

CONCLUSIONS

Results revealed several candidate genes involved in Berkshire meat quality; most of these genes are involved in lipid metabolism and intramuscular fat deposition. These results can provide a basis for future research on the genomic characteristics of Berkshire pigs.

Effect of thyroglobulin gene polymorphisms on growth, carcass composition and meat quality traits in Chinese beef cattle

The thyroglobulin (TG) gene has been studied as an important gene related to fat deposition, since not only does TG gene locate in a quantitative trait locus with an effect on fat deposition, but also it encodes the precursor of thyroid hormones which have crucial biological functions in energy metabolism. In the present study, we identified four novel SNPs at the 5' flanking region of the bovine TG gene. Association analysis indicated that the G275A, G277C, G280A and C281G SNPs were significantly associated with average daily gain (ADG, P < 0.01 for G275A and G277C, P < 0.05 for G280A and C281G). Five haplotypes for the 4 SNPs were constructed and their effects on growth, carcass composition and meat quality traits were evaluated. The results showed no significant effect of haplotype on ADG. Meanwhile, no significant association was found between 4 SNPs and other growth, carcass composition and meat quality traits including intramuscular fat. Bioinformatics analysis showed that 4 SNPs may results in potential transcription factor binding site changes. Results of this study suggest that TG gene-specific SNPs may be a useful marker for growth traits in marker assisted selection programs in beef cattle.

Association analysis of polymorphism in thyroglobulin gene promoter with milk production traits in riverine buffalo (Bubalus bubalis)

Polymorphism within the promoter region of bovine thyroglobulin has been reported to be associated with milk and meat quality. In this study, we investigated the genetic variation within thyroglobulin promoter region of swamp and riverine buffaloes using PCR-SSCP technique and sequencing, and also analyzing association of polymorphism with the milk production traits. The study revealed four conformational patterns, A, B, C, and D among 323 buffaloes of two riverine breeds and different swamp populations. The frequency of SSCP variant 'A' was found to be invariably high among all buffalo populations. Variant 'C' was found to be absent in pure swamp population and present with higher frequency among riverine dairy breeds Mehsana and Nili Ravi. Frequency of D variant was observed to be highest in buffalo population, representing riverine and hybrid types. Sequencing of three representative PCR products of each of the SSCP variants, revealed three polymorphic sites responsible, 33C > T, 176G > T and 221C > T, in the buffalo TG promoter region. Further, association studies of SSCP variants with various milk production and milk quality traits indicated significant effect on fat percentage in buffaloes belonging to Mehsana and Nili Ravi dairy breeds. The preliminary results also showed the substantial variations in the distribution of SSCP variants' frequencies across swamp and riverine buffaloes, two distinct populations being reared for meat and milk production, respectively.

Effect of polymorphisms in candidate genes on carcass and meat quality traits in double muscled Piemontese cattle

The aim of this study was to investigate the association between 10 candidate genes and carcass weight and conformation, carcass daily gain, and meat quality (pH, color, cooking loss, drip loss and shear force) in 990 double-muscled Piemontese young bulls. Animals were genotyped at each of the following genes: growth hormone, growth hormone receptor, pro-opiomelanocortin, pro-opiomelanocortin class 1 homeobox 1, melanocortin-4 receptor, corticotrophin-releasing hormone, diacylglycerol O-acyltransferase-1, thyroglobulin, carboxypeptidase E and gamma-3 regulatory subunit of the AMP-activated protein kinase. All the investigated SNPs had additive effects which were relevant for at least one of the traits. Relevant associations between the investigated SNPs and carcass weight, carcass daily gain and carcass conformation were detected, whereas associations of SNPs with meat quality were moderate. Results confirmed some of previously reported associations, but diverged for others. Validation in other cattle breeds is required to use these SNPs in gene-assisted selection programs for enhancement of carcass traits and meat quality.

Selection for genetic markers in beef cattle reveals complex associations of thyroglobulin and casein1-s1 with carcass and meat traits

Genetic markers in casein (CSN1S1) and thyroglobulin (TG) genes have previously been associated with fat distribution in cattle. Determining the nature of these genetic associations (additive, recessive, or dominant) has been difficult, because both markers have small minor allele frequencies in most beef cattle populations. This results in few animals homozygous for the minor alleles. selection to increase the frequencies of the minor alleles for 2 SNP markers in these genes was undertaken in a composite population. The objective was to obtain better estimates of genetic effects associated with these markers and determine if there were epistatic interactions. Selection increased the frequencies of minor alleles for both SNP from <0.30 to 0.45. Bulls (n = 24) heterozygous for both SNP were used in 3 yr to produce 204 steer progeny harvested at an average age of 474 d. The combined effect of the 9 CSN1S1 × TG genotypes was associated with carcass-adjusted fat thickness (P < 0.06) and meat tenderness predicted at the abattoir by visible and near-infrared reflectance spectroscopy (P < 0.04). Genotype did not affect BW from birth through harvest, ribeye area, marbling score, slice shear force, or image-based yield grade (P > 0.10). Additive, dominance, and epistatic SNP association effects were estimated from genotypic effects for adjusted fat thickness and predicted meat tenderness. Adjusted fat thickness showed a dominance association with TG SNP (P < 0.06) and an epistatic additive CSN1S1 × additive TG association (P < 0.03). For predicted meat tenderness, heterozygous TG meat was more tender than meat from either homozygote (P < 0.002). Dominance and epistatic associations can result in different SNP allele substitution effects in populations where SNP have the same linkage disequilibrium with causal mutations but have different frequencies. Although the complex associations estimated in this study would contribute little to within-population selection response, they could be important for marker-assisted management or reciprocal selection schemes.

[Artificial selection for cattle based on high-density SNP markers]

With the implementation of genetic improvement in recent years, artificial selection has greatly improved beef cattle production performance and its genetic basis has been dramatically changed. In this study, based on the Illumina BovineSNP50 (54K) and BovineHD (770K) BeadChip and the FST value, we analyzed the genetic differentiation of cattle and screened the imprints of selection in bovine genome. Finally, we found 47104 OUTLIER SNP loci and 3064 candidate genes, for example, CLIC5, TG, CACNA2D1, and FSHR etc. The biological processes and molecular functions of genes were analyzed through gene annotation.The results of this study established a genome-wide map of selection footprints in beef cattle genome and a clue for in-depth study of artificial selection and understanding of biological evolution.Our results indicate that artificial selection has played an important role in cattle breed genetic improvement.

Haplotype analysis improved evidence for candidate genes for intramuscular fat percentage from a genome wide association study of cattle

In genome wide association studies (GWAS), haplotype analyses of SNP data are neglected in favour of single point analysis of associations. In a recent GWAS, we found that none of the known candidate genes for intramuscular fat (IMF) had been identified. In this study, data from the GWAS for these candidate genes were re-analysed as haplotypes. First, we confirmed that the methodology would find evidence for association between haplotypes in candidate genes of the calpain-calpastatin complex and musculus longissimus lumborum peak force (LLPF), because these genes had been confirmed through single point analysis in the GWAS. Then, for intramuscular fat percent (IMF), we found significant partial haplotype substitution effects for the genes ADIPOQ and CXCR4, as well as suggestive associations to the genes CEBPA, FASN, and CAPN1. Haplotypes for these genes explained 80% more of the phenotypic variance compared to the best single SNP. For some genes the analyses suggested that there was more than one causative mutation in some genes, or confirmed that some causative mutations are limited to particular subgroups of a species. Fitting the SNPs and their interactions simultaneously explained a similar amount of the phenotypic variance compared to haplotype analyses. Haplotype analysis is a neglected part of the suite of tools used to analyse GWAS data, would be a useful method to extract more information from these data sets, and may contribute to reducing the missing heritability problem.

Genetic association of phosphodiesterase 1B (PDE1B) with carcass traits in Korean cattle

Quantitative trait loci for fat deposition and carcass traits have been identified in the vicinity of the gene encoding phosphodiesterase 1B (PDE1B) on bovine chromosome 5. Therefore, the PDE1B gene can be considered as a positional and functional candidate gene for carcass traits in beef cattle. This study aimed to identify single nucleotide polymorphisms (SNPs) in the PDE1B gene and to evaluate their associations with carcass traits in Korean cattle. Eight SNPs, g.440T>G, g.17122A>G, g.17507A>C, g.17575A>G, g.17607T>C, g.17609C>A, g.17692C>T, and g.17707C>G, were identified in the region ranged from exon 1 to intron 6. Five of them were used for association analysis because of their availability of restriction fragment length polymorphisms. As a result, g.17122A>G in intron 3 was significantly associated with backfat thickness (BFT), and g.17507A>C in exon 5 was associated with longissimus dorsi muscle area (LMA, P < 0.05). Animals with the AG genotype of g.17122 had thicker BFT than those with the AA genotype. Animals with the AA or AC genotype of g.17507A>C had larger LMA than those with the CC genotype. We suggested the PDE1B gene as a candidate gene for carcass traits of beef cattle. Fine mapping would be required for application to marker-assisted selection.

A genome scan for quantitative trait loci influencing carcass, post-natal growth and reproductive traits in commercial Angus cattle

To gain insight into the number of loci of large effect that underlie variation in cattle, a quantitative trait locus (QTL) scan for 14 economically important traits was performed in two commercial Angus populations using 390 microsatellites, 11 single nucleotide polymorphisms (SNPs) and one duplication loci. The first population comprised 1769 registered Angus bulls born between 1955 and 2003, with Expected Progeny Differences computed by the American Angus Association. The second comprised 38 half-sib families containing 1622 steers with six post-natal growth and carcass phenotypes. Linkage analysis was performed by half-sib least squares regression with gridqtl or Bayesian Markov chain Monte Carlo analysis of complex pedigrees with loki. Of the 673 detected QTL, only 118 have previously been reported, reflecting both the conservative approach to QTL reporting in the literature, and the more liberal approach taken in this study. From 33 to 71% of the genetic variance and 35 to 56% of the phenotypic variance in each trait was explained by the detected QTL. To analyse the effects of 11 SNPs and one duplication locus within candidate genes on each trait, a single marker analysis was performed by fitting an additive allele substitution model in both mapping populations. There were 53 associations detected between the SNP/duplication loci and traits with -log(10) P(nominal) ≥ 4.0, where each association explained 0.92% to 4.4% of the genetic variance and 0.01% to 1.86% of the phenotypic variance. Of these associations, only six SNP/duplication loci were located within 8 cM of a QTL peak for the trait, with two being located at the QTL peak: SST_DG156121:c.362A>G for ribeye muscle area and TG_X05380:c.422C>T for calving ease. Strong associations between several SNP/duplication loci and trait variation were obtained in the absence of any detected linked QTL. However, we reject the causality of several commercialized DNA tests, including an association between TG_X05380:c.422C>T and marbling in Angus cattle.

Association of CAPN1 316, CAPN1 4751 and TG5 markers with bovine meat quality traits in Mexico

We examined allele and genotype frequencies for the molecular markers CAPN1 316, CAPN1 4751 and TG5, and determined whether they are associated with beef quality traits in Mexican cattle. One hundred and twenty-four longissimus dorsi muscle samples were collected from cattle from north, central and southern Mexico. CAPN1 316 and CAPN1 4751 frequencies were determined using the allelic discrimination assay and the TG5 marker was typed by PCR-RFLP. Meat quality traits included intramuscular fat content (IMF) and tenderness determined by Warner-Bratzler shear force (WBSF) at 24 h postmortem. The association test was made using a mixed model, including genotypes, genetic group, and sampling location as fixed effects. Least squares means and significant interactions were compared using least significant differences based on the mixed procedure. CAPN1 316 CC was found at a low frequency (0.03) and has been reported as a favorable genotype associated with tenderness meat. Genotype frequencies for CAPN1 4751 were similar in favorable (CC) and unfavorable (TT) genotypes (0.26 and 0.28, respectively). The TG5 CC genotype had a frequency of 0.73, while the TT genotype frequency was 0.01. The means for WBSF and IMF were 4.08 ± 1.35 kg and 5.23 ± 2.14%, respectively. Sampling site and the CAPN1 316 genotypes significantly affected WBSF (P < 0.05). Samples collected from Hermosillo, Sonora, had the lowest WBSF (P < 0.05), while those collected in Veracruz were toughest (WBSF = 5.267 kg). The effect of GG and TG5 genotypes on IMF was significant (P < 0.05). CAPN1 316 and TG5 markers were found to be significantly associated with beef quality traits and thus will be useful for Mexican beef characterization.

Association of thyroglobulin gene variants with carcass and meat quality traits in beef cattle

Thyroid hormones play an important role in regulating metabolism and can affect homeostasis of fat deposition. The gene encoding thyroglobulin (TG), producing the precursor for thyroid hormones, has been proposed as a positional and functional candidate gene for a QTL with an effect on fat deposition. In the present study, we identified 6 novel SNPs at the 3' flanking region of the TG gene. The SNP marker association analysis indicated that the T354C, G392A, A430G and T433G SNP markers were significantly associated with marbling score (P < 0.05). Animals with the new homozygote genotype had higher marbling score than those with the other genotypes. Otherwise, the linkage disequilibrium analysis indicated that these four SNPs were completely linked (r (2) = 1). Results from this study suggest that TG gene-specific SNP may be a useful marker for meat quality traits in future marker assisted selection programs in beef cattle.

Association analysis of single nucleotide polymorphisms in DGAT1, TG and FABP4 genes and intramuscular fat in crossbred Bos taurus cattle

Previous studies have indicated that single nucleotide polymorphisms (SNPs) in the diacylglycerol-O-acyltransferase1 (DGAT1), thyroglobulin (TG) and adipose fatty acid binding protein (FABP4) genes are associated with intramuscular fat (IMF) levels or marbling scores in beef. The objectives were to estimate the frequency of SNPs in these candidate genes in purebred Irish cattle (n=459) and to determine if individual SNPs are associated with IMF values of longissimus thoracis et lumborum (LTL) and semimembranosus (SM) muscle of crossbred animals (n=138). Results indicated no significant association between the SNPs and IMF, despite the power of this study being sufficient to detect an association with SNPs in the DGAT1 and FABP4 genes. The results confirm the lack of an association found by many other studies and suggest that these SNPs are not influential on the divergent IMF levels in the crossbred population tested.

Hyaluronan: is it a biomarker for adipose development within bovine muscle?

Based on an association with extracellular matrix remodelling, mitosis, proliferation and adipogenic differentiation, the glycosaminoglycan hyaluronan (HA) was assessed as a marker for intramuscular fat (IMF) development (marbling) in bovine loin muscle (longissimus dorsi, LD). Loin samples collected from the quartering site of feedlot-finished Wagyu–Angus and Jersey–Limousin steers were assayed for percentage IMF (IMF%) and HA after assignment of AUS-MEAT marbling scores. There was a moderate phenotypic correlation (r2 = 0.69) between IMF% and marbling score but little variance was explained by HA concentration. Breed was not a significant factor in marbling score or IMF% but did influence the HA concentration of the LD, with Wagyu–Angus steers having 2-fold more HA than Jersey–Limousin steers at the same marbling score. The non-linear decline in fat-adjusted HA levels as marbling score increased suggests that HA concentration was associated with lean growth potential of the muscle rather than adipogenesis. Using a different experimental approach, differences in distribution and amount of HA could not be discerned in histological sections of LD from age-matched Wagyu–Hereford heifers allocated to a low (score 1) or medium (score 3) marbling score group. These findings were consistent with the absence of differences between the two groups for other indicators of fatness (IMF% and P8 fat depth), maturity and myofibre characteristics despite an increase in oxidative capacity of the muscle with age. The data support the conclusion that the concentration of HA in the LD alone was not predictive of development of intramuscular fat.

Comparison of skeletal muscle transcriptional profiles in dairy and beef breeds bulls

A cDNA microarray (18 263 probes) was used for transcriptome analysis of bovine skeletal muscle (m. semitendinosus) in 12-month-old bulls of the beef breed Limousin (LIM) and the typical dairy breed Holstein-Friesian (HF, used as a reference). We aimed to identify the genes whose expression may reflect the muscle phenotype of beef bulls. A comparison of muscle transcriptional profiles revealed significant differences in expression of 393 genes between HF and LIM. We classified biological functions of 117 genes with over 2-fold differences in expression between the examined breeds. Among them, 72 genes were up-regulated and 45 genes were down-regulated in LIM vs. HF. The genes were involved in protein metabolism and modifications (22 genes), signal transduction (15), nucleoside, nucleotide and nucleic acid metabolism (13), cell cycle (9), cell structure and motility (9), developmental processes (9), intracellular protein traffic (7), cell proliferation and differentiation (6), cell adhesion (6), lipid, fatty acid and steroid metabolism (5), transport (5), and other processes. For the purpose of microarray data validation, we randomly selected 4 genes: trip12, mrps30, pycrl, and c-erbb3. Real-time RT-PCR results showed similar trends in gene expression changes as those observed in microarray studies. Basing on results of the present study, we proposed a model of the regulation of skeletal muscle growth and differentiation, with a principal role of the somatotropic pathway. It may explain at least in part the development of muscle phenotype in LIM bulls. We assume that the growth hormone directly or indirectly (through IGF-1) activates the calcium-signaling pathway with calcineurin, which stimulates myogenic regulatory factors (MRFs) and inhibits early growth response gene. The inhibition results in indirect activation of MRFs and impaired activation of TGF-beta1 and myostatin, which finally facilitates terminal muscle differentiation.

A splice site single nucleotide polymorphism of the fatty acid binding protein 4 gene appears to be associated with intramuscular fat deposition in longissimus muscle in Australian cattle

Fatty acid binding protein 4 (FABP4) is a candidate gene affecting fatness traits of mammals. However, its association with fatness traits in cattle and other livestock species is not consistent from one study to another. Here, we sequenced the coding sequence of FABP4 looking for non-synonymous variants. We identified a splice site mutation between the third exon and the third intron of bovine FABP4. We genotyped this SNP, FABP4:g.2502C>G, in 1409 cattle with intramuscular fat measurements from seven breeds. The average allele frequency of the C allele was 0.66 with a range of 0.45 to 0.85. A regression on the number of G alleles shows a statistically significant effect of alpha = 0.11, P = 0.044. This appears to confirm an association between IMF and variation at FABP4, with an effect of 0.3% of the variation in our sample when using this SNP.

Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds

The imprints of domestication and breed development on the genomes of livestock likely differ from those of companion animals. A deep draft sequence assembly of shotgun reads from a single Hereford female and comparative sequences sampled from six additional breeds were used to develop probes to interrogate 37,470 single-nucleotide polymorphisms (SNPs) in 497 cattle from 19 geographically and biologically diverse breeds. These data show that cattle have undergone a rapid recent decrease in effective population size from a very large ancestral population, possibly due to bottlenecks associated with domestication, selection, and breed formation. Domestication and artificial selection appear to have left detectable signatures of selection within the cattle genome, yet the current levels of diversity within breeds are at least as great as exists within humans.

Variation at CPE but not CEBPA appears to be associated with intramuscular fat deposition in the longissimus muscle of cattle

An important step in the localisation of quantitative trait loci is the confirmation of trait-marker associations in independent studies. In this report, we test three single nucleotide polymorphisms (SNP) of two genes for associations to intramuscular fat (IMF) measurements in cattle. We genotyped SNP of carboxypeptidase E (CPE) and ccaat/enhancer binding protein, α (CEBPA) in a sample of a total of 813 cattle of taurine, composite and indicine breeds. All three polymorphisms showed significant differences between breeds, with the widest range found in CEBPA:g.271A > C where the A allele frequency ranged from P = 0.07 in Brahman to 0.88 in Shorthorn. The taurine breeds showed high linkage disequilibrium between the pair of CPE SNP, with all four breeds showing r2 = 1.0. The Brahman and Santa Gertrudis showed r2 ≤ 0.17. Both CPE:g.445C > T and CPE:g.601C > T SNP showed significant allele substitution effects to IMF in animals of taurine ancestry, with an allele substitution effect of α = 0.22, P = 0.020 for CPE:g.445C > T, explaining 0.4% of the phenotypic variance.

Bovine gene polymorphisms related to fat deposition and meat tenderness

Leptin, thyroglobulin and diacylglycerol O-acyltransferase play important roles in fat metabolism. Fat deposition has an influence on meat quality and consumers' choice. The aim of this study was to determine allele and genotype frequencies of polymorphisms of the bovine genes, which encode leptin (LEP), thyroglobulin (TG) and diacylglycerol O-acyltransferase (DGAT1). A further objective was to establish the effects of these polymorphisms on meat characteristics. We genotyped 147 animals belonging to the Nelore (Bos indicus), Canchim (5/8 Bos taurus + 3/8 Bos indicus), Rubia Gallega X Nelore (1/2 Bos taurus + 1/2 Bos indicus), Brangus Three-way cross (9/16 Bos taurus + 7/16 Bos indicus) and Braunvieh Three-way cross (3/4 Bos taurus + 1/4 Bos indicus) breeds. Backfat thickness, total lipids, marbling score, ribeye area and shear force were fitted, using the General Linear Model (GLM) procedure of the SAS software. The least square means of genotypes and genetic groups were compared using Tukey's test. Allele frequencies vary among the genetic groups, depending on Bos indicus versus Bos taurus influence. The LEP polymorphism segregates in pure Bos indicus Nelore animals, which is a new finding. The T allele of TG is fixed in Nelore, and DGAT1 segregates in all groups, but the frequency of allele A is lower in Nelore animals. The results showed no association between the genotypes and traits studied, but a genetic group effect on these traits was found. So, the genetic background remains relevant for fat deposition and meat tenderness, but the gene markers developed for Bos taurus may be insufficient for Bos indicus.

A genome map of divergent artificial selection between Bos taurus dairy cattle and Bos taurus beef cattle

A number of cattle breeds have become highly specialized for milk or beef production, following strong artificial selection for these traits. In this paper, we compare allele frequencies from 9323 single nucleotide polymorphism (SNP) markers genotyped in dairy and beef cattle breeds averaged in sliding windows across the genome, with the aim of identifying divergently selected regions of the genome between the production types. The value of the method for identifying selection signatures was validated by four sources of evidence. First, differences in allele frequencies between dairy and beef cattle at individual SNPs were correlated with the effects of those SNPs on production traits. Secondly, large differences in allele frequencies generally occurred in the same location for two independent data sets (correlation 0.45) between sliding window averages. Thirdly, the largest differences in sliding window average difference in allele frequencies were found on chromosome 20 in the region of the growth hormone receptor gene, which carries a mutation known to have an effect on milk production traits in a number of dairy populations. Finally, for the chromosome tested, the location of selection signatures between dairy and beef cattle was correlated with the location of selection signatures within dairy cattle.

Association analysis of thyroglobulin gene variants with carcass and meat quality traits in beef cattle

Thyroid hormones play an important role in regulating the metabolism and can affect the homeostasis of fat deposition. The gene encoding thyroglobulin (TG), producing the precursor for thyroid hormones, has been proposed as a positional and functional candidate gene for a QTL with an effect on fat deposition. In the present study, we identified 6 novel SNPs at the 3' flanking region of the TG gene. The SNP marker association analysis indicated that the SNP markers G133C, G156A, C220T and A506C were significantly associated with marbling score (P<0.05, N=271). Animals with the new homozygote genotype had higher marbling scores than those with the other genotypes. Besides, the linkage disequilibrium analysis indicated that these 4 SNPs were completely linked (r2 = 1). Results of this study suggest that the TG-gene-specific SNP may be a useful marker for meat quality traits in future marker-assisted selection programmes in beef cattle.

Identification of genetic markers for fat deposition and meat tenderness on bovine chromosome 5: development of a low-density single nucleotide polymorphism map

As genetic markers, SNP are well suited for the development of genetic tests for production traits in livestock. They are stable through many generations and can provide direct assessment of individual animal's genetic merit if they are in linkage disequilibrium and phase with functional genetic variation. Bovine chromosome 5 has been shown to harbor genetic variation affecting production traits in multiple cattle populations; thus, this chromosome was targeted for SNP-based marker development and subsequent association analysis with carcass and growth phenotypes. Discovery of SNP was performed in a panel of 16 sires representing two sires from each of seven beef breeds and two Holstein sires by PCR amplification and sequencing using primers designed from genomic sequence obtained by low-coverage sequencing of bacterial artificial chromosome (BAC) clones. From 550 SNP, 296 (54%) were tentatively identified as having a minor allele frequency >10%. Forty-five SNP derived from 15 BAC were chosen based on minor allele frequency and were genotyped in 564 steers and their sires. Production and carcass data were collected on the steers as a part of the Germplasm Evaluation (GPE), Cycle VII Project at the U.S. Meat Animal Research Center (Clay Center, NE), which involves of the evaluation of sires from seven of the most popular U.S. breeds. Haplotypes based on seven SNP derived from a BAC containing the bovine genes HEM1 and PDE1B were associated with traits related to carcass fat. Steers homozygous for the major haplotype had 0.15 +/- 0.04 cm less subcutaneous fat, 0.57 +/- 0.18 kg less rib fat, 0.18 +/- 0.07 lower yield grade, 1.11 +/- 0.35% less predicted fat yield, and 0.79 +/- 0.3% greater predicted retail product yield than heterozygotes. The frequency of the major haplotype was 0.70 in the steers, and it ranged from 0.44 (Limousin) to 0.98 (Simmental and Gelbvieh) in a panel consisting of an average of 20 purebred sires from each of the seven breeds. A second set of haplotypes based on four SNP derived from a BAC containing the genes NOL1 and CHD4 was associated with Warner-Bratzler shear force. Steers homozygous for the major haplotype had 0.27 +/- 0.11 kg greater shear force than those heterozygous for the major haplotype and one of two minor haplotypes. The frequency of the major haplotype was 0.59 in the steers and ranged from 0.27 (Hereford) to approximately 0.95 (Angus and Red Angus) in the panel of purebred sires. These results demonstrate the feasibility of targeting QTL regions for SNP-based marker development and that a low level of coverage can identify markers associated with phenotypic traits.

Assessing the association of single nucleotide polymorphisms at the thyroglobulin gene with carcass traits in beef cattle

The objective of this study was to assess the association of SNP in the thyroglobulin gene, including a previously reported marker in current industry use, with marbling score in beef cattle. Three populations, designated GPE6, GPE7, and GPE8, were studied. The GPE6 population sampled breeds that could be used as alternative germplasm sources in beef cattle production, including Wagyu, Swedish Red and White, Friesian, and Norwegian Red. The GPE7 population sampled 7 popular beef cattle breeds used in temperate climates of the United States: Angus, Charolais, Gelbvieh, Hereford, Limousin, Red Angus, and Simmental. The GPE8 population sampled Bos indicus-influenced breeds used in subtropical regions of the country and subtropical and tropical regions of the world, including Beefmaster, Bonsmara, Brangus, and Romosinuano. Evaluation of 6 SNP in the thyroglobulin gene, including 5 newly described variations, showed no association (P > 0.10) with marbling score in these populations, except a tendency (P < 0.10) for an association with the previously described marker in GPE6. Closer examination of the GPE6 data revealed that the source of the tendency was an association (P < 0.02) with marbling in animals of Wagyu inheritance. Animals having Wagyu background and inheriting the TT genotype had a greater marbling score (599 +/- 20) than those inheriting the CC (540 +/- 10) or the CT (541 +/- 11) genotype. No association was detected with any other carcass trait for this marker in the 3 populations. Furthermore, none of the 5 newly described markers in the gene displayed an association with marbling score. The data indicate that markers at the thyroglobulin gene may be a useful predictor of marbling performance for producers raising Wagyu-based cattle. Although associations with marbling score in the remaining populations were not large or significant, the TT genotype had the numerically greatest marbling score in each population.