Identification of diagnostic genes and vital microRNAs involved in rheumatoid arthritis: based on data mining and experimental verification

Chemokines and chemokine receptors as promising targets in rheumatoid arthritis

Rheumatoid arthritis (RA) is an autoimmune disease that commonly causes inflammation and bone destruction in multiple joints. Inflammatory cytokines, such as IL-6 and TNF-α, play important roles in RA development and pathogenesis. Biological therapies targeting these cytokines have revolutionized RA therapy. However, approximately 50% of the patients are non-responders to these therapies. Therefore, there is an ongoing need to identify new therapeutic targets and therapies for patients with RA. In this review, we focus on the pathogenic roles of chemokines and their G-protein-coupled receptors (GPCRs) in RA. Inflamed tissues in RA, such as the synovium, highly express various chemokines to promote leukocyte migration, tightly controlled by chemokine ligand-receptor interactions. Because the inhibition of these signaling pathways results in inflammatory response regulation, chemokines and their receptors could be promising targets for RA therapy. The blockade of various chemokines and/or their receptors has yielded prospective results in preclinical trials using animal models of inflammatory arthritis. However, some of these strategies have failed in clinical trials. Nonetheless, some blockades showed promising results in early-phase clinical trials, suggesting that chemokine ligand-receptor interactions remain a promising therapeutic target for RA and other autoimmune diseases.

Rheumatoid arthritis fibroblast-like synoviocytes maintain tumor-like biological characteristics through ciRS-7-dependent regulation of miR-7

BACKGROUND

Altered phenotype of Fibroblast-like synoviocyte(FLS) is an important cause of the pathogenesis and progression of rheumatoid arthritis(RA), but the specific mechanism causing this change has not yet been fully explained. The exact mechanism by which the biological properties of FLS change in RA is still unclear. microRNAs (miRNAs) have been shown to affect changes in the biological properties of RA-FLS, but the critical miRNAs remain to be discovered. Thus, we first used miRNA microarray and WGCNA to confirm the RA-FLS miRNA landscape and establish their biological functions via network analyses at the system level, as well as to provide a platform for modulating the overall phenotypic effects of RA-FLS.

METHODS

We enrolled a total of 3 patients with RA and 3 healthy participants, constructed a network analysis of via miRNA microarray and RNA-sequencing. Furthermore, the coexpression analyses of miR-7 and ciRS-7 were verified by siRNA transfection, overexpression and qPCR analyses. Finally, we evaluated the effects of adjusting the expression levels of miR-7 and ciRS-7 on RA-FLS, respectively.

RESULTS

We identified distinct miRNA features in RA-FLS, including miR-7, which was significantly lower expressed. Furthermore, we discovered the negative regulatory relationship between ciRS-7 and miR-7 in RA-FLS. Finally, we overexpressed miR-7 in RA-FLS and discovered that miR-7 inhibited RA-FLS hyperproliferation, migration, invasion, and apoptosis, whereas ciRS-7 overexpression reversed these effects.

CONCLUSIONS

The results indicate that the dysregulation of miR-7 in FLS may be involved in the pathological processes of RA and that ciRS-7 induced the suppression of tumor-like biological characters of RA-FLS via modulation of miR-7. These findings help us understand the essential roles of a regulatory interaction between ciRS-7 and miR-7 mediating disease activity of RA, and will facilitate to develop potential intervention target for RA.

Emerging landscape of molecular interaction networks: Opportunities, challenges and prospects

Network biology finds application in interpreting molecular interaction networks and providing insightful inferences using graph theoretical analysis of biological systems. The integration of computational bio-modelling approaches with different hybrid network-based techniques provides additional information about the behaviour of complex systems. With increasing advances in high-throughput technologies in biological research, attempts have been made to incorporate this information into network structures, which has led to a continuous update of network biology approaches over time. The newly minted centrality measures accommodate the details of omics data and regulatory network structure information. The unification of graph network properties with classical mathematical and computational modelling approaches and technologically advanced approaches like machine-learning- and artificial intelligence-based algorithms leverages the potential application of these techniques. These computational advances prove beneficial and serve various applications such as essential gene prediction, identification of drug–disease interaction and gene prioritization. Hence, in this review, we have provided a comprehensive overview of the emerging landscape of molecular interaction networks using graph theoretical approaches. With the aim to provide information on the wide range of applications of network biology approaches in understanding the interaction and regulation of genes, proteins, enzymes and metabolites at different molecular levels, we have reviewed the methods that utilize network topological properties, emerging hybrid network-based approaches and applications that integrate machine learning techniques to analyse molecular interaction networks. Further, we have discussed the applications of these approaches in biomedical research with a note on future prospects.

