Bioinformatics analysis of transcriptome dynamics during growth in angus cattle longissimus muscle

Antioxidant Activity of Different Tissues from Rabbits Fed Dietary Bovine Colostrum Supplementation

Recent advances in animal nutrition have indicated that bovine colostrum (BC), due to its content of macronutrients, micronutrients and bioactive compounds, is an excellent health supplement. To the best of our knowledge, no studies on the effect of BC on antioxidant status have been performed in rabbits. This study aimed to investigate the effect of two BC concentrations on antioxidant status and gene expression of antioxidant enzymes in some tissues of rabbits. Thirty New Zealand White male rabbits were randomly divided into three experimental diets, containing 0% (CON), 2.5%, and 5% of BC (BC-2.5 and BC-5, respectively). The activity of antioxidant enzymes in plasma (catalase: CAT; glutathione peroxidase: GPx; superoxide dismutase: SOD), and the enzymes' gene expression in the liver and longissimus dorsi muscle, were determined. Results showed no significant differences, neither in plasma nor in tissues. A significant tissue-related effect has been observed regarding the mRNA levels of SOD and GPx, which were higher in the LD (p = 0.022) and liver (p = 0.001), respectively. Further studies, considering modifications of the length and dosage of dietary BC supplementation, are required to update the current state of knowledge in rabbits, as well as to fully understand the potential value of BC for possible application in farming use.

Comparative transcriptome profiles of four sexually size dimorphic fish

Sexual size dimorphism is widespread in fish species. Although sex growth differences in multiple species have been studied successively, the commonalities of regulatory mechanisms across sexually dimorphic species are unknown. In this study, we performed RNA-seq analysis of four representative fish (loach, half-smooth tongue sole, yellow catfish, and Nile tilapia) with significant growth differences between females and males. Clean reads were identified from four fish species, ranging from 45,718,052 to 57,733,120. Following comparison transcriptome analysis, there were 1,132 and 1,108, 1,290 and 1,102, 4,732 and 4,266, 748 and 192 differentially expressed genes (DEGs) in the brain and muscle of loach, half-smooth tongue sole, yellow catfish, and Nile tilapia, respectively. Furthermore, the expression levels were validated by quantitative real-time PCR (qRT-PCR). Comparative transcriptome profiles of four fish described here will provide fundamental information for further studies on the commonalities of sexually size dimorphic fish in regulating growth differences between females and males.

Embryonic transcriptome unravels mechanisms and pathways underlying embryonic development with respect to muscle growth, egg production, and plumage formation in native and broiler chickens

Background: Muscle development, egg production, and plumage colors are different between native and broiler chickens. The study was designed to investigate why improved Aseel (PD4) is colorful, stronger, and grew slowly compared with the control broiler (CB).

Methods: A microarray was conducted using the 7th-day embryo (7EB) and 18th-day thigh muscle (18TM) of improved Aseel and broiler, respectively. Also, we have selected 24 Gallus gallus candidate reference genes from NCBI, and total RNA was isolated from the broiler, improved Aseel embryo tissues, and their expression profiles were studied by real-time quantitative PCR (qPCR). Furthermore, microarray data were validated with qPCR using improved Aseel and broiler embryo tissues.

Results: In the differential transcripts screening, all the transcripts obtained by microarray of slow and fast growth groups were screened by fold change ≥ 1 and false discovery rate (FDR) ≤ 0.05. In total, 8,069 transcripts were differentially expressed between the 7EB and 18TM of PD4 compared to the CB. A further analysis showed that a high number of transcripts are differentially regulated in the 7EB of PD4 (6,896) and fewer transcripts are differentially regulated (1,173) in the 18TM of PD4 compared to the CB. On the 7th^- and 18th-day PD4 embryos, 3,890, 3,006, 745, and 428 transcripts were up- and downregulated, respectively. The commonly up- and downregulated transcripts are 91 and 44 between the 7th- and 18th-day of embryos. In addition, the best housekeeping gene was identified. Furthermore, we validated the differentially expressed genes (DEGs) related to muscle growth, myostatin signaling and development, and fatty acid metabolism genes in PD4 and CB embryo tissues by qPCR, and the results correlated with microarray expression data.

