Search

ST6GAL2, a member of the sialoglycosyltransferase family, primarily localizes within the cellular Golgi apparatus. However, the role of the ST6GAL2 gene in skeletal muscle growth and development remains elusive. In this study, the impact of the ST6GAL2 gene on the proliferation, differentiation, and apoptosis of primary chicken myoblasts at the cellular level was investigated. Quantitative fluorescent PCR was used to measure the expression levels of genes. Subsequently, using gene knockout mice, we assessed its effects on skeletal muscle growth and development in vivo. Our findings reveal that the ST6GAL2 gene promotes the expression of cell cycle and proliferation-related genes, including CCNB2 and PCNA, and apoptosis-related genes, such as Fas and Caspase-9. At the individual level, double knockout of ST6GAL2 inhibited the formation of both fast and slow muscle fibers in the quadriceps, extensor digitorum longus, and tibial anterior muscle, while promoting their formation in the gastrocnemius and soleus. These results collectively demonstrate that the ST6GAL2 gene facilitates the proliferation, apoptosis, and fusion processes of primary chicken myoblasts. Additionally, it promotes the enlargement of cross-sectional muscle fiber areas and regulates the formation of fast and slow muscle fibers at the individual level, albeit inhibiting muscle fusion. This study provides valuable insights into the role of the ST6GAL2 gene in promoting proliferation of skeletal muscle.

IntroductionEchinococcosis is a chronic zoonotic disease caused by tapeworms of the genus Echinococcus. The World Health Organization (WHO) has identified encapsulated disease as one of 17 neglected diseases to be controlled or eliminated by 2050. There is no accurate, early, non-invasive molecular diagnostic method to detect echinococcosis. The feasibility of circulating free DNA as a diagnostic method for echinococcosis has yielded inconclusive results in a number of published studies. However, there has been no systematic evaluation to date assessing the overall performance of these assays. We report here the first meta-analysis assessing the diagnostic accuracy of cfDNA in plasma, serum, and urine for echinococcosis.MethodsWe systematically searched PubMed, Embase, Cochrane Library, China National Knowledge Infrastructure (CNKI), and WeiPu databases up to 17 January 2024, for relevant studies. All analyses were performed using RevMan 5.3, Meta-DiSc 1.4, Stata 17.0, and R 4.3.1 software. The sensitivity, specificity, and other accuracy indicators of circulating free DNA for the diagnosis of echinococcosis were summarized. Subgroup analyses and meta-regression were performed to identify sources of heterogeneity.ResultsA total of 7 studies included 218 patients with echinococcosis and 214 controls (156 healthy controls, 32 other disease controls (non-hydatid patients), and 26 non-study-targeted echinococcosis controls were included). Summary estimates of the diagnostic accuracy of cfDNA in the diagnosis of echinococcosis were as follows: sensitivity (SEN) of 0.51 (95% CI: 0.45–0.56); specificity (SPE) of 0.99 (95% CI: 0.97–0.99); positive likelihood ratio (PLR) of 11.82 (95% CI: 6.74–20.74); negative likelihood ratio (NLR) of 0.57 (95% CI: 0.41–0.80); diagnostic ratio (DOR) of 36.63 (95% CI: 13.75–97.59); and area under the curve (AUC) value of 0.98 (95% CI: 0.96–1.00).ConclusionExisting evidence indicates that the combined specificity of circulating cfDNA for echinococcosis is high. However, the combined sensitivity performance is unsatisfactory due to significant inter-study heterogeneity. To strengthen the validity and accuracy of our findings, further large-scale prospective studies are required.Systematic review registrationThe systematic review was registered in the International Prospective Register of Systematic Reviews PROSPERO [CRD42023454158]. https://www.crd.york.ac.uk/PROSPERO/.

Commensal microbial–host interaction is crucial for host metabolism, growth, development, and immunity. However, research on microbial–host immunity in large animal models has been limited. This study was conducted to investigate the effects of the commensal microbiota on immune function in two model groups: germ-free (GF) and specific-pathogen-free (SPF) piglets. The weight and organ index of the spleen of the GF piglet were larger than those in the SPF piglet (P < 0.05). The histological structure of the red pulp area and mean area of germinal centers were larger in the SPF piglet than in the GF piglet (P < 0.05), whereas the areas of staining of B cells and T cells in the spleen and mesenteric lymph nodes (MLNs) were lower in the GF piglet (P < 0.05). We identified immune-related genes in the spleen and MLNs using RNA sequencing, and used real-time quantitative PCR to analyze the expression of core genes identified in gene set enrichment analysis. The expression levels of genes in the transforming growth factor-β/SMAD3 signaling pathway, Toll-like receptor 2/MyD88/nuclear factor-κB signaling pathway, and pro-inflammatory factor genes IL-6 and TNF-α in the spleen and MLNs were higher in the SPF piglet and in splenic lymphocytes compared with those in the GF and control group, respectively, under treatment with acetic acid, propionic acid, butyric acid, lipopolysaccharide (LPS), or concanavalin A (ConA). The abundances of plasma cells, CD8++ T cells, follicular helper T cells, and resting natural killer cells in the spleen and MLNs were significantly greater in the SPF piglet than in the GF piglet (P < 0.05). In conclusion, the commensal microbiota influenced the immune tissue structure, abundances of immune cells, and expression of immune-related pathways, indicating the importance of the commensal microbiota for spleen and MLNs development and function. In our study, GF piglet was used as the research model, eliminating the interference of microbiota in the experiment, and providing a suitable and efficient large animal research model for exploring the mechanism of “microbial-host” interactions.

Pigs are widely used as animal models in various studies related to humans. The interaction between the gut microbiota and the host has significant effects on the host’s health and disease status. However, although there have been many studies investigating the pig gut microbiota, the findings have been inconsistent due to variations in rearing conditions. Interactions between the gut microbiota and host have not been fully explored in pigs. Specific pathogen-free (SPF) pigs are ideal non-primate large animals to study the interactions between the gut microbiota and the host. In this study, we performed high-throughput sequencing analysis of the gut microbiota and the gut tissue transcriptome of six SPF pigs to provide a systematic understanding of the composition, function, and spatial distribution of gut microbiota in SPF pigs. We identified significant differences in microbial diversity and functionality among different gastrointestinal tract sites. Metagenomics data analysis revealed significant differences in alpha diversity and beta diversity of microbiota in different gastrointestinal sites of SPF pigs. Additionally, transcriptomic data indicated significant differences in gene expression as well as KEGG and GO functional enrichment between the small intestine and large intestine. Furthermore, by combining microbial metagenomics and host transcriptomics analyses, specific correlations were found between gut microbiota and host genes. These included a negative correlation between the TCN1 gene and Prevotella dentalis, possibly related to bacterial metabolic pathways involving vitamin B12, and a positive correlation between the BDH1 gene and Roseburia hominis, possibly because both are involved in fatty acid metabolism. These findings lay the groundwork for further exploration of the co-evolution between the microbiota and the host, specifically in relation to nutrition, metabolism, and immunity. In conclusion, we have elucidated the diversity of the gut microbiota in SPF pigs and conducted a detailed investigation into the interactions between the gut microbiota and host gene expression. These results contribute to our understanding of the intricate dynamics between the gut microbiota and the host, offering important references for advancements in life science research, bioproduct production, and sustainable development in animal husbandry.

IntroductionPigs are often used to study the intestinal development of newborns, particularly as preterm pig models that mimic the intestinal growth of human preterm infants. Neonatology’s study of delivery mode’s impact on neonatal development is crucial.MethodsWe established 14 newborn pigs delivered via cesarean sections (C-section, at 113 days of gestational age, CS group) and 8 naturally born pigs were used as controls (at 114 days of gestational age, NF group). The impact of two alternative delivery procedures (C-section and natural birth) on the levels of short-chain fatty acids (SCFAs) and organic acids in the hepatic and intestines of newborn pigs were compared using metabolomics. The underlying molecular pathways are examined at the “protein-metabolite” level by integrating proteomic data.ResultsThe findings demonstrated that the mode of delivery changed the metabolism of SCFAs in newborn pigs, perhaps by affecting the physiology levels of cyclic intermediates such as lactate and malate in the pyruvate metabolic pathway. Additionally, by participating in the fatty acid metabolism pathway, two distinct proteins (FASN and HSD17B4) may impact the physiological concentration of these tiny metabolites.DiscussionIn conclusion, this study provided reliable animal model data for understanding the physiological SCFA metabolic information and its affecting mechanism of large-gestational age preterm infants.

