Search

Obesity and diabetes mellitus have become major health problems worldwide. In recent years, genistein has been found to be capable of inhibiting obesity and alleviating insulin resistance. However, the molecular mechanism of genistein against obesity is still not fully understood. In this study, we used 3T3-L1 preadipocytes and obese mice as models to explore the molecular mechanism of genistein against obesity. We found that genistein can inhibit obesity and downregulate the expression of miR-222 in mouse adipose tissue. In 3T3-L1 preadipocytes, treatment with miR-222 inhibitor or genistein reduced the expression of miR-222 and promoted lipid decomposition, while miR-222 treatment increased the expression of miR-222 and inhibited lipolysis. Moreover, the dual-luciferase reporter assay system confirmed that BTG2 and adipor1 are the target genes of miR-222. Experiments conducted in vitro and in vivo suggest that genistein may regulate lipid metabolism in the adipose tissue of obese mice by regulating the expression of miR-222 and its target genes, BTG2 and adipor1. Our findings provide a new epigenetic mechanism underpinning the ability of genistein to resist obesity. These results may provide a reference point for the dietary treatment of obesity and type 2 diabetes mellitus.

Swift and complete spindle disassembly is essential for cell survival, yet how it happens is largely unknown. Here we used real-time live-cell microscopy and biochemical assays to show that a cysteine-rich protein CRIPT dictates the spindle disassembly in a redox-dependent manner in human cells. This previously reported cytoplasmic protein was found to have a confined nuclear localization during interphase but was distributed to spindles and underwent redox modifications to form disulfides within CXXC pairs during mitosis. Then, it interacts with and transfers redox response to tubulin subunits to induce microtubule depolymerization. The mutants with any of cysteine substitution completely block the spindle disassembly generating two cell populations with long-lasting metaphase spindles or spindle remnants. The live cell recordings of a disease-relevant mutant (CRIPTC3Y) revealed that microtubule depolymerization at spindle ends during anaphase and the entire spindle dissolution during telophase may share a common CRIPT-bearing redox-controlled mechanism.

Increasing evidence indicates that muscular dysfunction or alterations in skeletal muscle fiber-type composition not only are involved in muscle metabolism and function but also can limit functional capacity. Therefore, understanding the mechanisms regulating key events during skeletal myogenesis is necessary. Betaine is a naturally occurring component of commonly eaten foods. Here, we showed that 10 mM betaine supplementation in vitro significantly repressed myoblast proliferation and enhanced myoblast differentiation. This effect can be mediated by regulation of miR-29b-3p. Further analysis showed that betaine supplementation in vitro regulated skeletal muscle fiber-type composition through the induction of NFATc1 and the negative regulation of MyoD expression. Furthermore, mice fed with 10 mM betaine in water for 133 days showed no impairment in overall health. Consistently, betaine supplementation increased muscle mass, promoted muscle formation, and modulated the ratio of fiber types in skeletal muscle in vivo. These findings shed light on the diverse biological functions of betaine and indicate that betaine supplementation may lead to new therapies for diseases such as muscular dystrophy or other diseases related to muscle dysfunction. KEY MESSAGES: Betaine supplementation inhibits proliferation and promotes differentiation of C2C12 myoblasts. Betaine supplementation regulates fast to slow muscle fiber-type conversion and is associated with NFATc1/MyoD. Betaine supplementation enhances skeletal myogenesis in vivo. Betaine supplementation does not impair health of mice.

Background The Liangshan breed of pig is well known for its excellent meat quality attributes, but poor carcass quality and growth attributes have limited its economic value. By crossing the Liangshan with western breeds that have strong growth attributes may improve the general performance of the breed. Objectives The objective of this research was to study the effects of slaughter age and cross-breeding on reproduction, growth, carcass, meat quality, and to estimate the heterosis values of the various traits. Materials Samples were taken from 636 pigs of three different genotypes [Liangshan pure breed (LS), Duroc-Liangshan cross (DL), and Berkshire-Liangshan cross (BL)]. Results DL and BL F1 generations showed improved reproductive and growth traits than purebred Liangshan pigs. DL and BL had better carcass quality than LS, with a significant reduction in backfat thickness and sebum rate. Crossbreeds had higher L24 h and shear force, and lower marbling scores than LS. Drip loss in DL was markedly lower than in LS and BL, but BL had the lowest cooking loss, higher L45 min and pH45 min. By estimating heterosis, reproduction traits were significantly improved in both DL and BL. Growth performance showed no significant difference among the breeds. Backfat thickness and loin muscle area was significantly improved in the DL and BL breeds, although there was no difference in lean percentage and sebum rate. BL showed improved heterosis, with excellent intramuscular fat, higher pH45 min, lower L45 min and cooking loss. Conclusions The results of this study suggest that crossing Berkshire with Liangshan pigs could be an alternative genotype in the commercial pork production system.

Myogenic differentiation, which occurs in the process of muscle development, is a highly ordered process. Increasing evidence indicates that microRNAs (miRNAs) are important regulators in myogenic processes. In this study, we found that miR-204-5p expression gradually decreased when myoblasts were induced to differentiate. Our results suggested that miR-204-5p blunted myoblast differentiation, which was accompanied with a decreased proportion of myosin heavy chain (MyHC)-positive cells in myoblasts with augmented expression of miR-204-5p. Furthermore, overexpression of miR-204-5p significantly decreased the MyHC composition of slow-twitch fibers in myoblasts. Luciferase activity assays confirmed that miR-204-5p directly targeted the 3'-untranslated region (3'-UTR) of myocyte enhancer factor 2C (MEF2C) and estrogen-related receptor gamma (ERRγ). Small interfering RNA (siRNA) technology successfully inhibited the expression of MEF2C and ERRγ. Interference with MEF2C or ERRγ inhibited myoblast differentiation and the formation of slow-twitch fibers. Meanwhile, co-transfection of either si-MEF2C or si-ERRγ with miR-204-5p mimics resulted in a more severe attenuation of myogenic differentiation. In summary, this study demonstrates that miR-204-5p inhibits myoblast differentiation by targeting MEF2C and ERRγ. Our findings suggest that miR-204-5p is a potential regulator that could influence myogenesis.

