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

Genome-wide Association Study for Carcass Primal Cut Yields Using Single-step Bayesian Approach in Hanwoo Cattle

The importance of meat and carcass quality is growing in beef cattle production to meet both producer and consumer demands. Primal cut yields, which reflect the body compositions of carcass, could determine the carcass grade and, consequently, command premium prices. Despite its importance, there have been few genome-wide association studies on these traits. This study aimed to identify genomic regions and putative candidate genes related to 10 primal cut traits, including tenderloin, sirloin, striploin, chuck, brisket, top round, bottom round, shank, flank, and rib in Hanwoo cattle using a single-step Bayesian regression (ssBR) approach. After genomic data quality control, 43,987 SNPs from 3,745 genotyped animals were available, of which 3,467 had phenotypic records for the analyzed traits. A total of 16 significant genomic regions (1-Mb window) were identified, of which five large-effect quantitative trait loci (QTLs) located on chromosomes 6 at 38–39 Mb, 11 at 21–22 Mb, 14 at 6–7 Mb and 26–27 Mb, and 19 at 26–27 Mb were associated with more than one trait, while the remaining 11 QTLs were trait-specific. These significant regions were harbored by 154 genes, among which TOX, FAM184B, SPP1, IBSP, PKD2, SDCBP, PIGY, LCORL, NCAPG, and ABCG2 were noteworthy. Enrichment analysis revealed biological processes and functional terms involved in growth and lipid metabolism, such as growth (GO:0040007), muscle structure development (GO:0061061), skeletal system development (GO:0001501), animal organ development (GO:0048513), lipid metabolic process (GO:0006629), response to lipid (GO:0033993), metabolic pathways (bta01100), focal adhesion (bta04510), ECM–receptor interaction (bta04512), fat digestion and absorption (bta04975), and Rap1 signaling pathway (bta04015) being the most significant for the carcass primal cut traits. Thus, identification of quantitative trait loci regions and plausible candidate genes will aid in a better understanding of the genetic and biological mechanisms regulating carcass primal cut yields.

Genes and proteins associated with ribeye area and meat tenderness in a commercial Nellore cattle population

Despite several studies on genetic markers and differently expressed genes related to ribeye area (REA) and tenderness traits in beef cattle, there is divergence in the results regarding the genes associated with these traits. Thirteen genes that had been associated or have biological functions that may influence such phenotypes were included in this study. A total of five genes for REA (IGF-1, IGF-2, MSTN, NEDD4, and UBE4A) and eight genes for meat tenderness (CAPN1, CAPN2, CAST, HSPB1, DNAJA1, FABP4, SCD, and PRKAG3) were selected from previously studies in beef cattle. Genes and its respective proteins expression were validated in a commercial population of Nellore cattle using quantitative real-time PCR (RT-qPCR) and advanced mass spectrometry (LC / MS-MS) techniques, respectively. MSTN gene was upregulated in animals with low REA. CAPN1, CAPN2, CAST, HSPB1, and DNAJA1 genes were upregulated in animals with tougher meat. The proteins translated by these genes were not differentially expressed. Our results could confirm the potential of some studied genes as biomarkers for carcass and meat quality in Nellore cattle.

Effect of dietary restriction and compensatory growth on performance, carcass characteristics, and metabolic hormone concentrations in Angus and Belgian Blue steers

