Ame-miR-1-3p of bee venom reduced cell viability through the AZIN1/OAZ1-ODC1-polyamines pathway and enhanced the defense ability of honeybee (Apis mellifera L.).

Bee venom serves as an essential defensive weapon for bees and also finds application as a medicinal drug. MicroRNAs (miRNAs) serve as critical regulators and have been demonstrated to perform a variety of biological functions. However, the presence of miRNAs in bee venom needs to be confirmed. Therefore, we conducted small RNA sequencing and identified 158 known miRNAs, 15 conserved miRNAs and 4 novel miRNAs. It is noteworthy that ame-miR-1-3p, the most abundant among them, accounted for over a quarter of all miRNA reads. To validate the function of ame-miR-1-3p, we screened 28 candidate target genes using transcriptome sequencing and three target gene prediction software (miRanda, PITA and TargetScan) for ame-miR-1-3p. Subsequently, we employed real-time quantitative reverse transcription PCR (qRT-PCR), Western blot and other technologies to confirm that ame-miR-1-3p inhibits the relative expression of antizyme inhibitor 1 (AZIN1) by targeting the 3' untranslated region (UTR) of AZIN1. This, in turn, caused ODC antizyme 1 (OAZ1) to bind to ornithine decarboxylase 1 (ODC1) and mark ODC1 for proteolytic destruction. The reduction in functional ODC1 ultimately resulted in a decrease in polyamine biosynthesis. Furthermore, we determined that ame-miR-1-3p accelerates cell death through the AZIN1/OAZ1-ODC1-polyamines pathway. Our studies demonstrate that ame-miR-1-3p diminishes cell viability and it may collaborate with sPLA2 to enhance the defence capabilities of honeybees (Apis mellifera L.). Collectively, these data further elucidate the defence mechanism of bee venom and expand the potential applications of bee venom in medical treatment.

Haplotype-resolved 3D chromatin architecture of the hybrid pig.

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Knockdown of RGMA improves ischemic stroke via Reprogramming of Neuronal Metabolism.

Neuronal energy metabolism dysregulation is involved in various pathologies of Ischemia-reperfusion (I/R), yet the role of RGMA in neuronal metabolic reprogramming has not been reported. In this study, we found that RGMA expression significantly increased after I/R, and compared to control mice, mice with MCAO/R showed an increase in glycolytic metabolic products and the expression of glycolytic pathway proteins. Furthermore, RGMA levels are closely related to neuronal energy metabolism. We discovered that knockdown of RGMA can shift neuronal energy metabolism towards oxidative phosphorylation and the pentose phosphate pathway, thereby protecting mice from ischemic reperfusion injury. Mechanistically, knockdown of RGMA can downregulate PGK1 expression, reducing the increase in glycolytic flux following ischemia reperfusion. Moreover, we found that knockdown of RGMA can reduce the interaction between USP10 and PGK1, thus affecting the ubiquitination degradation of PGK1. In summary, our data suggest that RGMA may regulate neuronal energy metabolism by inhibiting the USP10-mediated deubiquitination of PGK1, thus protecting it from I/R injury. This study provides new ideas for clarifying the intrinsic mechanism of neuronal damage after I/R.

An intronic enhancer of Cebpa regulates adipocyte differentiation and adipose tissue development via long-range loop formation.

