Rumen and fecal microbiomes are related to diet and production traits in Bos indicus beef cattle.

Background: Ruminants harbor a complex microbial community within their gastrointestinal tract, which plays major roles in their health and physiology. Brazil is one of the largest producers of beef in the world and more than 90% of the beef cattle herds are composed of pure and crossbred Nelore (Bos indicus). Despite its importance to the Brazilian economy and human feeding, few studies have characterized the Nelore microbiome. Therefore, using shotgun metagenomics, we investigated the impact of diet on the composition and functionality of the Nelore microbiome, and explored the associations between specific microbial taxa and their functionality with feed efficiency and methane emission. Results: The ruminal microbiome exhibited significantly higher microbial diversity, distinctive taxonomic profile and variations in microbial functionality compared to the fecal microbiome, highlighting the distinct contributions of the microbiomes of these environments. Animals subjected to different dietary treatments exhibited significant differences in their microbiomes' archaeal diversity and in the abundance of 89 genera, as well as in the functions associated with the metabolism of components of each diet. Moreover, depending on the diet, feed-efficient animals and low methane emitters displayed higher microbial diversity in their fecal microbiome. Multiple genera were associated with an increase or decrease of the phenotypes. Upon analyzing the functions attributed to these taxa, we observed significant differences on the ruminal taxa associated with feed efficient and inefficient cattle. The ruminal taxa that characterized feed efficient cattle stood out for having significantly more functions related to carbohydrate metabolism, such as monosaccharides, di-/oligosaccharides and amino acids. The taxa associated with methane emission had functions associated with methanogenesis and the production of substrates that may influence methane production, such as hydrogen and formate. Conclusion: Our findings highlight the significant role of diet in shaping Nelore microbiomes and how its composition and functionality may affect production traits such as feed efficiency and methane emission. These insights provide valuable support for the implementation of novel feeding and biotechnological strategies.

Prune homolog 2 with BCH domain (PRUNE2) gene expression is associated with feed efficiency-related traits in Nelore steers.

Animal feeding is a critical factor in increasing producer profitability. Improving feed efficiency can help reduce feeding costs and reduce the environmental impact of beef production. Candidate genes previously identified for this trait in differential gene expression studies (e.g., case-control studies) have not examined continuous gene-phenotype variation, which is a limitation. The aim of this study was to investigate the association between the expression of five candidate genes in the liver, measured by quantitative real-time PCR and feed-related traits. We adopted a linear mixed model to associate liver gene expression from 52 Nelore steers with the following production traits: average daily gain (ADG), body weight (BW), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), Kleiber index (KI), metabolic body weight (MBW), residual feed intake (RFI), and relative growth ratio (RGR). The total expression of the prune homolog 2 (PRUNE2) gene was significantly associated with DMI, FCR, FE, and RFI (P < 0.05). Furthermore, we have identified a new transcript of PRUNE2 (TCONS_00027692, GenBank MZ041267) that was inversely correlated with FCR and FE (P < 0.05), in contrast to the originally identified PRUNE2 transcript. The cytochrome P450 subfamily 2B (CYP2B6), early growth response protein 1 (EGR1), collagen type I alpha 1 chain (COL1A1), and connective tissue growth factor (CTGF) genes were not associated with any feed efficiency-related traits (P > 0.05). The findings reported herein suggest that PRUNE2 expression levels affects feed efficiency-related traits variation in Nelore steers.

DNA methylation may affect beef tenderness through signal transduction in Bos indicus.

