Identification of novel mRNA isoforms associated with meat tenderness using RNA sequencing data in beef cattle

Genetic diversity of SNPs associated with candidate genes for heat stress in Coreño Creole cattle in Mexico

Mexican Coreño Creole cattle are an important genetic resource adapted to local environmental conditions, so the study of their genetic diversity is essential to know their status and implement conservation programs and their use for crossbreeding. This study evaluated the genetic diversity of heat stress tolerance characteristics of Coreño Creole cattle, and a gene ontology enrichment was performed to know the biological processes in which candidate genes are involved. A total of 48 samples from three localities of Nayarit were genotyped using 777 K Illumina BovineHD BeadChip and 34 single nucleotide polymorphisms associated with candidate genes were selected. Genetic diversity was analyzed using allelic frequencies, expected heterozygosity (He), and Wright's fixation index (FST) using PLINK v1.9 software. Candidate genes were uploaded to the open-source GOnet for pathway analysis and linkage to biological processes. Coreño Creole cattle showed low genetic diversity (He = 0.35), the average FST obtained was 0.044, and only eight markers had allele frequencies higher than 0.80 in the three locations. We found that the genes GOT1 and NCAD are related in the biological processes of stress response, cell differentiation, and homeostatic process. The results revealed that Coreño Creole cattle have low genetic diversity; this could be due to the isolation of these populations.

Novel lncRNA regulatory elements in milk somatic cells of Holstein dairy cows associated with mastitis

Despite regulatory elements such as long non - coding RNAs representing most of the transcriptome, the functional understanding of long non - coding RNAs in relation to major health conditions including bovine mastitis is limited. This study examined the milk somatic cell transcriptome from udder quarters of 6 Holstein dairy cows to identify differentially expressed long non - coding RNAs using RNA - Sequencing. Ninety - four differentially expressed long non - coding RNAs are identified, 5 of which are previously annotated for gene name and length, 11 are annotated for gene name and 78 are novel, having no gene name or length previously annotated. Significant inflammatory response and regulation of immune response pathways (false discovery rate < 0.05) are associated with the differentially expressed long non - coding RNAs. QTL annotation analysis revealed 31 QTL previously annotated in the genomic regions of the 94 differentially expressed long non - coding RNAs, and the majority are associated with milk traits. This research provides a better understanding of long non - coding RNAs regulatory elements in milk somatic cells, which may enhance current breeding strategies for more adaptable or high mastitis resistant cattle.

Physiological and transcriptomic analysis dissects the molecular mechanism governing meat quality during postmortem aging in Hu sheep (Ovis aries)

Introduction

Hu sheep, known for its high quality and productivity, lack fundamental scientific research in China.

Methods

This study focused on the effects of 24 h postmortem aging on the meat physiological and transcriptomic alteration in Hu sheep.

Results

The results showed that the 24 h aging process exerts a substantial influence on the mutton color, texture, and water content as compared to untreated group. Transcriptomic analysis identified 1,668 differentially expressed genes. Functional enrichment analysis highlighted the importance of glycolysis metabolism, protein processing in endoplasmic reticulum, and the FcγR-mediated phagocytosis pathway in mediating meat quality modification following postmortem aging. Furthermore, protein-protein interaction analysis uncovered complex regulatory networks involving glycolysis, the MAPK signaling pathway, protein metabolism, and the immune response.

Discussion

Collectively, these findings offer valuable insights into the molecular mechanisms underlying meat quality changes during postmortem aging in Hu sheep, emphasizing the potential for improving quality control strategies in mutton production.

Identification of Novel mRNA Isoforms Associated with Acute Heat Stress Response Using RNA Sequencing Data in Sprague Dawley Rats

