Gene co-expression networks associated with carcass traits reveal new pathways for muscle and fat deposition in Nelore cattle

Large-scale transcriptomic and genomic analyses reveal a novel functional gene SERPINB6 for chicken carcass traits

BACKGROUND

Carcass traits are crucial indicators of meat production efficiency. However, the molecular regulatory mechanisms associated with these traits remain unclear.

RESULTS

In this study, we conducted comprehensive transcriptomic and genomic analyses on 399 Tiannong partridge chickens to identify key genes and variants associated with carcass traits and to elucidate the underlying regulatory mechanisms. Based on association analyses with the elastic net (EN) model, we identified 12 candidate genes (AMY1A, AP3B2, CEBPG, EEF2, EIF4EBP1, FGFR1, FOXD3, GOLM1, LOC107052698, PABPC1, SERPINB6 and TBC1D16) for 4 carcass-related traits, namely live weight, dressed weight, eviscerated weight, and breast muscle weight. SERPINB6 was identified as the only overlapping gene by 3 analyses, EN model analysis, weighted gene co-expression network analysis and differential expression analysis. Cell-level experiments confirmed that SERPINB6 promotes the proliferation of chicken DF1 cells and primary myoblasts. Further expression genome-wide association study and association analysis indicated that rs317934171 is the critical site that enhances SERPINB6 expression. Furthermore, a dual-luciferase reporter assay proved that gga-miR-1615 targets the 3'UTR of SERPINB6.

CONCLUSIONS

Collectively, our findings reveal that SERPINB6 serves as a novel gene for chicken carcass traits by promoting fibroblast and myoblast proliferation. Additionally, the downstream variant rs317934171 regulates SERPINB6 expression. These results identify a new target gene and molecular marker for the molecular mechanisms of chicken carcass traits.

A review of the role of transcription factors in regulating adipogenesis and lipogenesis in beef cattle

In the past few decades, genomic selection and other refined strategies have been used to increase the growth rate and lean meat production of beef cattle. Nevertheless, the fast growth rates of cattle breeds are often accompanied by a reduction in intramuscular fat (IMF) deposition, impairing meat quality. Transcription factors play vital roles in regulating adipogenesis and lipogenesis in beef cattle. Meanwhile, understanding the role of transcription factors in regulating adipogenesis and lipogenesis in beef cattle has gained significant attention to increase IMF deposition and meat quality. Therefore, the aim of this paper was to provide a comprehensive summary and valuable insight into the complex role of transcription factors in adipogenesis and lipogenesis in beef cattle. This review summarizes the contemporary studies in transcription factors in adipogenesis and lipogenesis, genome-wide analysis of transcription factors, epigenetic regulation of transcription factors, nutritional regulation of transcription factors, metabolic signalling pathways, functional genomics methods, transcriptomic profiling of adipose tissues, transcription factors and meat quality and comparative genomics with other livestock species. In conclusion, transcription factors play a crucial role in promoting adipocyte development and fatty acid biosynthesis in beef cattle. They control adipose tissue formation and metabolism, thereby improving meat quality and maintaining metabolic balance. Understanding the processes by which these transcription factors regulate adipose tissue deposition and lipid metabolism will simplify the development of marbling or IMF composition in beef cattle.

The Characterization of Subcutaneous Adipose Tissue in Sunit Sheep at Different Growth Stages: A Comprehensive Analysis of the Morphology, Fatty Acid Profile, and Metabolite Profile

Adipose tissue is a crucial economically significant trait that significantly influences the meat quality and growth performance of domestic animals. To reveal the changes in adipose tissue metabolism during the growth of naturally grazing sheep, we evaluated the thickness, adipocyte morphology, fatty acid profile, and metabolite profile of subcutaneous adipose tissue (SAT) from naturally grazing Sunit sheep at 6, 18, and 30 months of age (referred to as Mth-6, Mth-18, and Mth-30, respectively). The fat thickness and adipocyte number were significantly increased with the growth of the sheep (p < 0.05), and the increase of which from Mth-18 to Mth-30 was less than that from Mth-6 to Mth-18. Additionally, the alpha-linolenic acid metabolism was enhanced and fatty acid (FA) elongation increased with growth. The metabolomic analysis revealed 76 differentially expressed metabolites (DEMs) in the SAT in different growth stages. Interestingly, we observed elongation of FAs in lipids correlated with sheep growth. Furthermore, the expression of acylcarnitines was downregulated, and fatty acid amides, aspartic acid, acetic acid and phosphocholine were upregulated in Mth-18 and Mth-30 compared to Mth-6. Altogether, the study found that the difference in SAT in Mth-6 was great compared to Mth-18 and Mth-30. An increase in fat deposition via adipocyte proliferation with the growth of the sheep in naturally grazing. The DEMs of acylcarnitines, fatty acid amides, aspartic acid, acetic acid, and phosphocholine emerged as potential key regulators of adipose tissue metabolism. These findings illustrate the variation in and metabolic mechanism of sheep adipose tissue development under natural grazing, thus providing valuable insights into improving the edible quality of sheep meat and developing the mutton sheep industry.

The Identification of Functional Genes Affecting Fat-Related Meat Traits in Meat-Type Pigeons Using Double-Digest Restriction-Associated DNA Sequencing and Molecular Docking Analysis

The Chinese indigenous Shiqi (SQ) pigeon and the imported White King (WK) pigeon are two meat-type pigeon breeds of economical and nutritional importance in China. They displayed significant differences in such meat quality traits as intramuscular fat (IMF) content and fatty acid (FA) compositions in the breast muscles. In this study, we aimed to screen candidate genes that could affect fat-related meat quality traits in meat-type pigeons. We investigated the polymorphic variations at the genomic level using double-digest restriction-associated DNA (ddRAD) sequencing in 12 squabs of SQ and WK pigeons that exhibited significant inter-breed differences in IMF content as well as FA and amino acid compositions in the breast muscles, and screened candidate genes influencing fat-related traits in squabs through gene ontology analysis and pathway analysis. By focusing on 6019 SNPs, which were located in genes with correct annotations and had the potential to induce changes in the encoded proteins, we identified 19 genes (ACAA1, ACAA2, ACACB, ACADS, ACAT1, ACOX3, ACSBG1, ACSBG2, ACSL1, ACSL4, ELOVL6, FADS1, FADS2, HACD4, HADH, HADHA, HADHB, MECR, OXSM) as candidate genes that could affect fat-related traits in squabs. They were significantly enriched in the pathways of FA metabolism, degradation, and biosynthesis (p < 0.05). Results from molecular docking analysis further revealed that three non-synonymous amino acid alterations, ACAA1(S357N), ACAA2(T234I), and ACACB(H1418N), could alter the non-bonding interactions between the enzymatic proteins and their substrates. Since ACAA1, ACAA2, and ACACB encode rate-limiting enzymes in FA synthesis and degradation, alterations in the enzyme-substrate binding affinity may subsequently affect the catalytic efficiency of enzymes. We suggested that SNPs in these three genes were worthy of further investigation into their roles in explaining the disparities in fat-related traits in squabs.