Conclusion: Our study identified DEGs that regulate the myostatin signaling and differentiation pathway; glycolysis and gluconeogenesis; fatty acid metabolism; Jak-STAT, mTOR, and TGF-β signaling pathways; tryptophan metabolism; and PI3K-Akt signaling pathways in PD4. The results revealed that the gene expression architecture is present in the improved Aseel exhibiting embryo growth that will help improve muscle development, differentiation, egg production, protein synthesis, and plumage formation in PD4 native chickens. Our findings may be used as a model for improving the growth in Aseel as well as optimizing the growth in the broiler.

Identification of differentially expressed genes at different post-natal development stages of longissimus dorsi muscle in Tianzhu white yak

The regulatory mechanisms controlling post-natal muscle development in the yak (Bos grunniens) are still largely unknown, yet the growth and development of muscle is a complex process that plays a crucial role in determining the yield and quality of an animal's meat. In this study, we performed a transcriptome analysis based on the RNA sequencing (RNA-Seq) of yak longissimus dorsi muscle tissue obtained from calves (6 months of age; 6 M), young adults (30 months of age; 30 M) and adult (54 months of age; 54 M) to identify which genes are differentially expressed and to investigate their temporal expression profiles. In total, 1788 differentially expressed genes (DEGs) (|log2FC| ≥ 1, P-adjusted < 0.05) were detected by pairwise comparisons between the different age groups. The expression levels of 10 of the DEGs were confirmed using reverse transcription-quantitative PCR (RT-qPCR), and the results were consistent with the transcriptome profile. A time-series expression profile analysis clustered the DEGs into four groups that could be divided into two classes (P < 0.05): class 1 profiles, which had up-regulated patterns of gene expression and class 2 profiles, which featured down-regulated patterns. Based on that cluster analysis, GO enrichment analysis revealed 1073, 127, and 184 terms as significantly enriched in biological process (BP), cellular component (CC), and molecular function (MF) categories in the class 1 profiles, while 714, 66, and 206 terms were significantly enriched in BP, CC, and MF in the class 2 profiles. A KEGG pathway analysis revealed that DEGs from the class 1 profiles were enriched in 62 pathways, with the most enriched being the phosphoinositide 3-kinase (PI3K) - protein kinase B (Akt)-signaling pathway. The DEGs from the class 2 profiles were enriched in 16 pathways, of which forkhead box protein O (FoxO) - signaling was the most enriched. Taken together, these results provide insight into the mechanisms of skeletal muscle development, as well suggesting some potential genes of importance for yak meat production.

Profiling and Functional Analysis of Long Noncoding RNAs and mRNAs during Porcine Skeletal Muscle Development

Gene transcripts or mRNAs and long noncoding RNAs (lncRNAs) are differentially expressed during porcine skeletal muscle development. However, only a few studies have been conducted on skeletal muscle transcriptome in pigs based on timepoints according to the growth curve for porcine. Here, we investigated gene expression in Qingyu pigs at three different growth stages: the inflection point with the maximum growth rate (MGI), the inflection point of the gradually increasing stage to the rapidly increasing stage (GRI), and the inflection point of the rapidly increasing stage to the slowly increasing stage (RSI). Subsequently, we explored gene expression profiles during muscle development at the MGI, GRI and RSI stages by Ribo-Zero RNA sequencing. Qingyu pigs reached the MGI, GRI and RSI stages at 156.40, 23.82 and 288.97 days of age with 51.73, 3.14 and 107.03 kg body weight, respectively. A total of 14,530 mRNAs and 11,970 lncRNAs were identified at the three stages, and 645, 323 differentially expressed genes (DEGs) and 696, 760 differentially expressed lncRNAs (DELs) were identified in the GRI vs. MGI, and RSI vs. MGI, comparisons. Functional enrichment analysis revealed that genes involved in immune system development and energy metabolism (mainly relate to amino acid, carbohydrate and lipid) were enriched at the GRI and MGI stages, respectively, whereas genes involved in lipid metabolism were enriched at the RSI stage. We further characterized G1430, an abundant lncRNA. The full-length sequence (316 nt) of lncRNA G1430 was determined by rapid amplification of cDNA ends (RACE). Subcellular distribution analysis by quantitative real-time PCR (qRT-PCR) revealed that G1430 is a cytoplasmic lncRNA. Binding site prediction and dual luciferase assay showed that lncRNA G1430 directly binds to microRNA 133a (miR-133a). Our findings provide the basis for further investigation of the regulatory mechanisms and molecular genetics of muscle development in pigs.

