Porcine muscle sensory attributes associate with major changes in gene networks involving CAPZB, ANKRD1, and CTBP2

Genome-wide detection of runs of homozygosity in Ding'an pigs revealed candidate genes relating to meat quality traits

BACKGROUND

Ding'an (DA) pig, a native Chinese breed, is renowned for its excellent meat quality, disease resistance, high reproductive performance, and adaptability. Its meat quality traits hold significant economic value. However, its conservation population has been declining due to the impact of commercialized breeds and African swine fever, which is not conducive to its development and utilization.

RESULTS

This study utilized whole-genome resequencing data from 15 DA pigs to reveal their genetic characteristics and current resource status. We analyzed the length, number, and distribution patterns of Runs of Homozygosity (ROH) in DA pigs, as well as high-frequency ROH regions. The results identified 23,208,098 single nucleotide polymorphisms (SNPs), 4,497,242 insertion and deletion (InDels), 13,622 copy number variation (CNVs), and 399,934 structure variation (SVs). Further analysis revealed relatively high genetic diversity and low inbreeding levels in DA pigs. Through functional gene enrichment analysis of high-frequency ROH regions, we identified multiple candidate genes associated with specific traits in DA pigs, including meat quality (ANKRD1, CPNE5, MYOM1), fat deposition (OBSCN, MAPK4, PNPLA1, PACSIN1, GRM4), and skeletal muscle development (LRPPRC, WNT9A).

CONCLUSIONS

This study conducted whole-genome sequencing and ROH analysis on DA pigs, revealing high genetic diversity and low inbreeding levels within the population. Through functional gene enrichment analysis of high-frequency ROH regions, we identified multiple candidate genes associated with meat quality, fat deposition, and skeletal muscle development. These findings not only enhance our understanding of the genetic mechanisms underlying the unique traits of DA pigs but also provide valuable insights for practical applications. Specifically, the identified candidate genes and genomic regions can guide conservation efforts to maintain genetic diversity and mitigate inbreeding risks. Meanwhile, these genetic insights can be integrated into breeding programs to improve meat quality and other economically important traits, thereby supporting the sustainable development and utilization of DA pigs.

Genome-Wide Association Study of Conformation Traits in Brazilian Holstein Cattle

The linear conformation of animals exerts an influence on health, reproduction, production, and welfare, in addition to longevity, which directly affects the profitability of milk-producing farms. The objectives of this study were (1) to perform genome-wide association studies (GWASs) of conformation traits, namely the Rump, Feet and Legs, Mammary System, Dairy Strength, and Final Classification traits, and (2) to identify genes and related pathways involved in physiological processes associated with conformation traits in Brazilian Holstein cattle. Phenotypic and genotypic data from 2339 Holstein animals distributed across the states of Rio Grande do Sul, Paraná, São Paulo, and Minas Gerais were used. The genotypic data were obtained with a 100 K SNP marker panel. The single-step genome-wide association study (ssGWAS) method was employed in the analyses. Genes close to a significant SNP were identified in an interval of 100 kb up- and downstream using the Ensembl database available in the BioMart tool. The DAVID database was used to identify the main metabolic pathways and the STRING program was employed to create the gene regulatory network. In total, 36 significant SNPs were found on 15 chromosomes; 27 of these SNPs were linked to genes that may influence the traits studied. Fourteen genes most closely related to the studied traits were identified, as well as four genes that showed interactions in important metabolic pathways such as myogenesis, adipogenesis, and angiogenesis. Among the total genes, four were associated with myogenesis (TMOD2, TMOD3, CCND2, and CTBP2), three with angiogenesis (FGF23, FGF1, and SCG3), and four with adipogenesis and body size and development (C5H12orf4, CCND2, EMILIN1, and FGF6). These results contribute to a better understanding of the biological mechanisms underlying phenotypic variability in conformation traits in Brazilian Holstein cattle.

A Whole-Genome Scan Revealed Genomic Features and Selection Footprints of Mengshan Cattle

(1) Background: Mengshan cattle from the Yimeng mountainous region in China stand out as a unique genetic resource, known for their adaptive traits and environmental resilience. However, these cattle are currently endangered and comprehensive genomic characterization remains largely unexplored. This study aims to address this gap by investigating the genomic features and selection signals in Mengshan cattle. (2) Methods: Utilizing whole-genome resequencing data from 122 cattle, including 37 newly sequenced Mengshan cattle, we investigated population structure, genetic diversity, and selection signals. (3) Results: Our analyses revealed that current Mengshan cattle primarily exhibit European taurine cattle ancestry, with distinct genetic characteristics indicative of adaptive traits. We identified candidate genes associated with immune response, growth traits, meat quality, and neurodevelopment, shedding light on the genomic features underlying the unique attributes of Mengshan cattle. Enrichment analysis highlighted pathways related to insulin secretion, calcium signaling, and dopamine synapse, further elucidating the genetic basis of their phenotypic traits. (4) Conclusions: Our results provide valuable insights for further research and conservation efforts aimed at preserving this endangered genetic resource. This study enhances the understanding of population genetics and underscores the importance of genomic research in informing genetic resources and conservation initiatives for indigenous cattle breeds.

A network-based approach to understanding gene-biological processes affecting economically important traits of Nelore cattle

This study aimed to build gene-biological process networks with differentially expressed genes associated with economically important traits of Nelore cattle from 17 previous studies. The genes were clustered into three groups by evaluated traits: group 1, production traits; group 2, carcass traits; and group 3, meat quality traits. For each group, a gene-biological process network analysis was performed with the differentially expressed genes in common. For production traits, 37 genes were found in common, of which 13 genes were enriched for six Gene Ontology (GO) terms; these terms were not functionally grouped. However, the enriched GO terms were related to homeostasis, the development of muscles and the immune system. For carcass traits, four genes were found in common. Thus, it was not possible to functionally group these genes into a network. For meat quality traits, the analysis revealed 222 genes in common. CSRP3 was the only gene differentially expressed in all three groups. Non-redundant biological terms for clusters of genes were functionally grouped networks, reflecting the cross-talk between all biological processes and genes involved. Many biological processes and pathways related to muscles, the immune system and lipid metabolism were enriched, such as striated muscle cell development and triglyceride metabolic processes. This study provides insights into the genetic mechanisms of production, carcass and meat quality traits of Nelore cattle. This information is fundamental for a better understanding of the complex traits and could help in planning strategies for the production and selection systems of Nelore cattle.

