Expression patterns and genetic variation of the ovine skeletal muscle transcriptome of sheep from five Spanish meat breeds

Transcriptomic and metabolomic dissection of skeletal muscle of crossbred Chongming white goats with different meat production performance

BACKGROUND

The transcriptome and metabolome dissection of the skeletal muscle of high- and low- growing individuals from a crossbred population of the indigenous Chongming white goat and the Boer goat were performed to discover the potential functional differentially expressed genes (DEGs) and differential expression metabolites (DEMs).

RESULTS

A total of 2812 DEGs were detected in 6 groups at three time stages (3,6,12 Month) in skeletal muscle using the RNA-seq method. A DEGs set containing seven muscle function related genes (TNNT1, TNNC1, TNNI1, MYBPC2, MYL2, MHY7, and CSRP3) was discovered, and their expression tended to increase as goat muscle development progressed. Seven DEGs (TNNT1, FABP3, TPM3, DES, PPP1R27, RCAN1, LMOD2) in the skeletal muscle of goats in the fast-growing and slow-growing groups was verified their expression difference by reverse transcription-quantitative polymerase chain reaction. Further, through the Liquid chromatography-mass spectrometry (LC-MS) approach, a total of 183 DEMs in various groups of the muscle samples and these DEMs such as Queuine and Keto-PGF1α, which demonstrated different abundance between the goat fast-growing group and slow-growing group. Through weighted correlation network analysis (WGCNA), the study correlated the DEGs with the DEMs and identified 4 DEGs modules associated with 18 metabolites.

CONCLUSION

This study benefits to dissection candidate genes and regulatory networks related to goat meat production performance, and the joint analysis of transcriptomic and metabolomic data provided insights into the study of goat muscle development.

Genome-wide association study reveals candidate genes for traits related to meat quality in Colombian Creole hair sheep

Genome-wide association studies (GWAS) allow identifying genomic regions related to traits of economic importance in animals of zootechnical interest. The objective of this research was to conduct a genome-wide association study on meat quality traits using the Illumina OvineSNPs50 BeadChip array. The animals were sampled in the departments of Córdoba, Cesar, and Valle del Cauca. The genotypes obtained with the Illumina OvineSNP50 BeadChip microarray were analyzed SNP (single-nucleotide polymorphism) data to conduct a GWAS for pH and water-holding capacity (WHC) traits measured after 7 days of maturation, in the Longissimus dorsi (LD) muscle, in 167 Creole hair sheep of 12 months old belonging to Pelibuey (CHSP, n = 60), Ethiopian (CHSE, n = 44), and Sudan (CHSS, n = 63) breeds. The GWAS was done using a mixed linear model (MLMA) and based on the Ovis aries v3.1 genome. The CHSE showed the lowest meat juice release and, consequently, the highest water-holding capacity (WHC = 30.6 ± 0.1), suggesting that this breed has better performance in the meat industry compared with CHSS (WHC = 41.7 ± 0.1) and CHSP (WHC = 36.8 ± 0.1), since there is a relationship between WHC and juiciness. For the character pH, it was not possible to annotate genes related to meat quality, while, for the WHC, they have obtained 11 candidate genes associated (ELOVL2, ARAP2, LOC101102527, SHOC2, AIPL1, CSRNP3, IFRD, KDM8, NANS, DAPK1, IBN2, TPM2). Particularly, ELOVL2, ARAP2, IBN2, and TPM2 genes are involved in muscle contraction and fatty acid composition in sheep. In this study, we generated a baseline for GWAS related to meat quality traits in Colombian Creole hair sheep that can be used for future genomic selection plans.

Hepatic Transcriptome Analysis Reveals Genes, Polymorphisms, and Molecules Related to Lamb Tenderness

Tenderness is a key meat quality trait that determines the public acceptance of lamb consumption, so genetic improvement toward lamb with higher tenderness is pivotal for a sustainable sheep industry. However, unravelling the genomics controlling the tenderness is the first step. Therefore, this study aimed to identify the transcriptome signatures and polymorphisms related to divergent lamb tenderness using RNA deep sequencing. Since the molecules and enzymes that control muscle growth and tenderness are metabolized and synthesized in the liver, hepatic tissues of ten sheep with divergent phenotypes: five high- and five low-lamb tenderness samples were applied for deep sequencing. Sequence analysis identified the number of reads ranged from 21.37 to 25.37 million bases with a mean value of 22.90 million bases. In total, 328 genes are detected as differentially expressed (DEGs) including 110 and 218 genes that were up- and down-regulated, respectively. Pathway analysis showed steroid hormone biosynthesis as the dominant pathway behind the lamb tenderness. Gene expression analysis identified the top high (such as TP53INP1, CYP2E1, HSD17B13, ADH1C, and LPIN1) and low (such as ANGPTL2, IGFBP7, FABP5, OLFML3, and THOC5) expressed candidate genes. Polymorphism and association analysis revealed that mutation in OLFML3, ANGPTL2, and THOC5 genes could be potential candidate markers for tenderness in sheep. The genes and pathways identified in this study cause variation in tenderness, thus could be potential genetic markers to improve meat quality in sheep. However, further validation is needed to confirm the effect of these markers in different sheep populations so that these could be used in a selection program for lamb with high tenderness.