Background: Angiogenesis is essential for various physiological and pathological processes, such as embryonic development and cancer cell proliferation, migration, and invasion. Long noncoding RNAs (lncRNAs) play pivotal roles in normal homeostasis and disease processes by regulating gene expression through various mechanisms, including competing endogenous RNAs (ceRNAs) of target microRNAs (miRNAs). The lncRNA MYU is known to promote prostate cancer proliferation via the miR-184/c-Myc regulatory axis and to be upregulated in vascular endothelial cells under hypoxic conditions, which often occurs in solid tumors. In the present study, we investigated whether MYU might affect cancer growth by regulating angiogenesis in vascular endothelial cells under hypoxia. Methods: The expression of MYU-regulated miR-23a-3p and interleukin-8 (IL-8) in HUVEC cell lines was examined using qRT-PCR. The CCK-8 assay, EdU assay, wound-healing assay, and tube-formation assay were used to assess the effects of MYU on cell proliferation, migration, and tube formation of HUVEC cells in vitro. The dual-luciferase reporter assay was performed to examine the effects of miR-23a-3p on MYU and IL-8 expression. Results: We found that the overexpression of MYU and knockdown of miR-23a-3p in human umbilical vein endothelial cells (HUVECs) under hypoxia promoted cell proliferation, migration, and tube formation. Mechanistically, MYU was shown to bind competitively to miR-23a-3p, thereby preventing miR-23a-3p binding to the 3′ untranslated region of IL-8 mRNA. In turn, increased production of pro-angiogenic IL-8 promoted HUVEC proliferation, migration, and tube formation under hypoxia. Conclusion: This study identified a new role for lncRNA MYU as a ceRNA for miR-23a-3p and uncovered a novel MYU–miR-23a-3p–IL-8 regulatory axis for angiogenesis. MYU and/or miR-23a-3p may thus represent new targets for the treatment of hypoxia-related diseases by promoting angiogenesis.

The normal growth and development of skeletal muscle is essential for the health of the body. The regulation of skeletal muscle by intestinal microorganisms and their metabolites has been continuously demonstrated. Acetate is the predominant short-chain fatty acids synthesized by gut microbiota through the fermentation of dietary fiber; however, the underlying molecular mechanisms governing the interaction between acetate and skeletal muscle during the rapid growth stage remains to be further elucidated. Herein, specific pathogen-free (SPF) mice, germ-free (GF) mice, and germ-free mice supplemented with sodium acetate (GS) were used to evaluate the effects of acetate on the skeletal muscle growth and development of young mice with gut microbiota deficiency. We found that the concentration of serum acetate, body mass gain, succinate dehydrogenase activity, and expression of the myogenesis maker gene of skeletal muscle in the GS group were higher than those in the GF group, following sodium acetate supplementation. Furthermore, the transcriptome analysis revealed that acetate activated the biological processes that regulate skeletal muscle growth and development in the GF group, which are otherwise inhibited due to a gut microbiota deficiency. The in vitro experiment showed that acetate up-regulated Gm16062 to promote skeletal muscle cell differentiation. Overall, our findings proved that acetate promotes skeletal muscle growth and development in young mice via increasing Gm16062 expression.

Abstract Using an adult female miniature pig model with diet-induced weight gain/weight loss, we investigated the regulatory mechanisms of three-dimensional (3D) genome architecture in adipose tissues (ATs) associated with obesity. We generated 249 high-resolution in situ Hi-C chromatin contact maps of subcutaneous AT and three visceral ATs, analyzing transcriptomic and chromatin architectural changes under different nutritional treatments. We find that chromatin architecture remodeling underpins transcriptomic divergence in ATs, potentially linked to metabolic risks in obesity development. Analysis of chromatin architecture among subcutaneous ATs of different mammals suggests the presence of transcriptional regulatory divergence that could explain phenotypic, physiological, and functional differences in ATs. Regulatory element conservation analysis in pigs and humans reveals similarities in the regulatory circuitry of genes responsible for the obesity phenotype and identified non-conserved elements in species-specific gene sets that underpin AT specialization. This work provides a data-rich tool for discovering obesity-related regulatory elements in humans and pigs.

Abstract Background Skeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined. Results In this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers. Conclusions This study provides new insight into the chromatin interactions that regulate Myh genes expression.

Abstract Liver development is a complex process that is regulated by a series of signaling pathways. Three-dimensional (3D) chromatin architecture plays an important role in transcriptional regulation; nonetheless, its dynamics and role in the rapid transition of core liver functions during development and obesity-induced metabolic stress remain largely unexplored. To investigate the dynamic chromatin architecture during liver development and under metabolic stress, we generated high-resolution maps of chromatin architecture for porcine livers across six major developmental stages (from embryonic day 38 to the adult stage) and under a high-fat diet-induced obesity. The characteristically loose chromatin architecture supports a highly plastic genome organization during early liver development, which fundamentally contributes to the rapid functional transitions in the liver after birth. We reveal the multi-scale reorganization of chromatin architecture and its influence on transcriptional regulation of critical signaling processes during liver development, and show its close association with transition in hepatic functions (i.e., from hematopoiesis in the fetus to metabolism and immunity after birth). The limited changes in chromatin structure help explain the observed metabolic adaptation to excessive energy intake in pigs. These results provide a global overview of chromatin architecture dynamics associated with the transition of physiological liver functions between prenatal development and postnatal maturation, and a foundational resource that allows for future in-depth functional characterization.

Abstract Background The three-dimensional (3D) architecture of the genome has a highly ordered and hierarchical nature, which influences the regulation of essential nuclear processes at the basis of gene expression, such as gene transcription. While the hierarchical organization of heterochromatin and euchromatin can underlie differences in gene expression that determine evolutionary differences among species, the way 3D genome architecture is affected by evolutionary forces within major lineages remains unclear. Here, we report a comprehensive comparison of 3D genomes, using high resolution Hi-C data in fibroblast cells of fish, chickens, and 10 mammalian species. Results This analysis shows a correlation between genome size and chromosome length that affects chromosome territory (CT) organization in the upper hierarchy of genome architecture, whereas lower hierarchical features, including local transcriptional availability of DNA, are selected through the evolution of vertebrates. Furthermore, conservation of topologically associating domains (TADs) appears strongly associated with the modularity of expression profiles across species. Additionally, LINE and SINE transposable elements likely contribute to heterochromatin and euchromatin organization, respectively, during the evolution of genome architecture. Conclusions Our analysis uncovers organizational features that appear to determine the conservation and transcriptional regulation of functional genes across species. These findings can guide ongoing investigations of genome evolution by extending our understanding of the mechanisms shaping genome architecture.

The gastrointestinal tract contains a complex microbial community. Peyer’s patches (PPs) play an important role in inducing mucosal immune responses in the gastrointestinal tract. However, little is known about the effect of commensal microbiota on the host’s PPs. Here, we analyzed the phenotypic-to-transcriptome changes in the intestine PPs of specific pathogen-free (SPF) and germ-free (GF) piglets (fed in an environment with and without commensal microbiota, respectively) to elucidate the role of commensal microbiota in host intestine mucosal immunity. Analyses of anatomical and histological characteristics showed that commensal microbiota deficiency led to PP hypoplasia, especially regarding B and T cells. A total of 12,444 mRNAs were expressed in 12 libraries; 2,156 and 425 differentially expressed (DE) mRNAs were detected in the jejunal PP (JPP) and ileal PP (IPP), respectively (SPF vs. GF). The shared DE mRNAs of the JPP and IPP were mainly involved in basic physiological and metabolic processes, while the specific DE mRNAs were enriched in regulating immune cells in the JPP and microbial responses and cellular immunity in the IPP. Commensal microbiota significantly modulated the expression of genes related to B-cell functions, including activation, proliferation, differentiation, apoptosis, receptor signaling, germinal center formation, and IgA isotype class switching, particularly in the JPP. TLR4 pathway-related genes were induced in response to microbial colonization and in LPS/SCFA-treated B cells. We also detected 69 and 21 DE lncRNAs in the JPP and IPP, respectively, and four one-to-one lncRNA-mRNA pairs were identified. These findings might represent key regulatory axes for host intestine mucosal immunity development during microbial colonization. Overall, the findings of this study revealed that commensal microbiota modulated phenotypic characteristics and gene expression in the piglet intestine PPs and underscored the importance of early microbial colonization for host mucosal immunity development.

