Identification of differentially expressed and methylated genes associated with rheumatoid arthritis based on network

The biomarkers associated with epithelial-mesenchymal transition in human keloids

INTRODUCTION

A keloid is a type of benign fibrotic disease with similar features to malignancies, including anti-apoptosis, over-proliferation, and invasion. Epithelial-mesenchymal transition (EMT) is a crucial mechanism that regulates the metastatic behavior of tumors. Thus, identifying EMT biomarkers is paramount in comprehensively understanding keloid pathogenesis.

METHODS

To identify the differentially expressed genes (DEGs) GSE92566 dataset, with 3 normal skin and 4 keloid tissues, was downloaded from GEO databases to identify the differentially expressed genes (DEGs). Further, EMT-related genes were downloaded from dbEMT 2.0 databases and intersected with GSE92566 DEGs to identify EMT-related-DEGs (ERDEGs). Subsequently, the ERDEGs were used for GO, KEGG, gene set enrichment analysis (GSEA), protein-protein interaction (PPI), and miRNAs-mRNAs network analysis. To predict small molecules for EMT inhibition, the ERDEGs were imported to cMAP databases, whereas hub genes were imported to DGidb databases. Finally, we carried out qRT-PCR and in vitro experiments to validate our findings.

RESULTS

A total of 122 ERDEGs were identified, including 59 upregulated and 63 down-regulated genes. Moreover, enrichment analysis revealed that focal adhesion, AMPK signal pathway, Wnt signal pathway, and EMT biological process were significantly enriched. STRING databases and Cytoscape software were used to construct the PPI network and EMT-related hub genes. Further, 3 modules were explored from the PPI network using the Molecular Complex Detection (MCODE) plugin. In the Cytohubba plugin, 10 hub genes were explored, including FN1, EGF, SOX9, CDH2, PROM1, EPCAM, KRT19, ITGB1, CD24, and KRT18. These genes were then enriched for the focal adhesion pathway. We constructed a microRNA (miRNA)-mRNA network, which predicted hsa-miR-155-5p (8 edges), hsa-miR-124-3p (7 edges), hsa-miR-145-5p (5 edges), hsa-miR-20a-5p (5 edges) and hsa-let-7b-5p (4 edges) as the most connected miRNAs regulating EMT. Based on the ERDEGs and 10 hub genes mentioned above, ribavirin demonstrated high drug-targeting relevance. Subsequently, qRT-PCR confirmed that the expression of FN1, ITGB1, CDH2, and EPCAM corroborated with previous findings. qRT-PCR also showed that the expression levels of hsa-miR-124-3p and hsa-miR-145-5p were significantly lower in keloids and hsa-miR-155-5p was upregulated in keloids. Finally, by treating human keloid fibroblasts (HKFs) with ribavirin in vitro, we confirmed that ribavirin could inhibit HKFs proliferation and EMT.

CONCLUSION

In summary, this work provides novel EMT biomarkers in keloids and predicts new small target molecules for keloid therapy. Our findings improve the understanding of keloid pathogenesis, providing new treatment options.

Cell-Specific Gene Networks and Drivers in Rheumatoid Arthritis Synovial Tissues

Rheumatoid arthritis (RA) is a common autoimmune and inflammatory disease characterized by inflammation and hyperplasia of the synovial tissues. RA pathogenesis involves multiple cell types, genes, transcription factors (TFs) and networks. Yet, little is known about the TFs, and key drivers and networks regulating cell function and disease at the synovial tissue level, which is the site of disease. In the present study, we used available RNA-seq databases generated from synovial tissues and developed a novel approach to elucidate cell type-specific regulatory networks on synovial tissue genes in RA. We leverage established computational methodologies to infer sample-specific gene regulatory networks and applied statistical methods to compare network properties across phenotypic groups (RA versus osteoarthritis). We developed computational approaches to rank TFs based on their contribution to the observed phenotypic differences between RA and controls across different cell types. We identified 18,16,19,11 key regulators of fibroblast-like synoviocyte (FLS), T cells, B cells, and monocyte signatures and networks, respectively, in RA synovial tissues. Interestingly, FLS and B cells were driven by multiple independent co-regulatory TF clusters that included MITF, HLX, BACH1 (FLS) and KLF13, FOSB, FOSL1 (synovial B cells). However, monocytes were collectively governed by a single cluster of TF drivers, responsible for the main phenotypic differences between RA and controls, which included RFX5, IRF9, CREB5. Among several cell subset and pathway changes, we also detected reduced presence of NKT cell and eosinophils in RA synovial tissues. Overall, our novel approach identified new and previously unsuspected KDG, TF and networks and should help better understanding individual cell regulation and co-regulatory networks in RA pathogenesis, as well as potentially generate new targets for treatment.

