Symmorphosis and livestock bioenergetics: production animal muscle has low mitochondrial volume fractions

Genetic and phenotypic associations of mitochondrial DNA copy number, SNP, and haplogroups with growth and carcass traits in beef cattle

Mitochondrial DNA copy number (mtDNA CN) is heritable and easily obtained from low-pass sequencing (LPS). This study investigated the genetic correlation of mtDNA CN with growth and carcass traits in a multi-breed and crossbred beef cattle population. Blood, leucocyte, and semen samples were obtained from 2,371 animals and subjected to LPS that resulted in nuclear DNA (nuDNA) and mtDNA sequence reads. Mitochondrial DNA CN was estimated as the ratio of mtDNA to nuDNA coverages. Variant calling was performed from mtDNA, and 11 single nucleotide polymorphisms (SNP) were identified in the population. Samples were classified in taurine haplogroups. Haplogroup and mtDNA type were further classified based on the 11 segregating SNP. Growth and carcass traits were available for between 7,249 and 60,989 individuals. Associations of mtDNA CN, mtDNA haplogroups, mtDNA types, and mtDNA SNP with growth and carcass traits were estimated with univariate animal models, and genetic correlations were estimated with a bivariate animal model based on pedigree. Mitochondrial DNA CN tended (P-value ≤0.08) to be associated with birth weight and weaning weight. There was no association (P-value >0.10) between mtDNA SNP, haplogroups, or types with growth and carcass traits. Genetic correlation estimates of mtDNA CN were -0.30 ± 0.16 with birth weight, -0.31 ± 0.16 with weaning weight, -0.15 ± 0.14 with post-weaning gain, -0.11 ± 0.19 with average daily dry-matter intake, -0.04 ± 0.22 with average daily gain, -0.29 ± 0.13 with mature cow weight, -0.11 ± 0.13 with slaughter weight, -0.14 ± 0.13 with carcass weight, -0.07 ± 0.14 with carcass backfat, 0.14 ± 0.14 with carcass marbling, and -0.06 ± 0.14 with ribeye area. In conclusion, mtDNA CN was negatively correlated with most traits investigated, and the genetic correlation was stronger with growth traits than with carcass traits.

Molecular regulation of high muscle mass in developing Blonde d'Aquitaine cattle foetuses

The Blonde d'Aquitaine (BA) is a French cattle breed with enhanced muscularity, partly attributable to a MSTN mutation. The BA m. Semitendinosus has a faster muscle fibre isoform phenotype comprising a higher proportion of fast type IIX fibres compared to age-matched Charolais (CH). To better understand the molecular network of modifications in BA compared to CH muscle, we assayed the transcriptomes of the m. Semitendinosus at 110, 180, 210 and 260 days postconception (dpc). We used a combination of differential expression (DE) and regulatory impact factors (RIF) to compare and contrast muscle gene expression between the breeds. Prominently developmentally regulated genes in both breeds reflected the replacement of embryonic myosin isoforms (MYL4, MYH3) with adult isoforms (MYH1) and the upregulation of mitochondrial metabolism (CKMT2, AGXT2L1) in preparation for birth. However, the transition to a fast, glycolytic muscle phenotype in the MSTN mutant BA is detectable through downregulation of various slow twitch subunits (TNNC1, MYH7, TPM3, CSRP3) beyond 210 dpc, and a small but consistent genome-wide reduction in mRNA encoding the mitoproteome. Across the breeds, NRIP2 is the regulatory gene possessing a network change most similar to that of MSTN.

Summary: We explored the molecular basis of high muscle mass in Blonde d’Aquitaine cattle, and found that mRNA encoding slow fibres and mitochondrial proteins were downregulated in this breed.

Progesterone signalling in broiler skeletal muscle is associated with divergent feed efficiency

Background

We contrast the pectoralis muscle transcriptomes of broilers selected from within a single genetic line expressing divergent feed efficiency (FE) in an effort to improve our understanding of the mechanistic basis of FE.

Results