Genomic Regions Associated with Wool, Growth and Reproduction Traits in Uruguayan Merino Sheep

The aim of this study was to identify genomic regions and genes associated with the fiber diameter (FD), clean fleece weight (CFW), live weight (LW), body condition score (BCS), pregnancy rate (PR) and lambing potential (LP) of Uruguayan Merino sheep. Phenotypic records of approximately 2000 mixed-age ewes were obtained from a Merino nucleus flock. Genome-wide association studies were performed utilizing single-step Bayesian analysis. For wool traits, a total of 35 genomic windows surpassed the significance threshold (PVE ≥ 0.25%). The proportion of the total additive genetic variance explained by those windows was 4.85 and 9.06% for FD and CFW, respectively. There were 42 windows significantly associated with LWM, which collectively explained 43.2% of the additive genetic variance. For BCS, 22 relevant windows accounted for more than 40% of the additive genetic variance, whereas for the reproduction traits, 53 genomic windows (24 and 29 for PR and LP, respectively) reached the suggestive threshold of 0.25% of the PVE. Within the top 10 windows for each trait, we identified several genes showing potential associations with the wool (e.g., IGF-1, TGFB2R, PRKCA), live weight (e.g., CAST, LAP3, MED28, HERC6), body condition score (e.g., CDH10, TMC2, SIRPA, CPXM1) or reproduction traits (e.g., ADCY1, LEPR, GHR, LPAR2) of the mixed-age ewes.

Consequences of a low protein diet on the liver and longissimus dorsi transcriptome of Duroc × Iberian crossbred pigs

Low protein diets supplied during the growing period of pigs can diminish their growth rate and increase the intramuscular fat (IMF) content which affects the sensorial and technological characteristics of the products. In the present study, the effects of a low protein diet supplied during the growing period of Duroc × Iberian crossbred pigs on several phenotypic traits and on liver and longissimus dorsi transcriptome were analysed at the beginning (EARLY) and at the end (LATE) of the growing period. Two experimental groups of 10 crossbred pigs each were fed two isocaloric diets with different protein content: control diet (C) with 16.5% protein and 0.8% lysine and low protein diet (LP) with 11% CP and 0.6% lysine. Animals fed LP diet have a slower growth than those fed C diet, but no effect of LP diet was observed on the IMF content. The transcriptomes of liver and longissimus dorsi were characterised and quantified through RNA-sequencing (RNA-seq). In liver, 134 and 480 differentially expressed annotated genes and new isoforms (DEGs) were detected between C and LP diets for EARLY and LATE animals, respectively. In muscle, 128 and 68 DEGs were detected at EARLY and LATE time-points. Functional interpretation revealed that LP diet may inhibit immune system molecules and processes in both tissues at EARLY stage. In liver, the DEGs mainly affect lipid and cholesterol metabolic processes, while in muscle, the expression changes would be involved in growth, development and meat quality. In conclusion, a low protein diet supplied during the growing period seems to slow down the growth of Duroc × Iberian crossbred pigs, but it also seems to affect multiple biological processes that could compromise the immune system of Duroc × Iberian crossbred pigs. Therefore, these results question the adequacy of this type of regime in Duroc × Iberian pigs that must be studied in greater depth before being implemented.

Differences in the Loin Tenderness of Iberian Pigs Explained through Dissimilarities in Their Transcriptome Expression Profile

Simple Summary

The Iberian pig is the most representative autochthonous breed of the Mediterranean region with unique genetic and phenotypic characteristics. The breed has been successfully preserved by its high-quality meat and high-priced products. Tenderness is one of the most relevant meat quality traits, and meat tenderization is influenced by genetic and environmental effects such as pre-slaughter handling and post-mortem conditions. Tenderness could be included in Iberian pig breeding programs, mainly focused on the improvement of premium-cuts percentage, in order to avoid the meat quality decline. A better biological understanding of this trait is needed. In the current study, we analyze the transcriptome of pigs divergent for Warner–Bratzler shear force through RNA-seq technique for the identification, characterization and quantification of candidate genes involved in biological pathways, networks and functions affecting meat tenderness.

Abstract

Tenderness is one of the most important meat quality traits and it can be measured through shear force with the Warner–Bratzler test. In the current study, we use the RNA-seq technique to analyze the transcriptome of Longissimus dorsi (LD) muscle in two groups of Iberian pigs (Tough and Tender) divergent for shear force breeding values. We identified 200 annotated differentially expressed genes (DEGs) and 245 newly predicted isoforms. The RNAseq expression results of 10 genes were validated with quantitative PCR (qPCR). Functional analyses showed an enrichment of DE genes in biological processes related to proteolysis (CTSC, RHOD, MYH8, ACTC1, GADD45B, CASQ2, CHRNA9 and ANKRD1), skeletal muscle tissue development (ANKRD1, DMD, FOS and MSTN), lipid metabolism (FABP3 and PPARGC1A) and collagen metabolism (COL14A1). The upstream analysis revealed a total of 11 transcription regulatory factors that could regulate the expression of some DEGs. Among them, IGF1, VGLL3 and PPARG can be highlighted since they regulate the expression of genes involved in biological pathways that could affect tenderness. The experiment revealed a set of candidate genes and regulatory factors suggestive to search polymorphisms that could be incorporated in a breeding program for improving meat tenderness.

Cloning and expression profiling of muscle regulator ANKRD2 in domestic chicken Gallus gallus

Striated muscle signaling protein and transcriptional regulator ANKRD2 participates in myogenesis, myogenic differentiation, muscle adaptation and stress response. It is preferentially expressed in slow, oxidative fibers of mammalian skeletal muscle. In this study, we report on characterization of chicken ANKRD2. The chicken ANKRD2 coding region contains 1002 bp and encodes a 334-amino acid protein which shares approximately 58% identity with human and mouse orthologs, mostly in the conserved region of ankyrin repeats. Comprehensive analysis of the ANKRD2 gene and protein expression in adult chicken demonstrated its predominant expression in red muscles of thigh and drumstick, compared to white muscle. It was not detected in heart and white pectoral muscle. Uneven expression of ANKRD2 in chicken skeletal muscles, observed by immunohistochemistry, was attributed to its selective expression in slow, oxidative, type I and fast, oxidative-glycolytic, type IIA myofibers. Association of chicken ANKRD2 with phenotypic differences between red and white muscles points to its potential role in the process of myofiber-type specification. In addition to expression in slow oxidative myofibers, as demonstrated for mammalian protein, chicken ANKRD2 was also detected in fast fibers with mixed oxidative and glycolytic metabolism. This finding suggests that ANKRD2 is responsive to metabolic differences between types of avian myofibers and orientates future studies towards investigation of its role in molecular mechanisms of myofiber-type-specific gene expression.

