Kyoto Encyclopedia Of Genes And Genomes Pathway Analysis Research Articles

Artesunate inhibits airway remodeling in asthma via the MAPK signaling pathway.

Background: Artesunate (ART), is a semi-synthetic water-soluble artemisinin derivative extracted from the plant Artemisia annua, which is often used to treating malaria. In vivo and in vitro studies suggested it may help decrease inflammation and attenuate airway remodeling in asthma. However, its underlying mechanism of action is not elucidated yet. Herein, an attempt is made to investigate the ART molecular mechanism in treating asthma. Methods: The BALB/c female mice sensitized via ovalbumin (OVA) have been utilized to establish the asthma model, followed by carrying out ART interventions. Lung inflammation scores by Haematoxylin and Eosin (H&E), goblet cell hyperplasia grade by Periodic Acid-Schiff (PAS), and collagen fibers deposition by Masson trichrome staining have been utilized for evaluating how ART affected asthma. RNA-sequencing (RNA-seq) analyses were performed to identify differentially expressed genes (DEGs). The DEGs were analyzed by Gene Ontology (GO) terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and Protein-Protein interaction (PPI) function analyses. Hub clusters were found by Cytoscape MCODE. Subsequently, Real-Time quantitative PCR (RT-qPCR) verified the mRNA expression profiles of DEGs. Finally, immunohistochemistry (IHC) and western blots have validated the relevant genes and potential pathways. Results: ART considerably attenuated inflammatory cell infiltration, mucus secretion, and collagen fibers deposition. KEGG pathway analysis revealed that the ART played a protective role via various pathways including the mitogen-activated protein kinase (MAPK) pathway as one of them. Moreover, ART could alleviate the overexpression of found in inflammatory zone 1(FIZZ1) as revealed by IHC and Western blot analyses. ART attenuated OVA-induced asthma by downregulating phosphorylated p38 MAPK. Conclusion: ART exerted a protective function in a multitarget and multi-pathway on asthma. FIZZ1 was a possible target for asthma airway remodeling. The MARK pathway was one of the key pathways by which ART protected against asthma.

The NF-κB pathway plays a vital role in rat salivary gland atrophy model

ObjectiveThe aim of this study was to explore the histopathological and genetic changes in the submandibular glands after duct ligation and provide important clues to functional regeneration. DesignWe established a rat salivary gland duct ligation model and observed pathological changes in the rat submandibular gland on day 1 and weeks 1, 2, 3, and 4 using hematoxylin and eosin staining, Alcian blue–periodic acid Schiff staining, Masson staining, terminal deoxynucleotidyl transferase dUTP nick-end labeling (TUNEL), and immunohistochemical staining. RNA sequencing was performed on normal salivary glands and those from the ligation model after 1 week. Significantly differentially expressed genes were selected, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed. ResultsApoptosis levels and histological and functional KEGG pathway analyses showed that injury to the salivary gland after ligation gradually increased. The TGF-β pathway was activated and promoted fibrosis. RNA sequencing results and further verification of samples at week 1 showed that the NF-κB pathway plays a vital role in salivary gland atrophy. ConclusionsOur results detailed the pathological changes in the submandibular gland after ligation and the important functions of the NF-κB pathway.

Abstract P2-03-20: Clinical profiling and comprehensive analysis of candidate genes related to breast cancer estrogen receptor intratumour heterogeneity

