Effects of the F94L myostatin gene mutation in beef × dairy crossed cattle on muscle fiber type, live performance, carcass characteristics, and boxed beef and retail cut yields

Producer live performance data and carcasses from steers (n = 116) resulting from the mating of four Limousin/Angus sires heterozygous for the F94L myostatin mutation to Jersey/Holstein dams were utilized to evaluate the effects of one copy of the F94L allele on live performance, carcass traits and USDA grades, and boxed beef and retail yields. Slaughter data were collected at time of harvest and carcass data were collected 48 hours postmortem. One side from each of the 58 carcasses was fabricated into boxed beef and retail cuts by experienced lab personnel 5-8 d postmortem. One copy of the F94L allele did not affect gestation length, birth weight, percentage of unassisted births, feedlot average daily gain, live weight at harvest, hot carcass weight, or dressing percentage (P > 0.05). Muscle fiber analysis indicated that the increase in muscularity by the F94L allele in the semitendinosus and longissimus was likely due to hyperplasia as there was a 19% increase in the quantity of myosin heavy chain type IIA and IIX fibers in the semitendinosus (P < 0.05) with no effect on muscle fiber size (P > 0.05). Carcasses from steers with one F94L allele had larger ribeye areas (99.2 vs. 92.3 sq.cm.), greater ribeye width:length ratios (0.498 vs. 0.479), lower USDA yield grades (2.21 vs. 2.66), and lower marbling scores (438 vs. 480) (P < 0.05). Additionally, for boxed beef yields, one F94L allele, vs. zero F94L alleles, increased (P < 0.05) 85/15 trimmings (+0.59%), top round (+0.28%), strip loin (+0.12%), eye round (+0.11%), tenderloin (+0.07%), boneless foreshank (+0.07%), cap/wedge (+0.06%), and tri-tip (+0.04%). Overall, carcasses from steers with one F94L allele had a greater boxed beef yield (+1.06%), boxed beef plus 85/15 trimmings yield (+1.65%), and total retail cuts plus ground beef 85/15 yield (+1.78%) than carcasses from steers with zero F94L alleles (P < 0.05). One copy of the F94L allele utilized in beef-on-dairy breeding system had no significant impact on live performance traits but resulted in lower marbling scores and increased muscularity as evidenced through larger, more beef-shaped ribeyes, lower USDA yield grades, and greater carcass cutout yields (both boxed beef and retail yields).

The Blonde d'Aquitaine T3811>G3811 mutation in the myostatin gene: association with growth, carcass, and muscle phenotypes in veal calves

The mutation T3811 → G3811 (TG3811) discovered in the myostatin gene of the Blonde d'Aquitaine breed is suspected of contributing to the outstanding muscularity of this breed. An experiment was designed to estimate the effect of this mutation in an F2 and back-cross Blonde d'Aquitaine × Holstein population. By genotyping all known mutations in the myostatin gene, it was ensured that the TG3811 mutation was indeed the only known mutation segregating in this population. Fifty-six calves (43 F2, 13 back-cross) were intensively fattened and slaughtered at 24.0 ± 1.4 wk of age. The effects of the mutation were estimated by comparing the calves with the [T/T] (n = 18), [T/G] (n = 30), and [G/G] (n = 8) genotypes. Highly significant substitution effects (P < 0.001), above + 1.2 phenotypic SD, were shown on carcass yield and muscularity scores. Birth weight (P < 0.001) was positively affected by the mutation (+0.8 SD) but not growth rate (P = 0.97), while carcass length (P = 0.03), and fatness (P ≤ 0.03) were negatively affected (-0.5 to -0.7 SD). The characteristics of the Triceps brachii muscle were affected by the mutation (P < 0.001), with lower ICDH activity (oxidative) and a higher proportion of myosin type 2X muscle fibers (fast twitch). The effects of the TG3811 mutation were similar to those of other known myostatin mutations, although the Blonde d'Aquitaine animals, which are predominantly [G/G] homozygous, do not exhibit extreme double muscling.

Detection of selection signatures in Limousin cattle using whole-genome resequencing

