Integration of multi-omics data reveals cis-regulatory variants that are associated with phenotypic differentiation of eastern from western pigs.

BACKGROUND: The genetic mechanisms that underlie phenotypic differentiation in breeding animals have important implications in evolutionary biology and agriculture. However, the contribution of cis-regulatory variants to pig phenotypes is poorly understood. Therefore, our aim was to elucidate the molecular mechanisms by which non-coding variants cause phenotypic differences in pigs by combining evolutionary biology analyses and functional genomics. RESULTS: We obtained a high-resolution phased chromosome-scale reference genome with a contig N50 of 18.03 Mb for the Luchuan pig breed (a representative eastern breed) and profiled potential selective sweeps in eastern and western pigs by resequencing the genomes of 234 pigs. Multi-tissue transcriptome and chromatin accessibility analyses of these regions suggest that tissue-specific selection pressure is mediated by promoters and distal cis-regulatory elements. Promoter variants that are associated with increased expression of the lysozyme (LYZ) gene in the small intestine might enhance the immunity of the gastrointestinal tract and roughage tolerance in pigs. In skeletal muscle, an enhancer-modulating single-nucleotide polymorphism that is associated with up-regulation of the expression of the troponin C1, slow skeletal and cardiac type (TNNC1) gene might increase the proportion of slow muscle fibers and affect meat quality. CONCLUSIONS: Our work sheds light on the molecular mechanisms by which non-coding variants shape phenotypic differences in pigs and provides valuable resources and novel perspectives to dissect the role of gene regulatory evolution in animal domestication and breeding.

Exploring Multi-Tissue Alternative Splicing and Skeletal Muscle Metabolism Regulation in Obese- and Lean-Type Pigs.

Alternative splicing (AS) is a crucial mechanism in post-transcriptional regulation, contributing significantly to the diversity of the transcriptome and proteome. In this study, we performed a comprehensive AS profile in nine tissues obtained from Duroc (lean-type) and Luchuan (obese-type) pigs. Notably, 94,990 AS events from 14,393 genes were identified. Among these AS events, it was observed that 80% belonged to the skipped exon (SE) type. Functional enrichment analysis showed that genes with more than ten AS events were closely associated with tissue-specific functions. Additionally, the analysis of overlap between differentially alternative splicing genes (DSGs) and differentially expressed genes (DEGs) revealed the highest number of overlapped genes in the heart and skeletal muscle. The novelty of our study is that it identified and validated three genes (PYGM, MAPK11 and CAMK2B) in the glucagon signaling pathway, and their alternative splicing differences were highly significant across two pig breeds. In conclusion, our study offers novel insights into the molecular regulation of diverse tissue physiologies and the phenotypic differences between obese- and lean-type pigs, which are helpful for pig breeding.

Regulation of myo-miR-24-3p on the Myogenesis and Fiber Type Transformation of Skeletal Muscle.

Skeletal muscle plays critical roles in providing a protein source and contributing to meat production. It is well known that microRNAs (miRNAs) exert important effects on various biological processes in muscle, including cell fate determination, muscle fiber morphology, and structure development. However, the role of miRNA in skeletal muscle development remains incompletely understood. In this study, we observed a critical miRNA, miR-24-3p, which exhibited higher expression levels in Tongcheng (obese-type) pigs compared to Landrace (lean-type) pigs. Furthermore, we found that miR-24-3p was highly expressed in the dorsal muscle of pigs and the quadriceps muscle of mice. Functionally, miR-24-3p was found to inhibit proliferation and promote differentiation in muscle cells. Additionally, miR-24-3p was shown to facilitate the conversion of slow muscle fibers to fast muscle fibers and influence the expression of GLUT4, a glucose transporter. Moreover, in a mouse model of skeletal muscle injury, we demonstrated that overexpression of miR-24-3p promoted rapid myogenesis and contributed to skeletal muscle regeneration. Furthermore, miR-24-3p was found to regulate the expression of target genes, including Nek4, Pim1, Nlk, Pskh1, and Mapk14. Collectively, our findings provide evidence that miR-24-3p plays a regulatory role in myogenesis and fiber type conversion. These findings contribute to our understanding of human muscle health and have implications for improving meat production traits in livestock.

