Toll-like receptor (TLR) 9 is an endosomal receptor recognizing bacterial DNA/CpG-containing oligodeoxynucleotides (CpG ODN). Blocking CpG ODN/TLR9 activity represents a strategy for therapeutic prevention of immune system overactivation. Herein, we report that a synthetic peptide (SP) representing the leucine-rich repeat 11 subdomain of the human TLR9 extracellular domain could attenuate CpG ODN/TLR9 activity in RAW264.7 cells by binding to CpG ODN and decreasing its internalization. Our results demonstrate that preincubation with SP specifically inhibited CpG ODN- but not lipopolysaccharide (LPS)- and lipopeptide (PAM3CSK4)-stimulated TNF-α and IL-6 release. Preincubation of SP with CpG ODN dose-dependently decreased TLR9-driven phosphorylation of IκBα and ERK and activation of NF-κB/p65. Moreover, SP dose-dependently decreased FAM-labeled CpG ODN internalization, whereas non-labeled CpG ODN reversed the inhibition. The KD value of SP-CpG ODN binding was within the micromolar range. Our results demonstrated that SP was a specific inhibitor of CpG ODN/TLR9 activity via binding to CpG ODN, leading to reduced ODN internalization and decreased activation of subsequent pathways within cells. Thus, SP could be used as a potential CpG ODN antagonist to block TLR9 signaling.

Cypermethrin is one of the most highly effective synthetic pyrethroid insecticides. The toxicity of cypermethrin to the reproductive and nervous systems has been well studied. However, little is known about the toxic effect of cypermethrin on immune cells such as macrophages. Here, we investigated the cytotoxicity of cypermethrin on macrophages and the underlying molecular mechanisms. We found that cypermethrin reduced cell viability and induced apoptosis in RAW 264.7 cells. Cypermethrin also increased reactive oxygen species (ROS) production and DNA damage in a dose-dependent manner. Moreover, cypermethrin-induced G1 cell cycle arrest was associated with an enhanced expression of p21, wild-type p53, and down-regulation of cyclin D1, cyclin E and CDK4. In addition, cypermethrin treatment activated MAPK signal pathways by inducing c-Jun N-terminal kinase (JNK) and extracellular regulated protein kinases 1/2 ERK1/2 phosphorylation, and increased the cleaved poly ADP-ribose polymerase (PARP). Further, pretreatment with antioxidant N-acetylcysteine (NAC) effectively abrogated cypermethrin-induced cell cytotoxicity, G1 cell cycle arrest, DNA damage, PARP activity, and JNK and ERK1/2 activation. The specific JNK inhibitor (SP600125) and ERK1/2 inhibitor (PD98059) effectively reversed the phosphorylation level of JNK and ERK1/2, and attenuated the apoptosis. Taken together, these data suggested that cypermethrin caused immune cell death via inducing cell cycle arrest and apoptosis regulated by ROS-mediated JNK/ERK pathway.

Post-translational modifications have been identified to be of great importance in cancers and lysine acetylation, which can attract the multifunctional transcription factor BRD4, has been identified as a potential therapeutic target. In this paper, we identify that BRD4 has an important role in colorectal cancer; and that its inhibition substantially wipes out tumor cells. Treatment with inhibitor MS417 potently affects cancer cells, although such effects were not always outright necrosis or apoptosis. We report that BRD4 inhibition also limits distal metastasis by regulating several key proteins in the progression of epithelial-to-mesenchymal transition (EMT). This effect of BRD4 inhibitor is demonstrated via liver metastasis in animal model as well as migration and invasion experiments in vitro. Together, our results demonstrate a new application of BRD4 inhibitor that may be of clinical use by virtue of its ability to limit metastasis while also being tumorcidal.