In this study, data genotyping by sequence (GBS) was used to perform single step GWAS (ssGWAS) to identify SNPs associated with the litter traits in domestic pigs and search for candidate genes in the region of significant SNPs. After quality control, 167,355 high-quality SNPs from 532 pigs were obtained. Phenotypic traits on 2112 gilt litters from 532 pigs were recorded including total number born (TNB), number born alive (NBA), and litter weight born alive (LWB). A single-step genomic BLUP approach (ssGBLUP) was used to implement the genome-wide association analysis at a 5% genome-wide significance level. A total of 8, 23 and 20 significant SNPs were associated with TNB, NBA, and LWB, respectively, and these significant SNPs accounted for 62.78%, 79.75%, and 58.79% of genetic variance. Furthermore, 1 (SSC14: 16314857), 4 (SSC1: 81986236, SSC1: 66599775, SSC1: 161999013, and SSC1: 267883107), and 5 (SSC9: 29030061, SSC2: 32368561, SSC5: 110375350, SSC13: 45619882 and SSC13: 45647829) significant SNPs for TNB, NBA, and LWB were inferred to be novel loci. At SSC1, the AIM1 and FOXO3 genes were found to be associated with NBA; these genes increase ovarian reproductive capacity and follicle number and decrease gonadotropin levels. The genes SLC36A4 and INTU are involved in cell growth, cytogenesis and development were found to be associated with LWB. These significant SNPs can be used as an indication for regions in the Sus scrofa genome for variability in litter traits, but further studies are expected to confirm causative mutations.

The development of skeletal muscle is a complex process involving the proliferation, differentiation, apoptosis, and changing of muscle fiber types in myoblasts. Many reports have described the involvement of microRNAs in the myogenesis of myoblasts. In this study, we found that the expression of miR-152 was gradually down-regulated during myoblast proliferation, but gradually up-regulated during the differentiation of myoblasts. Transfection with miR-152 mimics restrained cell proliferation and decreased the expression levels of cyclin E, CDK4, and cyclin D1, but promoted myotube formation and significantly increased the mRNA expression levels of MyHC, MyoD, MRF4, and MyoG in C2C12 myoblasts. However, treatment with miR-152 inhibitors promoted cell proliferation and restrained differentiation. Moreover, over-expression of miR-152 significantly decreased E2F3 production in C2C12 myoblasts. A luciferase assay confirmed that miR-152 could bind to the 3' UTR of E2F3. In conclusion, this study showed that miR-152 inhibited proliferation and promoted myoblast differentiation by targeting E2F3.

Aging significantly affects the cardiac muscle (CM) and skeletal muscles (SM). Since the aging process of CM and SM may be different, high throughput RNA sequencing was performed using CM and SM in different age conditions to evaluate the expression profiles of messenger RNA (mRNA), long non-coding RNA (lncRNA), micro RNA (miRNA), and circular (circRNA). Several mRNAs, lncRNAs, and miRNAs were highly expressed and consistently appeared in both ages in one of the two muscle tissues. Gene ontology (GO) annotation described that these genes were required for maintaining normal biological functions of CM and SM tissues. Furthermore, 26 mRNAs, 4 lncRNAs, 22 miRNAs, and 26 circRNAs were differentially expressed during cardiac muscle aging. Moreover, 81 mRNAs, 5 lncRNAs, 79 miRNAs, and 62 circRNAs were differentially expressed during aging of skeletal muscle. When comparing the expression profiles of CM and SM during aging, the senescence process in CM and SM was found to be fundamentally different. In addition, we assessed multi-group cooperative control relationships and constructed circRNA-miRNA-mRNA co-expression networks in muscular aging. In conclusion, our findings will contribute to the understanding of muscular aging and provide a foundation for future studies on the molecular mechanisms underlying muscular aging.

Background Species living at high altitude are subject to strong selective pressures due to inhospitable environments (e.g., hypoxia, low temperature, high solar radiation, and lack of biological production), making these species valuable models for comparative analyses of local adaptation. Studies that have examined high-altitude adaptation have identified a vast array of rapidly evolving genes that characterize the dramatic phenotypic changes in high-altitude animals. However, how high-altitude environment shapes gene expression programs remains largely unknown. Findings We generated a total of 910 Gb of high-quality RNA-seq data for 180 samples derived from 6 tissues of 5 agriculturally important high-altitude vertebrates (Tibetan chicken, Tibetan pig, Tibetan sheep, Tibetan goat, and yak) and their cross-fertile relatives living in geographically neighboring low-altitude regions. Of these, ∼75% reads could be aligned to their respective reference genomes, and on average ∼60% of annotated protein coding genes in each organism showed FPKM expression values greater than 0.5. We observed a general concordance in topological relationships between the nucleotide alignments and gene expression–based trees. Tissue and species accounted for markedly more variance than altitude based on either the expression or the alternative splicing patterns. Cross-species clustering analyses showed a tissue-dominated pattern of gene expression and a species-dominated pattern for alternative splicing. We also identified numerous differentially expressed genes that could potentially be involved in phenotypic divergence shaped by high-altitude adaptation. Conclusions These data serve as a valuable resource for examining the convergence and divergence of gene expression changes between species as they adapt or acclimatize to high-altitude environments.