Identification of Differentially Expressed Genes in the Longissimus Dorsi Muscle of Luchuan and Duroc Pigs by Transcriptome Sequencing.

The Duroc pig originated in the United States and is a typical lean-meat pig. The breed grows fast, and the body size is large, but the meat quality is poor. The Luchuan pig is one of eight excellent local breeds in China; it has tender meat but is small in size. To study the factors that determine growth, we selected the longissimus dorsi muscle of Luchuan and Duroc pigs for transcriptome sequencing. The results of the transcriptome showed that 3682 genes were differentially expressed (DEGs) in the longissimus dorsi muscle of Duroc and Luchuan pigs. We screened out genes related to muscle development and selected the MYL2 (Myosin light chain-2) gene to perform preliminary research. Gene Ontology (GO) enrichment of biological functions and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the gene products were mainly involved in the Akt/FoxO signaling pathway, fatty acid metabolism, arachidonic acid metabolism and glycine, serine and threonine metabolism. Such pathways contributed to skeletal muscle growth, fatty acid metabolism and intramuscular fat deposition. These results provide insight into the mechanisms underlying the formation of skeletal muscle and provide candidate genes to improve growth traits, as well as contribute to improving the growth and development traits of pigs through molecular breeding.

Integrative metabolomic and transcriptomic analysis reveals difference in glucose and lipid metabolism in the longissimus muscle of Luchuan and Duroc pigs.

Luchuan pig, an obese indigenous Chinese porcine breed, has a desirable meat quality and reproductive capacity. Duroc, a traditional western breed, shows a faster growth rate, high feed efficiency and high lean meat rate. Given the unique features these two porcine breeds have, it is of interest to investigate the underlying molecular mechanisms behind their distinctive nature. In this study, the metabolic and transcriptomic profiles of longissimus dorsi muscle from Duroc and Luchuan pigs were compared. A total of 609 metabolites were identified, 77 of which were significantly decreased in Luchuan compared to Duroc, and 71 of which were significantly elevated. Most differentially accumulated metabolites (DAMs) upregulated in Luchuan were glycerophospholipids, fatty acids, oxidized lipids, alcohols, and amines, while metabolites downregulated in Luchuan were mostly amino acids, organic acids and nucleic acids, bile acids and hormones. From our RNA-sequencing (RNA-seq) data we identified a total of 3638 differentially expressed genes (DEGs), 1802 upregulated and 1836 downregulated in Luchuan skeletal muscle compared to Duroc. Combined multivariate and pathway enrichment analyses of metabolome and transcriptome results revealed that many of the DEGs and DAMs are associated with critical energy metabolic pathways, especially those related to glucose and lipid metabolism. We examined the expression of important DEGs in two pathways, AMP-activated protein kinase (AMPK) signaling pathway and fructose and mannose metabolism, using Real-Time Quantitative Reverse Transcription PCR (qRT-PCR). Genes related to glucose uptake, glycolysis, glycogen synthesis, fatty acid synthesis (PFKFB1, PFKFB4, MPI, TPI1, GYS1, SLC2A4, FASN, IRS1, ULK1) are more activated in Luchuan, while genes related to fatty acid oxidation, cholesterol synthesis (CPT1A, HMGCR, FOXO3) are more suppressed. Energy utilization can be a decisive factor to the distinctive metabolic, physiological and nutritional characteristics in skeletal muscle of the two breeds we studied. Our research may facilitate future porcine breeding projects and can be used to reveal the potential molecular basis of differences in complex traits between various breeds.

Genome-Wide Association Studies, Runs of Homozygosity Analysis, and Copy Number Variation Detection to Identify Reproduction-Related Genes in Bama Xiang Pigs.