Gene coexpression network analysis reveals perirenal adipose tissue as an important target of prenatal malnutrition in sheep

We have previously demonstrated that pre- and early postnatal malnutrition in sheep induced depot- and sex-specific changes in adipose morphological features, metabolic outcomes, and transcriptome in adulthood, with perirenal (PER) as the major target followed by subcutaneous (SUB) adipose tissue. We aimed to identify coexpressed and hub genes in SUB and PER to identify the underlying molecular mechanisms contributing to the early nutritional programming of adipose-related phenotypic outcomes. Transcriptomes of SUB and PER of male and female adult sheep with different pre- and early postnatal nutrition histories were used to construct networks of coexpressed genes likely to be functionally associated with pre- and early postnatal nutrition histories and phenotypic traits using weighted gene coexpression network analysis. The modules from PER showed enrichment of cell cycle regulation, gene expression, transmembrane transport, and metabolic processes associated with both sexes' prenatal nutrition. In SUB (only males), a module of enriched adenosine diphosphate metabolism and development correlated with prenatal nutrition. Sex-specific module enrichments were found in PER, such as chromatin modification in the male network but histone modification and mitochondria- and oxidative phosphorylation-related functions in the female network. These sex-specific modules correlated with prenatal nutrition and adipocyte size distribution patterns. Our results point to PER as a primary target of prenatal malnutrition compared to SUB, which played only a minor role. The prenatal programming of gene expression and cell cycle, potentially through epigenetic modifications, might be underlying mechanisms responsible for observed changes in PER expandability and adipocyte-size distribution patterns in adulthood in both sexes.

The Effect of the Feeding System on Fat Deposition in Yak Subcutaneous Fat

Fat deposition is very important to the growth and reproduction of yaks. In this study, the effect of the feeding system on fat deposition in yaks was explored by transcriptomics and lipidomics. The thickness of the subcutaneous fat in yaks under stall (SF) and graze feeding (GF) was evaluated. The transcriptomes and lipidomes of the subcutaneous fat in yaks under different feeding systems were detected by RNA-sequencing (RNA-Seq) and non-targeted lipidomics based on ultrahigh-phase liquid chromatography tandem mass spectrometry (UHPLC-MS), respectively. The differences in lipid metabolism were explored, and the function of differentially expressed genes (DEGs) was evaluated by gene ontology (GO) and Kyoto encyclopedia of genes and genome (KEGG) analysis. Compared with GF yaks, SF yaks possessed stronger fat deposition capacity. The abundance of 12 triglycerides (TGs), 3 phosphatidylethanolamines (PEs), 3 diglycerides (DGs), 2 sphingomyelins (SMs) and 1 phosphatidylcholine (PC) in the subcutaneous fat of SF and GF yaks was significantly different. Under the mediation of the cGMP-PKG signaling pathway, the blood volume of SF and GF yaks may be different, which resulted in the different concentrations of precursors for fat deposition, including non-esterified fatty acid (NEFA), glucose (GLU), TG and cholesterol (CH). The metabolism of C16:0, C16:1, C17:0, C18:0, C18:1, C18:2 and C18:3 in yak subcutaneous fat was mainly realized under the regulation of the INSIG1, ACACA, FASN, ELOVL6 and SCD genes, and TG synthesis was regulated by the AGPAT2 and DGAT2 genes. This study will provide a theoretical basis for yak genetic breeding and healthy feeding.

Bovine HOXA11 Gene Identified from RNA-Seq: mRNA Profile Analysis and Genetic Variation Detection Using ME Method and Their Associations with Carcass Traits

The Homeobox A11 (HOXA11) gene regulates limb skeletal development and muscle growth, thus, it was selected as a candidate gene for bovine carcass traits. In this study, we analyzed the mRNA expression level of HOXA11 in various tissues and cells, and determined the genetic variations in the HOXA11 gene, which might be used as molecular markers for cattle breeding. The mRNA expression profiles of HOXA11 in bovine different tissues showed that HOXA11 was highly expressed in both fat and muscle. The gene expression trend of HOXA11 in myoblasts and adipocytes indicated that HOXA11 might be involved in the differentiation of bovine myoblasts and adipocytes. The data in the Ensembl database showed that there are two putative insertion/deletion (InDel) polymorphisms in the bovine HOXA11 gene. The insertion site (rs515880802) was located in the upstream region (NC_037331.1: g. 68853364-68853365) and named as P1-Ins-4-bp, and the deletion site (rs517582703) was located in the intronic region (NC_037331.1: g. 68859510-68859517) and named as P2-Del-8-bp. These polymorphisms within the HOXA11 gene were identified and genotyped by PCR amplification, agarose gel electrophoresis and DNA sequencing in the 640 Shandong Black Cattle Genetic Resource (SDBCGR) population. Moreover, the mutation frequency was very low after detection, so the mathematical expectation (ME) method was used for detection. Statistical analysis demonstrated that P1-Ins-4-bp was significantly correlated with the beef shoulder (p = 0.012) and tongue root (p = 0.004). Meanwhile, P2-Del-8-bp displayed a significant correlation with the back tendon (p = 0.008), money tendon (p = 2.84 × 10^-4), thick flank (p = 0.034), beef shin (p = 9.09 × 10^-7), triangle thick flank (p = 0.04), triangle flank (p = 1.00 × 10^-6), rump (p = 0.018) and small tenderloin (p = 0.043) in the female SDBCGR population. In summary, these outcomes may provide a new perspective for accelerating the molecular breeding of cattle through marker-assisted selection (MAS) strategies.

Genomic and transcriptomic analyses enable the identification of important genes associated with subcutaneous fat deposition in Holstein cows