Application of a virtual muscle model to gene expression data pointed to a coordinated reduction in slow twitch muscle isoforms of the contractile apparatus (MYH15, TPM3, MYOZ2, TNNI1, MYL2, MYOM3, CSRP3, TNNT2), consistent with diminishment in associated slow machinery (myoglobin and phospholamban) in the high FE animals. These data are in line with the repeated transition from red slow to white fast muscle fibres observed in agricultural species selected on mass and FE. Surprisingly, we found that the expression of 699 genes encoding the broiler mitoproteome is modestly–but significantly–biased towards the high FE group, suggesting a slightly elevated mitochondrial content. This is contrary to expectation based on the slow muscle isoform data and theoretical physiological capacity arguments. Reassuringly, the extreme 40 most DE genes can successfully cluster the 12 individuals into the appropriate FE treatment group. Functional groups contained in this DE gene list include metabolic proteins (including opposing patterns of CA3 and CA4), mitochondrial proteins (CKMT1A), oxidative status (SEPP1, HIG2A) and cholesterol homeostasis (APOA1, INSIG1). We applied a differential network method (Regulatory Impact Factors) whose aim is to use patterns of differential co-expression to detect regulatory molecules transcriptionally rewired between the groups. This analysis clearly points to alterations in progesterone signalling (via the receptor PGR) as the major driver. We show the progesterone receptor localises to the mitochondria in a quail muscle cell line.

Conclusions

Progesterone is sometimes used in the cattle industry in exogenous hormone mixes that lead to a ~20% increase in FE. Because the progesterone receptor can localise to avian mitochondria, our data continue to point to muscle mitochondrial metabolism as an important component of the phenotypic expression of variation in broiler FE.

Electronic supplementary material

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

Chicken muscle mitochondrial content appears co-ordinately regulated and is associated with performance phenotypes

Mitochondrial content is a fundamental cellular bioenergetic phenotype. Previous work has hypothesised possible links between variation in muscle mitochondrial content and animal performance. However, no population screens have been performed in any production species. Here, we have designed a high throughput molecular approach to estimate mitochondrial content in commercial broilers. Technical validity was established using several approaches, including its performance in monoclonal DF-1 cells, cross-tissue comparisons in tissues with differing metabolic demands (white fat<breast muscle<drumstick muscle<heart muscle) and, as a negative control, a near absence of mtDNA amplification from whole blood. We screened breast muscle and thigh muscle in 80 birds individually phenotyped for 11 growth and development traits. Substantial individual variation (fivefold) was discovered in both breast and thigh muscle mitochondrial content. Interestingly, across birds we detected a very strong positive relationship between breast and thigh content (correlation coefficient 0.61; P<0.0001), consistent with coordinate regulatory control across the musculature. Further, breast muscle mitochondrial content is negatively correlated with breast muscle yield (−0.27; P=0.037), abdominal fat content (−0.31; P=0.017) and carcass yield (−0.26; P=0.045). Therefore, low breast muscle mitochondrial content is associated with more muscular birds possessing higher abdominal fat, the latter being in line with biomedical models of obesity. Finally, thigh mitochondrial content is negatively correlated with the bow out leg defect (−0.30; P=0.011). Overall, our data point to mitochondrial content as a promising consideration in predictive modelling of production traits.

Summary: We have detected variation in muscle mitochondrial content across chickens. This variation relates to phenotypes and also indicates coordinate regulation of mitochondrial biogenesis across the musculature.

RNA sequencing for global gene expression associated with muscle growth in a single male modern broiler line compared to a foundational Barred Plymouth Rock chicken line

Background

Modern broiler chickens exhibit very rapid growth and high feed efficiency compared to unselected chicken breeds. The improved production efficiency in modern broiler chickens was achieved by the intensive genetic selection for meat production. This study was designed to investigate the genetic alterations accumulated in modern broiler breeder lines during selective breeding conducted over several decades.

Methods

To identify genes important in determining muscle growth and feed efficiency in broilers, RNA sequencing (RNAseq) was conducted with breast muscle in modern pedigree male (PeM) broilers (n = 6 per group), and with an unselected foundation broiler line (Barred Plymouth Rock; BPR). The RNAseq analysis was carried out using Ilumina Hiseq (2 x 100 bp paired end read) and raw reads were assembled with the galgal4 reference chicken genome. With normalized RPM values, genes showing >10 average read counts were chosen and genes showing <0.05 p-value and >1.3 fold change were considered as differentially expressed (DE) between PeM and BPR. DE genes were subjected to Ingenuity Pathway Analysis (IPA) for bioinformatic functional interpretation.

