Discovery of novel and differentially expressed microRNAs between fetal and adult backfat in cattle

Transcriptome analysis of adipose tissue in grazing cattle: Identifying key regulators of fat metabolism

The taste and tenderness of meat are the main determinants of carcass quality in many countries. This study aimed to discuss the mechanisms of intramuscular fat deposition in grazing and house-breeding cattle. We performed transcriptome analysis to characterize messenger RNA and microRNA (miRNA) expression profiles. A total of 456 and 66 differentially expressed genes (DEGs) and differentially expressed (DE) miRNAs were identified in the adipose tissue of grazing and house-breeding cattle. Kyoto Encyclopedia of Genes and Genomes pathway analysis identified the association of DEGs with fatty acid metabolism, fatty acid degradation, peroxisome proliferator-activated receptors signaling pathway, adenosine monophosphate-activated protein kinase signaling pathway, adipocytokine signaling pathway, and the association of DE miRNAs with mitogen-activated protein kinase signaling pathway. Apolipoprotein L domain containing 1, pyruvate dehydrogenase kinase 4, and sphingosine-1-phosphate lyase 1 genes may be the key regulators of fat metabolism in grazing cattle. Finally, we found that miR-211 and miR-331-5p were negatively correlated with the elongation of very long-chain fatty acids protein 6 (ELOVL6), and miR-331-5p might be the new regulator involved in fat metabolism. The results indicated that ELOVL6 participated in various functions and pathways related to fat metabolism. Meanwhile, miR-331-5p, as a new regulator, might play an essential role in this process. Our findings laid a more in-depth and systematic research foundation for the formation mechanism and characteristics of adipose tissue in grazing cattle.

MicroRNA breed and parent-of-origin effects provide insights into biological pathways differentiating cattle subspecies in fetal liver

Introduction: MicroRNAs (miRNAs) play a crucial role in regulating gene expression during key developmental processes, including fetal development. Brahman (Bos taurus indicus) and Angus (Bos taurus taurus) cattle breeds represent two major cattle subspecies with strikingly different phenotypes. Methods: We analyzed miRNA expression in liver samples of purebred and reciprocal crosses of Angus and Brahman to investigate breed and parent-of-origin effects at the onset of accelerated fetal growth. Results: We identified eight novel miRNAs in fetal liver samples and 14 differentially expressed miRNAs (DEMs) between purebred samples. Correlation of gene expression modules and miRNAs by breed and parent-of-origin effects revealed an enrichment of genes associated with breed-specific differences in traits such as heat tolerance (Brahman) and fat deposition (Angus). We demonstrate that genes predicted to be targets of DEMs were more likely to be differentially expressed than non-targets (p-value < 0.05). We identified several miRNAs (bta-miR-187, bta-miR-216b, bta-miR-2284c, bta-miR-2285c, bta-miR-2285cp, bta-miR-2419-3p, bta-miR-2419-5p, and bta-miR-11984) that showed similar correlation patterns as bta-miR-2355-3p, which has been associated with the glutamatergic synapse pathway, a key facilitator of heat tolerance. Furthermore, we report Angus-breed-specific miRNAs (bta-miR-2313-5p, btamiR-490, bta-miR-2316, and bta-miR-11990) that may be involved in fat deposition. Finally, we showed that the DEMs identified in fetal liver are involved in Rap1, MAPK, and Ras signalling pathways, which are important for fetal development, muscle development and metabolic traits such as fat metabolism. Conclusion: Our work sheds light on the miRNA expression patterns that contribute to gene expression differences driving phenotypic differences in indicine and taurine cattle.

Bta-miR-493 Inhibits Bovine Preadipocytes Differentiation by Targeting BMPR1A via the TGFβ/BMP and p38MAPK Signaling Pathways

Fat deposition significantly impacts the meat yield and the meat quality of beef cattle, which is closely associated with the preadipocyte proliferation and differentiation. Bta-miR-493 is expressed differentially in the backfat of Qinchuan cattle of various ages, and it may be involved in the fat accumulation of beef cattle. However, the role and functional mechanism of bta-miR-493 in fat deposition are still unclear. Tissue-specific expression analysis showed that the level of bta-miR-493 was moderately abundant in the adipose tissues of beef cattle. Moreover, the expression of bta-miR-493 in perirenal fat, subcutaneous fat, and longissimus dorsi muscle of Japanese black cattle was significantly higher than that in Qinchuan cattle. EdU staining, cell counting assay, and Oil Red O staining indicated that bta-miR-493 promoted the proliferation of bovine preadipocytes but inhibited their differentiation. The dual luciferase assay and transcriptomic analysis confirmed that bone morphogenetic protein receptor 1A (BMPR1A) is a target gene of bta-miR-493. Moreover, rescue experiments showed that bta-miR-493 could promote bovine preadipocyte proliferation but inhibit their differentiation by targeting BMPR1A via the TGFβ/BMP and p38MAPK signaling pathways. Overall, our findings revealed a bta-miR-493-BMPR1A-TGFβ/BMP/p38MAPK regulatory axis that will serve as a theoretical foundation for the molecular breeding of beef cattle with high intramuscular fat deposition.