Subcutaneous fat deposition has many important roles in dairy cattle, including immunological defense and mechanical protection. The main objectives of this study are to identify key candidate genes regulating subcutaneous fat deposition in high-producing dairy cows by integrating genomic and transcriptomic datasets. A total of 1,654 genotyped Holstein cows are used to perform a genome-wide association study (GWAS) aiming to identify genes associated with subcutaneous fat deposition. Subsequently, weighted gene co-expression network analyses (WGCNA) are conducted based on RNA-sequencing data of 34 cows and de-regressed estimated breeding values of subcutaneous fat deposition. Lastly, differentially expressed (DE) mRNA, lncRNA, and differentially alternative splicing genes are obtained for 12 Holstein cows with extreme and divergent phenotypes for subcutaneous fat deposition. Forty-six protein-coding genes are identified as candidate genes regulating subcutaneous fat deposition in Holstein cattle based on the GWAS. Eleven overlapping genes are identified based on the analyses of DE genes and WGCNA. Furthermore, the candidate genes identified based on the GWAS, WGCNA, and analyses of DE genes are significantly enriched for pathways involved in metabolism, oxidative phosphorylation, thermogenesis, fatty acid degradation, and glycolysis/gluconeogenesis pathways. Integrating all findings, the NID2, STARD3, UFC1, DEDD, PPP1R1B, and USP21 genes are considered to be the most important candidate genes influencing subcutaneous fat deposition traits in Holstein cows. This study provides novel insights into the regulation mechanism underlying fat deposition in high-producing dairy cows, which will be useful when designing management and breeding strategies.

Identification of the hub genes related to adipose tissue metabolism of bovine

Due to the demand for high-quality animal protein, there has been consistent interest in how to obtain more high-quality beef. As well-known, the adipose content of beef has a close connection with the taste and quality of beef, and cattle with different energy or protein diet have corresponding effects on the lipid metabolism of beef. Thus, we performed weighted gene co-expression network analysis (WGCNA) with subcutaneous adipose genes from Norwegian red heifers fed different diets to identify hub genes regulating bovine lipid metabolism. For this purpose, the RNA sequencing data of subcutaneous adipose tissue of 12-month-old Norwegian red heifers (n = 48) with different energy or protein levels were selected from the GEO database, and 7,630 genes with the largest variation were selected for WGCNA analysis. Then, three modules were selected as hub genes candidate modules according to the correlation between modules and phenotypes, including pink, magenta and grey60 modules. GO and KEGG enrichment analysis showed that genes were related to metabolism, and participated in Rap, MAPK, AMPK, VEGF signaling pathways, and so forth. Combined gene interaction network analysis using Cytoscape software, eight hub genes of lipid metabolism were identified, including TIA1, LOC516108, SNAPC4, CPSF2, ZNF574, CLASRP, MED15 and U2AF2. Further, the expression levels of hub genes in the cattle tissue were also measured to verify the results, and we found hub genes in higher expression in muscle and adipose tissue in adult cattle. In summary, we predicted the key genes of lipid metabolism in the subcutaneous adipose tissue that were affected by the intake of various energy diets to find the hub genes that coordinate lipid metabolism, which provide a theoretical basis for regulating beef quality.

Comparative proteomics reveals genetic mechanisms of body weight in Hu sheep and Dorper sheep

Body weight (BW) is a critical economic trait for meat production in sheep, and it is a complex trait affected by numerous elements. The aim of this study was to investigate the genetic mechanisms of sheep BW by a label-free proteomics approach. The result showed, a total of 27, 14, 61, and 65 differentially abundant proteins (DAPs) were identified in the Hu_HBW vs. Hu_LBW, DP_HBW vs. DP_LBW, Hu_HBW vs. DP_HBW, and Hu_LBW vs. DP_LBW comparisons, respectively. Five proteins (including ILK, AHCYL2, MLIP, CYB5A, and SMTNL1) related to fat synthesis and muscle development were detected in the Hu sheep group. In the Dorper sheep group, the screened DAPs strictly related to muscle development and fat synthesis were significantly enriched in MAP kinase activity (MAPK12), Arachidonic acid metabolism, and Steroid hormone biosynthesis (PGFS, LOC101107119) pathways. Several DAPs related to immune responses (SERPINA1, FGG, SERPINC1, and LOC101108131), fat deposition (APOH, GC, AHSG, SKP1, ACSL1, ACAT1, and ACADS), and muscle development (LMOD3 and LRRC39) were detected in the Hu vs. Dorper sheep comparison. These analyses indicated that the BW of sheep is regulated via a variety of pathways, and these DAPs can be further investigated as candidate markers for predicting the BW of sheep. SIGNIFICANCE: Body weight is one of the key traits in sheep and involves multiple coordinated regulatory mechanisms, but the genetic mechanism of BW is still unclear in sheep. In the current study, the label-free method was used to identify the proteins and pathways related to BW using LT muscle of Hu sheep and Dorper sheep with different BW. These findings will provide new candidate proteins and vital pathways into the molecular mechanisms involved growth traits in sheep.

Integrative Analysis Between Genome-Wide Association Study and Expression Quantitative Trait Loci Reveals Bovine Muscle Gene Expression Regulatory Polymorphisms Associated With Intramuscular Fat and Backfat Thickness

Understanding the architecture of gene expression is fundamental to unravel the molecular mechanisms regulating complex traits in bovine, such as intramuscular fat content (IMF) and backfat thickness (BFT). These traits are economically important for the beef industry since they affect carcass and meat quality. Our main goal was to identify gene expression regulatory polymorphisms within genomic regions (QTL) associated with IMF and BFT in Nellore cattle. For that, we used RNA-Seq data from 193 Nellore steers to perform SNP calling analysis. Then, we combined the RNA-Seq SNP and a high-density SNP panel to obtain a new dataset for further genome-wide association analysis (GWAS), totaling 534,928 SNPs. GWAS was performed using the Bayes B model. Twenty-one relevant QTL were associated with our target traits. The expression quantitative trait loci (eQTL) analysis was performed using Matrix eQTL with the complete SNP dataset and 12,991 genes, revealing a total of 71,033 cis and 36,497 trans-eQTL (FDR &lt; 0.05). Intersecting with QTL for IMF, we found 231 eQTL regulating the expression levels of 117 genes. Within those eQTL, three predicted deleterious SNPs were identified. We also identified 109 eQTL associated with BFT and affecting the expression of 54 genes. This study revealed genomic regions and regulatory SNPs associated with fat deposition in Nellore cattle. We highlight the transcription factors FOXP4, FOXO3, ZSCAN2, and EBF4, involved in lipid metabolism-related pathways. These results helped us to improve our knowledge about the genetic architecture behind important traits in cattle.