The molecular mechanisms underlying heat stress tolerance in animals to high temperatures remain unclear. This study identified the differentially expressed mRNA isoforms which narrowed down the most reliable DEG markers and molecular pathways that underlie the mechanisms of thermoregulation. This experiment was performed on Sprague Dawley rats housed at 22 °C (control group; CT), and three acute heat-stressed groups housed at 42 °C for 30 min (H30), 60 min (H60), and 120 min (H120). Earlier, we demonstrated that acute heat stress increased the rectal temperature of rats, caused abnormal changes in the blood biochemical parameters, as well as induced dramatic changes in the expression levels of genes through epigenetics and post-transcriptional regulation. Transcriptomic analysis using RNA-Sequencing (RNA-Seq) data obtained previously from blood (CT and H120), liver (CT, H30, H60, and H120), and adrenal glands (CT, H30, H60, and H120) was performed. The differentially expressed mRNA isoforms (DEIs) were identified and annotated by the CLC Genomics Workbench. Biological process and metabolic pathway analyses were performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 225, 5764, and 4988 DEIs in the blood, liver, and adrenal glands were observed. Furthermore, the number of novel differentially expressed transcript lengths with annotated genes and novel differentially expressed transcript with non-annotated genes were 136 and 8 in blood, 3549 and 120 in the liver, as well as 3078 and 220 in adrenal glands, respectively. About 35 genes were involved in the heat stress response, out of which, Dnaja1, LOC680121, Chordc1, AABR07011951.1, Hsp90aa1, Hspa1b, Cdkn1a, Hmox1, Bag3, and Dnaja4 were commonly identified in the liver and adrenal glands, suggesting that these genes may regulate heat stress response through interactions between the liver and adrenal glands. In conclusion, this study would enhance our understanding of the complex underlying mechanisms of acute heat stress, and the identified mRNA isoforms and genes can be used as potential candidates for thermotolerance selection in mammals.

Bibliometric Analysis of Research on the Main Genes Involved in Meat Tenderness

Simple Summary

A bibliometric analysis was carried out to know the evolution of research on genes associated with meat tenderness, of interest for the development of selection programs. Since 1993, studies have been limited to a few researchers in high-income countries due to the costs associated with the techniques. The main findings showed that the scientific production had a discontinuous growth because science experienced a significant change since approximately 2010. Marker-assisted selection was widely used, evaluating mainly CAPN (calpain) and CAST (calpastatin) genes for their contribution to meat tenderness, especially in cattle. However, the effects are small; therefore, genomic selection was implemented by genotyping thousands of single nucleotide polymorphisms (SNPs) for further explanation of genetic variation. The results shown are important for scholars to identify emerging methodologies and gaps in the literature and to know who the prolific authors and institutions in the field for possible collaborations, etc., are.

Abstract

Tenderness is one of the main characteristics of meat because it determines its price and acceptability. This is the first bibliometric study on the trend of research on the role of genes in meat tenderness. A total of 175 original and English-language articles published up to 2021 were retrieved from Scopus. The bibliometric analysis was carried out with VOSviewer (version 1.6.18, Eck and Waltman, Leiden, Netherlands) and complemented with the Analyze search results service from Scopus. Erroneous and duplicate data were eliminated, and incomplete information was added to standardize the results. Scientific production was evaluated by means of quantity, quality and structure indicators. As a first glance, 8.816% of authors have published more than 50% of papers mainly related to genes encoding the calpain (CAPN)-calpastatin (CAST) system and single nucleotide polymorphisms (SNPs). Among other findings, a strong link was found between the contribution of the main countries (led by the United States with) and their institutions (led by the USDA Agricultural Research Service with) to their gross domestic product. Most studies on the topic are published in the Journal of Animal Science, and other journals with high impact according to the number of citations and different metrics. Finally, when evaluating the most cited articles, the occurrence and association of the main keywords, it was confirmed that research is focused on the role of CAPN and CAST genes and of SNPs in beef tenderness. The change in science was emphasized; although marker-assisted selection is still used, genes have an infinitesimal effect on complex traits. Therefore, since about 2010, new research groups adopted genomic selection to evaluate dense panels of SNPs and better explain genetic variation in meat tenderness.

