Integrated analysis of miRNA and mRNA paired expression profiling of prenatal skeletal muscle development in three genotype pigs

Circular RNA transcriptome across multiple tissues reveal skeletal muscle-specific circPSME4 regulating myogenesis

There are significant differences in meat production, growth rate and other traits between Western commercial pigs and Chinese local pigs. Comparative transcriptome approaches have identified many coding and non-coding candidate genes associated various traits. However, the expression and function of circular RNAs (circRNAs) in different pig tissues are largely unknown. In this study, we conducted a comprehensive analysis of the genome-wide circRNA expression profile across ten tissues in Luchuan (a Chinese local breed) and Duroc (a Western commercial breed) pigs. We identified a total of 56,254 circRNAs, of which 42.9 % were not previously annotated. We found that 33.7 % of these circRNAs were differentially expressed. Enrichment analysis revealed that differentially expressed circRNAs might contribute to the phenotypic differentiation between Luchuan and Duroc pigs. We identified 538 tissue-specific circRNAs, most of which were specifically expressed in the brain and skeletal muscle. Competitive endogenous RNA network analysis suggested that skeletal muscle-specific circPSME4 was co-expressed with MYOD1 and targeted by ssc-miR-181d-3p. Functional analysis revealed that circPSME4 knockdown could promote the proliferation and differentiation of myoblasts. Together, our findings provide valuable resources of circRNAs for animal breeding and biomedical research. We demonstrated that circPSME4 is a novel regulator of skeletal muscle development.

Role of Actin-Binding Proteins in Skeletal Myogenesis

Maintenance of skeletal muscle quantity and quality is essential to ensure various vital functions of the body. Muscle homeostasis is regulated by multiple cytoskeletal proteins and myogenic transcriptional programs responding to endogenous and exogenous signals influencing cell structure and function. Since actin is an essential component in cytoskeleton dynamics, actin-binding proteins (ABPs) have been recognized as crucial players in skeletal muscle health and diseases. Hence, dysregulation of ABPs leads to muscle atrophy characterized by loss of mass, strength, quality, and capacity for regeneration. This comprehensive review summarizes the recent studies that have unveiled the role of ABPs in actin cytoskeletal dynamics, with a particular focus on skeletal myogenesis and diseases. This provides insight into the molecular mechanisms that regulate skeletal myogenesis via ABPs as well as research avenues to identify potential therapeutic targets. Moreover, this review explores the implications of non-coding RNAs (ncRNAs) targeting ABPs in skeletal myogenesis and disorders based on recent achievements in ncRNA research. The studies presented here will enhance our understanding of the functional significance of ABPs and mechanotransduction-derived myogenic regulatory mechanisms. Furthermore, revealing how ncRNAs regulate ABPs will allow diverse therapeutic approaches for skeletal muscle disorders to be developed.

Identification and functional analysis of miRNAs in skeletal muscle of juvenile and adult largemouth bass, Micropterus salmoides

MicroRNAs (miRNAs) are considered key regulators to post-transcriptionally regulate gene expression affecting multiple biological activities. However, the developmental process of fish skeletal muscles is regulated by complicated molecular mechanism that has not been completely well-described. In this study, two small RNAs libraries from skeletal muscle of juvenile as well as adult largemouth bass (LMB) were obtained and sequenced using deep sequencing to investigate the development-related miRNAs. We identified an overall number of 486 already recognized miRNAs in addition to 43 novel miRNAs. Comparison of two different skeletal muscle development stages led to the identification of 220 differently expressed miRNAs between juvenile and adult LMB containing 116 up-regulated as well as 104 down-regulated miRNAs. Of them, confirmation of some differently expressed miRNAs was performed via a stem-loop qRT-PCR, which exhibited differently expressed level in juvenile and adult LMB. Furthermore, GO and KEGG enrichment analyses of targets of differently-expressed miRNAs were carried out. Additionally, the analysis of miRNAs-targets interaction network showed that miR-181b-5p_R-1, miR-725 and miR-103 as the nodal miRNAs has over 20 target genes. Moreover, miR-103 could bind the 3'-UTR of actr8, which was validated via dual-luciferase reporter assay. It has been reasonably hypothesized that miR-103 may play a crucial role, which regulate skeletal muscle development of LMB. The present study provides the first identification of miRNA expression profiles at two different skeletal muscle development stages in LMB. Results may be valuable in interpreting the regulatory role miRNAs plays in the growth and developmental process of skeletal muscle and its possible use in LMB breeding.

Profiling and Functional Analysis of mRNAs during Skeletal Muscle Differentiation in Goats

Skeletal myogenesis is a complicated biological event that involves a succession of tightly controlled gene expressions. In order to identify novel regulators of this process, we performed mRNA-Seq studies of goat skeletal muscle satellite cells (MuSCs) cultured under proliferation (GM) and differentiation (DM1/DM5) conditions. A total of 19,871 goat genes were expressed during these stages, 198 of which represented novel transcripts. Notably, in pairwise comparisons at the different stages, 2551 differentially expressed genes (DEGs) were identified (p < 0.05), including 1560 in GM vs. DM1, 1597 in GM vs. DM5, and 959 in DM1 vs. DM5 DEGs. The time-series expression profile analysis clustered the DEGs into eight gene groups, three of which had significantly upregulated and downregulated patterns (p < 0.05). Functional enrichment analysis showed that DEGs were enriched for essential biological processes such as muscle structure development, muscle contraction, muscle cell development, striated muscle cell differentiation, and myofibril assembly, and were involved in pathways such as the MAPK, Wnt and PPAR signaling pathways. Moreover, the expression of eight DEGs (MYL2, DES, MYOG, FAP, PLK2, ADAM, WWC1, and PRDX1) was validated. These findings offer novel insights into the transcriptional regulation of skeletal myogenesis in goats.

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.

Chemerin effect on transcriptome of the porcine endometrium during implantation determined by RNA-sequencing†

It is well known that the body's metabolism and reproduction are closely related. Chemerin is one of many biologically active proteins secreted by the adipose tissue involved in the regulation of the energy homeostasis of the organism. In the present study, RNA-Sequencing (RNA-Seq) was performed to investigate the differentially expressed genes (DEGs), long non-coding RNAs (lncRNAs) and alternatively spliced (AS) transcripts in the cultured in vitro porcine endometrium exposed to chemerin for 24 hours (CHEM; 400 ng/ml) collected during the implantation period (15 to 16 days of gestation). High-throughput sequencing of transcriptomes was performed on the Illumina NovaSeq 6000 platform (Illumina, USA). In the current study, among all 130 DEGs, 58 were up-regulated and were 72 down-regulated in the CHEM-treated group. DEGs were assigned to 73 functional annotations. Twelve identified lncRNAs indicated a difference in the expression profile after CHEM administration. Additionally, we detected 386 differentially AS events encompassed 274 protein-coding genes and 2 lncRNAs. All AS events were divided into 5 alternative splicing types: alternative 3' splice site (A3SS), 5' splice site (A5SS), mutually exclusive exons (MXE), retention intron (RI), and skipping exon (SE). Within all AS events, we identified 42 A3SS, 43 A5SS, 53 MXE, 9 RI, and 239 SE. In summary, CHEM affects the transcriptomic profile of the porcine endometrium, controlling the expression of numerous genes, including those involved in the cell migration and adhesion, angiogenesis, inflammation, and steroidogenesis. It can be assumed that CHEM may be an important factor for a proper course of gestation and embryo development.

A Novel miRNA Y-56 Targeting IGF-1R Mediates the Proliferation of Porcine Skeletal Muscle Satellite Cells Through AKT and ERK Pathways

As a key regulator of gene transcription and post-transcriptional modification, miRNAs play a wide range of roles in skeletal muscle development. Skeletal muscle satellite cells contribute to postnatal growing muscle fibers. Thus, the goal of this study was to explore the effects of novel miRNA Y-56 on porcine skeletal muscle satellite cells (PSCs). We found that Y-56 was highly expressed in porcine muscle tissues, and its expression was higher in Bama Xiang pigs than in Landrace pigs. The EdU assay, cell counting kit-8, and flow cytometry results showed that Y-56 overexpression suppressed cell proliferation and cell cycle, whereas Y-56 inhibition resulted in the opposite consequences. The results of qRT-PCR and Western blot showed that Y-56 remarkably inhibited the expression levels of cyclin-dependent kinase 4 (CDK4), proliferating cell nuclear antigen (PCNA), and cyclin D1. We identified that IGF-1R was a direct target of Y-56 by dual-luciferase reporter assay. Moreover, IGF-1R overexpression promoted the proliferation and cell cycle process of PSCs and upregulated the expression of CDK4, PCNA, and cyclin D1. Conversely, IGF-1R knockdown had the opposite effect. Furthermore, IGF-1R overexpression partially reversed the inhibition of the cell proliferation and cell cycle process of PSCs and the downregulation of the expression of CDK4, PCNA, and Cyclin D1 caused by Y-56 overexpression. Finally, Y-56 inhibited the protein expression levels of p-AKT and p-ERK. Collectively, our findings suggested that Y-56 represses the proliferation and cell cycle process of PSCs by targeting IGF-1R-mediated AKT and ERK pathways.