Transcriptome study digs out BMP2 involved in adipogenesis in sheep tails

Background

Hu sheep and Tibetan sheep in China are characterized by fat tails and thin tails, respectively. Several transcriptomes have been conducted in different sheep breeds to identify the differentially expressed genes (DEGs) underlying this trait. However, these studies identified different DEGs in different sheep breeds.

Results

Hence, RNA sequencing was performed on Hu sheep and Tibetan sheep. We obtained a total of 45.57 and 43.82 million sequencing reads, respectively. Two libraries mapped reads from 36.93 and 38.55 million reads after alignment to the reference sequences. 2108 DEGs were identified, including 1247 downregulated and 861 upregulated DEGs. GO and KEGG analyses of all DEGs demonstrated that pathways were enriched in the regulation of lipolysis in adipocytes and terms related to the chemokine signalling pathway, lysosomes, and glycosaminoglycan degradation. Eight genes were selected for validation by RT–qPCR. In addition, the transfection of BMP2 overexpression into preadipocytes resulted in increased PPAR-γ expression and expression. BMP2 potentially induces adipogenesis through LOX in preadipocytes. The number of lipid drops in BMP2 overexpression detected by oil red O staining was also greater than that in the negative control.

Conclusion

In summary, these results showed that significant genes (BMP2, HOXA11, PPP1CC and LPIN1) are involved in the regulation of adipogenesis metabolism and suggested novel insights into metabolic molecules in sheep fat tails.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08657-8.

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.

Comparative transcriptomic analysis reveals region-specific expression patterns in different beef cuts

BACKGROUND

Beef cuts in different regions of the carcass have different meat quality due to their distinct physiological function. The objective of this study was to characterize the region-specific expression differences using comparative transcriptomics analysis among five representative beef cuts (tenderloin, longissimus lumborum, rump, neck, chuck).

RESULTS

We obtained 15,701 expressed genes in 30 muscle samples across five regions from carcass meat. We identified a total of 80 region-specific genes (RSGs), ranging from three (identified in the rump cut) to thirty (identified in the longissimus lumborum cut), and detected 25 transcription factors (TFs) for RSGs. Using a co-expression network analysis, we detected seven region-specific modules, including three positively correlated modules and four negatively correlated modules. We finally obtained 91 candidate genes related to meat quality, and the functional enrichment analyses showed that these genes were mainly involved in muscle fiber structure (e.g., TNNI1, TNNT1), fatty acids (e.g., SCD, LPL), amino acids (ALDH2, IVD, ACADS), ion channel binding (PHPT1, SNTA1, SUMO1, CNBP), protein processing (e.g., CDC37, GAPDH, NRBP1), as well as energy production and conversion (e.g., ATP8, COX8B, NDUFB6). Moreover, four candidate genes (ALDH2, CANX, IVD, PHPT1) were validated using RT-qPCR analyses which further supported our RNA-seq results.

CONCLUSIONS

Our results provide valuable insights into understanding the transcriptome regulation of meat quality in different beef cuts, and these findings may further help to improve the selection for health-beneficial meat in beef cattle.

Transcriptomic profile of longissimus thoracis associated with fatty acid content in Nellore beef cattle

The beef fatty acid (FA) profile has the potential to impact human health, and displays polygenic and complex features. This study aimed to identify the transcriptomic FA profile in the longissimus thoracis muscle in Nellore beef cattle finished in feedlot. Forty-four young bulls were sampled to assess the beef FA profile by considering 14 phenotypes and including differentially expressed genes (DEG), co-expressed (COE), and differentially co-expressed genes (DCO) analyses. All samples (n = 44) were used for COE analysis, whereas 30 samples with extreme phenotypes for the beef FA profile were used for DEG and DCO. A total of 912 DEG were identified, and the polyunsaturated (n = 563) and unsaturated ω-3 (n = 346) FA sums groups were the most frequently observed. The COE analyses identified three modules, of which the blue module (n = 1776) was correlated with eight of 14 FA phenotypes. Also, 759 DCO genes were listed, and the oleic acid (n = 358) and monounsaturated fatty acids sum (n = 120) were the most frequent. Furthermore, 243 and 13, 319 and seven, and 173 and 12 gene ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways were enriched respectively for the DEG, COE, and DCO analyses. Combining the results, we highlight the unexplored GIPC2, ASB5, and PPP5C genes in cattle. Besides LIPE and INSIG2 genes in COE modules, the ACSL3, ECI1, DECR2, FITM1, and SDHB genes were signaled in at least two analyses. These findings contribute to understand the genetic mechanisms underlying the beef FA profile in Nellore beef cattle finished in feedlot.

Transcriptome Analysis of Bovine Rumen Tissue in Three Developmental Stages

Rumen development is a crucial physiological challenge for ruminants. However, the molecular mechanism regulating rumen development has not been clearly elucidated. In this study, we investigated genes involved in rumen development in 13 rumen tissues from three developmental stages (birth, youth, and adult) using RNA sequencing. We identified that 6,048 genes were differentially expressed among three developmental stages. Using weighted correlation network analysis, we found that 12 modules were significantly associated with developmental stages. Functional annotation and protein–protein interaction (PPI) network analysis revealed that CCNB1, CCNB2, IGF1, IGF2, HMGCL, BDH1, ACAT1, HMGCS2, and CREBBP involved in rumen development. Integrated transcriptome with GWAS information of carcass weight (CW), stomach weight (SW), marbling score (MS), backfat thickness (BFT), ribeye area (REA), and lean meat weight (LMW), we found that upregulated DEGs (fold change 0∼1) in birth–youth comparison were significantly enriched with GWAS signals of MS, downregulated DEGs (fold change >3) were significantly enriched with GWAS signals of SW, and fold change 0∼1 up/downregulated DEGs in birth–adult comparison were significantly enriched with GWAS signals of CW, LMW, REA, and BFT. Furthermore, we found that GWAS signals for CW, LMW, and REA were enriched in turquoise module, and GWAS signals for CW was enriched in lightgreen module. Our study provides novel insights into the molecular mechanism underlying rumen development in cattle and highlights an integrative analysis for illustrating the genetic architecture of beef complex traits.