Cebpa is a master transcription factor gene for adipogenesis. However, the mechanisms of enhancer-promoter chromatin interactions controlling Cebpa transcriptional regulation during adipogenic differentiation remain largely unknown. To reveal how the three-dimensional structure of Cebpa changes during adipogenesis, we generated high-resolution chromatin interactions of Cebpa in 3T3-L1 preadipocytes and 3T3-L1 adipocytes using circularized chromosome conformation capture sequencing (4C-seq). We revealed dramatic changes in chromatin interactions and chromatin status at interaction sites during adipogenic differentiation. Based on this, we identified five active enhancers of Cebpa in 3T3-L1 adipocytes through epigenomic data and luciferase reporter assays. Next, epigenetic repression of Cebpa-L1-AD-En2 or -En3 by the dCas9-KRAB system significantly down-regulated Cebpa expression and inhibited adipocyte differentiation. Furthermore, experimental depletion of cohesin decreased the interaction intensity between Cebpa-L1-AD-En2 and the Cebpa promoter and down-regulated Cebpa expression, indicating that long-range chromatin loop formation was mediated by cohesin. Two transcription factors, RXRA and PPARG, synergistically regulate the activity of Cebpa-L1-AD-En2. To test whether Cebpa-L1-AD-En2 plays a role in adipose tissue development, we injected dCas9-KRAB-En2 lentivirus into the inguinal white adipose tissue (iWAT) of mice to suppress the activity of Cebpa-L1-AD-En2. Repression of Cebpa-L1-AD-En2 significantly decreased Cebpa expression and adipocyte size, altered iWAT transcriptome, and affected iWAT development. We identified functional enhancers regulating Cebpa expression and clarified the crucial roles of Cebpa-L1-AD-En2 and Cebpa promoter interaction in adipocyte differentiation and adipose tissue development.

Inhibition of GPR17/ID2 Axis Improve Remyelination and Cognitive Recovery after SAH by Mediating OPC Differentiation in Rat Model.

Myelin sheath injury contributes to cognitive deficits following subarachnoid hemorrhage (SAH). G protein-coupled receptor 17 (GPR17), a membrane receptor, negatively regulates oligodendrocyte precursor cell (OPC) differentiation in both developmental and pathological contexts. Nonetheless, GPR17's role in modulating OPC differentiation, facilitating remyelination post SAH, and its interaction with downstream molecules remain elusive. In a rat SAH model induced by arterial puncture, OPCs expressing GPR17 proliferated prominently by day 14 post-onset, coinciding with compromised myelin sheath integrity and cognitive deficits. Selective Gpr17 knockdown in oligodendrocytes (OLs) via adeno-associated virus (AAV) administration revealed that reduced GPR17 levels promoted OPC differentiation, restored myelin sheath integrity, and improved cognitive deficits by day 14 post-SAH. Moreover, GPR17 knockdown attenuated the elevated expression of the inhibitor of DNA binding 2 (ID2) post-SAH, suggesting a GPR17-ID2 regulatory axis. Bi-directional modulation of ID2 expression in OLs using AAV unveiled that elevated ID2 counteracted the restorative effects of GPR17 knockdown. This resulted in hindered differentiation, exacerbated myelin sheath impairment, and worsened cognitive deficits. These findings highlight the pivotal roles of GPR17 and ID2 in governing OPC differentiation and axonal remyelination post-SAH. This study positions GPR17 as a potential therapeutic target for SAH intervention. The interplay between GPR17 and ID2 introduces a novel avenue for ameliorating cognitive deficits post-SAH.

Research Note: Spermatozoa proteins identification in domesticated pigeons by proteomic analysis.

Proteins are considered major effectors of sperm function. However, the proteins expressed in pigeon sperm have not been explored. Here, we collected semen from meat and racing pigeons using the electroejaculation method and identified proteins in pigeon sperm using the proteomics approach. A total of 1,641 proteins were identified in the sperm of domesticated pigeons. Of which, 1,541 proteins were reliably quantified, and gene ontology (GO) and associated bioinformatics analyses indicated that annotated proteins were linked to the oxidation-reduction process, integral component membrane, and protein binding, etc. Among quantified proteins, 1,515 and 1,507 proteins were respectively presented in White King pigeons and racing pigeons, and 1,481 proteins were shared between these 2 types of pigeons, including axonemal dynein, solute carrier, cilia- and flagella-associated protein, outer dense fiber protein, etc. Proteins in our constructed protein-protein interaction (PPI) network are involved in oxidative phosphorylation, sperm axoneme assembly, cilium-dependent cell motility, axonemal dynein complex assembly, flagellated sperm motility, etc. In conclusion, this study characterized the sperm proteome of pigeons and provided a foundation for the subsequent research screening markers for fertility evaluation of pigeons.

Integrated Analysis of Immune Infiltration and Hub Pyroptosis-Related Genes for Multiple Sclerosis.