Gene Network Analyses Identify Co-regulated Transcription Factors and BACH1 as a Key Driver in Rheumatoid Arthritis Fibroblast-like Synoviocytes

RNA-sequencing and differential gene expression studies have significantly advanced our understanding of pathogenic pathways underlying Rheumatoid Arthritis (RA). Yet, little is known about cell-specific regulatory networks and their contributions to disease. In this study, we focused on fibroblast-like synoviocytes (FLS), a cell type central to disease pathogenesis and joint damage in RA. We used a strategy that computed sample-specific gene regulatory networks (GRNs) to compare network properties between RA and osteoarthritis FLS. We identified 28 transcription factors (TFs) as key regulators central to the signatures of RA FLS. Six of these TFs are new and have not been previously implicated in RA, and included BACH1, HLX, and TGIF1. Several of these TFs were found to be co-regulated, and BACH1 emerged as the most significant TF and regulator. The main BACH1 targets included those implicated in fatty acid metabolism and ferroptosis. The discovery of BACH1 was validated in experiments with RA FLS. Knockdown of BACH1 in RA FLS significantly affected the gene expression signatures, reduced cell adhesion and mobility, interfered with the formation of thick actin fibers, and prevented the polarized formation of lamellipodia, all required for the RA destructive behavior of FLS. This is the first time that BACH1 is shown to have a central role in the regulation of FLS phenotypes, and gene expression signatures, as well as in ferroptosis and fatty acid metabolism. These new discoveries have the potential to become new targets for treatments aimed at selectively targeting the RA FLS.

LCK, FOXC1 and hsa-miR-146a-5p as potential immune effector molecules associated with rheumatoid arthritis

Rheumatoid arthritis (RA) is a type of systemic immune disease characterized by chronic inflammatory disease of the joints. However, the aetiology and underlying molecular events of RA are unclear. Here, we applied bioinformatics analysis to identify potential immune effector molecules involved in RA. The three microarray datasets were downloaded from the Gene Expression Omnibus (GEO) database. We used the R software screen 115 overlapping differentially expressed genes (DEGs). Subsequently, we constructed a protein-protein interaction (PPI) network encoded by these DEGs and identified 10 genes closely associated with RA - LCK, GZMA, GZMB, CD2, LAG3, IL-15, TNFRSF4, CD247, CCR5 and CCR7. Furthermore, in the miRNA-hub gene networks, we screened out hsa-miR-146a-5p, which is the miRNA controlling the largest number of hub genes. Finally, we found some transcription factors that closely interact with hub genes, such as FOXC1, GATA2, YY1, RUNX2, SREBF1, CEBPB and NFIC. This study successfully predicted that LCK, FOXC1 and hsa-miR-146a-5p can be used as potential immune effector molecules of RA. Our study may have potential implications for future prediction of disease progression in patients with symptomatic RA, and has important significance for the pathogenesis and targeted therapy of RA.

Predictive plasma biomarker for gestational diabetes: A case-control study in China

OBJECTIVE

This study aims to find new plasma biomarkers in early pregnancy.

DESIGN

The original study enrolled 1219 pregnant women. We investigated protein expression profiles of placental tissues from women with GDM (n = 89) and normal glucose tolerance (NGT) (n = 83). Maternal plasma samples between two groups in early and middle pregnancy were used for validation of candidate biomarkers.

METHODS

Differentially expressed proteins (DEPs) were identified by label-free quantitative proteomics from human placenta samples between two groups. Several DEPs were validated in plasma by Luminex assays. An automatic biochemical analyzer was used to detect blood lipid indexes. The associations of GAL-3BP with biochemical indicators were demonstrated by Pearson's correlation analysis. Binary logistic regression was used to model potential predictive indicators in early pregnancy of GDM. Receiver operating characteristic (ROC) curve was used to evaluate the diagnostic accuracy of the predictive model and the value of GAL-3BP.

RESULTS

123 DEPs were found in placenta involved in ribosomal function, pancreatic secretion, oxidative phosphorylation, and inflammatory signaling pathway. Plasma GAL-3BP are significantly higher in women with GDM than NGT in the first (p = 0.008) and second (p = 0.026) trimester, but C9 and VWF have no difference. The predictive value of GAL-3BP in the first trimester of pregnancy (AUC 0.64) is better than that in the second trimester (AUC 0.61), and combined predictive model of TG and GAL-3BP at early pregnancy has greater predictive and diagnostic value for GDM (AUC 0.69) than individual GAL-3BP (AUC 0.64).