The trehalose glycolipid C18Brar promotes antibody and T-cell immune responses to Mannheimia haemolytica and Mycoplasma ovipneumoniae whole cell antigens in sheep

Bronchopneumonia is a common respiratory disease in livestock. Mannheimia haemolytica is considered the main causative pathogen leading to lung damage in sheep, with Mycoplasma ovipneumoniae and ParaInfluenza virus type 3, combined with adverse physical and physiological stress, being predisposing factors. A balance of humoral and cellular immunity is thought to be important for protection against developing respiratory disease. In the current study, we compared the ability of the trehalose glycolipid adjuvant C18Brar (C18-alkylated brartemicin analogue) and three commercially available adjuvant systems i.e., Quil-A, Emulsigen-D, and a combination of Quil-A and aluminium hydroxide gel, to stimulate antibody and cellular immune responses to antigens from inactivated whole cells of M. haemolytica and M. ovipneumoniae in sheep. C18Brar and Emulsigen-D induced the strongest antigen-specific antibody responses to both M. haemolytica and M. ovipneumoniae, while C18Brar and Quil-A promoted the strongest antigen-specific IL-17A responses. The expression of genes with known immune functions was determined in antigen-stimulated blood cultures using Nanostring nCounter technology. The expression levels of CD40, IL22, TGFB1, and IL2RA were upregulated in antigen-stimulated blood cultures from animals vaccinated with C18Brar, which is consistent with T-cell activation. Collectively, the results demonstrate that C18Brar can promote both antibody and cellular responses, notably Th17 immune responses in a ruminant species.

lncRNA-miRNA-mRNA ceRNA Network Involved in Sheep Prolificacy: An Integrated Approach

Understanding the molecular pattern of fertility is considered as an important step in breeding of different species, and despite the high importance of the fertility, little success has been achieved in dissecting the interactome basis of sheep fertility. However, the complex mechanisms associated with prolificacy in sheep have not been fully understood. Therefore, this study aimed to use competitive endogenous RNA (ceRNA) networks to evaluate this trait to better understand the molecular mechanisms responsible for fertility. A competitive endogenous RNA (ceRNA) network of the corpus luteum was constructed between Romanov and Baluchi sheep breeds with either good or poor genetic merit for prolificacy using whole-transcriptome analysis. First, the main list of lncRNAs, miRNAs, and mRNA related to the corpus luteum that alter with the breed were extracted, then miRNA−mRNA and lncRNA−mRNA interactions were predicted, and the ceRNA network was constructed by integrating these interactions with the other gene regulatory networks and the protein−protein interaction (PPI). A total of 264 mRNAs, 14 lncRNAs, and 34 miRNAs were identified by combining the GO and KEGG enrichment analyses. In total, 44, 7, 7, and 6 mRNAs, lncRNAs, miRNAs, and crucial modules, respectively, were disclosed through clustering for the corpus luteum ceRNA network. All these RNAs involved in biological processes, namely proteolysis, actin cytoskeleton organization, immune system process, cell adhesion, cell differentiation, and lipid metabolic process, have an overexpression pattern (Padj < 0.01). This study increases our understanding of the contribution of different breed transcriptomes to phenotypic fertility differences and constructed a ceRNA network in sheep (Ovis aries) to provide insights into further research on the molecular mechanism and identify new biomarkers for genetic improvement.

Circular RNA, microRNA and Protein Profiles of the Longissimus Dorsi of Germany ZIKA and Sichuan White Rabbits

Due to the dietetic properties and remarkable nutritive value of rabbit meat, its industry is increasing rapidly. However, the association between circular RNAs, microRNAs, and proteins and muscle fiber type, and meat quality of rabbit is still unknown. Here, using deep sequencing and iTRAQ proteomics technologies we first identified 3159 circRNAs, 356 miRNAs, and 755 proteins in the longissimus dorsi tissues from Sichuan white (SCWrabs) and Germany great line ZIKA rabbits (ZIKArabs). Next, we identified 267 circRNAs, 3 miRNAs, and 29 proteins differentially expressed in the muscle tissues of SCWrabs and ZIKArabs. Interaction network analysis revealed some key regulation relationships between noncoding RNAs and proteins that might be associated with the muscle fiber type and meat quality of rabbit. Further, miRNA isoforms and gene variants identified in SCWrabs and ZIKArabs revealed some pathways and biological processes related to the muscle development. This is the first study of noncoding RNA and protein profiles for the two rabbit breeds. It provides a valuable resource for future studies in rabbits and will improve our understanding of the molecular regulation mechanisms in the muscle development of livestock. More importantly, the output of our study will benefit the researchers and producers in the rabbit breeding program.