Purpose: Studies on overall immune infiltration and pyroptosis in patients with multiple sclerosis (MS) are limited. This study explored immune cell infiltration and pyroptosis in MS using bioinformatics and experimental validation. Methods: The GSE131282 and GSE135511 microarray datasets including brain autopsy tissues from controls and MS patients were downloaded for bioinformatic analysis. The gene expression-based deconvolution method, CIBERSORT, was used to determine immune infiltration. Differentially expressed genes (DEGs) and functional enrichments were analyzed. We then extracted pyroptosis-related genes (PRGs) from the DEGs by using machine learning strategies. Their diagnostic ability for MS was evaluated in both the training set (GSE131282 dataset) and validation set (GSE135511 dataset). In addition, messenger RNA (mRNA) expression of PRGs was validated using quantitative real-time polymerase chain reaction (qRT-PCR) in cortical tissue from an experimental autoimmune encephalomyelitis (EAE) model of MS. Moreover, the functional enrichment pathways of each hub PRG were estimated. Finally, co-expressed competitive endogenous RNA (ceRNA) networks of PRGs in MS were constructed. Results: Among the infiltrating cells, naive CD4+ T cells (P=0.006), resting NK cells (P=0.002), activated mast cells (P=0.022), and neutrophils (P=0.002) were significantly higher in patients with MS than in controls. The DEGs of MS were screened. Analysis of enrichment pathways showed that the pathways of transcriptional regulatory mechanisms and ion channels associating with pyroptosis. Four PRGs genes CASP4, PLCG1, CASP9 and NLRC4 were identified. They were validated in both the GSE135511 dataset and the EAE model by using qRT-PCR. CASP4 and NLRC4 were ultimately identified as stable hub PRGs for MS. Single-gene Gene Set Enrichment Analysis showed that they mainly participated in biosynthesis, metabolism, and organism resistance. ceRNA networks containing CASP4 and NLRC4 were constructed. Conclusion: MS was associated with immune infiltration. CASP4 and NLRC4 were key biomarkers of pyroptosis in MS.

Corrigendum: Case report and literature analysis: Autoimmune cerebellar ataxia associated with homer-3 antibodies.

[This corrects the article DOI: 10.3389/fneur.2022.951659.].

Cooked pork-derived exosome nanovesicles mediate metabolic disorder-microRNA could be the culprit.

In this study, exosomes from cooked meat were extracted by ultra-high-speed centrifugation. Approximately 80% of exosome vesicles were within 20-200 nm. In addition, the surface biomarkers of isolated exosomes were evaluated using flow cytometry. Further studies showed the exosomal microRNA profiles were different among cooked porcine muscle, fat and liver. Cooked pork-derived exosomes were chronically administered to ICR mice by drinking for 80 days. The mice plasma levels of miR-1, miR-133a-3p, miR-206 and miR-99a were increased to varying degrees after drinking exosome enriched water. Furthermore, GTT and ITT results confirmed an abnormal glucose metabolism and insulin resistance in mice. Moreover, the lipid droplets were significantly increased in the mice liver. A transcriptome analysis performed with mice liver samples identified 446 differentially expressed genes (DEGs). Functional enrichment analysis found that DEGs were enriched in metabolic pathways. Overall, the results suggest that microRNAs derived form cooked pork may function as a critical regulator of metabolic disorder in mice.

Whole-genome sequence association study identifies cyclin dependent kinase 8 as a key gene for the number of mummified piglets.