Limousin, a renowned beef breed originating from central France, has been selectively bred over the last 100 years to improve economically important traits. We used whole-genome sequencing data from 10 unrelated Limousin bull calves to detect polymorphisms and identify regions under selection. A total of 13 943 766 variants were identified. Moreover, 311 852 bi-allelic SNPs and 92 229 indels located on autosomes were fixed for the alternative allele in all sequenced animals, including the previously reported missense deleterious F94L mutation in MSTN. We performed a whole-genome screen to discover genomic regions with excess homozygosity, using the pooled heterozygosity score and identified 171 different candidate selective sweeps. In total, 68 candidate genes were found in only 57 of these regions, indicating that a large fraction of the genome under selection might lie in non-coding regions and suggesting that a majority of adaptive mutations might be regulatory in nature. Many QTL were found within candidate selective sweep regions, including QTL associated with shear force or carcass weight. Among the putative selective sweeps, we located genes (MSTN, NCKAP5, RUNX2) that potentially contribute to important phenotypes in Limousin. Several candidate regions and genes under selection were also found in previous genome-wide selection scans performed in Limousin. In addition, we were able to pinpoint candidate causative regulatory polymorphisms in GRIK3 and RUNX2 that might have been under selection. Our results will contribute to improved understanding of the mechanisms and targets of artificial selection and will facilitate the interpretation of GWASs performed in Limousin.

Enhanced estimates of carcass and meat quality effects for polymorphisms in myostatin and µ-calpain genes

The objective of this study was to enhance estimates of additive, dominance, and epistatic effects of marker polymorphisms on beef carcass and quality traits. Myostatin (MSTN) F94L SNP and the µ-calpain (CAPN1) 316 and 4751 SNP haplotype have previously been associated with fat and muscle traits in beef cattle. Multiyear selection in a composite population segregating these polymorphisms increased minor allele (F94L L) and chosen haplotype (CAPN1 CC and GT) frequencies to intermediate levels resulting in more precise estimates of additive and nonadditive genetic effects. During the 3 yr after selection, 176 steers were evaluated for growth, carcass, meat quality, tenderness (n = 103), and meat color traits. The statistical model included year, age of dam, age of the steer, and genotype in a random animal model. The 9 genotypes (3 CAPN1 diplotypes × 3 F94L genotypes) affected marbling score, ribeye area, adjusted fat thickness, vision yield grade (all P < 0.001), slice shear force (P = 0.03), and CIE L* reflectance (P = 0.01). Linear contrasts of the 9 genotypes estimated additive, recessive, and epistatic genetic effects. Significant additive effects of the F94L L allele decreased marbling score, adjusted fat thickness, vision yield grade, and slice shear force; and increased ribeye area and CIE L* reflectance. The homozygous F94L FF and LL genotypes differed by 1.3 to 1.9 phenotypic SD for most carcass traits and by 0.8 to 0.9 SD for slice shear force and CIE L* reflectance but carcass weight differed by only 3 kg (0.1 SD). The L allele was partially recessive to F for ribeye area (P = 0.02) and the heterozygous FL means tended to be closer to the FF genotype than the LL genotype for other carcass traits but differences from additive were not significant. The CAPN1 additive × F94L additive effect on slice shear force was the only significant epistatic estimate. The F94L L allele is prevalent in Limousin but nearly absent in other U.S. purebreds. This allele had about half of the effects on birth weight, muscle, and fat traits reported for severe MSTN mutations in Belgian Blue and Piedmontese breeds. The interaction between MSTN and CAPN1 genotypes may reflect the strong additive effects of MSTN F94L L allele on fat and muscle traits interfering with the phenotypic effect of CAPN1 genotype on meat tenderness.

The Possible Role of Complete Loss of Myostatin in Limiting Excessive Proliferation of Muscle Cells (C2C12) via Activation of MicroRNAs