BACKGROUND: Beef tenderness is a complex trait of economic importance for the beef industry. Understanding the epigenetic mechanisms underlying this trait may help improve the accuracy of breeding programs. However, little is known about epigenetic effects on Bos taurus muscle and their implications in tenderness, and no studies have been conducted in Bos indicus. RESULTS: Comparing methylation profile of Bos indicus skeletal muscle with contrasting beef tenderness at 14 days after slaughter, we identified differentially methylated cytosines and regions associated with this trait. Interestingly, muscle that became tender beef had higher levels of hypermethylation compared to the tough group. Enrichment analysis of predicted target genes suggested that differences in methylation between tender and tough beef may affect signal transduction pathways, among which G protein signaling was a key pathway. In addition, different methylation levels were found associated with expression levels of GNAS, PDE4B, EPCAM and EBF3 genes. The differentially methylated elements correlated with EBF3 and GNAS genes overlapped CpG islands and regulatory elements. GNAS, a complex imprinted gene, has a key role on G protein signaling pathways. Moreover, both G protein signaling pathway and the EBF3 gene regulate muscle homeostasis, relaxation, and muscle cell-specificity. CONCLUSIONS: We present differentially methylated loci that may be of interest to decipher the epigenetic mechanisms affecting tenderness. Supported by the previous knowledge about regulatory elements and gene function, the methylation data suggests EBF3 and GNAS as potential candidate genes and G protein signaling as potential candidate pathway associated with beef tenderness via methylation.

Multi-Omics Approach Reveals miR-SNPs Affecting Muscle Fatty Acids Profile in Nelore Cattle.

MicroRNAs (miRNAs) are key regulators of gene expression, potentially affecting several biological processes, whose function can be altered by sequence variation. Hence, the integration of single nucleotide polymorphisms (SNP) and miRNAs can explain individual differences in economic traits. To provide new insights into the effects of SNPs on miRNAs and their related target genes, we carried out a multi-omic analysis to identify SNPs in miRNA mature sequences (miR-SNPs) associated with fatty acid (FA) composition in the Nelore cattle. As a result, we identified 3 miR-SNPs in different miRNAs (bta-miR-2419-3p, bta-miR-193a-2, and bta-miR-1291) significantly associated with FA traits (p-value < 0.02, Bonferroni corrected). Among these, the rs110817643C>T, located in the seed sequence of the bta-miR-1291, was associated with different ω6 FAs, polyunsaturated FA, and polyunsaturated:saturated FA ratios. Concerning the other two miR-SNPs, the rs43400521T>C (located in the bta-miR-2419-3p) was associated with C12:0 and C18:1 cis-11 FA, whereas the rs516857374A>G (located in the bta-miR-193a-2) was associated with C18:3 ω6 and ratio of ω6/ω3 traits. Additionally, to identify potential biomarkers for FA composition, we described target genes affected by these miR-SNPs at the mRNA or protein level. Our multi-omics analysis outlines the effects of genetic polymorphism on miRNA, and it highlights miR-SNPs and target candidate genes that control beef fatty acid composition.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes.

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

Iron Content Affects Lipogenic Gene Expression in the Muscle of Nelore Beef Cattle.

Iron (Fe) is an essential mineral for metabolism and plays a central role in a range of biochemical processes. Therefore, this study aimed to identify differentially expressed (DE) genes and metabolic pathways in Longissimus dorsi (LD) muscle from cattle with divergent iron content, as well as to investigate the likely role of these DE genes in biological processes underlying beef quality parameters. Samples for RNA extraction for sequencing and iron, copper, manganese, and zinc determination were collected from LD muscles at slaughter. Eight Nelore steers, with extreme genomic estimated breeding values for iron content (Fe-GEBV), were selected from a reference population of 373 animals. From the 49 annotated DE genes (FDR<0.05) found between the two groups, 18 were up-regulated and 31 down-regulated for the animals in the low Fe-GEBV group. The functional enrichment analyses identified several biological processes, such as lipid transport and metabolism, and cell growth. Lipid metabolism was the main pathway observed in the analysis of metabolic and canonical signaling pathways for the genes identified as DE, including the genes FASN, FABP4, and THRSP, which are functional candidates for beef quality, suggesting reduced lipogenic activities with lower iron content. Our results indicate metabolic pathways that are partially influenced by iron, contributing to a better understanding of its participation in skeletal muscle physiology.