The STAT4 gene encodes a transcriptional factor that transmits signals induced by several key cytokines which play important roles in the development of autoimmune diseases. The aim of this study was to explore the association of STAT4 polymorphism with Graves' disease (GD) and Hashimoto's thyroiditis (HT). A total of 1048 autoimmune thyroid diseases (AITDs) patients (693 with GD and 355 with HT) and 909 age- and gender-matched controls were examined. STAT4 polymorphisms (rs7574865/rs10181656/ rs7572482) were genotyped by multiplex polymerase chain reaction (PCR) and ligase detection reaction (LDR). The results indicated that the frequencies of rs7574865 genotypes in patients with GD differed significantly from the controls (p=0.028), the T allele frequency of GD patients was also significantly higher than the controls (p=0.020). The genotypes of rs10181656 differed significantly in GD patients from controls (p=0.012); G allele frequencies were significantly higher in AITD patients than the controls (p=0.014 and 0.031, respectively). The frequencies of haplotype GC with GD and HT patients were significantly lower than their controls (p=0.015 and 0.030, respectively). In contrast, the frequencies of haplotype TG with GD and HT patients were significantly higher than their controls (p=0.016 and 0.048, respectively). These findings strongly suggest that STAT4 rs7574865/rs10181656 polymorphisms increase the risk of AITD in a Chinese population.

Congenital cataract is the leading cause of blindness among children worldwide. Patients with posterior subcapsular congenital cataract (PSC) in the central visual axis can result in worsening vision and stimulus deprivation amblyopia. However, the pathogenesis of PSC remains unclear. This study aims to explore the functional regulation and mechanism of HTRA1 in human lens epithelial cells (HLECs). HTRA1 was significantly downregulated in the lens capsules of children with PSC compared to normal controls. HTRA1 is a suppression factor of transforming growth factor-β (TGF-β) signalling pathway, which plays a key role in cataract formation. The results showed that the TGF-β/Smad signalling pathway was activated in the lens tissue of PSC. The effect of HTRA1 on cell proliferation, migration and apoptosis was measured in HLECs. In primary HLECs, the downregulation of HTRA1 can promote the proliferation and migration of HLECs by activating the TGF-β/Smad signalling pathway and can significantly upregulate the TGF-β/Smad downstream target genes FN1 and α-SMA. HTRA1 was also knocked out in the eyes of C57BL/6J mice via adeno-associated virus-mediated RNA interference. The results showed that HTRA1 knockout can significantly upregulate p-Smad2/3 and activate the TGF-β/Smad signalling pathway, resulting in abnormal proliferation and irregular arrangement of lens epithelial cells and leading to the occurrence of subcapsular cataract. To conclude, HTRA1 was significantly downregulated in children with PSC, and the downregulation of HTRA1 enhanced the proliferation and migration of HLECs by activating the TGF-β/Smad signalling pathway, which led to the occurrence of PSC.

Bis-(2,3-dibromo-4,5-dihydroxy-phenyl)-methane (BDDPM) is a bromophenol first isolated from Rhodomelaceae confervoides. Our previous studies showed that BDDPM exerts PTP1B-inhibiting activity and anti-cancer activity against a wide range of tumor cells while it also showed lower cytotoxicity against normal cells. In the present study, we found that BDDPM exhibits significant activities toward angiogenesis in vitro. BDDPM inhibits multiple angiogenesis processes, including endothelial cell sprouting, migration, proliferation, and tube formation. Further kinase assays investigations found that BDDPM is a potent selective, but multi-target, receptor tyrosine kinase (RTKs) inhibitor. BDDPM (10 μM) inhibits the activities of fibroblast growth factor receptor 2 and 3 (FGFR2, 3), vascular endothelial growth factor receptor 2 (VEGFR2) and platelet-derived growth factor receptor α (PDGFRα) (inhibition rate: 57.7%, 78.6%, 78.5% and 71.1%, respectively). Moreover, BDDPM also decreases the phosphorylation of protein kinase B (PKB/Akt) and endothelial nitric oxide synthase (eNOS), as well as nitric oxide (NO) production in a dose dependent manner. These results indicate that BDDPM can be exploited as an anti-angiogenic drug, or as a lead compound for the development of novel multi-target RTKs inhibitors.