The whole-genome sequencing (WGS) data can potentially discover all genetic variants. Studies have shown the power of WGS for genome-wide association study (GWAS) lies in the ability to identify quantitative trait loci and nucleotides (QTNs). However, the resequencing of thousands of target individuals is expensive. Genotype imputation is a powerful approach for WGS and to identify causal mutations. This study aimed to evaluate the imputation accuracy from genotyping-by-sequencing (GBS) to WGS in two pig breeds using a resequencing reference population and to detect single-nucleotide polymorphisms (SNPs) and candidate genes for farrowing interval (FI) of different parities using the data before and after imputation for GWAS. Six hundred target pigs, 300 Landrace and 300 Large White pigs, were genotyped by GBS, and 60 reference pigs, 20 Landrace and 40 Large White pigs, were sequenced by whole-genome resequencing. Imputation for pigs was conducted using Beagle software. The average imputation accuracy (allelic R2) from GBS to WGS was 0.42 for Landrace pigs and 0.45 for Large White pigs. For Landrace pigs (Large White pigs), 4,514,934 (5,533,290) SNPs had an accuracy >0.3, resulting an average accuracy of 0.73 (0.72), and 2,093,778 (2,468,645) SNPs had an accuracy >0.8, resulting an average accuracy of 0.94 (0.93). Association studies with data before and after imputation were performed for FI of different parities in two populations. Before imputation, 18 and 128 significant SNPs were detected for FI in Landrace and Large White pigs, respectively. After imputation, 125 and 27 significant SNPs were identified for dataset with an accuracy >0.3 and 0.8 in Large White pigs, and 113 and 18 SNPs were found among imputed sequence variants. Among these significant SNPs, six top SNPs were detected in both GBS data and imputed WGS data, namely, SSC2: 136127645, SSC5: 103426443, SSC6: 27811226, SSC10: 3609429, SSC14: 15199253, and SSC15: 150297519. Overall, many candidate genes could be involved in FI of different parities in pigs. Although imputation from GBS to WGS data resulted in a low imputation accuracy, association analyses with imputed WGS data were optimized to detect QTNs for complex trait. The obtained results provide new insight into genotype imputation, genetic architecture, and candidate genes for FI of different parities in Landrace and Large White pigs.

Background The number of animals born dead, which includes the number of mummified (NM) and stillborn (NS) animals, is the most important trait to directly quantify the reproductive loss in domestic pigs. In this study, 282 Landrace sows and 250 Large White sows were genotyped by sequencing (GBS). A total of 816 and 1068 litter records for NM and NS were collected from them. A genome-wide association study (GWAS) was conducted to reveal the genetic difference between NM and NS. Results A total of 248 and 10 genome-wide significant SNPs were detected for NM and NS across numerous parities in Landrace pigs. The corresponding numbers for Large White pigs were 175 and 6, respectively. All of the detected SNPs were parity specific for both NM and NS in two breeds. Based on significant SNPs, in total 242 (146 for Landrace pig, 96 for Large White pig) and 10 significant chromosome regions (8 for Landrace pigs, 2 for Large White pigs) were found for NM and NS, respectively. Among them, 237 (142 for Landrace pig, 95 for Large White pig) and 8 significant chromosome regions (6 for Landrace pigs, 2 for Large White pigs) for NM and NS were not reported in previous studies. A list of candidate genes at the identified loci was proposed, including HMGB1, SOX5, KCNJ8, ABCC9 and YY1 for NM, ASTN1 for NS. Conclusion This is the first time when GBS data was used to identify genetic regions affecting NM and NS in Landrace and Large White pigs. Many identified informative SNPs and candidate genes advance our understanding of the genetic architecture of NM and NS in pigs. However, further studies are needed to validate using larger populations with more breeds. Electronic supplementary material The online version of this article (10.1186/s12863-018-0692-x) contains supplementary material, which is available to authorized users.

Total number born (TNB), number born alive (NBA), and litter weight born alive (LWB) are critically important traits in pig production. The sow’s parity is one of the major factors influencing litter traits. Because of monogenic or polygenic contributions and the presence of temporal gene effects in different sows’ parities, it is difficult to clarify the biological and genetic background. To systematically explore the genetic mechanism of litter traits, we conducted 18 GWASs using single-step GWAS (ssGWAS) based on two breeds (908 Landrace and 1,130 Large White sow litter records) for each litter trait in different parities. A total of 300 Landrace and 300 Large White sows were genotyped by sequencing (GBS). ssGWAS was performed separately for each breed and each parity due to population stratification and temporal gene effect. In summary, we identified 80 (15 for Landrace and 65 for Large White), 227 (52 for Landrace, 175 for Large White), and 187 (34 for Landrace, 153 for Large White) single nucleotide polymorphisms (SNPs) affecting TNB, NBA, and LWB, respectively. Of them, we suggest that a total of 22 loci (SSC1: 125098202, SSC1: 117560058, SSC14: 147794697, SSC8: 84823302, SSC9: 143554876, and SSC9: 138766097 for Landrace; SSC1: 4023577, SSC1: 3859573, SSC1: 4891063, SSC16: 5197665, SSC10: 32050819, SSC13: 13552924, SSC13: 92819, SSC17: 3579607, SSC13: 196698221, SSC7: 30918403, SSC16: 46221484, SSC16: 46169204, SSC2: 41988642, SSC2: 44475457, SSC2: 42521875, and SSC7: 58411951 for Large White) are shared by TNB, NBA, and LWB. These results indicate the existence of gene temporal effect in each parity. Furthermore, our findings suggest four interesting candidate genes (FBXL7, ALDH1A2, LEPR, and DDX1) associated with litter traits in different parities that have a major effect on embryonic development progression. In conclusion, 22 crucial SNPs and four interesting candidate genes were identified for three litter traits across six parities. These findings advance our understanding of the genetic architecture of litter traits and confirm the presence of temporal gene effects in different parities. Importantly, functional validation studies for findings of particular interest are recommended in litter traits.