CONCLUSIONS

Plasma TG and GAL-3BP has good predictive and diagnostic value at early pregnancy, suggesting that these two indicators may be used as biomarkers for early prediction and diagnosis of GDM.

SIGNIFICANCE

The advantage of this study is that circulating TG and GAL-3BP might differentiate the progress of women with GDM and normal glucose tolerance (NGT) at the early stage of pregnancy. It is the first study to consider the role of GAL-3BP as an early predictive biomarker in the development of GDM during the whole pregnancy. Another advantage is that volunteers in this study were recruited from two provinces in China to eliminate the impacts of environmental confounders. The similar changes of blood glucose/lipid indicators for women with GDM and NGT in both regions was found in the first and second trimester of pregnancy, which added to the reliability of analytical results.

Identifying potential signatures for atherosclerosis in the context of predictive, preventive, and personalized medicine using integrative bioinformatics approaches and machine-learning strategies

Background

Atherosclerosis is a major contributor to morbidity and mortality worldwide. Although several molecular markers associated with atherosclerosis have been developed in recent years, the lack of robust evidence hinders their clinical applications. For these reasons, identification of novel and robust biomarkers will directly contribute to atherosclerosis management in the context of predictive, preventive, and personalized medicine (PPPM). This integrative analysis aimed to identify critical genetic markers of atherosclerosis and further explore the underlying molecular immune mechanism attributing to the altered biomarkers.

Methods

Gene Expression Omnibus (GEO) series datasets were downloaded from GEO. Firstly, differential expression analysis and functional analysis were conducted. Multiple machine-learning strategies were then employed to screen and determine key genetic markers, and receiver operating characteristic (ROC) analysis was used to assess diagnostic value. Subsequently, cell-type identification by estimating relative subsets of RNA transcript (CIBERSORT) and a single-cell RNA sequencing (scRNA-seq) data were performed to explore relationships between signatures and immune cells. Lastly, we validated the biomarkers' expression in human and mice experiments.

Results

A total of 611 overlapping differentially expressed genes (DEGs) included 361 upregulated and 250 downregulated genes. Based on the enrichment analysis, DEGs were mapped in terms related to immune cell involvements, immune activating process, and inflaming signals. After using multiple machine-learning strategies, dehydrogenase/reductase 9 (DHRS9) and protein tyrosine phosphatase receptor type J (PTPRJ) were identified as critical biomarkers and presented their high diagnostic accuracy for atherosclerosis. From CIBERSORT analysis, both DHRS9 and PTPRJ were significantly related to diverse immune cells, such as macrophages and mast cells. Further scRNA-seq analysis indicated DHRS9 was specifically upregulated in macrophages of atherosclerotic lesions, which was confirmed in atherosclerotic patients and mice.

Conclusions

Our findings are the first to report the involvement of DHRS9 in the atherogenesis, and the proatherogenic effect of DHRS9 is mediated by immune mechanism. In addition, we confirm that DHRS9 is localized in macrophages within atherosclerotic plaques. Therefore, upregulated DHRS9 could be a novel potential target for the future predictive diagnostics, targeted prevention, patient stratification, and personalization of medical services in atherosclerosis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-022-00289-y.

Identification of SLAMF1 as an immune-related key gene associated with rheumatoid arthritis and verified in mice collagen-induced arthritis model

Background

Rheumatoid arthritis (RA) is the most common inflammatory arthropathy. Immune dysregulation was implicated in the pathogenesis of RA. Thus, the aim of the research was to determine the immune related biomarkers in RA.

Methods

We downloaded the gene expression data of RA in GSE89408 and GSE45291 from Gene Expression Omnibus public database (GEO). Differentially expressed genes (DEGs) were identified between RA and control groups. Infiltrating immune cells related genes were obtained by ssGSEA and weighted gene co-expression network analysis (WGCNA). We performed functional enrichment analysis of differentially expressed immunity-related genes (DEIRGs) by “clusterProfiler” R package, key genes screening by protein-protein interaction (PPI) network of DEIRGs. And mice collagen-induced arthritis (CIA) model was employed to verify these key genes.

Results

A total of 1,885 up-regulated and 1,899 down-regulated DEGs were identified in RA samples. The ssGSEA analysis showed that the infiltration of 25 cells was significantly different. 603 immune related genes were obtained by WGCNA, and 270 DEIRGs were obtained by taking the intersection of DEGs and immune related genes. Enrichment analyses indicated that DEIRGs were associated with immunity related biological processes. 4 candidate biomarkers (CCR7, KLRK1, TIGIT and SLAMF1) were identified from the PPI network of DEIRGs and literature research.