Compensatory growth (CG) is the ability of an animal to undergo accelerated growth after a period of restricted feeding. However, there is a dearth of information in relation to the effect of genotype on CG response, thus the objective of this study was to evaluate CG response in two contrasting breed types, namely Aberdeen Angus (AN) and Belgian Blue (BB). Crossbred AN × Holstein-Friesian or BB × Holstein-Friesian steers were assigned to one of two treatment groups in a two (genotypes) × two (diets) factorial design. For 99 days, one group (11 AN and 12 BB) was offered a high energy control diet (H-H) whereas the second group (11 AN and 12 BB) was offered an energy restricted diet (L-H). At the end of the differential feeding period (99 days), both groups of animals were then offered a high energy control diet for a further 200 days. All animals were then slaughtered on day-299 of the study. During feed restriction, L-H had lower DM intake (DMI), had greater feed conversion ratio (FCR) and lower plasma concentrations of insulin, IGF-1, leptin, glucose, urea, betahydroxybutyrate and smaller M. longissimus thoracis or lumborum muscle and fat depths compared to H-H steers. During realimentation, there was no difference in DMI between diets; however, L-H had greater live weight gain compared to H-H steers. Overall, H-H consumed greater quantities on a DM basis, however, had a higher FCR compared to L-H steers. By the end of the realimentation period, there was no difference in plasma metabolite or hormone concentrations, linear body measurements, ultrasonically scanned fat depths, carcass conformation, dressing percentage or fat class between H-H and L-H steers. At slaughter, carcass weights were affected by diet with greater values for H-H compared to L-H steers. Genotype affected measures associated with body composition including pelvic width and both muscle and fat depths (P < 0.05). Overall, L-H had a CG (or recovery) index of 0.52 and did not make up for the loss of gains during the differential feeding period; however, M. longissimus thoracis et lumborum, a tissue of high economic value, recovered completely making it a target of interest for further investigation.

Transcriptomic diversity in longissimus thoracis muscles of Barbari and Changthangi goat breeds of India

The present study is an attempt to examine the differential expression of genes in longissimus thoracis muscles between meat and wool type Indian goat breeds. Barbari goat is considered the best meat breed while Changthangi is famous for its fine fibre quality. RNA sequencing data was generated from four biological replicates of longissimus thoracis muscles of Barbari and Changthangi goats. A clear demarcation could be observed between the breeds in terms of expression of genes associated with lipid metabolism (FASN, SCD, THRSP, DGAT2 and FABP3). Most significant genes with high connectivity identified by gene co-expression network analysis were associated with triacylglycerol biosynthesis pathway in Barbari goat. Highly interactive genes identified in Changthangi goat were mainly associated with muscle fibre type. This study provides an insight into the differential expression of genes in longissimus thoracis muscles between Barbari and Changthangi goats that are adapted to and reared in different agro-climatic regions.

Few Fixed Variants between Trophic Specialist Pupfish Species Reveal Candidate Cis-Regulatory Alleles Underlying Rapid Craniofacial Divergence

Investigating closely related species that rapidly evolved divergent feeding morphology is a powerful approach to identify genetic variation underlying variation in complex traits. This can also lead to the discovery of novel candidate genes influencing natural and clinical variation in human craniofacial phenotypes. We combined whole-genome resequencing of 258 individuals with 50 transcriptomes to identify candidate cis-acting genetic variation underlying rapidly evolving craniofacial phenotypes within an adaptive radiation of Cyprinodon pupfishes. This radiation consists of a dietary generalist species and two derived trophic niche specialists-a molluscivore and a scale-eating species. Despite extensive morphological divergence, these species only diverged 10 kya and produce fertile hybrids in the laboratory. Out of 9.3 million genome-wide SNPs and 80,012 structural variants, we found very few alleles fixed between species-only 157 SNPs and 87 deletions. Comparing gene expression across 38 purebred F1 offspring sampled at three early developmental stages, we identified 17 fixed variants within 10 kb of 12 genes that were highly differentially expressed between species. By measuring allele-specific expression in F1 hybrids from multiple crosses, we found that the majority of expression divergence between species was explained by trans-regulatory mechanisms. We also found strong evidence for two cis-regulatory alleles affecting expression divergence of two genes with putative effects on skeletal development (dync2li1 and pycr3). These results suggest that SNPs and structural variants contribute to the evolution of novel traits and highlight the utility of the San Salvador Island pupfish system as an evolutionary model for craniofacial development.

E3 Ubiquitin Ligase TRIP12: Regulation, Structure, and Physiopathological Functions

The Thyroid hormone Receptor Interacting Protein 12 (TRIP12) protein belongs to the 28-member Homologous to the E6-AP C-Terminus (HECT) E3 ubiquitin ligase family. First described as an interactor of the thyroid hormone receptor, TRIP12’s biological importance was revealed by the embryonic lethality of a murine model bearing an inactivating mutation in the TRIP12 gene. Further studies showed the participation of TRIP12 in the regulation of major biological processes such as cell cycle progression, DNA damage repair, chromatin remodeling, and cell differentiation by an ubiquitination-mediated degradation of key protein substrates. Moreover, alterations of TRIP12 expression have been reported in cancers that can serve as predictive markers of therapeutic response. The TRIP12 gene is also referenced as a causative gene associated to intellectual disorders such as Clark–Baraitser syndrome and is clearly implicated in Autism Spectrum Disorder. The aim of the review is to provide an exhaustive and integrated overview of the different aspects of TRIP12 ranging from its regulation, molecular functions and physio-pathological implications.