Uncovering genetic variation through resequencing is limited by the fact that only sequences with similarity to the reference genome are examined. Reference genomes are often incomplete and cannot represent the full range of genetic diversity as a result of geographical divergence and independent demographic events. To more comprehensively characterize genetic variation of pigs (Sus scrofa), we generated de novo assemblies of nine geographically and phenotypically representative pigs from Eurasia. By comparing them to the reference pig assembly, we uncovered a substantial number of novel SNPs and structural variants, as well as 137.02-Mb sequences harboring 1737 protein-coding genes that were absent in the reference assembly, revealing variants left by selection. Our results illustrate the power of whole-genome de novo sequencing relative to resequencing and provide valuable genetic resources that enable effective use of pigs in both agricultural production and biomedical research.

Domestication and subsequent selective pressures have produced a large variety of pig coat colors in different regions and breeds. The melanocortin 1 receptor (MC1R) gene plays a crucial role in determining coat color of mammals. Here, we investigated genetic diversity and selection at the coding region of the porcine melanocortin receptor 1 (MC1R) in Tibetan pigs and Landrace pigs. By contrast, genetic variability was much lower in Landrace pigs than in Tibetan pigs. Meanwhile, haplotype analysis showed that Tibetan pigs possessed shared haplotypes, suggesting a possibility of recent introgression event by way of crossbreeding with neighboring domestic pigs or shared ancestral polymorphism. Additionally, we detected positive selection at the MC1R in both Tibetan pigs and Landrace pigs through the dN/dS analysis. These findings suggested that novel phenotypic change (dark coat color) caused by novel mutations may help Tibetan pigs against intensive solar ultraviolet (UV) radiation and camouflage in wild environment, whereas white coat color in Landrace were intentionally selected by human after domestication. Furthermore, both the phylogenetic analysis and the network analysis provided clues that MC1R in Asian and European wild boars may have initially experienced different selective pressures, and MC1R alleles diversified in modern domesticated pigs.

Growth and fat deposition are important economic traits due to the influence on production in pigs. In this study, a dataset of 1200 pigs with 345,570 SNPs genotyped by sequencing (GBS) was used to conduct a GWAS with single-marker regression method to identify SNPs associated with body weight and backfat thickness (BFT) and to search for candidate genes in Landrace and Yorkshire pigs . A total of 27 and 13 significant SNPs were associated with body weight and BFT, respectively. In the region of 149.85–149.89 Mb on SSC6, the SNP (SSC6: 149876737) for body weight and the SNP (SSC6: 149876507) for BFT were in the same locus region (a gap of 230 bp). Two SNPs were located in the DOCK7 gene, which is a protein-coding gene that plays an important role in pigmentation . Two SNPs located on SSC8: 54567459 and SSC11: 33043081 were found to overlap weight and BFT; however, no candidate gene was found in these regions. In addition, based on other significant SNPs, two positional candidate genes, NSRP1 and CADPS, were proposed to influence weight. In conclusion, this is the first study report using GBS data to identify the significant SNPs for weight and BFT. A total of four particularly interesting SNPs and one potential candidate genes (DOCK7) were found for these traits in domestic pigs. This study improves our knowledge to better understand the complex genetic architecture of weight and BFT, but further validation studies of these candidate loci and genes are recommended in pigs

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.

Excessive accumulation of adipose tissue is a main cause of obesity or overweight, which is significantly involved in increasing the risk of diseases. Recently, numerous studies have proved that microRNAs (miRNAs) play important roles in adipogenesis by negatively regulating gene expression at posttranscriptional levels. In this study, we showed that miR-125a-5p was expressed at lower levels in the adipose tissues of high-fat diet (HFD)-fed mice than the normal chow (NCW)-fed mice. MiR-125a-5p expression were strongly up-regulated by nearly five-fold, when 3T3-L1 preadipocyte were induced and differentiated into mature adipocytes. Functional analysis indicated that overexpression of miR-125a-5p promoted 3T3-L1 preadipocyte proliferation and inhibited its differentiation. By contrast, inhibition of miR-125a-5p repressed 3T3-L1 preadipocyte proliferation and accelerated its differentiation. Furthermore, a dual-luciferase reporter assay demonstrated that signal transducer and activator of transcription 3 (STAT3) is a direct target gene of miR-125a-5p during 3T3-L1 preadipocyte differentiation. Further analysis confirmed that the process of miR-125a-5p inhibiting 3T3-L1 preadipocyte differentiation might be associated with the regulation of fatty acid metabolism related genes. Taken together, our results indicated that miR-125a-5p might promote 3T3-L1 preadipocyte proliferation, whereas inhibiting 3T3-L1 preadipocyte differentiation by negatively regulating STAT3.

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.

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.

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.