Molecular signatures of muscle growth and composition deciphered by the meta-analysis of age-related public transcriptomics data

The lean-to-fat ratio is a major issue in the beef meat industry from both carcass and meat production perspectives. This industrial perspective has motivated meat physiologists to use transcriptomics technologies to decipher mechanisms behind fat deposition within muscle during the time course of muscle growth. However, synthetic biological information from this volume of data remains to be produced to identify mechanisms found in various breeds and rearing practices. We conducted a meta-analysis on 10 transcriptomic data sets stored in public databases, from the longissimus thoracis of five different bovine breeds divergent by age. We updated gene identifiers on the last version of the bovine genome (UCD1.2), and the 715 genes common to the 10 studies were subjected to the meta-analysis. Of the 238 genes differentially expressed (DEG), we identified a transcriptional signature of the dynamic regulation of glycolytic and oxidative metabolisms that agrees with a known shift between those two pathways from the animal puberty. We proposed some master genes of the myogenesis, namely MYOG and MAPK14, as probable regulators of the glycolytic and oxidative metabolisms. We also identified overexpressed genes related to lipid metabolism (APOE, LDLR, MXRA8, and HSP90AA1) that may contribute to the expected enhanced marbling as age increases. Lastly, we proposed a transcriptional signature related to the induction (YBX1) or repression (MAPK14, YWAH, ERBB2) of the commitment of myogenic progenitors into the adipogenic lineage. The relationships between the abundance of the identified mRNA and marbling values remain to be analyzed in a marbling biomarkers discovery perspectives.

Comparative transcriptomic analysis of hypothalamus-pituitary-liver axis in bighead carp (Hypophthalmichthys nobilis) with differential growth rate

Background

Growth rate is one of the most important features for aquaculture species and deciphering its regulation mechanism has great significance both in genetics and in economics. Hypothalamus-pituitary growth axis (HP growth axis) or neuro-endocrine axis plays a vital role in growth regulation in different aquaculture animals.

Results

In this study, the HP and liver transcriptomes of two female groups (H and L) with phenotypically extreme growth rate were sequenced using RNA-Seq. A total of 30,524 and 22,341 genes were found expressed in the two tissues, respectively. The average expression levels for the two tissues were almost the same, but the median differed significantly. A differential expression analysis between H and L groups identified 173 and 204 differentially expressed genes (DEGs) in HP and liver tissue, respectively. Pathway analysis revealed that DEGs in HP tissue were enriched in regulation of cell proliferation and angiogenesis while in liver tissue these genes were overrepresented in sterol biosynthesis and transportation. Genomic overlapping analyses found that 4 and 5 DEGs were within growth-related QTL in HP and liver tissue respectively. A deeper analysis of these 9 genes revealed 3 genes were functionally linked to the trait of interest. The expression of 2075 lncRNAs in HP tissue and 1490 in liver tissue were also detected, and some of lncRNAs were highly expressed in the two tissues.

Conclusions