A Comparative Analysis of Metabolic Profiles of Embryonic Skeletal Muscle from Lantang and Landrace Pigs

Simple Summary

The pig is one of the most important domesticated meat animals. Some studies have revealed that pigs with low meat production show more intense myogenesis at the early stage of embryonic muscle development than pigs with high meat production. Here, by gas chromatography–mass spectrometry GC–MS based metabolomics, we concluded that the nucleotide metabolism and energy metabolism of the longissimuslumborum (LL) were increased in Lantang pigs compared with Landrace pigs, indicating rapid synthesis of nucleic acids and ATP to meet the material and energy requirements of rapid cell proliferation and differentiation in Lantang pigs.

Abstract

Elucidation of the complex regulation of porcine muscle development is key to increasing pork output and improving pork quality. However, the molecular mechanisms involved in early porcine embryonic muscle development in different pig breeds remain largely unknown. Here, GC–MS based metabolomics and metabolomic profiling was used to examine the longissimus lumborum (LL) of the Lantang (LT) and the Landrace (LR) pig at embryonic day 35 (E35). Metabolites showed clear separation between LT and LR, with 40 metabolites having higher abundances in LT and 14 metabolites having lower abundances in LT compared with LR. In addition, these metabolic changes were mainly associated with nucleotide metabolism and energy metabolism, such as purine metabolism, pyrimidine metabolism, the pentose phosphate pathway, and the TCA cycle. More interestingly, the contents of DNA, RNA, and ATP per unit mass of LL tissues were higher in LT, indicating rapid synthesis of nucleic acids and ATP, to meet both the material and energy requirements of rapid cell proliferation and differentiation. Furthermore, enzyme activity associated with the TCA cycle and pentose phosphate pathway, including α-ketoglutaric dehydrogenase (KGDH), malate dehydrogenase (MDH), pyruvate dehydrogenase (PDH), succinate dehydrogenase (SDH), and glucose-6-phosphate dehydrogenase (G6PDH), were higher in LT. Based on these results, we conclude that there are significant differences in nucleotide metabolism and energy metabolism of LL between LT and LR, and we speculate that the enhanced nucleic acid metabolism and energy metabolism in LT can meet the material and energy requirements of rapid cell proliferation and differentiation, making myogenesis more intense in LT compared to LR which might be the metabolic mechanism underlying the distinct skeletal muscle development in the two breeds.

Role of MiR-325-3p in the Regulation of CFL2 and Myogenic Differentiation of C2C12 Myoblasts

Skeletal myogenesis is required to maintain muscle mass and integrity, and impaired myogenesis is causally linked to the etiology of muscle wasting. Recently, it was shown that excessive uptake of saturated fatty acids (SFA) plays a significant role in the pathogenesis of muscle wasting. Although microRNA (miRNA) is implicated in the regulation of myogenesis, the molecular mechanism whereby SFA-induced miRNAs impair myogenic differentiation remains largely unknown. Here, we investigated the regulatory roles of miR-325-3p on CFL2 expression and myogenic differentiation in C2C12 myoblasts. PA impeded myogenic differentiation, concomitantly suppressed CFL2 and induced miR-325-3p. Dual-luciferase analysis revealed that miR-325-3p directly targets the 3'UTR of CFL2, thereby suppressing the expression of CFL2, a crucial factor for actin dynamics. Transfection with miR-325-3p mimic resulted in the accumulation of actin filaments (F-actin) and nuclear Yes-associated protein (YAP) in myoblasts and promoted myoblast proliferation and cell cycle progression. Consequently, miR-325-3p mimic significantly attenuated the expressions of myogenic factors and thereby impaired the myogenic differentiation of myoblasts. The roles of miR-325-3p on CFL2 expression, F-actin modulation, and myogenic differentiation suggest a novel miRNA-mediated regulatory mechanism of myogenesis and PA-inducible miR-325-3p may be a critical mediator between obesity and muscle wasting.

Dynamic transcriptome profiles of postnatal porcine skeletal muscle growth and development

BACKGROUND

Skeletal muscle growth and development are closely associated with the quantity and quality of pork production. We performed a transcriptomic analysis of 12 Longissimus dorsi muscle samples from Tibetan piglets at four postnatal stages of 0, 14, 30, and 60 days using RNA sequencing.

RESULTS

According to the pairwise comparisons between the libraries of the muscle samples at the four postnatal stages, a total of 4115 differentially expressed genes (DEGs) were identified in terms of |log2(fold change)| ≥ 1 and an adjusted P value < 0.01. Short-time series expression miner (STEM) analysis of the DEGs identified eight significantly different expression profiles, which were divided into two clusters based on the expression pattern. DEGs in cluster I displayed a pattern of decreasing to a nadir, and then a rise, and the significantly enriched gene ontology (GO) terms detected using them were involved in multiple processes, of which the cell cycle, immunocyte activation and proliferation, as well as actin cytoskeleton organization, were the top three overrepresented processes based on the GO terms functional classification. DEGs in cluster II displayed a pattern of increasing to a peak, then declining, which mainly contributed to protein metabolism. Furthermore, besides the pathways related to immune system, a few diseases, and protein metabolism, the DEGs in clusters I and II were significantly enriched in pathways related to muscle growth and development, such as the Rap1, PI3K-Akt, AMPK, and mTOR signaling pathways.

CONCLUSIONS

This study revealed GO terms and pathways that could affect the postnatal muscle growth and development in piglets. In addition, this study provides crucial information concerning the molecular mechanisms of muscle growth and development as well as an overview of the piglet transcriptome dynamics throughout the postnatal period in terms of growth and development.

MicroRNA-183-5p protects human derived cell line SH-SY5Y cells from mepivacaine-induced injury

With the gradual recognition of the side effects of local anesthetics, the nerve injury caused by local anesthetics has received growing attention. This research intended to delve into miR-183-5p changes in mepivacaine-mediated SH-SY5Y cell injury, as well as its modulatory mechanism on cell apoptosis. RT-qPCR was adopted for assaying miR-183-5p and PDCD4 mRNA expression. Our team respectively transfected miR-183-5p mimic and inhibitor to enhance or inhibit miR-183-5p function. We employed Western blot for detecting PDCD4 protein levels, as well as flow cytometry and Hoechst 33342/PI double staining for determining cell apoptosis rate. Additionally, our crew applied an ELISA kit for measuring TNF-α, IL-1β, IL-6, and IL-8 contents. The level of reactive oxygen species (ROS) production was examined by the Image-iT LIVE Green ROS detection Kit. As well as dual-luciferase reporter experiment for verifying the targeting link of miR-183-5p with PDCD4. In mepivacaine-induced cell apoptosis in SH-SY5Y cells, miR-183-5p expression was down-regulated. TNF-α, IL-1β, IL-6, and IL-8 contents were elevated. The rate of apoptosis increased visibly, cleaved caspase-3 and Bax levels waxed, whereas Bcl-2 level waned. MiR-183-5p could alleviate the damaging impact of mepivacaine. Dual-luciferase reporter experiments demonstrated that miR-183-5p directly targeted PDCD4. Collectively, we concluded that a high concentration of mepivacaine can cause SH-SY5Y cell damage, miR-183-5p functions crucially in mepivacaine-mediated cell damage. This study provides a theoretical basis for elucidating the mechanism of mepivacaine-induced nerve cell damage, and overexpressed miR-183-5p likely become a novel strategy to combat mepivacaine-induced nerve damage.Abbreviations:miRNA: Micro RNA; PDCD4: Programmed Cell Death 4; MDA: Malondialdehyde; SOD: Superoxide Dismutase; ROS: Reactive Oxygen Species; WT: Wild Type; Mut: Mutant; UTR: Untranslated Region; IL-6: Interleukin-6; IL-1β: Interleukin-1β; TNF-α: Tumor Necrosis Factor-α; IL-8: Interleukin-8; COX-2: Cyclooxygenase-2; iNOS: inducible NOS; MEP: Mepivacaine.

Biological differentiation of traditional Chinese medicine from excessive to deficient syndromes in AIDS: Comparative microRNA microarray profiling and syndrome-specific biomarker identification