Abstract 1. Background and Aim Breast cancer has the highest incidence among the world’s population. One third of patients diagnosed with early breast cancer progressed into metastatic outcomes due to the inherent heterogeneity and evolutionary features of tumors. The main clinical manifestation is the discordance of receptor expression in metastases. ESR1 the estrogen receptor (ER) encoding gene, has been proved affected ER expression. Nevertheless, the correlation between ER transformation and the somatic mutation profiles in breast cancer metastatic lesions remains unclear. Availability of advanced high-throughput technologies, and the development of bioinformatics tools has greatly accelerated our cognition of the molecular basis of cancer. In this study, we aimed to evaluate the frequency, clinical characteristics and prognostic value of ER receptor conversion between primary and first metastatic lesions in 115 patients whose primary ER status are positive. We also compared the results of NGS of 42 patients and explored latent genes related to ER discordance in expression. 2. Method Qualified metastatic tumor tissue sequencing information was available for 42 patients. Samples were accessed using the FoundationOne CDx assay (Foundation Medicine, Cambridge, MA, USA) panel including 733 tumor related genes. High frequency mutant genes were defined as the genes count more than 4 times in each group. All data were analyzed using SPSS 25.0 software (Chicago, IL, USA). Pearson’s chi-squared (χ2) test, Fisher’s exact test, were used to test the associations between different variables. The Cox proportional hazard models was used to evaluate the relationship between KMT2C mutation and Disease free survival (DFS) in the patients cohorts. The results with p < 0.05 were considered statistically significant. 3. Results and Conclusions All 115 patients (primary ER status is positive) were divided into two groups, 70(60.87%) patients ER status remained positive while 45(39.14%) patients ER status transformed into negative. The histological type (p=0.008) and DFS (p=0.004) were significantly different in two subgroups. There were no significant correlations between the age at the time of diagnosis, BMI classification, menopausal status, surgery, histological grade, tumor size, lymph node status, metastasis, Ki67, adjuvant radiotherapy, neoadjuvant, metastasis site before second biopsy, first treatment after recurrent and first PFS after rebiopsy. Univariable analysis of DFS proved metastasis to be distinct risk factor for poor survival (hazard ratio [HR] > 1, p = 0.001). By contrast, neoadjuvant proved to be protective factor for better survival (hazard ratio [HR] < 1, p < 0.001). The surgery, histological type, histological grade, tumorsize, metastasis, neoadjuvant, first treatment after recurrent, first PFS after rebiopsy were subsequently analyzed with the multivariable Cox analysis. Metastasis (HR>1, p=0.007) and neoadjuvant (HR<1, p<0.001) remained independent prognostic factors of DFS. We compared the differentially mutated genes in the ER differential expression BC patient groups. Mutations in TP53 (n=24) and PIK3CA (n=19) occurred most frequently in both groups and account for 57.14% and 49.23% separately. Notably, mutations in KMT2C (n=6) was detected in ER+ to – subgroup and accounted for 37.5%. According to the results of NGS, the distribution of KMT2C mutations was different among the two cohorts. We identified the mutation in codon 321 was the hot spot of KMT2C in our patient cohort. Next, we used the Kyoto Encyclopedia of Genes and Genomes (KEGG) database to analyze the enriched pathways of KMT2C. KEGG pathway analysis revealed that patients with KMT2C mutations harbored significantly more mutations in genes involved in the Ubiquitin mediated proteolysis and Lysine degradation signaling pathway. Table 1. Patients’ basic clinicopathological characteristics among two groups Table 2. Univariate and Multivariate Analysis Between Clinicopathological Characteristics and DFS of 115 BC patients Table 3. Mutated Genes Rank and patients number Citation Format: Xi Wang, Yongmei Yin, Ziyi Fu. Clinical profiling and comprehensive analysis of candidate genes related to breast cancer estrogen receptor intratumour heterogeneity [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P2-03-20.

Tumor necrosis factor-related lncRNAs predict prognosis and immunotherapy response for patients with lung adenocarcinoma.

Lung adenocarcinoma (LUAD) has become one of the most lethal cancers, for which the recurrence and survival rates remain unfavorable. The tumor necrosis factor (TNF) family is involved in tumorigenesis and tumor progression. Various long non-coding RNAs (lncRNAs) play important roles by mediating the TNF family in cancer. Therefore, this study aimed to construct a TNF-related lncRNA signature to predict prognosis and immunotherapy response in LUAD. The expression of TNF family members and their related lncRNAs in a total of 500 enrolled LUAD patients was collected from The Cancer Genome Atlas (TCGA). Univariate Cox and the least absolute shrinkage and selection operator (LASSO)-Cox analysis was used to construct a TNF family-related lncRNA prognostic signature. Kaplan-Meier (KM) survival analysis was used to evaluate survival status. The time-dependent area under the receiver operating characteristic (ROC) curve (AUC) values were used to assess the predictive value of the signature to 1-, 2-, and 3-year overall survival (OS). Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were applied to identify the signature-related biological pathways. Furthermore, tumor immune dysfunction and exclusion (TIDE) analysis was employed to evaluate immunotherapy response. A total of 8 TNF-related lncRNAs significantly associated with OS of LUAD patients were used to construct a TNF family-related lncRNA prognostic signature. According to risk score, these patients were divided into high- and low-risk subgroups. The KM survival analysis indicated that patients in the high-risk group showed significantly less favorable OS than that of low-risk group. The AUC values in predicting 1-, 2-, and 3-year OS were 0.740, 0.738, and 0.758, respectively. Moreover, the GO and KEGG pathway analyses demonstrated that these lncRNAs were closely involved in immune-related signaling pathways. The further TIDE analysis indicated that high-risk patients had a lower TIDE score than that of low-risk patients, indicating that high-risk patients may be appropriate candidates for immunotherapy. For the first time, this study constructed and validated a prognostic predictive signature of LUAD patients based on TNF-related lncRNAs, and the signature showed good performance to predict immunotherapy response. Therefore, this signature may provide new strategies for individualized treatment of LUAD patients.

Differential expression profile of microRNAs in the lung tissues of coal workers with pneumoconiosis and patients with silicosis.

The purpose of this study was to characterize the microRNA (miRNA) profile of the lung tissues from coal workers' pneumoconiosis (CWP) and silicosis and to analyze the changes in downstream genes, biological processes, and signaling pathways based on the differently expressed miRNAs. Lung tissues from three CWP patients, eight silicosis patients, and four healthy controls were collected and analyzed for their miRNA profiles using Affymetrix® GeneChip® miRNA Arrays. Differentially expressed miRNAs (DEMs) were identified between the different groups. The miRanda and TargetScan databases were used to predict the putative target genes, and volcano and heat maps were drawn. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analyses were then performed to screen the DEMs-associated biological process and signaling pathways, respectively. Further identification with a comprehensive literature research involving particle exposure, fibrosis, inflammation and lung cancer were used to further screen DEMs of CWP and silicosis. Microarray data showed that 375 and 88 miRNAs were differentially expressed in CWP and silicosis lung tissues compared with healthy lung tissues, while 34 miRNAs were differentially expressed in CWP compared with silicosis lung tissues. The GO and KEGG pathway analyses showed that, the target genes were mainly enriched in the TGF-β, MAPK, p53 and other signal pathways. These results provided insight into the miRNA-related underlying mechanisms of CWP and silicosis, and they provided new clues for miRNAs as biomarkers for the diagnosis and differential diagnosis of these two diseases.

Potential correlation of oral flora with pemphigus vulgaris – A case control study

Background/purposeOral flora is related to various immune-related diseases. Herein we explored the characteristics of oral flora in patients with pemphigus vulgaris (PV) and analyzed the correlation between oral flora and PV. Materials and methodsTwenty-two untreated patients with PV and 12 healthy controls (HC) were included in this case-control study. The characteristics of salivary microbiome were assessed by high-throughput sequencing using the 16S rRNA Illumina MiSeq approach, and differences between the PV and HC groups were determined. The Kyoto Encyclopedia of Genes and Genomes (KEGG) database was applied to screen key metabolic pathways and preliminarily explore potential mechanisms underlying PV occurrence and development. ResultsThe abundance of oral flora in the PV group was significantly lower than that in the HC group, and there were characteristic changes. The relative abundance of Prevotella and Agrobacterium in the PV group was significantly higher than that in the HC group (P < 0.05) and that of Neisseria, Lautropia, and Fusobacterium was significantly lower (P < 0.05). There was a linear correlation between Prevotella and serum Dsg3 level in PV. KEGG pathway analyses indicated significant differences in nine metabolic pathways between the PV and HC groups (P < 0.05), namely carbohydrate metabolism, digestive system, neurodegenerative disease, glycan biosynthesis and metabolism, drug resistance: antimicrobial, infectious disease: viral, circulatory system, excretory system, and nervous system. ConclusionThe oral flora of patients with PV presented characteristic changes, and several metabolic pathways were affected, including N-glycan biosynthesis and metabolism. Prevotella spp. appear to require the most attention in PV. We believe that oral flora dysbacteriosis contributes to PV occurrence and development.

Five-hub genes identify potential mechanisms for the progression of asthma to lung cancer.

Previous studies have shown that asthma is a risk factor for lung cancer, while the mechanisms involved remain unclear. We attempted to further explore the association between asthma and non-small cell lung cancer (NSCLC) via bioinformatics analysis. We obtained GSE143303 and GSE18842 from the GEO database. Lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) groups were downloaded from the TCGA database. Based on the results of differentially expressed genes (DEGs) between asthma and NSCLC, we determined common DEGs by constructing a Venn diagram. Enrichment analysis was used to explore the common pathways of asthma and NSCLC. A protein-protein interaction (PPI) network was constructed to screen hub genes. KM survival analysis was performed to screen prognostic genes in the LUAD and LUSC groups. A Cox model was constructed based on hub genes and validated internally and externally. Tumor Immune Estimation Resource (TIMER) was used to evaluate the association of prognostic gene models with the tumor microenvironment (TME) and immune cell infiltration. Nomogram model was constructed by combining prognostic genes and clinical features. 114 common DEGs were obtained based on asthma and NSCLC data, and enrichment analysis showed that significant enrichment pathways mainly focused on inflammatory pathways. Screening of 5 hub genes as a key prognostic gene model for asthma progression to LUAD, and internal and external validation led to consistent conclusions. In addition, the risk score of the 5 hub genes could be used as a tool to assess the TME and immune cell infiltration. The nomogram model constructed by combining the 5 hub genes with clinical features was accurate for LUAD. Five-hub genes enrich our understanding of the potential mechanisms by which asthma contributes to the increased risk of lung cancer.