Multi-omic data integration for the study of production, carcass, and meat quality traits in Nellore cattle

Data integration using hierarchical analysis based on the central dogma or common pathway enrichment analysis may not reveal non-obvious relationships among omic data. Here, we applied factor analysis (FA) and Bayesian network (BN) modeling to integrate different omic data and complex traits by latent variables (production, carcass, and meat quality traits). A total of 14 latent variables were identified: five for phenotype, three for miRNA, four for protein, and two for mRNA data. Pearson correlation coefficients showed negative correlations between latent variables miRNA 1 (mirna1) and miRNA 2 (mirna2) (−0.47), ribeye area (REA) and protein 4 (prot4) (−0.33), REA and protein 2 (prot2) (−0.3), carcass and prot4 (−0.31), carcass and prot2 (−0.28), and backfat thickness (BFT) and miRNA 3 (mirna3) (−0.25). Positive correlations were observed among the four protein factors (0.45–0.83): between meat quality and fat content (0.71), fat content and carcass (0.74), fat content and REA (0.76), and REA and carcass (0.99). BN presented arcs from the carcass, meat quality, prot2, and prot4 latent variables to REA; from meat quality, REA, mirna2, and gene expression mRNA1 to fat content; from protein 1 (prot1) and mirna2 to protein 5 (prot5); and from prot5 and carcass to prot2. The relations of protein latent variables suggest new hypotheses about the impact of these proteins on REA. The network also showed relationships among miRNAs and nebulin proteins. REA seems to be the central node in the network, influencing carcass, prot2, prot4, mRNA1, and meat quality, suggesting that REA is a good indicator of meat quality. The connection among miRNA latent variables, BFT, and fat content relates to the influence of miRNAs on lipid metabolism. The relationship between mirna1 and prot5 composed of isoforms of nebulin needs further investigation. The FA identified latent variables, decreasing the dimensionality and complexity of the data. The BN was capable of generating interrelationships among latent variables from different types of data, allowing the integration of omics and complex traits and identifying conditional independencies. Our framework based on FA and BN is capable of generating new hypotheses for molecular research, by integrating different types of data and exploring non-obvious relationships.

Interaction of MyoD and MyoG with Myoz2 gene in bovine myoblast differentiation

This study aims to explore the functional role of Myoz2 in myoblast differentiation, and elucidate the potential factors interact with Myoz2 in promoter transcriptional regulation. The temporal-spatial expression results showed that the bovine Myoz2 gene was highest expressed in longissimus dorsi, and in individual growth stages and myoblast differentiation stages. Knockdown of Myoz2 inhibited the differentiation of myoblast, and negative effect of MyoD, MyoG, MyH and MEF2A expression on mRNA levels. Subsequently, the promoter region of bovine Myoz2 gene with 1.7 Kb sequence was extracted, and then it was set as eight series of deleted fragments, which were ligated into pGL3-basic to detect core promoter regions of Myoz2 gene in myoblasts and myotubes. Transcription factors MyoD and MyoG were identified as important cis-acting elements in the core promoter region (-159/+1). Also, it was highly conserved in different species based on dual-luciferase analysis and multiple sequence alignment analysis, respectively. Furthermore, a chromatin immunoprecipitation (ChIP) analysis combined with site-directed mutation and siRNA interference and overexpression confirmed that the combination of MyoD and MyoG occurred in region -159/+1, and played an important role in the regulation of bovine Myoz2 gene. These findings explored the regulatory network mechanism of Myoz2 gene during the development of bovine skeletal muscle.

Use of a graph neural network to the weighted gene co-expression network analysis of Korean native cattle