Traditional Chinese medicine (TCM) has been widely applied as a supplementary therapy of human immunodeficiency virus infection and acquired immunodeficiency syndrome (HIV/AIDS) in China. TCM has a positive effect on improving the quality of life, prolonging life, and ameliorating the symptoms of HIV/AIDS patients. Yang deficiency of spleen and kidney (YDSK) syndrome is a typical deficient TCM syndrome in AIDS patients, and accumulation of heat-toxicity (AHT) syndrome is a common excessive syndrome in the earlier stage of AIDS. Thus, accurate diagnosis of these two syndromes can improve the targeted treatment effect, and predict the prognosis of the disease. However, the scientific basis of TCM syndromes remains lacking, greatly hindering the accuracy of diagnosis and effectiveness of treatment. In this research, microRNA (miRNA) microarray and quantitative real-time polymerase chain reaction combined with bioinformatics were used for comparative analysis between YDSK and AHT patients. Significantly differential expressed miRNAs (SDE-miRNAs) of each TCM syndrome were identified, including hsa-miR-766-3p and hsa-miR-1260a and so on, as well hsa-miR-6124, hsa-let-7g-5p and so on, for YDSK and AHT, respectively. Biological differences were found between their SDE-miRNAs based on bioinformatics analyses, for example, ErbB signaling pathway mainly linked to AHT, while focal adhesion dominated in YDSK. Syndrome-specific SDE-miRNAs were further identified as potential biomarkers, including hsa-miR-30e-5p, hsa-miR-144-5p for YDSK and hsa-let-7g-5p, hsa-miR-126-3p for AHT, respectively. All of them have laid biological and clinical bases for TCM diagnosis and treatment of AIDS syndrome at the miRNA level, offering potential diagnostic indicators of immune reconstitution.

Skeletal Muscle Transcriptome Analysis of Hanzhong Ma Duck at Different Growth Stages Using RNA-Seq

As one of the most important poultry worldwide, ducks (Anas platyrhynchos) are raised mainly for meat and egg products, and muscle development in ducks is important for meat production. Therefore, an investigation of gene expression in duck skeletal muscle would significantly contribute to our understanding of muscle development. In this study, twenty-four cDNA libraries were constructed from breast and leg muscles of Hanzhong Ma ducks at day 17, 21, 27 of the embryo and postnatal at 6-month-old. High-throughput sequencing and bioinformatics were used to determine the abundances and characteristics of transcripts. A total of 632,172,628 (average 52,681,052) and 637,213,938 (average 53,101,162) reads were obtained from the sequencing data of breast and leg muscles, respectively. Over 71.63% and 77.36% of the reads could be mapped to the Anas platyrhynchos genome. In the skeletal muscle of Hanzhong duck, intron variant (INTRON), synonymous variant (SYNONYMOUS_CODING), and prime 3' UTR variant (UTR_3_PRIME) were the main single nucleotide polymorphisms (SNP) annotation information, and "INTRON", "UTR_3_PRIME", and downstream-gene variant (DOWNSTREAM) were the main insertion-deletion (InDel) annotation information. The predicted number of alternative splicing (AS) in all samples were mainly alternative 5' first exon (transcription start site)-the first exon splicing (TSS) and alternative 3' last exon (transcription terminal site)-the last exon splicing (TTS). Besides, there were 292 to 2801 annotated differentially expressed genes (DEGs) in breast muscle and 304 to 1950 annotated DEGs in leg muscle from different databases. It is worth noting that 75 DEGs in breast muscle and 49 DEGs in leg muscle were co-expressed at all developmental points of comparison, respectively. The RNA-Seq data were confirmed to be reliable by qPCR. The identified DEGs, such as CREBL2, RHEB, GDF6, SHISA2, MYLK2, ACTN3, RYR3, and STMN1, were specially highlighted, indicating their strong associations with muscle development in the Hanzhong Ma duck. KEGG pathway analysis suggested that regulation of actin cytoskeleton, oxidative phosphorylation, and focal adhesion were involved in the development of skeletal muscle. The findings from this study can contribute to future investigations of the growth and development mechanism in duck skeletal muscle.

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.

Role of microRNAs in myogenesis and their effects on meat quality in pig - A review

The demand for food is increasing day by day because of the increasing global population. Therefore, meat, the easiest and largely available source of protein, needs to be produced in large amounts with good quality. The pork industry is a significant shareholder in fulfilling the global meat demands. Notably, myogenesis- development of muscles during embryogenesis- is a complex mechanism which culminates in meat production. But the molecular mechanisms which govern the myogenesis are less known. The involvement of miRNAs in myogenesis and meat quality, which depends on factors such as myofiber composition and intramuscular fat contents which determine the meat color, flavor, juiciness, and water holding capacity, are being extrapolated to increase both the quantity and quality of pork. Various kinds of microRNAs (miRNAs), miR-1, miR-21, miR22, miR-27, miR-34, miR-127, miR-133, miR-143, miR-155, miR-199, miR-206, miR-208, miR-378, and miR-432 play important roles in pig skeletal muscle development. Further, the quality of meat also depends upon myofiber which is developed through the expression of different kinds of miRNAs at different stages. This review will focus on the mechanism of myogenesis, the role of miRNAs in myogenesis, and meat quality with a focus on the pig.

MiR-743a-5p regulates differentiation of myoblast by targeting Mob1b in skeletal muscle development and regeneration

The microRNAs (miRNAs) play an important role in regulating myogenesis by targeting mRNA. However, the understanding of miRNAs in skeletal muscle development and diseases is unclear. In this study, we firstly performed the transcriptome profiling in differentiating C2C12 myoblast cells. Totally, we identified 187 miRNAs and 4260 mRNAs significantly differentially expressed that were involved in myoblast differentiation. We carried out validation of microarray data based on 5 mRNAs and 5 miRNAs differentially expressed and got a consistent result. Then we constructed and validated the significantly up- and down-regulated mRNA-miRNA interaction networks. Four interaction pairs (miR-145a-5p-Fscn1, miR-200c-5p-Tmigd1, miR-27a-5p-Sln and miR-743a-5p-Mob1b) with targeted relationships in differentiated myoblast cells were demonstrated. They are all closely related to myoblast development. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated cell cycle signals important for exploring skeletal muscle development and disease. Functionally, we discovered that miR-743a targeting gene Mps One Binder Kinase Activator-Like 1B (Mob1b) gene in differentiated C2C12. The up-regulated miR-743a can promote the differentiation of C2C12 myoblast. While the down-regulated Mob1b plays a negative role in differentiation. In addition, the expression profile of miR-743a and Mob1b are consistent with skeletal muscle recovery after Cardiotoxin (CTX) injury. Our study revealed that miR-743a-5p regulates myoblast differentiation by targeting Mob1b involved in skeletal muscle development and regeneration. Our findings made a further exploration for mechanisms in myogenesis and might provide potential possible miRNA-based target therapies for skeletal muscle regeneration and disease in the near future.

Identification and Characterization of lncRNAs Related to the Muscle Growth and Development of Japanese Flounder (Paralichthys olivaceus)

Long noncoding RNAs (lncRNAs) play an important role in many life activities, but the expression pattern and function of lncRNAs in Japanese flounder skeletal muscle are unclear. In this study, 751 lncRNAs were selected from skeletal muscle in different development stages of the Japanese flounder [stage A: larval 7 days post hatching (dph); stage B: juvenile about 90 dph; stage C (female) and stage D (male): adult about 24 months] using coding potential analysis methods. In total, 232, 211, 194, 28, 29, and 14 differentially expressed lncRNAs and 9549, 8673, 9181, 1821, 1080, and 557 differentially expressed mRNAs were identified in comparisons of A versus B, A versus C, A versus D, B versus C, B versus D, and C versus D, respectively. We identified the cis- and trans-regulatory target genes of differentially expressed lncRNAs, and lncRNA-gene interaction networks were constructed using the Cytoscape program. In total, there were 200, 200, 200, 93, 47, and 11 cis-regulation relationships, and 29, 19, 24, 38, 8, and 47 trans-regulation relationships in the comparisons between A versus B, A versus C, A versus D, B versus C, B versus D, and C versus D, respectively. These results indicate that lncRNA may participate in the development of Japanese flounder skeletal muscle through cis- or trans-acting mechanisms, thus providing a scientific basis for further study of the biological function of lncRNA in Japanese flounder skeletal muscle. Based on these relationships, functional annotation of the related lncRNAs was performed by gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Differentially expressed genes associated with muscle development were enriched in multiple pairs of comparisons (e.g., differentially expressed genes LOC109640370, LOC109634180, LOC109643555, rusc1, and LOC109626999 were enriched in the actin-binding term, and differentially expressed genes LOC109640370, was, LOC109644970, LOC109643555, and itga9 were enriched in the regulation of the actin cytoskeleton pathway in the KEGG pathway analysis in the comparison of stages C and D). We predicted lncRNA-mRNA, miRNA-mRNA, and lncRNA-miRNA regulatory relationships and constructed interactive networks using Cytoscape software. Co-expression networks show that most lncRNAs interact with one or two predicted miRNAs involved in muscle growth and development. These results provide a basis for further study of the function of lncRNAs on skeletal muscle in different developmental stages of Japanese flounder.

