Reference gene selection and myosin heavy chain (MyHC) isoform expression in muscle tissues of domestic yak (Bos grunniens)

Myosin heavy chain isoform expression and meat quality characteristics of different muscles in yak (Bos grunniens)

Myosin heavy chain (MHC) isoforms and meat quality characteristics of different muscles were investigated to explore their potential relationships in yaks. Results showed that semitendinosus (ST), longissimus thoracis (LT), and infraspinatus (IS) have a greater ratio of MHC IIb (47.84%), MHC IIa (73.27%), and MHC I (24.26%), respectively, than the other two muscles. Compared with LT or ST, IS exhibited more intense color, greater water-holding capacity, and initial tenderness with higher intermuscular fat (IMF) and collagen (of lower cross-linking level), presenting overall better quality. Variations in MHC isoforms accounted for the muscle-specific meat quality. Specifically, MHC I was positively associated with redness, myoglobin, IMF, collagen, pH, and thermal stability and negatively associated with myofibril fragmentation index, fiber thickness, collagen cross-linking, and drip loss. These results provide insights into the relationships between MHC isoforms and meat quality in yaks and the MHC I isoform has an extensive influence on meat quality.

Effect of slaughter age and postmortem aging time on tenderness and water-holding capacity of yak (Bos grunniens) longissimus thoracis muscle

The present study investigated the effect of slaughter age (2.43 ± 0.20, 4.15 ± 0.19, 6.62 ± 0.18, 10.59 ± 0.74 years) and postmortem aging time (1, 24, and 72 h) on the tenderness and water-holding capacity (WHC) of yak longissimus thoracis muscles to determine the most suitable age for slaughter to ensure product consistency. Under conventional postmortem aging conditions (4 °C), muscles of each age group exhibited the effect of cold shortening. Once the cold shortening occurred, the age effect on thickening muscle fiber and developing cross-links of collagen, considered to intensify the meat toughness, became less important. Owing to greater carcass weight and intramuscular fat, muscles of the older carcass (over 6-year-old) were less influenced by the cold shortening effect during the chilling process and showed lessened sarcomere contraction, delayed formation of drip loss channels, and increased level of myofibril fragmentation index (MFI) and myofiber structural disintegration, resulting in greater tenderness and WHC, especially 6-7 years group. Aging of 72 h structurally disintegrated the collagen cross-linking and integrity of muscle fibers and elevated the MFI, improving the meat tenderness. Therefore, the suitable slaughter age for yak is 6-7 years old and after 72 h aging, improved quality of yak meat can be obtained.

Transcriptome-Based Evaluation of Optimal Reference Genes for Quantitative Real-Time PCR in Yak Stomach throughout the Growth Cycle

Efficient nutritional assimilation and energy metabolism in the stomachs of yaks contribute to their adaption to harsh environments. Accurate gene expression profile analysis will help further reveal the molecular mechanism of nutrient and energy metabolism in the yak stomach. RT-qPCR is regarded as an accurate and dependable method for analyzing gene expression. The selection of reference genes is essential to obtain meaningful RT-qPCR results, especially in longitudinal gene expression studies of tissues and organs. Our objective was to select and validate optimal reference genes from across the transcriptome as internal controls for longitudinal gene expression studies in the yak stomach. In this study, 15 candidate reference genes (CRGs) were determined according to transcriptome sequencing (RNA-seq) results and the previous literature. The expression levels of these 15 CRGs were quantified using RT-qPCR in the yak stomach, including the rumen, reticulum, omasum and abomasum at five stages: 0 days, 20 days, 60 days, 15 months and three years old (adult). Subsequently, the expression stabilities of these 15 CRGs were evaluated via four algorithms: geNorm, NormFinder, BestKeeper and the comparative C_T method. Furthermore, RefFinder was employed to obtain a comprehensive ranking of the stability of CRGs. The analysis results indicate that RPS15, MRPL39 and RPS23 are the most stable genes in the yak stomach throughout the growth cycle. In addition, to verify the reliability of the selected CRGs, the relative expression levels of HMGCS2 were quantified via RT-qPCR using the three most stable or the three least stable CRGs. Overall, we recommend combining RPS15, MRPL39 and RPS23 as reference genes for the normalization of RT-qPCR data in the yak stomach throughout the growth cycle.

Importance of Housekeeping Gene Optimization for the Analysis of mRNA Expression During Wound Healing in a Third-Degree Burn Injury Model

Wound healing evaluation methods in a third-degree burn injury model are categorized as histological (re-epithelialization and granulation tissue formation) and molecular (quantitative polymerase chain reaction). In general, mRNA expression is normalized to those of the housekeeping gene. Although the housekeeping gene expression is generally stable, it has been reported that the stability of these genes depends on the wound healing process and treatment method. In this study, we identified the most stable housekeeping gene (TATA-binding protein) for studying gene expression in a third-degree burn injury model, in which wound healing was promoted by grafting human amnion-derived mesenchymal cells. We investigated the wound healing effect of human amnion-derived mesenchymal cells in the injury model. The formation of granulation tissue, the differentiation from fibroblasts to myofibroblasts, and functional vascular structure were promoted in the full-thickness skin excision site by treatment with these cells. The expression of angiogenic, pro-inflammatory and anti-inflammatory related mRNA was measured and normalized to that of the housekeeping gene, showing that treatment with the cells promoted the infiltration of endothelial cells and differentiation of M1 and M2 macrophages. In conclusion, wound healing in a third-degree burn injury model can be accurately analyzed using the optimized housekeeping gene.