Glucose metabolism is a basic biological process that shows substantial variation within and between species. Using pig as a model organism, we investigated differences in glucose metabolic genes in seven tissues from domesticated pigs (Rongchang pig and Tibetan pig, meanwhile, the Tibetan pig just as a special case of the domesticated pig under plateau condition) and wild boar. We found large differences in the expression of genes involved in multiple aspects of glucose metabolism, including genes associated with glucose transport, gluconeogenesis, and glycolysis. In addition, we identified microRNAs (miRNAs) that may be involved in the divergence of glucose metabolism in pig. A combined analysis of mRNA and miRNA expression indicated that some miRNA:mRNA pairs showed ab facto function in it. Our results provide a valuable resource for further determination of miRNA regulatory roles in pig glucose metabolism and reveal the divergence of glucose metabolism in pigs under domestication.

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.

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. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3719-1) contains supplementary material, which is available to authorized users.

Recent evidence suggests that hypoxia caused by acute myocardial infarction can induce cardiomyocyte apoptosis. Exosomes are signalling mediators that contribute to intercellular communication by transporting cytosolic components including miRNAs, mRNAs, and proteins. However, the systemic regulation and function of exosomal miRNAs in hypoxic cardiomyocytes are currently not well understood. Here, we used small RNA sequencing to investigate the effects of hypoxia stress on miRNAome of rat cardiomyoblast cells (H9c2) and corresponding exosomes. We identified 92 and 62 miRNAs in cells and exosomes, respectively, that were differentially expressed between hypoxia and normoxia. Hypoxia strongly modulated expression of hypoxia-associated miRNAs in H9c2 cells, and altered the miRNAome of H9c2 cells-derived exosomes. Functional enrichment analysis revealed extensive roles of differentially expressed exosomal miRNAs in the HIF-1 signalling pathway and in apoptosis-related pathways including the TNF, MAPK, and mTOR pathways. Furthermore, gain- and loss-of-function analysis demonstrated potential anti-apoptotic effects of the hypoxia-induced exosomal miRNAs, including miR-21-5p, miR-378-3p, miR-152-3p, and let-7i-5p; luciferase reporter assay confirmed that Atg12 and Faslg are targets of miR-152-3p and let-7i-5p, respectively. To summarize, this study revealed that hypoxia-induced exosomes derived from H9c2 cells loaded cardioprotective miRNAs, which mitigate hypoxia-induced H9c2 cells apoptosis.

The conventional BLUP uses the phenotype and pedigree information to predict the estimated breeding values (EBV) of individuals in genetic evaluation. However, to obtain the phenotypes of interest in swine breeding program requires relatively a long time until the completion of performance testing, which produces an extra breeding cost for the culled pigs in swine industry. An alternative solution is to pre-select a predefined proportion of young replacement piglets for future performance testing through genomic selection, which could reduce the number of testing animals entering into performance testing program and hence reduce the breeding costs associated with the tests. In this study, four strategies of genomic selection applications in swine breeding program were compared through simulation to investigate the potential benefits in the different strategies in swine breeding program. For comparison purpose, the conventional BLUP selection was simulated as strategy 1 to form the benchmark basis for comparisons. Strategy 2 was an extreme case for applying genomic selection, where newborn piglets were selected directly based on their genomic enhanced breeding value (GEBV) only without further performance testing. Strategy 3 was to pre-select piglets, based on their GEBV, for entry to performance testing in early stage, and then the breeding stock were selected ultimately based on their EBVs predicted by BLUP method when the phenotypic records were available. Similar to strategy 3, strategies 4 and 5 also used the GEBVs to pre-select replacement piglets in an early stage; however, the breeding stocks in strategies 4 and 5 were selected based on the breeding values obtained using the bi-variable model and the conventional index method to combine GEBV and EBV information, respectively, when individuals had both GEBV and phenotypes available after the performance testing. Comparing these strategies, strategy 4 resulted in the highest accuracy in first three generations and achieved the best cumulative selection response in the last generation, followed by strategies 1, 3, 5, and 2. The proportion of pre-selection of boars and sows in the early stage affected the efficiency of genomic selection substantially.

Gut microbiota play important roles in host health and diseases. Several factors, in particular antibiotics, affect the gut microbiota of pigs. Cross-fostering has been applied as a regular practice to equalize litter size, reduce pre-weaning mortality and increase body weight. However, the effect of cross-fostering on cecal microbiota is unclear. In this study, we fostered three Jinhua pigs to a Yorkshire sow. The fostered Jinhua piglets grew significantly faster than their biological siblings. To explore whether the cecal microbiota of piglets will alter during fostering, we characterized the cecal microbiota of all piglets by examining the V3 hypervariable region of 16S rDNA. We observed altered cecal microbiota in these piglets using the Illumina HiSeq 2000 platform, and this was accompanied with an increase in growth rate after fostering. The relative abundance of Bacteroidetes in the fostered Jinhua piglets was decreased compared with their biological siblings, although still higher in comparison with their new littermates. Beta-diversity analysis also showed that the cecal microbiota of the adopted Jinhua piglets differed from their biological siblings with a shift toward their step-siblings. Our data show that cecal microbiota of piglets were altered after cross-fostering while the growth rate increased.