Discovery of significant porcine SNPs for swine breed identification by a hybrid of information gain, genetic algorithm, and frequency feature selection technique

Background

The number of porcine Single Nucleotide Polymorphisms (SNPs) used in genetic association studies is very large, suitable for statistical testing. However, in breed classification problem, one needs to have a much smaller porcine-classifying SNPs (PCSNPs) set that could accurately classify pigs into different breeds. This study attempted to find such PCSNPs by using several combinations of feature selection and classification methods. We experimented with different combinations of feature selection methods including information gain, conventional as well as modified genetic algorithms, and our developed frequency feature selection method in combination with a common classification method, Support Vector Machine, to evaluate the method’s performance. Experiments were conducted on a comprehensive data set containing SNPs from native pigs from America, Europe, Africa, and Asia including Chinese breeds, Vietnamese breeds, and hybrid breeds from Thailand.

Results

The best combination of feature selection methods—information gain, modified genetic algorithm, and frequency feature selection hybrid—was able to reduce the number of possible PCSNPs to only 1.62% (164 PCSNPs) of the total number of SNPs (10,210 SNPs) while maintaining a high classification accuracy (95.12%). Moreover, the near-identical performance of this PCSNPs set to those of bigger data sets as well as even the entire data set. Moreover, most PCSNPs were well-matched to a set of 94 genes in the PANTHER pathway, conforming to a suggestion by the Porcine Genomic Sequencing Initiative.

Conclusions

The best hybrid method truly provided a sufficiently small number of porcine SNPs that accurately classified swine breeds.

Transcriptome profiling of muscle in Nelore cattle phenotypically divergent for the ribeye muscle area

This study aimed to use RNA-Seq to identify differentially expressed genes (DEGs) in muscle of uncastrated Nelore males phenotypically divergent for ribeye muscle area (REA). A total of 80 animals were phenotyped for REA, and 15 animals each with the highest REA and the lowest REA were selected for analyses. DEGs found (N = 288) belonging to families related to muscle cell growth, development, motility and proteolysis, such as actin, myosin, collagen, integrin, solute carrier, ubiquitin and kelch-like. Functional analysis showed that many of the significantly enriched gene ontology terms were closely associated with muscle development, growth, and degradation. Through co-expression network analysis, we predicted three hub genes (PPP3R1, FAM129B and UBE2G1), these genes are involved in muscle growth, proteolysis and immune system. The genes expression levels and its biological process found this study may result in differences in muscle deposition, and therefore, Nelore animals with different REA proportions.

Sequencing and characterization of lncRNAs in the breast muscle of Gushi and Arbor Acres chickens

Chicken muscle quality is one of the most important factors determining the economic value of poultry, and muscle development and growth are affected by genetics, environment, and nutrition. However, little is known about the molecular regulatory mechanisms of long non-coding RNAs (lncRNAs) in chicken skeletal muscle development. Our study aimed to better understand muscle development in chickens and thereby improve meat quality. In this study, Ribo-Zero RNA-Seq was used to investigate differences in the expression profiles of muscle development related genes and associated pathways between Gushi (GS) and Arbor Acres (AA) chickens. We identified two muscle tissue specific expression lncRNAs. In addition, the target genes of these lncRNAs were significantly enriched in certain biological processes and molecular functions, as demonstrated by Gene Ontology (GO) analysis, and these target genes participate in five signaling pathway, as revealed by an analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Taken together, these data suggest that different lncRNAs might be involved in regulating chicken muscle development and growth and provide new insight into the molecular mechanisms of lncRNAs.

Genome-wide association studies, meta-analyses and derived gene network for meat quality and carcass traits in pigs

A large number of quantitative trait loci (QTL) for meat quality and carcass traits has been reported in pigs over the past 20 years. However, few QTL have been validated and the biological meaning of the genes associated to these QTL has been underexploited. In this context, a meta-analysis was performed to compare the significant markers with meta-QTL previously reported in literature. Genome association studies were performed for 12 traits, from which 144 SNPs were found out to be significant (P < 0.05). They were validated in the meta-analysis and used to build the Association Weight Matrix, a matrix framework employed to investigate co-association of pairwise SNP across phenotypes enabling to derive a gene network. A total of 45 genes were selected from the Association Weight Matrix analysis, from which 25 significant transcription factors were identified and used to construct the networks associated to meat quality and carcass traits. These networks allowed the identification of key transcription factors, such as SOX5 and NKX2–5, gene–gene interactions (e.g. ATP5A1, JPH1, DPT and NEDD4) and pathways related to the regulation of adipose tissue metabolism and skeletal muscle development. Validated SNPs and knowledge of key genes driving these important industry traits might assist future strategies in pig breeding.

Identification of differentially expressed genes in longissimus dorsi muscle between Wei and Yorkshire pigs using RNA sequencing

Intramuscular fat (IMF) content is an important trait closely related to meat quality, which is highly variable among pig breeds from diverse genetic backgrounds. High-throughput sequencing has become a powerful technique for analyzing the whole transcription profiles of organisms. In order to elucidate the molecular mechanism underlying porcine meat quality, we adopted RNA sequencing to detect transcriptome in the longissimus dorsi muscle of Wei pigs (a Chinese indigenous breed) and Yorkshire pigs (a Western lean-type breed) with different IMF content. For the Wei and Yorkshire pig libraries, over 57 and 64 million clean reads were generated by transcriptome sequencing, respectively. A total of 717 differentially expressed genes (DEGs) were identified in our study (false discovery rate < 0.05 and fold change > 2), with 323 up-regulated and 394 down-regulated genes in Wei pigs compared with Yorkshire pigs. Gene Ontology analysis showed that DEGs significantly related to skeletal muscle cell differentiation, phospholipid catabolic process, and extracellular matrix structural constituent. Pathway analysis revealed that DEGs were involved in fatty acid metabolism, steroid biosynthesis, glycerophospholipid metabolism, and protein digestion and absorption. Quantitative real time PCR confirmed the differential expression of 11 selected DEGs in both pig breeds. The results provide useful information to investigate the transcriptional profiling in skeletal muscle of different pig breeds with divergent phenotypes, and several DEGs can be taken as functional candidate genes related to lipid metabolism (ACSL1, FABP3, UCP3 and PDK4) and skeletal muscle development (ASB2, MSTN, ANKRD1 and ANKRD2).