Cisatracurium besylate is an ideal non-depolarizing muscle relaxant which is widely used in clinical application. However, some studies have suggested that cisatracurium besylate can affect cell proliferation. Moreover, its specific mechanism of action remains unclear. Here, we found that the number of GFP-LC3 (green fluoresent protein-light chain 3) positive autophagosomes and the rate of mitochondria fracture both increased significantly in drug-treated GFP-LC3 and MitoDsRed stable HeLa cells. Moreover, cisatracurium promoted the co-localization of LC3 and mitochondria and induced formation of autolysosomes. Levels of mitochondrial proteins decreased, which were reversed by the lysosome inhibitor Bafinomycin A1. Similar results with evidence of dose-dependent effects were found in both HeLa and Human Umbilical Vein Endothelial Cells (HUVECs). Cisatracurium lowered HUVEC viability to 0.16 (OD490) at 100 µM and to 0.05 (OD490) after 48 h in vitro; it increased the cell death rate to 56% at 100 µM and to 60% after 24 h in a concentration- and time-dependent manner (p < 0.01). Cell proliferation decreased significantly by four fold in Atg5 WT (wildtype) MEF (mouse embryonic fibroblast) (p < 0.01) but was unaffected in Atg5 KO (Knockout) MEF, even upon treatment with a high dose of cisatracurium. Cisatracurium induced significant increase in cell death of wild-type MEFs even in the presence of the apoptosis inhibitor zVAD. Thus, we conclude that activation of both the autophagic cell death and cell apoptosis pathways contributes to cisatracurium-mediated cell injury.

Selenoproteins that contain selenocysteine (Sec) are found in all kingdoms of life. Although they constitute a small proportion of the proteome, selenoproteins play essential roles in many organisms. In photosynthetic eukaryotes, selenoproteins have been found in algae but are missing in land plants (embryophytes). In this study, we explored the evolutionary dynamics of Sec incorporation by conveying a genomic search for the Sec machinery and selenoproteins across Archaeplastida. We identified a complete Sec machinery and variable sizes of selenoproteomes in the main algal lineages. However, the entire Sec machinery was missing in the Bangiophyceae-Florideophyceae clade (BV) of Rhodoplantae (red algae) and only partial machinery was found in three species of Archaeplastida, indicating parallel loss of Sec incorporation in different groups of algae. Further analysis of genome and transcriptome data suggests that all major lineages of streptophyte algae display a complete Sec machinery, although the number of selenoproteins is low in this group, especially in subaerial taxa. We conclude that selenoproteins tend to be lost in Archaeplastida upon adaptation to a subaerial or acidic environment. The high number of redox-active selenoproteins found in some bloom-forming marine microalgae may be related to defense against viral infections. Some of the selenoproteins in these organisms may have been gained by horizontal gene transfer from bacteria.

As one of the largest transcription factor families in plants, bZIP transcription factors play important regulatory roles in different biological processes, especially in the process of stress response. Salt stress inhibits the growth and yield of sugar beet. However, bZIP-related studies in sugar beet (Beta vulgaris L.) have not been reported. This study aimed to identify the bZIP transcription factors in sugar beet and analyze their biological functions and response patterns to salt stress. Using bioinformatics, 48 BvbZIP genes were identified in the genome of sugar beet, encoding 77 proteins with large structural differences. Collinearity analysis showed that three pairs of BvbZIP genes were fragment replication genes. The BvbZIP genes were grouped according to the phylogenetic tree topology and conserved structures, and the results are consistent with those reported in Arabidopsis. Under salt stress, the expression levels of most BvbZIP genes were decreased, and only eight genes were up-regulated. GO analysis showed that the BvbZIP genes were mainly negatively regulated in stress response. Protein interaction prediction showed that the BvbZIP genes were mainly involved in light signaling and ABA signal transduction, and also played a certain role in stress responses. In this study, the structures and biological functions of the BvbZIP genes were analyzed to provide foundational data for further mechanistic studies and for facilitating the efforts toward the molecular breeding of stress-resilient sugar beet.

