Project description:Commercial cuts yield is an important trait for beef production, which affects the final value of the products, but its direct determination is a challenging procedure to be implemented in practice. The measurement of ribeye area (REA) and backfat thickness (BFT) can be used as indirect measures of meat yield. REA and BFT are important traits studied in beef cattle due to their strong implication in technological (carcass yield) and nutritional characteristics of meat products, like the degree of muscularity and total body fat. Thus, the aim of this work was to study the Longissimus dorsi muscle transcriptome of Nellore cattle, associated with REA and BFT, to find differentially expressed (DE) genes, metabolic pathways, and biological processes that may regulate these traits.By comparing the gene expression level between groups with extreme genomic estimated breeding values (GEBV), 101 DE genes for REA and 18 for BFT (false discovery rate, FDR 10%) were identified. Functional enrichment analysis for REA identified two KEGG pathways, MAPK (Mitogen-Activated Protein Kinase) signaling pathway and endocytosis pathway, and three biological processes, response to endoplasmic reticulum stress, cellular protein modification process, and macromolecule modification. The MAPK pathway is responsible for fundamental cellular processes, such as growth, differentiation, and hypertrophy. For BFT, 18 biological processes were found to be altered and grouped into 8 clusters of semantically similar terms. The DE genes identified in the biological processes for BFT were ACHE, SRD5A1, RSAD2 and RSPO3. RSAD2 has been previously shown to be associated with lipid droplet content and lipid biosynthesis.In this study, we identified genes, metabolic pathways, and biological processes, involved in differentiation, proliferation, protein turnover, hypertrophy, as well as adipogenesis and lipid biosynthesis related to REA and BFT. These results enlighten some of the molecular processes involved in muscle and fat deposition, which are economically important carcass traits for beef production.
Project description:The meat market has enormous importance for the world economy, and the quality of the product offered to the consumer is fundamental for the success of the sector. In this study, we analyzed a database which contained information on 2470 animals from a commercial farm in the state of São Paulo, Brazil. Of this total, 2181 animals were genotyped, using 777,962 single-nucleotide polymorphisms (SNPs). After quality control analysis, 468,321 SNPs provided information on the number of genotyped animals. Genome-wide association analyses (GWAS) were performed for the characteristics of the rib eye area (REA), subcutaneous fat thickness (SFT), shear force at 7 days' ageing (SF7), and intramuscular fat (IMF), with the aid of the single-step genomic best linear unbiased prediction (ssGBLUP) method, with the purpose of identifying possible genomic windows (~1 Mb) responsible for explaining at least 0.5% of the genetic variance of the traits under analysis (≥0.5%). These genomic regions were used in a gene search and enrichment analyses using MeSH terms. The distributed heritability coefficients were 0.14, 0.20, 0.18, and 0.21 for REA, SFT, SF7, and IMF, respectively. The GWAS results indicated significant genomic windows for the traits of interest in a total of 17 chromosomes. Enrichment analyses showed the following significant terms (FDR ≤ 0.05) associated with the characteristics under study: for the REA, heat stress disorders and life cycle stages; for SFT, insulin and nonesterified fatty acids; for SF7, apoptosis and heat shock proteins (HSP27); and for IMF, metalloproteinase 2. In addition, KEGG (Kyoto encyclopedia of genes and genomes) enrichment analysis allowed us to highlight important metabolic pathways related to the studied phenotypes, such as the growth hormone synthesis, insulin-signaling, fatty acid metabolism, and ABC transporter pathways. The results obtained provide a better understanding of the molecular processes involved in the expression of the studied characteristics and may contribute to the design of selection strategies and future studies aimed at improving the productivity of Nellore cattle.
Project description:Two complementary methods, namely Multi-Trait Meta-Analysis and Versatile Gene-Based Test for Genome-wide Association Studies (VEGAS), were used to identify putative pleiotropic genes affecting carcass traits in Bos indicus (Nellore) cattle. The genotypic data comprised over 777,000 single-nucleotide polymorphism markers scored in 995 bulls, and the phenotypic data included deregressed breeding values (dEBV) for weight measurements at birth, weaning and yearling, as well visual scores taken at weaning and yearling for carcass finishing precocity, conformation and muscling. Both analyses pointed to the pleomorphic adenoma gene 1 (PLAG1) as a major pleiotropic gene. VEGAS analysis revealed 224 additional candidates. From these, 57 participated, together with PLAG1, in a network involved in the modulation of the function and expression of IGF1 (insulin like growth factor 1), IGF2 (insulin like growth factor 2), GH1 (growth hormone 1), IGF1R (insulin like growth factor 1 receptor) and GHR (growth hormone receptor), suggesting that those pleiotropic genes operate as satellite regulators of the growth pathway.
Project description:This study evaluated phenotypic (rph) and genetic correlations (rg) between 8 feed efficiency traits and other traits of economic interest including weight at selection (WS), loin-eye area (LEA), backfat thickness (BF), and rump fat thickness (RF) in Nellore cattle. Feed efficiency traits were gain:feed, residual feed intake (RFI), residual feed intake adjusted for backfat thickness (RFIb) and for backfat and rump fat thickness (RFIsf), residual body weight gain (RG), residual intake and body weight gain (RIG), and residual intake and body weight gain using RFIb (RIGb) and RFIsf (RIGsf). The variance components were estimated by the restricted maximum likelihood method using a two-trait animal model. The heritability estimates (h2) were 0.14, 0.24, 0.20, 0.22, 0.19, 0.15, 0.11 and 0.11 for gain:feed, RFI, RFIb, RFIsf, RG, RIG, RIGb and RIGsf, respectively. All rph values between traits were close to zero, except for the correlation of feed efficiency traits with dry matter intake and average daily gain. High rg values were observed for the correlation of dry matter intake, average daily gain and metabolic weight with WS and hip height (>0.61) and low to medium values (0.15 to 0.48) with the carcass traits (LEA, BF, RF). Among the feed efficiency traits, RG showed the highest rg with WS and hip height (0.34 and 0.25) and the lowest rg with subcutaneous fat thickness (-0.17 to 0.18). The rg values of RFI, RFIb and RFIsf with WS (0.17, 0.23 and 0.22), BF (0.37, 0.33 and 0.33) and RF (0.30, 0.31 and 0.32) were unfavorable. The rg values of gain:feed, RIG, RIGb and RIGsf with WS were low and favorable (0.07 to 0.22), while medium and unfavorable (-0.22 to -0.45) correlations were observed with fat thickness. The inclusion of subcutaneous fat thickness in the models used to calculate RFI did not reduce the rg between these traits. Selecting animals for higher feed efficiency will result in little or no genetic change in growth and will decrease subcutaneous fat thickness in the carcass.
Project description:Data integration using hierarchical analysis based on the central dogma or common pathway enrichment analysis may not reveal non-obvious relationships among omic data. Here, we applied factor analysis (FA) and Bayesian network (BN) modeling to integrate different omic data and complex traits by latent variables (production, carcass, and meat quality traits). A total of 14 latent variables were identified: five for phenotype, three for miRNA, four for protein, and two for mRNA data. Pearson correlation coefficients showed negative correlations between latent variables miRNA 1 (mirna1) and miRNA 2 (mirna2) (-0.47), ribeye area (REA) and protein 4 (prot4) (-0.33), REA and protein 2 (prot2) (-0.3), carcass and prot4 (-0.31), carcass and prot2 (-0.28), and backfat thickness (BFT) and miRNA 3 (mirna3) (-0.25). Positive correlations were observed among the four protein factors (0.45-0.83): between meat quality and fat content (0.71), fat content and carcass (0.74), fat content and REA (0.76), and REA and carcass (0.99). BN presented arcs from the carcass, meat quality, prot2, and prot4 latent variables to REA; from meat quality, REA, mirna2, and gene expression mRNA1 to fat content; from protein 1 (prot1) and mirna2 to protein 5 (prot5); and from prot5 and carcass to prot2. The relations of protein latent variables suggest new hypotheses about the impact of these proteins on REA. The network also showed relationships among miRNAs and nebulin proteins. REA seems to be the central node in the network, influencing carcass, prot2, prot4, mRNA1, and meat quality, suggesting that REA is a good indicator of meat quality. The connection among miRNA latent variables, BFT, and fat content relates to the influence of miRNAs on lipid metabolism. The relationship between mirna1 and prot5 composed of isoforms of nebulin needs further investigation. The FA identified latent variables, decreasing the dimensionality and complexity of the data. The BN was capable of generating interrelationships among latent variables from different types of data, allowing the integration of omics and complex traits and identifying conditional independencies. Our framework based on FA and BN is capable of generating new hypotheses for molecular research, by integrating different types of data and exploring non-obvious relationships.
Project description:Sex condition has been demonstrated to alter meat quality and sex is a major factor that affects the fatty acid composition of lipids of carcass dissectible or intramuscular depot fats. But the possible genetic molecular mechanism of gender causing meat quality differences is not well defined. Qinchuan cattle, Qinghai yak and Guangxi buffalo are three typical indigenous species of cattle in China. Obivious differences of meat quality exist among the three species of cattle. Few studies have been conducted to elucidate the muscle tissue expression of genes involved in pathways and mechanisms leading to meat quality differences beyond the phenotype properties of beef. Bovine Genome Arrays were used to construct muscle expression profiles of the longuissimus dorsi from Qinchuan cattle at 36 months and screen differentially expressed genes in the longuissimus dorsi muscle tissues among different genders of Qinchuan cattle, between Qinchuan cattle and Qinghai yak, and between Qinchuan cattle and Guangxi buffalo.
Project description:MicroRNAs (miRNAs) regulate gene expression. The aim of this study was to identify circulating miRNAs that are involved with meat yield and connect exosomes and longissimus dorsi (LD) muscle in Korean cattle steer. Thus, we performed analyses of the carcass characteristics, exosome extraction, bovine miRNA array, mature miRNA qPCR, bioinformatics, and qPCR. Our analysis of the carcass characteristics relative to the yield grade (YG) showed that the yield index (YI) and rib eye area were the highest, whereas the backfat thickness was the lowest for YG A cattle among the three YGs. We performed a miRNA microarray and a Venn diagram analysis to sort the circulating miRNAs that connect exosomes and LD muscle. Mature miRNA qPCR showed that miR-15a (r = 0.84), miR-26b (r = 0.91), and miR-29c (r = 0.92) had positive relationships with exosomes and LD muscle. In YG A cattle, miR-26b was considered to be a circulating miRNA connecting exosomes and LD muscle because the expression pattern of miR-26b was similar in the microarray and miRNA qPCR for miR-26b-targeted genes, where the expression levels of DIAPH3 and YOD1 were lower than those in YG C cattle. Our results suggest that circulating miR-26b may participate in cell maintenance by regulating DIAPH3 and YOD1 in the LD muscle of Korean cattle. 29 or 30 month Korean cattle steers were classified by yield grade(YG) A, B, and C. miRNA expression patterns for exosome and LD muscle were compared among YGs.
Project description:Cattle plays an important role in providing essential nutrients through meat production. Thus, we focused on epigenetic factors associated with meat yield. To investigating circulating miRNAs that are involved with meat yield and connect biofluids and longissimus dorsi (LD) muscle in Korean cattle, we performed analyses of the carcass characteristics, miRNA array, qPCR, and bioinformatics. Carcass characteristics relative to the yield grade (YG) showed that the yield index and rib eye area were the highest, whereas the backfat thickness was the lowest for YG A (equal to high yield grade) cattle among the three YGs. miRNA array sorted the circulating miRNAs that connect biofluids and LD muscle. miRNA qPCR showed that miR-15a (r = 0.84), miR-26b (r = 0.91), and miR-29c (r = 0.92) had positive relationships with biofluids and LD muscle. In YG A cattle, miR-26b was considered to be a circulating miRNA connecting biofluids and LD muscle because the target genes of miR-26b was more involved with myogenesis. Then, miR-26b targeted genes, DIAPH3 and YOD1 were downregulated in YG A cattle. Our results suggest that miR-15a, miR-26b, and miR-29c are upregulated in biofluids and LD muscle whereas, downregulation of DIAPH3 and YOD1 in the LD muscle of finishing cattle steers.
Project description:Heat stress has a detrimental impact on cattle health, welfare and productivity by affecting gene expression, metabolism and immune response, but little is known on the epigenetic mechanisms mediating the effect of temperatureat the cellular and organism level. In this study, we investigated genome-wide DNA methylation in blood samples collected from 5 bulls of the heat stress resilient Nellore breed and 5 bulls of the Angus that are more heat stress susceptible, exposed to the sun and high temperature-high humidity during the summer season of the Brazilian South-East region. The methylomes were analyzed during and after the exposure by Reduced Representation Bisulfite Sequencing, which provided genome-wide single-base resolution methylation profiles. Significant methylation changes between stressful and recovery periods were observed in 822 genes. Among these, 352 were only seen in Angus, 367 were specific to Nellore, and 103 showed significant changes in methylation patterns in both breeds. KEGG and Gene Ontology (GO) enrichment analyses showed that responses were breed-specific. Interestingly, in Nellore significant genes and pathways were mainly involved in stress responses and cellular defense and were undermethylated during heat stress. Whereas in Angus the response was less focused. These preliminary results suggest that heat challenge induces changes in methylation patterns in specific loci, which should be further scrutinized to assess their role in heat tolerance.