Myostatin (MSTN) is a member of the TGF-β superfamily that negatively regulates skeletal muscle growth and differentiation. However, the mechanism by which complete MSTN deletion limits excessive proliferation of muscle cells remains unclear. In this study, we knocked out MSTN in mouse myoblast lines using a Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR/Cas9) system and sequenced the mRNA and miRNA transcriptomes. The results show that complete loss of MSTN upregulates seven miRNAs targeting an interaction network composed of 28 downregulated genes, including TGFB1, FOS and RB1. These genes are closely associated with tumorigenesis and cell proliferation. Our study suggests that complete loss of MSTN may limit excessive cell proliferation via activation of miRNAs. These data will contribute to the treatment of rhabdomyosarcoma (RMS).

Characterization of a molt-related myostatin gene (FmMstn) from the banana shrimp Fenneropenaeus merguiensis

Myostatin is an important member of the transforming growth factor (TGF) family that functions to regulate muscle growth in animals. In this study, the myostatin gene (FmMstn) and two slightly different (short and long forms) cDNAs of the banana shrimp Fenneropenaeus merguiensis were cloned and characterized. Similar to Mstn gene of the scallop, fish and mammal, FmMstn gene consists of 3 exons and 2 introns. The 2kb upstream promoter region of the FmMstn gene consists of putative response elements for myocyte enhancing factor (MEF2) and E-box factors. The longest open reading frame of the short Mstn consists of 1260bp encoding for a protein with 420 amino acid residues. The long FmMstn is almost identical to the short FmMstn with the exception of 8 amino acid insertions. FmMstn is most similar to the Mstn of Litopenaeus vannamei and Penaeus monodon sharing >92-98% amino acid sequence identity. Multiple sequence alignment results revealed high degree of amino acid conservation of the cysteine residues and mature peptide of the FmMstn with Mstn from other animals. FmMstn transcript was detected in the heart, muscle, optic nerve and thoracic ganglion. FmMstn transcript level in muscle is higher in early postmolt, decreases in intermolt and increases again towards ecdysis. Higher expression level of FmMstn is also observed in smaller shrimp of the same age. Knock-down of FmMstn gene by RNAi can cause a significant increase in molt cycle duration and failure of some shrimp to undergo ecdysis. Direct DNA sequencing results revealed that FmMstn gene is highly polymorphic and several potential SNPs have been identified. Some SNPs are associated with the size difference of the shrimp. In summary, the result of this study indicates that shrimp FmMstn gene is molt/growth-related and the presence of SNP suggests that it could be a candidate gene for shrimp genetic improvement research.

Comparison of breeding value prediction for two traits in a Nellore-Angus crossbred population using different Bayesian modeling methodologies

The objectives of this study were to 1) compare four models for breeding value prediction using genomic or pedigree information and 2) evaluate the impact of fixed effects that account for family structure. Comparisons were made in a Nellore-Angus population comprising F₂, F₃ and half-siblings to embryo transfer F₂ calves with records for overall temperament at weaning (TEMP; n = 769) and Warner-Bratzler shear force (WBSF; n = 387). After quality control, there were 34,913 whole genome SNP markers remaining. Bayesian methods employed were BayesB (π̃ = 0.995 or 0.997 for WBSF or TEMP, respectively) and BayesC (π = 0 and π̃), where π̃ is the ideal proportion of markers not included. Direct genomic values (DGV) from single trait Bayesian analyses were compared to conventional pedigree-based animal model breeding values. Numerically, BayesC procedures (using π̃) had the highest accuracy of all models for WBSF and TEMP (ρ̂_gĝ = 0.843 and 0.923, respectively), but BayesB had the least bias (regression of performance on prediction closest to 1, β̂_y,x = 2.886 and 1.755, respectively). Accounting for family structure decreased accuracy and increased bias in prediction of DGV indicating a detrimental impact when used in these prediction methods that simultaneously fit many markers.

Large-effect pleiotropic or closely linked QTL segregate within and across ten US cattle breeds

Background

The availability of high-density SNP assays including the BovineSNP50 (50 K) enables the identification of novel quantitative trait loci (QTL) and improvement of the resolution of the locations of previously mapped QTL. We performed a series of genome-wide association studies (GWAS) using 50 K genotypes scored in 18,274 animals from 10 US beef cattle breeds with observations for twelve body weights, calving ease and carcass traits.