Identifying Key Genes and Functionally Enriched Pathways of Diverse Adipose Tissue Types in Cattle

Background: Fat is a tissue that not just stores energy and plays a protective role; it is also a vital endocrine organ that generates and integrates signals to influence metabolism. Meanwhile, the excessive accumulation of lipids in adipose tissue can lead to metabolic disturbance and diseases. To date, the complicated molecular mechanisms of bovine adipose tissue are still unknown. This study aimed to identify key genes and functionally enriched pathways in various adipose tissue types.

Results: The RNAseq data of 264 samples were downloaded from Gene Expression Omnibus (GEO) and analyzed by weighted gene co-expression network analysis (WGCNA). We identified 19 modules that significantly associated with at least one adipose tissue type. The brown module from GSE39618 was most closely associated with intramuscular fat tissue, which contained 550 genes. These genes were significantly enriched in pathways that related to inflammation and disease, such as TNF signaling pathway, IL-17 signaling pathway, and NF-kappa B signaling pathway. The pink module (GSE39618) that contained 58 genes was most closely associated with omental fat tissue. The turquoise (GSE39618), blue (GSE116775), and yellow (GSE65125) module were most closely associated with subcutaneous fat tissue. Genes in these modules were significantly enriched in pathways related to fat metabolism, such as the PPAR signaling pathway, fatty acid metabolism and PI3K-Akt signaling pathway. At last, key genes for intramuscular fat (PTGS2 and IL6), omental fat (ARHGEF5 and WT1), and subcutaneous fat (KIT, QR6Q1, PKD2L1, etc.) were obtained and verified. In addition, it was found that IL10 and VCAM1 might be potential genes to distinguish adipose and muscle.

Conclusion: The study applied WGCNA to generate a landscape of adipose tissue and provide a basis for identifying potential pathways and hub genes of different adipose tissue types.

Integrated Analysis of mRNA and MicroRNA Co-expressed Network for the Differentiation of Bovine Skeletal Muscle Cells After Polyphenol Resveratrol Treatment

Resveratrol (RSV) has been confirmed to benefit human health. Resveratrol supplemented in the feeds of animals improved pork, chicken, and duck meat qualities. In this study, we identified differentially expressed (DE) messenger RNAs (mRNAs) (n = 3,856) and microRNAs (miRNAs) (n = 93) for the weighted gene co-expression network analysis (WGCNA) to investigate the co-expressed DE mRNAs and DE miRNAs in the primary bovine myoblasts after RSV treatment. The mRNA results indicated that RSV treatments had high correlations with turquoise module (0.91, P-value = 0.01) and blue module (0.93, P-value < 0.01), while only the turquoise module (0.96, P-value < 0.01) was highly correlated with the treatment status using miRNA data. After biological enrichment analysis, the 2,579 DE genes in the turquoise module were significantly enriched in the Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The top two GO terms were actin filament-based process (GO:0030029) and actin cytoskeleton organization (GO:0030036). The top two KEGG pathways were regulation of actin cytoskeleton (bta04810) and tight junction (bta04530). Then, we constructed the DE mRNA co-expression and DE miRNA co-expression networks in the turquoise module and the mRNA–miRNA targeting networks based on their co-expressions in the key module. In summary, the RSV-induced miRNAs participated in the co-expression networks that could affect mRNA expressions to regulate the primary myoblast differentiation. Our study provided a better understanding of the roles of RSV in inducing miRNA and of the characteristics of DE miRNAs in the key co-expressed module in regulation of mRNAs and revealed new candidate regulatory miRNAs and genes for the beef quality traits.

Transcriptome analysis in comparing carcass and meat quality traits of Jiaxing Black Pig and Duroc × Duroc × Berkshire × Jiaxing Black Pig crosses

This study compares two typical strains: Chinese local excellent meat quality of Jiaxing Black (JXB) Pig and quadratic crossbred pig strain Duroc × Duroc × Berkshire × Jiaxing Black (DDBJ). It was found that between the two pig strains, carcass traits and meat quality traits differed significantly. This is exemplified by the leanness and dressing out percent of DDBJ that were significantly higher than JXB pigs of the same age (P < 0.05) and the better growth rate of DDBJ pigs as to JXB pigs was shown by quantifying muscle proliferation and differentiation of longissimus dorsi muscle employing Hematoxylin and Eosin staining of longissimus dorsi muscle. Nutrients such as inosinic acid, intramuscular fat, and free amino acids in the longissimus dorsi muscle were significantly higher in JXB pigs than DDBJ pigs (p < 0.0001); saturated fatty acids were higher in JXB than in DDBJ pigs (p = 0.0097); essential amino acids and fresh taste amino acids (serine, glutamic acid, proline, glycine, alanine) of JXB pigs was higher than that of DDBJ pigs (p < 0.0001) and amino acids in longissimus dorsi muscle of JXB pigs surpasses the amino acid concentration of DDBJ pigs (p < 0.0001), thus showing the superiority of JXB in terms of meat quality. However, the content of polyunsaturated fatty acids, which is responsible for poor meat quality, was significantly higher in the longissimus dorsi muscle of DDBJ pig than JXB pigs (p < 0.0001); RNA-seq analysis of 5 biological replicates from two of the strains was performed. The screening of 164 up-regulated genes and 183 down-regulated genes found in longissimus dorsi muscle of DDBJ was done and the results identified differentially expressed genes related to muscle development, adipogenesis, amino acid metabolism, fatty acid metabolism and inosine synthesis. In conclusion, the study identified functional genes, elucidated the mechanisms associated with carcass quality traits, meat quality traits and other related traits, and provided means of genetic enhancement to improve meat quality traits and carcass traits in Chinese commercial pigs.

Genetic parameters for traditional and novel ultrasound carcass traits in Nellore cattle