Above all, the results of the present study greatly contributed to the knowledge of the regulation of growth and then assisted the design of new selection strategies for bighead carp with improved growth-related traits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5691-4) contains supplementary material, which is available to authorized users.

Glutathione content and expression of proteins involved with glutathione metabolism differs in longissimus dorsi, subcutaneous adipose, and liver tissues of finished vs. growing beef steers

Increased tissue redox state may result in sub-optimal growth. Our goal was to determine if glutathione (GSH) content and expression of proteins involved with GSH metabolism change in longissimus dorsi (LD), subcutaneous adipose (SA), and liver tissues of growing vs. finishing steer phenotypes. Tissues were taken from 16 Angus steers (BW = 209 ± 29.4 kg) randomly assigned (n = 8) to develop through Growing (final BW = 301 ± 7.06 kg) vs. Finished (final BW = 576 ± 36.9 kg) growth phases, and at slaughter had achieved different rib-eye area (REA) (53.2, 76.8 cm2), marbling scores (296, 668), and 12th rib adipose thickness (0.54, 1.73 cm), respectively (Amino Acids, doi:10.1007/s00726-018-2540-8). GSH content (mg/g wet tissue) was determined by a commercial assay and the relative content of target protein and mRNA in tissue homogenates was determined by Western blot and reverse-transcribed PCR analyses, respectively. The effect of growth phase (Finished vs. Growing) was assessed by ANOVA using the GLM procedure of SAS. The LD of Finished steers had more (P < 0.04) GSH (42%) and GSH synthesizing (GCLC, 61%; GCLM, 21%) and metabolizing (GPX1, 42%; GPX3, 73%; GGT1, 56%) enzymes, and less (P < 0.02) GPX2 (46%), EAAC1 (30%) and glutamine synthetase (GS) (28%), whereas GTRAP3-18 and ARL6IP1 did not differ (P > 0.57). Principal component analysis found that GSH content of LD was associated with REA and marbling score. The SA of Finished steers had less (P < 0.04) GSH (38%), GSH metabolizing (GPX4, 52%; GGT1, 71%) enzyme mRNA, and GTRAP3-18 (123%) and ARL6IP1 (43%), whereas the mRNA content of GSH-synthesizing enzymes and content of EAAC1 and GS did not differ (P > 0.32). The liver of Finished steers had less (P < 0.02) mRNA content of GSH synthesizing (GCLC, 39%; GSS 29%) and metabolizing (GPX1, 30%) enzymes, and more (P < 0.01) GSTM1 metabolizing enzyme (114%). The change in GSH content as steers fattened indicate an increased antioxidant capacity in the LD of Finished steers, and a decreased antioxidant capacity in SA, consistent with changes in enzyme and transporter expression. Changes in liver enzyme and transporter expression were consistent with no change in GSH content. The relationship of EAAC1 regulatory proteins (GTRAP3-18, ARL6IP1) to GSH, EAAC1, and GS content differs and changes as Growing steers develop into Finished phenotypes. These findings provide mechanistic insight into how antioxidant capacity occurs in tissues of economic and metabolic importance as cattle fatten.

Skeletal muscle and liver gene expression profiles in finishing steers supplemented with Amaize

Our main objective was to evaluate the effects of feeding α-amylase (Amaize, Alltech Inc., Nicholasville, KY, USA) for 140 days on skeletal muscle and liver gene transcription in beef steers. Steers fed Amaize had lower average daily gain (p = .03) and gain:feed ratio (p = .05). No differences (p > .10) in serum metabolites or carcass traits were detected between the two groups but Amaize steers tended (p < .15) to have increased 12th rib fat depth. Microarray analysis of skeletal muscle revealed 21 differentially expressed genes (DEG), where 14 were up-regulated and seven were down-regulated in Amaize-fed steers. The bioinformatics analysis indicated that metabolic pathways involved in fat formation and deposition, stress response, and muscle function were activated, while myogenesis was inhibited in Amaize-fed steers. The quantitative PCR results for liver revealed a decrease (p < .01) in expression of fatty acid binding protein 1 (FABP1) and 3-hydroxybutyrate dehydrogenase 1 (BDH1) with Amaize. Because these genes are key for intracellular fatty acid transport, oxidation and ketone body production, data suggest a reduction in hepatic lipid catabolism. Future work to investigate potential positive effects of Amaize on cellular stress response, muscle function, and liver function in beef cattle appears warranted.