In the general framework of the weighted gene co-expression network analysis (WGCNA), a hierarchical clustering algorithm is commonly used to module definition. However, hierarchical clustering depends strongly on the topological overlap measure. In other words, this algorithm may assign two genes with low topological overlap to different modules even though their expression patterns are similar. Here, a novel gene module clustering algorithm for WGCNA is proposed. We develop a gene module clustering network (gmcNet), which simultaneously addresses single-level expression and topological overlap measure. The proposed gmcNet includes a “co-expression pattern recognizer” (CEPR) and “module classifier”. The CEPR incorporates expression features of single genes into the topological features of co-expressed ones. Given this CEPR-embedded feature, the module classifier computes module assignment probabilities. We validated gmcNet performance using 4,976 genes from 20 native Korean cattle. We observed that the CEPR generates more robust features than single-level expression or topological overlap measure. Given the CEPR-embedded feature, gmcNet achieved the best performance in terms of modularity (0.261) and the differentially expressed signal (27.739) compared with other clustering methods tested. Furthermore, gmcNet detected some interesting biological functionalities for carcass weight, backfat thickness, intramuscular fat, and beef tenderness of Korean native cattle. Therefore, gmcNet is a useful framework for WGCNA module clustering.

Preliminary Results about Lamb Meat Tenderness Based on the Study of Novel Isoforms and Alternative Splicing Regulation Pathways Using Iso-seq, RNA-seq and CTCF ChIP-seq Data

Tenderness is an important indicator of meat quality. Novel isoforms associated with meat tenderness and the role of the CCCTC-binding factor (CTCF) in regulating alternative splicing to produce isoforms in sheep are largely unknown. The current project studied six sheep from two crossbred populations (Dorper × Hu × Hu, DHH and Dorper × Dorper × Hu, DDH) with divergent meat tenderness. Pooled Iso-seq data were used to annotate the sheep genomes. Then, the updated genome annotation and six RNA-seq data were combined to identify differentially expressed isoforms (DEIs) in muscles between DHH and DDH. These data were also combined with peaks detected from CTCF ChIP-seq data to investigate the regulatory role of CTCF for the alternative splicing. As a result, a total of 624 DEIs were identified between DDH and DHH. For example, isoform 7.524.18 transcribed from CAPN3 may be associated with meat tenderness. In addition, a total of 86 genes were overlapped between genes with transcribed DEIs and genes in differential peaks identified by CTCF ChIP-seq. Among these overlapped genes, ANKRD23 produces different isoforms which may be regulated by CTCF via methylation. As preliminary research, our results identified novel isoforms associated with meat tenderness and revealed the possible regulating mechanisms of alternative splicing to produce isoforms.

Identification of novel alternative splicing associated with mastitis disease in Holstein dairy cows using large gap read mapping

Background

Mastitis is a very common disease in the dairy industry that producers encounter daily. Transcriptomics, using RNA-Sequencing (RNA - Seq) technology, can be used to study the functional aspect of mastitis resistance to identify animals that have a better immune response to mastitis. When the cow has mastitis, not only genes but also specific mRNA isoforms generated via alternative splicing (AS) could be differentially expressed (DE), leading to the phenotypic variation observed. Therefore, the objective of this study was to use large gap read mapping to identify mRNA isoforms DE between healthy and mastitic milk somatic cell samples (N = 12). These mRNA isoforms were then categorized based on being 1) annotated mRNA isoforms for gene name and length, 2) annotated mRNA isoforms with different transcript length and 3) novel mRNA isoforms of non - annotated genes.

Results

Analysis identified 333 DE transcripts (with at least 2 mRNA isoforms annotated, with at least one being DE) between healthy and mastitic samples corresponding to 303 unique genes. Of these 333 DE transcripts between healthy and mastitic samples, 68 mRNA isoforms are annotated in the bovine genome reference (ARS.UCD.1.2), 249 mRNA isoforms had novel transcript lengths of known genes and 16 were novel transcript lengths of non - annotated genes in the bovine genome reference (ARS.UCD.1.2). Functional analysis including gene ontology, gene network and metabolic pathway analysis was performed on the list of 288 annotated and unique DE mRNA isoforms. In total, 67 significant metabolic pathways were identified including positive regulation of cytokine secretion and immune response. Additionally, numerous DE novel mRNA isoforms showed potential involvement with the immune system or mastitis. Lastly, QTL annotation analysis was performed on coding regions of the DE mRNA isoforms, identifying overlapping QTLs associated with clinical mastitis and somatic cell score.