Longitudinal epitranscriptome profiling reveals the crucial role of N6-methyladenosine methylation in porcine prenatal skeletal muscle development

N⁶-methyladenosine (m⁶A) represents the most abundantly occurring mRNA modification and is involved in the regulation of skeletal muscle development. However, the status and function of m⁶A methylation in prenatal myogenesis remains unclear. In this study, we first demonstrated that knockdown of METTL14, an m⁶A methyltransferase, inhibited the differentiation and promoted the proliferation of C2C12 myoblast cells. Then, using a refined m⁶A-specific methylated RNA immunoprecipitation (RIP) with next generation sequencing (MeRIP-seq) method that is optimal for use with samples containing small amounts of RNA, we performed transcriptome-wide m⁶A profiling for six prenatal skeletal muscle developmental stages spanning two important waves of porcine myogenesis. The results revealed that, along with a continuous decrease in the mRNA expression of the m⁶A reader protein insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1), the m⁶A methylome underwent highly dynamic changes across different development stages, with most of the affected genes being enriched in pathways related to skeletal muscle development. RNA immunoprecipitation confirmed that IGF2BP1 targets 76 genes involved in pathways associated with muscle development, including the key marker genes MYH2 and MyoG. Moreover, small interfering RNA (siRNA)-mediated knockdown of IGF2BP1 induced phenotypic changes in C2C12 myoblasts similar to those observed with knockdown of METTL14. In conclusion, we clarified the dynamics of m⁶A methylation and identified key genes involved in the regulatory network of porcine skeletal muscle development.

Analysis of miRNAs and their target genes associated with mucosal damage caused by transport stress in the mallard duck intestine

Bowel health is an important factor for duck rearing that has been linked to feed uptake and growth and death rates. Because the regulatory networks associated with acute stress-mediated injury in the duck gastrointestinal tract have not clearly elucidated, we aimed to explore potential miRNA-mRNA pairs and their regulatory roles in oxidative stress injury caused by transport stress. Here, 1-day-old mallard ducklings from the same breeder flock were collected and transported for 8 h, whereas the control group was not being transported. Various parameters reflecting oxidative stress and the tissue appearance of the intestine were assessed. The data showed that the plasma T-AOC and SOD concentrations were decreased in the transported ducklings. The intestine of the transported ducklings also displayed significant damage. High-throughput sequencing of the intestine revealed 44 differentially expressed miRNAs and 75 differentially expressed genes, which constituted 344 miRNA-mRNA pairs. KEGG pathway analysis revealed that the metabolic, FoxO signaling, influenza A and TGF-β signaling pathways were mainly involved in the mechanism underlying the induction of intestinal damage induced by simulated transport stress in ducks. A miRNA-mRNA pair, miR-217-5p/CHRDL1, was selected to validate the miRNA-mRNA negative relationship, and the results showed that miR-217-5p could influence CHRDL1 expression. This study provides new useful information for future research on the regulatory network associated with mucosal damage in the duck intestine.

miR-206 knockout shows it is critical for myogenesis and directly regulates newly identified target mRNAs

The muscle specific miRNA, miR-206, is important for the process of myogenesis; however, studying the function of miR-206 in muscle development and differentiation still proves challenging because the complement of mRNA targets it regulates remains undefined. In addition, miR-206 shares close sequence similarity to miR-1, another muscle specific miRNA, making it hard to study the impact of miR-206 alone in cell culture models. Here we used CRISPR/Cas9 technology to knockout miR-206 in C2C12 muscle cells. We show that knocking out miR-206 significantly impairs and delays differentiation and myotube formation, revealing that miR-206 alone is important for myogenesis. In addition, we use an experimental affinity purification technique to identify new mRNA targets of miR-206 in C2C12 cells. We identified over one hundred mRNAs as putative miR-206 targets. Functional experiments on six of these targets indicate that Adam19, Bgn, Cbx5, Smarce1, and Spg20 are direct miR-206 targets in C2C12 cells. Our data show a unique and important role for miR-206 in myogenesis.

Biological Differentiation of Dampness-Heat Syndromes in Chronic Hepatitis B: From Comparative MicroRNA Microarray Profiling to Biomarker Identification

Increasing interest is aroused by traditional Chinese medicine (TCM) treatment of chronic hepatitis B (CHB) based on specific TCM syndrome. As the most common CHB syndromes, spleen-stomach dampness-heat (SSDH) syndrome and liver-gallbladder dampness-heat (LGDH) syndrome are still apt to be confused in TCM diagnosis, greatly hindering the stable exertion of TCM effectiveness. It is urgently needed to provide objective and biological evidences for differentiation and identification of the two significant syndromes. In this study, microRNA (miRNA) microarray analyses coupled with bioinformatics were employed for comparative miRNA profiling of SSDH and LGDH patients. It was found that the two syndromes had both the same and different significantly differentially expressed miRNAs (SDE-miRNAs). Commonness and specificity were also both found between their SDE-miRNA-based bioinformatics analyses, including Hierarchical Clustering, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and miRNA-GO/pathway networks. Furthermore, syndrome-specific SDE-miRNAs were identified as the potential biomarkers, including hsa-miR-1273g-3p and hsa-miR-4419b for SSDH as well as hsa-miR-129-1-3p and hsa-miR-129-2-3p for LGDH. All these laid biological and clinical bases for classification and diagnosis of the two significant CHB dampness-heat syndromes including SSDH and LGDH, providing more opportunities for better application of TCM efficacy and superiority in CHB treatment.

Transcriptome analysis of miRNA and mRNA in the livers of pigs with highly diverged backfat thickness

Fat deposition is very important in pig production, and its mechanism is not clearly understood. MicroRNAs (miRNAs) play critical roles in fat deposition and energy metabolism. In the current study, we investigated the mRNA and miRNA transcriptome in the livers of Landrace pigs with extreme backfat thickness to explore miRNA-mRNA regulatory networks related to lipid deposition and metabolism. A comparative analysis of liver mRNA and miRNA transcriptomes from pigs (four pigs per group) with extreme backfat thickness was performed. We identified differentially expressed genes from RNA-seq data using a Cufflinks pipeline. Seventy-one differentially expressed genes (DEGs), including twenty-eight well annotated on the porcine reference genome genes, were found. The upregulation genes in pigs with higher backfat thickness were mainly involved in fatty acid synthesis, and included fatty acid synthase (FASN), glucokinase (GCK), phosphoglycerate dehydrogenase (PHGDH), and apolipoprotein A4 (APOA4). Cytochrome P450, family 2, subfamily J, polypeptide 34 (CYP2J34) was lower expressed in pigs with high backfat thickness, and is involved in the oxidation of arachidonic acid. Moreover, 13 differentially expressed miRNAs were identified. Seven miRNAs were associated with fatty acid synthesis, lipid metabolism, and adipogenic differentiation. Based on comprehensive analysis of the transcriptome of both mRNAs and miRNAs, an important regulatory network, in which six DEGs could be regulated by differentially expressed miRNAs, was established for fat deposition. The negative correlate in the regulatory network including, miR-545-5p and GRAMD3, miR-338 and FASN, and miR-127, miR-146b, miR-34c, miR-144 and THBS1 indicate that direct suppressive regulation may be involved in lipid deposition and energy metabolism. Based on liver mRNA and miRNA transcriptomes from pigs with extreme backfat thickness, we identified 28 differentially expressed genes and 13 differentially expressed miRNAs, and established an important miRNA-mRNA regulatory network. This study provides new insights into the molecular mechanisms that determine fat deposition in pigs.

Comprehensive Profiles of mRNAs and miRNAs Reveal Molecular Characteristics of Multiple Organ Physiologies and Development in Pigs

The pig (Sus scrofa) is not only an important livestock animal but also widely used as a biomedical model. However, the understanding of the molecular characteristics of organs and of the developmental skeletal muscle of the pig is severely limited. Here, we performed a comprehensive transcriptome profiling of mRNAs and miRNAs across nine tissues and three skeletal muscle developmental stages in the Guizhou miniature pig. The reproductive organs (ovary and testis) had greater transcriptome complexity and activity than other tissues, and the highest transcriptome similarity was between skeletal muscle and heart (R = 0.79). We identified 1,819 mRNAs and 96 miRNAs to be tissue-specific in nine organs. Testis had the largest number of tissue-specific mRNAs (992) and miRNAs (40). Only 15 genes and two miRNAs were specifically expressed in skeletal muscle and fat, respectively. During postnatal skeletal muscle development, the mRNAs associated with focal adhesion, Notch signaling, protein digestion, and absorption pathways were up-regulated from D0 to D30 and then down-regulated from D30 and D240, while genes with opposing expression patterns were significantly enriched in the oxidative phosphorylation and proteasome pathways. The miRNAs mainly regulated genes associated with insulin, Wnt, fatty acid biosynthesis, Notch, MAPK, TGF-beta, insulin secretion, ECM-receptor interaction, focal adhesion, and calcium signaling pathways. We also identified 37 new miRNA-mRNA interaction pairs involved in skeletal muscle development. Overall, our data not only provide a rich resource for understanding pig organ physiology and development but also aid the study of the molecular functions of mRNA and miRNA in mammals.