The objective of this study was to contribute to genetic improvement for carcass quality in Nellore cattle by estimating genetic parameters for carcass characteristics measured by traditional and novel ultrasound adjusted for body weight. Records of 19,910 Nellore breed animals born between 2001 and 2015 as part of a genetic improvement program were evaluated for this purpose. The (co)variance components and genetic parameters for yearling weight (YW), slaughter conformation at yearling age (SCY), ribeye area (REA), REA adjusted for 100 kg of body weight (REA100), REA adjusted for 450 kg of body weight (REA450), subcutaneous fat thickness over the loin (SFTL) and the rump (SFTR), intramuscular fat (IF), and height × width ratio of the sirloin steak (RATIO) were estimated using the restricted maximum likelihood (REML) method and the three-trait animal model. The estimated heritability values were 0.39 (YW), 0.33 (SCY), 0.31 (REA), 0.42 (REA100), 0.38 (REA450), 0.31 (RATIO), 0.17 (SFTL), 0.27 (SFTR), and 0.41 (MAR). The genetic correlations between REA and YW (0.43) and SCY and YW (0.78), respectively, suggested that selecting for REA or SCY could result in an increase in body weight. Because of its moderate and negative correlation with YW, REA100 was ineffective as a muscling indicator (- 0.38). The indicators REA450 and RATIO had a - 0.12 and 0.11 low correlation with YW, respectively. Thus, it is concluded that REA450 and RATIO were the best indicators of muscling, with heritability of medium magnitude, making the selection viable, as well as having a low correlation with body weight.

Identification of the Key Genes Associated with the Yak Hair Follicle Cycle

The development of hair follicles in yak shows significant seasonal cycles. In our previous research, transcriptome data including mRNAs and lncRNAs in five stages during the yak hair follicles (HFs) cycle were detected, but their regulation network and the hub genes in different periods are yet to be explored. This study aimed to screen and identify the hub genes during yak HFs cycle by constructing a mRNA-lncRNA co-expression network. A total of 5000 differently expressed mRNA (DEMs) and 729 differently expressed long noncoding RNA (DELs) were used to construct the co-expression network, based on weighted genes co-expression network analysis (WGCNA). Four temporally specific modules were considered to be significantly associated with the HFs cycle of yak. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the modules are enriched into Wnt, EMC-receptor interaction, PI3K-Akt, focal adhesion pathways, and so on. The hub genes, such as FER, ELMO1, PCOLCE, and HOXC13, were screened in different modules. Five hub genes (WNT5A, HOXC13, DLX3, FOXN1, and OVOL1) and part of key lncRNAs were identified for specific expression in skin tissue. Furthermore, immunofluorescence staining and Western blotting results showed that the expression location and abundance of DLX3 and OVOL1 are changed following the process of the HFs cycle, which further demonstrated that these two hub genes may play important roles in HFs development.

Comparative Signatures of Selection Analyses Identify Loci Under Positive Selection in the Murrah Buffalo of India

India is home to a large and diverse buffalo population. The Murrah breed of North India is known for its milk production, and it has been used in breeding programs in several countries. Selection signature analysis yield valuable information about how the natural and artificial selective pressures have shaped the genomic landscape of modern-day livestock species. Genotype information was generated on six buffalo breeds of India, namely, Murrah, Bhadawari, Mehsana, Pandharpuri, Surti, and Toda using ddRAD sequencing protocol. Initially, the genotypes were used to carry out population diversity and structure analysis among the six breeds, followed by pair-wise comparisons of Murrah with the other five breeds through XP-EHH and F _ST methodologies to identify regions under selection in Murrah. Admixture results showed significant levels of Murrah inheritance in all the breeds except Pandharpuri. The selection signature analysis revealed six regions in Murrah, which were identified in more than one pair-wise comparison through both XP-EHH and F _ST analyses. The significant regions overlapped with QTLs for milk production, immunity, and body development traits. Genes present in these regions included SLC37A1, PDE9A, PPBP, CXCL6, RASSF6, AFM, AFP, ALB, ANKRD17, CNTNAP2, GPC5, MYLK3, and GPT2. These genes emerged as candidates for future polymorphism studies of adaptability and performance traits in buffaloes. The results also suggested ddRAD sequencing as a useful cost-effective alternative for whole-genome sequencing to carry out diversity analysis and discover selection signatures in Indian buffalo breeds.

Interference with ACSL1 gene in bovine adipocytes: Transcriptome profiling of circRNA related to unsaturated fatty acid production

Long-chain acyl-CoA synthetase 1 (ACSL1) is a member of the acyl-CoA synthetase family that plays a vital role in lipid metabolism. We have previously shown that the ACSL1 gene regulates the composition of unsaturated fatty acids (UFAs) in bovine skeletal muscle, which in turn regulates the fatty acid synthesis and the generation of lipid droplets. Here, we used RNA-Seq to screen circRNAs that regulated the expression of ACSL1 gene and other UFA synthesis-related genes by RNA interference and noninterference in bovine adipocytes. The results of KEGG pathway analysis showed that the parental genes of differentially expressed (DE)-circRNAs were primarily enriched in the adipocytokine signaling pathway. The prediction results showed that novel_circ_0004855, novel_circ_0001507, novel_circ_0001731, novel_circ_0005276, novel_circ_0002060, novel_circ_0005405 and novel_circ_0004254 regulated UFA synthesis-related genes by interacting with the related miRNAs. These results could help expand our knowledge of the molecular mechanisms of circRNAs in the regulation of UFA synthesis in bovine adipocytes.

Comparative Analyses of Gene Co-expression Networks: Implementations and Applications in the Study of Evolution

Similarities and differences in the associations of biological entities among species can provide us with a better understanding of evolutionary relationships. Often the evolution of new phenotypes results from changes to interactions in pre-existing biological networks and comparing networks across species can identify evidence of conservation or adaptation. Gene co-expression networks (GCNs), constructed from high-throughput gene expression data, can be used to understand evolution and the rise of new phenotypes. The increasing abundance of gene expression data makes GCNs a valuable tool for the study of evolution in non-model organisms. In this paper, we cover motivations for why comparing these networks across species can be valuable for the study of evolution. We also review techniques for comparing GCNs in the context of evolution, including local and global methods of graph alignment. While some protein-protein interaction (PPI) bioinformatic methods can be used to compare co-expression networks, they often disregard highly relevant properties, including the existence of continuous and negative values for edge weights. Also, the lack of comparative datasets in non-model organisms has hindered the study of evolution using PPI networks. We also discuss limitations and challenges associated with cross-species comparison using GCNs, and provide suggestions for utilizing co-expression network alignments as an indispensable tool for evolutionary studies going forward.

Effects of Intensive Fattening With Total Mixed Rations on Carcass Characteristics, Meat Quality, and Meat Chemical Composition of Yak and Mechanism Based on Serum and Transcriptomic Profiles