In mice CIA model, the immunohistochemical stain showed SLAMF1 has a significantly high expression in diseased joints. And flow cytometry analysis shows the expression of SLAMF1 on CIA mice-derived CTL cells, Th, NK cells, NKT cells, classical dendritic cell (cDCs) and monocytes/macrophages was also significantly higher than corresponding immune cells from HC mice.

Conclusion

Our study identified SMLAF1 as a key biomarker in the development and progression of RA, which might provide new insight for exploring the pathogenesis of RA.

Kinetochore-associated protein 1 promotes the invasion and tumorigenicity of cervical cancer cells via matrix metalloproteinase-2 and matrix metalloproteinase-9

Cervical cancer, a common cancer in women, has become a serious social burden. Kinetochore-associated protein 1 (KNTC1) that regulates the cell cycle by regulating mitosis is related to the malignant behavior of different types of tumors. However, its role in the development of cervical cancer remains unclear. In this study, we initially explored the role of KNTC1 in cervical cancer. KNTC1 expression and relevant information were downloaded from The Cancer Genome Atlas (TCGA) and dataset GSE63514 in the Gene Expression Omnibus (GEO) database for bioinformatics analyses. Cell proliferation was detected by cell counting kit-8 (CCK8) and colony formation assays. Wound healing and Transwell assays were used to evaluate cell migration and invasion abilities. Protein expression levels of matrix metallopeptidase 2 (MMP2) and matrix metallopeptidase 9 (MMP9) were measured by western blotting. Nude mouse models of subcutaneous xenograft tumor were constructed to analyze tumor growth in vivo. CCK8 and colony formation assay results demonstrated that the proliferation rate of SiHa and C-33A cells decreased when KNTC1 was silenced. Western blot and Transwell assays indicated that KNTC1 knockdown weakened the invasion and migration abilities of SiHa and C-33A cells and decreased the expression of MMP-2 and MMP-9. In-vivo experiments suggested that the inhibition of KNTC1 reduced tumor growth. Taken together, our study showed that KNTC1 plays an important role in cervical cancer. Further, we verified the promotional effect of KNTC1 on cervical cancer through in-vivo and in-vitro experiments and speculated that KNTC1 might mediate tumor invasion via MMP9 and MMP2.

Dehydrogenase reductase 9 (SDR9C4) and related homologs recognize a broad spectrum of lipid mediator oxylipins as substrates

Bioactive oxylipins play multiple roles during inflammation and in the immune response, with termination of their actions partly dependent on the activity of yet-to-be characterized dehydrogenases. Here, we report that human microsomal dehydrogenase reductase 9 (DHRS9, also known as SDR9C4 of the short-chain dehydrogenase/reductase (SDR) superfamily) exhibits a robust oxidative activity toward oxylipins with hydroxyl groups located at carbons C9 and C13 of octadecanoids, C12 and C15 carbons of eicosanoids, and C14 carbon of docosanoids. DHRS9/SDR9C4 is also active toward lipid inflammatory mediator dihydroxylated Leukotriene B₄ and proresolving mediators such as tri-hydroxylated Resolvin D1 and Lipoxin A₄, although notably, with lack of activity on the 15-hydroxyl of prostaglandins. We also found that the SDR enzymes phylogenetically related to DHRS9, i.e., human SDR9C8 (or retinol dehydrogenase 16), the rat SDR9C family member known as retinol dehydrogenase 7, and the mouse ortholog of human DHRS9 display similar activity toward oxylipin substrates. Mice deficient in DHRS9 protein are viable, fertile, and display no apparent phenotype under normal conditions. However, the oxidative activity of microsomal membranes from the skin, lung, and trachea of Dhrs9^−/− mice toward 1 μM Leukotriene B₄ is 1.7- to 6-fold lower than that of microsomes from wild-type littermates. In addition, the oxidative activity toward 1 μM Resolvin D1 is reduced by about 2.5-fold with DHRS9-null microsomes from the skin and trachea. These results strongly suggest that DHRS9 might play an important role in the metabolism of a wide range of bioactive oxylipins in vivo.

MiR-144-3p induced by SP1 promotes IL-1β-induced pyroptosis in chondrocytes via PTEN/PINK1/Parkin axis