Transcriptional changes in mesenteric and subcutaneous adipose tissue from Holstein cows in response to plane of dietary energy

Background

Dairy cows can readily overconsume dietary energy during most of the prepartum period, often leading to higher prepartal concentrations of insulin and glucose and excessive body fat deposition. The end result of these physiologic changes is greater adipose tissue lipolysis post-partum coupled with excessive hepatic lipid accumulation and compromised health. Although transcriptional regulation of the adipose response to energy availability is well established in non-ruminants, such regulation in cow adipose tissue depots remains poorly characterized.

Results

Effects of ad-libitum access to high [HIGH; 1.62 Mcal/kg of dry matter (DM)] or adequate (CON; 1.35 Mcal/kg of DM) dietary energy for 8 wk on mesenteric (MAT) and subcutaneous (SAT) adipose tissue transcript profiles were assessed in non-pregnant non-lactating Holstein dairy cows using a 13,000-sequence annotated bovine oligonucleotide microarray. Statistical analysis revealed 409 and 310 differentially expressed genes (DEG) due to tissue and diet. Bioinformatics analysis was conducted using the Dynamic Impact Approach (DIA) with the KEGG pathway database. Compared with SAT, MAT had more active biological processes related to adipose tissue accumulation (adiponectin secretion) and signs of pro-inflammatory processes due to adipose tissue expansion and macrophage infiltration (generation of ceramides). Feeding the HIGH diet led to changes in mRNA expression of genes associated with cell hypertrophy (regucalcin), activation of adipogenesis (phospholipid phosphatase 1), insulin signaling activation (neuraminidase 1) and angiogenesis (semaphorin 4G, plexin B1). Further, inflammation due to HIGH was underscored by mRNA expression changes associated with oxidative stress response (coenzyme Q3, methyltransferase), ceramide synthesis (N-acylsphingosine amidohydrolase 1), and insulin signaling (interferon regulatory factor 1, phosphoinositide-3-kinase regulatory subunit 1, retinoic acid receptor alpha). Activation of ribosome in cows fed HIGH indicated the existence of greater adipocyte growth rate (M-phase phosphoprotein 10, NMD3 ribosome export adaptor).

Conclusions

The data indicate that long-term ad-libitum access to a higher-energy diet led to transcriptional changes in adipose tissue that stimulated hypertrophy and the activity of pathways associated with a slight but chronic inflammatory response. Further studies would be helpful in determining the extent to which mRNA results also occur at the protein level.

Transcriptome-Wide Analysis Reveals the Role of PPARγ Controlling the Lipid Metabolism in Goat Mammary Epithelial Cells

To explore the large-scale effect of peroxisome proliferator-activated receptor γ (PPARG) in goat mammary epithelial cells (GMEC), an oligonucleotide microarray platform was used for transcriptome profiling in cells overexpressing PPARG and incubated with or without rosiglitazone (ROSI, a PPARγ agonist). A total of 1143 differentially expressed genes (DEG) due to treatment were detected. The Dynamic Impact Approach (DIA) analysis uncovered the most impacted and induced pathways “fatty acid elongation in mitochondria,” “glycosaminoglycan biosynthesis-keratan sulfate,” and “pentose phosphate pathway.” The data highlights the central role of PPARG in milk fatty acid metabolism via controlling fatty acid elongation, biosynthesis of unsaturated fatty acid, lipid formation, and lipid secretion; furthermore, its role related to carbohydrate metabolism promotes the production of intermediates required for milk fat synthesis. Analysis of upstream regulators indicated that PPARG participates in multiple physiological processes via controlling or cross talking with other key transcription factors such as PPARD and NR1H3 (also known as liver-X-receptor-α). This transcriptome-wide analysis represents the first attempt to better understand the biological relevance of PPARG expression in ruminant mammary cells. Overall, the data underscored the importance of PPARG in mammary lipid metabolism and transcription factor control.