The objective of this study was to investigate the effects of intensive fattening with total mixed rations (TMR) on carcass characteristics, meat quality, and chemical composition of the yak meat. Theoretical data has been provided for evaluating the quality of yak meat during natural grazing and short-term fattening. Based on the analysis, we found that in fattening yak, the carcass weight (CWT) was increased by 106.43%, whereas the cooking loss, tenderness, and drop loss were significantly improved due to higher intramuscular fat content and lower moisture (P < 0.05). Protein, fat, calcium, and amino acids were also much higher (P < 0.01) in fattening yak compared with the grazing yak. The levels of albumin (ALB), lactate dehydrogenase (LDH), triglyceride (TRIG), and amylase (AMYL) in serum indicated better nutritional status for fattening yaks. The transcriptomics analysis showed that the high expression of ACSL1 and ACACB genes improved the synthesis and deposition of fat in fattening yak, whereas the regulation of SLC7A8, ATP1A4, ATP1A1, SLC3A2, and CPA3 gene expression weakened the proteolysis. These results indicated that fattening with TMR improves the yield and quality of the yak meat.

Selection signatures in two oldest Russian native cattle breeds revealed using high-density single nucleotide polymorphism analysis

Native cattle breeds can carry specific signatures of selection reflecting their adaptation to the local environmental conditions and response to the breeding strategy used. In this study, we comprehensively analysed high-density single nucleotide polymorphism (SNP) genotypes to characterise the population structure and detect the selection signatures in Russian native Yaroslavl and Kholmogor dairy cattle breeds, which have been little influenced by introgression with transboundary breeds. Fifty-six samples of pedigree-recorded purebred animals, originating from different breeding farms and representing different sire lines, of the two studied breeds were genotyped using a genome-wide bovine genotyping array (Bovine HD BeadChip). Three statistical analyses—calculation of fixation index (F_ST) for each SNP for the comparison of the pairs of breeds, hapFLK analysis, and estimation of the runs of homozygosity (ROH) islands shared in more than 50% of animals—were combined for detecting the selection signatures in the genome of the studied cattle breeds. We confirmed nine and six known regions under putative selection in the genomes of Yaroslavl and Kholmogor cattle, respectively; the flanking positions of most of these regions were elucidated. Only two of the selected regions (localised on BTA 14 at 24.4–25.1 Mbp and on BTA 16 at 42.5–43.5 Mb) overlapped in Yaroslavl, Kholmogor and Holstein breeds. In addition, we detected three novel selection sweeps in the genome of Yaroslavl (BTA 4 at 4.74–5.36 Mbp, BTA 15 at 17.80–18.77 Mbp, and BTA 17 at 45.59–45.61 Mbp) and Kholmogor breeds (BTA 12 at 82.40–81.69 Mbp, BTA 15 at 16.04–16.62 Mbp, and BTA 18 at 0.19–1.46 Mbp) by using at least two of the above-mentioned methods. We expanded the list of candidate genes associated with the selected genomic regions and performed their functional annotation. We discussed the possible involvement of the identified candidate genes in artificial selection in connection with the origin and development of the breeds. Our findings on the Yaroslavl and Kholmogor breeds obtained using high-density SNP genotyping and three different statistical methods allowed the detection of novel putative genomic regions and candidate genes that might be under selection. These results might be useful for the sustainable development and conservation of these two oldest Russian native cattle breeds.

Genome-Wide Identification of RNA Editing Sites Affecting Intramuscular Fat in Pigs

Simple Summary

In order to investigate RNA editing sites affecting IMF (which is associated with pork quality and human insulin resistance.), we analyzed the transcriptome and genome sequencing data of a high- and low- groups composed of full-sib pairs pig with opposite IMF phenotypes. Finally, a total of 36 nonredundant RNA editing sites in the longissimus dorsi muscle, which may reveal the potential importance of RNA editing in IMF were identified. Four were selected as candidate sites associated with IMF. Our findings provide some new insights of RNA editing function in pig longissimus dorsi muscle.

Abstract

Intramuscular fat (IMF) is essential for improving the palatability and flavor of meat, and it is strongly associated with human insulin resistance. RNA editing is a widespread regulating event in different tissues. Here, we investigated the global RNA editing difference of two groups of pig with different IMF contents to find the potential editing sites affecting IMF. In this research, RES-Scanner and REDItools were used to identify RNA editing sites based on the whole genome and transcriptome sequencing data of the high and low groups composed of three full-sib pairs with opposite IMF phenotypes. A total of 295 RNA editing sites were investigated in the longissimus dorsi muscle, and 90.17% of these sites caused A to G conversion. After annotation, most editing sites were located in noncoding regions (including five sites located on the 3′ UTR regions). Five editing sites (including two sites that could lead to nonsynonymous amino acid changes) were located in the exons of genes. A total of 36 intergroup (high and low IMF) differential RNA editing sites were found in 33 genes. Some candidate editing sites, such as sites in acyl-coenzymeA: cholesterol acyltransferase 1 (ACAT1), coatomer protein, subunit alpha (COPA), and nuclear receptor coactivator 3 (NCOA3), were selected as candidate RNA editing sites associated with IMF. One site located on the 3′ UTR region of growth hormone secretagogue receptor (GHSR) may regulate GHSR expression by affecting the interaction of miRNA and mRNA. In conclusion, we identified a total of 36 nonredundant RNA editing sites in the longissimus dorsi muscle, which may reveal the potential importance of RNA editing in IMF. Four were selected as candidate sites associated with IMF. Our findings provide some new insights of RNA editing function in pig longissimus dorsi muscle which useful for pig IMF breeding or human insulin resistances research.

The application of omics in ruminant production: a review in the tropical and sub-tropical animal production context

The demand for animal products (e.g. dairy and beef) in tropical regions is expected to increase in parallel with the public demand for sustainable practices, due to factors such as population growth and climate change. The necessity to increase animal production output must be achieved with better management and production technologies. For this to happen, novel research methodologies, animal selection and postgenomic tools play a pivotal role. Indeed, improving breeder selection programs, the quality of meat and dairy products as well as animal health will contribute to higher sustainability and productivity. This would surely benefit regions where resource quality and quantity are increasingly unstable, and research is still very incipient, which is the case of many regions in the tropics. The purpose of this review is to demonstrate how omics-based approaches play a major role in animal science, particularly concerning ruminant production systems and research associated to the tropics and developing countries. SIGNIFICANCE: Environmental conditions in the tropics make livestock production harder, compared to temperate regions. Due to global warming, the sustainability of livestock production will become increasingly problematic. The use of novel omics technologies could generate useful information to understand adaptation mechanisms of resilient breeds and/or species. The application of omics to tropical animal production is still residual in the currently available literature. With this review, we aim to summarize the most notable results in the field whilst encouraging further research to deal with the future challenges that animal production in the tropics will need to face.