OBJECTIVE: Pigs, an ideal biomedical model for human diseases, suffer from about 50% early embryonic and fetal death, a major cause of fertility loss worldwide. However, identifying the causal variant remains a huge challenge. This study aimed to detect single nucleotide polymorphisms (SNPs) and candidate genes for the number of mummified (NM) piglets using the imputed whole-genome sequence (WGS) and validate the potential candidate genes. METHODS: The imputed WGS was introduced from genotyping-by-sequencing (GBS) using a multi-breed reference population. We performed genome-wide association studies (GWAS) for NM piglets at birth from a Landrace pig populatiGWAS peak located on SSC11: 0.10 to 7.11 Mbp (Top SNP, SSC11:1,889,658 bp; p = 9.98E-13) was identified in cyclin dependent kinase on. A total of 300 Landrace pigs were genotyped by GBS. The whole-genome variants were imputed, and 4,252,858 SNPs were obtained. Various molecular experiments were conducted to determine how the genes affected NM in pigs. RESULTS: A strong GWAS peak located on SSC11: 0.10 to 7.11 Mbp (Top SNP, SSC11:1,889,658 bp; p = 9.98E-13) was identified in cyclin dependent kinase 8 (CDK8) gene, which plays a crucial role in embryonic retardation and lethality. Based on the molecular experiments, we found that Y-box binding protein 1 (YBX1) was a crucial transcription factor for CDK8, which mediated the effect of CDK8 in the proliferation of porcine ovarian granulosa cells via transforming growth factor beta/small mother against decapentaplegic signaling pathway, and, as a consequence, affected embryo quality, indicating that this pathway may be contributing to mummified fetal in pigs. CONCLUSION: A powerful imputation-based association study was performed to identify genes associated with NM in pigs. CDK8 was suggested as a functional gene for the proliferation of porcine ovarian granulosa cells, but further studies are required to determine causative mutations and the effect of loci on NM in pigs.

Differential Expression Analysis of tRNA-Derived Small RNAs from Subcutaneous Adipose Tissue of Obese and Lean Pigs.

Epigenetic factors, including non-coding RNA regulation, play a vital role in the development of obesity and have been well researched. Transfer RNA-derived small RNA (tsRNA) is a class of non-coding RNA proven to be involved in various aspects of mammalian biology. Here we take pigs as a model for obesity research and use tsRNA-seq to investigate the difference in tsRNA expression in the subcutaneous adipose tissue of obese and lean pigs to elucidate the role of tsRNA in obesity development. A total of 482 tsRNAs were identified in pig adipose tissue, of which 123 were significantly differentially accumulated tsRNAs compared with the control group. The tRF-5c was the main type of these tsRNAs. The largest number of tsRNAs produced was the Gly-carrying tRNA, which produced 81 tsRNAs. Functional enrichment analysis revealed that differential tsRNAs indirectly participated in MAPK, AMPK, insulin resistance, the TNF signaling pathway, adipocytokine signaling pathway, and other signaling pathways by interacting with target genes. These are involved in bioenergetic metabolic regulatory processes, suggesting that tsRNAs may influence these pathways to mediate the regulation of energy metabolism in porcine adipocytes to promote lipid deposition, thus contributing to obesity. Our findings suggest a potential function of tsRNA in regulating obesity development.

Analysis of mRNA and lncRNA Expression Profiles of Breast Muscle during Pigeon (Columbalivia) Development.

The breast muscle is essential for flight and determines the meat yield and quality of the meat type in pigeons. At present, studies about long non-coding RNA (lncRNA) expression profiles in skeletal muscles across the postnatal development of pigeons have not been reported. Here, we used transcriptome sequencing to examine the White-King pigeon breast muscle at four different ages (1 day, 14 days, 28 days, and 2 years old). We identified 12,918 mRNAs and 9158 lncRNAs (5492 known lncRNAs and 3666 novel lncRNAs) in the breast muscle, and 7352 mRNAs and 4494 lncRNAs were differentially expressed in the process of development. We found that highly expressed mRNAs were mainly related to cell-basic and muscle-specific functions. Differential expression and time-series analysis showed that differentially expressed genes were primarily associated with muscle development and functions, blood vessel development, cell cycle, and energy metabolism. To further predict the possible role of lncRNAs, we also conducted the WGCNA and trans/cis analyses. We found that differentially expressed lncRNAs such as lncRNA-LOC102093252, lncRNA-G12653, lncRNA-LOC110357465, lncRNA-G14790, and lncRNA-LOC110360188 might respectively target UBE2B, Pax7, AGTR2, HDAC1, Sox8 and participate in the development of the muscle. Our study provides a valuable resource for studying the lncRNAs and mRNAs of pigeon muscles and for improving the understanding of molecular mechanisms in muscle development.