Dental caries, particularly secondary caries, which is the main contributor to dental repair failure, has been the subject of extensive research due to its biofilm-mediated, sugar-driven, multifactorial, and dynamic characteristics. The clinical utility of restorations is improved by cleaning bacteria nearby and remineralizing marginal crevices. In this study, a novel multifunctional dental resin composite (DRC) composed of Sr-N-co-doped titanium dioxide (Sr-N-TiO2) nanoparticles and nano-hydroxyapatite (n-HA) reinforcing fillers with improved antibacterial and mineralization properties is proposed. The experimental results showed that the anatase-phase Sr-N-TiO2 nanoparticles were synthesized successfully. After this, the curing depth (CD) of the DRC was measured from 4.36 ± 0.18 mm to 5.10 ± 0.19 mm, which met the clinical treatment needs. The maximum antibacterial rate against Streptococcus mutans (S. mutans) was 98.96%, showing significant inhibition effects (p < 0.0001), which was experimentally verified to be derived from reactive oxygen species (ROS). Meanwhile, the resin exhibited excellent self-remineralization behavior in an SBF solution, and the molar ratio of Ca/P was close to that of HA. Moreover, the relative growth rate (RGR) of mouse fibroblast L929 indicated a high biocompatibility, with the cytotoxicity level being 0 or I. Therefore, our research provides a suitable approach for improving the antibacterial and mineralization properties of DRCs.

Tomato is a crop that requires a sufficient supply of potassium (K) for optimal productivity and quality. K⁺-deficiency stress decreases tomato yield and quality. To further delve into the mechanism of the response to K⁺-deficiency and to screen out low-K⁺ tolerant genes in tomatoes, BGISEQ-500-based RNA sequencing was performed using two tomato genotypes (low-K⁺ tolerant JZ34 and low-K⁺ sensitive JZ18). We identified 1936 differentially expressed genes (DEGs) in JZ18 and JZ34 at 12 and 24 h after K⁺-deficiency treatment. According to the Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analyses, the DEGs that changed significantly primarily included transcription factors, transporters, kinases, oxidative stress proteins, and hormone signaling-and glycometabolism-related genes. The experimental results confirmed the induced expression of the responsive genes in the low-K⁺ signaling pathway. The largest group of DEGs comprised up to 110 oxidative stress-related genes. In total, 19 ethylene response factors (ERFs) demonstrated differential expression between JZ18 and JZ34 in response to K⁺-deficiency. Furthermore, we confirmed 20 DEGs closely related to K⁺-deficiency stress by quantitative RT-PCR (qRT-PCR), some of which affected the root configuration, these DEGs could be further studied for use as molecular targets to explore novel approaches, and to acquire more effective K acquisition efficiencies for tomatoes. A hypothesis involving possible cross-talk between phytohormone signaling cues and reactive oxygen species (ROS) leading to root growth in JZ34 is proposed. The results provide a comprehensive foundation for the molecular mechanisms involved in the response of tomatoes to low K⁺ stress.

Amur ide (Leuciscus waleckii), an important aquaculture species, inhabits neutral freshwater but can tolerate high salinity or alkalinity. As an extreme example, the population in Dali Nor lake inhabits alkalized soda water permanently, and migrates from alkaline water to neutral freshwater to spawn. In this study, we performed comparative transcriptome profiling study on the livers of Amur ide to interrogate the expression differences between the population that permanently inhabit freshwater in Ganggeng Nor lake (FW) and the spawning population that recently migrated from alkaline water into freshwater (SM). A total of 637,234,880 reads were generated, resulting in 53,440 assembled contigs that were used as reference sequences. Comparisons of these transcriptome files revealed 444 unigenes with significant differential expression (p-value ≤ 0.01, fold-change ≥ 2), including 246 genes that were up-regulated in SM and 198 genes that were up-regulated in FW. The gene ontology (GO) enrichment analysis and KEGG pathway analysis indicated that the mTOR signaling pathway, Janus kinase-signal transducer and activator of transcription (JAK-STAT) signaling pathway, and oxidative phosphorylation were highly likely to affect physiological changes during spawning migration. Overall, this study demonstrates that transcriptome changes played a role in Amur ide spawning migration. These results provide a foundation for further analyses on the physiological and molecular mechanisms underlying Amur ide spawning migration.