Comparison of microRNAs in adipose and muscle tissue from seven indigenous Chinese breeds and Yorkshire pigs

Elucidation of the pig microRNAome is essential for interpreting functional elements of the genome and understanding the genetic architecture of complex traits. Here, we extracted small RNAs from skeletal muscle and adipose tissue, and we compared their expression levels between one Western breed (Yorkshire) and seven indigenous Chinese breeds. We detected the expression of 172 known porcine microRNAs (miRNAs) and 181 novel miRNAs. Differential expression analysis found 92 and 12 differentially expressed miRNAs in adipose and muscle tissue respectively. We found that different Chinese breeds shared common directional miRNA expression changes compared to Yorkshire pigs. Some miRNAs differentially expressed across multiple Chinese breeds, including ssc-miR-129-5p, ssc-miR-30 and ssc-miR-150, are involved in adipose tissue function. Functional enrichment analysis revealed that the target genes of the differentially expressed miRNAs are associated mainly with signaling pathways rather than metabolic and biosynthetic processes. The miRNA-target gene and miRNA-phenotypic traits networks identified many hub miRNAs that regulate a large number of target genes or phenotypic traits. Specifically, we found that intramuscular fat content is regulated by the greatest number of miRNAs in muscle tissue. This study provides valuable new candidate miRNAs that will aid in the improvement of meat quality and production.

Differentially expressed mRNAs, proteins and miRNAs associated to energy metabolism in skeletal muscle of beef cattle identified for low and high residual feed intake

Background

Feed efficiency is one of the most important parameters that affect beef production costs. The energy metabolism of skeletal muscle greatly contributes to variations in feed efficiency. However, information regarding differences in proteins involved in the energy metabolism of the skeletal muscle in beef cattle divergently identified for feed efficiency is scarce. In this study, we aimed to investigate energy metabolism of skeletal muscle of Nellore beef cattle, identified for low and high residual feed intake using a proteomics approach. We further assessed the expression of candidate microRNAs as a one of the possible mechanisms controlling the biosynthesis of the proteins involved in energy metabolism that were differentially abundant between high and low residual feed intake animals.

Results

A greater abundance of 14–3-3 protein epsilon (P = 0.01) was observed in skeletal muscle of residual feed intake (RFI) high animals (RFI-High). Conversely, a greater abundance of Heat Shock Protein Beta 1 (P < 0.01) was observed in the skeletal muscle of RFI-Low cattle. A greater mRNA expression of YWHAE, which encodes the 14–3-3 protein epsilon, was also observed in the skeletal muscle of RFI-High animals (P = 0.01). A lower mRNA expression of HSPB1, which encodes the Heat Shock Protein Beta 1, was observed in the skeletal muscle of RFI-High animals (P = 0.01). The miR-665 was identified as a potential regulator of the 14–3-3 protein epsilon, and its expression was greater in RFI-Low animals (P < .001). A greater expression of miR-34a (P = 0.01) and miR-2899 (P < .001) was observed in the skeletal muscle of RFI-High animals, as both miRNAs were identified as potential regulators of HSPB1 expression.

Conclusion

Our results show that Nellore cattle divergently identified for feed efficiency by RFI present changes in the abundance of proteins involved in energy expenditure in skeletal muscle. Moreover, our data point towards that miR-665, miR34a and miR-2899 are likely involved in controlling both 14-3-3 epsilon and HSPB1 proteins identified as differentially abundant in the skeletal muscle of RFI-High and RFI-Low Nellore cattle.

Electronic supplementary material

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

Integrative microRNA-mRNA Analysis of Muscle Tissues in Qianhua Mutton Merino and Small Tail Han Sheep Reveals Key Roles for oar-miR-655-3p and oar-miR-381-5p

The Qianhua Mutton Merino (QHMM) is a new variety of sheep (Ovis aries) with improved meat performance compared with the traditional Small Tail Han (STH) sheep variety. We recently reported the transcriptome profiling of longissimus muscle tissues between QHMM and STH sheep. In the present study, we aimed to evaluate key micro (mi)RNA-mRNA networks associated with sheep muscle growth and development. We used miRNA sequencing to obtain longissimus muscle miRNA profiles from QHMM and STH sheep. We identified a total of 153 known sheep miRNAs, of which 4 were differentially expressed (DE) between the 2 sheep varieties. We combined these results with mRNA library data to build an miRNA-mRNA network, including 26 target genes of the 4 DE miRNAs. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that 26 target genes were significantly enriched in 86 biological processes, including muscle organogenesis, myoblast migration, cell proliferation, and adipose tissue development, and in 9 metabolic pathways, including carbohydrate, nucleotide, and amino acid metabolic pathways. oar-miR-655-3p and its target gene ACSM3 and oar-miR-381-5p and its target gene ABAT were selected for subsequent analysis based on GO and KEGG analyses. The binding sites of oar-miR-655-3p with ACSM3 and oar-miR-381-5p with ABAT were validated by a dual-luciferase reporter gene detection system. This represents the first integrative analysis of miRNA-mRNA networks in QHMM and STH muscles and suggests that DE miRNAs, especially oar-miR-655-3p and oar-miR-381-5p, play crucial roles in muscle growth and development.

Dysregulated miRNAome and Proteome of PPRV Infected Goat PBMCs Reveal a Coordinated Immune Response

In this study, the miRNAome and proteome of virulent Peste des petits ruminants virus (PPRV) infected goat peripheral blood mononuclear cells (PBMCs) were analyzed. The identified differentially expressed miRNAs (DEmiRNAs) were found to govern genes that modulate immune response based on the proteome data. The top 10 significantly enriched immune response processes were found to be governed by 98 genes. The top 10 DEmiRNAs governing these 98 genes were identified based on the number of genes governed by them. Out of these 10 DEmiRNAs, 7 were upregulated, and 3 were downregulated. These include miR-664, miR-2311, miR-2897, miR-484, miR-2440, miR-3533, miR-574, miR-210, miR-21-5p, and miR-30. miR-664 and miR-484 with proviral and antiviral activities, respectively, were upregulated in PPRV infected PBMCs. miR-210 that inhibits apoptosis was downregulated. miR-21-5p that decreases the sensitivity of cells to the antiviral activity of IFNs and miR-30b that inhibits antigen processing and presentation by primary macrophages were downregulated, indicative of a strong host response to PPRV infection. miR-21-5p was found to be inhibited on IPA upstream regulatory analysis of RNA-sequencing data. This miRNA that was also highly downregulated and was found to govern 16 immune response genes in the proteome data was selected for functional validation vis-a-vis TGFBR2 (TGF-beta receptor type-2). TGFBR2 that regulates cell differentiation and is involved in several immune response pathways was found to be governed by most of the identified immune modulating DEmiRNAs. The decreased luciferase activity in Dual Luciferase Reporter Assay indicated specific binding of miR-21-5p and miR-484 to their target thus establishing specific binding of the miRNAs to their targets.This is the first report on the miRNAome and proteome of virulent PPRV infected goat PBMCs.

Comprehensive Analysis of MicroRNA⁻Messenger RNA from White Yak Testis Reveals the Differentially Expressed Molecules Involved in Development and Reproduction

Testis development is a vital and tightly regulated process in mammals. Understanding the biological mechanisms underlying testis development will benefit the animal reproduction industry. Expression changes in microRNA and messenger RNA in response to dynamic regulation effects have been associated with this process. However, very little is known about the roles of these molecules in yak development. Using whole-genome small RNA and messenger RNA sequencing, we performed a comprehensive analysis of the microRNA–messenger RNA interaction network expression in the testicles of Tianzhu white yaks during three developmental stages. Using Short Time-series Expression Miner analysis we identified 589 differentially expressed microRNAs (DERs) and 3383 differentially expressed messenger RNAs (DEGs) in the three age groups. A total of 93 unique DEGs are primarily involved in reproduction and testis development. Subsequently, four integration networks were constructed according to the DEGs and DERs in three biological processes. Nineteen DEGs were potentially regulated by 60 DERs, of which miR-574 and target gene AURKA played a crucial role in yak testis development and reproduction. The results of this study provide a basis for further exploration of the microRNA–messenger RNA interactions in testis development and reproduction and aid in uncovering the molecular mechanisms of spermatogenesis in male mammals.