Results

A total of 159 large-effects QTL (defined as 1-Mb genome windows explaining more than 1% of additive genetic variance) were identified. In general, more QTL were identified in analyses with bigger sample sizes. Four large-effect pleiotropic or closely linked QTLs located on BTA6 at 37–42 Mb (primarily at 38 Mb), on BTA7 at 93 Mb, on BTA14 at 23–26 Mb (primarily at 25 Mb) and on BTA20 at 4 Mb were identified in more than one breed. Several breed-specific large-effect pleiotropic or closely linked QTL were also identified. Some identified QTL regions harbor genes known to have large effects on a variety of traits in cattle such as PLAG1 and MSTN and others harbor promising candidate genes including NCAPG, ARRDC3, ERGIC1, SH3PXD2B, HMGA2, MSRB3, LEMD3, TIGAR, SEPT7, and KIRREL3. Gene ontology analysis revealed that genes involved in ossification and in adipose tissue development were over-represented in the identified pleiotropic QTL. Also, the MAPK signaling pathway was identified as a common pathway affected by the genes located near the pleiotropic QTL.

Conclusions

This largest GWAS ever performed in beef cattle, led us to discover several novel across-breed and breed-specific large-effect pleiotropic QTL that cumulatively account for a significant percentage of additive genetic variance (e.g. more than a third of additive genetic variance of birth and mature weights; and calving ease direct in Hereford). These results will improve our understanding of the biology of growth and body composition in cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-442) contains supplementary material, which is available to authorized users.

Hepatic transcriptome profiling identifies differences in expression of genes associated with changes in metabolism and postnatal growth between Hereford and Holstein-Friesian bulls

This study examined liver transcriptomic profiles of cattle distinctly different in meat and milk production capacity. It was performed on bulls of two different genetic backgrounds: Herefords (H), a meat breed, and Holstein-Friesians (HF), a dairy breed. Using bovine long oligo-microarrays and qPCR, we identified 128 genes that are differentially expressed between the two breeds. In H bulls, we observed up-regulation of genes involved in fatty acid biosynthesis and lipid metabolism (CD36, CAT, HSD3B1, FABP1, ACAA1) and involved in insulin signaling (INSR, INSIG2, NR4A1) and down-regulation of genes involved in somatotropic axis signaling (IGF1, GHR, IGFBP3) as compared to HF. Transcriptome profiling of these two breeds allowed us to pinpoint the transcriptional differences between Holstein and Hereford bulls at hepatic level associated with changes in metabolism and postnatal growth.

Gene-based single nucleotide polymorphism discovery in bovine muscle using next-generation transcriptomic sequencing

Background

Genetic information based on molecular markers has increasingly being used in cattle breeding improvement programmes, as a mean to improve conventionally phenotypic selection. Advances in molecular genetics have led to the identification of several genetic markers associated with genes affecting economic traits. Until recently, the identification of the causative genetic variants involved in the phenotypes of interest has remained a difficult task. The advent of novel sequencing technologies now offers a new opportunity for the identification of such variants. Despite sequencing costs plummeting, sequencing whole-genomes or large targeted regions is still too expensive for most laboratories. A transcriptomic-based sequencing approach offers a cheaper alternative to identify a large number of polymorphisms and possibly to discover causative variants. In the present study, we performed a gene-based single nucleotide polymorphism (SNP) discovery analysis in bovine Longissimus thoraci, using RNA-Seq. To our knowledge, this represents the first study done in bovine muscle.

Results

Messenger RNAs from Longissimus thoraci from three Limousin bull calves were subjected to high-throughput sequencing. Approximately 36–46 million paired-end reads were obtained per library. A total of 19,752 transcripts were identified and 34,376 different SNPs were detected. Fifty-five percent of the SNPs were found in coding regions and ~22% resulted in an amino acid change. Applying a very stringent SNP quality threshold, we detected 8,407 different high-confidence SNPs, 18% of which are non synonymous coding SNPs. To analyse the accuracy of RNA-Seq technology for SNP detection, 48 SNPs were selected for validation by genotyping. No discrepancies were observed when using the highest SNP probability threshold. To test the usefulness of the identified SNPs, the 48 selected SNPs were assessed by genotyping 93 bovine samples, representing mostly the nine major breeds used in France. Principal component analysis indicates a clear separation between the nine populations.

