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.

Integrated transcriptome analysis reveals roles of long non-coding RNAs (lncRNAs) in caprine skeletal muscle mass and meat quality

Long non-coding RNAs (lncRNAs) play important roles in the growth and development of skeletal muscle. However, there is limited information on goats. In this study, expression profiles of lncRNAs in Longissimus dorsi muscle from Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with divergent meat yield and meat quality were compared using RNA-sequencing. Based on our previous microRNA (miRNA) and mRNA profiles obtained from the same tissues, the target genes and binding miRNAs of differentially expressed lncRNAs were obtained. Subsequently, lncRNA-mRNA interaction networks and a ceRNA network of lncRNA-miRNA-mRNA were constructed. A total of 136 differentially expressed lncRNAs were identified between the two breeds. Fifteen cis target genes and 143 trans target genes were found for differentially expressed lncRNAs, and they were enriched in muscle contraction, muscle system process, muscle cell differentiation, and p53 signaling pathway. A total of 69 lncRNA-trans target gene pairs were constructed, with close relationship with muscle development, intramuscular fat deposition, and meat tenderness. A total of 16 lncRNA-miRNA-mRNA ceRNA pairs were identified, of which some reportedly associated with skeletal muscle development and fat deposition were found. The study will provide an improved understanding of the roles of lncRNAs in caprine meat yield and meat quality.

A conserved haplotype in Wagyu cattle contains RAB4A whose encoded protein regulates glucose trafficking in muscle and fat cells

The Wagyu breed of taurine cattle possess favourable genetics for intramuscular fat (IMF) but genomic loci associated with the trait remain under characterised. Here, we report the identification of a previously unidentified genomic region possessing a particular haplotype structure in Wagyu. Through deployment of a genome-wide haplotype detection analysis that captures regions conserved in a target population but not other populations we screened 100 individual Wagyu and contrasted them with 100 individuals from two independent comparison breeds, Charolais and Angus, using high-density SNPs. An extreme level of Wagyu conservation was assigned to a single genomic window (spanning genomic coordinates BTA28:41 088-300 265 bp). In fact, a five-SNP region spanning 27 096 bp is almost perfectly conserved among the 100 Wagyu individuals assayed and partially overlaps RAB4A. Focussing in, two consecutive SNPs (genomic coordinates 236 949 and 239 950) are apparently fixed within the Wagyu (BB and AA respectively), but at mixed frequencies in the other two breeds. These SNPs are located in the two introns straddling exon 7. In a separate analysis using the 1000 Bulls database, we found that, coincident with exon 7 of RAB4A first allele frequencies were highest in the high IMF Japanese Native (Wagyu) breeds (0.78) and lowest in the low IMF indicine breeds (Nelore and Brahman), with intermediate marbling breeds (Angus and Charolais) assigned intermediate rankings (0.42). RAB4A is known to encode a protein that regulates intracellular trafficking of the insulin-regulated glucose transporter GLUT4. RAB4A can be considered an attractive new positional candidate for IMF development.

Identification of Genomic Regions for Carcass Quality Traits within the American Simmental Association Carcass Merit Program

USDA quality and yield grade are primary driving forces for carcass value in the United States. Carcass improvements can be achieved by making selection decisions based on the results of genetic evaluations in the form of expected progeny differences (EPD), real-time ultrasound imaging, and physical evaluation of candidate breeding animals. In an effort to advance their ability to accurately predict the breeding value of potential sires for carcass traits, the American Simmental Association launched the Carcass Merit Program as a means to collect progeny sire group carcass information. All records were extracted from the American Simmental Association database. Progeny data were organized by sire family and progeny performance phenotypes were constructed. Sire genotypes were filtered, and a multi-locus mixed linear model was used to perform an association analysis on the genotype data, while correcting for cryptic relatedness and pedigree structure. Three chromosomes were found to have genome-wide significance and this conservative approach identified putative QTL in those regions. Three hundred ninety-three novel regions were identified across all traits, as well as 290 novel positional candidate genes. Correlations between carcass characteristics and maternal traits were less unfavorable than those previously reported.

Fine mapping of quantitative trait loci underlying sensory meat quality traits in three French beef cattle breeds

Improving the traits that underlie meat quality is a major challenge in the beef industry. The objective of this paper was to detect QTL linked to sensory meat quality traits in 3 French beef cattle breeds. We genotyped 1,059, 1,219, and 947 young bulls and their sires belonging to the Charolais, Limousin, and Blonde d'Aquitaine breeds, respectively, using the Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA). After estimating relevant genetic parameters using VCE software, we performed a linkage disequilibrium and linkage analysis on 4 meat traits: intramuscular fat content, muscle lightness, shear force, and tenderness score. Heritability coefficients largely ranged between 0.10 and 0.24; however, they reached a maximum of 0.44 and 0.50 for intramuscular fat content and tenderness score, respectively, in the Charolais breed. The 2 meat texture traits, shear force and tenderness score, were strongly genetically correlated (-0.91 in the Charolais and Limousin breed and -0.86 in the Blonde d'Aquitaine breed), indicating that they are 2 different measures of approximately the same trait. The genetic correlation between tenderness and intramuscular fat content differed across breeds. Using a significance threshold of 5 × 10(-4) for QTL detection, we found more than 200 significant positions across the 29 autosomal chromosomes for the 4 traits in the Charolais and Blonde d'Aquitaine breeds; in contrast, there were only 78 significant positions in the Limousin breed. Few QTL were common across breeds. We detected QTL for intramuscular fat content located near the myostatin gene in the Charolais and Blonde d'Aquitaine breeds. No mutation in this gene has been reported for the Blonde d'Aquitaine breed; therefore, it suggests that an unknown mutation could be segregating in this breed. We confirmed that, in certain breeds, markers in the calpastatin and calpain 1 gene regions affect tenderness. We also found new QTL as several QTL on chromosome 3 that are significantly associated with meat tenderness in the Blonde d'Aquitaine breed. Overall, these results greatly contribute to the goal of building a panel of markers that can be used to select animals of high meat quality.