Whole transcriptome analyses and comparison reveal the metabolic differences between oxidative and glycolytic skeletal muscles of yak

Mammalian skeletal muscle is composed of various muscle fibers that exhibit different physiological and metabolic features. Muscle fiber type composition has significant influences on the meat quality of livestock. In this study, we comprehensively analyzed the whole transcriptome profiles of the oxidative muscle biceps femoris (BF) and the glycolytic muscle obliquus externus abdominis (OEA) of yak. A total of 1436 mRNAs, 1172 lncRNAs, and 218 circRNAs were differentially expressed in the oxidative muscles compared with the glycolytic muscles. KEGG annotation showed that differentially expressed mRNAs regulated by lncRNA and circRNA were mainly involved in PPAR signaling pathway, citrate cycle (TCA cycle), and PI3K-Akt signaling pathway, which reflect the different metabolic properties between oxidative and glycolytic muscles. In addition, regulatory networks associated with muscle fiber type conversion and mitochondria energy metabolism in muscles were constructed. Our study provides new evidence for a better understanding of the molecular mechanisms underlying skeletal muscle fiber determination and meat quality traits of yak.

Mutation of the MYH3 gene causes recessive cleft palate in Limousine cattle

Background

The palate is a structure separating the oral and nasal cavities and its integrity is essential for feeding and breathing. The total or partial opening of the palate is called a cleft palate and is a common malformation in mammals with environmental or hereditary aetiologies. Generally, it compromises life expectancy in the absence of surgical repair. A new form of non-syndromic cleft palate arose recently in Limousine cattle, with animals referred to the French National Observatory of Bovine Abnormalities since 2012. Since the number of affected animals has increased steadily, this study was undertaken to identify the cause of this disease.

Results

Based on pedigree analysis, occurrence of cleft palate in Limousine cattle was concordant with an autosomal recessive mode of inheritance. Genotyping of 16 affected animals and homozygosity mapping led to the identification of a single disease-associated haplotype on Bos taurus chromosome (BTA)19. The genome of two affected animals was sequenced, and their sequences were compared to the ARS-UCD1.2 reference genome to identify variants. The likely causal variants were compared to the variant database of the 1000 bull genome project and two fully linked mutations in exon 24 of the MYH3 (myosin heavy chain) gene were detected: a 1-bp non-synonymous substitution (BTA19:g.29609623A>G) and a 11-bp frameshift deletion (BTA19:g.29609605-29609615del). These two mutations were specific to the Limousine breed, with an estimated allele frequency of 2.4% and are predicted to be deleterious. The frameshift leads to a premature termination codon. Accordingly, mRNA and protein analyses in muscles from wild-type and affected animals revealed a decrease in MYH3 expression in affected animals, probably due to mRNA decay, as well as an absence of the MYH3 protein in these animals. MYH3 is mostly expressed in muscles, including craniofacial muscles, during embryogenesis, and its absence may impair palate formation.

Conclusions

We describe a new form of hereditary cleft palate in Limousine cattle. We identified two fully linked and deleterious mutations, ultimately leading to the loss-of-function of the MYH3 protein. The mutations were included on the Illumina EuroG10k v8 and EuroGMD v1 SNP chips and are used to set up a reliable eradication strategy in the French Limousine breed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00762-2.

Reference gene selection in bovine caruncular epithelial cells under pregnancy-associated hormones exposure

Examination of transcriptional regulation occurring during pregnancy establishment and maintenance requires the identification of endogenous reference genes characterized by high expression stability. Since the expression of some reference genes may be modulated by pregnancy-associated hormones, the goal of our study was to identify suitable reference genes unaffected by hormonal treatment. In our study bovine caruncular epithelial cells were subjected to progesterone, estrogen and prostaglandin F_2α treatment. Ten candidate reference genes (ACTR1A, CNOT11, HDAC1, HPRT1, RPL19, RPS9, SDHA, SUZ12, UXT and ZNF131) were evaluated with the use of four approaches (geNorm, NormFinder, BestKeeper, delta Ct). We found that RPS9 and SUZ12 displayed the highest expression stability in the tested material. Moreover, HPRT1 and SDHA were found inappropriate for RT-qPCR data normalization as they demonstrated the highest expression variability out of all candidates analysed. Hence geNorm calculations shown that the use of just two best-performing genes would be sufficient for obtaining reliable results, we propose that RPS9 and SUZ12 be used as suitable endogenous controls in future studies investigating gene expression in normal and compromised pregnancies.

Genome-Wide Identification of Reference Genes for Reverse-Transcription Quantitative PCR in Goat Rumen

Simple Summary