Conclusions

The RNA-Seq data and the collection of newly discovered coding SNPs improve the genomic resources available for cattle, especially for beef breeds. The large amount of variation present in genes expressed in Limousin Longissimus thoracis, especially the large number of non synonymous coding SNPs, may prove useful to study the mechanisms underlying the genetic variability of meat quality traits.

Association of an ACSL1 gene variant with polyunsaturated fatty acids in bovine skeletal muscle

Background

The intramuscular fat deposition and the fatty acid profiles of beef affect meat quality. High proportions of unsaturated fatty acids are related to beef flavor and are beneficial for the nutritional value of meat. Moreover, a variety of clinical and epidemiologic studies showed that particularly long-chain omega-3 fatty acids from animal sources have a positive impact on human health and disease.

Results

To screen for genetic factors affecting fatty acid profiles in beef, we initially performed a microsatellite-based genome scan in a F₂ Charolais × German Holstein resource population and identified a quantitative trait locus (QTL) for fatty acid composition in a region on bovine chromosome 27 where previously QTL affecting marbling score had been detected in beef cattle populations. The long-chain acyl-CoA synthetase 1 (ACSL1) gene was identified as the most plausible functional and positional candidate gene in the QTL interval due to its direct impact on fatty acid metabolism and its position in the QTL interval. ACSL1 is necessary for synthesis of long-chain acyl-CoA esters, fatty acid degradation and phospholipid remodeling. We validated the genomic annotation of the bovine ACSL1 gene by in silico comparative sequence analysis and experimental verification. Re-sequencing of the complete coding, exon-flanking intronic sequences, 3' untranslated region (3'UTR) and partial promoter region of the ACSL1 gene revealed three synonymous mutations in exons 6, 7, and 20, six noncoding intronic gene variants, six polymorphisms in the promoter region, and four variants in the 3' UTR region. The association analysis identified the gene variant in intron 5 of the ACSL1 gene (c.481-233A>G) to be significantly associated with the relative content of distinct fractions and ratios of fatty acids (e.g., n-3 fatty acids, polyunsaturated, n-3 long-chain polyunsaturated fatty acids, trans vaccenic acid) in skeletal muscle. A tentative association of the ACSL1 gene variant with intramuscular fat content indicated that an indirect effect on fatty acid composition via modulation of total fat content of skeletal muscle cannot be excluded.

Conclusions

The initial QTL analysis suggested the ACSL1 gene as a positional and functional candidate gene for fatty acid composition in bovine skeletal muscle. The findings of subsequent association analyses indicate that ACSL1 or a separate gene in close proximity might play a functional role in mediating the lipid composition of beef.

Quantitative trait loci for organ weights and adipose fat composition in Jersey and Limousin back-cross cattle finished on pasture or feedlot

A QTL study of live animal and carcass traits in beef cattle was carried out in New Zealand and Australia. Back-cross calves (385 heifers and 398 steers) were generated, with Jersey and Limousin backgrounds. This paper reports on weights of eight organs (heart, liver, lungs, kidneys, spleen, gastro-intestinal tract, fat, and rumen contents) and 12 fat composition traits (fatty acid (FA) percentages, saturated and monounsaturated FA subtotals, and fat melting point). The New Zealand cattle were reared and finished on pasture, whilst Australian cattle were reared on grass and finished on grain for at least 180 days. For organ weights and fat composition traits, 10 and 12 significant QTL locations (P<0.05), respectively, were detected on a genome-wide basis, in combined-sire or within-sire analyses. Seven QTL significant for organ weights were found at the proximal end of chromosome 2. This chromosome carries a variant myostatin allele (F94L), segregating from the Limousin ancestry, and this is a positional candidate for the QTL. Ten significant QTL for fat composition were found on chromosomes 19 and 26. Fatty acid synthase and stearoyl-CoA desaturase (SCD1), respectively, are positional candidate genes for these QTL. Two FA QTL found to be common to sire groups in both populations were for percentages of C14:0 and C14:1 (relative to all FAs) on chromosome 26, near the SCD1 candidate gene.