Transcriptomic profile of semitendinosus muscle of bulls of different breed and performance

The aim of the study was to compare the transcriptomic profiles of fully differentiated skeletal muscle derived from bulls belonging to different breeds of varying performance. Microarray analyses were performed to determine the differences in the expression profiles of genes between semitendinosus muscles of 15-month-old beef-breed bulls (Limousin—LIM and Hereford—HER) and dairy-breed bulls (Holstein Friesian—HF). These analyses allowed for the identification of those genes the expression of which is similar and characteristic of fully differentiated muscle in beef breeds, but differs in skeletal muscle of a typical dairy breed. The analysis revealed 463 transcripts showing similar expression in the semitendinosus muscle of beef breeds (LIM/HER), in comparison with the dairy breed (HF). Among the identified genes, 227 were upregulated and 236 were downregulated in beef breeds. The ontological analyses revealed that the largest group of genes similarly expressed in LIM and HER was involved in the processes of protein metabolism and development of muscle organ. In beef breeds, some genes involved in protein synthesis and proteolysis showed an upregulation, including ctsd, ctsf, fhl2, fhl3, fst, sirt1, and trim63, whereas some were downregulated, including bmpr1a, bmpr2, mstn, smad2, hspa8, gsk3β, and tgfβ2. The expression of the chosen genes was confirmed by RT-qPCR technique. Thus, it can be assumed that the identified genes involved in the regulation of growth and development of muscle tissue and the processes of protein metabolism in the examined cattle breeds may be responsible for the greater gain of muscle mass in beef-breed bulls.

Transcriptomic Profile of Primary Culture of Skeletal Muscle Cells Isolated from Semitendinosus Muscle of Beef and Dairy Bulls

The aim of the study was to identify differences in the transcriptomic profiles of primary muscle cell cultures derived from the semitendinosus muscle of bulls of beef breeds (Limousin (LIM) and Hereford (HER)) and a dairy breed (Holstein-Friesian (HF)) (n = 4 for each breed). Finding a common expression pattern for proliferating cells may point to such an early orientation of the cattle beef phenotype at the transcriptome level of unfused myogenic cells. To check this hypothesis, microarray analyses were performed. The analysis revealed 825 upregulated and 1300 downregulated transcripts similar in both beef breeds (LIM and HER) and significantly different when compared with the dairy breed (HF) used as a reference. Ontological analyses showed that the largest group of genes were involved in muscle organ development. Muscle cells of beef breeds showed higher expression of genes involved in myogenesis (including erbb-3, myf5, myog, des, igf-1, tgfb2) and those encoding proteins comprising the contractile apparatus (acta1, actc1, myh3, myh11, myl1, myl2, myl4, tpm1, tnnt2, tnnc1). The obtained results confirmed our hypothesis that the expression profile of several groups of genes is common in beef breeds at the level of proliferating satellite cells but differs from that observed in typical dairy breeds.

A Gene-Set Enrichment and Protein-Protein Interaction Network-Based GWAS with Regulatory SNPs Identifies Candidate Genes and Pathways Associated with Carcass Traits in Hanwoo Cattle

Non-synonymous SNPs and protein coding SNPs within the promoter region of genes (regulatory SNPs) might have a significant effect on carcass traits. Imputed sequence level data of 10,215 Hanwoo bulls, annotated and filtered to include only regulatory SNPs (450,062 SNPs), were used in a genome-wide association study (GWAS) to identify loci associated with backfat thickness (BFT), carcass weight (CWT), eye muscle area (EMA), and marbling score (MS). A total of 15, 176, and 1 SNPs were found to be significantly associated (p < 1.11 × 10-7) with BFT, CWT, and EMA, respectively. The significant loci were BTA4 (CWT), BTA6 (CWT), BTA14 (CWT and EMA), and BTA19 (BFT). BayesR estimated that 1.1%~1.9% of the SNPs contributed to more than 0.01% of the phenotypic variance. So, the GWAS was complemented by a gene-set enrichment (GSEA) and protein-protein interaction network (PPIN) analysis in identifying the pathways affecting carcass traits. At p < 0.005 (~2,261 SNPs), 25 GO and 18 KEGG categories, including calcium signaling, cell proliferation, and folate biosynthesis, were found to be enriched through GSEA. The PPIN analysis showed enrichment for 81 candidate genes involved in various pathways, including the PI3K-AKT, calcium, and FoxO signaling pathways. Our finding provides insight into the effects of regulatory SNPs on carcass traits.