Identification of potential biomarkers for colorectal cancer by clinical database analysis and Kaplan-Meier curves analysis.

This study aimed to explore critical genes as potential biomarkers for the diagnosis and prognosis of colorectal cancer (CRC) for clinical utility. To identify and screen candidate genes involved in CRC carcinogenesis and disease progression, we downloaded microarray datasets GSE89076, GSE73360, and GSE32323 from the GEO database identified differentially expressed genes (DEGs), and performed a functional enrichment analysis. A protein-protein interaction network was constructed, and correlated module analysis was performed using STRING and Cytoscape. The Kaplan-Meier survival curve shows the survival of the hub genes. The expression of cyclin-dependent kinase (CDK1), cyclin B1 (CCNB1), and PCNA in tissues and changes in tumor grade were analyzed. A total of 329 DEGs were identified, including 264 upregulated and 65 downregulated genes. The functions and pathways of DEGs include the mitotic cell cycle, poly(A) RNA binding replication, ATP binding, DNA replication, ribosome biogenesis in eukaryotes, and RNA transport. Forty-seven Hub genes were identified, and biological process analysis showed that these genes were mainly enriched in cell cycle and DNA replication. Patients with mutations in CDK1, PCNA, and CCNB1 had poorer survival rates. CDK1, PCNA, and CCNB1 were significantly overexpressed in the tumor tissues. The expression of CDK1 and CCNB1 gradually decreased with increasing tumor grade. CDK1, CCNB1, and PCNA can be used as potential markers for the diagnosis and prognosis of CRC. These genes are overexpressed in colon cancer tissues and are associated with low survival rates in CRC patients.