Adipose triglyceride lipase (ATGL), catalyzing the initial step of hydrolysis of triacylglycerol (TAG) in adipocytes, has been known to be inhibited by G0/G1 switch protein 2 (G0S2). In this study, we determined tissue expression pattern and polymorphism of G0S2 gene in porcine. The results showed that the G0S2 transcript levels were very high in the liver and, to a lesser degree, in adipose tissues of greater omentum and suet fat; and low G0S2 transcript levels were observed in other tissues. A comparative study on the transcript levels between ATGL and G0S2 genes showed that ATGL transcript levels were high in all six adipose tissues, but negligible in the liver. Higher transcript levels were obtained for sows in adipose tissues of the inner layer of subcutaneous fat and suet fat, but higher expression values were found for boars in the liver, spleen, and stomach. 19 single nucleotide polymorphisms (SNPs), including 4 nonsynonymous SNPs (g.-307A>T, g.-394C>G, g.-565G>A, and g.-566T>C), were found in porcine G0S2 genomic DNA. Association analyses showed that the g.-565G>A and g.-742T>A SNPs were associated with back fat thickness (BFT). In conclusion, G0S2 mRNAs are abundantly expressed in porcine liver and adipose tissues of greater omentum and suet fat, and sex affects porcine G0S2 tissue transcript levels; meanwhile, the genetic diversity of porcine G0S2 gene is abundant and 2 SNPs are a genetic factor affecting BFT.

Here, we analyzed the distribution of H-FABP/(HinfI, MspI, and HaeIII) and ACSL4/RsaI polymorphisms, and the associations of these 4 polymorphic loci with intramuscular fat (IMF) content and backfat thickness (BFT) in Yanan, Jinhua, Duroc, Landrace, Yorkshire, and Duroc x (Landrace x Yorkshire) (DLY) pigs. H-FABP/HinfI polymorphisms were present in all the 6 populations. At the ACSL4/RsaI locus, sows had 3 genotypes, whereas boars only had haplotype A or G, in Duroc, Landrace, Yorkshire, and DLY pigs. H-FABP/(MspI and HaeIII) and ACSL4/RsaI polymorphisms were absent in Yanan and Jinhua pigs. Linkage disequilibrium analysis indicated that the 3 loci (HinfI, MspI, and HaeIII) were separated. Association analysis showed that the H-FABP/HinfI locus significantly affected IMF content in DLY (P0.05) and Yanan (P0.001) pigs. The highest IMF content was recorded in the adH haplotype of the 3 H-FABP polymorphic loci (2.59%, P0.05) in DLY pigs. At the ACSL4/RsaI locus, higher IMF content was recorded for sows with a GG genotype or boars with a G haplotype compared to those with an AA genotype (2.53 vs 2.10%, P0.05) or A haplotype (2.48 vs 1.73%, P0.01) in DLY pigs. Significant differences were not obtained among these 4 polymorphic loci and BFT (P0.05). The results indicate that H-FABP and ACSL4 genes might serve as markers to improve IMF content (but not BFT) in the pig breeding system.

The balance between biosynthesis and oxidation of fatty acids determines adipose deposition in mammals. Obese and lean pigs show obvious differences in total adipose mass and therefore offer an attractive model for comparative studies. We found that obese Rongchang pigs, when compared with lean Landrace pigs, exhibited significantly higher mRNA levels for five genes encoding acyl-CoA dehydrogenases involved in mitochondrial fatty-acid β-oxidation in eight different adipose tissues. These changes in gene expression were positively correlated with adipocyte volume in the eight adipose tissues. Based on these results, we hypothesize that acyl-CoA dehydrogenase genes participate in the regulation of fat mass in pigs.

We report the sequencing at 131× coverage, de novo assembly and analyses of the genome of a female Tibetan wild boar. We also resequenced the whole genomes of 30 Tibetan wild boars from six major distributed locations and 18 geographically related pigs in China. We characterized genetic diversity, population structure and patterns of evolution. We searched for genomic regions under selection, which includes genes that are involved in hypoxia, olfaction, energy metabolism and drug response. Comparing the genome of Tibetan wild boar with those of neighboring Chinese domestic pigs further showed the impact of thousands of years of artificial selection and different signatures of selection in wild boar and domestic pig. We also report genetic adaptations in Tibetan wild boar that are associated with high altitudes and characterize the genetic basis of increased salivation in domestic pig.

The accumulation of inbreeding and the loss of genetic diversity is a potential problem in the modern swine breeds in China. Therefore, the purpose of this study was to analyze the pedigrees of Chinese Duroc (CD), Landrace (CL) and Yorkshire (CY) swine to estimate the past and current rates of inbreeding, and to identify the main causes of genetic diversity loss. Pedigree files from CD, CL and CY containing, 4529, 16,776 and 22,600 records, respectively, were analyzed. Pedigree completeness indexes of the three breeds, accounting for one generation back, were 83.72, 93.93 and 93.59%, respectively. The estimated average annual inbreeding rates for CD, CL and CY in recent three years were 0.21, 0.19 and 0.13%, respectively. The estimated average percentage of genetic diversity loss within each breed in recent three years was about 8.92, 2.19, and 3.36%, respectively. The average relative proportion of genetic diversity loss due to unequal contributions of founders in CD, CL and CY was 69.09, 57.95 and 60.57%, and due to random genetic drift was 30.91, 42.05 and 39.43%, respectively. The estimated current effective population size for CD, CL and CY was 76, 117 and 202, respectively. Therefore, CD has been found to have lost considerable genetic diversity, demanding priority for optimizing the selection and mating to control future coancestry and inbreeding. Unequal contribution of founders was a major cause of genetic diversity loss in Chinese swine breeds and random genetic drift also showed substantial impact on the loss of diversity.

Functional differences in the different types of adipose tissue and the impact of their dysfunction on metabolism are associated with the regional distribution of adipose depots. Here we show a genome-wide comparison between the transcriptomes of one source of subcutaneous and two sources of visceral adipose tissue in the pig using an RNA-seq approach. We obtained ~32.3 million unique mapped reads which covered ~80.2% of the current annotated transcripts across these three sources of adipose tissue. We identified various genes differentially expressed between subcutaneous and visceral adipose tissue, which are potentially associated with the inflammatory features of visceral adipose tissue. These results are of benefit for understanding the phenotypic, metabolic and functional differences between different types of adipose tissue that are deposited in different body sites.