Characterization of genes for beef marbling based on applying gene coexpression network

Marbling is an important trait in characterization beef quality and a major factor for determining the price of beef in the Korean beef market. In particular, marbling is a complex trait and needs a system-level approach for identifying candidate genes related to the trait. To find the candidate gene associated with marbling, we used a weighted gene coexpression network analysis from the expression value of bovine genes. Hub genes were identified; they were topologically centered with large degree and BC values in the global network. We performed gene expression analysis to detect candidate genes in M. longissimus with divergent marbling phenotype (marbling scores 2 to 7) using qRT-PCR. The results demonstrate that transmembrane protein 60 (TMEM60) and dihydropyrimidine dehydrogenase (DPYD) are associated with increasing marbling fat. We suggest that the network-based approach in livestock may be an important method for analyzing the complex effects of candidate genes associated with complex traits like marbling or tenderness.

Effects of new polymorphisms in the bovine myocyte enhancer factor 2D (MEF2D) gene on the expression rates of the longissimus dorsi muscle

Myocyte enhancer factor 2D (MEF2D), a product of the MEF2D gene, belongs to the myocyte enhancer factor 2 (MEF2) protein family which is involved in vertebrate skeletal muscle development and differentiation during myogenesis. The aim of the present study was to search for polymorphisms in the bovine MEF2D gene and to analyze their effect on MEF2D mRNA and on protein expression levels in the longissimus dorsi muscle of Polish Holstein–Friesian cattle. Overall, three novel variations, namely, insertion/deletion g.−818_−814AGCCG and g.−211C<A transversion in the promoter region as well as g.7C<T transition in the 5′untranslated region (5′UTR), were identified by DNA sequencing. A total, 375 unrelated bulls belonging to six different cattle breeds were genotyped, and three combined genotypes (Ins-C-C/Ins-C-C, Del-A-T/Del-A-T and Ins-C-C/Del-A-T) were determined. The frequency of the combined genotype Ins-C-C/Ins-C-C and Del-A-T/Del-A-T was varied between the breeds and the average frequency was 0.521 and 0.037, respectively. Expression analysis showed that the MEF2D variants were highly correlated with MEF2D mRNA and protein levels in the longissimus dorsi muscle of Polish Holstein–Friesian bulls carrying the three different combined genotypes. The highest MEF2D mRNA and protein levels were estimated in the muscle of bulls with the Ins-C-C/Ins-C-C homozygous genotype as compared to the Del-A-T/Del-A-T homozygotes (P < 0.01) and Ins-C-C/Del-A-T heterozygotes (P < 0.05). A preliminary association study showed no significant differences in the carcass quality traits between bulls with various MEF2D combined genotypes in the investigated population of Polish Holstein–Friesian cattle.

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.

Identification of candidate genes related to bovine marbling using protein-protein interaction networks

Complex traits are determined by the combined effects of many loci and are affected by gene networks or biological pathways. Systems biology approaches have an important role in the identification of candidate genes related to complex diseases or traits at the system level. The present study systemically analyzed genes associated with bovine marbling score and identified their relationships. The candidate nodes were obtained using MedScan text-mining tools and linked by protein-protein interaction (PPI) from the Human Protein Reference Database (HPRD). To determine key node of marbling, the degree and betweenness centrality (BC) were used. The hub nodes and biological pathways of our network are consistent with the previous reports about marbling traits, and also suggest unknown candidate genes associated with intramuscular fat. Five nodes were identified as hub genes, which was consistent with the network analysis using quantitative reverse-transcription PCR (qRT-PCR). Key nodes of the PPI network have positive roles (PPARγ, C/EBPα, and RUNX1T1) and negative roles (RXRA, CAMK2A) in the development of intramuscular fat by several adipogenesis-related pathways. This study provides genetic information for identifying candidate genes for the marbling trait in bovine.

The effect of measurement error of phenotypes on genome wide association studies

Background

There is an unspoken assumption that imprecision of measurement of phenotypes will not have large systematic effects on the location of significant associations in a genome wide association study (GWAS). In this report, the effects of two independent measurements of the same trait, subcutaneous fat thickness, were examined in GWAS of 940 individuals.

Results

The trait values obtained by two independent groups working to the same trait definition were correlated with r = 0.72. The allele effects obtained from the two analyses were only moderately correlated, with r = 0.53, and there was one significant (P < 0.0001) association in common to the two measurements. The correlation between allele effects was approximately equal to the square of the correlation between the trait measurements. An important quantitative trait locus (QTL) on BTA14 appeared to be shifted distally by 1 Mb along the chromosome. The divergence in GWAS was stronger with data coded into two discrete classes. Univariate trimming of the top and bottom 5% of data, a method used to control for erroneous trait values, decreased the similarity between the GWAS and increased the apparent shift of the QTL on BTA14. Stringent bivariate trimming of data, using only trait values that were similar to each other in the two data sets, substantially improved the correlation of trait values and allele effects in the GWAS, and showed evidence for two QTL on BTA14 separated by 1 Mb. Despite the reduction in sample size due to trimming, more SNP were significant. Using the mean of the two measurements of the trait was not as efficient as bivariate trimming.

Conclusions

It is recommended that trait values in GWAS experiments be examined for repeatability before the experiment is performed. For traits that do not have high repeatability (r < 0.95), two or more independent measurements of the same trait should be obtained for all samples, and individuals genotyped that have highly correlated trait measurements.