Novel Genes Associated with Dairy Traits in Sarda Sheep

The aim of the present research was to analyze the variability of 45 SNPs from different genes involved in metabolism and innate immunity to perform an association analysis with the milk yield, composition and milk coagulation traits. A population of 1112 Sarda breed sheep was sampled. Genotyping was generated by a TaqMan Open Array^TM. Thirty out of the 45 SNPs were polymorphic, and 12 displayed a minor allele frequency higher than 0.05. An association analysis showed that the variability at genes PRKAG3 and CD14 was significantly associated with the daily milk yield. The variability at PRKAG3 was also associated with the protein and casein content, somatic cell score and bacterial score. The variation at the PRKAA2 gene was associated with the milk lactose concentration. The SNPs at CD14 were also associated with the traditional milk coagulation properties, while the SNPs at GHR and GHRHR were associated with k_SR, a derived coagulation parameter related to the rate of syneresis. The information provided here is new and increases our knowledge of genotype-phenotype interactions in sheep. Our findings might be useful in appropriate breeding schemes to be set up for the Sarda sheep breed, but these should be confirmed by further studies, possibly performed on independent populations.

Key Genes Regulating Skeletal Muscle Development and Growth in Farm Animals

Simple Summary

Skeletal muscle mass is an important economic trait, and muscle development and growth is a crucial factor to supply enough meat for human consumption. Thus, understanding (candidate) genes regulating skeletal muscle development is crucial for understanding molecular genetic regulation of muscle growth and can be benefit the meat industry toward the goal of increasing meat yields. During the past years, significant progress has been made for understanding these mechanisms, and thus, we decided to write a comprehensive review covering regulators and (candidate) genes crucial for muscle development and growth in farm animals. Detection of these genes and factors increases our understanding of muscle growth and development and is a great help for breeders to satisfy demands for meat production on a global scale.

Abstract

Farm-animal species play crucial roles in satisfying demands for meat on a global scale, and they are genetically being developed to enhance the efficiency of meat production. In particular, one of the important breeders’ aims is to increase skeletal muscle growth in farm animals. The enhancement of muscle development and growth is crucial to meet consumers’ demands regarding meat quality. Fetal skeletal muscle development involves myogenesis (with myoblast proliferation, differentiation, and fusion), fibrogenesis, and adipogenesis. Typically, myogenesis is regulated by a convoluted network of intrinsic and extrinsic factors monitored by myogenic regulatory factor genes in two or three phases, as well as genes that code for kinases. Marker-assisted selection relies on candidate genes related positively or negatively to muscle development and can be a strong supplement to classical selection strategies in farm animals. This comprehensive review covers important (candidate) genes that regulate muscle development and growth in farm animals (cattle, sheep, chicken, and pig). The identification of these genes is an important step toward the goal of increasing meat yields and improves meat quality.

Expression patterns and genetic variation of the ovine skeletal muscle transcriptome of sheep from five Spanish meat breeds

The goal of the current study is to analyse the gene expression profile of the ovine skeletal muscle as well as to characterize the genetic variation of transcripts expressed in such tissue. This aim has been achieved by sequencing the longissimus dorsi transcriptomes of 50 sheep distributed in five pools representing the Canaria de Pelo, Roja Mallorquina, Gallega, Xisqueta and Ripollesa Spanish autochthonous breeds. Approximately, 363 million reads per pool have been produced and 71.9–82.9% have been successfully mapped to the ovine genome in a paired-end mode (2 × 75 bp). The 200 most expressed muscle transcripts (≈1% of the total transcript count) account for 51% (Canaria de Pelo) to 67% (Gallega) of the total ovine skeletal muscle mRNA expression. These highly expressed genes play key roles in pathways related with striated muscle contraction, gluconeogenesis, glycolysis, citric acid cycle and respiratory electron transport. RNA-Sequencing of muscle transcripts has also revealed that ~72% of the SNPs detected with this approach are shared by at least two pools, and 10% of them segregate in the five pools under analysis. Most of the substitutions detected by RNA-Seq are synonymous or missense and only a minority are predicted to have consequences on protein function.