Conclusion

This study identified novel mRNA isoforms generated via AS that could lead to differences in the immune response of Holstein dairy cows and be potentially implemented in future breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08430-x.

Characterization of novel lncRNA muscle expression profiles associated with meat quality in beef cattle

The aim of this study was to identify novel lncRNA differentially expressed (DE) between divergent animals for beef tenderness and marbling traits in Nellore cattle. Longissimus thoracis muscle samples from the 20 most extreme bulls (of 80 bulls set) for tenderness, tender (n = 10) and tough (n = 10) groups, and marbling trait, high (n = 10) and low (n = 10) groups were used to perform transcriptomic analysis using RNA‐Sequencing. For tenderness, 29 lncRNA were DE (p‐value ≤ 0.01) in tough beef animals in relation to tender beef animals. We observed that genic lncRNAs, for example, lncRNA_595.1, were overlapping exonic part of the PICK gene, while lncRNA_3097.2 and lncRNA_3129.5 overlapped intronic part of the genes GADL1 and PSMD6. The lncRNA associated with PICK1, GADL1, and PMD6 genes were enriched in the pathways associated with the ionotropic glutamate receptor, gamma‐aminobutyric acid synthesis, and the ubiquitin–proteasome pathway. For marbling, 50 lncRNA were DE (p‐value ≤ 0.01) in high marbling group compared with low marbling animals. The genic lncRNAs, such as lncRNA_3191.1, were overlapped exonic part of the ITGAL gene, and the lncRNA_512.1, lncRNA_3721.1, and lncRNA_41.4 overlapped intronic parts of the KRAS and MASP1 genes. The KRAS and ITGAL genes were enriched in pathways associated with integrin signaling, which is involved in intracellular signals in response to the extracellular matrix, including cell form, mobility, and mediates progression through the cell cycle. In addition, the lincRNAs identified to marbling trait were associated with several genes related to calcium binding, muscle hypertrophy, skeletal muscle, lipase, and oxidative stress response pathways that seem to play a role important in the physiological processes related to meat quality. These findings bring new insights to better understand the biology mechanisms involved in the gene regulation of these traits, which will be valuable for a further investigation of the interactions between lncRNA and mRNAs, and of how these interactions may affect meat quality traits.

Effects of Neck-Arm Restraint Suspension of Beef Carcasses on Meat Quality and Proteome of Different Muscles During Post-mortem Aging

Beef quality is the first deciding factor for consumers to consider before purchasing. The aim of this study was to evaluate the effects of suspension and aging time on beef quality. We compared the differences in pH, drip loss, cooking loss, color, shear force, myofibril fragmentation index (MFI), and electron microscope of three muscle tissues between Achilles tendon (AT) and neck-arm restraint (NR) suspensions during seven aging periods (days 0, 1, 2, 3, 7, 14, and 21) after slaughter using the carcasses of six Xinjiang brown cattle. We found that NR suspension could significantly increase the water loss rate and MFI, as well as reduce the shear force compared to AT suspension. The muscle fiber structure with NR suspension was more severely damaged. The proteomics of longissimus dorsi was checked for the post-mortem days 1, 7, and 14. We detected 50, 26, and 29 differentially expressed proteins between NR and AT suspension at post-mortem days 1, 7, and 14, respectively. These proteins were involved in metabolic and muscle structure associated pathways and contributed to a comprehensive understanding of suspension-dependent meat quality regulation by proteins in beef cattle. To conclude, NR suspension can accelerate the aging time of beef carcasses, which will reduce the cost of carcass suspension and bring more benefits in the beef industry.

