1.Number of crawled literature
A total of 3,352 literature related to photobiomodulation (PBM group), 20,442 literature related to photodynamic therapy (PDT group), 5,874 LED-related literature (LED group), and 326,426 articles on laser-related were crawled. (laser group).
2.Gene identification
A total of 146 human genes related to PBM, 268 human genes related to PDT were identified, and the frequency of top10 genes was shown in Figure 2a, 2b. The 142 human genes related to LED and 3188 human genes related to laser were identified.
3.Confirm the genes in LED/laser-related phototherapy
In this experiment, the intersected genes of the LED group genes and (PBM + PDT) group genes were LED-related phototherapy. In the same way, the genes in laser-related phototherapy were confirmed. By using the Venn diagram, as shown in Figure 3, it was clear that there were 161 human genes related to LED, of which 50 genes were intersected with PBM, 54 genes were intersected with PDT, and 29 genes were repeated intersecting. The 3188 human genes related to laser light therapy, of which 138 genes were intersected with PBM, 141 genes were intersected with PDT, and 59 genes were repeated intersecting. In summary, by crawling the literature and identifying genes, a total of 75 phototherapy-related genes with LED and 220 phototherapy-related genes with laser were found. Moreover, the genes were deserved to be attention in the public intersection of LED, laser, PBM, and PDT. We proposed that the genes in the public intersection may be the key genes that play a role in phototherapy. For this reason, enrichment analysis was performed on these genes.
4.Comparison of PBM and PDT
In the R language, the clusterprofiler package was used to enrich functions and pathways. Among them, a total of 57 pathways in the PBM group were enriched, and the results of enrichment were shown in Figure 4a, b; a total of 109 pathways in the PDT group were enriched. The results of pathway and function enrichment were shown in Figure 5a, b; Compare the two results, a total of 45 pathways were in the intersection, which intersection is 78.9% of the KEGG result of PBM.
5.Comparison of LED and laser
By crawling the literature and identifying genes, a total of 75 phototherapy-related genes with LED and 220 phototherapy-related genes with laser were found. No KEGG result in the LED-related phototherapy, which was caused by the number of identified genes was not enough. A total of 104 pathways were enriched in the laser group, as shown in Figure 6. However, it can be known that 74 genes were intersected with laser in 75 LED-related phototherapy genes from the Venn diagram. The intersected genes account for up to 98.7% of the 91 LED genes. Therefore, it can be concluded that the genes affected by the LED in phototherapy are the same as the laser light source.
6.The 29 genes that were located in the common intersection of PBM, PDT, LED, and laser
Known from the Venn diagram that a total of 29 genes were a common intersection in these four groups, and so we speculated that these genes may be the keys in phototherapy. Therefore, we also performed KEGG and GO analysis for these 29 genes. A total of 7 pathways were enriched. The results were shown in Figure 7a. At the same time, we also performed a GO enrichment analysis for these genes, and a total of 365 results were enriched, as shown in Figure 7b. It was obvious that these genes were related to tyrosine, ultraviolet, light, and oxygen stress, etc.
7.Sequencing results of LED stimulate-MH7A
To explore the mechanisms of PBM for RA disease, we sequenced LED-MH7A cells. First, MH7A cells were stimulated with TNF-α. The experimental groups were irradiated with 630nm LED for 15min at room temperature, while control groups were wrapped with aluminum foil at room temperature. A total of 2950 differentially expressed genes were identified. 1885 and 1065 differentially expressed genes were up-regulated,down-regulated, respectively. The volcano is shown in Figure 8. KEGG enriched to 29 signal pathways, as shown in Figure 9. Among the 29 signaling pathways, 12 pathways were matched with those obtained by PBM analysis in this paper, which are Cytokine-cytokine receptor interaction, Legionellosis, Prion diseases, Viral protein interaction with cytokine and cytokine receptor, Complement and coagulation cascades, Bladder cancer, Fluid shear stress, and atherosclerosis, Prostate cancer, MAPK signaling pathway, TNF signaling pathway, IL-17 signaling pathway, Hepatitis B. All data, including sequencing reads and single-nucleus expression matrices have been deposited in NCBI’s Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/). The accession number for the mRNA-seq data reported in this study is GSE172223.
8.PPI network and cluster analysis
For further analysis, we got 12 pathways that were matched in the KEGG results of MH7A and the KEGG results of PBM. The pathways were Cytokine-cytokine receptor interaction, Legionellosis, Prion diseases, Viral protein interaction with cytokine and cytokine receptor, Complement and coagulation cascades, Bladder cancer, Fluid shear stress and atherosclerosis, Prostate cancer, MAPK signaling pathway, TNF signaling pathway, IL-17 signaling pathway, Hepatitis B. The genes were constructed a PPI network in 12 pathways. The outcomes of the string profile obtained from STRING showed that contained 129 nodes and 1235 edges. The top 10 outstanding proteins were identified as hub genes, which might play critical roles in the PBM through PPI analysis of STRING software including IL6, TNF, VEGFA, IL10, EGFR, EGF, CCL2, MMP9, IL1B, FOS. Among these significant proteins, IL6 was the most important protein and connected with 84 nodes in the network. A total of 6 clusters were generated in MCODE, and the top one cluster was selected as hub modules by the scores evaluated in MCODE. MCODE 1, which had 32 nodes and 369 edges, had the highest score in these clusters, such as Figure 10. In this cluster, the 9 nodes, including IL6, VEGFA, CXCL2, CSF1P, ITGAM, TNF, FOS, MMP9, PGF were in both MH7A cells and the PBM group. Maybe, these 9 nodes are key genes for PBM. On the one hand, in this cluster, we can know that some protein(red circles) of PBM results obtained by the method proposed in this paper already coincide with the parts of bio-MH7A cell sequencing. It proves the feasibility of our method. On the other hand, in LED stimulate-MH7A cells, the yellow circles, including CSF3, CXCL3, SOCS3, HMOX1, HBEGF, IL1B, IL1A, PLAU, and CCL20, represent the differentially expressed genes and have not been studied. Therefore, these genes have more research value and potential.
9. Effects of 630-nm LED radiation on mRNA gene expressions of inflammatory cytokines
For the accuracy of the PPI network and the effect of 640nm LED red light on rheumatoid disease, we performed RT-qPCR to analysis mRNA levels of inflammatory.The result was shown in Figure 11; there was a significant decrease in the gene expressions of inflammatory cytokines IL-6, IL-8, and IL-1Βand MMP-3 expression; meanwhile, the expression of antiinflammatory factor IL-10 was promoted with the radiation doses increased. Collectively, these results demonstrate that the gene expression and production of inflammatory cytokines by MH7A cells can be inhibited by 630-nm LED radiation.