RNA-Seq-based transcriptomics analysis during the photodynamic therapy of primary cells in secondary hyperparathyroidism.

The aim of this study was to identify changes in gene expression before and after 5-aminolevulinic acid-mediated photodynamic therapy (5-ALA-PDT) and to investigate the potential mechanism of 5-ALA-PDT based on ribonucleic acid sequencing (RNA-Seq) analysis. Secondary hyperparathyroidism (SHPT) primary cells were isolated from surgically excised specimens and exposed to laser light. The transcription profiles of SHPT primary cells were identified through RNA-Seq. Differentially expressed genes (DEGs) were identified. Enrichment of functions and signaling pathway analysis were performed based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Quantitative real-time polymerase chain reaction (RT-qPCR) and western blot analysis were used to validate genes based on RNA-Seq results. In total, 1320 DEGs were identified, of which 1019 genes were upregulated and 301 genes were downregulated. GO and KEGG pathway analyses identified significantly enriched pathways in DEGs, including TGF beta in extracellular matrix (ECM), negative regulation of triglyceride biosynthetic process, protein heterodimerization activity, systemic lupus erythematosus, ECM-receptor interaction, focal adhesion and protein digestion and absorption. Protein-protein interaction (PPI) network analyses identified potential heat shock protein (HSP) interactions among the DEGs. Eight HSP genes were also identified that were most likely involved in 5-ALA-PDT, which were further validated by RT-qPCR and western blotting. The findings of this descriptive study reveal changes in the transcriptome profile during 5-ALA-PDT, suggesting that gene expression and mutation, signaling pathways, and the molecular network are altered in SHPT primary cells. The above findings provide new insight for further studies on the mechanisms underlying 5-ALA-PDT in SHPT.

FNBP4 is a Potential Biomarker Associated with Cuproptosis and Promotes Tumor Progression in Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors that lacks an efficient therapeutic approach because of its elusive molecular mechanisms. This study aimed to investigate the biological function and potential mechanism of formin-binding protein 4 (FNBP4) in HCC. FNBP4 expression in tissues and cells were detected by quantitative real-time PCR (qRT‒PCR), Western blot, and immunohistochemistry (IHC). The Kaplan-Meier method was used to explore the correlation between the FNBP4 expression and clinical survival. MTT, EdU incorporation, colony formation, and Transwell assays were performed to evaluate the function of FNBP4 in cell proliferation and migration in vitro. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was used to explore the potential mechanism of FNBP4. The prognostic risk signature and nomogram were constructed to demonstrate the prognostic value of FNBP4. We found that FNBP4 was upregulated in patients with HCC and associated with poor overall survival (OS). Furthermore, knockdown of FNBP4 inhibited the proliferation and migration in HCC cells. Then, we performed a KEGG pathway analysis of the coexpressed genes associated with FNBP4 and found that FNBP4 may be associated with tumor-related signaling pathways and cuproptosis. We verified that FNBP4 could cause cell cycle progression and inactivation of the hippo signaling pathway. A prognostic risk signature containing three FNBP4-related differentially expressed cuproptosis regulators (DECRs) was established and can be used as an independent risk factor to evaluate the prognosis of patients with HCC. In addition, a nomogram including a risk score and clinicopathological factors was used to predict patient survival probabilities. FNBP4, as a potential biomarker associated with cuproptosis, promotes HCC cell proliferation and metastasis. We provide a new potential strategy for HCC treatment by targeting FNBP4.

Association of gut microbiota and SCFAs with finishing weight of Diannan small ear pigs.

Finishing weight is a key economic trait in the domestic pig industry. Evidence has linked the gut microbiota and SCFAs to health and production performance in pigs. Nevertheless, for Diannan small ear (DSE) pigs, a specific pig breed in China, the potential effect of gut microbiota and SCFAs on their finishing weight remains unclear. Herein, based on the data of the 16S ribosomal RNA gene and metagenomic sequencing analysis, we found that 13 OTUs could be potential biomarkers and 19 microbial species were associated with finishing weight. Among these, carbohydrate-decomposing bacteria of the families Streptococcaceae, Lactobacillaceae, and Prevotellaceae were positively related to finishing weight, whereas the microbial taxa associated with intestinal inflammation and damage exhibited opposite effects. In addition, interactions of these microbial species were found to be linked with finishing weight for the first time. Gut microbial functional annotation analysis indicated that CAZymes, such as glucosidase and glucanase could significantly affect finishing weight, given their roles in increasing nutrient absorption efficiency. Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthologies (KOs) and KEGG pathways analysis indicated that glycolysis/gluconeogenesis, phosphotransferase system (PTS), secondary bile acid biosynthesis, ABC transporters, sulfur metabolism, and one carbon pool by folate could act as key factors in regulating finishing weight. Additionally, SCFA levels, especially acetate and butyrate, had pivotal impacts on finishing weight. Finishing weight-associated species Prevotella sp. RS2, Ruminococcus sp. AF31-14BH and Lactobacillus pontis showed positive associations with butyrate concentration, and Paraprevotella xylaniphila and Bacteroides sp. OF04-15BH were positively related to acetate level. Taken together, our study provides essential knowledge for manipulating gut microbiomes to improve finishing weight. The underlying mechanisms of how gut microbiome and SCFAs modulate pigs' finishing weight required further elucidation.