Genetic mapping of quantitative trait loci for meat quality and muscle metabolic traits in cattle

A whole-genome scan was carried out in New Zealand and Australia to detect quantitative trait loci (QTL) for live animal and carcass composition traits and meat quality attributes in cattle. Backcross calves (385 heifers and 398 steers) were generated, with Jersey and Limousin backgrounds. The New Zealand cattle were reared and finished on pasture, whilst Australian cattle were reared on grass and finished on grain for at least 180 days. This paper reports on meat quality traits (tenderness measured as shear force at 4-5 ages on two muscles as well as associated traits of meat colour, pH and cooking loss) and a number of metabolic traits. For meat quality traits, 18 significant QTL (P < 0.05), located in nine linkage groups, were detected on a genome-wise basis, in combined-sire (seven QTL) or within-sire analyses (11 QTL). For metabolic traits, 11 significant QTL (P < 0.05), located in eight linkage groups, were detected on a genome-wise basis, in combined-sire (five QTL) or within-sire analyses (six QTL). BTA2 and BTA3 had QTL for both metabolic traits and meat quality traits. Six significant QTL for meat quality and metabolic traits were found at the proximal end of chromosome 2. BTA2 and BTA29 were the most common chromosomes harbouring QTL for meat quality traits; QTL for improved tenderness were associated with Limousin-derived and Jersey-derived alleles on these two chromosomes, respectively.

Mutations in the myostatin gene leading to hypermuscularity in mammals: indications for a similar mechanism in fish?

The transforming growth factor β (TGF-β) superfamily encodes secreted factors that are important in regulating embryonic development and tissue homeostatis in adults. Myostatin (MSTN, encoded by MSTN) or 'growth and differentiation factor 8', a member of this superfamily, is a negative regulator of skeletal muscle growth and is highly conserved among animal species. In 1997, a mutation associated with the so-called double-muscling phenotype in cattle was found in the MSTN gene. During the years following the discovery of the first MSTN mutation, other mutations were found in cattle and other mammalian species, and MSTN became one of the most thoroughly studied genes in animals. The aim of this review is mainly to describe the functional mutations located in the MSTN genes of several mammalian species, leading to double muscling in these animals. Furthermore, in light of the increasing importance of fish genetics, the possibility of functional mutations in piscine MSTN with a similar effect as in mammals, and a genetic model for MSTN research in fish, will also be discussed.

Discovery of novel genetic networks associated with 19 economically important traits in beef cattle

Quantitative or complex traits are determined by the combined effects of many loci, and are affected by genetic networks or molecular pathways. In the present study, we genotyped a total of 138 mutations, mainly single nucleotide polymorphisms derived from 71 functional genes on a Wagyu x Limousin reference population. Two hundred forty six F₂ animals were measured for 5 carcass, 6 eating quality and 8 fatty acid composition traits. A total of 2,280 single marker-trait association runs with 120 tagged mutations selected based on the HAPLOVIEW analysis revealed 144 significant associations (P < 0.05), but 50 of them were removed from the analysis due to the small number of animals (≤ 9) in one genotype group or absence of one genotype among three genotypes. The remaining 94 single-trait associations were then placed into three groups of quantitative trait modes (QTMs) with additive, dominant and overdominant effects. All significant markers and their QTMs associated with each of these 19 traits were involved in a linear regression model analysis, which confirmed single-gene associations for 4 traits, but revealed two-gene networks for 8 traits and three-gene networks for 5 traits. Such genetic networks involving both genotypes and QTMs resulted in high correlations between predicted and actual values of performance, thus providing evidence that the classical Mendelian principles of inheritance can be applied in understanding genetic complexity of complex phenotypes. Our present study also indicated that carcass, eating quality and fatty acid composition traits rarely share genetic networks. Therefore, marker-assisted selection for improvement of one category of these traits would not interfere with improvement of another.