Genomic analysis of bovine respiratory disease and lung consolidation in preweaned Holstein calves using clinical scoring and lung ultrasound

Bovine respiratory disease (BRD) is a leading cause of morbidity and mortality in dairy calves, with detrimental long-term effects that include stunted growth, increased age at first calving, and decreased milk production in first lactation. The objectives of this study were to establish a protocol for objective and efficient assessment of BRD phenotypes in preweaned dairy calves, develop a genomic reference population with well-defined clinical and subclinical phenotypes, identify chromosomal regions associated with BRD in a genome-wide association study, estimate genetic parameters of BRD, and predict genomic breeding values of dairy calves. A total of 1,107 Holstein calves from 6 dairy farms in southern Wisconsin were examined using clinical respiratory scoring and lung ultrasound at 3 and 6 wk of age. The clinical respiratory score was based on visual appraisal of eyes, nose, ears, cough, and temperature. Lung ultrasound scores were assigned based on the amount of consolidation present. Calves were genotyped with a commercially available SNP array and after quality control and imputation to higher density, 690,291 SNP markers and 1,014 individuals remained. Single-step genome-wide association study and single-step genomic best linear unbiased prediction were applied to binary 3- and 6-wk phenotypes considered as overall respiratory healthy versus affected (RESP3, RESP6) or as presence or absence of lung consolidation (CON3, CON6). Lung ultrasound combined with a clinical scoring system allowed for efficient and objective assessment for the prevalence of BRD. Proportions of variance attributed to 1-Mb non-overlapping windows suggested genomic regions that may contain putative candidate genes, most notably regions on Bos taurus autosomes 1, 6, 7, 10, 11, 12, 15, 17, 18, 27, and 28 that explained 0.70 to 1.45% of the genetic variance. Heritability estimates were higher at 3 wk (0.214 and 0.241 for CON3 and RESP3, respectively) than 6 wk (0.084 and 0.111 for CON6 and RESP6, respectively), and mean reliabilities of genomic estimated breeding vales for calves with genotypes and phenotypes ranged from 0.12 for CON6 to 0.30 for RESP3.

Muscle transcriptome signature and gene regulatory network analysis in two divergent lines of a hilly bovine species Mithun (Bos frontalis)

A massive bovine, Bos frontalis, also known as Mithun or Gayal, found at higher altitude is very promising meat and milk animal. For candidate gene and marker discovery, RNA-seq data was generated from longissimus dorsi muscle tissues with Illumina-HiSeq. Such markers can be used in future for genetic gain of traits like feed conversion efficiency (FCE) and average daily gain (ADG). Analysis revealed 297differentially expressed genes (DEGs) having 173 up and 124 down-regulated unigenes. Extensive conservation was found in genic region while comparing with Bos taurus. Analysis revealed 57 pathways having 112 enzymes, 72 transcriptional factors and cofactors, 212 miRNAs regulating 71 DEGs, 25,855 SSRs, mithun-specific 104,822 variants and 7288 indels, gene regulatory network (GRN) having 24 hub-genes and transcriptional factors regulating cell proliferation, immune tolerance and myogenesis. This is first report of muscle transcriptome depicting candidate genes with GRN controlling FCE and ADG. Reported putative molecular markers, candidate genes and hub proteins can be valuable genomic resources for association studies in genetic improvement programme.

Breed-dependent microRNA expression in the primary culture of skeletal muscle cells subjected to myogenic differentiation

Background