RNA sequencing to determine the contribution of kinase receptor transactivation to G protein coupled receptor signalling in vascular smooth muscle cells

G protein coupled receptor (GPCR) signalling covers three major mechanisms. GPCR agonist engagement allows for the G proteins to bind to the receptor leading to a classical downstream signalling cascade. The second mechanism is via the utilization of the β-arrestin signalling molecule and thirdly via transactivation dependent signalling. GPCRs can transactivate protein tyrosine kinase receptors (PTKR) to activate respective downstream signalling intermediates. In the past decade GPCR transactivation dependent signalling was expanded to show transactivation of serine/threonine kinase receptors (S/TKR). Kinase receptor transactivation enormously broadens the GPCR signalling paradigm. This work utilizes next generation RNA-sequencing to study the contribution of transactivation dependent signalling to total protease activated receptor (PAR)-1 signalling. Transactivation, assessed as gene expression, accounted for 50 percent of the total genes regulated by thrombin acting through PAR-1 in human coronary artery smooth muscle cells. GPCR transactivation of PTKRs is approximately equally important as the transactivation of the S/TKR with 209 and 177 genes regulated respectively, via either signalling pathway. This work shows that genome wide studies can provide powerful insights into GPCR mediated signalling pathways.

Whole genome scan reveals the genetic signature of African Ankole cattle breed and potential for higher quality beef

BACKGROUND

Africa is home to numerous cattle breeds whose diversity has been shaped by subtle combinations of human and natural selection. African Sanga cattle are an intermediate type of cattle resulting from interbreeding between Bos taurus and Bos indicus subspecies. Recently, research has asserted the potential of Sanga breeds for commercial beef production with better meat quality as compared to Bos indicus breeds. Here, we identified meat quality related gene regions that are positively selected in Ankole (Sanga) cattle breeds as compared to indicus (Boran, Ogaden, and Kenana) breeds using cross-population (XP-EHH and XP-CLR) statistical methods.

RESULTS

We identified 238 (XP-EHH) and 213 (XP-CLR) positively selected genes, of which 97 were detected from both statistics. Among the genes obtained, we primarily reported those involved in different biological process and pathways associated with meat quality traits. Genes (CAPZB, COL9A2, PDGFRA, MAP3K5, ZNF410, and PKM2) involved in muscle structure and metabolism affect meat tenderness. Genes (PLA2G2A, PARK2, ZNF410, MAP2K3, PLCD3, PLCD1, and ROCK1) related to intramuscular fat (IMF) are involved in adipose metabolism and adipogenesis. MB and SLC48A1 affect meat color. In addition, we identified genes (TIMP2, PKM2, PRKG1, MAP3K5, and ATP8A1) related to feeding efficiency. Among the enriched Gene Ontology Biological Process (GO BP) terms, actin cytoskeleton organization, actin filament-based process, and protein ubiquitination are associated with meat tenderness whereas cellular component organization, negative regulation of actin filament depolymerization and negative regulation of protein complex disassembly are involved in adipocyte regulation. The MAPK pathway is responsible for cell proliferation and plays an important role in hyperplastic growth, which has a positive effect on meat tenderness.

CONCLUSION

Results revealed several candidate genes positively selected in Ankole cattle in relation to meat quality characteristics. The genes identified are involved in muscle structure and metabolism, and adipose metabolism and adipogenesis. These genes help in the understanding of the biological mechanisms controlling beef quality characteristics in African Ankole cattle. These results provide a basis for further research on the genomic characteristics of Ankole and other Sanga cattle breeds for quality beef.

Genetic association of marbling score with intragenic nucleotide variants at selection signals of the bovine genome

Selection signals of Korean cattle might be attributed largely to artificial selection for meat quality. Rapidly increased intragenic markers of newly annotated genes in the bovine genome would help overcome limited findings of genetic markers associated with meat quality at the selection signals in a previous study. The present study examined genetic associations of marbling score (MS) with intragenic nucleotide variants at selection signals of Korean cattle. A total of 39 092 nucleotide variants of 407 Korean cattle were utilized in the association analysis. A total of 129 variants were selected within newly annotated genes in the bovine genome. Their genetic associations were analyzed using the mixed model with random polygenic effects based on identical-by-state genetic relationships among animals in order to control for spurious associations produced by population structure. Genetic associations of MS were found (P<3.88×10-4) with six intragenic nucleotide variants on bovine autosomes 3 (cache domain containing 1, CACHD1), 5 (like-glycosyltransferase, LARGE), 16 (cell division cycle 42 binding protein kinase alpha, CDC42BPA) and 21 (snurportin 1, SNUPN; protein tyrosine phosphatase, non-receptor type 9, PTPN9; chondroitin sulfate proteoglycan 4, CSPG4). In particular, the genetic associations with CDC42BPA and LARGE were confirmed using an independent data set of Korean cattle. The results implied that allele frequencies of functional variants and their proximity variants have been augmented by directional selection for greater MS and remain selection signals in the bovine genome. Further studies of fine mapping would be useful to incorporate favorable alleles in marker-assisted selection for MS of Korean cattle.

Recent advances in omic technologies for meat quality management

The knowledge of the molecular organization of living organisms evolved considerably during the last years. The methodologies associated also progressed with the development of the high-throughput sequencing (SNP array, RNAseq, etc.) and of genomic tools allowing the simultaneous analysis of hundreds or thousands of genes, proteins or metabolites. In farm animals, some proteins, mRNAs or metabolites whose abundance has been associated with meat quality traits have been detected in pig, cattle, chicken. They constitute biomarkers for the assessment and prediction of qualities of interest in each species, with potential biomarkers across species. The ongoing development of rapid methods will allow their use for decision-making and management tools in slaughterhouses, to better allocate carcasses or cuts to the appropriate markets. Besides, their application on living animals will help to improve genetic selection and to adapt a breeding system to fulfill expected quality level. The ultimate goal is to propose effective molecular tools for the management of product quality in meat production chains.