Litter size and teat number are economically important traits in the porcine industry. However, the genetic mechanisms influencing these traits remain unknown. In this study, we analyzed the genetic basis of litter size and teat number in Bama Xiang pigs and evaluated the genomic inbreeding coefficients of this breed. We conducted a genome-wide association study to identify runs of homozygosity (ROH), and copy number variation (CNV) using the novel Illumina PorcineSNP50 BeadChip array in Bama Xiang pigs and annotated the related genes in significant single nucleotide polymorphisms and common copy number variation region (CCNVR). We calculated the ROH-based genomic inbreeding coefficients (F ROH) and the Spearman coefficient between F ROH and reproduction traits. We completed a mixed linear model association analysis to identify the effect of high-frequency copy number variation (HCNVR; over 5%) on Bama Xiang pig reproductive traits using TASSEL software. Across eight chromosomes, we identified 29 significant single nucleotide polymorphisms, and 12 genes were considered important candidates for litter-size traits based on their vital roles in sperm structure, spermatogenesis, sperm function, ovarian or follicular function, and male/female infertility. We identified 9,322 ROHs; the litter-size traits had a significant negative correlation to F ROH. A total of 3,317 CNVs, 24 CCNVR, and 50 HCNVR were identified using cnvPartition and PennCNV. Eleven genes related to reproduction were identified in CCNVRs, including seven genes related to the testis and sperm function in CCNVR1 (chr1 from 311585283 to 315307620). Two candidate genes (NEURL1 and SH3PXD2A) related to reproduction traits were identified in HCNVR34. The result suggests that these genes may improve the litter size of Bama Xiang by marker-assisted selection. However, attention should be paid to deter inbreeding in Bama Xiang pigs to conserve their genetic diversity.

Effects of the Probiotic, Lactobacillus delbrueckii subsp. bulgaricus, as a Substitute for Antibiotics on the Gastrointestinal Tract Microbiota and Metabolomics Profile of Female Growing-Finishing Pigs.

Lactobacillus delbrueckii subsp. bulgaricus (LDB) is an approved feed additive on the Chinese 'Approved Feed Additives' list. However, the possibility of LDB as an antibiotic replacement remains unclear. Particularly, the effect of LDB on microbiota and metabolites in the gastrointestinal tract (GIT) requires further explanation. This study aimed to identify the microbiota and metabolites present in fecal samples and investigate the relationship between the microbiota and metabolites to evaluate the potential of LDB as an antibiotic replacement in pig production. A total of 42 female growing-finishing pigs were randomly allocated into the antibiotic group (basal diet + 75 mg/kg aureomycin) and LDB (basal diet + 3.0 × 109 cfu/kg LDB) groups. Fecal samples were collected on days 0 and 30. Growth performance was recorded and assessed. 16S rRNA sequencing and liquid chromatography-mass spectrometry-based non-targeted metabolomics approaches were used to analyze the differences in microbiota and metabolites. Associations between the differences were calculated using Spearman correlations with the Benjamini−Hochberg adjustment. The LDB diet had no adverse effect on feed efficiency but slightly enhanced the average daily weight gain and average daily feed intake (p > 0.05). The diet supplemented with LDB increased Lactobacillus abundance and decreased that of Prevotellaceae_NK3B31_group spp. Dietary-supplemented LDB enhanced the concentrations of pyridoxine, tyramine, D-(+)-pyroglutamic acid, hypoxanthine, putrescine and 5-hydroxyindole-3-acetic acid and decreased the lithocholic acid concentration. The Lactobacillus networks (Lactobacillus, Peptococcus, Ruminococcaceae_UCG-004, Escherichia-Shigella, acetophenone, tyramine, putrescine, N-methylisopelletierine, N1-acetylspermine) and Prevotellaceae_NK3B31_group networks (Prevotellaceae_NK3B31_group, Treponema_2, monolaurin, penciclovir, N-(5-acetamidopentyl)acetamide, glycerol 3-phosphate) were the most important in the LDB effect on pig GIT health in our study. These findings indicate that LDB may regulate GIT function through the Lactobacillus and Prevotellaceae_NK3B31_group networks. However, our results were restrained to fecal samples of female growing-finishing pigs; gender, growth stages, breeds and other factors should be considered to comprehensively assess LDB as an antibiotic replacement in pig production.