Objective Investigation of muscle growth at different developmental stages is an appropriate strategy for studying the mechanisms underlying muscle development and differences in phenotypes. In particular, the muscle development mechanisms and the difference between the fastest and slowest growth rates. Methods In this study, we used a growth curve model to fit the growth inflection point (IP) of QingYu pigs and compared differences in the long non-coding RNA (lncRNA) transcriptome of muscle both at the growth IP and plateau phase (PP). Results The growth curve of the QingYu pig had a good fit (R2 = 0.974) relative to a typical S-curve and reached the IP at day 177.96. At the PP, marbling, intramuscular fat, and monounsaturated fatty acids had increased significantly and the percentage of lean muscle and polyunsaturated fatty acids had decreased. A total of 1,199 mRNAs and 62 lncRNAs were differentially expressed at the IP compared with the PP. Additional to gene ontology and Kyoto encyclopedia of genes and genomes pathway analyses, these differentially expressed protein coding genes were principally related to muscle growth and lipid metabolism. Conclusion Our results suggest that the identified differentially expressed lncRNAs, could play roles in muscle growth, fat deposition and regulation of fatty acid composition at the IP and PP.

Hypoxia is a common hallmark of healthy tissues in physiological states or chronically inflamed tissues in pathological states. Mammalian cells sense and adapt to hypoxia mainly through hypoxia-inducible factor (HIF) signaling. Many studies have shown that hypoxia and HIF signaling play an important regulatory role in development and function of innate immune cells and T cells, but their role in B cell biology is still controversial. B cells experience a complex life cycle (including hematopoietic stem cells, pro-B cells, pre-B cells, immature B cells, mature naïve B cells, activated B cells, plasma cells, and memory B cells), and the partial pressure of oxygen (PO2) in the corresponding developmental niche of stage-specific B cells is highly dynamic, which suggests that hypoxia and HIF signaling may play an indispensable role in B cell biology. Based on the fact that hypoxia niches exist in the B cell life cycle, this review focuses on recent discoveries about how hypoxia and HIF signaling regulate the development, metabolism, and function of B cells, to facilitate a deep understanding of the role of hypoxia in B cell-mediated adaptive immunity and to provide novel strategies for vaccine adjuvant research and the treatment of immunity-related or infectious diseases.

A comprehensive transcriptomic survey of the pig could enable mechanistic understanding of tissue specialization and accelerate its use as a biomedical model. Here the authors characterize four distinct transcript types in 31 adult pig tissues to dissect their distinct structural and transcriptional features and uncover transcriptomic variability related to tissue physiology.

The bursa of Fabricius (BF) is the critical humoral immune organ to birds, playing an essential role in B lymphocyte differentiation. However, unlike other poultries, surgical removal of pigeon BF did not limit humoral immune responsiveness. To investigate the expression profiles and the potential role of mRNA and long non-coding RNA (LncRNA) in squab BFs, transcriptome analysis was performed by RNA-Sequencing (RNA-Seq) over three developmental stages (1-day, 13 and 26 days old). We identified 13,072 mRNAs and 19,129 lncRNAs, of which 2,752 mRNAs and 1,515 lncRNAs were differential expressed (DE) in pigeon BFs over three developmental stages. Cluster analysis presented different expression patterns in DE mRNAs and lncRNAs. Functional enrichment analysis revealed that DE lncRNAs and mRNAs with distinct expression patterns might play crucial roles in the immune system process and tissue morphogenesis. In particular, some DE genes and lncRNAs with higher expression levels in 13D or 26D are related to lymphocyte activation and differentiation, adaptive immune response, positive regulation of immune response, leukocyte migration, etc. Protein-protein interaction (PPI) network and Molecular Complex Detection (MCODE) analysis sreened six significant modules containing 37 genes from immune-related DE gene cluster, which is closely linked in B cell activation, lymphocyte differentiation, B cell receptor signaling pathway, etc. Our study characterizes mRNA and lncRNA transcriptomic variability in pigeon BFs over different developmental stages and enhances understanding of the mechanisms underlying physiological functions of pigeon BF.

Skeletal muscle formation is an extremely important step in animal growth and development. Recent studies have found that TMEM8c (also known as Myomaker, MYMK), a muscle-specific transmembrane protein, can promote myoblast fusion and plays a key role in the normal development of skeletal muscle. However, the effect of Myomaker on porcine (Sus scrofa) myoblast fusion and the underlying regulatory mechanisms remain largely unknown. Therefore, in this study, we focused on the role and corresponding regulatory mechanism of the Myomaker gene during skeletal muscle development, cell differentiation, and muscle injury repair in pigs. We obtained the entire 3′ UTR sequence of porcine Myomaker using the 3′ RACE approach and found that miR-205 inhibited porcine myoblast fusion by targeting the 3′ UTR of Myomaker. In addition, based on a constructed porcine acute muscle injury model, we discovered that both the mRNA and protein expression of Myomaker were activated in the injured muscle, while miR-205 expression was significantly inhibited during skeletal muscle regeneration. The negative regulatory relationship between miR-205 and Myomaker was further confirmed in vivo. Taken together, the present study reveals that Myomaker plays a role during porcine myoblast fusion and skeletal muscle regeneration and demonstrates that miR-205 inhibits myoblast fusion through targeted regulation of the expression of Myomaker.

The liver is the central organ for metabolism and influence the growth and development of the animals. To date, little is known about the microRNA (miRNA) in pigeon livers, particularly in different developmental stages. A comprehensive investigation into miRNA transcriptomes in livers across 3 pigeon developmental stages (1, 14, 28 d old) and an adult stage (2 y old) was performed by small RNA sequencing. We identified 312 known miRNA, 433 conserved miRNA, and 192 novel miRNA in pigeon livers. A set of differentially expressed (DE) miRNA in livers were screened out during pigeon development. This set of miRNA might be involved in hepatospecific phenotype and liver development. A Short Time-series Expression Miner analysis indicated significant expression variations in DE miRNA during liver development of pigeons. These DE miRNA with different expression patterns might play essential roles in response to growth factor, cell morphogenesis, and gland development, etc. Protein–protein interaction network and Molecular Complex Detection analysis identified several vital target genes (e.g., TNRC6B, FRS2, PTCH1, etc.) of DE miRNA, which is closely linked in liver development and enriched in PI3K cascade and regulation of growth. Our results expanded the repertoire of pigeon miRNA and may be of help in better understanding the mechanism of squab's rapid development from the perspective of liver development.