Results

The results indicate that 2,464 DE genes were identified in the comparison between PeM and BPR. Interestingly, the expression of genes encoding mitochondrial proteins in chicken are significantly biased towards the BPR group, suggesting a lowered mitochondrial content in PeM chicken muscles compared to BPR chicken. This result is inconsistent with more slow muscle fibers bearing a lower mitochondrial content in the PeM. The molecular, cellular and physiological functions of DE genes in the comparison between PeM and BPR include organismal injury, carbohydrate metabolism, cell growth/proliferation, and skeletal muscle system development, indicating that cellular mechanisms in modern broiler lines are tightly associated with rapid growth and differential muscle fiber contents compared to the unselected BPR line. Particularly, PDGF (platelet derived growth factor) signaling and NFE2L2 (nuclear factor, erythroid 2-like 2; also known as NRF2) mediated oxidative stress response pathways appear to be activated in modern broiler compared to the foundational BPR line. Upstream and network analyses revealed that the MSTN (myostatin) –FST (follistatin) interactions and inhibition of AR (androgen receptor) were predicted to be effective regulatory factors for DE genes in modern broiler line. PRKAG3 (protein kinase, AMP-activated, gamma 3 non-catalytic subunit) and LIPE (lipase E) are predicted as core regulatory factors for myogenic development, nutrient and lipid metabolism.

Conclusion

The highly upregulated genes in PeM may represent phenotypes of subclinical myopathy commonly observed in the commercial broiler breast tissue, that can lead to muscle hardening, named as woody breast. By investigating global gene expression in a highly selected pedigree broiler line and a foundational breed (Barred Plymouth Rock), the results provide insight into cellular mechanisms that regulate muscle growth, fiber composition and feed efficiency.

Electronic supplementary material

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

Mitochondrial metabolism: a driver of energy utilisation and product quality?

High feed efficiency is a very desirable production trait as it positively influences resource utilisation, profitability and environmental considerations, albeit at the possible expense of product quality. The modern broiler is arguably the most illustrative model species as it has been transformed over the past half century into an elite feed converter. Some producers are currently reporting that 42-day-old birds gain 1 kg of wet weight for every 1.35 kg of dry weight consumed. Its large breast muscle is exclusively composed of large, low mitochondrial-content Type IIB fibres, which may contribute to low maintenance costs and high efficiency. In an effort to gain a better understanding of individual variation in chicken feed efficiency, our group has been exploring the biology of the mitochondrion at multiple levels of organisation. The mitochondrion is the organelle where much biochemical energy transformation occurs in the cell. Using Cobb-Vantress industrial birds as our primary experimental resource, we have explored the tissue content, structure and function of the mitochondrion and its relationship to growth, development, efficiency and genetic background. While much remains to be understood, recent highlights include (1) variation in muscle mitochondrial content that is associated with performance phenotypes, (2) altered muscle mitochondrial gene and protein expression in birds differing in feed efficiency, (3) variation in isolated mitochondrial function in birds differing in feed efficiency and (4) evidence for an unexpected role for the mitochondrially localised progesterone receptor in altering bird muscle metabolism. Mitochondrial function is largely conserved across the vertebrates, so the same metabolic principles appear to apply to the major production species, whether monogastric or ruminant. A speculative role for the mitochondria in aspects of meat quality and in influencing postmortem anaerobic metabolism will conclude the manuscript.

RNF14 is a regulator of mitochondrial and immune function in muscle

BACKGROUND

Muscle development and remodelling, mitochondrial physiology and inflammation are thought to be inter-related and to have implications for metabolism in both health and disease. However, our understanding of their molecular control is incomplete.

RESULTS