SQLE Promotes Differentiation and Apoptosis of Bovine Skeletal Muscle-Derived Mesenchymal Stem Cells.

In this study, Squalene epoxidase (SQLE) overexpression vector was transfected into bovine skeletal muscle-derived mesenchymal stem/stromal cells (MSCs) to study the molecular mechanism of SQLE regulating meat quality through myogenesis. We initially profiled the expression of SQLE in cattle embryos and adults, in the muscle tissue of four different cattle varieties, and in 11 different tissues/organs of Guangxi cattle variety. Subsequently, we isolated and cultured bovine skeletal muscle-derived MSCs and detected the expression of SQLE during cell proliferation and differentiation. Then, we constructed a bovine SQLE overexpression vector and transfected it into bovine skeletal muscle-derived MSCs by liposome transfection. Cell Counting Kit-8 (CCK-8), 5-ethynyl-20-deoxyuridine (EdU), flow cytometry, immunofluorescence, and quantitative polymerase chain reaction assays were used to characterize cell proliferation and differentiation in detail. The results showed that the relative expression level of bovine SQLE gene in brain tissue was the highest, and in adult muscle tissue was significantly higher than that in embryonic stage. Especially, the expression of SQLE was significantly upregulated in cell differentiation stage. Furthermore, the proliferation, cell cycle, apoptosis, and myoblast differentiation assays indicated that SQLE significantly promoted the differentiation and apoptosis of bovine skeletal muscle-derived MSCs, but inhibited their proliferation. In conclusion, our study reveals the role of SQLE in myoblast differentiation. These results will provide new clues for the regulation network of bovine muscle development.

Two RNAs or DNAs May Artificially Fuse Together at a Short Homologous Sequence (SHS) during Reverse Transcription or Polymerase Chain Reactions, and Thus Reporting an SHS-Containing Chimeric RNA Requires Extra Caution.

Tens of thousands of chimeric RNAs have been reported. Most of them contain a short homologous sequence (SHS) at the joining site of the two partner genes but are not associated with a fusion gene. We hypothesize that many of these chimeras may be technical artifacts derived from SHS-caused mis-priming in reverse transcription (RT) or polymerase chain reactions (PCR). We cloned six chimeric complementary DNAs (cDNAs) formed by human mitochondrial (mt) 16S rRNA sequences at an SHS, which were similar to several expression sequence tags (ESTs).These chimeras, which could not be detected with cDNA protection assay, were likely formed because some regions of the 16S rRNA are reversely complementary to another region to form an SHS, which allows the downstream sequence to loop back and anneal at the SHS to prime the synthesis of its complementary strand, yielding a palindromic sequence that can form a hairpin-like structure.We identified a 16S rRNA that ended at the 4th nucleotide(nt) of the mt-tRNA-leu was dominant and thus should be the wild type. We also cloned a mouse Bcl2-Nek9 chimeric cDNA that contained a 5-nt unmatchable sequence between the two partners, contained two copies of the reverse primer in the same direction but did not contain the forward primer, making it unclear how this Bcl2-Nek9 was formed and amplified. Moreover, a cDNA was amplified because one primer has 4 nts matched to the template, suggesting that there may be many more artificial cDNAs than we have realized, because the nuclear and mt genomes have many more 4-nt than 5-nt or longer homologues. Altogether, the chimeric cDNAs we cloned are good examples suggesting that many cDNAs may be artifacts due to SHS-caused mis-priming and thus greater caution should be taken when new sequence is obtained from a technique involving DNA polymerization.

Cloning of Porcine Pituitary Tumor Transforming Gene 1 and Its Expression in Porcine Oocytes and Embryos.