The intestine is a tubular organ with multiple functions such as digestion absorption and immunity, but the functions of each intestinal segments are different. Intestinal regionalization is necessary for normal physiological function, but it also means the research results obtained at specific sites may not be applicable to other intestinal segments. In order to comprehensively describe the functional changes in the intestine, different intestinal segments and their contents (duodenum, jejunum, ileum, cecum, colon, and rectum) of guinea pigs were collected for RNA seq and 16S rRNA seq, respectively. The results showed differential genes of each intestinal segment mainly involve mucosa, digestion, absorption, and immunity. The gene sets related to fat, bill salts, vitamins, aggregates, amino acids, and water absorption were highly expressed in the small intestine, and the gene sets related to metal ions, nucleotides, and SCFAs were highly expressed in the large intestine. In terms of immunity, the CD8+ T, Th1, eosinophils, pDCs, and natural killer (NK) T cells in the small intestine showed higher scores than those in the large intestine, while the pattern-recognition receptor signaling pathway-related genes are highly expressed in the large intestine. In terms of microbial composition, Proteobacteria and Actinobacteria are abundant in the small intestine, while Firmicutes and Spirochaete are abundant in large intestine. The correlation analysis showed a high correlation between intestinal microorganisms and gene modules related to digestion and absorption. In addition, cross-species analysis showed the SCFA metabolism gene expression trends in human and rodent intestine were different. In conclusion, we analyzed the changes in substance transport, immune and microbial composition between different intestinal segments of guinea pigs, and explored the relationship between intestinal transcriptome and microorganisms, our research will provides a reference for subsequent intestinal-related research.

MicroRNAs (miRNAs) are a group of crucial regulators in the process of animal growth and development. However, little is known about the expression and function of miRNAs in pigeon muscles. To identify the miRNAs participating in the rapid development of pigeon pectoral muscles and quantitate their expression levels of pectoral muscles in different age stages, we performed miRNA transcriptome analysis in pigeon pectoral muscles by sequencing small RNAs over three different age stages (1-day old, 28 days old, and 2 years old). Dual-luciferase reporter assay was applied to validate the interaction between miRNA and its target gene. We identified 304 known miRNAs, 201 conserved miRNAs, and 86 novel miRNAs in pigeon pectoral muscles. 189 differentially expressed (DE) miRNAs were screened out during pigeon development. A short time-series expression miner (STEM) analysis indicated 89 DE miRNAs were significantly clustered in a progressively decreasing expression profile, and mainly enriched in biosynthesis-related GO categories and signaling pathways for MAPK and TGF-β. Dual-luciferase reporter assay indicated that a progressively down-regulated miRNA (miR-20b-5p) could directly target Krüppel-like factor 3 (KLF3) gene. To sum-up, our data expand the repertoire of pigeon miRNAs and enhance understanding of the mechanisms underlying rapid development in squabs.

Abstract Background Pigeon crop has the unique ability to produce a nutrient rich substance termed pigeon ‘milk’ (PM), which has functional resemblance with the mammalian milk. Previous researches have demonstrated that a large number of exosomes and exosomal miRNAs exist in mammalian milk, and many of them are associated with immunity, growth and development. However, to date, little is known about the exosomes and exosomal miRNAs in PM. Results In this study, we isolated the exosomes from PM and used small RNA sequencing to investigate the distribution and expression profiles of exosomal miRNAs. A total of 301 mature miRNAs including 248 conserved and 53 novel miRNAs were identified in five lactation stages i.e. 1d, 5d, 10d, 15d, and 20d. From these, four top 10 conserved miRNAs (cli-miR-21-5p, cli-miR-148a-3p, cli-miR-10a-5p and cli-miR-26a-5p) were co-expressed in all five stages. We speculate that these miRNAs may have important role in the biosynthesis and metabolism of PM. Moreover, similar to the mammalian milk, a significant proportion of immune and growth-related miRNAs were also present and enriched in PM exosomes. Furthermore, we also identified 41 orthologous miRNAs group (giving rise to 81 mature miRNA) commonly shared with PM, human, bovine and porcine breast milk. Additionally, functional enrichment analysis revealed the role of exosomal miRNAs in organ development and in growth-related pathways including the MAPK, Wnt and insulin pathways. Conclusions To sum-up, this comprehensive analysis will contribute to a better understanding of the underlying functions and regulatory mechanisms of PM in squabs.

Betaine is a natural compound present in commonly consumed foods and may have a potential role in the regulation of glucose and lipids metabolism. However, the underlying molecular mechanism of its action remains largely unknown. Here, we show that supplementation with betaine contributes to improved high-fat diet (HFD)-induced gut microbiota dysbiosis and increases anti-obesity strains such as Akkermansia muciniphila, Lactobacillus, and Bifidobacterium. In mice lacking gut microbiota, the functional role of betaine in preventing HFD-induced obesity, metabolic syndrome, and inactivation of brown adipose tissues are significantly reduced. Akkermansia muciniphila is an important regulator of betaine in improving microbiome ecology and increasing strains that produce short-chain fatty acids (SCFAs). Increasing two main members of SCFAs including acetate and butyrate can significantly regulate the levels of DNA methylation at host miR-378a promoter, thus preventing the development of obesity and glucose intolerance. However, these beneficial effects are partially abolished by Yin yang (YY1), a common target gene of the miR-378a family. Taken together, our findings demonstrate that betaine can improve obesity and associated MS via the gut microbiota-derived miR-378a/YY1 regulatory axis, and reveal a novel mechanism by which gut microbiota improve host health.

Local hypoxia has recently been reported to occur in the white adipose tissue (WAT) microenvironment during obesity. Adipocytes have a unique life cycle that reflects the different stages of adipogenesis in the WAT niche. Long non-coding RNAs (lncRNAs) play an important role in the cellular response to hypoxia. However, the differentially hypoxic responses of preadipocytes during adipogenesis and the potential role of lncRNAs in this process remain to be elucidated. Here, we evaluated the differentially hypoxic responses of primary hamster preadipocytes during adipogenesis and analyzed mRNA and lncRNA expression in same Ribo-Zero RNA-seq libraries. Hypoxia induced HIF-1α protein during adipogenesis and caused divergent changes of cell phenotypes. A total of 10,318 mRNAs were identified to be expressed in twenty libraries (five timepoints), and 3,198 differentially expressed mRNAs (DE mRNAs) were detected at five timepoints (hypoxia vs. normoxia). Functional enrichment analysis revealed the shared and specific hypoxia response pathways in the different stages of adipogenesis. Hypoxia differentially modulated the expression profile of adipose-associated genes, including adipokines, lipogenesis, lipolysis, hyperplasia, hypertrophy, inflammatory, and extracellular matrix. We also identified 4,296 lncRNAs that were expressed substantially and detected 1,431 DE lncRNAs at five timepoints. Two, 3, 5, 13, and 50 DE mRNAs at D0, D1, D3, D7, and D11, respectively, were highly correlated and locus-nearby DE lncRNAs and mainly involved in the cell cycle, vesicle-mediated transport, and mitochondrion organization. We identified 28 one-to-one lncRNA-mRNA pairs that might be closely related to adipocyte functions, such as ENSCGRT00015041780-Hilpda, TU2105-Cdsn, and TU17588-Ltbp3. These lncRNAs may represent the crucial regulation axis in the cellular response to hypoxia during adipogenesis. This study dissected the effects of hypoxia in the cell during adipogenesis, uncovered novel regulators potentially associated with WAT function, and may provide a new viewpoint for interpretation and treatment of obesity.

High-altitude acclimatization is a representative example of vertebrates’ acclimatization to harsh and extreme environments. Previous studies reported sufficient evidence for a molecular genetic basis of high-altitude acclimatization, and genomic patterns of genetic variation among populations and species have been widely elucidated in recent years. However, understanding of the miRNA role in high-altitude acclimatization have lagged behind, especially in non-model species. To investigate miRNA expression alterations of goats that were induced by high-altitude stress, we performed comparative miRNA transcriptome analysis on six hypoxia-sensitive tissues (heart, kidney, liver, lung, skeletal muscle, and spleen) in two goat populations from distinct altitudes (600 and 3000 m). We obtained the expression value of 1391 mature miRNAs and identified 138 differentially expressed (DE) miRNAs between high and low altitudes. Combined with tissue specificity analysis, we illustrated alterations of expression levels among altitudes and tissues, and found that there were coexisting tissue-specific and -conserved mechanisms for hypoxia acclimatization. Notably, the interplay between DE miRNA and DE target genes strongly indicated post-transcriptional regulation in the hypoxia inducible factor 1, insulin, and p53 signaling pathways, which might play significant roles in high-altitude acclimatization in domestic goats. It’s also worth noting that we experimentally confirmed miR-106a-5p to have a negative regulation effect on angiogenesis by directly targeting FLT-1. These results provide insight into the complicated miRNA expression patterns and regulatory mechanisms of high-altitude acclimatization in domestic goats.