Genome-wide association study for meat tenderness in beef cattle identifies patterns of the genetic contribution in different post-mortem stages

The beef tenderization process during the post-mortem period is one of the most important sensorial attributes and it is well-established. The aim of this study was to identify the genetic contribution pattern to meat tenderness at 7-(LMD7), 14-(LMD14), and 21-(LMD21) days post-mortem. The heritabilities for LMD7 (0.194), LMD14 (0.142) and LMD21 (0.048) are well established in the population evaluated here. However, its genetic contribution in terms of genomic candidate regions is still poorly understood. Tenderness was measured in the Longissiums thoracis using Warner-Bratzler shear force in the three post-mortem periods. A total of 4323 crossbred beef cattle were phenotyped and genotyped using the Illumina BovineSNP50K. The percentage of the total genetic variance was estimated using the weighted single-step genomic best linear unbiased prediction method. The main candidate windows for LMD7 were associated with proteolysis of myofibrillar structures and the weakening endomysium and perimysium. Candidate windows for LMD14 and LMD21 were mapped in bovine QTLs for body composition, height and growth. Results presented herein highlight, the largest contribution of proteolysis related processes before 14-days post-mortem and body composition characteristics in later stages for meat tenderness.

RNA Profiles of the Korat Chicken Breast Muscle with Increased Carnosine Content Produced through Dietary Supplementation with β-Alanine or L-Histidine

Simple Summary

Carnosine is a bioactive food component with several potential health benefits for humans due to its physiological functions. Dietary supplementation with β-alanine or L-histidine can increase the carnosine content of skeletal muscles in chickens. Dietary supplementation with β-alanine or L-histidine has produced a slow-growing chicken variety with high carnosine content in the breast meat; however, the supplementation with L-histidine alone softens the meat toughness, which may affect consumers’ willingness to buy the meat. Gene expression is a key factor that influences meat quality. Understanding the molecular mechanisms that affect carnosine content and meat toughness would allow the production of more value-added slow-growing chickens. We compared global gene expression in chicken breast muscles with differing carnosine contents and meat toughness produced through dietary supplementation with β-alanine or L-histidine. We identified differentially expressed genes involved in regulating myosin, collagen, intramuscular fat, and calpain—factors that may affect meat tenderness. Pathway enrichment analysis indicated that the insulin-related and adipocytokine signaling pathways were altered by dietary supplementation with β-alanine or L-histidine. These data will be useful for future studies on carnosine content and meat toughness in slow-growing chickens.

Abstract

Korat chicken (KRC) is a slow-growing chicken bred in Thailand, whose meat exhibits a unique toughness. A previous study produced KRC breast meat containing high carnosine content through dietary supplementation with β-alanine or L-histidine; however, the KRC that were fed an L-histidine-supplemented diet produced meat that was significantly more tender. Herein, we performed RNA-Seq to identify candidate genes involved in the regulation of carnosine content and meat toughness. Total RNA was isolated from five female KRC breast muscles in each treatment group that KRC fed diets without supplementation, supplemented with β-alanine or L-histidine. Compared to the non-supplemented group, we identified 118 and 198 differentially expressed genes (DEGs) in the β-alanine or L-histidine supplementation groups, respectively. Genes potentially related to meat tenderness—i.e., those regulating myosin, collagen, intramuscular fat, and calpain—were upregulated (LOC107051274, ACSBG1, and CAPNS2) and downregulated (MYO7B, MYBPH, SERPINH1, and PGAM1). However, carnosine synthase gene was not identified. Functional enrichment analysis identified pathways affected by dietary supplementation, including the insulin signaling pathway (β-alanine supplementation) and the insulin resistance and adipocytokine signaling pathways (L-histidine supplementation). The FoxO signaling pathway was identified as a regulatory network for both supplementation groups. The identified genes can be used as molecular markers of meat tenderness in slow-growing chickens.