Identification of TRPC6 as a Novel Diagnostic Biomarker of PM-Induced Chronic Obstructive Pulmonary Disease Using Machine Learning Models.

Chronic obstructive pulmonary disease (COPD) was the third most prevalent cause of mortality worldwide in 2010; it results from a progressive and fatal deterioration of lung function because of cigarette smoking and particulate matter (PM). Therefore, it is important to identify molecular biomarkers that can diagnose the COPD phenotype to plan therapeutic efficacy. To identify potential novel biomarkers of COPD, we first obtained COPD and the normal lung tissue gene expression dataset GSE151052 from the NCBI Gene Expression Omnibus (GEO). A total of 250 differentially expressed genes (DEGs) were investigated and analyzed using GEO2R, gene ontology (GO) functional annotation, and Kyoto Encyclopedia of Genes and Genomes (KEGG) identification. The GEO2R analysis revealed that TRPC6 was the sixth most highly expressed gene in patients with COPD. The GO analysis indicated that the upregulated DEGs were mainly concentrated in the plasma membrane, transcription, and DNA binding. The KEGG pathway analysis indicated that the upregulated DEGs were mainly involved in pathways related to cancer and axon guidance. TRPC6, one of the most abundant genes among the top 10 differentially expressed total RNAs (fold change ≥ 1.5) between the COPD and normal groups, was selected as a novel COPD biomarker based on the results of the GEO dataset and analysis using machine learning models. The upregulation of TRPC6 was verified in PM-stimulated RAW264.7 cells, which mimicked COPD conditions, compared to untreated RAW264.7 cells by a quantitative reverse transcription polymerase chain reaction. In conclusion, our study suggests that TRPC6 can be regarded as a potential novel biomarker for COPD pathogenesis.

Identification and analysis of key microRNAs derived from osteoarthritis synovial fluid exosomes.

Identification and analysis of key microRNAs derived from osteoarthritis synovial fluid exosomes.

Influence of reactive oxygen species responsive self-assembled nanomicelle loaded with pyroptosis inhibitor on full-thickness skin defects in diabetic rats

