Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 25 February 2022
Sec. Computational Genomics
This article is part of the Research Topic Advanced computational strategies to tackle human complex diseases View all 28 articles

Drug Discovery in Spinal Cord Injury With Ankylosing Spondylitis Identified by Text Mining and Biomedical Databases

  • Department of Orthopaedics, Shanghai Changzheng Hospital, Naval Medical University, Shanghai, China

Spinal cord injury (SCI) and ankylosing spondylitis (AS) are common inflammatory diseases in spine surgery. However, it is a project where the relationship between the two diseases is ambiguous and the efficiency of drug discovery is limited. Therefore, the study aimed to investigate new drug therapies for SCI and AS. First, text mining was used to obtain the interacting genes related to SCI and AS, and then, the functional analysis was conducted. Protein–protein interaction (PPI) networks were constructed by STRING online and Cytoscape software to identify hub genes. Last, hub genes and potential drugs were performed after undergoing drug–gene interaction analysis, and MicroRNA and transcription factors regulatory networks were also analyzed. Two hundred five genes common to “SCI” and “AS” identified by text mining were enriched in inflammatory responses. PPI network analysis showed that 30 genes constructed two significant modules. Ultimately, nine (SST, VWF, IL1B, IL6, CXCR4, VEGFA, SERPINE1, FN1, and PROS1) out of 30 genes could be targetable by a total of 13 drugs. In conclusion, the novel core genes contribute to a novel insight for latent functional mechanisms and present potential prognostic indicators and therapeutic targets in SCI and AS.

Introduction

Spinal cord injury (SCI) is a serious complication of traumatic diseases such as spinal fracture or dislocation, burdening families, economics, and society (GBD 2017 US Neurological Disorders Collaborators et al., 2021). On the other hand, ankylosing spondylitis (AS) is a chronic disease that mainly affects the spine and accumulates the sacroiliac and surrounding joints (Braun and Sieper, 2007). However, the relationship between the two diseases is unclear and the efficiency of drug discovery is limited. In addition, if patients with AS happen to have SCI, then the question whether AS aggregates the injury degree of the spinal cord needs to be answered. Therefore, mining the underlying pathomechanism is a necessary way to further understand the relationship between the two diseases and discover potential therapies for destroying their interaction.

There is possibly an inevitable connection between the SCI and AS. The incidence of SCI appears to be higher in patients with AS than in the general population, which, of course, is directly related to the increasing incidence of vertebral column fracture in patients with AS (Jacobs and Fehlings, 2008). Under the pathologic mechanism of the disease, combined with an unstable spine, the AS population is predisposed to highly distracting injuries and spinal epidural hematoma formation further increases the severity of SCI. Hence, they likely share similar pathogenesis, genes, and antigens.

In recent years, with the development of bioinformatics methods and computational technology, the genetic mechanism behind diseases can be interpreted to investigate deleterious variants associated with target drugs, such as cancers, degenerative diseases, and COVID-19 (Cheng et al., 2016; Woolston et al., 2019; DiNardo et al., 2020; Li et al., 2021). Having been applied in exploring biomarkers, searching for drug targets, and studying molecular regulation, bioinformatics analysis acts as a catalyst to clarify the potential pathogenic agents (Li et al., 2019a; Tolios et al., 2020). Text mining of biomedical literature has been understood as a valid way to reveal the relationship between genes and pathologies while combining with other bioinformatical tools is an effective protocol to verify the existing therapy.

Despite extensive research on SCI and AS, there still is a gap in understanding intersecting genes and latent targets for the treatment of SCI and AS. This study aimed to mine the correlations between the SCI and AS, further obtaining new drugs to control the progression of the two, including text mining, genetic functional analysis, protein–protein interaction (PPI) analysis, drug–gene interaction, and related regulatory molecular analysis. First, common genes were obtained from the SCI and AS preliminary gene lists. Second, interacting genes conducted functional analysis for clarifying related functions and pathways. Third, PPI analysis is to further screen out core genes and clusters, further purposing drug selection. Last, the hub genes obtained from drug–gene interaction were undergone MicroRNA analysis and transcriptional factors (TFs) regulatory network to explore the potential mechanisms (Figure 1).

FIGURE 1
www.frontiersin.org

FIGURE 1. Summary of overall data mining result. 1) Text mining: 860 genes were found by using the searching term “spinal cord injury” and 423 genes were found via the term “ankylosing spondylitis” in pubmed2ensemble, ultimately interacting 205 genes. 2) Gene set enrichment: DAVID functional enrichment analysis was performed using biological process, cellular component, molecular function, and signaling pathways analysis. Next, 30 genes were screened out by using the STRING and Cytoscape software. 3) Drug–gene interaction and functional analysis; 30 genes were imported into the DGIdb and 13 drugs were regarded as the potential medical therapy, whereas nine genes were selected as the final genes that completed the functional analysis. 4) Nine hub genes were participating in the analysis of their MicroRNA and TFs regulatory network.