Skeletal muscle in livestock develops into meat, an important source of protein and other nutrients for human consumption. The muscle is largely composed of a fixed number of multinucleated myofibers determined during late gestation and remains constant postnatally. A population of postnatal muscle stem cells, called satellite cells, gives rise to myoblast cells that can fuse with the existing myofibers, thus increasing their size. This requires a delicate balance of transcription and growth factors and specific microRNA (miRNA) expressed by satellite cells and their supporting cells from the muscle stem cell niche. The role of transcription and growth factors in bovine myogenesis is well-characterized; however, very little is known about the miRNA activity during this process. We have hypothesized that the expression of miRNA can vary between primary cultures of skeletal muscle cells isolated from the semitendinosus muscles of different cattle breeds and subjected to myogenic differentiation.

Results

After a 6-day myogenic differentiation of cells isolated from the muscles of the examined cattle breeds, we found statistically significant differences in the number of myotubes between Hereford (HER)/Limousine (LIM) beef breeds and the Holstein-Friesian (HF) dairy breed (p ≤ 0.001). The microarray analysis revealed differences in the expression of 23 miRNA among the aforementioned primary cultures. On the basis of a functional analysis, we assigned 9 miRNA as molecules responsible for differentiation progression (miR-1, -128a, -133a, -133b, -139, -206, -222, -486, and -503). The target gene prediction and functional analysis revealed 59 miRNA-related genes belonging to the muscle organ development process.

Conclusion

The number of myotubes and the miRNA expression in the primary cultures of skeletal muscle cells derived from the semitendinosus muscles of the HER/LIM beef cattle breeds and the HF dairy breed vary when cells are subjected to myogenic differentiation. The net effect of the identified miRNA and their target gene action should be considered the result of the breed-dependent activity of satellite cells and muscle stem cell niche cells and their mutual interactions, which putatively can be engaged in the formation of a larger number of myotubes in beef cattle-related cells (HER/LIM) during in vitro myogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4492-5) contains supplementary material, which is available to authorized users.

Expression of genes related to the regulation of muscle protein turnover in Angus and Nellore bulls

We aimed to evaluate the expression of genes related to the regulation of muscle protein turnover in the longissimus dorsi (LD) muscle of Angus and Nellore bulls and to estimate the within-breed correlations of gene expression and performance traits. Thirteen genes related to the IGF-1 and myostatin pathways were studied. Thirteen animals, with an initial average BW of 381.2 ± 11.8 kg, from each breed were used in a completely randomized 2 × 2 factorial design (2 breeds and 2 feeding levels). The diet consisted of corn silage and a corn-soybean meal concentrate in a roughage-to-concentrate ratio of 30:70. Cattle were fed ad libitum (with 9 animals from each breed) or feed restricted (a 55% restriction of total DMI of ad libitum-fed animals, calculated as percentage of metabolic BW, with 4 animals of each breed). The experimental period lasted for 82 d and it was preceded by a 28-d adaptation period. The performance traits evaluated were slaughter body weight, total ADG (from d 1 to 82 of the trial), initial ADG (from d 1 to 41 of the trial), final ADG (from d 42 to 82 of the trial), total DMI (from d 1 to 82 of the trial), initial DMI (from d 1 to 41 of the trial), final DMI (from d 42 to 82 of the trial), HCW, LD weight (LDW), and rib eye area (REA). After slaughter, samples were taken from the LD muscle between the 12th and 13th ribs for gene expression analysis by quantitative reverse transcription PCR. There was no difference ( > 0.05) in the expression of any of the genes studied between ad libitum-fed Angus and ad libitum-fed Nellore, whereas feed restriction increased the expression of (; < 0.001), (; = 0.05), and (; = 0.04) and decreased the expression of ( < 0.01). The REA was negatively correlated to (; = 0.01), (; = 0.02), and ( = 0.05). The HCW was negatively correlated to ( = 0.01) and ( = 0.01) and tended to be negatively correlated to ( = 0.07), whereas the LDW tended to be negatively correlated to ( = 0.08). The genes , , and seem to be important for muscle growth and may be worthy of further investigation as future strategies for increasing muscle in livestock.

Transcriptomic profile adaptations following exposure of equine satellite cells to nutriactive phytochemical gamma-oryzanol

BACKGROUND

Adult skeletal muscle myogenesis depends on the activation of satellite cells that have the potential to differentiate into new fibers. Gamma-oryzanol (GO), a commercially available nutriactive phytochemical, has gained global interest on account of its muscle-building and regenerating effects. Here, we investigated GO for its potential influence on myogenesis, using equine satellite cell culture model, since the horse is a unique animal, bred and exercised for competitive sport. To our knowledge, this is the first report where the global gene expression in cultured equine satellite cells has been described.