CircRNA Profiling Reveals an Abundant circBDP1 that Regulates Bovine Fat Development by Sponging miR-181b/miR-204 Targeting Sirt1/TRARG1

The proliferation and differentiation of preadipocytes is an important factor determining bovine fat development, which is closely related to the feed conversion ratio, carcass traits, and beef quality. The purpose of this study was to identify the effects of candidate circRNA and miRNA on the proliferation and differentiation of bovine preadipocytes in order to provide basic materials for molecular breeding in cattle. circRNA sequencing was performed on bovine adipocyte samples at different differentiation time points, and a total of 1830 differentially expressed circRNAs were identified. Among them, circBDP1, derived from the bovine BDP1 gene, has potential binding sites for miR-204 (known as a regulator of bovine fat development) and miR-181b, which gives us a hint that circBDP1 may regulate bovine fat development by adsorbing miR-204 and miR-181b. Here, our results revealed that circBDP1 overexpression promoted the proliferation and differentiation of bovine preadipocytes. The miRNA profile of bovine adipocytes at different differentiation time points was also analyzed using the small RNA sequencing method, and a total of 89 differentially expressed miRNAs were identified, including miR-204 and miR-181b. As expected, dual-luciferase reporter results showed that circBDP1 competitively adsorbed miR-181b and miR-204. Overexpression and interference of miR-181b in bovine preadipocytes and 3T3-L1 showed that miR-181b promoted the proliferation and differentiation of preadipocytes. Further results displayed that miR-181b and miR-204 simultaneously targeted the SIRT1 gene, and miR-204 also targeted the 3' UTR region of the TRARG1 gene. In summary, this study found that miR-181b and miR-204 were involved in fat development by targeting SIRT1 and TRARG1. The results of this study will lay a foundation for the research of fat development and beef cattle industry.

Integrated analysis of expression profiles with meat quality traits in cattle

MicroRNAs (miRNAs) play a vital role in improving meat quality by binding to messenger RNAs (mRNAs). We performed an integrated analysis of miRNA and mRNA expression profiling between bulls and steers based on the differences in meat quality traits. Fat and fatty acids are the major phenotypic indices of meat quality traits to estimate between-group variance. In the present study, 90 differentially expressed mRNAs (DEGs) and 18 differentially expressed miRNAs (DEMs) were identified. Eighty-three potential DEG targets and 18 DEMs were used to structure a negative interaction network, and 75 matching target genes were shown in this network. Twenty-six target genes were designated as intersection genes, screened from 18 DEMs, and overlapped with the DEGs. Seventeen of these genes enriched to 19 terms involved in lipid metabolism. Subsequently, 13 DEGs and nine DEMs were validated using quantitative real-time PCR, and seven critical genes were selected to explore the influence of fat and fatty acids through hub genes and predict functional association. A dual-luciferase reporter and Western blot assays confirmed a predicted miRNA target (bta-miR-409a and PLIN5). These findings provide substantial evidence for molecular genetic controls and interaction among genes in cattle.

A bovine miRNA, bta-miR-154c, withstands in vitro human digestion but does not affect cell viability of colorectal human cell lines after transfection

Colorectal cancer (CRC) is the third most frequent human cancer with over 1.3 million new cases globally. CRC is a complex disease caused by interactions between genetic and environmental factors; in particular, high consumption of red meat, including beef, is considered a risk factor for CRC initiation and progression. Recent data demonstrate that exogenous microRNAs (miRNAs) entering the body via ingestion could pose an effect on the consumer. In this study, we focused on bovine miRNAs that do not share a seed sequence with humans and mice. We identified bta-miR-154c, a bovine miRNA found in edible parts of beef and predicted via cross-species bioinformatic analysis to affect cancer-related pathways in human cells. When bovine tissue was subjected to cooking and a simulation of human digestion, bta-miR-154c was still detected after all procedures, albeit at reduced concentrations. However, lipofection of bta-miR-154c in three different colorectal human cell lines did not affect their viability as evaluated at various time points and concentrations. These data indicate that bta-miR-154c (a) may affect cancer-related pathways in human cells, (b) can withstand digestion and be detected after all stages of an in vitro digestion protocol, but (c) it does not appear to alter epithelial cell viability after entering human enterocytes, even at supraphysiological amounts. Further experiments will elucidate whether bta-miR-154c exerts a different functional effect on the human gut epithelium, which may cause it to contribute to CRC progression through its consumption.