Identification of common regulators of genes in co-expression networks affecting muscle and meat properties

Understanding the genetic contributions behind skeletal muscle composition and metabolism is of great interest in medicine and agriculture. Attempts to dissect these complex traits combine genome-wide genotyping, expression data analyses and network analyses. Weighted gene co-expression network analysis (WGCNA) groups genes into modules based on patterns of co-expression, which can be linked to phenotypes by correlation analysis of trait values and the module eigengenes, i.e. the first principal component of a given module. Network hub genes and regulators of the genes in the modules are likely to play an important role in the emergence of respective traits. In order to detect common regulators of genes in modules showing association with meat quality traits, we identified eQTL for each of these genes, including the highly connected hub genes. Additionally, the module eigengene values were used for association analyses in order to derive a joint eQTL for the respective module. Thereby major sites of orchestrated regulation of genes within trait-associated modules were detected as hotspots of eQTL of many genes of a module and of its eigengene. These sites harbor likely common regulators of genes in the modules. We exemplarily showed the consistent impact of candidate common regulators on the expression of members of respective modules by RNAi knockdown experiments. In fact, Cxcr7 was identified and validated as a regulator of genes in a module, which is involved in the function of defense response in muscle cells. Zfp36l2 was confirmed as a regulator of genes of a module related to cell death or apoptosis pathways. The integration of eQTL in module networks enabled to interpret the differentially-regulated genes from a systems perspective. By integrating genome-wide genomic and transcriptomic data, employing co-expression and eQTL analyses, the study revealed likely regulators that are involved in the fine-tuning and synchronization of genes with trait-associated expression.

Identification and characterization of a differentially expressed protein (CAPZB) in skeletal muscle between Meishan and Large White pigs

Actin capping protein beta (CAPZB) protein was identified with considerable differences in the longissimus dorsi muscle between Large White and Meishan pigs using proteomics approach. However, in pigs, the information on CAPZB is very limited. In this study, we cloned and characterized the porcine actin capping protein beta (CAPZB) gene. In addition, we present two novel porcine CAPZB splice variants CAPZB1 and CAPZB2. CAPZB1 was expressed in all twenty tissues. However, CAPZB2 was predominantly expressed in the skeletal muscle and heart. In addition, the two isoforms had different expression profiles during the skeletal muscle development and between breeds. Moreover, the SNP T394G was identified in the coding region of the CAPZB gene, which was significantly associated with the carcass traits including the LFW, CFW, SFT and LEA. Data presented in our study suggests that the CAPZB gene may be a candidate gene of meat production trait and provides useful information for further studies on its roles in porcine skeletal muscle.

A structural approach to understanding the interactions between colour, water-holding capacity and tenderness

The colour, water-holding capacity (WHC) and tenderness of meat are primary determinants of visual and sensory appeal. Although there are many factors which influence these quality traits, the end-results of their influence is often through key changes to the structure of muscle proteins and their spatial arrangement. Water acts as a plasticiser of muscle proteins and water is lost from the myofibrillar lattice structure as a result of protein denaturation and consequent reductions in the muscle fibre volume with increasing cooking temperature. Changes in the myofilament lattice arrangement also impact the light scattering properties and the perceived paleness of the meat. Causes of variation in the quality traits of raw meat do not generally correspond to variations in cooked meat and the differences observed between the raw muscle and cooked or further processed meat are discussed. The review will also identify the gaps in our knowledge and where further investigation would beneficial.

Discovery of candidate genes for muscle traits based on GWAS supported by eQTL-analysis

Biochemical and biophysical processes that take place in muscle under relaxed and stressed conditions depend on the abundance and activity of gene products of metabolic and structural pathways. In livestock at post-mortem, these muscle properties determine aspects of meat quality and are measurable. The conversion of muscle to meat mimics pathological processes associated with muscle ischemia, injury or damage in humans and it is an economic factor in pork production. Linkage, association, and expression analyses independently contributed to the identification of trait-associated molecular pathways and genes. We aim at providing multiple evidences for the role of specific genes in meat quality by integrating a genome-wide association study (GWAS) for meat quality traits and the detection of eQTL based on trait-correlated expressed genes and trait-associated markers. The GWAS revealed 51 and 200 SNPs significantly associated with meat quality in a crossbred Pietrain×(German Landrace×Large White) (Pi×(GL×LW)) and a purebred German Landrace (GL) population, respectively. Most significant SNPs in Pi×(GL×LW) were located on chromosomes (SSC) 4 and 6. The data of 47,836 eQTLs at a significance level of p<10^-5 were used to scale down the number candidate genes located in these regions. These SNPs on SSC4 showed association with expression levels of ZNF704, IMPA1, and OXSR1; SSC6 SNPs were associated with expression of SIGLEC10 and PIH1D1. Most significant SNPs in GL were located on SSC6 and associated with expression levels of PIH1D1, SIGLEC10, TBCB, LOC100518735, KIF1B, LOC100514845, and two unknown genes. The abundance of transcripts of these genes in muscle, in turn, is significantly correlated with meat quality traits. We identified several genes with evidence for their candidacy for meat quality arising from the integrative approach of a genome-wide association study and eQTL analysis.

Correlated mRNAs and miRNAs from co-expression and regulatory networks affect porcine muscle and finally meat properties

Background

Physiological processes aiding the conversion of muscle to meat involve many genes associated with muscle structure and metabolic processes. MicroRNAs regulate networks of genes to orchestrate cellular functions, in turn regulating phenotypes.

Results

We applied weighted gene co-expression network analysis to identify co-expression modules that correlated to meat quality phenotypes and were highly enriched for genes involved in glucose metabolism, response to wounding, mitochondrial ribosome, mitochondrion, and extracellular matrix. Negative correlation of miRNA with mRNA and target prediction were used to select transcripts out of the modules of trait-associated mRNAs to further identify those genes that are correlated with post mortem traits.

Conclusions

Porcine muscle co-expression transcript networks that correlated to post mortem traits were identified. The integration of miRNA and mRNA expression analyses, as well as network analysis, enabled us to interpret the differentially-regulated genes from a systems perspective. Linking co-expression networks of transcripts and hierarchically organized pairs of miRNAs and mRNAs to meat properties yields new insight into several biological pathways underlying phenotype differences. These pathways may also be diagnostic for many myopathies, which are accompanied by deficient nutrient and oxygen supply of muscle fibers.