METHODS

Equine satellite cells were isolated from semitendinosus muscle and cultured until the second day of differentiation. Differentiating cells were incubated with GO for the next 24 h. Subsequently, total RNA from GO-treated and control cells was isolated, amplified, labeled, and hybridized to two-color Horse Gene Expression Microarray slides. Quantitative PCR was used for the validation of microarray data.

RESULTS

Our results revealed 58 genes with changed expression in GO-treated vs. control cells. Analysis of expression changes suggests that various processes are reinforced by GO in differentiating equine satellite cells, including inhibition of myoblast differentiation, increased proliferation and differentiation, stress response, and increased myogenic lineage commitment.

CONCLUSIONS

The present study may confirm putative muscle-enhancing abilities of GO; however, the collective role of GO in skeletal myogenesis remains equivocal. The diversity of these changes is likely due to heterogenous growth rate of cells in primary culture. Genes identified in our study, modulated by the presence of GO, may become potential targets of future research investigating impact of this supplement in skeletal muscle on proteomic and biochemical level.

Genetic dissection of milk yield traits and mastitis resistance quantitative trait loci on chromosome 20 in dairy cattle

Intense selection to increase milk yield has had negative consequences for mastitis incidence in dairy cattle. Due to low heritability of mastitis resistance and an unfavorable genetic correlation with milk yield, a reduction in mastitis through traditional breeding has been difficult to achieve. Here, we examined quantitative trait loci (QTL) that segregate for clinical mastitis and milk yield on Bos taurus autosome 20 (BTA20) to determine whether both traits are affected by a single polymorphism (pleiotropy) or by multiple closely linked polymorphisms. In the latter but not the former situation, undesirable genetic correlation could potentially be broken by selecting animals that have favorable variants for both traits. First, we performed a within-breed association study using a haplotype-based method in Danish Holstein cattle (HOL). Next, we analyzed Nordic Red dairy cattle (RDC) and Danish Jersey cattle (JER) with the goal of determining whether these QTL identified in Holsteins were segregating across breeds. Genotypes for 12,566 animals (5,966 HOL, 5,458 RDC, and 1,142 JER) were determined by using the Illumina Bovine SNP50 BeadChip (50K; Illumina, San Diego, CA), which identifies 1,568 single nucleotide polymorphisms on BTA20. Data were combined, phased, and clustered into haplotype states, followed by within- and across-breed haplotype-based association analyses using a linear mixed model. Association signals for both clinical mastitis and milk yield peaked in the 26- to 40-Mb region on BTA20 in HOL. Single-variant association analyses were carried out in the QTL region using whole sequence level variants imputed from references of 2,036 HD genotypes (BovineHD BeadChip; Illumina) and 242 whole-genome sequences. The milk QTL were also segregating in RDC and JER on the BTA20-targeted region; however, an indication of differences in the causal factor(s) was observed across breeds. A previously reported F279Y mutation (rs385640152) within the growth hormone receptor gene showed strong association with milk, fat, and protein yields. In HOL, the highest peaks for milk yield and susceptibility to mastitis were separated by over 3.5 Mb (3.8 Mb by haplotype analysis, 3.6 Mb by single nucleotide polymorphism analysis), suggesting separate genetic variants for the traits. Further analysis yielded 2 candidate mutations for the mastitis QTL, at 33,642,072 bp (rs378947583) in an intronic region of the caspase recruitment domain protein 6 gene and 35,969,994 bp (rs133596506) in an intronic region of the leukemia-inhibitory factor receptor gene. These findings suggest that it may be possible to separate these beneficial and detrimental genetic factors through targeted selective breeding.

Deleted copy number variation of Hanwoo and Holstein using next generation sequencing at the population level

BACKGROUND

Copy number variation (CNV), a source of genetic diversity in mammals, has been shown to underlie biological functions related to production traits. Notwithstanding, there have been few studies conducted on CNVs using next generation sequencing at the population level.

RESULTS