The maternal-to-embryonic transition (MET) is a complex process that occurs during early mammalian embryogenesis and is characterized by activation of the zygotic genome, initiation of embryonic transcription, and replacement of maternal mRNA with embryonic mRNA. The objective of this study was to reveal the temporal expression and localization patterns of PTTG1 during early porcine embryonic development and to establish a relationship between PTTG1 and the MET. To achieve this goal, reverse transcription-polymerase chain reaction (RT-PCR) was performed to clone porcine PTTG1. Subsequently, germinal vesicle (GV)- and metaphase II (MII)-stage oocytes, zygotes, 2-, 4-, and 8-cell-stage embryos, morulas, and blastocysts were produced in vitro and their gene expression was analyzed. The results revealed that the coding sequence of porcine PTTG1 is 609-bp in length and that it encodes a 202-aa polypeptide. Using qRT-PCR, PTTG1 mRNA expression was observed to be maintained at high levels in GV- and MII-stage oocytes. The transcript levels in oocytes were also significantly higher than those in embryos from the zygote to blastocyst stages. Immunohistochemical analyses revealed that porcine PTTG1 was primarily localized to the cytoplasm and partially localized to the nucleus. Furthermore, the PTTG1 protein levels in MII-stage oocytes and zygotes were significantly higher than those in embryos from the 2-cell to blastocyst stage. After fertilization, the level of this protein began to decrease gradually until the blastocyst stage. The results of our study suggest that porcine PTTG1 is a new candidate maternal effect gene (MEG) that may participate in the processes of oocyte maturation and zygotic genome activation during porcine embryogenesis.

Factors affecting efficiency of introducing foreign DNA and RNA into parthenogenetic or in vitro-fertilized porcine eggs by cytoplasmic microinjection.

Cytoplasmic microinjection (CI) of foreign gene into in vivo fertilized zygotes has emerged as a useful tool for transgenic pig production. In the current study, we investigated factors affecting transgenic efficiency and developmental potential of parthenogenetic (PA) and in vitro-fertilized (IVF) porcine embryos produced by CI. These factors included adding of RNase inhibitor, DNA or RNA concentration, injection time, and different structures of plasmids. Our results showed that adding of 1-4 U/µL of RNase inhibitor did not have negative effect on development potential of CI-PA embryos, and RNase inhibitor injection significantly increased EGFP expressing rate of CI-PA embryos. High injection DNA concentration and long injection interval after PA significantly reduced blastocyst formation. Different molecular structures such as DNA or RNA affected CI-PA embryos development, and RNA had little harmful effect on pig's early embryonic development. EGFP expression rate of CI-IVF embryos was improved following the increase of foreign DNA concentration, but blastocyst formation rate was decreased. Injection time after IVF did not show any significant difference on embryonic development, but longer interval resulted in a significantly lower EGFP expressing rate. Cas9 mRNA and myostatin (GDF-8) sgRNA co-injection indicated that the mutation rate of CI-IVF group was significantly higher than that of CI-PA. The CI-IVF-generated embryos were then transferred to six recipient pigs, but no live piglets were obtained. The following pronuclear formation assessment showed more than 76.1% IVF zygotes were polyspermy. These results demonstrate that CI-PA and CI-IVF were effective methods for production of transgenic pig embryos. However, polyspermic fertilization and poor quality of porcine IVF blastocysts are still the main problem of resulting in pregnancy failure.

Successful cloning of an adult breeding boar from the novel Chinese Guike No. 1 swine specialized strain.