Identification and Validation of ATF3 Serving as a Potential Biomarker and Correlating With Pharmacotherapy Response and Immune Infiltration Characteristics in Rheumatoid Arthritis

Background: Although disease-modifying antirheumatic drugs (DMARDs) have significantly improved the prognosis of patients with rheumatoid arthritis (RA), approximately 40% of RA patients have limited response. Therefore, it was essential to explore new biomarkers to improve the therapeutic effects on RA. This study aimed to develop a new biomarker and validate it by an in vitro study.

Methods: The RNA-seq and the clinicopathologic data of RA patients were downloaded from Gene Expression Omnibus (GEO) databases. Differentially expressed genes were screened in the GPL96 and GPL570 databases. Then, weighted gene co-expression network analysis (WGCNA) was used to explore the most correlated gene modules to normal and RA synovium in the GPL96 and GPL570 databases. After that, the differentially expressed genes were intersected with the correlated gene modules to find the potential biomarkers. The CIBERSORT tool was applied to investigate the relationship between activated transcription factor 3 (ATF3) expression and the immune cell infiltration, and Gene Set Enrichment Analysis (GSEA) was used to investigate the related signaling pathways of differentially expressed genes in the high and low ATF3 groups. Furthermore, the relationships between ATF3 expression and clinical parameters were also explored in the GEO database. Finally, the role of ATF3 was verified by in vitro cell experiments.

Results: We intersected the differentially expressed genes and the most correlated gene modules in the GPL570 and GPL96 databases and identified that ATF3 is a significant potential biomarker and correlates with some clinical–pharmacological variables. Immune infiltration analysis showed that activated mast cells had a significant infiltration in the high ATF3 group in the two databases. GSEA showed that metabolism-associated pathways belonged to the high ATF3 groups and that inflammation and immunoregulation pathways were enriched in the low ATF3 group. Finally, we validated that ATF3 could promote the proliferation, migration, and invasion of RA fibroblast-like synoviocyte (FLS) and MH7A. Flow cytometry showed that ATF3 expression could decrease the proportion of apoptotic cells and increase the proportion of S and G2/M phase cells.

Conclusion: We successfully identified and validated that ATF3 could serve as a novel biomarker in RA, which correlated with pharmacotherapy response and immune cell infiltration.

Identification of Hub mRNAs and lncRNAs in Atrial Fibrillation Using Weighted Co-expression Network Analysis With RNA-Seq Data

Atrial fibrillation (AF)/paroxysmal AF (PAF) is the main cause of cardiogenic embolism. In recent years, the progression from paroxysmal AF to persistent AF has attracted more and more attention. However, the molecular mechanism of the progression of AF is unclear. In this study, we performed RNA sequencing for normal samples, paroxysmal AF and persistent AF samples to identify differentially expressed gene (DEG) and explore the roles of these DEGs in AF. Totally, 272 differently expressed mRNAs (DEmRNAs) and 286 differentially expressed lncRNAs (DElncRNAs) were identified in paroxysmal AF compared to normal samples; 324 DEmRNAs and 258 DElncRNAs were found in persistent atrial fibrillation compared with normal samples; and 520 DEmRNAs and 414 DElncRNAs were identified in persistent AF compared to paroxysmal AF samples. Interestingly, among the DEGs, approximately 50% were coding genes and around 50% were non-coding RNAs, suggesting that lncRNAs may also have a crucial role in the progression of AF. Bioinformatics analysis demonstrated that these DEGs were significantly related to regulating multiple AF associated pathways, such as the regulation of vascular endothelial growth factor production and binding to the CXCR chemokine receptor. Furthermore, weighted gene co-expression network analysis (WGCNA) was conducted to identify key modules and hub RNAs and lncRNAs to determine their potential associations with AF. Five hub modules were identified in the progression of AF, including blue, brown, gray, turquoise and yellow modules. Interestingly, blue module and turquoise module were significantly negatively and positively correlated to the progression of AF respectively, indicating that they may have a more important role in the AF. Moreover, the hub protein-protein interaction (PPI) networks and lncRNA-mRNA regulatory network were constructed. Bioinformatics analysis on the hub PPI network in turquoise was involved in regulating immune response related signaling, such as leukocyte chemotaxis, macrophage activation, and positive regulation of α-β T cell activation. Our findings could clarify the underlying molecular changes associated fibrillation, and provide a useful resource for identifying AF marker.