Objective: To investigate the influence of reactive oxygen species (ROS) responsive self-assembled nanomicelle loaded with pyroptosis inhibitor on full-thickness skin defects in diabetic rats. Methods: Experimental research methods were employed. A nucleotide-binding oligomerization domain (NOD) 1/2 inhibitor (NOD-IN-1) was encapsulated with nanomicelle polyethylene glycol-block-polypropylene sulfide (PEG-b-PPS), and the resulting product was called PEPS@NOD-IN-1. The morphology and hydration particle size of PEG-b-PPS and PEPS@NOD-IN-1 were observed by transmission electron microscope and particle size analyzer, respectively, and the encapsulation rate and drug loading rate of PEPS@NOD-IN-1 to NOD-IN-1 and the cumulative release rate of NOD-IN-1 by PEPS@NOD-IN-1 in phosphate buffer solution (PBS) alone and hydrogen peroxide-containing PBS within 40 h were measured and calculated by microplate reader, and the sample number was 3. Twenty-four male Sprague-Dawley rats aged 6-7 weeks were injected with streptozotocin to induce type 1 diabetes mellitus. Six full-thickness skin defect wounds were made on the back of each rat. The injured rats were divided into PBS group, NOD-IN-1 group, PEG-b-PPS group, and PEPS@NOD-IN-1 group with corresponding treatment according to the random number table, with 6 rats in each group. The wound healing was observed on post injury day (PID) 3, 7, and 12, and the wound healing rate was calculated. The ROS levels in wound tissue were detected by immunofluorescence method on PID 3. On PID 7, the granulation tissue thickness in wound was assessed by hematoxylin-eosin staining, the mRNA expressions of NOD1 and NOD2 were detected by real-time fluorescence quantitative reverse transcription polymerase chain reaction, and the protein expressions of NOD1, NOD2, and GSDMD-N terminals were detected by Western blotting. Six wounds from different rats in each group were taken for detection of the above indicators. Wound tissue (3 samples per group) was taken from rats in PBS group and PEPS@NOD-IN-1 group on PID 7, and transcriptome sequencing was performed using high-throughput sequencing technology platform. Differentially expressed genes (DEGs) significantly down-regulated in PEPS@NOD-IN-1 group as compared with PBS group were screened, and the enrichment analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) was performed. The DEG heatmap of the NOD-like receptor pathway, a pyroptosis-related pathway, was made. Protein-protein interaction (PPI) analysis of DEGs in heatmap was performed through the STRING database to screen key genes of PEPS@NOD-IN-1 regulating the NOD-like receptor pathway. Data were statistically analyzed with analysis of variance for repeated measurement, one-way analysis of variance, and Tukey test. Results: PEG-b-PPS and PEPS@NOD-IN-1 were in spherical structures of uniform size, with hydration particle sizes of (134.2±3.3) and (143.1±2.3) nm, respectively. The encapsulation rate of PEPS@NOD-IN-1 to NOD-IN-1 was (60±5)%, and the drug loading rate was (15±3)%. The release of NOD-IN-1 from PEPS@NOD-IN-1 in PBS alone was slow, and the cumulative release rate at 40 h was only (12.4±2.3)%. The release of NOD-IN-1 from PEPS@NOD-IN-1 in hydrogen peroxide-containing PBS within 10 h was very rapid, and the cumulative release rate at 10 h reached (90.1±3.6)%. On PID 3 and 7, the wounds of rats in the four groups were gradually healed, and the healing in PEPS@NOD-IN-1 group was better than that in the other three groups. On PID 12, the wound scab area in PBS group was large, the wound epithelialization in NOD-IN-1 group and PEG-b-PPS group was obvious, and the wound in PEPS@NOD-IN-1 group was close to complete epithelialization. Compared with those in PBS group, NOD-IN-1 group, and PEG-b-PPS group, the wound healing rates on PID 7 and 12 in PEPS@NOD-IN-1 group were significantly increased (P<0.05), the level of ROS in wound tissue on PID 3 was significantly decreased (P<0.05), the thickness of granulation tissue in wound on PID 7 was significantly thickened (P<0.05), and the mRNA expressions of NOD1 and NOD2 and the protein expressions of NOD1, NOD2, and GSDMD-N terminals in wound tissue on PID 7 were significantly decreased (P<0.05). KEGG pathway analysis showed that DEGs significantly down-regulated in PEPS@NOD-IN-1 group as compared with PBS group were significantly enriched in NOD-like receptors, hypoxia-inducible factors, mitogen-activated protein kinases, and tumor necrosis factor (TNF) pathways. In the DEG heatmap of NOD-like receptor pathway, the genes regulating pyroptosis mainly involved NOD1, NOD2, NOD-like receptor thermoprotein domain-related protein 3, Jun, signal transduction and transcriptional activator 1 (STAT1), TNF-α-induced protein 3. The PPI results showed that NOD1, NOD2, and STAT1 were the key genes of PEPS@NOD-IN-1 regulating the NOD-like receptor pathway. Conclusions: PEPS@NOD-IN-1 can down-regulate the level of local ROS in wounds and the expression of NOD1, NOD2, and GSDMD-N terminals, the key regulators of pyroptosis, thereby promoting the repair of full-thickness skin defect wounds in diabetic rats. PEPS@NOD-IN-1 can also significantly down-regulate the pyroptosis, inflammation, and hypoxia-related pathways of wounds, and regulate NOD-like receptor pathways by down-regulating key genes NOD1, NOD2, and STAT1.

Identification and analysis of lncRNA, miRNA and mRNA related to subcutaneous and intramuscular fat in Laiwu pigs.