Illumina NGS data was obtained for ten Holsteins, a dairy cattle, and 22 Hanwoo, a beef cattle. The sequence data for each of the 32 animals varied from 13.58-fold to almost 20-fold coverage. We detected a total of 6,811 deleted CNVs across the analyzed individuals (average length = 2732.2 bp) corresponding to 0.74% of the cattle genome (18.6 Mbp of variable sequence). By examining the overlap between CNV deletion regions and genes, we selected 30 genes with the highest deletion scores. These genes were found to be related to the nervous system, more specifically with nervous transmission, neuron motion, and neurogenesis. We regarded these genes as having been effected by the domestication process. Further analysis of the CNV genotyping information revealed 94 putative selected CNVs and 954 breed-specific CNVs.

CONCLUSIONS

This study provides useful information for assessing the impact of CNVs on cattle traits using NGS at the population level.

Transcriptional background of beef marbling - novel genes implicated in intramuscular fat deposition

The purpose of this study was to identify novel marbling-related genes by comparison of the global gene expression in semitendinosus muscle of 15-month-old Limousin (LIM), Holstein-Friesian (HF) and Hereford (HER) bulls. Muscle of LIM was lean with low intramuscular fat (IMF) content (0.53%) unlike the marbled muscles of HER and HF characterized by higher amounts of IMF (1.10 and 0.81%, respectively). The comparison of muscle transcriptional profile between marbled and lean beef revealed significant differences in expression of 144 genes, presumably involved in consecutive stages of adipose tissue development, such as preadipocyte proliferation and differentiation, adipocyte maturation, lipid filling and lipid metabolism leading to increased IMF deposition and marbling development. Correlation coefficients and regression analysis for nine of them (gadd45a, pias3, ccrn4l, diras3, pou5f1, hoxa9, atp2a2 and pim1) validated by real-time qPCR confirmed their moderate-high correlation with IMF% and explained up to 70.5% of the total variability in IMF deposition in the bulls.

Proteomic profiling of bovine M. longissimus lumborum from Crossbred Aberdeen Angus and Belgian Blue sired steers varying in genetic merit for carcass weight

Bovine skeletal muscle is a tissue of significant value to the beef industry and global economy. Proteomic analyses offer the opportunity to detect molecular mechanisms regulating muscle growth and intramuscular fat accumulation. The current study aimed to investigate differences in protein abundance in skeletal muscle tissue of cattle from two breeds of contrasting maturity (early vs. late maturing), adiposity, and muscle growth potential, namely, Belgian Blue (BB) × Holstein Friesian and Aberdeen Angus (AA) × Holstein Friesian. Twenty AA (n = 10) and BB (n = 10) sired steers, the progeny of sires of either high or low genetic merit, expressed as expected progeny difference for carcass weight (EPDcwt), and bred through AI, were evaluated as 4 genetic groups, BB-High, BB-Low, AA-High, and AA-Low (n = 5 per treatment). Chemical composition analysis of M. longissimus lumborum showed greater protein and moisture and decreased lipid concentrations for BB-sired compared with AA-sired steers. To investigate the effects of both sire breed and EPDcwt on M. longissimus lumborum, proteomic analysis was performed using 2-dimensional difference gel electrophoresis followed by mass spectrometry. Proteins were identified from their peptide sequences, using the National Center for Biotechnology Information (NCBI) and Swiss-prot databases. Metabolic enzymes involved in glycolysis (glycogen phosphorylase, phosphoglycerate mutase) and the citric acid cycle (aconitase 2, oxoglutarate dehydrogenase) were increased in AA- vs. BB-sired steers. Expression of proteins involved in cell structure, such as myosin light chain isoforms and troponins I and T, were also altered due to sire breed. Furthermore, heat shock protein β-1 and peroxiredoxin 6, involved in cell defense, had increased abundance in muscle of AA-sired relative to BB-sired steers. Protein abundance of glucose-6-phosphate isomerase, enolase-3, and pyruvate kinase was greater in AA-sired animals of High compared with Low EPDcwt. Changes in the expression of these proteins were supported by gene expression analysis using quantitative real-time PCR. This information will aid in our understanding of genetic influences controlling muscle growth and fat accumulation and could contribute to future breeding programs to increase lean tissue gain of beef cattle.

The transcriptomic signature of myostatin inhibitory influence on the differentiation of mouse C2C12 myoblasts

