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.

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

The Warner-Bratzler shear force (WBSF) and myofibrillar fragmentation index (MFI) are complementary methodologies used to measure beef tenderness. Longissimus thoracis samples from the 20 most extreme bulls (out of 80 bulls set) for WBSF (tender (n = 10) and tough (n = 10)) and MFI (high (n = 10) and low (n = 10)) traits were collected to perform transcriptomic analysis using RNA-Sequencing. All analysis were performed through CLC Genomics Workbench. A total of 39 and 27 transcripts for WBSF and MFI phenotypes were DE, respectively. The possible DE novel mRNA isoforms, for WBSF and MFI traits, are myosin encoders (e.g. MYL1 and MYL6). In addition, we identified potential mRNA isoforms related to genes affecting the speed fibers degradation during the meat aging process. The DE novel transcripts are transcripted by genes with biological functions related to oxidative process, energy production and striated muscle contraction. The results suggest that the identified mRNA isoforms could be used as potential candidate to select animals in order to improve meat tenderness.

Use of gene expression profile to identify potentially relevant transcripts to myofibrillar fragmentation index trait.

The myofibrillar fragmentation index (MFI) is an indicative trait for meat tenderness. Longissimus thoracis muscle samples from the 20 most extreme bulls (out of 80 bulls set) for MFI (high (n = 10) and low (n = 10) groups) trait were used to perform transcriptomic analysis, using RNA Sequencing (RNA-Seq). An average of 24.616 genes was expressed in the Nellore muscle transcriptome analysis. A total of 96 genes were differentially expressed (p value ≤ 0.001) between the two groups of divergent bulls for MFI. The HEBP2 and BDH1 genes were overexpressed in animals with high MFI. The MYBPH and MYL6, myosin encoders, were identified. The differentially expressed genes were related to increase mitochondria efficiency, especially in cells under oxidative stress conditions, and these also were related to zinc and calcium binding, membrane transport, and muscle constituent proteins, such as actin and myosin. Most of those genes were involved in metabolic pathways of oxidation-reduction, transport of lactate in the plasma membrane, and muscle contraction. This is the first study applying MFI phenotypes in transcriptomic studies to identify and understand differentially expressed genes for beef tenderness. These results suggest that differences detected in gene expression between high and low MFI animals are related to reactive mechanisms and structural components of oxidative fibers under the condition of cellular stress. Some genes may be selected as positional candidate genes to beef tenderness, MYL6, MYBPH, TRIM63, TRIM55, TRIOBP, and CHRNG genes. The use of MFI phenotypes could enhance results of meat tenderness studies.