Long non-coding RNAs (lncRNAs) regulate adipocyte differentiation and metabolism, However, their function on subcutaneous and intramuscular adipose tissues in pigs is unclear. Intramuscular fat (IMF) is an important indicator for evaluating meat quality. Breeds with high IMF content are often accompanied by high subcutaneous fat (SCF), which severely affects the meat rate of pigs. It is of great significance for porcine breeding to study the mechanism of lncRNA related to adipogenesis and lipid metabolism. We identified differentially expressed lncRNAs, miRNAs and mRNAs in subcutaneous and intramuscular adipose tissues in three female Laiwu pigs by deep RNA-sequencing(|log2foldchange|≥1, P_value ≤ 0.05). The gene expression profiles of IMF and SCF in Laiwu pigs were comparatively analyzed by Bioinformatics methods to identify key lncRNAs, miRNAs, and mRNAs associated with lipid metabolism and adipogenesis. A total of 1209 lncRNAs (DElncRNAs), 286 miRNAs (DEmiRNAs), and 1597 mRNAs (DEgenes) were differentially expressed between two types of adipose. Among them, 17 DElncRNAs and 103 target genes play a role in the co-expression network, as well as 59 DElncRNAs, 44 DEmiRNAs, and 88 DEgenes involved in ceRNA network. In GO(Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis of DElncRNAs their target genes involved in many adipogenesis and lipid metabolism biological processes and signaling pathways, such as PPAR signaling pathway, Wnt signaling pathway, MAPK signaling pathway. By constructing co-expression and ceRNAs network we found that Wnt signaling pathway play a critical regulatory role in intramuscular adipogenesis and lipid accumulation in Laiwu pigs. TCONS_00006525, TCONS_00046551 and TCONS_00000528 may target WNT5A, WNT10B and FDZ3 in co-expression network, TCONS_00026517 and other lncRNAs regulate the expression of PPARG, RXRG and SCD in ceRNA network, and were involved in Wnt signaling pathway. This study provides a theoretical basis for further understanding the post-transcriptional regulation mechanism of meat quality formation, predicting and treating diseases caused by ectopic fat.

Differential expression profiles and functional prediction of circRNA in mice with traumatic heterotopic ossification.

Traumatic heterotopic ossification (HO) is an intractable sequela incited by inflammatory insult. To date, the exact molecular mechanisms of traumatic HO formation remain unclear. Recent studies have indicated that circular RNAs (circRNAs) participate in various human skeletal diseases. Although the formation of HO recapitulates many programs during bone development and remodeling, few data are available concerning whether circRNAs could participate in this pathological osteogenesis. To investigate the differentially expressed circRNAs (DE-circRNAs) in HO formation, microarray assay was performed to analyze the circRNA expression profile in four pairs of mice HO tissues and normal tissues. Then, qRT-PCR was applied to verify the microarray data. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed the biological functions of the differentially expressed circRNAs target genes. Cytoscape software was used to construct the circRNA-miRNA-mRNA network for circRNAs with different expression levels as well as the target genes. We demonstrated that 491 circRNAs were significantly differentially expressed in mouse HO tissues by a fold-change ≥ 2 and p-value ≤ 0.05. Among them, the expressions of 168 circRNAs were increased, while 323 were decreased. The expression levels of 10 selected circRNAs were verified successfully by qRT-PCR. GO analysis exhibited that these DE-circRNAs participated in a series of cellular processes. KEGG pathway analysis revealed that multiple upregulated and downregulated pathways were closely related to the DE-circRNAs in HO mice. The circRNA-miRNA-mRNA networks demonstrated that DE-circRNAs may be involved in the pathological osteogenesis of HO through the circRNA-targeted miRNA-mRNA axis. Our study first demonstrated the expression profiles and predicted the potential functions of DE-circRNAs in mice traumatic HO, which may shed new light on the elucidation of mechanisms as well as provide novel potential peripheral biological diagnostic markers and therapeutic targets for traumatic HO.

A five-pseudouridylation-associated-LncRNA classifier for primary prostate cancer prognosis prediction.

Background: Prostate cancer (PCa) is one of the most common cancers in males around the globe, and about one-third of patients with localized PCa will experience biochemical recurrence (BCR) after radical prostatectomy or radiation therapy. Reportedly, a proportion of patients with BCR had a poor prognosis. Cumulative studies have shown that RNA modifications participate in the cancer-related transcriptome, but the role of pseudouridylation occurring in lncRNAs in PCa remains opaque. Methods: Spearman correlation analysis and univariate Cox regression were utilized to determine pseudouridylation-related lncRNAs with prognostic value in PCa. Prognostic pseudouridylation-related lncRNAs were included in the LASSO (least absolute shrinkage and selection operator) regression algorithm to develop a predictive model. KM (Kaplan-Meier) survival analysis and ROC (receiver operating characteristic) curves were applied to validate the constructed model. A battery of biological cell assays was conducted to confirm the cancer-promoting effects of RP11-468E2.5 in the model. Results: A classifier containing five pseudouridine-related lncRNAs was developed to stratify PCa patients on BCR and named the "ψ-lnc score." KM survival analysis showed patients in the high ψ-lnc score group experienced BCR more than those in the low ψ-lnc score group. ROC curves demonstrated that ψ-lnc score outperformed other clinical indicators in BCR prediction. An external dataset, GSE54460, was utilized to validate the predictive model's efficacy and authenticity. A ceRNA (competitive endogenous RNA) network was constructed to explore the model's potential molecular functions and was annotated through GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analyses. RP11-468E2.5 was picked for further investigation, including pan-cancer analysis and experimental validation. Preliminarily, RP11-468E2.5 was confirmed as a tumor promoter. Conclusion: We provide some evidence that pseudouridylation in lncRNA played a role in the development of PCa and propose a novel prognostic classifier for clinical practice.

Titanium Dioxide Nanoparticles Altered the lncRNA Expression Profile in Human Lung Cells.

Respiration is considered to be the main occupational or environmental exposure pathway of titanium dioxide nanoparticles (TiO2 NPs), and the lung is considered to be the target organ of respiratory exposure; however, the mechanism of respiratory toxicity is not fully understood. In this study, the effect of TiO2 NPs on the expression profile of long non-coding RNA (lncRNA) in bronchial epithelial cells (BEAS-2B) was investigated to understand their potential toxic mechanism. BEAS-2B cells were treated with 100 μg/mL TiO2 NPs for 48 h, then RNA sequencing was performed to screen the differential lncRNAs compared with the control group, and the enrichment pathways of the differentially expressed lncRNAs were further analyzed using the Kyoto Encyclopedia of Genes and Genomes (KEGG). The results identified a total of 45,769 lncRNAs, and 277 different lncRNAs were screened. KEGG pathway analysis showed that the targeted mRNAs of these different lncRNAs were enriched in the pyrimidine metabolism pathway. This work demonstrates that TiO2 NPs could alter the lncRNA expression profile in BEAS-2B cells, and epigenetics may play a role in the mechanism of respiratory toxicity induced by TiO2 NPs.

System Analysis Based on Pancreatic Cancer Progression Identifies BRINP2 as a Novel Prognostic Biomarker.

Pancreatic adenocarcinoma (PAAD) is a malignant tumor of the digestive system, which develops rapidly and has no obvious early symptoms. This study aims to discover the biomarkers associated with PAAD development. We obtained RNA expression of PAAD patient samples and corresponding clinical data from The cancer genome atlas (TCGA), and screened out BMP/RA-inducible neural-specific protein 2 (BRINP2) gene which is highly associated with PAAD severity. Then, gene ontology (GO) enrichment, Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis and single-sample gene set enrichment analysis (ssGSEA) analysis were performed to explore the biological functions of BRINP2. Subsequently, long non-coding RNA (lncRNAs) associated with BRINP2 were screened out via correlation analysis, and Cox regression analysis and least absolute shrinkage selection operator (LASSO) regression analysis were used to construct the risk prediction model. We further validated the expression level of BRINP2 and its associated lncRNAs in BRINP2-associated lncRNAs prognostic model in vitro. We proposed that BRINP2 might be correlated to the tumor immune microenvironment and could also be used as a biomarker for PAAD progression. GO enrichment analysis and KEGG pathway analysis showed that the prognostic model was highly correlated to immune microenvironment-related pathways. Additionally, we established a BRINP2-associated lncRNAs prognostic model consisting of three lncRNAs. We validated the expression trends of BRINP2 and its associated lncRNAs in BRINP2-associated lncRNAs prognostic model in PAAD cells with various severity of metastatic potential using the quantitative real-time PCR (qRT-PCR). Meanwhile, pRRophetic R package was employed to predict potential therapeutic drugs for BRINP2-associated lncRNAs prognostic model of PAAD. The results suggest that BRINP2 can be used as a novel prognostic biomarker for PAAD.

Mechanism of Xiezhuo Huayu Yiqi Tongluo Formula in the Treatment of Uric Acid Nephropathy Based on Network Pharmacology, Molecular Docking, and In Vivo Experiments.

Xiezhuo Huayu Yiqi Tongluo Formula (XHYTF) consists of 14 Chinese herbal medicines. In this study, we investigated the potential mechanism of XHYTF in the treatment of uric acid nephropathy (UAN) through network pharmacology, molecular docking, and in vivo methods. Using various pharmacological databases and analysis platforms, information on the active ingredients and targets of Chinese herbal medicine was collected, and UAN disease targets were retrieved using OMIM, Gene Cards, and NCBI. Then common target proteins were integrated. A Drug-Component-Target (D-C-T) map was constructed to screen core compounds and build a protein-protein interaction (PPI) network. Further, Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed for common targets, and a Drug-Component-Target-Pathway (D-C-T-P) network diagram was constructed. The molecular docking simulation was performed to verify the binding affinity between core components and hub targets. Subsequently, the UAN rat model was established, followed by the collection of serum and renal tissues. The expression levels of indicators in the serum were determined using an enzyme-linked immunosorbent assay. The pathological changes of renal tissues were detected using H & E staining and Masson staining. The expression of related proteins in renal tissue was detected by western blot. In the study, 216 active ingredients and 439 targets in XHYTF were screened, and 868 targets were identified as being related to UAN. Among them, 115 were common targets. Based on the D-C-T network, quercetin, luteolin, β-sitosterol, and stigmasterol were observed to be the key active ingredients of XHYTF that were effective against UAN. The analysis of the PPI network revealed TNF, IL6, AKT1, PPARG, and IL1β as the 5 key targets. GO enrichment analysis revealed that the pathways were mainly concentrated in cell killing, regulation of signaling receptor activity, and other activities. Subsequently, KEGG pathway analysis revealed that multiple signaling pathways, including the HIF-1, PI3K-Akt, IL-17, and other signaling pathways, were closely related to the action of XHYTF. All 5 key targets were confirmed to interact with all core active ingredients. In vivo experiments indicated that XHYTF significantly reduced blood uric acid and creatinine levels, alleviated inflammatory cell infiltration in kidney tissues, reduced the levels of serum inflammatory factors such as TNF-α and IL1β, and ameliorated renal fibrosis in rats with UAN. Finally, western blot revealed decreased levels of PI3K and AKT1 proteins in the kidney, which confirmed the hypothesis. Collectively, our observations demonstrated that XHYTF significantly protects kidney function, including alleviation of inflammation and renal fibrosis via multiple pathways. This study provided novel insights into the treatment of UAN using traditional Chinese medicines.