Jing Ye1†
Bingbing Xu1†
Baoshi Fan2Jiying Zhang1
Fuzhen Yuan1
Yourong Chen1Zewen Sun1
Xin Yan1Yifan Song1Shitang Song1Meng Yang2
Jia-Kuo Yu1*- 1Knee Surgery Department of the Institution of Sports Medicine, Peking University Third Hospital, Beijing Key Laboratory of Sports Injuries, Beijing, China
- 2School of Clinical Medicine, Weifang Medical University, Weifang, China
Background: Unlike bone tissue, little progress has been made regarding cartilage regeneration, and many challenges remain. Furthermore, the key roles of cartilage lesion caused by traumas, focal lesion, or articular overstress remain unclear. Traumatic injuries to the meniscus as well as its degeneration are important risk factors for long-term joint dysfunction, degenerative joint lesions, and knee osteoarthritis (OA) a chronic joint disease characterized by degeneration of articular cartilage and hyperosteogeny. Nearly 50% of the individuals with meniscus injuries develop OA over time. Due to the limited inherent self-repair capacity of cartilage lesion, the Biomaterial drug-nanomedicine is considered to be a promising alternative. Therefore, it is important to elucidate the gene potential regeneration mechanisms and discover novel precise medication, which are identified through this study to investigate their function and role in pathogenesis.
Methods: We downloaded the mRNA microarray statistics GSE117999, involving paired cartilage lesion tissue samples from 12 OA patients and 12 patients from a control group. First, we analyzed these statistics to recognize the differentially expressed genes (DEGs). We then exposed the gene ontology (GO) annotation and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway enrichment analyses for these DEGs. Protein-protein interaction (PPI) networks were then constructed, from which we attained eight significant genes after a functional interaction analysis. Finally, we identified a potential nanomedicine attained from this assay set, using a wide range of inhibitor information archived in the Search Tool for the Retrieval of Interacting Genes (STRING) database.
Results: Sixty-six DEGs were identified with our standards for meaning (adjusted P-value < 0.01, |log2 - FC| ≥1.2). Furthermore, we identified eight hub genes and one potential nanomedicine - Selenocysteine based on these integrative data.
Conclusion: We identified eight hub genes that could work as prospective biomarkers for the diagnostic and biomaterial drug treatment of cartilage lesion, involving the novel genes CAMP, DEFA3, TOLLIP, HLA-DQA2, SLC38A6, SLC3A1, FAM20A, and ANO8. Meanwhile, these genes were mainly associated with immune response, immune mediator induction, and cell chemotaxis. Significant support is provided for obtaining a series of novel gene targets, and we identify potential mechanisms for cartilage regeneration and final nanomedicine immunotherapy in regenerative medicine.
Introduction
Cartilage lesion occur as a result of destructive joint diseases, such as osteoarthritis (OA) (Hu et al., 2018). They can cause disability, joint pain, movement limitation, and function impairment (Renders et al., 2014). Currently, there are no efficient actions for cartilage regeneration, due to the substandard inherent repair capacity of the damaged part of the cartilage lesion (Alcidi et al., 2007). Articular cartilage has the unique effect of conducting stress and reducing friction, and its damage can lead to joint dysfunction, and even disability (Gao et al., 2019). The articular cartilage itself has no blood supply, nerves, or lymphoid tissues. It also lacks chondrocytes, and has relatively low sculpting capacity (Cheng et al., 2019). The chondrocytes in the matrix pits have low metabolic activity (Lawrence et al., 2017). They mainly obtain necessary nutrients and excrete metabolites through penetration. The inability or difficulty to repair itself has been a major challenge for the orthopaedic community (Pang et al., 2013). Clinical treatments of articular cartilage lesion can be divided into being either reparative or non-reparative (Uehara et al., 2015). There are two types of reparative surgery: biological and non-biological. Non-reparative operations are debridement and joint irrigation. Non-biological methods, such as artificial joint prosthesis replacement and articular surface shaping, have achieved very good results and functional recovery (Yao et al., 2015; Ragni et al., 2019). Additionally, extreme and burning knee cartilage lesion often need total knee arthroplasty (TKA), which is a surgical option for easing the pain and facilitating knee reconstruction (Math et al., 2006).
The range of methods extensively used for cartilage regeneration in clinical practice involves distinctive types of scaffolds that imitate the native atmosphere, following: (1) mosaicplasty—the replacement of the missing cartilage with an autologous transplant-collected cartilage, such as the double layer collagen type I/III scaffold (MACI) (Liu et al., 2018); (2) microfracture—the disorder of the collagen scaffold, where double layer type I collagen sponge hols chondroitin sulfate to support recruitment of bone marrow stromal cells to the cartilage defect site (Wang et al., 2019); (3) ACI—the vitro culture of autologous chondrocytes that are collected from a special area of the cartilage and then injected into the defect site, covering them with a recyclable collagen membrane (Zhang et al., 2019); and (4) MACI—the transplantation of a viable scaffold surrounding the previous autologous chondrocytes culture (Li et al., 2016). However, currently, the effectiveness of these methods regarding cartilage regeneration remains far from satisfactory (Chen et al., 2020). Consequently, there is a vital need to identify potential mechanisms for cartilage regeneration, and to mine efficacious nanomedicines for regenerative medicine.
Bioinformatics is an emerging interdisciplinary field that manages the storing, repossession, sharing and best use of data and skills, for problem-solving and decision-making purposes (Hasman et al., 2011). The development and renewal of bioinformatics has provided us with the opportunity to mine large databases and uncover more meaningful solutions (Cheng, 2018). Massive databases have increased from cartilage lesion samples during these years, and a great deal of the differentially expressed genes (DEGs) have been determined using the gene ontology (GO) annotation and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway enrichment analyses (Huang et al., 2018). At present, some basic bioinformatics tools have been applied in the clinic, providing powerful weapons for clinical diagnosis, prevention, treatment, and clinical efficacy evaluation. For example, the GeneBank database and OMIM database are widely used by clinicians to search for pathogen or human disease-related gene information. They can then design specific primers and probes through biological software for the diagnosis, typing, quantification, and the identification of drug resistance genes of clinical genes. This plays an important role in the prevention of infectious diseases, genetic diseases, tumors, early diagnosis, treatment, and prognosis.
In our analysis, we took the GSE117999 mRNA expression statistics with the Gene Expression Omnibus (GEO). DEG analyses were made from the cartilage tissue from 12 osteoarthritis patients and 12 patients without osteoarthritis (marked control), using the R software (Jalal et al., 2017) by means of the R-limma et al. packages (Ritchie et al., 2015). Subsequently, GO and KEGG analyses were accomplished using the website of DAVID (Database for Annotation, Visualization and Integrated Discovery). The Protein-protein interaction networks (PPI) networks were constructed, and eight HUB genes were identified. Finally, we discovered a cathelicidin antimicrobial peptide (CAMP) inhibitor, Selenocysteine, using the DGIdb database (Cotto et al., 2018). Using this approach (Figure 1/Figure 2), we recognized numerous potentially important cartilage regeneration - associated genes and pathways.
Figure 1 General text mining approach. The data mining was expanded to distinguish genes connected with osteoarthritis (OA) and patients without osteoarthritis (marked control) (N-OA), using the Gene Expression Omnibus. The obtained genes were investigated for their function and gene pathway using GO and KEGG. The enrichment was attained by PPI with STRING. The final enriched gene list was identified using the DGIdb.
Figure 2 Work flow of statistics mining. (A) Text mining: on the pursuit concepts osteoarthritis (OA) and patients without osteoarthritis (N-OA), text digging was accomplished with GEO, and 62,976 genes were established. (B) Gene enrichment: gene enrichment analysis was accomplished with GO and KEGG to enrich 66 genes. Meanwhile, by uploading Search Tool for the STRING, 66 hub genes were enriched, and among them eight significant genes were targetable by drugs using DGIdb. (C) Drug-gene interaction: these eight genes were established in DGIdb and the inhibitor Selenocysteine drugs were recognized as having a prospective impact on cartilage regeneration.
Materials and Methods
mRNA Microarray Statistics and Database Acquisition
The mRNA statistics text microarray GSE117999 (Rai MF et al., 2018) (the Patients/tissues details are given in the Supplementary Material) were download on the community GEO database (Barrett et al., 2013) (https://www.ncbi.nlm.nih.gov/geo), and executed on the GPL20844 platform. GSE117999 contains 12 patients with osteoarthritis, and 12 patients without osteoarthritis (arthroscopic partial meniscectomy). The normalized log2 ratio expressive OA/N-OA of the GSE117999 dataset, normalized by the L and Q process (Sbarrato et al., 2017), was taken. Analysis interested gene IDs were changed to gene symbols on the Agilent-072363 SurePrint G3 Human GE v3 8x60K Microarray 039494 (Feature Number Version). Institutional Review Board approved the study protocol. Prior to participation, a written informed consent was obtained from each patient (Brophy et al., 2018).
DEG Analysis
DEG analysis is the genes expressed at meaningfully levels by numerous methods of analysis (Rahmatallah et al., 2016). We used limma and R to recognize the DEGs in the tissues of OA patients, and compared them to the patients without osteoarthritis (arthroscopic partial meniscectomy). Genes with |log2 FC| ≥1.2 and adjusted P values < 0.01 (moderated t-statistics, corrected by B and H method) were taken in the analysis period (Solari et al., 2017).
GO and KEGG Analysis
GO (Carbon et al., 2017; Thomas, 2017) comprises a set of terms defining gene produces of biological process (BP), molecular function (MF), and cellular component (CC). KEGG (Kanehisa et al., 2017) offers statistics of recognized biological pathways. We used DAVID (Dennis et al., 2003) to envision the biological function and pathways enrichment of DEGs (The noteworthy was determined as p-value < 0.05).
PPI Networks Construction
The STRING (Version 11.0) (Wei et al., 2019) catalogue was used. This databank includes over 25 million proteins and interactions, complexed within 5,000 organisms. The interaction score was set to ≥ 0.900, and the PPI (Protein-Protein Interaction) networks were created.
Drug-Gene Interactions
Drugs were selected based on the hub genes that served as promising targets by using the Drug-Gene Interaction Database (DGIdb; http://www.dgidb.org/search_interactions). In this study, the final drug was approved by the Food and Drug Administration (FDA) (Yang et al., 2020). This study was prepared so that data on drug-gene interactions and gene target ability could be acquired. Furthermore, PubChem (https://www.ncbi.nlm.nih.gov/pccompound/) was used to confirm whether the medicines identified in our enquiry could target the genes that we identified.
Quantitative Real-Time PCR (qRT-PCR)
After total RNA was extracted and mRNA purified (RNeasy Mini kit, Qiagen, Hilden, Germany), mRNA was converted to cDNA using the Trans-Scriptor First-Strand cDNA Synthesis SuperMix (TransScript, #AT301, Beijing, China). The assays-on-demand primers and probes and TaqMan Universal Master Mix were used to examine gene expression by the MiniOpticonTM RT-PCR system (Bio-Rad, Hercules, CA, USA) according to the user’s manual and amplified detection was conducted using SYBR-Green RealMastcrMix (Bio-Rad). The expression values of mRNAs were normalized against glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and relatively calculated using 2−ΔΔCt method (Ye et al., 2019).
Results
Identification of DEGs
In our analysis, we recognized DEGs from 12 paired cartilage tissues (OA), using R-limma, compared with N-OA (Figure 1/Figure 2). Using both |log2 FC| ≥ 1.2 and adjusted p value < 0.01, 66 genes were recognized, including 35 UP-GENEs and 31 Down-GENEs (Table 1) (Figure 3) (Supplementary Table).
Table 1 Sixty-six DEG were recognized in GSE117999.
Figure 3 Volcano Plot image displaying the 66 genes that were significantly up-regulated or down-regulated (P-value < 0.05) in OA cartilage, compared to normal controls.
GO and KEGG Analysis
To discover further prospective targets of those DEGs in cartilage regeneration, we performed GO and KEGG analyses on OA, using a p value of <0.05 (Figure 4). As revealed in Figure 4, it demonstrations the all noteworthy terms for each of the following: the BP, CC, MF, and KEGG pathways of DEGs.
Figure 4 GO and KEGG pathways of DEGs meaningfully enriched in OA.
The marks of the UP-GENEs and the DOWN-GENEs are shown. As revealed in Table 2, in the BP group, the UP-GENEs were primarily enriched for genes complicated in the intracellular oestrogen receptor signalling pathway, innate immune response in mucosa and antibacterial humoral response. The DOWN-GENEs were enriched genes in the antigen processing and the presentation of peptide or polysaccharide antigen, via MHC class II. They were also enriched in antigen processing and presentation, and the interferon-gamma-mediated signalling pathway. In the CC group, the UP-GENEs were chiefly enriched for genes associated with the vital constituent of the plasma membrane. The DOWN-GENEs were enriched for genes associated with the MHC class II protein complex, the integral component of lumenal side of endoplasmic reticulum membrane, and nuclear bodies. In the MF group, the UP-GENEs were chiefly enriched from genes associated with transcription factor activity and sequence-specific DNA binding. The downregulated DEGs were enriched for genes with MHC class II receptor activity. In the KEGG pathway group, the DOWN-GENEs were enriched for genes in Asthma, Graft-versus-host disease, and Allograft rejection.
Table 2 The GO and KEGG Enrichment Analysis of DEGs.
PPI Networks and Analysis
The DEGs were uploaded in STRING, and the results were explored. Eight HUB Genes with scores of >0.900 (maximum confidence) were chosen to create the PPI networks: CAMP and TOLLIP, DEFA3, HLA-DQA2, SLC38A6, SLC3A1, FAM20A and ANO8 (Figure 5). CAMP translates the set of an antimicrobial peptide family, characterized by chemoattractant, immune mediator induction, and immunoreaction regulation (Snoussi et al., 2018; Uysal et al., 2019). Toll interacting protein (TOLLIP) encodes a ubiquitin-binding protein that regulates inflammatory signalling (Diao et al., 2016; Shah et al., 2017). Defensin alpha 3 (DEFA3) belongs to antimicrobial and cytotoxic peptides involved in host defence (Okazaki et al., 2007; Froy and Sthoeger, 2009). Major histocompatibility complex, class II, DQ alpha 2 (HLA-DQA2) are a family of the HLA class II alpha family. Many investigators suggest this involved in the release of CLIP molecule (Rudy and Lew, 1994). Solute carrier family 38 member 6 (SLC38A6) has possibly a relation with the glutamate-glutamine cycle regulation, responsible for preventing excitotoxicity (Schioth et al., 2013; Bagchi et al., 2014). Solute carrier family 3 member 1 (SLC3A1) translates a type II membrane glycoprotein that encodes neutral amino acids associated with cystinuria (Ma et al., 2018). FAM20A golgi is associated with the secretory pathway pseudokinase (FAM20A), which encodes a protein that might function in haematopoiesis and is associated with amelogenesis imperfecta and gingival hyperplasia syndrome (Beres et al., 2018; Koruyucu et al., 2018). Anoctamin 8 (ANO8) is associated with a human disorder that is often overexpressed in diverse cancers (Katoh and Katoh, 2005; Ousingsawat et al., 2011).
Figure 5 PPI networks of DEGs. Sixty-six genes were clarified into the PPI network using the STRING databank. ↑ indicates up-Genes; ↓ indicates down-Genes.
Drug-Gene Interactions
The functional enrichment investigation was conducted to identify the last prospective HUB genes with the sifter stricture, p value <1X10^10, as the edge in the GenCLiP 2.0 website (Wang et al., 2019) (Figure 6). DGIdb database was used to hunt the potential targets CAMP and its small organic compounds Selenocysteine. Selenocysteine is the main form of selenium in proteins. The study of selenocysteine biosynthesis and the mechanism of protein participation is a classic protein biochemical. This important supplement of molecular biology is also the basis for further research into the biological functions and applications of selenoproteins (Stadtman, 1996). The structure of Selenocysteine is shown in Figure 7.
Figure 6 The functional enrichment analysis of eight applicant cartilage regeneration - associated genes (green shows that the gene is enhanced in the related pathway, black shows not on the contrary).
Figure 7 The schematic structure of Selenocysteine.
The Expression of CAMP was Inverse in OA Cartilage Tissues
To examine the expression of mRNAs in the articular cartilage, qRT-PCR was performed in 3 samples. We found that the expression of CAMP was significantly upregulated in the OA group compared with the normal group (Figure 8).
Figure 8 The expression of CAMP transcripts differentially expressed between normal and osteoarthritis was validated by real-time PCR. The data is presented as mean with standard error of the mean.
Discussion
This analysis was designed to recognize prospective cartilage regeneration - associated genes, by comparing cartilage tissues in patients without osteoarthritis (arthroscopic partial meniscectomy). Thirty-five UP and 31 DOWN DEGs were recognized. We accomplished GO and KEGG comment examinations. Next, a PPI network was created, and eight significant genes were recognized. Finally, eight genes, CAMP and TOLLIP, DEFA3, HLA-DQA2, SLC38A6, SLC3A1, FAM20A, and ANO8 were identified as being significantly associated with immune response, immune mediator induction, and cell chemotaxis. The CAMP inhibitor Selenocysteine may be a nanomedicine potential candidate for cartilage regeneration.
In summary, only the CAMP gene has its commercial inhibitor of the eight hub genes, there have been many studies into the role of transcription factors in regulating the promoter region of the CAMP gene (Horibe et al., 2013). These studies have revealed that many biomolecules and factors can regulate the expression of the human antibacterial peptide CAMP gene, and that different signalling pathways can also affect the expression of the gene (Hancock and Diamond, 2000). Due to the important role that the CAMP protein plays in infectious diseases, tumors, and other diseases, it may be a target for the diagnosis and treatment of these diseases (Kovach et al., 2012). It is expected that the molecular regulatory mechanism of the CAMP gene during the occurrence of diseases, especially infectious diseases, will become a topic of significant further research (Rosenberger et al., 2004).
Selenocysteine is the main form of selenium in proteins. The determination of its codon UGA increases the number of amino acids that make up the biosynthetic protein, from 20 to 21 (Mousa et al., 2017). It is also the only amino acid that contains a metalloid element. It is mostly located in the active center of selenoprotein or selenoenzyme (especially in antioxidant enzymes). At the same time, as an essential trace element of the human body, it has a very significant role in anti-oxidation, immune regulation, and anti-tumor roles (Hatfield et al., 2014; Varlamova and Cheremushkina, 2017). This important supplement of molecular biology is also the basis for further research into the biological functions and applications of selenoproteins (Stadtman, 1996). It is one of the hotspots of research in the fields of protein biochemistry and molecular biology (Ren et al., 2018).
Interestingly, Se deficiency has been proposed as an underlying contributing factor for the chronic osteochondral disease Kashin-Beck disease, an important but neglected disease in parts of China. It was previously shown that Se(IV), as well as superoxide dismutase, can prevent damage done to cultured human embryonic cartilage cells caused by various etiological environmental substrates, and increase the activity of GSHpx while decreasing the production of lipid peroxides (Peng and Yang, 1991). In another report, the disease was associated with the incidence of Se deficiency in regions where the disease is prevalent (Peng et al., 1992). Finally, a mixture containing glycosaminoglycans, selenium, and vitamin E was shown to be exceptionally capable of promoting osteochondral repair in a rabbit model of knee osteochondral defect after 6 weeks of treatment (Handl et al., 2007). These past studies provide some evidence to support our current finding that selenocysteine may be effective drug to support cartilage regeneration.
To conclude, 66 prospective candidate cartilage regeneration - associated genes, which have earlier been involved in numerous pathways related to pathogenesis. All of these DEG applicant cartilage regeneration -associated genes should be further established through biological trials. Moreover, CAMP, TOLLIP, DEFA3, HLA-DQA2, SLC38A6, SLC3A1, FAM20A, and ANO8, as prospective markers for cartilage regeneration; they have not been linked previously, either in diagnosis or in research. Thus, the eight targets were considered to be potential therapeutic targets of cartilage regeneration. Thus, the CAMP inhibitor selenocysteine is considered to be a potential nanomedicine candidate for regenerative medicine.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
All authors contributed to the data analysis. JYe and BX completed the analysis using bioinformatics instruments. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (ID. 5192010500) and the National Key R&D Program of China (No.2017YFB1303000).
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.
Acknowledgments
We would like to thank Editage (www.editage.com) for English language editing.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2020.01138/full#supplementary-material
References
Alcidi, L., Beneforti, E., Maresca, M., Santosuosso, U., Zoppi, M. (2007). Low power radiofrequency electromagnetic radiation for the treatment of pain due to osteoarthritis of the knee. Reumatismo 59, 140–145. doi: 10.4081/reumatismo.2007.140
Bagchi, S., Baomar, H. A., Al-Walai, S., Al-Sadi, S., Fredriksson, R. (2014). Histological analysis of SLC38A6 (SNAT6) expression in mouse brain shows selective expression in excitatory neurons with high expression in the synapses. PloS One 9, e95438. doi: 10.1371/journal.pone.0095438
Barrett, T., Wilhite, S. E., Ledoux, P., Evangelista, C., Kim, I. F., Tomashevsky, M., et al. (2013). NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res. 41, D991–D995. doi: 10.1093/nar/gks1193
Beres, F., Lignon, G., Rouziere, S., Mauprivez, C., Simon, S., Berdal, A., et al. (2018). Physicochemical analysis of human pulpal mineralization secondary to FAM20A mutations. Connect. Tissue Res. 59, 46–51. doi: 10.1080/03008207.2018.1435644
Brophy, R. H., Zhang, B., Cai, L., Wright, R. W., Sandell, L. J., Rai, M. F. (2018). Transcriptome comparison of meniscus from patients with and without osteoarthritis. Osteoarthritis Cartilage 26, 422–432. doi: 10.1016/j.joca.2017.12.004
Carbon, S., Dietze, H., Lewis, S. E., Mungall, C. J., Munoz-Torres, M. C., Basu, S. (2017). Expansion of the Gene Ontology knowledgebase and resources. Nucleic Acids Res. 45, D331–d338. doi: 10.1093/nar/gkw1108
Chen, H., Sun, T., Yan, Y., Ji, X., Sun, Y., Zhao, X., et al. (2020). Cartilage matrix-inspired biomimetic superlubricated nanospheres for treatment of osteoarthritis. Biomaterials 242, 119931. doi: 10.1016/j.biomaterials.2020.119931
Cheng, B., Liu, Y., Zhao, Y., Li, Q., Liu, Y., Wang, J., et al. (2019). The role of anthrax toxin protein receptor 1 as a new mechanosensor molecule and its mechanotransduction in BMSCs under hydrostatic pressure. Sci. Rep. 9, 12642–12642. doi: 10.1038/s41598-019-49100-5
Cheng, P. F. (2018). Medical bioinformatics in melanoma. Curr. Opin. Oncol. 30, 113–117. doi: 10.1097/CCO.0000000000000428
Cotto, K. C., Wagner, A. H., Feng, Y. Y., Kiwala, S., Coffman, A. C., Spies, G., et al. (2018). DGIdb 3.0: a redesign and expansion of the drug-gene interaction database. Nucleic Acids Res. 46, D1068–d1073. doi: 10.1093/nar/gkx1143
Dennis, G., Jr., Sherman, B. T., Hosack, D. A., Yang, J., Gao, W., Lane, H. C., et al. (2003). DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 4, P3. doi: 10.1186/gb-2003-4-5-p3
Diao, N., Zhang, Y., Chen, K., Yuan, R., Lee, C., Geng, S., et al. (2016). Deficiency in Toll-interacting protein (Tollip) skews inflamed yet incompetent innate leukocytes in vivo during DSS-induced septic colitis. Sci. Rep. 6, 34672. doi: 10.1038/srep34672
Froy, O., Sthoeger, Z. M. (2009). Defensins in systemic lupus erythematosus. Ann. N. Y. Acad. Sci. 1173, 365–369. doi: 10.1111/j.1749-6632.2009.04622.x
Gao, F., Xu, Z., Liang, Q., Li, H., Peng, L., Wu, M., et al. (2019). Osteochondral Regeneration with 3D-Printed Biodegradable High-Strength Supramolecular Polymer Reinforced-Gelatin Hydrogel Scaffolds. Adv. Sci. (Weinh) 6, 1900867–1900867. doi: 10.1002/advs.201900867
Hancock, R. E., Diamond, G. (2000). The role of cationic antimicrobial peptides in innate host defences. Trends Microbiol. 8, 402–410. doi: 10.1016/S0966-842X(00)01823-0
Handl, M., Amler, E., Bräun, K., Holzheu, J., Trc, T., Imhoff, A. B., et al. (2007). Positive effect of oral supplementation with glycosaminoglycans and antioxidants on the regeneration of osteochondral defects in the knee joint. Physiol. Res. 56, 243–249. doi: 10.1007/s000130050125
Hasman, A., Ammenwerth, E., Dickhaus, H., Knaup, P., Lovis, C., Mantas, J., et al. (2011). Biomedical informatics–a confluence of disciplines? Methods Inf. Med. 50, 508–524. doi: 10.3414/ME11-06-0003
Hatfield, D. L., Tsuji, P. A., Carlson, B. A., Gladyshev, V. N. (2014). Selenium and selenocysteine: roles in cancer, health, and development. Trends Biochem. Sci. 39, 112–120. doi: 10.1016/j.tibs.2013.12.007
Horibe, K., Nakamichi, Y., Uehara, S., Nakamura, M., Koide, M., Kobayashi, Y., et al. (2013). Roles of cathelicidin-related antimicrobial peptide in murine osteoclastogenesis. Immunology 140, 344–351. doi: 10.1111/imm.12146
Hu, P.-F., Sun, F.-F., Jiang, L.-F., Bao, J.-P., Wu, L.-D. (2018). Paeoniflorin inhibits IL-1β-induced MMP secretion via the NF-κB pathway in chondrocytes. Exp. Ther. Med. 16, 1513–1519. doi: 10.3892/etm.2018.6325
Huang, H., Zheng, J., Shen, N., Wang, G., Zhou, G., Fang, Y., et al. (2018). Identification of pathways and genes associated with synovitis in osteoarthritis using bioinformatics analyses. Sci. Rep. 8, 10050. doi: 10.1038/s41598-018-28280-6
Jalal, H., Pechlivanoglou, P., Krijkamp, E., Alarid-Escudero, F., Enns, E., Hunink, M. G. M. (2017). An Overview of R in Health Decision Sciences. Med. Decis. Making 37, 735–746. doi: 10.1177/0272989X16686559
Kanehisa, M., Furumichi, M., Tanabe, M., Sato, Y., Morishima, K. (2017). KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 45, D353–d361. doi: 10.1093/nar/gkw1092
Katoh, M., Katoh, M. (2005). Identification and characterization of TMEM16H gene in silico. Int. J. Mol. Med. 15, 353–358. doi: 10.3892/ijmm.15.2.353
Koruyucu, M., Seymen, F., Gencay, G., Gencay, K., Tuna, E. B., Shin, T. J., et al. (2018). Nephrocalcinosis in Amelogenesis Imperfecta Caused by the FAM20A Mutation. Nephron 139, 189–196. doi: 10.1159/000486607
Kovach, M. A., Ballinger, M. N., Newstead, M. W., Zeng, X., Bhan, U., Yu, F. S., et al. (2012). Cathelicidin-related antimicrobial peptide is required for effective lung mucosal immunity in Gram-negative bacterial pneumonia. J. Immunol. 189, 304–311. doi: 10.4049/jimmunol.1103196
Lawrence, K. M., Jones, R. C., Jackson, T. R., Baylie, R. L., Abbott, B., Bruhn-Olszewska, B., et al. (2017). Chondroprotection by urocortin involves blockade of the mechanosensitive ion channel Piezo1. Sci. Rep. 7, 5147–5147. doi: 10.1038/s41598-017-04367-4
Li, X., Ding, J., Zhang, Z., Yang, M., Yu, J., Wang, J., et al. (2016). Kartogenin-Incorporated Thermogel Supports Stem Cells for Significant Cartilage Regeneration. ACS Appl. Mater. Interfaces 8, 5148–5159. doi: 10.1021/acsami.5b12212
Liu, H., Cheng, Y., Chen, J., Chang, F., Wang, J., Ding, J., et al. (2018). Component effect of stem cell-loaded thermosensitive polypeptide hydrogels on cartilage repair. Acta Biomater. 73, 103–111. doi: 10.1016/j.actbio.2018.04.035
Ma, Y. Y., Liu, Y. P., Li, D., Li, X. Y., Song, J. Q., Yang, Y. L. (2018). Clinical, Biochemical, and Genetic Findings of Cystinuria in Chinese Children. Clin. Lab. 64, 1145–1151. doi: 10.7754/Clin.Lab.2018.180110
Math, K. R., Zaidi, S. F., Petchprapa, C., Harwin, S. F. (2006). Imaging of total knee arthroplasty. Semin. Musculoskelet. Radiol. 10, 47–63. doi: 10.1055/s-2006-934216
Mousa, R., Notis Dardashti, R., Metanis, N. (2017). Selenium and Selenocysteine in Protein Chemistry. Angew. Chem. Int. Ed. Engl. 56, 15818–15827. doi: 10.1002/anie.201706876
Okazaki, T., Ota, Y., Yuki, N., Hayashida, A., Shoda, A., Nishikawa, M., et al. (2007). Plasma levels of alpha-defensins 1-3 are an indicator of neutrophil activation in pregnant and post-partum women. J. Obstet. Gynaecol Res. 33, 645–650. doi: 10.1111/j.1447-0756.2007.00626.x
Ousingsawat, J., Kongsuphol, P., Schreiber, R., Kunzelmann, K. (2011). CFTR and TMEM16A are separate but functionally related Cl- channels. Cell Physiol. Biochem. 28, 715–724. doi: 10.1159/000335765
Pang, H., Wu, X.-H., Fu, S.-L., Luo, F., Zhang, Z.-H., Hou, T.-Y., et al. (2013). Prevascularisation with endothelial progenitor cells improved restoration of the architectural and functional properties of newly formed bone for bone reconstruction. Int. Orthop 37, 753–759. doi: 10.1007/s00264-012-1751-y
Peng, A., Yang, C. L. (1991). Examination of the roles of selenium in the Kaschin-Beck disease. Cartilage cell test and model studies. Biol. Trace Elem Res. 28, 1–9. doi: 10.1007/BF02990457
Peng, A., Yang, C., Rui, H., Li, H. (1992). Study on the pathogenic factors of Kashin-Beck disease. J. Toxicol. Environ. Health 35, 79–90. doi: 10.1080/15287399209531597
Ragni, E., Perucca Orfei, C., De Luca, P., Colombini, A., Viganò, M., Lugano, G., et al. (2019). Identification of miRNA Reference Genes in Extracellular Vesicles from Adipose Derived Mesenchymal Stem Cells for Studying Osteoarthritis. Int. J. Mol. Sci. 20, 1108. doi: 10.3390/ijms20051108
Rahmatallah, Y., Emmert-Streib, F., Glazko, G. (2016). Gene set analysis approaches for RNA-seq data: performance evaluation and application guideline. Brief Bioinform. 17, 393–407. doi: 10.1093/bib/bbv069
Ren, X., Zou, L., Lu, J., Holmgren, A. (2018). Selenocysteine in mammalian thioredoxin reductase and application of ebselen as a therapeutic. Free Radic. Biol. Med. 127, 238–247. doi: 10.1016/j.freeradbiomed.2018.05.081
Renders, G. A. P., Mulder, L., Lin, A. S., Langenbach, G. E. J., Koolstra, J. H., Guldberg, R. E., et al. (2014). Contrast-enhanced microCT (EPIC-μCT) ex vivo applied to the mouse and human jaw joint. Dentomaxillofac. Radiol. 43, 20130098–20130098. doi: 10.1259/dmfr.20130098
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43, e47. doi: 10.1093/nar/gkv007
Rosenberger, C. M., Gallo, R. L., Finlay, B. B. (2004). Interplay between antibacterial effectors: a macrophage antimicrobial peptide impairs intracellular Salmonella replication. Proc. Natl. Acad. Sci. U.S.A. 101, 2422–2427. doi: 10.1073/pnas.0304455101
Rudy, G., Lew, A. M. (1994). Limited polymorphism of the HLA-DQA2 promoter and identification of a variant octamer. Hum. Immunol. 39, 225–229. doi: 10.1016/0198-8859(94)90264-X
Sbarrato, T., Spriggs, R. V., Wilson, L., Jones, C., Dudek, K., Bastide, A., et al. (2017). An improved analysis methodology for translational profiling by microarray. Rna 23, 1601–1613. doi: 10.1261/rna.060525.116
Schioth, H. B., Roshanbin, S., Hagglund, M. G., Fredriksson, R. (2013). Evolutionary origin of amino acid transporter families SLC32, SLC36 and SLC38 and physiological, pathological and therapeutic aspects. Mol. Aspects Med. 34, 571–585. doi: 10.1016/j.mam.2012.07.012
Shah, J. A., Musvosvi, M., Shey, M., Horne, D. J., Wells, R. D., Peterson, G. J., et al. (2017). A Functional Toll-Interacting Protein Variant Is Associated with Bacillus Calmette-Guerin-Specific Immune Responses and Tuberculosis. Am. J. Respir. Crit. Care Med. 196, 502–511. doi: 10.1164/rccm.201611-2346OC
Snoussi, M., Talledo, J. P., Del Rosario, N. A., Mohammadi, S., Ha, B. Y., Kosmrlj, A., et al. (2018). Heterogeneous absorption of antimicrobial peptide LL37 in Escherichia coli cells enhances population survivability. Elife 7, 1–21. doi: 10.7554/eLife.38174
Solari, A., Goeman, J. J., Minimally adaptive, B. H. (2017). A tiny but uniform improvement of the procedure of Benjamini and Hochberg. Biom. J. 59, 776–780. doi: 10.1002/bimj.201500253
Stadtman, T. C. (1996). Selenocysteine. Annu. Rev. Biochem. 65, 83–100. doi: 10.1146/annurev.bi.65.070196.000503
Thomas, P. D. (2017). The Gene Ontology and the Meaning of Biological Function. Methods Mol. Biol. 1446 , 15–24. doi: 10.1007/978-1-4939-3743-1_2
Uehara, K., Hara, M., Matsuo, T., Namiki, G., Watanabe, M., Nomura, Y. (2015). Hyaluronic acid secretion by synoviocytes alters under cyclic compressive load in contracted collagen gels. Cytotechnology 67, 19–26. doi: 10.1007/s10616-013-9669-9
Uysal, P., Simsek, G., Durmus, S., Sozer, V., Aksan, H., Yurt, S., et al. (2019). Evaluation of plasma antimicrobial peptide LL-37 and nuclear factor-kappaB levels in stable chronic obstructive pulmonary disease. Int. J. Chron. Obstruct. Pulmon. Dis. 14, 321–330. doi: 10.2147/COPD.S185602
Varlamova, E. G., Cheremushkina, I. V. (2017). Contribution of mammalian selenocysteine-containing proteins to carcinogenesis. J. Trace Elem Med. Biol. 39, 76–85. doi: 10.1016/j.jtemb.2016.08.004
Wang, C., Feng, N., Chang, F., Wang, J., Yuan, B., Cheng, Y., et al. (2019). Injectable Cholesterol-Enhanced Stereocomplex Polylactide Thermogel Loading Chondrocytes for Optimized Cartilage Regeneration. Adv. Healthc. Mater. 8, e1900312. doi: 10.1002/adhm.201900312
Wang, J. H., Zhao, L. F., Wang, H. F., Wen, Y. T., Jiang, K. K., Mao, X. M., et al. (2019). GenCLiP 3: mining human genes’ functions and regulatory networks from PubMed based on co-occurrences and natural language processing. Bioinformatics. 36. (6), 1973–1975. doi: 10.1093/bioinformatics/btz807
Wei, S., Yun, X., Ruan, X., Wei, X., Zheng, X., Gao, M. (2019). Identification of potential pathogenic candidates or diagnostic biomarkers in papillary thyroid carcinoma using expression and methylation profiles. Oncol. Lett. 18, 6670–6678. doi: 10.3892/ol.2019.11059
Yang, D., He, Y., Wu, B., Deng, Y., Wang, N., Li, M., et al. (2020). Integrated bioinformatics analysis for the screening of hub genes and therapeutic drugs in ovarian cancer. J. Ovarian Res. 13, 10. doi: 10.1186/s13048-020-0613-2
Yao, Q., Wei, B., Liu, N., Li, C., Guo, Y., Shamie, A. N., et al. (2015). Chondrogenic regeneration using bone marrow clots and a porous polycaprolactone-hydroxyapatite scaffold by three-dimensional printing. Tissue Eng. Part A. 21, 1388–1397. doi: 10.1089/ten.tea.2014.0280
Ye, J., Zha, J., Shi, Y., Li, Y., Yuan, D., Chen, Q., et al. (2019). Co-inhibition of HDAC and MLL-menin interaction targets MLL-rearranged acute myeloid leukemia cells via disruption of DNA damage checkpoint and DNA repair. Clin. Epigenet. 11, 137. doi: 10.1186/s13148-019-0723-0
Keywords: cartilage, biomaterial drug, nanomedicine, immune response, clinical translation, text mining
Citation: Ye J, Xu B, Fan B, Zhang J, Yuan F, Chen Y, Sun Z, Yan X, Song Y, Song S, Yang M and Yu J-K (2020) Discovery of Selenocysteine as a Potential Nanomedicine Promotes Cartilage Regeneration With Enhanced Immune Response by Text Mining and Biomedical Databases. Front. Pharmacol. 11:1138. doi: 10.3389/fphar.2020.01138
Received: 24 March 2020; Accepted: 13 July 2020;
Published: 24 July 2020.
Edited by:
Salvatore Salomone, University of Catania, ItalyReviewed by:
Daniel Torres, University of Hawaii at Manoa, United StatesGuangdong Zhou, Shanghai Jiao-Tong University School of Medicine, China
Copyright © 2020 Ye, Xu, Fan, Zhang, Yuan, Chen, Sun, Yan, Song, Song, Yang and Yu. 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: Jia-Kuo Yu, yujiakuo@126.com
†These authors have contributed equally to this work