Effects of nutrient restriction during early or mid-gestation in bovine on placental development and miRNA expression in the cotyledon

The objective of this study was to determine effects of maternal nutrient restriction (NR) during early or mid-gestation on uterine composition and miRNA expression in cotyledons. Primiparous Angus-cross cows (n = 38) were synchronized and inseminated using male sexed semen, blocked by body condition score and body weight (BW), and assigned to treatments. Animals were fed either: control (CON; gain 1 kg/week) or NR (55% maintenance energy and crude protein requirements) based on BW. An initial set of animals were fed either NR (n = 8) or CON (n = 8) from day 30-110 of gestation. A second set of animals were fed CON (n = 8) d 30-190 (CON/CON); NR (n = 7) day 30-110 followed by CON day 110-190 (NR/CON); or CON (n = 7) day 30-110 followed by NR day 110-190 (CON/NR). Cows were harvested on day 110 or 190 of gestation to collect placental tissues. RNA was isolated from cotyledon samples (3 animals/group) prior to microarray analysis using known Bos taurus microRNA sequences. Relative microRNA abundance was analyzed via ANOVA. Maternal NR increased (P < 0.05) cotyledon weight and total placentome surface area irrespective of gestational day. At day 110 of gestation, 51 microRNAs were reduced while 91 microRNAs observed greater abundance (P < 0.05) in NR verses CON cotyledons. At day 190 of gestation, 40 microRNAs were reduced and 26 microRNAs were increased (P < 0.05) in both NR/CON and CON/NR verses CON cotyledons. Top KEGG pathway analysis included: axon guidance, endocytosis, neuroactive ligand receptor interaction, and MAPK signaling pathway. Early-gestation maternal NR altered microRNA abundance to a greater extent than mid-gestation NR.

Yak miR-2285o-3p attenuates hypoxia-induced apoptosis by targeting caspase-3

miRNAs are a class of hairpin‐derived RNAs, 21–24 nucleotides in length, which are involved in a range of biological processes. The bta‐miR‐2285 family has over 40 members spanning the entire bovine genome. We previously found that bta‐miR‐2285o‐3p was highly expressed in yak heart and lung when compared with cattle, which prompted us to investigate its potential function in high‐altitude adaptation of yaks. In this study, we detected wide‐spread high expression of bta‐miR‐2285o‐3p in yak tissues. Further experiments revealed that the protein tyrosine phosphatase receptor type M (PTPRM) gene was the host gene of bta‐miR‐2285o‐3p and that two linked SNPs in bta‐mir‐2285o precursor affected the biogenesis of mature miRNA (bta‐miR‐2285o‐3p). Functional analysis in vitro indicated that bta‐miR‐2285o‐3p attenuated hypoxia‐induced apoptosis by targeting very low‐density lipoprotein receptor (VLDLR), phosphatase and tensin homolog (PTEN) and caspase‐3. Expression level analysis in vivo revealed the high negative Pearson’s correlation between bta‐miR‐2285o‐3p and caspase3 in yak, highlighting the potential important roles of bta‐miR‐2285o‐3p in yak high‐altitude adaptation. Our study provides a typical model for deciphering the function of miRNAs in environmental adaptation.

miR-424 Promotes Bovine Adipogenesis Through an Unconventional Post-Transcriptional Regulation of STK11

Adipose tissue is the largest energy reservoir and secretory organ in the animal body, and is essential for maintaining normal physiological functions and metabolic balance. MicroRNAs regulate the process of adipogenic differentiation through post-transcriptional regulatory mechanisms. In the present study, miR-424 was upregulated during bovine adipocyte differentiation both in vivo and in vitro. The overexpression and interference of miR-424 exhibited the positive regulatory role in the differentiation of bovine adipocytes. Furthermore, miR-424 directly binds to the three prime untranslated region (3' UTR) of serine/threonine kinase 11 (STK11, also called LKB1), a master upstream gene in the AMP-activated protein kinase (AMPK) cascade, and up-regulates its expression. Functional studies showed that the knockdown of STK11 attenuated the pro-adipogenic effect of miR-424. Post-transcriptional regulation of STK11 by miR-424 was mediated potentially in an RNA binding protein (RBP) binding site-dependent manner. In conclusion, our study shows that miR-424 promotes bovine adipogenesis through an unconventional post-transcriptional regulation of STK11, which may serve as a potential target for the regulation of bovine adipogenesis and the improvement of livestock breeding efficiency.

Comparative muscle transcriptome associated with carcass traits of Nellore cattle

Background

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.

Results

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.

Conclusion

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3897-x) contains supplementary material, which is available to authorized users.