UBE3B and ZRANB1 polymorphisms and transcript abundance are associated with water holding capacity of porcine M. longissimus dorsi

The degradation of myofibrillar proteins during meat maturation affects the water holding capacity (WHC) of meat. Our study sought to identify polymorphisms in UBE3B and ZRANB1, genes encoding proteins involved in ubiquitination, and to evaluate the relationship between genotype, transcript abundance, and WHC of pork. A single SNP of ZRANB1, c.552A>G (p.Ile153Val), and two silent SNPs of UBE3B, c.1921A>T and c.4292C>T, were associated with muscle pH, conductivity, meat colour, or drip loss in German Landrace (GL, n=266) and Pietrain×(Large White×German Landrace) (PiF1, n=316). Further, carriers of the minor alleles at the SNPs tended to have increased transcript abundance. Consistent with the protein degradation promoting and inhibiting effects of UBE3B and ZRANB1, respectively, and the expected impact on WHC, their expressions were positively and negatively associated with WHC. The results implicate that the SNPs in both genes are in linkage with a causal site that affects transcript abundance and WHC.

Cloning, expression, and bioinformatics analysis of the sheep CARP gene

The cardiac ankyrin repeat protein (CARP) is a multifunctional protein that is expressed specifically in mammalian cardiac muscle and plays important roles in stress responses, transcriptional regulation, myofibrillar assembly, and the development of cardiac and skeletal muscle. In this study, the sheep homolog of the CARP gene was cloned and characterized. The coding region of the gene consists of 960 bp and encodes 319 amino acids with molecular weight 36.2 KD. Bioinformatics analysis demonstrated that the 3' untranslated region (3'-UTR) of the gene contains many AU-rich elements that are associated with mRNA stability and a potential regulatory site for miRNA binding. The protein was predicted to contain 14 potential phosphorylation sites and an O-GlcNAc glycosylation site and to be expressed in both the nucleus and cytoplasm. The evolutionary analysis revealed that the sheep CARP exhibited a high level of homology with the mammalian counterparts; however, the protein exhibited an increased evolutionary distance from the chicken, frog, and fish homologs. RT-PCR revealed that in addition to its high mRNA expression level in cardiac muscle, trace amounts of the sheep CARP mRNA were expressed in the skeletal muscle, stomach, and small intestine. However, western blot analysis demonstrated that the CARP protein was expressed only in cardiac muscle. The coding sequence was cloned into the pET30a-TEV-LIC vector, and the soluble CARP-MBP (maltose-binding protein) fusion protein was expressed in a prokaryotic host and purified by affinity chromatography. Our data provide the basis for future studies of the structure and function of sheep CARP.

Associations between muscle gene expression pattern and technological and sensory meat traits highlight new biomarkers for pork quality assessment

Meat quality (MQ) results from complex phenomenon and despite improved knowledge on MQ development, its variability remains high. The identification of biomarkers and the further development of rapid tests would thus be helpful to evaluate MQ in pork industries. Using transcriptomics, the present study aimed at identifying biomarkers of eight pork quality traits: ultimate pH, drip loss, lightness, redness, hue angle, intramuscular fat, shear force and tenderness, based on an experimental design inducing a high variability in MQ. Associations between microarray gene expression and pork traits (n=50 pigs) highlighted numerous potential biomarkers of MQ. Using quantitative RT-PCR, 113 transcript-trait correlations including 40 of these genes were confirmed (P<0.05, |r|≤0.73), out of which 60 were validated (P<0.05, |r|≤0.68) on complementary experimental data (n=50). Multiple regression models including 3 to 5 genes explained up to 59% of MQ trait variability. Moreover, functional analysis of correlated-trait genes provided information on the biological phenomena underlying MQ.

Elucidating molecular networks that either affect or respond to plasma cortisol concentration in target tissues of liver and muscle

Cortisol is a steroid hormone with important roles in regulating immune and metabolic functions and organismal responses to external stimuli are mediated by the glucocorticoid system. Dysregulation of the afferent and efferent axis of glucocorticoid signaling have adverse effects on growth, health status, and well-being. Glucocorticoid secretion and signaling show large interindividual variation that has a considerable genetic component; however, little is known about the underlying genetic variants. Here, we used trait-correlated expression analysis, screening for expression quantitative trait loci (eQTL), genome-wide association (GWA) studies, and causality modeling to identify candidate genes in porcine liver and muscle that affect or respond to plasma cortisol levels. Through trait-correlated expression, we characterized transcript activities in many biological functions in liver and muscle. Candidates from the list of trait-correlated expressed genes were narrowed using only those genes with an eQTL, and these were further prioritized by determining whether their expression was predicted to be related to variation in plasma cortisol levels. Using network edge orienting (NEO), a causality modeling algorithm, 26 of 990 candidates in liver were predicted to affect and 70 to respond to plasma cortisol levels. Of 593 candidates in muscle that were correlated with cortisol levels and were regulated by eQTL, 2 and 25 were predicted as effective and responsive, respectively, to plasma cortisol levels. Comprehensive data integration has helped to elucidate the complex molecular networks contributing to cortisol levels and thus its subsequent metabolic effects. The discrimination of up- and downstream effects of transcripts affecting or responding to plasma cortisol concentrations improves the understanding of the biology of complex traits related to growth, health, and well-being.

Gene expression profiling of soft and firm Atlantic salmon fillet

Texture of salmon fillets is an important quality trait for consumer acceptance as well as for the suitability for processing. In the present work we measured fillet firmness in a population of farmed Atlantic salmon with known pedigree and investigated the relationship between this trait and gene expression. Transcriptomic analyses performed with a 21 K oligonucleotide microarray revealed strong correlations between firmness and a large number of genes. Highly similar expression profiles were observed in several functional groups. Positive regression was found between firmness and genes encoding proteasome components (41 genes) and mitochondrial proteins (129 genes), proteins involved in stress responses (12 genes), and lipid metabolism (30 genes). Coefficients of determination (R(2)) were in the range of 0.64-0.74. A weaker though highly significant negative regression was seen in sugar metabolism (26 genes, R(2) = 0.66) and myofiber proteins (42 genes, R(2) = 0.54). Among individual genes that showed a strong association with firmness, there were extracellular matrix proteins (negative correlation), immune genes, and intracellular proteases (positive correlation). Several genes can be regarded as candidate markers of flesh quality (coiled-coil transcriptional coactivator b, AMP deaminase 3, and oligopeptide transporter 15) though their functional roles are unclear. To conclude, fillet firmness of Atlantic salmon depends largely on metabolic properties of the skeletal muscle; where aerobic metabolism using lipids as fuel, and the rapid removal of damaged proteins, appear to play a major role.