Accumulating evidence has indicated that microRNAs (miRNAs) and endoplasmic reticulum (ER) stress play critical roles in myoblast differentiation. However, the regulation roles of miRNAs and ER stress in myogenic differentiation have not been fully revealed and need to be further studied. Here, we discovered that the expression levels of miR-181a-5p were strongly upregulated during C2C12 cell differentiation. miR-181a-5p overexpression promoted ER stress and differentiation of C2C12 cells, which was accompanied by increasing expression levels of marker genes related to ER stress-mediated apoptosis and myogenic differentiation. Opposite results were observed after inhibition of the miR-181a-5p expression. The gain- and loss-of-function experiments on C2C12 cells showed that miR-181a-5p affected the development of muscle fiber type, but had no significant influence on C2C12 cell proliferation. In the ER-stressed C2C12 cells induced by thapsigargin (Tg), the expression levels of both miR-181a-5p and marker genes related to ER stress and myogenesis were upregulated. In the ER-stressed C2C12 cells and porcine muscle fibroblast (PMF) cells pretreated with Tg, we found that miR-181a-5p targeted glucose-regulated protein, 78kDa/binding immunoglobulin protein (GRP78/BIP), and influenced cell apoptosis. In conclusion, these results indicate that miR-181a-5p and ER stress have positive synergistic effects on myogenic differentiation by increasing the expression levels of myogenic differentiation key genes and activating the ER stress-mediated apoptosis signaling pathway.

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,SCDandTUBA1B) 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.

Adipose tissue is not only a storage organ involved in fuel metabolism, but also an endocrine organ involved in the regulation of insulin sensitivity, thermogenesis, immunity, and inflammation. There are anatomical, cellular, molecular and physiological differences among adipose tissues deposited in different body sites. However, current understanding of the intrinsic differences between the sub-compartments of the subcutaneous adipose tissue remains rudimentary. Here, we analyzed the genome-wide DNA methylation differences between the porcine superficial and deep backfat tissues using methylated DNA immunoprecipitation combined with high-throughput sequencing. We show that the genes with differentially methylated regions in their promoter are mainly involved in the processes of “lipid metabolism” and “regulation of immune-related cytokines”. Compared with the deep backfat tissue, the promoters of genes related to the ‘positive regulation of cytokine production’ were significantly hypermethylated in the superficial backfat tissue, which reflects the intrinsic functional and metabolic differences between the sub-compartments of the subcutaneous adipose tissue. This study provides epigenetic evidence for functionally relevant methylation differences between different layers of porcine backfat tissues.

It is well established that the metabolic risk factors of obesity and its comorbidities are more attributed to adipose tissue distribution rather than total adipose mass. Since emerging evidence suggests that epigenetic regulation plays an important role in the aetiology of obesity, we conducted a genome-wide methylation analysis on eight different adipose depots of three pig breeds living within comparable environments but displaying distinct fat level using methylated DNA immunoprecipitation sequencing. We aimed to investigate the systematic association between anatomical location-specific DNA methylation status of different adipose depots and obesity-related phenotypes. We show here that compared to subcutaneous adipose tissues which primarily modulate metabolic indicators, visceral adipose tissues and intermuscular adipose tissue, which are the metabolic risk factors of obesity, are primarily associated with impaired inflammatory and immune responses. This study presents epigenetic evidence for functionally relevant methylation differences between different adipose depots.

It is well established that the metabolic risk factors of obesity and its comorbidities are more attributed to adipose tissue distribution rather than total adipose mass. Since emerging evidence suggests that epigenetic regulation plays an important role in the aetiology of obesity, we conducted a genome-wide methylation analysis on eight different adipose depots of three pig breeds living within comparable environments but displaying distinct fat level using methylated DNA immunoprecipitation sequencing. We aimed to investigate the systematic association between anatomical location-specific DNA methylation status of different adipose depots and obesity-related phenotypes. We show here that compared to subcutaneous adipose tissues which primarily modulate metabolic indicators, visceral adipose tissues and intermuscular adipose tissue, which are the metabolic risk factors of obesity, are primarily associated with impaired inflammatory and immune responses. This study presents epigenetic evidence for functionally relevant methylation differences between different adipose depots.

Determination of an optimal number/set of endogenous control (EC) microRNA (miRNA) genes is a critical but often an underappreciated aspect of quantitative gene expression analysis. In this study, the expression stabilities of 13 selected porcine EC miRNA genes were compared in all 47 tissue-specific normal tissues, 10 types of adipose tissue, and four types of muscle tissue using an EvaGreen quantitative PCR approach. Seven, 12, and 11 genes exhibited credible stability in the three groups, respectively. Our analysis clearly showed that three optimal EC genes are adequate for an accurate normalization, which correlated well with the theoretical optimal number (r ≥ 0.841). In terms of economical and experimental feasibility, we recommend the use of the three most stable EC miRNA genes for calculating the normalization factor, that is, ssc-miR-17, -103 and -107 for all 47 different tissues. We also suggest that two sets of EC miRNA genes are appropriate for 10 types of adipose tissue (ssc-miR-17, -107 and -24) and four types of muscle tissue (ssc-miR-17, -23a and -103), respectively. We envision that these results will serve as a valuable reference for other studies aimed at measuring tissue-specific miRNA abundance in porcine samples.