Methods and Materials

Text Mining

Pubmed2ensemble, an available website, could be utilized to conduct text mining (http://pubmed2ensembl.ls.manchester.ac.uk/). Relevant kinds of literature and databases are mined in pubmed2ensemble to obtain significant genes. In this study, two queries were performed with “spinal cord injury” and “ankylosing spondylitis”. Two gene lists were generated, and each gene hits were extracted to acquire the common genes for the next analysis.

Gene Ontology and Signal Pathways Analysis

The DAVID (https://david.ncifcrf.gov/) was used to conduct gene ontology (GO) and signal pathways annotations for the genes corresponding to the SCI and AS intersection. Biological process, cellular component, and molecular function genes participated in explained functional information while performing pathways analysis with the Kyoto Encyclopedia of Genes and Genomes (KEGG). False discovery rate (FDR) < 0.05 was set as the cutoff.

Protein–Protein Interaction Analysis and Gene Module Analysis

First, the STRING database (https://www.string-db.org/), a useful tool to complete functional protein association networks, was selected to analyze the PPI. Second, we uploaded the genes to the site, and the maximum interaction score >0.9 (high confidence) was regarded as an important screening indicator. Third, the TSV file downloaded from the STRING was imported into the Cytoscape software for the Molecular Complex Detection (MCODE), finally building up the gene modules (cluster) to own the hub genes, which were applied to the subsequent steps of drug–gene interactions.

Drug–Gene Interaction and Functional Analysis

It was realizable for us to explore the interaction between drugs and genes, which would help us to search for the potential therapy target of SCI and AS. The drug–gene interaction database (https://www.dgidb.org/) is the same as the abovementioned website, open to everyone eager to find effective medicines to treat SCI and AS. The module genes would match with existing drugs whose criteria were obvious interaction type and approved by countries. Besides, final genes after interacting were continued to take part in the process of functional analysis, and FDR of consequences <0.05 was recognized as statistically significant.

MicroRNA Analysis and TFs Regulatory Network of Hub Genes

Targetscan (http://www.targetscan.org/) could be utilized to search for the microRNA of hub genes. In addition, TFs of core genes were relying on another online database named NetworkAnalyst (https://www.networkanalyst.ca/), a comprehensive network visual analytics platform for gene expression analysis.

Result

Identification of Common Genes

After eliminating the duplication of gene symbols, 605 genes were related to SCI, 218 genes were related to AS, and 205 common genes were found from text mining with conceptions of “spinal cord injury” and “ankylosing spondylitis”.

Gene Ontology and Signal Pathways Analysis

Common 205 genes were participating in the process of GO and KEGG analysis. Consequently, the number of GO terms for the biological process is 1,750, that of the cellular component is 37, and that of molecular function is 31. The five most enriched annotations of the biological process were 1) “response to organic substance” (FDR = 2.85E-35), 2) “response to external stimulus” (FDR = 7.78E-32), 3) “positive regulation of multicellular organismal process” (FDR = 1.62E-31), 4) “inflammatory response” (FDR = 1.81E-31), and 5) “positive regulation of response to stimulus” (FDR = 1.70E-30), respectively, containing 115, 96, 51, 58, and 93 genes from the query set. As for cellular component, the five most enriched annotation were 1) “extracellular space” (FDR = 1.71E-32), 2) “extracellular region part” (FDR = 2.19E-25), 3) “extracellular region” (FDR = 1.07E-24), 4) “cell surface” (FDR = 1.25E-22), and 5) “external side of plasma membrane” (FDR = 6.27E-19), containing 85, 123, 133, 54, and 31 genes. There are also five most enriched annotation of molecular function, 1) receptor binding (FDR = 4.41E-28), 2) cytokine receptor binding (FDR = 1.27E-18), 3) cytokine activity (FDR = 5.97E-17), 4) growth factor activity (FDR = 1.98E-12), and 5) identical protein binding (FDR = 3.12E-8).

In addition, the most enriched annotations of signaling pathways were that 1) Cytokine-cytokine receptor interaction (FDR = 3.17E-20), 2) HIF-1 signaling pathway (FDR = 2.22E-12), 3) NF-kappa B signaling pathway (FDR = 2.46E-11), 4) TNF signaling pathway (FDR = 5.86E-11), 5) Jak-STAT signaling pathway (FDR = 7.66E-8), respectively containing 39, 16, 15, 15, and 17 genes from the query set (Figure 2).

FIGURE 2
www.frontiersin.org