Gene expression analysis of mammary tissue during fetal bud formation and growth in two pig breeds--indications of prenatal initiation of postnatal phenotypic differences

BACKGROUND

The mammary gland is key to all mammal species; in particular in multiparous species like pigs the number and the shape of functional mammary gland complexes are major determinants of fitness. Accordingly, we aimed to catalog the genes relevant to mammogenesis in pigs. Moreover, we aimed to address the hypothesis that the extent and timing of proliferation, differentiation, and maturation processes during prenatal development contribute to postnatal numerical, morphological and functional properties of the mammary gland. Thus we focused on differentially expressed genes and networks relevant to mammary complex development in two breeds that are subject to different selection pressure on number, shape and function of teats and show largely different prevalence of non-functional inverted teats. The expression patterns of fetal mammary complexes obtained at 63 and 91 days post conception (dpc) from German Landrace (GL) and Pietrain (PI) were analyzed by Affymetrix GeneChip Porcine Genome Arrays.

RESULTS

The expression of 11,731 probe sets was analysed between the two stages within and among breeds. The analysis showed the largest distinction of samples of the breed GL at 63 dpc from all other samples. According to Ingenuity Pathways Analysis transcripts with abundance at the four comparisons made (GL63-GL91, PI63-PI93, GL63-PI63 and GL91-PI91) were predominantly assigned to biofunctions relevant to 'cell maintenance, proliferation, differentiation and replacement', 'organismal, organ and tissue development' and 'genetic information and nucleic acid processing'. Moreover, these transcripts almost exclusively belong to canonical pathways related to signaling rather than metabolic pathways. The accumulation of transcripts that are up-regulated in GL compared to PI indicate a higher proliferating activity in GL, whereas processes related to differentiation, maturation and maintenance of cells are more prominent in PI. Differential expression was validated by quantitative RT-PCR of five genes (GAB1, MAPK9, PIK3C2B, PIK3C3 and PRKCH) that are involved in several relevant signaling pathways.

CONCLUSIONS

The results indicate that mammary complex development in PI precedes GL. The differential expression between the two breeds at fetal stages likely reflects the prenatal initiation of postnatal phenotypes concerning the number and shape as well as functionality of teats.

Comparison of muscle transcriptome between pigs with divergent meat quality phenotypes identifies genes related to muscle metabolism and structure

BACKGROUND

Meat quality depends on physiological processes taking place in muscle tissue, which could involve a large pattern of genes associated with both muscle structural and metabolic features. Understanding the biological phenomena underlying muscle phenotype at slaughter is necessary to uncover meat quality development. Therefore, a muscle transcriptome analysis was undertaken to compare gene expression profiles between two highly contrasted pig breeds, Large White (LW) and Basque (B), reared in two different housing systems themselves influencing meat quality. LW is the most predominant breed used in pig industry, which exhibits standard meat quality attributes. B is an indigenous breed with low lean meat and high fat contents, high meat quality characteristics, and is genetically distant from other European pig breeds.

METHODOLOGY/PRINCIPAL FINDINGS

Transcriptome analysis undertaken using a custom 15 K microarray, highlighted 1233 genes differentially expressed between breeds (multiple-test adjusted P-value<0.05), out of which 635 were highly expressed in the B and 598 highly expressed in the LW pigs. No difference in gene expression was found between housing systems. Besides, expression level of 12 differentially expressed genes quantified by real-time RT-PCR validated microarray data. Functional annotation clustering emphasized four main clusters associated to transcriptome breed differences: metabolic processes, skeletal muscle structure and organization, extracellular matrix, lysosome, and proteolysis, thereby highlighting many genes involved in muscle physiology and meat quality development.

CONCLUSIONS/SIGNIFICANCE

Altogether, these results will contribute to a better understanding of muscle physiology and of the biological and molecular processes underlying meat quality. Besides, this study is a first step towards the identification of molecular markers of pork quality and the subsequent development of control tools.

Investigation of four porcine candidate genes (H-FABP, MYOD1, UCP3 and MASTR) for meat quality traits in Large White pigs

Meat quality traits are economically important traits of swine, and are controlled by multiple genes as complex quantitative traits. In the present study four genes, H-FABP (heart fatty acid-binding protein), MASTR (MEF2 activating motif and SAP domain containing transcriptional regulator), UCP3 (uncoupling protein 3) and MYOD1 (myogenic differentiation 1) were researched in Large White pigs. The polymorphisms H-FABP T/C of 5'UTR, MYOD1 g.257 A>C, UCP3 g.1406 G>A in exon 3 and MASTR c.187 C>T have been reported to be associated with meat quality traits in pigs. The aim of this study was to analyze the effect of single and multiple markers for single traits in Large White pigs. The single marker association analysis showed that the H-FABP and MASTR genes were associated with IMF (intramuscular fat content) (P < 0.05), and that the g.257 A>C of MYOD1 gene was most significantly related to muscle pH value (P < 0.01). The multiple markers for IMF were analyzed by combining the markers and quantitative trait modes into the linear regression. The results revealed that H-FABP and MASTR integrate gene networks for IMF. Thus, our study results suggested that H-FABP and MASTR polymorphisms could be used as genetic markers in the marker-assisted selection towards the improvement of IMF in Large White pigs.

Molecular characterization and different expression patterns of the muscle ankyrin repeat protein (MARP) family during porcine skeletal muscle development in vitro and in vivo

CARP, ANKRD2, and DARP belong to the ankyrin repeat protein (MARP) family and play a critical role in the integration of cytoskeletal architecture, stress response, and transcriptional regulation. In this study, we cloned the cDNA and promoter sequences of porcine ankyrin repeat protein (MARP) gene family. RT-PCR analysis revealed that porcine CARP gene was predominantly expressed in heart. ANKRD2 was widely expressed in many tissues, a high expression level was observed in the skeletal muscle and heart. DARP gene was expressed specifically in skeletal muscle and heart. Moreover, the expression of CARP and ANKRD2 was significantly different in porcine skeletal muscle among different developmental stages and between the two breeds. Expression analysis in porcine satellite cells showed that CARP and ANKRD2 were induced in differentiated porcine satellite cells, suggesting a role of them in myogenic differentiation. This result suggests that the MARP gene family may be important genes for skeletal muscle growth and provides useful information for further studies on their roles in porcine skeletal muscle.