Analysis of the genes controlling three quantitative traits in three diverse plant species reveals the molecular basis of quantitative traits

Most traits of agricultural importance are quantitative traits controlled by numerous genes. However, it remains unclear about the molecular mechanisms underpinning quantitative traits. Here, we report the molecular characteristics of the genes controlling three quantitative traits randomly selected from three diverse plant species, including ginsenoside biosynthesis in ginseng (Panax ginseng C.A. Meyer), fiber length in cotton (Gossypium hirsutum L. and G. barbadense L.) and grain yield in maize (Zea mays L.). We found that a vast majority of the genes controlling a quantitative trait were significantly more likely spliced into multiple transcripts while they expressed. Nevertheless, only one to four, but not all, of the transcripts spliced from each of the genes were significantly correlated with the phenotype of the trait. The genes controlling a quantitative trait were multiple times more likely to form a co-expression network than other genes expressed in an organ. The network varied substantially among genotypes of a species and was associated with their phenotypes. These findings indicate that the genes controlling a quantitative trait are more likely pleiotropic and functionally correlated, thus providing new insights into the molecular basis underpinning quantitative traits and knowledge necessary to develop technologies for efficient manipulation of quantitative traits.

Proteome alterations associated with the oleic acid and cis-9, trans-11 conjugated linoleic acid content in bovine skeletal muscle

Oleic acid (OA) and cis-9, trans-11 conjugated linoleic acid (c9t11-CLA) are fatty acids found in beef with beneficial effects in human health. This study investigated differentially abundant proteins (DAPs) in skeletal muscle of bovines with extreme values of OA, and c9t11-CLA. For each one of the fatty acids, twenty muscle samples were divided into two groups (N = 10_High; N = 10_Low) and analyzed by high definition mass spectrometry. We identified 103 and 133 DAPs between the groups for each fatty acid. We found 64 and 45 up-regulated and 39 and 68 down-regulated proteins for OA and c9t11-CLA, respectively. Comparative analysis between proteomic and transcriptomic data revealed eight and ten genes with a consistent between mRNA expression levels and protein abundance for OA and c9t11-CLA, respectively. Unconventional myosin-Id (MYO1D), mineralocorticoid receptor (NR3C2), geranylgeranyl transferase type-2 subunit-alpha (RABGGTA), and uveal autoantigen with coiled-coil domains and ankyrin repeats (UACA) were found as putative candidate proteins for OA content. Fatty acid synthase (FASN), tubulin alpha-4A chain (TUBA4A), vinculin (VCL), NADH dehydrogenase 1 alpha subcomplex 5 (NDUFA5), and prefoldin subunit 6 (PFDN6) for c9t11-CLA. Our findings contribute to a deeper understanding of the molecular mechanisms behind the regulation of the OA and c9t11-CLA content in cattle skeletal muscle. SIGNIFICANCE: Questions about the association between meat intake and disease incidence in humans has driven animal scientist to pursue a better understanding of the biological processes associated with differences in the intramuscular fat composition. The beneficial effects of oleic acid and conjugated linoleic acid in human health have been demonstrated by improving the immune system and preventing atherosclerosis, different types of cancers, hypertension, and diabetes. Previous genome-wide association and gene expression studies identified genomic regions and differentially expressed genes associated with the fatty acid profile in skeletal muscle. In this work, differences were evaluated at the protein level. The use of a label-free quantitative proteomic approach, compared with muscle transcriptome results obtained by RNA-sequencing, allowed us to earn new insights into the variability in fatty acid deposition in skeletal muscle of farm animals. This study opens new avenues to explore the effect of the fatty acids in the skeletal muscle of livestock animals, which is associated with nutritional values of the meat, and perhaps to understand the mechanisms correlated with metabolic diseases in other species.

Genetic Architecture of Carcass and Meat Quality Traits in Montana Tropical® Composite Beef Cattle

The Montana Tropical^® Composite is a recently developed beef cattle population that is rapidly expanding in Brazil and other tropical countries. This is mainly due to its improved meat quality and adaptation to tropical climate conditions compared to Zebu and Taurine cattle breeds, respectively. This study aimed to investigate the genetic architecture of ultrasound-based carcass and meat quality traits in Montana Tropical^® Composite beef cattle. Therefore, we estimated variance components and genetic parameters and performed genome-wide association studies using the weighted single-step Genomic Best Linear Unbiased Prediction (GBLUP) approach. A pedigree dataset containing 28,480 animals was used, in which 1,436 were genotyped using a moderate-density Single Nucleotide Polymorphism panel (30K; 30,105 SNPs). A total of 9,358, 5,768, 7,996, and 1,972 phenotypic records for the traits Longissimus muscle area (LMA), backfat thickness (BFT), rump fat thickness (RFT), and for marbling score (MARB), respectively, were used for the analyses. Moderate to high heritability estimates were obtained and ranged from 0.16 ± 0.03 (RFT) to 0.33 ± 0.05 (MARB). A high genetic correlation was observed between BFT and RFT (0.97 ± 0.02), suggesting that a similar set of genes affects both traits. The most relevant genomic regions associated with LMA, BFT, RFT, and MARB were found on BTA10 (5.4–5.8 Mb), BTA27 (25.2–25.5 Mb), BTA18 (60.6–61.0 Mb), and BTA21 (14.8–15.4 Mb). Two overlapping genomic regions were identified for RFT and MARB (BTA13:47.9–48.1 Mb) and for BFT and RFT (BTA13:61.5–62.3 Mb). Candidate genes identified in this study, including PLAG1, LYN, WWOX, and PLAGL2, were previously reported to be associated with growth, stature, skeletal muscle growth, fat thickness, and fatty acid composition. Our results indicate that ultrasound-based carcass and meat quality traits in the Montana Tropical^® Composite beef cattle are heritable, and therefore, can be improved through selective breeding. In addition, various novel and already known genomic regions related to these traits were identified, which contribute to a better understanding of the underlying genetic background of LMA, BFT, RFT, and MARB in the Montana Tropical Composite population.