Crossbreeding is widely used aimed at improving crossbred performance for poultry and livestock. Alleles that are specific to different purebreds will yield a large number of heterozygous single-nucleotide polymorphisms (SNPs) in crossbred individuals, which are supposed to have the power to alter gene function or regulate gene expression. For pork production, a classic three-way crossbreeding system of Duroc × (Landrace × Yorkshire) (DLY) is generally used to produce terminal crossbred pigs with stable and prominent performance. Nonetheless, little is known about the breed-of-origin effects from purebreds on DLY pigs. In this study, we first estimated the distribution of heterozygous SNPs in three kinds of three-way crossbred pigs via whole genome sequencing data originated from three purebreds. The result suggested that DLY is a more effective strategy for three-way crossbreeding as it could yield more stably inherited heterozygous SNPs. We then sequenced a DLY pig family and identified 95, 79, 132 and 42 allele-specific expression (ASE) genes in adipose, heart, liver and skeletal muscle, respectively. Principal component analysis and unrestricted clustering analyses revealed the tissue-specific pattern of ASE genes, indicating the potential roles of ASE genes for development of DLY pigs. In summary, our findings provided a lot of candidate SNP markers and ASE genes for DLY three-way crossbreeding system, which may be valuable for pig breeding and production in the future.

High-altitude inhospitable environments impose a formidable life challenge for the local animals. Training and exposure to high-altitude environments produce both distinct physiological and phenotypic characteristics. The mitochondrion, an organelle crucial for the energy production, plays an important role in hypoxia adaptation. In this study, we investigated the mitochondrial DNA (mtDNA) polymorphism and copy number variation between the population pairs from distinct altitudes across the multi-species. Higher mitochondrial DNA control region’s genetic diversity is conspicuous in high-altitude animals versus low-altitude relatives. We also found an accordant decrease of mtDNA copy number in most of the tissues from high-altitude animals. Compared to mammals, chickens have significantly distinct mitogenomic characteristics, and more significant changes in the skeletal muscle mtDNA copy number between high- and low-altitude individuals. Our study catches a snapshot of the biological similarities and differences in the mitochondrial high-altitude acclimation across the species.

Objective Hemicastration is a unilateral orchiectomy to remove an injured testis, which can induce hormonal changes and compensatory hypertrophy of the remaining testis, and may influence spermatogenesis. However, the underlying molecular mechanisms are poorly understood. Here, we investigated the impact of hemicastration on remaining testicular function. Methods Prepubertal mice (age 24 days) were hemicastrated, and their growth was monitored until they reached physical maturity (age 72 days). Subsequently, we determined testis DNA methylation patterns using reduced representation bisulfite sequencing of normal and hemicastrated mice. Moreover, we profiled the testicular gene expression patterns by RNA sequencing (RNA-seq) to examine whether methylation changes affected gene expression in hemicastrated mice. Results Hemicastration did not significantly affect growth or testosterone (p>0.05) compared with control. The genome-wide DNA methylation pattern of remaining testis suggested that substantial genes harbored differentially methylated regions (1,139) in gene bodies, which were enriched in process of protein binding and cell adhesion. Moreover, RNA-seq results indicated that 46 differentially expressed genes (DEGs) involved in meiotic cell cycle, synaptonemal complex assembly and spermatogenesis were upregulated in the hemicastration group, while 197 DEGs were downregulated, which were related to arachidonic acid metabolism. Integrative analysis revealed that proteasome 26S subunit ATPase 3 interacting protein gene, which encodes a protein crucial for homologous recombination in spermatocytes, exhibited promoter hypomethylation and higher expression level in hemicastrated mice. Conclusion Global profiling of DNA methylation and gene expression demonstrated that hemicastration-induced compensatory response maintained normal growth and testicular morphological structure in mice.

Abstract The physiological role of miRNAs is widely understood to include fine-tuning the post-transcriptional regulation of a wide array of biological processes. Extensive studies have indicated that exosomal miRNAs in the bodily fluids of various organisms can be transferred between living cells for the delivery of gene silencing signals. Here, we illustrated the expression characteristics of exosomal miRNAs in giant panda breast milk during distinct lactation periods and highlighted the enrichment of immune- and development-related endogenous miRNAs in colostral and mature giant panda milk. These miRNAs are stable, even under certain harsh conditions, via the protection of extracellular vesicles. These findings indicate that breast milk may facilitate the dietary intake of maternal miRNAs by infants for the regulation of postnatal development. We also detected exogenous plant miRNAs from the primary food source of the giant panda (bamboo) in the exosomes of giant panda breast milk that were associated with regulatory roles in basic metabolism and neuron development. This result suggested that dietary plant miRNAs are absorbed by host cells and subsequently secreted into bodily fluids as potential cross-kingdom regulators. In conclusion, exosomal miRNAs in giant panda breast milk may be crucial maternal regulators for the development of intrinsic ‘slink’ newborn cubs.

Abstract MicroRNAs are a class of small RNAs that are important in post-transcriptional gene regulation in animals and plants. These single-stranded molecules are widely distributed in organisms and influence fundamental biological processes. Interestingly, recent studies have reported that diet-derived plant miRNAs could regulate mammalian gene expression, and these studies have broadened our view of cross-kingdom communication. In the present study, we evaluated miRNA levels in cooked maize-containing chow diets, and found that plant miRNAs were resistant to the harsh cooking conditions to a certain extent. After feeding fresh maize to pigs (7 days), maize-derived miRNAs could be detected in porcine tissues and serum, and the authenticity of these plant miRNAs was confirmed by using oxidization reactions. Furthermore, in vivo and in vitro experiments demonstrated that dietary maize miRNAs could cross the gastrointestinal tract and enter the porcine bloodstream. In the porcine cells, we found that plant miRNAs are very likely to specifically target their endogenous porcine mRNAs and influence gene expression in a fashion similar to that of mammalian miRNAs. Our results indicate that maize-derived miRNAs can cross the gastrointestinal tract and present in pigs, and these exogenous miRNAs have the potential to regulate mammalian gene expression.

Abstract Background N 6 -methyladenosine (m6A) is the most prevalent internal form of modification in messenger RNA in higher eukaryotes and potential regulatory functions of reversible m6A methylation on mRNA have been revealed by mapping of m6A methylomes in several species. m6A modification in active gene regulation manifests itself as altered methylation profiles in a tissue-specific manner or in response to changing cellular or species living environment. However, up to date, there has no data on m6A porcine transcriptome-wide map and its potential biological roles in adipose deposition and muscle growth. Methods In this work, we used methylated RNA immunoprecipitation with next-generation sequencing (MeRIP-Seq) technique to acquire the first ever m6A porcine transcriptome-wide map. Transcriptomes of muscle and adipose tissues from three different pig breeds, the wild boar, Landrace, and Rongchang pig, were used to generate these maps. Results Our findings show that there were 5,872 and 2,826 m6A peaks respectively, in the porcine muscle and adipose tissue transcriptomes. Stop codons, 3′-untranslated regions, and coding regions were found to be mainly enriched for m6A peaks. Gene ontology analysis revealed that common m6A peaks in nuclear genes are associated with transcriptional factors, suggestive of a relationship between m6A mRNA methylation and nuclear genome transcription. Some genes showed tissue- and breed-differential methylation, and have novel biological functions. We also found a relationship between the m6A methylation extent and the transcript level, suggesting a regulatory role for m6A in gene expression. Conclusion This comprehensive map provides a solid basis for the determination of potential functional roles for RNA m6A modification in adipose deposition and muscle growth.