The rumen plays an essential role as a digestive organ and serves as the primary site of energy substrate absorption for the productive ruminants. Understanding gene expression profiles is necessary to explore the intrinsic regulatory mechanisms of rumen development in goats. The selection of suitable reference genes (RGs) was the primary assay before the real-time quantitative PCR (RT-qPCR). We identified sixteen genome-wide candidate RGs for normalization of gene expression assessments in goat rumen tissues. We demonstrate that the RGs selected (RPS4X and RPS6) were more stably expressed than the commonly used HKGs (ACTB and GAPDH) in goat rumen tissues, suggesting that the ribosomal protein gene family may be another source for the RG pool.

Abstract

As the largest chamber of the ruminant stomach, the rumen not only serves as the principal absorptive surface and nutrient transport pathway from the lumen into the animal, but also plays an important short-chain fatty acid (SCFA) metabolic role in addition to protective functions. Accurate characterization of the gene expression profiles of genes of interest is essential to the exploration of the intrinsic regulatory mechanisms of rumen development in goats. Thus, the selection of suitable reference genes (RGs) is an important prerequisite for real-time quantitative PCR (RT-qPCR). In the present study, 16 candidate RGs were identified from our previous transcriptome sequencing of caprine rumen tissues. The quantitative expressions of the candidate RGs were measured using the RT-qPCR method, and the expression stability of the RGs was assessed using the geNorm, NormFinder, and BestKeeper programs. GeNorm analysis showed that the M values were less than 0.5 for all the RGs except GAPT4, indicating that they were stably expressed in the rumen tissues throughout development. RPS4X and RPS6 were the two most stable RGs. Furthermore, the expressions of two randomly selected target genes (IGF1 and TOP2A), normalized by the selected most stable RGs (RPS4X and RPS6), were consistent with the results of RNA sequencing, while the use of GAPDH and ACTB as RGs resulted in altered profiles. Overall, RPS4X and RPS6 showed the highest expression stability and the lowest coefficients of variation, and could be used as the optimal reference combination for quantifying gene expression in rumen tissues via RT-qPCR analysis.

Bovine Satellite Cells Isolated after 2 and 5 Days of Tissue Storage Maintain the Proliferative and Myogenic Capacity Needed for Cultured Meat Production

Cultured meat is an emerging alternative food technology which aims to deliver a more ethical, sustainable, and healthy muscle-tissue-derived food item compared to conventional meat. As start-up companies are rapidly forming and accelerating this technology, many aspects of this multi-faceted science have still not been investigated in academia. In this study, we investigated if bovine satellite cells with the ability to proliferate and undergo myogenic differentiation could be isolated after extended tissue storage, for the purpose of increasing the practicality for cultured meat production. Proliferation of bovine satellite cells isolated on the day of arrival or after 2 and 5 days of tissue storage were analyzed by metabolic and DNA-based assays, while their myogenic characteristics were investigated using RT-qPCR and immunofluorescence. Extended tissue storage up to 5 days did not negatively affect proliferation nor the ability to undergo fusion and create myosin heavy chain-positive myotubes. The expression patterns of myogenic and muscle-specific genes were also not affected after tissue storage. In fact, the data indicated a positive trend in terms of myogenic potential after tissue storage, although it was non-significant. These results suggest that the timeframe of which viable myogenic satellite cells can be isolated and used for cultured meat production can be greatly extended by proper tissue storage.

Rational, Unbiased Selection of Reference Genes for Pluripotent Stem Cell-Derived Cardiomyocytes

Reverse transcription, quantitative polymerase chain reaction (RT-qPCR) is a powerful technique to quantify gene expression by transcript abundance. Expression of target genes is normalized to expression of stable reference genes to account for sample preparation variability. Thus, the identification and validation of stably expressed reference genes is crucial for making accurate, quantitative, statistical conclusions in gene expression studies. Traditional housekeeping genes identified decades ago based on high and relatively stable expression are often used, although many have shown these to not be valid, particularly in highly dynamic systems such as stem cell differentiation. In this study we outline a rational approach to identify stable reference genes valid throughout human pluripotent stem cell (hPSC) differentiation to hPSC-derived cardiomyocytes (hPSC-CMs). Several publicly available transcriptomic data sets were analyzed to identify genes with low variability in expression throughout differentiation. These putative novel reference genes were subsequently validated in RT-qPCR analyses to assess their stability under various perturbations, including maturation during extended culture, lactate purification, and various differentiation efficiencies. Expression in hPSC-CMs was also compared with whole human heart tissue. A core set of three novel reference genes (EDF1, DDB1, and ZNF384) exhibited robust stability across the conditions tested, whereas expression of the traditional housekeeping genes tested (ACTB, B2M, GAPDH, and RPL13A) varied significantly under these conditions. Impact statement This article presents an unbiased method for the selection and validation of novel reference genes for real-time quantitative polymerase chain reaction normalization using data from RNA sequencing datasets. This method identified more robust and stable reference genes for gene expression studies during human pluripotent stem cell differentiation to cardiomyocytes than commonly used reference genes. This study also provides a roadmap for identifying reference genes for assessing gene expression during other dynamic cellular processes, including stem cell differentiation to other cell types.