Integrated data mining of transcriptomic and proteomic datasets to predict the secretome of adipose tissue and muscle in ruminants

Adipose tissue and muscle are endocrine organs releasing signalling and mediator proteins termed adipokines and myokines. The identification of the complete set of proteins secreted by adipose tissue and muscle is a challenge to understand the molecular cross-talk between these tissues and to reveal potential targets to control body or muscle composition and metabolism.

Maternal Plane of Nutrition during Late Gestation and Weaning Age Alter Angus × Simmental Offspring Longissimus Muscle Transcriptome and Intramuscular Fat

In model organisms both the nutrition of the mother and the young offspring could induce long-lasting transcriptional changes in tissues. In livestock, such changes could have important roles in determining nutrient use and meat quality. The main objective was to evaluate if plane of maternal nutrition during late-gestation and weaning age alter the offspring's Longissimus muscle (LM) transcriptome, animal performance, and metabolic hormones. Whole-transcriptome microarray analysis was performed on LM samples of early (EW) and normal weaned (NW) Angus × Simmental calves born to grazing cows receiving no supplement [low plane of nutrition (LPN)] or 2.3 kg high-grain mix/day [medium plane of nutrition (MPN)] during the last 105 days of gestation. Biopsies of LM were harvested at 78 (EW), 187 (NW) and 354 (before slaughter) days of age. Despite greater feed intake in MPN offspring, blood insulin was greater in LPN offspring. Carcass intramuscular fat content was greater in EW offspring. Bioinformatics analysis of the transcriptome highlighted a modest overall response to maternal plane of nutrition, resulting in only 35 differentially expressed genes (DEG). However, weaning age and a high-grain diet (EW) strongly impacted the transcriptome (DEG = 167), especially causing a lipogenic program activation. In addition, between 78 and 187 days of age, EW steers had an activation of the innate immune system due presumably to macrophage infiltration of intramuscular fat. Between 187 and 354 days of age (the "finishing" phase), NW steers had an activation of the lipogenic transcriptome machinery, while EW steers had a clear inhibition through the epigenetic control of histone acetylases. Results underscored the need to conduct further studies to understand better the functional outcome of transcriptome changes induced in the offspring by pre- and post-natal nutrition. Additional knowledge on molecular and functional outcomes would help produce more efficient beef cattle.

Effects of capsaicin on adipogenic differentiation in bovine bone marrow mesenchymal stem cell

Capsaicin is a major constituent of hot chili peppers that influences lipid metabolism in animals. In this study, we explored the effects of capsaicin on adipogenic differentiation of bovine bone marrow mesenchymal stem cells (BMSCs) in a dose- and time-dependent manner. The BMSCs were treated with various concentrations of capsaicin (0, 0.1, 1, 5, and 10 μM) for 2, 4, and 6 days. Capsaicin suppressed fat deposition significantly during adipogenic differentiation. Peroxisome proliferator-activated receptor gamma, cytosine-cytosine-adenosine-adenosine-thymidine/enhancer binding protein alpha, fatty acid binding protein 4, and stearoyl-CoA desaturase expression decreased after capsaicin treatment. We showed that the number of apoptotic cells increased in dose- and time-dependent manners. Furthermore, we found that capsaicin increased the expression levels of apoptotic genes, such as B-cell lymphoma 2-associated X protein and caspase 3. Overall, capsaicin inhibits fat deposition by triggering apoptosis.