Mammalian mitochondrial biogenesis is a complex process involving mitochondrial proliferation and differentiation. Mitochondrial DNA transcription factor A (TFAM), which encodes a major component of a protein-mitochondrial DNA (mtDNA) complex, is regulated by peroxisome proliferator-activated receptor γ coactivator 1α (PGC1α). Testosterone is the primary male sex hormone and plays an increasingly important role in mammalian development through its interaction with androgen receptor (AR). However, the function of AR in mitochondrial biogenesis induced by testosterone deficiency has not been investigated. Here, we explored the molecular mechanism underlying the effect of testosterone deficiency on mitochondrial biogenesis using a Yorkshire boar model. Testosterone deficiency caused by castration induced changes in mtDNA copy numbers in various tissues, and AR showed the opposite tendency to that of mtDNA copy number, particularly in adipose tissues and muscle tissues. In addition, castration weakened the correlation of PGC1α and mtDNA copy number, while AR and TFAM showed a relatively high correlation in both control and castrated pigs. Furthermore, luciferase assays revealed that AR binds to potential AR elements in the TFAM promoter to promote TFAM expression. Taken together, testosterone may be involved in the pathway linking PGC1α to mitochondrial biogenesis through the interaction between AR and TFAM.

Epidemiological data have indicated that intrauterine growth retardation (IUGR) is a risk factor for the adult metabolic syndrome in pigs. However, the causative genetic mechanism leading to the phenotype in adulthood has not been well characterized. In the present study, both normal and IUGR adult pigs were used as models to survey the differences in global gene expression in livers through transcriptome sequencing. The transcriptome libraries generated 104.54 gb of data. In normal and IUGR pigs, 16,948 and 17,078 genes were expressed, respectively. A total of 1,322 differentially expressed genes (DEGs) were identified. Enrichment analysis of the DEGs revealed that the top overrepresented gene ontology (GO) terms and pathways were related to oxidoreductase activity, ATPase activity, amino catabolic process, glucose metabolism, and insulin signaling pathway. The increased gluconeogenesis (GNG) and decreased glycogen synthesis in the liver contributed to the glucose intolerance observed in IUGR. The reduced expression of insulin signaling genes (such as PI3K and AKT) indicated an elevated risk of diabetes in adulthood. Together, these findings provide a comprehensive understanding of the molecular mechanisms of adult IUGR pigs and valuable information for future studies of therapeutic intervention in IUGR metabolic syndrome.

Edible plant-derived exosome-like nanoparticles (EPDELNs) are novel naturally occurring plant ultrastructures that are structurally similar to exosomes. Many EPDELNs have anti-inflammatory properties. MicroRNAs (miRNAs) play a critical role in mediating physiological and pathological processes in animals and plants. Although miRNAs can be selectively encapsulated in extracellular vesicles, little is known about their expression and function in EPDELNs. In this study, we isolated nanovesicles from 11 edible fruits and vegetables and subjected the corresponding EPDELN small RNA libraries to Illumina sequencing. We identified a total of 418 miRNAs—32 to 127 per species—from the 11 EPDELN samples. Target prediction and functional analyses revealed that highly expressed miRNAs were closely associated with the inflammatory response and cancer-related pathways. The 418 miRNAs could be divided into three classes according to their EPDELN distributions: 26 “frequent” miRNAs (FMs), 39 “moderately present” miRNAs (MPMs), and 353 “rare” miRNAs (RMs). FMs were represented by fewer miRNA species than RMs but had a significantly higher cumulative expression level. Taken together, our in vitro results indicate that miRNAs in EPDELNs have the potential to regulate human mRNA.

The domestic pigs have been undergone intense selection pressures for these development of interested traits following domestication and modern breeding. This has altered many traits in most of pig breeds, such as growth rate, body weight, fertility, and immunity. Thus, the objectives of this study were to (1) detect these selection signatures and identify the candidate genes that show evidences of recent artificial selection at the level of whole genome, (2) be beneficial to understand the relationship between genomic structure and phenotypic diversity, and (3) highlight the key roles of these candidate genes in growth and development in the two breeds. The data consisted of total raw number of 345570 single nucleotide polymorphisms (SNPs) in 1200 individuals from the Chinese Landrace pigs (L, n = 600) and Yorkshire pigs (Y, n = 600). Based on these SNPs data, two complementary methods, population differentiation (Fst) and composite likelihood ratio test (CLR), were carried out to detect the selection signatures in this study. A total of 540 potential selection regions (50 kb) which contained 111 candidate genes were detected for Landrace-Yorkshire pair (L-Y) by Fst. In addition, 73 and 125 candidate genes were found for Landrace pigs and Yorkshire pigs by CLR test based on 321 and 628 potential selection regions, respectively. Some candidate genes are associated with important traits and signaling pathways including the ACACA, MECR, COL11A1, GHR, IGF1R, IGF2R, IFNG, and MTOR gene. The ACACA and MECR gene are related to fatty acid biosynthesis. The COL11A1 gene is essential for the development of the normal differentiation. The GHR, IGF1R, and IGF2R gene are significant candidate genes which play major roles in the growth and development in animals. The IFNG gene is associated with some aspects of immune response. The MTOR gene regulates many signaling pathways and signaling transduction pathway.

Genome-wide transcriptomic studies in humans and mice have become extensive and mature. However, a comprehensive and systematic understanding of protein-coding genes and long non-coding RNAs (lncRNAs) expressed during pig spleen development has not been achieved. LncRNAs are known to participate in regulatory networks for an array of biological processes. Here, we constructed 18 RNA libraries from developing fetal pig spleen (55 days before birth), postnatal pig spleens (0, 30, 180 days and 2 years after birth), and the samples from the 2-year-old Wild Boar. A total of 15,040 lncRNA transcripts were identified among these samples. We found that the temporal expression pattern of lncRNAs was more restricted than observed for protein-coding genes. Time-series analysis showed two large modules for protein-coding genes and lncRNAs. The up-regulated module was enriched for genes related to immune and inflammatory function, while the down-regulated module was enriched for cell proliferation processes such as cell division and DNA replication. Co-expression networks indicated the functional relatedness between protein-coding genes and lncRNAs, which were enriched for similar functions over the series of time points examined. We identified numerous differentially expressed protein-coding genes and lncRNAs in all five developmental stages. Notably, ceruloplasmin precursor (CP), a protein-coding gene participating in antioxidant and iron transport processes, was differentially expressed in all stages. This study provides the first catalog of the developing pig spleen, and contributes to a fuller understanding of the molecular mechanisms underpinning mammalian spleen development.

Extensive and in-depth investigations of high-altitude adaptation have been carried out at the level of morphology, anatomy, physiology and genomics, but few investigations focused on the roles of microRNA (miRNA) in high-altitude adaptation. We examined the differences in the miRNA transcriptomes of two representative hypoxia-sensitive tissues (heart and lung) between yak and cattle, two closely related species that live in high and low altitudes, respectively. In this study, we identified a total of 808 mature miRNAs, which corresponded to 715 pre-miRNAs in the two species. The further analysis revealed that both tissues showed relatively high correlation coefficient between yak and cattle, but a greater differentiation was present in lung than heart between the two species. In addition, miRNAs with significantly differentiated patterns of expression in two tissues exhibited co-operation effect in high altitude adaptation based on miRNA family and cluster. Functional analysis revealed that differentially expressed miRNAs were enriched in hypoxia-related pathways, such as the HIF-1α signaling pathway, the insulin signaling pathway, the PI3K-Akt signaling pathway, nucleotide excision repair, cell cycle, apoptosis and fatty acid metabolism, which indicated the important roles of miRNAs in high altitude adaptation. These results suggested the diverse degrees of miRNA transcriptome variation in different tissues between yak and cattle, and suggested extensive roles of miRNAs in high altitude adaptation.

MicroRNAs (miRNAs) are small non-coding RNAs that can regulate their target genes at the post-transcriptional level. Skeletal muscle comprises different fiber types that can be broadly classified as red, intermediate, and white. Recently, a set of miRNAs was found expressed in a fiber type-specific manner in red and white fiber types. However, an in-depth analysis of the miRNA transcriptome differences between all three fiber types has not been undertaken. Herein, we collected 15 porcine skeletal muscles from different anatomical locations, which were then clearly divided into red, white, and intermediate fiber type based on the ratios of myosin heavy chain isoforms. We further illustrated that three muscles, which typically represented each muscle fiber type (i.e., red: peroneal longus (PL), intermediate: psoas major muscle (PMM), white: longissimus dorsi muscle (LDM)), have distinct metabolic patterns of mitochondrial and glycolytic enzyme levels. Furthermore, we constructed small RNA libraries for PL, PMM, and LDM using a deep sequencing approach. Results showed that the differentially expressed miRNAs were mainly enriched in PL and played a vital role in myogenesis and energy metabolism. Overall, this comprehensive analysis will contribute to a better understanding of the miRNA regulatory mechanism that achieves the phenotypic diversity of skeletal muscles.

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.