A new approach of gene co-expression network inference reveals significant biological processes involved in porcine muscle development in late gestation

The integration of genetic information in the cellular and nuclear environments is crucial for deciphering the way in which the genome functions under different physiological conditions. Experimental techniques of 3D nuclear mapping, a high-flow approach such as transcriptomic data analyses, and statistical methods for the development of co-expressed gene networks, can be combined to develop an integrated approach for depicting the regulation of gene expression. Our work focused more specifically on the mechanisms involved in the transcriptional regulation of genes expressed in muscle during late foetal development in pig. The data generated by a transcriptomic analysis carried out on muscle of foetuses from two extreme genetic lines for birth mortality are used to construct networks of differentially expressed and co-regulated genes. We developed an innovative co-expression networking approach coupling, by means of an iterative process, a new statistical method for graph inference with data of gene spatial co-localization (3D DNA FISH) to construct a robust network grouping co-expressed genes. This enabled us to highlight relevant biological processes related to foetal muscle maturity and to discover unexpected gene associations between IGF2, MYH3 and DLK1/MEG3 in the nuclear space, genes that are up-regulated at this stage of muscle development.

Genome-wide DNA methylation profiles of porcine ovaries in estrus and proestrus

DNA methylation is an important epigenetic modification involved in the estrous cycle and the regulation of reproduction. Here, we investigated the genome-wide profiles of DNA methylation in porcine ovaries in proestrus and estrus using methylated DNA immunoprecipitation sequencing. The results showed that DNA methylation was enriched in intergenic and intron regions. The methylation levels of coding regions were higher than those of the 5'- and 3'-flanking regions of genes. There were 4,813 differentially methylated regions (DMRs) of CpG islands in the estrus vs. proestrus ovarian genomes. Additionally, 3,651 differentially methylated genes (DMGs) were identified in pigs in estrus and proestrus. The DMGs were significantly enriched in biological processes and pathways related to reproduction and hormone regulation. We identified 90 DMGs associated with regulating reproduction in pigs. Our findings can serve as resources for DNA methylome research focused on porcine ovaries and further our understanding of epigenetically regulated reproduction in mammals.

Integrated miRNA-mRNA analysis reveals regulatory pathways underlying the curly fleece trait in Chinese tan sheep

Background

Tan sheep is an indigenous Chinese breed well known for its beautiful curly fleece. One prominent breed characteristic of this sheep breed is that the degree of curliness differs markedly between lambs and adults, but the molecular mechanisms regulating the shift are still not well understood. In this study, we identified 49 differentially expressed (DE) microRNAs (miRNAs) between Tan sheep at the two stages through miRNA-seq, and combined the data with that in our earlier Suppression Subtractive Hybridization cDNA (SSH) library study to elucidate the mechanisms underlying curly fleece formation.

Results

Thirty-six potential miRNA-mRNA target pairs were identified using computational methods, including 25 DE miRNAs and 10 DE genes involved in the MAPK signaling pathway, steroid biosynthesis and metabolic pathways. With the differential expressions between lambs and adults confirmed by qRT-PCR, some miRNAs were already annotated in the genome, but some were novel miRNAs. Inhibition of KRT83 expression by miR-432 was confirmed by both gene knockdown with siRNA and overexpression, which was consistent with the miRNAs and targets prediction results.

Conclusion

Our study represents the comprehensive analysis of mRNA and miRNA in Tan sheep and offers detailed insight into the development of curly fleece as well as the potential mechanisms controlling curly hair formation in humans.

Electronic supplementary material

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

Examining the Genetic Background of Porcine Muscle Growth and Development Based on Transcriptome and miRNAome Data

Recently, selection in pigs has been focused on improving the lean meat content in carcasses; this focus has been most evident in breeds constituting a paternal component in breeding. Such sire-breeds are used to improve the meat quantity of cross-breed pig lines. However, even in one breed, a significant variation in the meatiness level can be observed. In the present study, the comprehensive analysis of genes and microRNA expression profiles in porcine muscle tissue was applied to identify the genetic background of meat content. The comparison was performed between whole gene expression and miRNA profiles of muscle tissue collected from two sire-line pig breeds (Pietrain, Hampshire). The RNA-seq approach allowed the identification of 627 and 416 differentially expressed genes (DEGs) between pig groups differing in terms of loin weight between Pietrain and Hampshire breeds, respectively. The comparison of miRNA profiles showed differential expression of 57 microRNAs for Hampshire and 34 miRNAs for Pietrain pigs. Next, 43 genes and 18 miRNAs were selected as differentially expressed in both breeds and potentially related to muscle development. According to Gene Ontology analysis, identified DEGs and microRNAs were involved in the regulation of the cell cycle, fatty acid biosynthesis and regulation of the actin cytoskeleton. The most deregulated pathways dependent on muscle mass were the Hippo signalling pathway connected with the TGF-β signalling pathway and controlling organ size via the regulation of ubiquitin-mediated proteolysis, cell proliferation and apoptosis. The identified target genes were also involved in pathways such as the FoxO signalling pathway, signalling pathways regulating pluripotency of stem cells and the PI3K-Akt signalling pathway. The obtained results indicate molecular mechanisms controlling porcine muscle growth and development. Identified genes (SOX2, SIRT1, KLF4, PAX6 and genes belonging to the transforming growth factor beta superfamily) could be considered candidate genes for determining muscle mass in pigs.

Comparing two approaches of miR-34a target identification, biotinylated-miRNA pulldown vs miRNA overexpression

microRNAs (miRNAs) are critical regulators of gene expression. For elucidating functional roles of miRNAs, it is critical to identify their direct targets. There are debates about whether pulldown of biotinylated miRNA mimics can be used to identify miRNA targets or not. Here we show that biotin-labelled miR-34a can be loaded to AGO2, and AGO2 immunoprecipitation can pulldown biotinylated miR-34a (Bio-miR pulldown). RNA-sequencing (RNA-seq) of the Bio-miR pulldown RNAs efficiently identified miR-34a mRNA targets, which could be verified with luciferase assays. In contrast to the approach of Bio-miR pulldown, RNA-seq of miR-34a overexpression samples had limited value in identifying direct targets of miR-34a. It seems that pulldown of 3'-Biotin-tagged miRNA can identify bona fide microRNA targets at least for miR-34a.

Genome-wide profiling of Sus scrofa circular RNAs across nine organs and three developmental stages

The spatio-temporal expression patterns of Circular RNA (circRNA) across organs and developmental stages are critical for its function and evolution analysis. However, they remain largely unclear in mammals. Here, we comprehensively analysed circRNAs in nine organs and three skeletal muscles of Guizhou miniature pig (S. scrofa), a widely used biomedical model animal. We identified 5,934 circRNAs and analysed their molecular properties, sequence conservation, spatio-temporal expression pattern, potential function, and interaction with miRNAs. S. scrofa circRNAs show modest sequence conservation with human and mouse circRNAs, are flanked by long introns, exhibit low abundance, and are expressed dynamically in a spatio-temporally specific manner. S. scrofa circRNAs show the greatest abundance and complexity in the testis. Notably, 31% of circRNAs harbour well-conserved canonical miRNA seed matches, suggesting that some circRNAs act as miRNAs sponges. We identified 149 circRNAs potentially associated with muscle growth and found that their host genes were significantly involved in muscle development, contraction, chromatin modification, cation homeostasis, and ATP hydrolysis-coupled proton transport; moreover, this set of genes was markedly enriched in genes involved in tight junctions and the calcium signalling pathway. Finally, we constructed the first public S. scrofa circRNA database, allowing researchers to query comprehensive annotation, expression, and regulatory networks of circRNAs.

Integrated analysis of mRNA and miRNA expression profiles in livers of Yimeng black pigs with extreme phenotypes for backfat thickness

Fat deposition is an important economic trait in farm animals as well as obesity related diseases in humans, and the liver is a central organ involved in regulating lipid synthesis and metabolism in mammals. In this study, the pig liver transcriptome of two groups (H and L) showing differences in backfat thickness were profiled using RNA-Seq and miRNA-Seq to further explore the molecular mechanism of fat deposition. A total of 238 differentially expressed genes (DEGs) and 58 differentially expressed miRNAs were identified between the H and L group. These genes and miRNAs were functionally related to lipid metabolism, including CYP1A1/2, HMGCS2, ACSS2, UBE2L6, miR-27a, and miR-31. Functional enrichment analysis revealed that genes associated with oxidative stress might be responsible for fat deposition in pigs. Two miRNA-mRNA interaction networks involved in lipid metabolism were identified, and these provided new insights into the molecular regulation that determines fat content in these pigs. Overall, our study furthers our understanding of the molecular mechanisms involved in fat deposition, and these results may help in the design of selection strategies to improve the quality of pork meat and to combat obesity in humans.