Plasma Metabolomic Profiling Reveals Preliminary Biomarkers of Pork Quality Based on pH Value.

This study aimed to identify biomarkers for pork quality evaluation. Firstly, the correlation between indicators of pork quality evaluation was investigated. The pH of pork meat at 45 min post slaughter showed a significant negative correlation with meat color indicators (r: -0.4868--0.3040). Subsequently, porcine plasma samples were further divided into low pH (pH = 6.16 ± 0.22) or high pH (pH = 6.75 ± 0.08) groups. Plasma metabolites in both sample groups were investigated using untargeted metabolomics. In total, 90 metabolites were recognized as differential metabolites using partial least squares discriminant analysis. Pathway enrichment analysis indicated these differential metabolites were enriched in amino acid metabolism and energy metabolism. Correlation analysis revealed that creatinine, L-carnitine, D-sphingosine, citraconic acid, and other metabolites may constitute novel plasma biomarkers with the pH value of pork meat. The current study provides important insights into plasma biomarkers for predicting pork quality based on pH value.

Transcriptome analysis revealed the roles of long non-coding RNA and mRNA in the bursa of Fabricius during pigeon (Columba livia) development.

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.

Case report and literature analysis: Autoimmune cerebellar ataxia associated with homer-3 antibodies.

Objective: We present a case of autoimmune cerebellar ataxia (ACA) associated with Homer protein homolog 3 (Homer-3) antibodies. Then, a review of the literature was conducted to summarize its clinical spectrum to improve clinicians' understanding of this rare entity. Case presentation: A 25-year-old man suffered from the subacute onset of cerebellar ataxia and psychiatric symptoms with abnormalities in the cerebellum on initial brain MRI and Homer-3 antibodies titers of 1:100 in the serum. His neurological symptoms did not improve after intravenous methylprednisolone but significantly improved following plasma exchange with a modified Rankin Scale (mRS) score of 1. However, 5 months later, he experienced relapse during oral prednisone tapering with enhanced cerebellar lesions and obvious cerebellar atrophy on repeated MRI. Various immunomodulatory approaches, including corticosteroids and plasma exchange, were utilized with no improvement. Then rituximab was given for the first time to treat Homer-3 autoimmunity with partial improvement of symptoms. However, the patient remained profoundly disabled with an mRS score of 4. Conclusion: ACA associated with Homer-3 antibodies may have a suboptimal response to corticosteroid therapy. More intense immunotherapy such as rituximab may contribute to the improvement of cerebellar syndrome. Relapsing courses and presentation of cerebellar atrophy may suggest a poor prognosis in this entity.

Detecting the selection signatures in Chinese Duroc,Landrace, Yorkshire, Liangshan, and Qingyu pigs.

Here we used two kinds of chips data from 5 pig breeds, Chinese Duroc (DD), Landrace (LL), Yorkshire (YY), Liangshan (LS), and Qingyu pigs (QY) in China to identify genes which show evidence of selection during domestication. Four breed pairs, LS-YY, QY-YY, DD-YY, and LL-YY pair, were performed to detect selection signatures using the Fst method. Then we identified a list of genes that played key roles in domestication and artificial selection. For example, the PTPRM gene was shared in LS-YY, QY-YY, and DD-YY pairs and it regulates a variety of cellular processes including cell growth, differentiation as signaling molecules. The HACD3 gene was shared in QY-YY and DD-YY pairs, and the HACD3 protein is involved in the production of very long-chain fatty acids of different chain lengths. Besides, the MYH11 gene that related to muscle contraction was found in LS-YY and LL-YY pair. These results suggested that genes related to immunity, disease resistance, and metabolism were subjected to strong selection pressure in Chinese domestic pigs in the progress of domestication and evolution; however, genes related to appearance, production performance, and reproduction were undergone strong artificial selection in commercial pig breeds.

Identification and expression pattern analysis of miRNAs in pectoral muscle during pigeon (Columba livia) development.

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.