Acute myocardial infarction (AMI) is an ischemic heart disease with high mortality worldwide. AMI triggers a hypoxic microenvironment and induces extensive myocardial injury, including autophagy and apoptosis. MiRNAs, which are a class of posttranscriptional regulators, have been shown to be involved in the development of ischemic heart diseases. We have previously reported that hypoxia significantly alters the miRNA transcriptome in rat cardiomyoblast cells (H9c2), including miR-27a-5p. In the present study, we further investigated the potential function of miR-27a-5p in the cardiomyocyte response to hypoxia, and showed that miR-27a-5p expression was downregulated in the H9c2 cells at different hypoxia-exposed timepoints and the myocardium of a rat AMI model. Follow-up experiments revealed that miR-27a-5p attenuated hypoxia-induced cardiomyocyte injury by regulating autophagy and apoptosis via Atg7, which partly elucidated the anti-hypoxic injury effects of miR-27a-5p. Taken together, this study shows that miR-27a-5p has a cardioprotective effect on hypoxia-induced H9c2 cell injury, suggesting it may be a novel target for the treatment of hypoxia-related heart diseases.

The gram-negative bacterium Actinobacillus pleuropneumoniae (APP) is an inhabitant of the porcine upper respiratory tract and the causative agent of porcine pleuropneumonia (PP). In recent years, knowledge about the proinflammatory cytokine and chemokine gene expression that occurs in lung and lymph node of the APP-infected swine has been advanced. However, systematic gene expression profiles on hilar nodes from pigs after infection with Actinobacillus pleuropneumoniae have not yet been reported. The transcriptional responses were studied in hilar nodes (HN) from swine experimentally infected with APP and the control groupusing Agilent Porcine Genechip, including 43,603 probe sets. 9,517 transcripts were identified as differentially expressed (DE) at the p ≤ 0.01 level by comparing the log2 (normalized signal) of the two groups named treatment group (TG) and controls (CG). Eight hundred and fifteen of these DE transcripts were annotated as pig genes in the GenBank database (DB). Two hundred and seventy-two biological process categories (BP), 75 cellular components and 171 molecular functions were substantially altered in the TG compared to CG. Many BP were involved in host immune responses (i.e., signaling, signal transmission, signal transduction, response to stimulus, oxidation reduction, response to stress, immune system process, signaling pathway, immune response, cell surface receptor linked signaling pathway). Seven DE gene pathways (VEGF signaling pathway, Long-term potentiation, Ribosome, Asthma, Allograft rejection, Type I diabetes mellitus and Cardiac muscle contraction) and statistically significant associations with host responses were affected. Many cytokines (including NRAS, PI3K, MAPK14, CaM, HSP27, protein phosphatase 3, catalytic subunit and alpha isoform), mediating the proliferation and migration of endothelial cells and promoting survival and vascular permeability, were activated in TG, whilst many immunomodulatory cytokines were suppressed. The significant changes in the expression patterns of the genes, GO terms, and pathways, led to a decrease of antigenic peptides with antigen presenting cells presented to T lymphocytes via the major histocompatibility complex, and alleviated immune response induced APP of HN. The immune response ability of HN in the APP-infected pigs was weakened; however, cell proliferation and migration ability was enhanced.

Porcine pleuropneumonia is a highly contagious respiratory disease that causes great economic losses worldwide. In this study, we aimed to explore the underlying relationship between infection and injury by investigation of the whole porcine genome expression profiles of swine lung tissues post-inoculated with experimentally Actinobacillus pleuropneumoniae. Expression profiling experiments of the control group and the treatment group were conducted using a commercially available Agilent Porcine Genechip including 43,603 probe sets. Microarray analysis was conducted on profiles of lung from challenged versus non-challenged swine. We found 11,929 transcripts, identified as differentially expressed at the p ≤0.01 level. There were 1188 genes annotated as swine genes in the GenBank Data Base. GO term analysis identified a total of 89 biological process categories, 82 cellular components and 182 molecular functions that were significantly affected, and at least 27 biological process categories that were related to the host immune response. Gene set enrichment analysis identified 13 pathways that were significantly associated with host response. Many proinflammatory-inflammatory cytokines were activated and involved in the regulation of the host defense response at the site of inflammation; while the cytokines involved in regulation of the host immune response were suppressed. All changes of genes and pathways of induced or repressed expression not only led to a decrease in antigenic peptides presented to T lymphocytes by APCs via the MHC and alleviated immune response injury induced by infection, but also stimulated stem cells to produce granulocytes (neutrophils, eosinophils, and basophils) and monocyte, and promote neutrophils and macrophages to phagocytose bacterial and foreign antigen at the site of inflammation. The defense function of swine infection with Actinobacillus pleuropneumoniae was improved, while its immune function was decreased.

Intermuscular adipose tissue is located between the muscle fiber bundles in skeletal muscles, and has similar metabolic features to visceral adipose tissue, which has been found to be related to a number of obesity-related diseases. Although various miRNAs are known to play crucial roles in adipose deposition and adipogenesis, the microRNA transcriptome of intermuscular adipose tissue has not, until now, been studied. Here, we sequenced the miRNA transcriptomes of porcine intermuscular adipose tissue by small RNA-sequencing and compared it to a representative subcutaneous adipose tissue. We found that the inflammation- and diabetes-related miRNAs were significantly enriched in the intermuscular rather than in the subcutaneous adipose tissue. A functional enrichment analysis of the genes predicted to be targeted by the enriched miRNAs also indicated that intermuscular adipose tissue was associated mainly with immune and inflammation responses. Our results suggest that the intermuscular adipose tissue should be recognized as a potential metabolic risk factor of obesity.

A better understanding of the control of lipogenesis is of critical importance for both human and animal physiology. This requires a better knowledge of the changes of gene expression during the process of adipose tissue development. Thus, the objective of the current study was to determine the effects of development on subcutaneous adipose tissue gene expression in growing and adult pigs. Here, we present a comprehensive investigation of mRNA transcriptomes in porcine subcutaneous adipose tissue across four developmental stages using digital gene expression profiling. We identified 3,274 differential expressed genes associated with oxidative stress, immune processes, apoptosis, energy metabolism, insulin stimulus, cell cycle, angiogenesis and translation. A set of universally abundant genes (ATP8, COX2, COX3, ND1, ND2, SCD and TUBA1B) was found across all four developmental stages. This set of genes may play important roles in lipogenesis and development. We also identified development-related gene expression patterns that are linked to the different adipose phenotypes. We showed that genes enriched in significantly up-regulated profiles were associated with phosphorylation and angiogenesis. In contrast, genes enriched in significantly down-regulated profiles were related to cell cycle and cytoskeleton organization, suggesting an important role for these biological processes in adipose growth and development. These results provide a resource for studying adipose development and promote the pig as a model organism for researching the development of human obesity, as well as being used in the pig industry.

MicroRNAs (miRNAs) play critical roles in many important biological processes, such as growth and development in mammals. Various studies of porcine muscle development have mainly focused on identifying miRNAs that are important for fetal and adult muscle development; however, little is known about the role of miRNAs in middle-aged muscle development. Here, we present a comprehensive investigation of miRNA transcriptomes across five porcine muscle development stages, including one prenatal and four postnatal stages. We identified 404 known porcine miRNAs, 118 novel miRNAs, and 101 miRNAs that are conserved in other mammals. A set of universally abundant miRNAs was found across the distinct muscle development stages. This set of miRNAs may play important housekeeping roles that are involved in myogenesis. A short time-series expression miner analysis indicated significant variations in miRNA expression across distinct muscle development stages. We also found enhanced differentiation- and morphogenesis-related miRNA levels in the embryonic stage; conversely, apoptosis-related miRNA levels increased relatively later in muscle development. These results provide integral insight into miRNA function throughout pig muscle development stages. Our findings will promote further development of the pig as a model organism for human age-related muscle disease research.