FIGURE 2. Gene ontology and signal pathways analysis. Gene ontology analysis classified common genes into biological processes, cellular components, and molecular functions. Green bar charts represent the biological process, blue bar charts represent the cellular component, red bar charts represent the molecular function, purple bar charts represent the signaling pathways, and orange line chart represents −log10 (FDR).

Module Screening of Protein–Protein Interaction Network

Meeting the limitation of high confidence and hiding non-interacting genes, candidate genes were selected to construct a PPI network via multiple protein functions of the STRING database. A “tsv” file was exported from the site and then imported into the Cytoscape software for the clustering analysis, illustrating that there are 142 nodes and 507 edges (Figure 3A). MCODE app, a necessary tool, was used to complete cluster analysis to further screening for the target genes, eventually acquiring various clusters. Cluster 1 constructed by 12 nodes and 86 edges and cluster 2 built up by 18 nodes and 58 edges were chosen to collect important genes (Figures 3B,C).

FIGURE 3
www.frontiersin.org

FIGURE 3. Protein–protein interaction analysis and gene module analysis. (A) Based on the STRING online database, 142 genes and 507 edges formed the network which was under the maximum interaction score >0.9 (high confidence). (B) Cluster 1: the first significant module from the PPI network, containing 12 nodes and 86 edges. (C) Cluster 2: the second significant module from the PPI network, containing 18 nodes and 54 edges.

Potential Therapeutics and Functional Analysis