Rheumatoid arthritis (RA) often leads to functional disabilities and deformities. MiRNA plays a vital role in cell pyroptosis. Nevertheless, the function and underlying mechanism of miR-144-3p in pyroptosis during the progression of RA remains unclear. In this study, N1511 cells were stimulated with IL-1β to construct a RA model. 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide (MTT) assay was performed to assess the cell viability. Cell pyroptosis was detected by flow cytometry. The levels of inflammatory cytokines (TNF-α, IL-6, and IL-18) were assessed by enzyme-linked immunosorbent assay (ELISA). The relationship among specific protein 1 (SP1), microRNA-144-3p (miR-144-3p), and phosphatase and tensin homolog (PTEN) was explored by dual-luciferase reporter assay, RNA immunoprecipitation (RIP), and chromatin immunoprecipitation (ChIP), respectively. The level of miR-144-3p in N1511 cells was upregulated by IL-1β. MiR-144-3p knockdown inhibited IL-1β-induced pyroptosis in N1511 cells, and the expressions of NOD-like receptor family pyrin domain containing 3 (NLRP3), Cleaved caspase-1, Gasdermin D (GSDMD), and Cleaved caspase-3 in IL-1β-stimulated N1511 cells were increased. The levels of inflammatory cytokines in N1511 cells were increased by IL-1β, which were restored by miR-144-3p knockdown. MiR-144-3p knockdown abolished IL-1β-induced inactivation of putative kinase 1 (PINK1)/Parkin RBR E3 ubiquitin-protein (Parkin) signalling. Moreover, transcription factor SP1 could upregulate miR-144-3p expression and miR-144-3p negatively regulated PTEN expression. In summary, MiR-144-3p induced by SP1 could promote IL-1β-induced chondrocyte pyroptosis via inhibiting PTEN expression and suppressing the activation of PINK1/Parkin signalling, which provided a new strategy against RA.

Sema4D correlates with tumour immune infiltration and is a prognostic biomarker in bladder cancer, renal clear cell carcinoma, melanoma and thymoma

Sema4D, a member of the immune semaphorin family, plays crucial roles in the immune regulation, bone resorption and nervous system. It is also involved in angiogenesis and tumour progression. However, systemic studies on the correlation between Sema4D expression and the immune infiltration or clinical outcomes in tumours are still limited. Here, we analysed the landscape of Sema4D expression and its prognostic value in the cancer genome atlas pan-cancer as well as the correlation between Sema4D and immune cell infiltration by Tumour Immune Estimation Resource and Gene Expression Profiling interactive analysis online tools. Results showed that a higher Sema4D expression was significantly correlated with a favourable overall survival in diverse solid tumours including bladder cancer (Hazards Ratio (HR)=0.68, p = .0095), kidney renal clear cell carcinoma (HR = 0.61, p = .0016), melanoma (HR = 0.58, p = 6.6e-05) and thymoma (HR = 0.1, p = .011). Interestingly, Sema4D expression has positive correlation with various tumour infiltrating immune cells and immune cell biomarkers in these tumours. These results suggest that Sema4D could be a prospective biomarker for calculating hazard ratio of tumour patients and their tumour immune infiltration levels.

Deer antler extract potentially facilitates xiphoid cartilage growth and regeneration and prevents inflammatory susceptibility by regulating multiple functional genes

BACKGROUND

Deer antler is a zoological exception due to its fantastic characteristics, including amazing growth rate and repeatable regeneration. Deer antler has been used as a key ingredient in traditional Chinese medicine relating to kidney and bone health for centuries. The aim of this study was to dissect the molecular regulation of deer antler extract (DAE) on xiphoid cartilage (XC).

METHODS

The DAE used in this experiment was same as the one that was prepared as previously described. The specific pathogen-free (SPF) grade Sprague-Dawley (SD) rats were randomly divided into blank group (n =10) and DAE group (n =10) after 1-week adaptive feeding. The DAE used in this experiment was same as the one that was prepared as previously described. The rats in DAE group were fed with DAE for 3 weeks at a dose of 0.2 g/kg per day according to the body surface area normalization method, and the rats in blank group were fed with drinking water. Total RNA was extracted from XC located in the most distal edge of the sternum. Illumina RNA sequencing (RNA-seq) in combination with quantitative real-time polymerase chain reaction (qRT-PCR) validation assay was carried out to dissect the molecular regulation of DAE on XC.

RESULTS

We demonstrated that DAE significantly increased the expression levels of DEGs involved in cartilage growth and regeneration, but decreased the expression levels of DEGs involved in inflammation, and mildly increased the expression levels of DEGs involved in chondrogenesis and chondrocyte proliferation.

CONCLUSIONS

Our findings suggest that DAE might serve as a complementary therapeutic regent for cartilage growth and regeneration to treat cartilage degenerative disease, such as osteoarthritis.