To normalize a set of quantitative real-time PCR (q-PCR) data, it is essential to determine an optimal number/set of housekeeping genes, as the abundance of housekeeping genes can vary across tissues or cells during different developmental stages, or even under certain environmental conditions. In this study, of the 20 commonly used endogenous control genes, 13, 18 and 17 genes exhibited credible stability in 56 different tissues, 10 types of adipose tissue and five types of muscle tissue, respectively. Our analysis clearly showed that three optimal housekeeping genes are adequate for an accurate normalization, which correlated well with the theoretical optimal number (r ≥ 0.94). In terms of economical and experimental feasibility, we recommend the use of the three most stable housekeeping genes for calculating the normalization factor. Based on our results, the three most stable housekeeping genes in all analysed samples (TOP2B, HSPCB and YWHAZ) are recommended for accurate normalization of q-PCR data. We also suggest that two different sets of housekeeping genes are appropriate for 10 types of adipose tissue (the HSPCB, ALDOA and GAPDH genes) and five types of muscle tissue (the TOP2B, HSPCB and YWHAZ genes), respectively. Our report will serve as a valuable reference for other studies aimed at measuring tissue-specific mRNA abundance in porcine samples.

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.

Lipoprotein lipase (LPL) is a key enzyme in lipid metabolism and is transported by glycosylphosphatidylinositol-anchored high-density lipoprotein-binding protein 1 (GPIHBP1) from the interstitial spaces to the capillary lumen. Here, we cloned a cDNA and the genomic locus of the porcine GPIHBP1 gene, and investigated its tissue expression pattern and its genetic effects on adipose traits. Porcine GPIHBP1 exhibits a four-exon/three-intron structure, including a 543 bp open reading frame that encodes 180 amino acids. The porcine GPIHBP1 protein shows 49%–65% homology and shares the major conserved structural domains of GPIHBP1 proteins in other mammals. Porcine GPIHBP1 mRNA levels were high in the adipose tissue, muscle and lung, and higher mRNA levels were observed in sows compared to boars in adipose tissues of the inner and outer layers of subcutaneous fat, abdominal fat, and suet fat. The mRNA expression pattern of porcine GPIHBP1 and LPL genes was similar in most tissues except for the lung. Thirty six single nucleotide polymorphisms (SNPs) were found in the porcine GPIHBP1 gene. Association analyses showed that the g.-255G>C and g.-626T>G SNPs are associated with intramuscular fat content, and that the g.-1557T>C and g.-1948G>A SNPs are associated with back fat thickness. In conclusion, porcine GPIHBP1 mRNA is abundantly expressed in the adipose tissue, muscle and lung, and gender affects GPIHBP1 mRNA expression levels; furthermore, four GPIHBP1 SNPs are genetic factors affecting adipose traits.

Domesticated organisms have experienced strong selective pressures directed at genes or genomic regions controlling traits of biological, agricultural or medical importance. The genome of native and domesticated pigs provide a unique opportunity for tracing the history of domestication and identifying signatures of artificial selection. Here we used whole-genome sequencing to explore the genetic relationships among the European native pig Berkshire and breeds that are distributed worldwide, and to identify genomic footprints left by selection during the domestication of Berkshire. Numerous nonsynonymous SNPs-containing genes fall into olfactory-related categories, which are part of a rapidly evolving superfamily in the mammalian genome. Phylogenetic analyses revealed a deep phylogenetic split between European and Asian pigs rather than between domestic and wild pigs. Admixture analysis exhibited higher portion of Chinese genetic material for the Berkshire pigs, which is consistent with the historical record regarding its origin. Selective sweep analyses revealed strong signatures of selection affecting genomic regions that harbor genes underlying economic traits such as disease resistance, pork yield, fertility, tameness and body length. These discoveries confirmed the history of origin of Berkshire pig by genome-wide analysis and illustrate how domestication has shaped the patterns of genetic variation.

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.

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.

We evaluated carcass and meat quality traits of two Chinese native crossbreeds Landrace x Meishan (LM) and Duroc x (Landrace x Meishan) (DLM) and two foreign crossbreeds Duroc x (Landrace x Yorkshire) (DLY) and PIC (an imported five-way crossbreed). One hundred and twenty weaned pigs (half castrated males and half females) were reared and slaughtered at a predestinated slaughter age. The general carcass and meat quality traits were measured and analyzed. The DLY and PIC crosses had significantly heavier live weights (93.39 and 96.33 kg, P < 0.01), significantly higher dressing percentages (80.65 and 79.39%, P < 0.05), significantly bigger loin areas (42.69 and 43.91 cm(2), P < 0.001), and significantly more lean carcasses (65.78 and 66.40%, P < 0.001) than LM and DLM. On the other hand, LM had a significantly lower live weight (70.29 kg, P < 0.01), significantly thicker back fat (3.54 cm, P < 0.001), significantly less lean carcasses (46.82%, P < 0.001), and significantly less ham and breech (26.53%, P < 0.05) than the other crossbreeds. Among meat quality parameters, LM had the highest intramuscular fat content (5.02%, P < 0.001) and the smallest fiber area (3126.45 μm(2), P < 0.01). However, PIC showed the lowest pH(1) (5.82, P < 0.01) and pH(2) (5.63, P < 0.01), the highest drip loss (2.89%, P < 0.01), and the lowest intramuscular fat (1.35%, P < 0.001). We concluded that LM and DLM had good meat quality traits but poorer carcass traits than DLY and PIC; DLY had good carcass and meat quality traits; PIC had good carcass traits, but it had less intramuscular fat, lower pH and higher drip loss.