MicroRNA-2400 promotes bovine preadipocyte proliferation

MicroRNAs (miRNAs) play critical roles in the proliferation of bovine preadipocytes. miR-2400 is a novel and unique miRNA from bovines. In the present study, we separated and identified preadipocytes from bovine samples. miR-2400 overexpression increased the rate of preadipocyte proliferation, which was analyzed with a combination of EdU and flow cytometry. Simultaneously, functional genes related to proliferation (PCNA, CCND2, CCNB1) were also increased, which was detected by real-time PCR. Furthermore, luciferase reporter assays showed that miR-2400 bound directly to the 3'untranslated regions (3'UTRs) of PRDM11 mRNA. These data suggested that miR-2400 could promote preadipocyte proliferation by targeting PRDM11.

Expression of microRNAs in Horse Plasma and Their Characteristic Nucleotide Composition

MicroRNAs (miRNAs) in blood plasma are stable under high levels of ribonuclease activity and could function in tissue-to-tissue communication, suggesting that they may have distinctive structural characteristics compared with non-circulating miRNAs. In this study, the expression of miRNAs in horse plasma and their characteristic nucleotide composition were examined and compared with non-plasma miRNAs. Highly expressed plasma miRNA species were not part of the abundant group of miRNAs in non-plasma tissues, except for the eca-let-7 family. eca-miR-486-5p, -92a, and -21 were among the most abundant plasma miRNAs, and their human orthologs also belong to the most abundant group of miRNAs in human plasma. Uracil and guanine were the most common nucleotides of both plasma and non-plasma miRNAs. Cytosine was the least common in plasma and non-plasma miRNAs, although levels were higher in plasma miRNAs. Plasma miRNAs also showed higher expression levels of miRNAs containing adenine and cytosine repeats, compared with non-plasma miRNAs. These observations indicate that miRNAs in the plasma have a unique nucleotide composition.

Genome-wide analysis of microRNAs identifies the lipid metabolism pathway to be a defining factor in adipose tissue from different sheep

MicroRNAs are short (17–24 nt) non-coding RNAs that are involved in post-transcriptional regulation of gene expression in multicellular organisms by affecting both the stability and translation of mRNAs. In recent years, deep sequencing of the transcriptome is increasingly being utilized with the promise of higher sensitivity for the identification of differential expression patterns as well as the opportunity to discover new transcripts, including new alternative isoforms and miRNAs. Here, we utilized RNA-seq technology to perform a genome-wide analysis of miRNAs from the adipose tissue of the two species of sheep to look for clues that might explain the fat deposition differences between the sheep. The RNA-seq analysis detected 3132 miRNAs from the adipose tissue of the Small-tail Han and Dorset sheep, of which 2893 were defined as potential new miRNAs. In addition, 54 miRNAs were differentially expressed between the two breeds of sheep. Gene ontology and pathway analyses of the predicted target genes that negatively associated with the differentially expressed miRNAs revealed that there was less active lipid metabolism in the adipose tissue of Small Tail Han sheep. This study can help understand the underling mechanisms responsible for the morphological differences related to fat deposition between two breeds of sheep.

Epigenetic marks: regulators of livestock phenotypes and conceivable sources of missing variation in livestock improvement programs

Improvement in animal productivity has been achieved over the years through careful breeding and selection programs. Today, variations in the genome are gaining increasing importance in livestock improvement strategies. Genomic information alone, however, explains only a part of the phenotypic variance in traits. It is likely that a portion of the unaccounted variance is embedded in the epigenome. The epigenome encompasses epigenetic marks such as DNA methylation, histone tail modifications, chromatin remodeling, and other molecules that can transmit epigenetic information such as non-coding RNA species. Epigenetic factors respond to external or internal environmental cues such as nutrition, pathogens, and climate, and have the ability to change gene expression leading to emergence of specific phenotypes. Accumulating evidence shows that epigenetic marks influence gene expression and phenotypic outcome in livestock species. This review examines available evidence of the influence of epigenetic marks on livestock (cattle, sheep, goat, and pig) traits and discusses the potential for consideration of epigenetic markers in livestock improvement programs. However, epigenetic research activities on farm animal species are currently limited partly due to lack of recognition, funding and a global network of researchers. Therefore, considerable less attention has been given to epigenetic research in livestock species in comparison to extensive work in humans and model organisms. Elucidating therefore the epigenetic determinants of animal diseases and complex traits may represent one of the principal challenges to use epigenetic markers for further improvement of animal productivity.