HtrA serine peptidase 2 (HtrA2), also named Omi, is a pro-apoptotic protein that exhibits dramatic changes in expression levels in a variety of disorders, including ischemia/reperfusion injury, cancer, and neurodegeneration. In our study, Omi/HtrA2 protein levels were high in the heart, brain, kidney and liver, with elevated heart/brain expression in aging mice. A similar expression pattern was observed at the mRNA level, which suggests that the regulation of Omi/HtrA2 is predominately transcriptional. Promoter binding by transcription factors is the main influencing factor of transcription, and to identify specific promoter elements that contribute to the differential expression of mouse Omi/HtrA2, we constructed truncated Omi/HtrA2 promoter/luciferase reporter vectors and analyzed their relative luciferase activity; it was greatest in the promoter regions at -1205~-838 bp and -146~+93 bp, with the -838~-649 bp region exhibiting negative regulatory activity. Bioinformatics analysis suggested that the Omi/HtrA2 gene promoter contains a CpG island at -709~+37 bp, and eight heat shock transcription factor 1 (HSF1) sites, two Sp1 transcription factor (SP1)sites, one activator protein (AP) site, seven p53 sites, and four YY1 transcription factor(YY1) sites were predicted in the core areas. Furthermore, we found that p53 and HSF1 specifically binds to the Omi/HtrA2 promoter using chromatin immunoprecipitation analysis. These results provide a foundation for understanding Omi/HtrA2 regulatory mechanisms, which could further understanding of HtrA-associated diseases.

Immortalized cell lines have been used in a wide range of applications in research on immune disorders and cellular metabolic regulation due to the stability and uniformity of their cellular characteristics. At present, the investigation into molecular functions and signaling pathways within bovine cells remains largely limited by the lack of immortalized model cells. Current methods for immortalizing bovine cells are mainly restricted to the ectopic expression of human telomerase reverse transcriptase (hTERT) through transient transfection or virus-mediated delivery, which have defects in efficiency and reliability. In this study, we identified bovine TERT (bTERT) as a novel potent biofactor for immortalizing bovine cells with great advantages over hTERT, and established an efficient and easily manipulated strategy for the immortalization of bovine primary cells. Through the homology-mediated end-joining-based insertion of bTERT at the ROSA26 locus, we successfully generated immortalized bovine fetal fibroblast cell lines with stable characteristics. The observed limitation of this strategy in immortalizing bovine bone marrow-derived macrophages was attributed to the post-translational modification of bTERT, causing inhibited nuclear localization and depressed activity of bTERT in this terminally differentiated cell. In summary, we constructed an innovative method to achieve the high-quality immortalization of bovine primary cells, thereby expanding the prospects for the future application of immortalized bovine model cell lines.