Somatic cloning, also known as somatic cell nuclear transfer (SCNT), is a promising technology which has been expected to rapidly extend the population of elaborately selected breeding boars with superior production performance. Chinese Guike No. 1 pig breed is a novel swine specialized strain incorporated with the pedigree background of Duroc and Chinese Luchuan pig breeds, thus inherits an excellent production performance. The present study was conducted to establish somatic cloning procedures of adult breeding boars from the Chinese Guike No. 1 specialized strain. Ear skin fibroblasts were first isolated from a three-year-old Chinese Guike No. 1 breeding boar, and following that, used as donor cell to produce nuclear transfer embryos. Such cloned embryos showed full in vitro development and with the blastocyst formation rate of 18.4 % (37/201, three independent replicates). Finally, after transferring of 1187 nuclear transfer derived embryos to four surrogate recipients, six live piglets with normal health and development were produced. The overall cloning efficiency was 0.5 % and the clonal provenance of such SCNT derived piglets was confirmed by DNA microsatellite analysis. All of the cloned piglets were clinically healthy and had a normal weight at 1 month of age. Collectively, the first successful cloning of an adult Chinese Guike No. 1 breeding boar may lay the foundation for future improving the pig production industry.

Culture at a Higher Temperature Mildly Inhibits Cancer Cell Growth but Enhances Chemotherapeutic Effects by Inhibiting Cell-Cell Collaboration.

Acute febrile infections have historically been used to treat cancer. To explore the underlying mechanism, we studied chronic effects of fever on cancer cell growth and chemotherapeutic efficacy in cell culture. We found that culturing cancer cells at 39°C mildly inhibited cell growth by arresting the cells at the G1 phase of the cell cycle. When cells were seeded in culture dishes at a lower density, e.g. about 1000-2000 cells per 35-mm dish, the growth inhibition was much greater, manifested as many fewer cell colonies in the 39°C dishes, compared with the results at a higher density seeding, e.g. 20,000 cells per dish, suggesting that cell-cell collaboration as the Allee effect in cell culture is inhibited at 39°C. Withdrawal of cells from serum enhanced the G1 arrest at 39°C and, for some cell lines such as A549 lung cancer cells, serum replenishment failed to quickly drive the cells from the G1 into the S and G2-M phases. Therapeutic effects of several chemotherapeutic agents, including clove bud extracts, on several cancer cell lines were more potent at 39°C than at 37°C, especially when the cells were seeded at a low density. For some cell lines and some agents, this enhancement is long-lasting, i.e. continuing after the cessation of the treatment. Collectively these results suggest that hyperthermia may inhibit cancer cell growth by G1 arrest and by inhibition of cell-cell collaboration, and may enhance the efficacy of several chemotherapeutic agents, an effect which may persist beyond the termination of chemotherapy.

Cyclin-dependent kinase 4 may be expressed as multiple proteins and have functions that are independent of binding to CCND and RB and occur at the S and G 2/M phases of the cell cycle.

Cyclin-dependent kinase 4 (CDK4) is known to be a 33 kD protein that drives G 1 phase progression of the cell cycle by binding to a CCND protein to phosphorylate RB proteins. Using different CDK4 antibodies in western blot, we detected 2 groups of proteins around 40 and 33 kD, respectively, in human and mouse cells; each group often appeared as a duplet or triplet of bands. Some CDK4 shRNAs could decrease the 33 kD wild-type (wt) CDK4 but increase some 40 kD proteins, whereas some other shRNAs had the opposite effects. Liquid chromatography-mass spectrometry/mass spectrometry analysis confirmed the existence of CDK4 isoforms smaller than 33 kD but failed to identify CDK4 at 40 kD. We cloned one CDK4 mRNA variant that lacks exon 2 and encodes a 26 kD protein without the first 74 amino acids of the wt CDK4, thus lacking the ATP binding sequence and the PISTVRE domain required for binding to CCND. Co-IP assay confirmed that this ΔE2 protein lost CCND1- and RB1-binding ability. Moreover, we found, surprisingly, that the wt CDK4 and the ΔE2 could inhibit G 1-S progression, accelerate S-G 2/M progression, and enhance or delay apoptosis in a cell line-specific manner in a situation where the cells were treated with a CDK4 inhibitor or the cells were serum-starved and then replenished. Hence, CDK4 seems to be expressed as multiple proteins that react differently to different CDK4 antibodies, respond differently to different shRNAs, and, in some situations, have previously unrecognized functions at the S-G 2/M phases of the cell cycle via mechanisms independent of binding to CCND and RB.