Dynamic transcriptomic analysis in hircine longissimus dorsi muscle from fetal to neonatal development stages

Muscle growth and development from fetal to neonatal stages consist of a series of delicately regulated and orchestrated changes in expression of genes. In this study, we performed whole transcriptome profiling based on RNA-Seq of caprine longissimus dorsi muscle tissue obtained from prenatal stages (days 45, 60, and 105 of gestation) and neonatal stage (the 3-day-old newborn) to identify genes that are differentially expressed and investigate their temporal expression profiles. A total of 3276 differentially expressed genes (DEGs) were identified (Q value < 0.01). Time-series expression profile clustering analysis indicated that DEGs were significantly clustered into eight clusters which can be divided into two classes (Q value < 0.05), class I profiles with downregulated patterns and class II profiles with upregulated patterns. Based on cluster analysis, GO enrichment analysis found that 75, 25, and 8 terms to be significantly enriched in biological process (BP), cellular component (CC), and molecular function (MF) categories in class I profiles, while 35, 21, and 8 terms to be significantly enriched in BP, CC, and MF in class II profiles. KEGG pathway analysis revealed that DEGs from class I profiles were significantly enriched in 22 pathways and the most enriched pathway was Rap1 signaling pathway. DEGs from class II profiles were significantly enriched in 17 pathways and the mainly enriched pathway was AMPK signaling pathway. Finally, six selected DEGs from our sequencing results were confirmed by qPCR. Our study provides a comprehensive understanding of the molecular mechanisms during goat skeletal muscle development from fetal to neonatal stages and valuable information for future studies of muscle development in goats.

Dynamic mRNA and miRNA expression analysis in response to hypoxia and reoxygenation in the blunt snout bream (Megalobrama amblycephala)

Adaptation to hypoxia is a complex process involving various pathways and regulation mechanisms. A better understanding of the genetic influence on these mechanisms could permit selection for hypoxia-sensitive fish. To aid this understanding, an integrated analysis of miRNA and mRNA expression was performed in Megalobrama amblycephala under four acute hypoxia and reoxygenation stages. A number of significantly differentially-expressed miRNAs and genes associated with oxidative stress were identified, and their functional characteristics were revealed by GO function and KEGG pathway analysis. They were found to be involved in HIF-1 pathways known to affect energy metabolism and apoptosis. MiRNA-mRNA interaction pairs were detected from comparison of expression between the four different stages. The function annotation results also showed that many miRNA-mRNA interaction pairs were likely to be involved in regulating hypoxia stress. As a unique resource for gene expression and regulation during hypoxia and reoxygenation, this study could provide a starting point for further studies to better understand the genetic background of hypoxia stress.

Transcriptomic profiling in muscle and adipose tissue identifies genes related to growth and lipid deposition

Growth performance and meat quality are important traits for the pig industry and consumers. Adipose tissue is the main site at which fat storage and fatty acid synthesis occur. Therefore, we combined high-throughput transcriptomic sequencing in adipose and muscle tissues with the quantification of corresponding phenotypic features using seven Chinese indigenous pig breeds and one Western commercial breed (Yorkshire). We obtained data on 101 phenotypic traits, from which principal component analysis distinguished two groups: one associated with the Chinese breeds and one with Yorkshire. The numbers of differentially expressed genes between all Chinese breeds and Yorkshire were shown to be 673 and 1056 in adipose and muscle tissues, respectively. Functional enrichment analysis revealed that these genes are associated with biological functions and canonical pathways related to oxidoreductase activity, immune response, and metabolic process. Weighted gene coexpression network analysis found more coexpression modules significantly correlated with the measured phenotypic traits in adipose than in muscle, indicating that adipose regulates meat and carcass quality. Using the combination of differential expression, QTL information, gene significance, and module hub genes, we identified a large number of candidate genes potentially related to economically important traits in pig, which should help us improve meat production and quality.

Integrated analysis of mRNA-seq and miRNA-seq reveals the potential roles of sex-biased miRNA-mRNA pairs in gonad tissue of dark sleeper (Odontobutis potamophila)

Background

The dark sleeper (Odontobutis potamophila) is an important commercial fish species which shows a sexually dimorphic growth pattern. However, the lack of sex transcriptomic data is hindering further research and genetically selective breeding of the dark sleeper. In this study, integrated analysis of mRNA and miRNA was performed on gonad tissue to elucidate the molecular mechanisms of sex determination and differentiation in the dark sleeper.

Results

A total of 143 differentially expressed miRNAs and 16,540 differentially expressed genes were identified. Of these, 8103 mRNAs and 75 miRNAs were upregulated in testes, and 8437 mRNAs and 68 miRNAs were upregulated in ovaries. Integrated analysis of miRNA and mRNA expression profiles predicted more than 50,000 miRNA-mRNA interaction sites, and among them 27,583 negative miRNA-mRNA interactions. A number of sex related genes were targeted by sex-biased miRNAs. The relationship between 15 sex-biased genes and 15 sex-biased miRNAs verified by using qRT-PCR were described. Additionally, a number of SNPs were revealed through the transcriptome data.

Conclusions

The overall results of this study facilitate our understanding of the molecular mechanism underlying sex determination and differentiation and provide valuable genomic information for selective breeding of the dark sleeper.

Electronic supplementary material

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

MicroRNA expression profiles differ between primary myofiber of lean and obese pig breeds

MicroRNAs (miRNAs) are non-coding small miRNAs ~22 nucleotides in length and play a vital role in muscle development by binding to messenger RNAs (mRNAs). Large White (LW, a lean type pig) and Meishan pigs (MS, a Chinese indigenous obese breed) have significant postnatal phenotype differences in growth rate, muscle mass and meat quality, and these differences are programmed during prenatal muscle development. Little research shed light directly on the miRNA transcriptome difference in prenatal muscles between these two distinct pig breeds. Myofiber phenotypes of LW and MS were measured at developmental stages of 35, 55 and 90 days post-conception (dpc), which revealed that the myogenesis process is more intense in MS than in LW at 35 dpc. To investigate the role of miRNAs involved in regulating muscle development at earlier stages of myogenesis and decipher the miRNAs transcriptome difference between LW and MS, here, the miRNAomes of longissimus dorsi muscle collected at 35 dpc from female LW and MS were analyzed by deep sequencing. Overall, 1147 unique miRNAs comprising 434 known miRNAs, 239 conserved miRNAs and 474 candidate miRNAs were identified. Expression analysis of the 10 most abundant miRNAs in every library indicated that functional miRNAome may be a small amount and tend to be greater expressed. These sets of miRNA may play house keeping roles that were involved in myogenesis. A total of 87 miRNAs were significantly differentially expressed between LW and MS (reads > 1000, P < 0.05). Gene ontology (GO) and KEGG pathway enrichment analysis revealed that the differentially expressed miRNAs (DE miRNAs) were associated mainly with muscle contraction, WNT, mTOR, and MAPK signaling pathways. Some myogenesis related miRNAs (miR-133, miR-1, miR-206 and miR-148a) are highly abundant in MS, while other miRNAs (let-7 family, miR-214, miR-181) highly expressed in LW. In addition, the expression patterns of miRNAs (miR-1, -133, -206) at three prenatal stages (35, 55 and 90 dpc) were determined using qRT-PCR. Notably, ssc-miR-133 was significantly more highly expressed in LW pigs skeletal muscle at all prenatal stages compared with its expression in LW pigs skeletal muscle. Taken together, the main functional miRNAs during muscle development are different between lean and obese pig breeds. The present study adds new information to existing data on porcine miRNAs and will be helpful to investigate the dominant (main functional) muscle-related miRNAs sets in different pig breeds.

Modulation of Host miRNAs Transcriptome in Lung and Spleen of Peste des Petits Ruminants Virus Infected Sheep and Goats