Guanidinoacetic acid (GAA), an amino acid derivative that is endogenous to animal tissues including muscle and nerve, has been reported to enhance muscular performance. MicroRNA (miRNA) is a post-transcriptional regulator that plays a key role in nutrient-mediated myogenesis. However, the effects of GAA on myogenic differentiation and skeletal muscle growth, and the potential regulatory mechanisms of miRNA in these processes have not been elucidated. In this study, we investigated the effects of GAA on proliferation, differentiation, and growth in C2C12 cells and mice. The results showed that GAA markedly inhibited the proliferation of myoblasts, along with the down-regulation of cyclin D1 (CCND1) and cyclin dependent kinase 4 (CDK4) mRNA expression, and the upregulation of cyclin dependent kinase inhibitor 1A (P21) mRNA expression. We also demonstrated that GAA treatment stimulated myogenic differentiation 1 (MyoD) and myogenin (MyoG) mRNA expression, resulting in an increase in the myotube fusion rate. Meanwhile, GAA supplementation promoted myotube growth through increase in total myosin heavy chain (MyHC) protein level, myotubes thickness and gastrocnemius muscle cross-sectional area. Furthermore, small RNA sequencing revealed that a total of eight miRNAs, including miR-133a-3p and miR-1a-3p cluster, showed differential expression after GAA supplementation. To further study the function of miR-133a-3p and miR-1a-3p in GAA-induced skeletal muscle growth, we transfected miR-133a-3p and miR-1a-3p mimics into myotube, which also induced muscle growth. Through bioinformatics and a dual-luciferase reporter system, the target genes of miR-133a-3p and miR-1a-3p were determined. These two miRNAs were shown to modulate the Akt/mTOR/S6K signaling pathway by restraining target gene expression. Taken together, these findings suggest that GAA supplementation can promote myoblast differentiation and skeletal muscle growth through miR-133a-3p- and miR-1a-3p-induced activation of the AKT/mTOR/S6K signaling pathway.

Intramuscular fat (IMF) content and composition are considered crucial indicators of porcine meat quality. However, the molecular mechanism of porcine IMF development is still mostly unclear. Recently, new evidence suggested that microRNA (miRNAs) play important roles in porcine intramuscular adipogenesis. Previously, microRNA-125a-5p (miR-125a-5p) was identified as an important regulator of adipogenesis. In the present study, we found that the expression of miR-125a-5p is dynamically regulated during porcine intramuscular preadipocytes differentiation and that its expression levels in different porcine muscle tissues were negatively involved with IMF content. To investigate the potential function role of miR-125a-5p in IMF development, porcine intramuscular preadipocytes were collected and transfected with miR-125a-5p mimics, inhibitors, or a negative control (NC), respectively. The results showed that overexpression of miR-125a-5p promoted proliferation and inhibited differentiation of porcine intramuscular preadipocytes while inhibition of miR-125a-5p had the opposite effects. Furthermore, a luciferase reporter assay demonstrated that porcine kruppel like factor 3 (KLF13) is a target gene of miR-125a-5p during porcine intramuscular preadipocytes differentiation. Interestingly, porcine ELOVL fatty acid elongase 6 (ELOVL6), a regulator of fatty acid composition, was also identified as a target gene of miR-125a-5p during porcine intramuscular adipogenesis. Further studies show that miR-125a-5p overexpression reduced total saturated fatty acids (SFA) content and monounsaturated fatty acids (MUFA)/SFA ratios while having no significant impact on polyunsaturated fatty acids (PUFA)/SFA and n-6/n-3 ratios. Taken together, our results identified that miR-125a-5p may be a novel regulator of porcine intramuscular adipogenesis and the fatty acid composition of porcine IMF.

Obesity is a major driver of metabolic diseases such as nonalcoholic fatty liver disease, certain cancers, and insulin resistance. However, there are no effective drugs to treat obesity. Betaine is a nontoxic, chemically stable and naturally occurring molecule. This study shows that dietary betaine supplementation significantly inhibits the white fat production in a high-fat diet (HFD)-induced obese mice. This might be due to betaine preventing the formation of new white fat (WAT), and guiding the original WAT to burn through stimulated mitochondrial biogenesis and promoting browning of WAT. Furthermore, dietary betaine supplementation decreases intramyocellular lipid accumulation in HFD-induced obese mice. Further analysis shows that betaine supplementation reduced intramyocellular lipid accumulation might be associated with increasing polyunsaturated fatty acids (PUFA), fatty acid oxidation, and the inhibition of fatty acid synthesis in muscle. Notably, by performing insulin-tolerance tests (ITTs) and glucose-tolerance tests (GTTs), dietary betaine supplementation could be observed for improvement of obesity and non-obesity induced insulin resistance. Together, these findings could suggest that inhibiting WAT production, intramyocellular lipid accumulation and inflammation, betaine supplementation limits HFD-induced obesity and improves insulin resistance.

MicroRNAs (miRNAs) are non-coding small RNA ∼22 nucleotides in length that can regulate the expression of a wide range of coding genes at the post-transcriptional level. Visceral adipose tissues (VATs) and subcutaneous adipose tissues (SATs), the two main fat compartments in mammals, are anatomically, physiologically, metabolically, and clinically distinct. Various studies of adipose tissues have focused mainly on DNA methylation, and mRNA and protein expression, nonetheless little research sheds directly light on the miRNA transcriptome differences between these two distinct adipose tissue types. Here, we present a comprehensive investigation of miRNA transcriptomes across six variant porcine adipose tissues by small RNA-sequencing. We identified 219 known porcine miRNAs, 97 novel miRNA*s, and 124 miRNAs that are conserved to other mammals. A set of universally abundant miRNAs (i.e., miR-148a-3p, miR-143-3p, miR-27b-3p, miR-let-7a-1-5p, and miR-let-7f-5p) across the distinct adipose tissues was found. This set of miRNAs may play important housekeeping roles that are involved in adipogenesis. Clustering analysis indicated significant variations in miRNA expression between the VATs and SATs, and highlighted the role of the greater omentum in responding to potential metabolic risk because of the observed enrichment in this tissue of the immune- and inflammation-related miRNAs, such as the members of miR-17-92 cluster and miR-181 family. Differential expression of the miRNAs between the VATs and SATs, and miRNA target prediction analysis revealed that the VATs-specific enriched miRNAs were associated mainly with immune and inflammation responses. In summary, the differences of miRNA expression between the VATs and SATs revealed some of their intrinsic differences and indicated that the VATs might be closely associated with increased risk of metabolic disorders.

Previous studies have indicated two main domestic pig dispersal routes in East Asia: one is from the Mekong region, through the upstream region of the Yangtze River (URYZ) to the middle and upstream regions of the Yellow River, the other is from the middle and downstream regions of the Yangtze River to the downstream region of the Yellow River, and then to northeast China. The URYZ was regarded as a passageway of the former dispersal route; however, this assumption remains to be further investigated. We therefore analyzed the hypervariable segements of mitochondrial DNA from 513 individual pigs mainly from Sichuan and the Tibet highlands and 1,394 publicly available sequences from domestic pigs and wild boars across Asia. From the phylogenetic tree, most of the samples fell into a mixed group that was difficult to distinguish by breed or geography. The total network analysis showed that the URYZ pigs possessed a dominant position in haplogroup A and domestic pigs shared the same core haplotype with the local wild boars, suggesting that pigs in group A were most likely derived from the URYZ pool. In addition, a region-wise network analysis determined that URYZ contains 42 haplotypes of which 22 are unique indicating the high diversity in this region. In conclusion, our findings confirmed that pigs from the URYZ were domesticated in situ.