Genes with expression levels correlating to drip loss prove association of their polymorphism with water holding capacity of pork

Six genes that were known to exhibit expression levels that are correlated to drip loss BVES, SLC3A2, ZDHHC5, CS, COQ9, and EGFR have been for candidate gene analysis. Based on in silico analysis SNPs were detected, confirmed by sequencing, and used for genotyping. The SNPs were genotyped in about 1,800 animals from six pig populations including commercial herds of Pietrain (PI) and German Landrace (DL), different commercial herds of Pietrain×(German Large White×German Landrace) (PIF1(a/b/c)), and one experimental F2-population Duroc×Pietrain (DUPI). Comparative and genetic mapping established the location of BVES on SSC1, of SLC3A2 and ZDHHC5 on SSC2, of CS on SSC5, of COQ9 on SSC6 and of EGFR on SSC9, respectively, coinciding with QTL regions for carcass and meat quality traits. BVES, SLC3A2, and CS revealed association at least with drip loss and with several other measures of water holding capacity (WHC). Moreover, COQ9 and EGFR were associated with several meat quality traits such as meat color and/or thawing loss. This study reveals statistic evidence in addition to the functional relationship of these genes to WHC previously evidenced by expression analysis. This study reveals positional and genetic statistical evidence for a link of genetic variation at these loci or close to them and promotes those six candidate genes as functional and/or positional candidate genes for meat quality traits.

Integrating expression profiling and whole-genome association for dissection of fat traits in a porcine model

Traits related to fatness, important as economic factors in pork production, are associated with serious diseases in humans. Genetical genomics is a useful approach for studying the effects of genetic variation at the molecular level in biological systems. Here we applied a whole-genome association analysis to hepatic gene expression traits, focusing on transcripts with expression levels that correlated with fatness traits in a porcine model. A total of 150 crossbred pigs [Pietrain × (German Large White × German Landrace)] were studied for transcript levels in the liver. The 24K Affymetrix expression microarrays and 60K Illumina single nucleotide polymorphism (SNP) chips were used for genotyping. A total of 663 genes, whose expression significantly correlated with the trait "fat area," were analyzed for enrichment of functional annotation groups as defined in the Ingenuity Pathways Knowledge Base (IPKB). Genes involved in metabolism of various macromolecules and nutrients as well as functions related to dynamic cellular processes correlated with fatness traits. Regions affecting the transcription levels of these genes were mapped and revealed 4,727 expression quantitative trait loci (eQTL) at P < 10⁻⁵, including 448 cis-eQTL. In this study, genome-wide association analysis of trait-correlated expression was successfully used in a porcine model to display molecular networks and list genes relevant to fatness traits.

Identification of expression QTL (eQTL) of genes expressed in porcine M. longissimus dorsi and associated with meat quality traits

BACKGROUND

Genetic analysis of transcriptional profiles is a promising approach for identifying and dissecting the genetics of complex traits like meat performance. Accordingly, expression levels obtained by microarray analysis were taken as phenotypes in a linkage analysis to map eQTL. Moreover, expression levels were correlated with traits related to meat quality and principle components with high loadings of these traits. By using an up-to-date annotation and localization of the respective probe-sets, the integration of eQTL mapping data and information of trait correlated expression finally served to point to candidate genes for meat quality traits.

RESULTS

Genome-wide transcriptional profiles of M. longissimus dorsi RNAs samples of 74 F2 animals of a pig resource population revealed 11,457 probe-sets representing genes expressed in the muscle. Linkage analysis of expression levels of these probe-sets provided 9,180 eQTL at the suggestive significance threshold of LOD > 2. We mapped 653 eQTL on the same chromosome as the corresponding gene and these were designated as 'putative cis-eQTL'. In order to link eQTL to the traits of interest, probe-sets were addressed with relative transcript abundances that showed correlation with meat quality traits at p ≤ 0.05. Out of the 653 'putative cis-eQTL', 262 transcripts were correlated with at least one meat quality trait. Furthermore, association of expression levels with composite traits with high loadings for meat quality traits generated by principle component analysis were taken into account leading to a list of 85 genes exhibiting cis-eQTL and trait dependent expression.

CONCLUSION

Holistic expression profiling was integrated with QTL analysis for meat quality traits. Correlations between transcript abundance and meat quality traits, combined with genetic positional information of eQTL allowed us to prioritise candidate genes for further study.

Functional genomics and genetical genomics approaches towards elucidating networks of genes affecting meat performance in pigs

The benefit of functional genomics is to identify key pathways and functional networks of genes and candidate genes underlying the genetic control of phenotypes. Genetical genomics, i.e. the integration of genetic analysis and expression phenotypes, has the potential to uncover regulatory networks controlling the coordinated expression of genes and to map variation on the level of DNA affecting the mRNA expression. Here we illustrate our own attempts to apply functional genomics and genetical genomics approaches in order to identify functional networks of genes relevant to traits related to meat performance. Expression data of 74 M longissimus dorsi samples obtained using Affymetrix GeneChips were correlated with drip loss and principal components (PCs) with high loadings of meat quality traits. Functional annotation analyses revealed that differences in water holding capacity, early pH decline and ultimate pH were related to the ubiquitin-proteasome system, mitochondrial metabolic pathways and muscle structural aspects. In particular, 1279 genes were correlated with drip loss (P <or= 0.001; q <or= 0.004). Negatively correlated transcripts were enriched in functional categories like extracellular matrix receptor interaction and Ca-signalling. Transcripts with a positive correlation represented oxidative phosphorylation, mitochondrial pathways and transporter activity. A linkage analysis revealed 897 expression QTL (eQTL) with 104 eQTL mapping in QTL regions for water holding capacity including 8 cis eQTL. The reduction of the multi-dimensional data sets of meat performance traits into lower dimensions of PC and the genetical genomics approach of eQTL analysis proved to be appropriate means to detect relevant biological pathways and to experimentally prioritize candidate genes.