Aging affects the structure and function of the liver. Hydroxysafflor yellow A (HSYA) effectively improves liver aging (LA) in mice, but the potential mechanisms require further exploration. In this study, an integrated approach combining network pharmacology and transcriptomics was used to elucidate the potential mechanisms of HSYA delay of LA. The targets of HSYA were predicted using the PharmMapper, SwissTargetPrediction, and CTD databases, and the targets of LA were collected from the GeneCards database. An ontology (GO) analysis and a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation of genes related to HSYA delay of LA were performed using the DAVID database, and Cytoscape software was used to construct an HSYA target pathway network. The BMKCloud platform was used to sequence mRNA from mouse liver tissue, screen differentially expressed genes (DEGs) that were altered by HSYA, and enrich their biological functions and signaling pathways through the OmicShare database. The results of the network pharmacology and transcriptomic analyses were combined. Then, quantitative real-time PCR (qRT-PCR) and Western blot experiments were used to further verify the prediction results. Finally, the interactions between HSYA and key targets were assessed by molecular docking. The results showed that 199 potentially targeted genes according to network pharmacology and 480 DEGs according to transcriptomics were involved in the effects of HSYA against LA. An integrated analysis revealed that four key targets, including HSP90AA1, ATP2A1, NOS1 and CRAT, as well as their three related pathways (the calcium signaling pathway, estrogen signaling pathway and cGMP-PKG signaling pathway), were closely related to the therapeutic effects of HSYA. A gene and protein expression analysis revealed that HSYA significantly inhibited the expressions of HSP90AA1, ATP2A1 and NOS1 in the liver tissue of aging mice. The molecular docking results showed that HSYA had high affinities with the HSP90AA1, ATP2A1 and NOS1 targets. Our data demonstrate that HSYA may delay LA in mice by inhibiting the expressions of HSP90AA1, ATP2A1 and NOS1 and regulating the calcium signaling pathway, the estrogen signaling pathway, and the cGMP-PKG signaling pathway.

The success of root canal therapy depends mainly on the complete elimination of the root canal bacterial biofilm. The validity and biocompatibility of root canal disinfectant materials are imperative for the success of root canal treatment. However, the insufficiency of the currently available root canal disinfectant materials highlights that more advanced materials are still needed. In this study, a nanozyme-loaded hydrogel (Fe3O4-CaO2-Hydrogel) was modified and analyzed as a root canal disinfectant material. Fe3O4-CaO2-Hydrogel was fabricated and examined for its release profile, biocompatibility, and antibacterial activity against E. faecalis and S. sanguis biofilms in vitro. Furthermore, its efficiency in eliminating the root canal bacterial biofilm removal in SD rat teeth was also evaluated. The results in vitro showed that Fe3O4-CaO2-Hydrogel could release reactive oxygen species (ROS). Moreover, it showed good biocompatibility, disrupting bacterial cell membranes, and inhibiting exopolysaccharide production (p < 0.0001). In addition, in vivo results showed that Fe3O4-CaO2-Hydrogel strongly scavenged on root canal biofilm infection and prevented further inflammation expansion (p < 0.05). Altogether, suggesting that Fe3O4-CaO2-Hydrogel can be used as a new effective biocompatible root canal disinfectant material. Our research provides a broad prospect for clinical root canal disinfection, even extended to other refractory infections in deep sites.

Sipunculus nudus, an old marine species, has great potential for use as functional seafood due to its various bioactivities. Its potential antithrombotic activity pushed us to isolate the bio-active components bio-guided by tracking fibrinolytic activity. As a result, a novel protease named as SK (the kinase obtained from S. nudus) was obtained, which possessed a molecular weight of 28,003.67 Da and 15 N-terminal amino acid sequences of PFPVPDPFVWDTSFQ. SK exerted inhibitory effects on thrombus formation through improving the coagulation system with dose-effect relationship within a certain range. Furthermore, in most cases SK got obviously better effect than that of urokinase. With the help of untargeted mass spectrometry-based metabolomics profiling, arachidonic acid, sphingolipid, and nicotinate and nicotinamide mechanism pathways were found to be important pathways. They revealed that the effect mechanism of SK on common carotid arterial thrombosis induced by FeCl₃ was achieved by inhibiting vessel contraction, platelet aggregation, adhesion, and release, correcting endothelial cell dysfunction and retarding process of thrombus formation. This study demonstrated SK was a promising thrombolytic agent on the basis of its comprehensive activities on thrombosis, and it should get further exploitation and utilization.