Peste des petits ruminants (PPR) is one of the highly contagious viral disease, characterized by fever, sore mouth, conjunctivitis, gastroenteritis, and pneumonia, primarily affecting sheep and goats. Reports suggested variable host response in goats and sheep and this host response vis-a-vis the expression of microRNAs (miRNAs) has not been investigated. Here, miRNAs were sequenced and proteomics data were generated to identify the role of differentially expressed miRNA (DEmiRNA) in PPR virus (PPRV) infected lung and spleen tissues of sheep and goats. In lungs, 67 and 37 DEmiRNAs have been identified in goats and sheep, respectively. Similarly, in spleen, 50 and 56 DEmiRNAs were identified in goats and sheep, respectively. A total of 20 and 11 miRNAs were found to be common differentially expressed in both the species in PPRV infected spleen and lung, respectively. Six DEmiRNAs—miR-21-3p, miR-1246, miR-27a-5p, miR-760-3p, miR-320a, and miR-363 were selected based on their role in viral infections, apoptosis, and fold change. The target prediction analysis of these six selected DEmiRNAs from the proteome data generated, revealed involvement of more number of genes in lung and spleen of goats than in sheep. On gene ontology analysis of host target genes these DEmiRNAs were found to regulate several immune response signaling pathways. It was observed that the pathways viz. T cell receptor signaling, Rap1 signaling, Toll-like receptor signaling, and B cell receptor signaling governed by DEmiRNAs were more perturbed in goats than in sheep. The data suggests that PPRV-induced miR-21-3p, miR-320a, and miR-363 might act cooperatively to enhance viral pathogenesis in the lung and spleen of sheep by downregulating several immune response genes. The study gives an important insight into the molecular pathogenesis of PPR by identifying that the PPRV—Izatnagar/94 isolate elicits a strong host response in goats than in sheep.

Comparative analysis of DNA methylome and transcriptome of skeletal muscle in lean-, obese-, and mini-type pigs

DNA methylation plays a pivotal role in biological processes by affecting gene expression. However, how DNA methylation mediates phenotype difference of skeletal muscle between lean-, obese-, and mini-type pigs remains unclear. We systematically carried out comparative analysis of skeletal muscle by integrating analysis of genome-wide DNA methylation, mRNA, lncRNA and miRNA profiles in three different pig breeds (obese-type Tongcheng, lean-type Landrace, and mini-type Wuzhishan pigs). We found that the differentially methylated genes (DMGs) were significantly associated with lipid metabolism, oxidative stress and muscle development. Among the identified DMGs, 253 genes were related to body-size and obesity. A set of lncRNAs and mRNAs including UCP3, FHL1, ANK1, HDAC4, and HDAC5 exhibited inversely changed DNA methylation and expression level; these genes were associated with oxidation reduction, fatty acid metabolism and cell proliferation. Gene regulatory networks involved in phenotypic variation of skeletal muscle were related to lipid metabolism, cellular movement, skeletal muscle development, and the p38 MAPK signaling pathway. DNA methylation potentially influences the propensity for obesity and body size by affecting gene expression in skeletal muscle. Our findings provide an abundant information of epigenome and transcriptome that will be useful for animal breeding and biomedical research.

MicroRNA-mRNA regulatory networking fine-tunes the porcine muscle fiber type, muscular mitochondrial respiratory and metabolic enzyme activities

Background

MicroRNAs (miRNAs) are small non-coding RNAs that play critical roles in diverse biological processes via regulation of gene expression including in skeletal muscles. In the current study, miRNA expression profile was investigated in longissimus muscle biopsies of malignant hyperthermia syndrome-negative Duroc and Pietrain pigs with distinct muscle metabolic properties in order to explore the regulatory role of miRNAs related to mitochondrial respiratory activity and metabolic enzyme activity in skeletal muscle.

Results

A comparative analysis of the miRNA expression profile between Duroc and Pietrain pigs was performed, followed by integration with mRNA profiles based on their pairwise correlation and computational target prediction. The identified target genes were enriched in protein ubiquitination pathway, stem cell pluripotency and geranylgeranyl diphosphate biosynthesis, as well as skeletal and muscular system development. Next, we analyzed the correlation between individual miRNAs and phenotypical traits including muscle fiber type, mitochondrial respiratory activity, metabolic enzyme activity and adenosine phosphate concentrations, and constructed the regulatory miRNA-mRNA networks associated with energy metabolism. It is noteworthy that miR-25 targeting BMPR2 and IRS1, miR-363 targeting USP24, miR-28 targeting HECW2 and miR-210 targeting ATP5I, ME3, MTCH1 and CPT2 were highly associated with slow-twitch oxidative fibers, fast-twitch oxidative fibers, ADP and ATP concentration suggesting an essential role of the miRNA-mRNA regulatory networking in modulating the mitochondrial energy expenditure in the porcine muscle. In the identified miRNA-mRNA network, a tight relationship between mitochondrial and ubiquitin proteasome system at the level of gene expression was observed. It revealed a link between these two systems contributing to energy metabolism of skeletal muscle under physiological conditions.

Conclusions

We assembled miRNA-mRNA regulatory networks based on divergent muscle properties between different pig breeds and further with the correlation analysis of expressed genes and phenotypic measurements. These complex networks relate to muscle fiber type, metabolic enzyme activity and ATP production and may contribute to divergent muscle phenotypes by fine-tuning the expression of genes. Altogether, the results provide an insight into a regulatory role of miRNAs in muscular energy metabolisms and may have an implication on meat quality and production.

Electronic supplementary material

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

Genome-wide analysis of DNA methylation in obese, lean, and miniature pig breeds

DNA methylation is a crucial epigenetic modification involved in diverse biological processes. There is significant phenotypic variance between Chinese indigenous and western pig breeds. Here, we surveyed the genome-wide DNA methylation profiles of blood leukocytes from three pig breeds (Tongcheng, Landrace, and Wuzhishan) by methylated DNA immunoprecipitation sequencing. The results showed that DNA methylation was enriched in gene body regions and repetitive sequences. LINE/L1 and SINE/tRNA-Glu were the predominant methylated repeats in pigs. The methylation level in the gene body regions was higher than in the 5' and 3' flanking regions of genes. About 15% of CpG islands were methylated in the pig genomes. Additionally, 2,807, 2,969, and 5,547 differentially methylated genes (DMGs) were identified in the Tongcheng vs. Landrace, Tongcheng vs. Wuzhishan, and Landrace vs. Wuzhishan comparisons, respectively. A total of 868 DMGs were shared by the three contrasts. The DMGs were significantly enriched in development- and metabolism-related biological processes and pathways. Finally, we identified 32 candidate DMGs associated with phenotype variance in pigs. Our research provides a DNA methylome resource for pigs and furthers understanding of epigenetically regulated phenotype variance in mammals.

Transcriptome analysis revealed chimeric RNAs, single nucleotide polymorphisms and allele-specific expression in porcine prenatal skeletal muscle

Prenatal skeletal muscle development genetically determines postnatal muscle characteristics such as growth and meat quality in pigs. However, the molecular mechanisms underlying prenatal skeletal muscle development remain unclear. Here, we performed the first genome-wide analysis of chimeric RNAs, single nuclear polymorphisms (SNPs) and allele-specific expression (ASE) in prenatal skeletal muscle in pigs. We identified 14,810 protein coding genes and 163 high-confidence chimeric RNAs expressed in prenatal skeletal muscle. More than 94.5% of the chimeric RNAs obeyed the canonical GT/AG splice rule and were trans-splicing events. Ten and two RNAs were aligned to human and mouse chimeric transcripts, respectively. We detected 106,457 high-quality SNPs (6,955 novel), which were mostly (89.09%) located within QTLs for production traits. The high proportion of non-exonic SNPs revealed the incomplete annotation status of the current swine reference genome. ASE analysis revealed that 11,300 heterozygous SNPs showed allelic imbalance, whereas 131 ASE variants were located in the chimeric RNAs. Moreover, 4 ASE variants were associated with various economically relevant traits of pigs. Taken together, our data provide a source for studies of chimeric RNAs and biomarkers for pig breeding, while illuminating the complex transcriptional events underlying prenatal skeletal muscle development in mammals.

Analysis of miRNAs and their target genes associated with lipid metabolism in duck liver

Fat character is an important index in duck culture that linked to local flavor, feed cost and fat intake for costumers. Since the regulation networks in duck lipid metabolism had not been reported very clearly, we aimed to explore the potential miRNA-mRNA pairs and their regulatory roles in duck lipid metabolism. Here, Cherry-Valley ducks were selected and treated with/without 5% oil added in feed for 2 weeks, and then fat content determination was performed on. The data showed that the fat contents and the fatty acid ratios of C17:1 and C18:2 were up-regulated in livers of oil-added ducks, while the C12:0 ratio was down-regulated. Then 21 differential miRNAs, including 10 novel miRNAs, were obtain from the livers by sequencing, and 73 target genes involved in lipid metabolic processes of these miRNAs were found, which constituted 316 miRNA-mRNA pairs. Two miRNA-mRNA pairs including one novel miRNA and one known miRNA, N-miR-16020-FASN and gga-miR-144-ELOVL6, were selected to validate the miRNA-mRNA negative relation. And the results showed that N-mir-16020 and gga-miR-144 could respectively bind the 3′-UTRs of FASN and ELOVL6 to control their expressions. This study provides new sights and useful information for future research on regulation network in duck lipid metabolism.