Identification of CDK2-Related Immune Forecast Model and ceRNA in Lung Adenocarcinoma, a Pan-Cancer Analysis

Background

Tumor microenvironment (TME) plays important roles in different cancers. Our study aimed to identify molecules with significant prognostic values and construct a relevant Nomogram, immune model, competing endogenous RNA (ceRNA) in lung adenocarcinoma (LUAD).

Methods

“GEO2R,” “limma” R packages were used to identify all differentially expressed mRNAs from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. Genes with P-value <0.01, LogFC>2 or <-2 were included for further analyses. The function analysis of 250 overlapping mRNAs was shown by DAVID and Metascape software. By UALCAN, Oncomine and R packages, we explored the expression levels, survival analyses of CDK2 in 33 cancers. “Survival,” “survminer,” “rms” R packages were used to construct a Nomogram model of age, gender, stage, T, M, N. Univariate and multivariate Cox regression were used to establish prognosis-related immune forecast model in LUAD. CeRNA network was constructed by various online databases. The Genomics of Drug Sensitivity in Cancer (GDSC) database was used to explore correlations between CDK2 expression and IC50 of anti-tumor drugs.

Results

A total of 250 differentially expressed genes (DEGs) were identified to participate in many cancer-related pathways, such as activation of immune response, cell adhesion, migration, P13K-AKT signaling pathway. The target molecule CDK2 had prognostic value for the survival of patients in LUAD (P = 5.8e-15). Through Oncomine, TIMER, UALCAN, PrognoScan databases, the expression level of CDK2 in LUAD was higher than normal tissues. Pan-cancer analysis revealed that the expression, stage and survival of CDK2 in 33 cancers, which were statistically significant. Through TISIDB database, we selected 13 immunodepressants, 21 immunostimulants associated with CDK2 and explored 48 genes related to these 34 immunomodulators in cBioProtal database (P < 0.05). Gene Set Enrichment Analysis (GSEA) and Metascape indicated that 49 mRNAs were involved in PUJANA ATM PCC NETWORK (ES = 0.557, P = 0, FDR = 0), SIGNAL TRANSDUCTION (ES = –0.459, P = 0, FDR = 0), immune system process, cell proliferation. Forest map and Nomogram model showed the prognosis of patients with LUAD (Log-Rank = 1.399e-08, Concordance Index = 0.7). Cox regression showed that four mRNAs (SIT1, SNAI3, ASB2, and CDK2) were used to construct the forecast model to predict the prognosis of patients (P < 0.05). LUAD patients were divided into two different risk groups (low and high) had a statistical significance (P = 6.223e-04). By “survival ROC” R package, the total risk score of this prognostic model was AUC = 0.729 (SIT1 = 0.484, SNAI3 = 0.485, ASB2 = 0.267, CDK2 = 0.579). CytoHubba selected ceRNA mechanism medicated by potential biomarkers, 6 lncRNAs-7miRNAs-CDK2. The expression of CDK2 was associated with IC50 of 89 antitumor drugs, and we showed the top 20 drugs with P < 0.05.

Conclusion

In conclusion, our study identified CDK2 related immune forecast model, Nomogram model, forest map, ceRNA network, IC50 of anti-tumor drugs, to predict the prognosis and guide targeted therapy for LUAD patients.