GDF8 (myostatin) is a unique cytokine strongly affecting the skeletal muscle phenotype in human and animals. The aim of the present study was to elucidate the molecular mechanism of myostatin influence on the differentiation of mouse C2C12 myoblasts, using the global-transcriptome analysis with the DNA microarray technique. Treatment with exogenous GDF8 strongly affected the growth and development of C2C12 mouse myoblasts. This was manifested by the inhibition of proliferation and differentiation as well as the impairment of cell fusion. DNA microarray analysis revealed 778 genes regulated by GDF8 in differentiating myoblasts (436 down-regulated and 235 up-regulated). Ontological analysis revealed their involvement in 17 types of biological processes, 10 types of molecular functions and 68 different signalling pathways. The effect of GDF8 was mainly mediated by the disruption of the cell cycle, calcium and insulin signalling pathways and expression of cytoskeletal and muscle specific proteins. The identified key-genes that could play a role as GDF8 targets in differentiating myoblasts are: Mef2, Hgf, Ilbl, Itgb1, Edn1, Ppargc1a.

Effect of sire breed and genetic merit for carcass weight on the transcriptional regulation of the somatotropic axis in longissimus dorsi of crossbred steers

The somatotropic axis plays an important role in postnatal growth, development, and differentiation of skeletal muscle. The aim of this study was to examine the effect of sire breed and sire EPD for carcass weight (EPD(cwt)) on the expression of components of the somatotropic axis in LM of beef cattle at slaughter. Crossbred Aberdeen Angus (AA; n = 17) and Belgian Blue (BB; n = 16) steers born to Holstein-Friesian dams and sired by bulls with either high (H) or low (L) EPD(cwt) were used in the study. Thus, there were 4 genetic groups [i.e., BBH (n = 8), BBL (n = 8), AAH (n = 8), and AAL (n = 9)]. Blood samples were collected via jugular venipuncture at regular intervals for analysis of plasma concentrations of IGF-1 and insulin. Total RNA was isolated from LM collected at slaughter, and the mRNA expression of IGF-1, IGF-2, their receptors (IGF-1R; IGF-2R), 6 IGFBP, acid labile subunit (ALS), and GH receptor (GHR) was measured by real-time reverse-transcription quantitative PCR. There was no effect of either sire breed or EPD(cwt) on concentrations of circulating IGF or insulin (P > 0.05). Gene expression of IGF-1R and IGFBP3 was upregulated in AA (P < 0.001) compared with BB, whereas IGF-1 was upregulated in H compared with L animals (P < 0.01). Correlation analysis indicated moderate positive associations between gene expression of IGFBP3 and IGF-1 (r = 0.54; P < 0.001) and IGF-1R (r = 0.48; P < 0.01). In addition, correlation analysis revealed that mRNA expression of IGFBP3 was moderately negatively associated with LM area per kilogram of carcass weight (r = -0.40; P < 0.05). Greater gene expression of IGF-1 and reduced transcript abundance of IGFBP3 in muscle may have a role in increased muscle growth potential in steers during the finishing period. These data will contribute to a better understanding of the molecular control of muscle growth at a tissue level in cattle.

Transcriptional pattern of TGF-β1 inhibitory effect on mouse C2C12 myoblasts differentiation

The aim of the present study was to define the effect of TGF-β1 on C2C12 myoblasts myogenesis. TGF-β1 together with its receptor is a negative auto-paracrine regulator of myogenesis, which influences the proliferation, differentiation, and functions of muscle cells. TGF-β1 exerts highly significant inhibitory effect on differentiation of C2C12 mouse myoblasts manifested by the impairment of cell fusion and very low expression of myosin heavy chain. The study of differentiating C2C12 mouse myoblasts treated with TGF-β1 revealed 502 genes (436 down-regulated and 66 up-regulated) with statistically different expression. TGF-β1-regulated genes were identified to be involved in 29 biological processes, 29 molecular functions groups and 59 pathways. The strongest inhibiting effect of TGF-β1 was observed in the cadherin and Wnt pathways. The key-genes that could play the role of TGF-β1 targets during myoblasts differentiation was identified such as: Max, Creb1, Ccna2, Bax, MdfI, Tef, Tubg1, Cxcl5, Rho, Calca and Lgals4.