Drug–gene interaction is the last step to search for the potential signatures and target medicines, investigating existing drugs to alleviate the damaging degree of SCI and AS. Among the PPI and MCODE analysis, a hub gene list that had been eliminated the duplicating part was a channel to the DGIdb (http://www.dgidb.org/), which was consist of 27 various clinical databases and literature, such as DrugBank and NCBI PubMed. With the limitations of presence or absence of approvals and definitive drug–gene interacting type, a sum of 13 drugs was shown in Table 1 and Figure 4, and the interacting score of Cysteamine was highest. Depending on the target drugs analysis result, nine genes (SST, VWF, IL1B, IL6, CXCR4, VEGFA, SERPINE1, FN1, and PROS1), which acted as the core to take part in the occurrence and development of SCI and AS underwent functional analysis (Table 2), which clarified HIF-1 signaling pathway and PI3K-Akt signaling pathway played an important role in the progression of the two diseases.

TABLE 1
www.frontiersin.org

TABLE 1. Potential drugs targeting genes with SCI and AS association.

FIGURE 4
www.frontiersin.org

FIGURE 4. The interactions of hub genes with drugs and MicroRNA. The interaction of nine hub genes with 13 drugs and MicroRNA.

TABLE 2
www.frontiersin.org

TABLE 2. Summary of gene set enrichment analysis.

MicroRNA and TF Regulatory Network

Nine hub genes were corresponding to one or two conserved MircoRNA shown in Figure 4. The TFs regulatory network consisted of nine genes and 122 TFs (Figure 5). CXCR4 was regulated by 47 TFs, FN1 was regulated by 46 TFs, VEGFA was found to be regulated by 39 TFs, and SERPINE1 was regulated by 36 TFs. The numbers of VWF, IL6, PROS1, SST, and IL1B were 29, 25, 25, 18, and 15, respectively. In addition, 11 TFs were found with an interacting degree >4 in the network, implying they may play a key role in the progress of diseases (Table 3).

FIGURE 5
www.frontiersin.org

FIGURE 5. The interactions of hub genes with TFs. The TFs regulatory network of nine hub genes. Red nodes represented nine hub genes and blue squares represented TFs.

TABLE 3
www.frontiersin.org

TABLE 3. Summary of potential TFs of hub genes.

Discussion

Leading to lifelong disability, SCI, a devastating disease, is damage to the spinal cord. As for AS, the autoimmune disease is characterized by inflammation of the sacroiliac joints and the attachment points of the spine. According to the study outcomes, patients with AS are more likely to fracture and develop SCI. However, there has been little or no research into drug discovery for SCI and AS diseases treatment. Therefore, we designed this study to identify the potential drugs of therapeutic value, further providing references for clinical diagnosis and treatment. These 13 drugs target the following nine genes: VEGFA (four drugs), IL1B (two drugs), SST, VWF, IL6, CXCR4, SERPINE1, FN1, and PROS1 (one drug each).

CXCR4 (C-X-C motif chemokine receptor 4), a member of the GPCR family, and CXCL12 constitute an axis to modulate the nerve system. The expression of the CXCL12/CXCR4 axis was increasing in the patients with SCI and AS (Knerlich-Lukoschus et al., 2011; Tysseling et al., 2011). Besides, a study found that CXCR4 could accelerate the ossification of fibroblasts in patients with AS (He et al., 2019). FN1 (fibronectin 1) is involved in the process of cell adhesion and migration (Cai et al., 2018) and could rely on the TGF-β/PI3K/Akt pathway to promote chondrocyte differentiation and collagen production in fractured bones (Zhang et al., 2021). As for PROS1 (protein s), one of its roles is to inhibit tumor metastasis by mediating inflammation and immunization (Maimon et al., 2021). Downregulated SERPINE1 (serpin family E member 1) could accelerate the progression of inflammation (Yang et al., 2017).

When it comes to the VWF (von Willebrand factor) and VEGFA (vascular endothelial growth factor), some reports have revealed the relationship between the two genes and the two diseases (Widenfalk et al., 2003; Divecha et al., 2005; Nesic et al., 2010; Taylan et al., 2012). Experiencing compressive injury, the levels of VWF and VEGFA expression in the rat spinal cord were decreased. Interestingly, the levels of two gene expressions were dramatically increased after decompression (Cheng et al., 2021), which was consensus with another study (Shen et al., 2020). Furthermore, VWF production in men with AS was more than normal people, and it is likely to be recognized as a biomarker of AS (Taylan et al., 2012). In addition, a study found that early VEGF delivery for SCI treatment attributed to its power of protecting vessels and inhibiting cell apoptosis was advantageous (Widenfalk et al., 2003), and it even serves as an important criterion for assessing SCI (Nesic et al., 2010).

GO analysis and KEGG enrichment analysis of DEGs showed that most DEGs were closely related to inflammatory response, response to external stimulus and response to organic substance. PI3k/Akt and HIF-1 are the major signaling pathways that have been identified to be significant in inflammatory diseases. PI3K/Akt signaling pathway is inhibited by endoplasmic reticulum stress after SCI, which is activated to realize motor recovery and nerve protection (Li et al., 2019b). In addition, owning the function of regulating spinal cord regeneration, PI3K/Akt is regarded as a potential target for SCI treatment (Lu et al., 2018). Concerning the HIF-1 signaling pathway, if combined with VEGF, then it is activated, secondary SCI will be inhibited, and local hypoxic ischemia environment will be improving, meanwhile promoting neurological recovery (Chen et al., 2017). Furthermore, the HIF-1 signaling pathway that can also be explored to prevent SCI involves the apoptosis of neuronal cells (Wang et al., 2016). On the basis of the previous studies, HIF-1 activity relies on the HIF-1α that is particularly controlled by the PI3K/Akt signaling pathway at the post-transcriptional protein level (DeBerardinis et al., 2008; Zhang et al., 2018).

The top five on the drug–gene interacting scoring system were Cysteamine (SST binder), Caplacizumab (VWF inhibitor), Canakinumab (IL1B inhibitor), Siltuximab (IL6 antagonist), and Plerixafor (CXCR4 inhibitor). Owning the highest interacting score in these drugs, Cysteamine deserves to be more explored for its therapeutic value for SCI patients with AS, which is mainly applied to cure Cystinosis (Hollywood et al., 2020), atherosclerosis (Wen et al., 2019), and Parkinson’s disease (Cicchetti et al., 2019). Caplacizumab is only used for acquired thrombotic thrombocytopenic purpura (Scully and Spero, 2019), whereas Siltuximab is mostly utilized for Castleman’s disease (van Rhee et al., 2014). Besides, Canakinumab has entered the clinical stage of lung cancer treatment, perhaps becoming the candidate for SCI or AS (Wong et al., 2020). On the basis of the existing researches, inhibition of CXCR4 could alleviate inflammatory pain induced by ischemia–reperfusion (Li et al., 2016).

More than a quarter of identified drugs target VEGFA, including Aflibercept, Bevacizumab, Pegaptanib sodium, and Ranibizumab. Bevacizumab approved by the US Food and Drug Administration (FDA) has been used in the management of neurosurgery diseases such as radiation-induced myelopathy and gliomas (Psimaras et al., 2016; Levenbaum et al., 2019). Surprisingly, some drugs (Aflibercept, Pegaptanib sodium, and Ranibizumab) for the therapy of age-related macular degeneration may also be a benefit for SCI patients with AS. As mentioned before, VEGFA owns the function of inducing angiogenesis, indicating that it plays a key role in the SCI and AS. Although the above drugs are not currently used for AS and SCI, they are approved by the FDA and may be possibly approved in the therapy of the two diseases in the future.

Interestingly, it is worth noting that STAT3 (signal transducer and activator of transcription 3), as the most connected TF, and MicroRNAs, respectively, corresponding to final genes play a significant role in the development of AS and SCI. The expression of STAT3 is high in patients with AS as well as SCI patients, so AS can further aggregate the degree of SCI and induce more severe inflammatory storms, finally influencing the prognosis of patients (Jo et al., 2020). Meantime, TFs and MicorRNA can also regulate the hub genes to mediate inflammations (Qiu et al., 2021). For example, MicroRNA-149 could suppress hepatic inflammatory response through antagonizing STAT3 signaling pathway (Zhang et al., 2017). Besides, MicroRNA-183 regulates LPS-induced oxidative stress by targeting FN1 (Xie et al., 2020) and MicroRNA-140-3p can ameliorate the progression of osteoarthritis via targeting CXCR4 (Ren et al., 2020). Therefore, MicorRNA and TFs may act as important targets for the treatment of SCI patients with AS, providing direction for later researches.

All consequences suggested that common genes are not only responsible for the SCI and AS but also cover a wide range of digestive system neoplasm and blood system via cytokine-cytokine interaction, PI3K/Akt signaling pathway, and HIF-1 signaling pathway. What calls for special attention is that those alternating genes can aggravate the severity of the SCI and AS, thus building up a vicious circle, and there is also a potential risk of tumor formation. Therefore, it is important to give early interventions for SCI patients with AS through potential drugs therapy. On the one hand, the progression of SCI and AS can be postponed by early treatment, finally breaking the vicious circle. On the other hand, tumors will not be given a chance to develop if intervention is used in the early phase, ultimately maintaining body physical activity to promote recovery.

Regrettably, the limitation of this study is that, with the update of the database we used, our study needs to catch up with the development to repeat the study. In addition, the criteria of screening out core genes are subjective, implicating objective criteria need to be constructed to identify the candidate.

Conclusion

In conclusion, our study obtained 205 genes common in both SCI and AS. GO and KEGG analyses of intersecting genes might reveal a novel prospective relationship between SCI and AS. Consequently, the final nine genes (SST, VWF, IL1B, IL6, CXCR4, VEGFA, SERPINE1, FN1, and PROS1) and 13 drugs may provide therapeutic value for SCI patients with AS.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author Contributions

CW designed the experiments, analyze the data, and wrote the manuscript. HM revised the manuscript. WW designed the experiments and analyze the data. XL devised the concept, designed the research, supervised the research, and wrote the manuscript. All authors read and approved the final manuscript.

Funding

This study was supported by the Research Projects of ShangHai ChangZheng Hospital (No. 0910 and No. 2020YCXYJ-ZD06).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Braun, J., and Sieper, J. (2007). Ankylosing Spondylitis. The Lancet 369 (9570), 1379–1390. doi:10.1016/S0140-6736(07)60635-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, X., Liu, C., Zhang, T. N., Zhu, Y. W., Dong, X., and Xue, P. (2018). Down‐regulation of FN1 Inhibits Colorectal Carcinogenesis by Suppressing Proliferation, Migration, and Invasion. J. Cel. Biochem. 119 (6), 4717–4728. doi:10.1002/jcb.26651

CrossRef Full Text | Google Scholar

Chen, H., Li, J., Liang, S., Lin, B., Peng, Q., Zhao, P., et al. (2017). Effect of Hypoxia-Inducible Factor-1/vascular Endothelial Growth Factor Signaling Pathway on Spinal Cord Injury in Rats. Exp. Ther. Med. 13 (3), 861–866. doi:10.3892/etm.2017.4049

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, F., Zhao, J., and Zhao, Z. (2016). Advances in Computational Approaches for Prioritizing Driver Mutations and Significantly Mutated Genes in Cancer Genomes. Brief Bioinform 17 (4), 642–656. doi:10.1093/bib/bbv068

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, X., Yu, Z., Xu, J., Quan, D., and Long, H. (2021). Pathophysiological Changes and the Role of Notch-1 Activation after Decompression in a Compressive Spinal Cord Injury Rat Model. Front. Neurosci. 15, 579431. doi:10.3389/fnins.2021.579431

PubMed Abstract | CrossRef Full Text | Google Scholar

Cicchetti, F., David, L. S., Siddu, A., and Denis, H. L. (2019). Cysteamine as a Novel Disease-Modifying Compound for Parkinson's Disease: Over a Decade of Research Supporting a Clinical Trial. Neurobiol. Dis. 130, 104530. doi:10.1016/j.nbd.2019.104530

PubMed Abstract | CrossRef Full Text | Google Scholar

DeBerardinis, R. J., Lum, J. J., Hatzivassiliou, G., and Thompson, C. B. (2008). The Biology of Cancer: Metabolic Reprogramming Fuels Cell Growth and Proliferation. Cel Metab. 7 (1), 11–20. doi:10.1016/j.cmet.2007.10.002

CrossRef Full Text | Google Scholar

DiNardo, C. D., Tiong, I. S., Quaglieri, A., MacRaild, S., Loghavi, S., Brown, F. C., et al. (2020). Molecular Patterns of Response and Treatment Failure after Frontline Venetoclax Combinations in Older Patients with AML. Blood 135 (11), 791–803. doi:10.1182/blood.2019003988

PubMed Abstract | CrossRef Full Text | Google Scholar

Divecha, H., Sattar, N., Rumley, A., Cherry, L., Lowe, G. D. O., Sturrock, R., et al. (20051979). Cardiovascular Risk Parameters in Men with Ankylosing Spondylitis in Comparison with Non-inflammatory Control Subjects: Relevance of Systemic Inflammation. Clin. Sci. (London, Engl. 109 (2), 171–176. doi:10.1042/cs20040326

CrossRef Full Text | Google Scholar

GBD 2017 US Neurological Disorders Collaborators Feigin, V. L., Vos, T., Alahdab, F., Amit, A. M. L., Bärnighausen, T. W., Beghi, E., et al. (2021). Burden of Neurological Disorders across the US from 1990-2017: A Global Burden of Disease Study. JAMA Neurol. 78 (2), 165–176. doi:10.1001/jamaneurol.2020.4152

PubMed Abstract | CrossRef Full Text | Google Scholar

He, C., Li, D., Gao, J., Li, J., Liu, Z., and Xu, W. (2019). Inhibition of CXCR4 Inhibits the Proliferation and Osteogenic Potential of Fibroblasts from Ankylosing Spondylitis via the Wnt/β-catenin P-athway. Mol. Med. Rep. 19 (4), 3237–3246. doi:10.3892/mmr.2019.9980

PubMed Abstract | CrossRef Full Text | Google Scholar

Hollywood, J. A., Przepiorski, A., D’Souza, R. F., Sreebhavan, S., Wolvetang, E. J., Harrison, P. T., et al. (2020). Use of Human Induced Pluripotent Stem Cells and Kidney Organoids to Develop a Cysteamine/mTOR Inhibition Combination Therapy for Cystinosis. Jasn 31 (5), 962–982. doi:10.1681/ASN.2019070712

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacobs, W. B., and Fehlings, M. G. (2008). Ankylosing Spondylitis and Spinal Cord Injury: Origin, Incidence, Management, and Avoidance. Foc 24 (1), E12. doi:10.3171/FOC/2008/24/1/E12

PubMed Abstract | CrossRef Full Text | Google Scholar

Jo, S., Won, E. J., Kim, M.-J., Lee, Y. J., Jin, S.-H., Park, P.-R., et al. (2020). STAT3 Phosphorylation Inhibition for Treating Inflammation and New Bone Formation in Ankylosing Spondylitis. Rheumatology (Oxford) 60, 3923. doi:10.1093/rheumatology/keaa846

PubMed Abstract | CrossRef Full Text | Google Scholar

Knerlich-Lukoschus, F., von der Ropp-Brenner, B., Lucius, R., Mehdorn, H. M., and Held-Feindt, J. (2011). Spatiotemporal CCR1, CCL3(MIP-1α), CXCR4, CXCL12(SDF-1α) Expression Patterns in a Rat Spinal Cord Injury Model of Posttraumatic Neuropathic Pain. J. Neurosurg. Spine 14 (5), 583–597. doi:10.3171/2010.12.SPINE10480

PubMed Abstract | CrossRef Full Text | Google Scholar

Levenbaum, E., Ellika, S., and Korones, D. N. (2019). Bevacizumab in Treating the Cystic Components of Pediatric Low‐grade Gliomas: A Report of Four Patients. Pediatr. Blood Cancer 66 (11), e27917. doi:10.1002/pbc.27917

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Zhang, X., Qi, X., Zhu, X., and Cheng, L. (2019a). Icariin Inhibits Endoplasmic Reticulum Stress-Induced Neuronal Apoptosis after Spinal Cord Injury through Modulating the PI3K/AKT Signaling Pathway. Int. J. Biol. Sci. 15 (2), 277–286. doi:10.7150/ijbs.30348

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Wu, Z., Zheng, D., Sun, Y., Wang, S., and Yan, Y. (2019b). Bioinformatics Analysis of the Regulatory lncRNA-miRNA-mRNA N-etwork and D-rug P-rediction in P-atients with H-ypertrophic C-ardiomyopathy. Mol. Med. Rep. 20 (1), 549–558. doi:10.3892/mmr.2019.10289

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, R., Li, Y., Liang, X., Yang, L., Su, M., and Lai, K. P. (2021). Network Pharmacology and Bioinformatics Analyses Identify Intersection Genes of Niacin and COVID-19 as Potential Therapeutic Targets. Brief Bioinform 22 (2), 1279–1290. doi:10.1093/bib/bbaa300

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X.-Q., Zhang, Z.-L., Tan, W.-F., Sun, X.-J., and Ma, H. (2016). Down-Regulation of CXCL12/CXCR4 Expression Alleviates Ischemia-Reperfusion-Induced Inflammatory Pain via Inhibiting Glial TLR4 Activation in the Spinal Cord. PloS one 11 (10), e0163807. doi:10.1371/journal.pone.0163807

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, H., Zhang, L.-H., Yang, L., and Tang, P-F. (2018). The PI3K/Akt/FOXO3a Pathway Regulates Regeneration Following Spinal Cord Injury in Adult Rats through TNF-α and P27kip1 Expression. Int. J. Mol. Med. 41 (5), 2832–2838. doi:10.3892/ijmm.2018.3459

CrossRef Full Text | Google Scholar

Maimon, A., Levi-Yahid, V., Ben-Meir, K., Halpern, A., Talmi, Z., Priya, S., et al. (2021). Myeloid Cell-Derived PROS1 Inhibits Tumor Metastasis by Regulating Inflammatory and Immune Responses via IL-10. J. Clin. Invest. 131 (10), e126089. doi:10.1172/JCI126089

CrossRef Full Text | Google Scholar

Nesic, O., Sundberg, L. M., Herrera, J. J., Mokkapati, V. U. L., Lee, J., and Narayana, P. A. (2010). Vascular Endothelial Growth Factor and Spinal Cord Injury Pain. J. Neurotrauma 27 (10), 1793–1803. doi:10.1089/neu.2010.1351

CrossRef Full Text | Google Scholar

Psimaras, D., Tafani, C., Ducray, F., Leclercq, D., Feuvret, L., Delattre, J. Y., et al. (2016). Bevacizumab in Late-Onset Radiation-Induced Myelopathy. Neurology 86 (5), 454–457. doi:10.1212/WNL.0000000000002345

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, Y., Zhao, Z., Chen, Q., Zhang, B., and Yang, C. (2021). MiR-495 Regulates Cell Proliferation and Apoptosis in H2O2 Stimulated Rat Spinal Cord Neurons through Targeting Signal Transducer and Activator of Transcription 3 (STAT3). Ann. Transl Med. 9 (6), 461. doi:10.21037/atm-21-102

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, T., Wei, P., Song, Q., Ye, Z., Wang, Y., and Huang, L. (2020). MiR-140-3p Ameliorates the Progression of Osteoarthritis via Targeting CXCR4. Biol. Pharm. Bull. 43 (5), 810–816. doi:10.1248/bpb.b19-00959

PubMed Abstract | CrossRef Full Text | Google Scholar

Scully, M., Cataland, S. R., Peyvandi, F., Coppo, P., Knöbl, P., Kremer Hovinga, J. A., et al. (2019). Cataland, Flora Peyvandi, Paul Coppo, Paul Knöbl, Johanna A. Kremer Hovinga, Ara Metjian, Javier de la Rubia, Katerina Pavenski, Filip Callewaert, Debjit Biswas, Hilde De Winter, and Robert K. ZeldinCaplacizumab Treatment for Acquired Thrombotic Thrombocytopenic Purpura. N. Engl. J. Med. 380 (4), 335–346. doi:10.1056/NEJMoa1806311

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, N., Wang, L., Wu, Y., Liu, Y., Pei, H., and Xiang, H. (2020). Adeno-associated Virus Packaged TRPC5 Gene Therapy Alleviated Spinal Cord Ischemic Reperfusion Injury in Rats. Neuroreport 31 (1), 29–36. doi:10.1097/WNR.0000000000001359

PubMed Abstract | CrossRef Full Text | Google Scholar

Taylan, A., Sari, I., Kozaci, D. L., Yildiz, Y., Bilge, S., Coker, I., et al. (2012). Evaluation of Various Endothelial Biomarkers in Ankylosing Spondylitis. Clin. Rheumatol. 31 (1), 23–28. doi:10.1007/s10067-011-1760-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Tolios, A., De Las Rivas, J., Hovig, E., Trouillas, P., Scorilas, A., and Mohr, T. (2020). Computational Approaches in Cancer Multidrug Resistance Research: Identification of Potential Biomarkers, Drug Targets and Drug-Target Interactions. Drug Resist. Updates 48, 100662. doi:10.1016/j.drup.2019.100662

PubMed Abstract | CrossRef Full Text | Google Scholar

Tysseling, V. M., Mithal, D., Sahni, V., Birch, D., Jung, H., Miller, R. J., et al. (2011). SDF1 in the Dorsal Corticospinal Tract Promotes CXCR4+ Cell Migration after Spinal Cord Injury. J. Neuroinflammation 8, 16. doi:10.1186/1742-2094-8-16

CrossRef Full Text | Google Scholar

van Rhee, F., Wong, R. S., Munshi, N., Rossi, J.-F., Ke, X.-Y., Fosså, A., et al. (2014). Siltuximab for Multicentric Castleman's Disease: a Randomised, Double-Blind, Placebo-Controlled Trial. Lancet Oncol. 15 (9), 966–974. doi:10.1016/S1470-2045(14)70319-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Li, J., Wu, D., Bu, X., and Qiao, Y. (2016). Hypoxia Promotes Apoptosis of Neuronal Cells through Hypoxia-Inducible Factor-1α-microRNA-204-B-Cell Lymphoma-2 Pathway. Exp. Biol. Med. (Maywood) 241 (2), 177–183. doi:10.1177/1535370215600548

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, Y., Ahmad, F., Mohri, Z., Weinberg, P. D., and Leake, D. S. (2019). Cysteamine Inhibits Lysosomal Oxidation of Low Density Lipoprotein in Human Macrophages and Reduces Atherosclerosis in Mice. Atherosclerosis 291, 9–18. doi:10.1016/j.atherosclerosis.2019.09.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Widenfalk, J., Lipson, A., Jubran, M., Hofstetter, C., Ebendal, T., Cao, Y., et al. (2003). Vascular Endothelial Growth Factor Improves Functional Outcome and Decreases Secondary Degeneration in Experimental Spinal Cord Contusion Injury. Neuroscience 120 (4), 951–960. doi:10.1016/s0306-4522(03)00399-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, C. C., Baum, J., Silvestro, A., Beste, M. T., Bharani-Dharan, B., Xu, S., et al. (2020). Inhibition of IL1β by Canakinumab May Be Effective against Diverse Molecular Subtypes of Lung Cancer: An Exploratory Analysis of the CANTOS Trial. Cancer Res. 80 (24), 5597–5605. doi:10.1158/0008-5472.CAN-19-3176

PubMed Abstract | CrossRef Full Text | Google Scholar

Woolston, A., Khan, K., Spain, G., Barber, L. J., Griffiths, B., Gonzalez-Exposito, R., et al. (2019). Genomic and Transcriptomic Determinants of Therapy Resistance and Immune Landscape Evolution during Anti-EGFR Treatment in Colorectal Cancer. Cancer Cell 36 (1), 35–50. doi:10.1016/j.ccell.2019.05.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, Q., Feng, L., Zhang, C., and Gan, L. (2020). MicroRNA-183 Regulates Lipopolysaccharide-Induced Oxidative Stress of Hippocampal Neurons by Targeting the Fibronectin 1 Gene. Neurosci. Lett. 725, 134902. doi:10.1016/j.neulet.2020.134902

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C.-H., Li, H.-C., Ku, T.-S., Wu, P.-C., Yeh, Y.-J., Cheng, J.-C., et al. (2017). Hepatitis C Virus Down-Regulates SERPINE1/PAI-1 Expression to Facilitate its Replication. J. Gen. Virol. 98 (9), 2274–2286. doi:10.1099/jgv.0.000901

CrossRef Full Text | Google Scholar

Zhang, H., Chen, X., Xue, P., Ma, X., Li, J., and Zhang, J. (2021). FN1 Promotes Chondrocyte Differentiation and Collagen Production via TGF-β/PI3K/Akt Pathway in Mice with Femoral Fracture. Gene 769, 145253. doi:10.1016/j.gene.2020.145253

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Su, J., Wang, Z., Qi, H., Ge, Z., Li, Z., et al. (2017). MicroRNA-149* Suppresses Hepatic Inflammatory Response through Antagonizing STAT3 Signaling Pathway. Oncotarget 8 (39), 65397–65406. doi:10.18632/oncotarget.18541

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Yao, L., Yang, J., Wang, Z., and Du, G. (2018). PI3K/Akt and HIF-1 S-ignaling P-athway in H-ypoxia-ischemia (Review). Mol. Med. Rep. 18 (4), 3547–3554. doi:10.3892/mmr.2018.9375

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: spinal cord injury, ankylosing spondylitis, text mining, drug discovery, bioinformatic analysis

Citation: Wang C, Ma H, Wu W and Lu X (2022) Drug Discovery in Spinal Cord Injury With Ankylosing Spondylitis Identified by Text Mining and Biomedical Databases. Front. Genet. 13:799970. doi: 10.3389/fgene.2022.799970

Received: 22 October 2021; Accepted: 19 January 2022;
Published: 25 February 2022.

Edited by:

Xiao Chang, Children’s Hospital of Philadelphia, United States

Reviewed by:

Xianhu Zhou, Ningbo University, China
Yichuan Liu, Children’s Hospital of Philadelphia (CHOP), United States

Copyright © 2022 Wang, Ma, Wu and Lu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xuhua Lu, xuhualu@hotmail.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.