Protein conjugations at post-translational levels are known to be essential to protein stability and function. Recently, it has been proven that the split protein CnaB2 (SpyTag/SpyCatcher, ST/SC) from Streptococcus pyogenes can induce covalent conjugation rapidly and efficiently under various conditions. The protein of interest fused with the split protein SC/ST could be assembled spontaneously. In light of this finding, we introduced the ST/SC protein coupling concept into the silkworm-bacmid protein expression system (SpyBEVS). As a proof of concept, we first examined and confirmed that a competent ligation occurred between ST/SC-fused protein partners in vitro in cultured silkworm cells and in vivo in silkworm larvae by co-infection of several recombinant baculoviruses. The protein conjugation could be also achieved sufficiently by a simple one-step mixture of purified ST/SC-tagged peptide-protein pairs in vitro. Given the flexibility and robustness of silkworm-BEVS, our results on SpyBEVS show an alternative method for enabling the production of protein decorations in vitro and inside of silkworms.

The Catharanthus roseus receptor-like kinase 1-like (CrRLK1L), which is a vital member of the plant receptor-like kinase family, plays versatile roles in plant growth, development, and stress response. Although the primary screening of tomato CrRLK1Ls has been reported previously, our knowledge of these proteins is still scarce. Using the latest genomic data annotations, a genome-wide re-identification and analysis of the CrRLK1Ls in tomatoes were conducted. In this study, 24 CrRLK1L members were identified in tomatoes and researched further. Subsequent gene structures, protein domains, Western blot analyses, and subcellular localization analyses all confirmed the accuracy of the newly identified SlCrRLK1L members. Phylogenetic analyses showed that the identified SlCrRLK1L proteins had homologs in Arabidopsis. Evolutionary analysis indicated that two pairs of the SlCrRLK1L genes had predicted segmental duplication events. Expression profiling analyses demonstrated that the SlCrRLK1L genes were expressed in various tissues, and most of them were up- or down-regulated by bacteria and PAMP treatments. Together, these results will lay the foundation for elaborating the biological roles of SlCrRLK1Ls in tomato growth, development, and stress response.

Pig diseases seriously threaten the health of pigs and the benefits of pig production. Previous research has indicated that Chinese native pigs, such as the Min (M) pig, has a better disease resistance ability than Large White (LW) pigs. However, the molecular mechanism of this resistance is still unclear. In our study, we used serum untargeted metabolomics and proteomics, interrogated to characterize differences in the molecular immunities between six resistant and six susceptible pigs raised in the same environment. A total of 62 metabolites were identified as being significantly exhibited in M and LW pigs. Ensemble feature selection (EFS) machine learning methods were used to predict biomarkers of metabolites and proteins, and the top 30 were selected and retained. Weighted gene co-expression network analysis (WGCNA) confirmed that four key metabolites, PC (18:1 (11 Z)/20:0), PC (14:0/P-18: 0), PC (18:3 (6 Z, 9 Z, 12 Z)/16:0), and PC (16:1 (9 Z)/22:2 (13 Z, 16 Z)), were significantly associated with phenotypes, such as cytokines, and different pig breeds. Correlation network analysis showed that 15 proteins were significantly correlated with the expression of both cytokines and unsaturated fatty acid metabolites. Quantitative trait locus (QTL) co-location analysis results showed that 13 of 15 proteins co-localized with immune or polyunsaturated fatty acid (PUFA)-related QTL. Moreover, seven of them co-localized with both immune and PUFA QTLs, including proteasome 20S subunit beta 8 (PSMB8), mannose binding lectin 1 (MBL1), and interleukin-1 receptor accessory protein (IL1RAP). These proteins may play important roles in regulating the production or metabolism of unsaturated fatty acids and immune factors. Most of the proteins could be validated with parallel reaction monitoring, which suggests that these proteins may play an essential role in producing or regulating unsaturated fatty acids and immune factors to cope with the adaptive immunity of different pig breeds. Our study provides a basis for further clarifying the disease resistance mechanism of pigs.