In this study we have confirmed that the ring finger 14 protein (RNF14), a poorly understood transcriptional regulator, influences the expression of both mitochondrial and immune-related genes. The prediction was based on a combination of network connectivity and differential connectivity in cattle (a non-model organism) and mice data sets, with a focus on skeletal muscle. They assigned similar probability to mammalian RNF14 playing a regulatory role in mitochondrial and immune gene expression. To try and resolve this apparent ambiguity we performed a genome-wide microarray expression analysis on mouse C2C12 myoblasts transiently transfected with two Rnf14 transcript variants that encode 2 naturally occurring but different RNF14 protein isoforms. The effect of both constructs was significantly different to the control samples (untransfected cells and cells transfected with an empty vector). Cluster analyses revealed that transfection with the two Rnf14 constructs yielded discrete expression signatures from each other, but in both cases a substantial set of genes annotated as encoding proteins related to immune function were perturbed. These included cytokines and interferon regulatory factors. Additionally, transfection of the longer transcript variant 1 coordinately increased the expression of 12 (of the total 13) mitochondrial proteins encoded by the mitochondrial genome, 3 of which were significant in isolated pair-wise comparisons (Mt-coxII, Mt-nd2 and mt-nd4l). This apparent additional mitochondrial function may be attributable to the RWD protein domain that is present only in the longer RNF14 isoform.

CONCLUSIONS

RNF14 influences the expression of both mitochondrial and immune related genes in a skeletal muscle context, and has likely implications for the inter-relationship between bioenergetic status and inflammation.

Inferring the in vivo cellular program of developing bovine skeletal muscle from expression data

We outline an in vivo cellular program of bovine longissimus muscle development inferred from expression data from 60 days post conception to 3months postnatal. Analytic challenges included changes in cellular composition, ambiguous 'diagnostic' markers of cell type and contrasts between cattle human and mouse myogenesis. Nevertheless, the expression profiles of the myosin isoforms support slow and fast muscle fibres emanating from primary and secondary myogenesis respectively, while expression of the prenatal myosin subunits is down regulated prior to birth. Of the canonical pro-myogenic transcription factors (TF), MYF6 and MYF5 are negatively co-expressed, with MYF6 displaying higher expression in the post-natal samples and MYF5, MYOG, HES6 and PAX7 displaying higher expression in early development. A set of TFs (SIX1, EYA2 and DACH2) considered important in undifferentiated murine cells were equally abundant in differentiated bovine cells. An examination of mammalian regulators of fibre composition, muscle mass and muscle metabolism, underscored the roles of PPARGC1A, TGFβ signalling and the NHR4 Nuclear Hormone Receptors on bovine muscle development. Enriched among the most variably expressed genes from the entire data set were molecules regulating mitochondrial metabolism of carbohydrate (PDK4), fat (UCP3), protein (AGXT2L1) and high energy phosphate (CKMT2). The dramatic increase in the expression of these transcripts, which may enable the peri-natal transition to metabolic independence critical for new-born herbivores, provides surprising evidence for substantial developmental remodelling of muscle mitochondria and reflects changes in nutrient availability. Overall, despite differences in size, metabolism and physiology, the muscle structural subunit expression program appears very similar in ruminants, rodents and humans.

mRNA expression of genes regulating oxidative phosphorylation in the muscle of beef cattle divergently ranked on residual feed intake

Our objective was to evaluate the effects of phenotypic ranking on residual feed intake (RFI) on the transcription of genes 1) involved in the respiratory chain complex and 2) coding for transcriptional factors regulating mitochondrial biogenesis, across two contrasting diet types. Beef heifers (n = 86) fed a diet comprising 70:30 concentrate-corn silage [low forage (LF)] over a 82-day period were ranked on RFI. The 10 highest (feed inefficient, high-RFI) and 10 lowest (feed efficient, low-RFI) ranking animals were selected for the current study. Biopsies of the M. longissimus dorsi were harvested following initial selection (LF diet) and again following a 6 wk period while the animals were offered a high-forage (HF) grass silage-only diet. Real-time PCR was used to quantify mRNA transcripts of 17 genes associated with cellular energetic efficiency. The mRNA expression of UCP3 tended to be upregulated (2.2-fold, P = 0.06) for the high-RFI compared with the low-RFI animals. mRNA transcripts coding for the transcription factor PGC-1α was 1.7-fold higher (P = 0.01) in low compared with high-RFI animals. A phenotype × diet interaction was evident for the abundance of ANT1 mRNA transcript, with greater (P = 0.04) expression levels detected in the low-RFI phenotype during the HF period, but no difference (P = 0.50) between phenotypes during the LF period. A phenotype × diet interaction was also evident for COX II with greater expression levels detected (P = 0.04) in the low compared with the high RFI phenotype while on LF but not the HF diet (P = 0.22). These data suggest an association between cellular energetic efficiency and RFI in cattle.