A combined GWAS approach reveals key loci for socially-affected traits in Yorkshire pigs.

Socially affected traits in pigs are controlled by direct genetic effects and social genetic effects, which can make elucidation of their genetic architecture challenging. We evaluated the genetic basis of direct genetic effects and social genetic effects by combining single-locus and haplotype-based GWAS on imputed whole-genome sequences. Nineteen SNPs and 25 haplotype loci are identified for direct genetic effects on four traits: average daily feed intake, average daily gain, days to 100 kg and time in feeder per day. Nineteen SNPs and 11 haplotype loci are identified for social genetic effects on average daily feed intake, average daily gain, days to 100 kg and feeding speed. Two significant SNPs from single-locus GWAS (SSC6:18,635,874 and SSC6:18,635,895) are shared by a significant haplotype locus with haplotype alleles 'GGG' for both direct genetic effects and social genetic effects in average daily feed intake. A candidate gene, MT3, which is involved in growth, nervous, and immune processes, is identified. We demonstrate the genetic differences between direct genetic effects and social genetic effects and provide an anchor for investigating the genetic architecture underlying direct genetic effects and social genetic effects on socially affected traits in pigs.

Genome-wide DNA methylation analysis in Chinese Chenghua and Yorkshire pigs.

BACKGROUND: The Chinese Chenghua pig (CHP) is a typical Chinese domestic fatty pig breed with superior meat quality characteristics, while the Yorkshire pig (YP) has the characteristics of fast growth and a high rate of lean meat. Long term natural selection and artificial selection resulted in great phenotypic differences between the two breeds, including growth, development, production performance, meat quality, and coat color. However, genome-wide DNA methylation differences between CHP and YP remain unclear. RESULTS: DNA methylation data were generated for muscle tissues of CHP and YP using reduced representation bisulfite sequencing (RRBS). In this study, a total of 2,416,211 CpG sites were identified. Besides, the genome-wide DNA methylation analysis revealed 722 differentially methylated regions (DMRs) and 466 differentially methylated genes (DMGs) in pairwise CHP vs. YP comparison. Six key genomic regions (Sus scrofa chromosome (SSC)1:253.47-274.23 Mb, SSC6:148.71-169.49 Mb, SSC7:0.25-9.86 Mb, SSC12:43.06-61.49 Mb, SSC14:126.43-140.95 Mb, and SSC18:49.17-54.54 Mb) containing multiple DMRs were identified, and differences of methylation patterns in these regions may be related to phenotypic differences between CHP and YP. Based on the functional analysis of DMGs, 8 DMGs (ADCY1, AGBL4, EXOC2, FUBP3, PAPPA2, PIK3R1, MGMT and MYH8) were considered as important candidate genes associated with muscle development and meat quality traits in pigs. CONCLUSIONS: This study explored the difference in meat quality between CHP and YP from the epigenetic point of view, which has important reference significance for the local pork industry and pork food processing.

Transcriptome analysis reveals key genes and pathways associated with piglet fetal mummification.

China has the largest pork consumption worldwide. However, the high incidence of piglet fetal mummification (3%-5%) is an important factor that causes the slow improvement of pig reproductive capacity, and the genetic mechanism is still unclear. This study aimed to identify candidate genes associated with piglet fetal mummification. RNA-seq technology was used to compare transcriptome profiling of blood from healthy and mummified piglets at different stages of pregnancy (35, 56, 77, and 98 days). A total of 137-420 differentially expressed genes (DEGs) were detected at each stage. Seven DEGs were significantly differentially expressed at various stages. IL-9R, TLR8, ABLIM3, FSH-α, ASCC1, PRKCZ, and GCK may play important roles in the course of piglet fetal mummification. The differential genes we identified between the groups were mainly enriched in immune and inflammation regulation, while others were mainly enriched in reproduction. Considering the function of candidate genes, IL-9R and TLR8 were suggested as the most promising candidate genes involved in mummified piglet traits. We speculate that during pregnancy, it may be the combined effects of the above-mentioned inflammation, immune response, and reproduction-related signaling pathways that affect the occurrence of mummified piglets and further affect pig reproduction.