Yue Feng
Zhao Wang
Nan Yang
Sijia Liu
Jiazhuo Yan
Jiayu Song
Shanshan Yang
Yunyan Zhang- Department of Gynecological Radiotherapy, Harbin Medical University Cancer Hospital, Harbin, China
Cervical cancer as a common gynecological malignancy threatens the health and lives of women. Resistance to radiotherapy is the primary cause of treatment failure and is mainly related to difference in the inherent vulnerability of tumors after radiotherapy. Here, we investigated signature genes associated with poor response to radiotherapy by analyzing an independent cervical cancer dataset from the Gene Expression Omnibus, including pre-irradiation and mid-irradiation information. A total of 316 differentially expressed genes were significantly identified. The correlations between these genes were investigated through the Pearson correlation analysis. Subsequently, random forest model was used in determining cancer-related genes, and all genes were ranked by random forest scoring. The top 30 candidate genes were selected for uncovering their biological functions. Functional enrichment analysis revealed that the biological functions chiefly enriched in tumor immune responses, such as cellular defense response, negative regulation of immune system process, T cell activation, neutrophil activation involved in immune response, regulation of antigen processing and presentation, and peptidyl-tyrosine autophosphorylation. Finally, the top 30 genes were screened and analyzed through literature verification. After validation, 10 genes (KLRK1, LCK, KIF20A, CD247, FASLG, CD163, ZAP70, CD8B, ZNF683, and F10) were to our objective. Overall, the present research confirmed that integrated bioinformatics methods can contribute to the understanding of the molecular mechanisms and potential therapeutic targets underlying radiotherapy resistance in cervical cancer.
Introduction
Cervical cancer is the fourth most common malignancy and the fourth leading cause of cancer-related mortality in women globally (Barker et al., 2015; Doja et al., 2020; Sahu and Pattanayak, 2020). Among women between the ages of 20 and 39, cervical cancer is the second leading cause of cancer deaths (Siegel et al., 2020). To date, approximately 80% of patients with cervical cancer need radiotherapy in the process of treatment (Jemal et al., 2011). Radiotherapy resistance, also known as radiosensitivity, remains a significant hurdle to cervical cancer therapeutics. Since 1999, the cisplatin-based concurrent chemoradiotherapy (CCRT) is the standard of treatment for advanced cervical cancer (Rose et al., 1999). Compared with radiotherapy alone, CCRT increases patients’ local control rates and can improve prognosis. However, cisplatin increases the incidence of acute hematological toxicity, which may lead to the interruption of CCRT treatment and poor prognosis in patients with cervical cancer (Green et al., 2001; Dueñas-González et al., 2011; Bazan et al., 2013; Liu et al., 2020). The 5 years survival rate for advanced cervical cancer is still only 5–15% (Chopra et al., 2018). Therefore, exploring the mechanisms underlying radiosensitivity at the molecular level is the essential strategy for increasing long-term survival of cervical cancer.
Resistance to radiotherapy can be intrinsic or acquired, but its specific mechanism is not yet clear. As a promising predictive application, bioinformatics has led to breakthroughs in the field of medical research (Kanehisa and Bork, 2003; Wadlow and Ramaswamy, 2005; Oliver et al., 2015; Tang et al., 2018; Gauthier et al., 2019; Yang and Deng, 2020; Zhang, 2020; Zhang et al., 2020b). For decades, many bioinformatics studies have been performed to investigate radiotherapy resistance-related mechanisms and detect molecular biomarkers (Zhang et al., 2020c,d). However, in radiotherapy, a tumor may gradually adapt to changes in physical and chemical environments and acquire radiotherapy resistance (Baskar et al., 2012; Stapleton et al., 2017). Reliable studies that focus on changes in gene expression during radiotherapy treatment based on bioinformatics are few. Hence, key genes for revealing molecular mechanisms to overcome cervical cancer radiotherapy resistance according to genomic changes should be identified.
In the current study, we acquired a cervical cancer-associated gene expression dataset, comparing cervix cancer tissue between pre-irradiation and mid-irradiation status to investigate the phenomenon of changes in gene expression during radiotherapy. Second, we performed gene expression profile analysis and identification of differentially expressed genes (DEGs). Correlations between these genes were discovered through the Pearson correlation analysis. Subsequently, random forest was used in optimizing risk genes. According to the random forest scoring results, the top 30 genes were selected as candidate genes for verification analysis. Finally, a total of 10 signature genes were verified through comprehensive literature searching as biomarkers for radiotherapy resistance in cervical cancer.
Materials and Methods
The whole designation was performed according to the flow chart (Figure 1).
Figure 1. The overall workflow.
Microarray Data Collection and Preprocessing
The gene expression profile was downloaded from the Gene Expression Omnibus (GEO) database1 which includes a variety of gene expression status of public genomics data (Dueñas-González et al., 2011). We retrieved the keyword “cervical cancer” in this functional database to obtain relevant information about radiotherapy resistance. The inclusion criteria were as follows: (a) the sample was treated with CCRT and without neoadjuvant chemotherapy; (b) number of samples were not less than 50; and (c) sample data were reliable.
The raw GEO dataset of GSE3578 was downloaded as the validation set. It had a total of 156 samples (78 pre-irradiation and 78 mid-irradiation samples). The dataset was based on the platform GPL2895, GE Healthcare/Amer sham Biosciences Code Link Human Whole Genome Bio array. Next, we preprocessed the data and converted each probe ID into a gene symbol. Probes without corresponding or mapping to multiple gene symbols were rejected to prevent ambiguity (Zhang et al., 2018).
Gene Expression Profile Analysis and Identification of Differentially Expressed Genes
The limma package was used in processing the downloaded expression profile and converting and filtering unqualified data to obtain DEGs (Ritchie et al., 2015; Liu et al., 2019b). The original data were chip data, and thus the limma package was used for our analysis. DEGs between the pre-irradiation and mid-irradiation samples were identified, and a P value of <0.05 and | log2fold change (FC)| of >0.5 were used as the cutoff thresholds.
Pearson Correlation Analysis
Pearson correlation analysis is a statistical method with powerful functions for accurately analyzing gene or protein co-expression correlations (Zhang et al., 2020a; Tang et al., 2021). In this study, it was used in assessing the relationships among the 316 DEGs obtained. According to the results of the analysis, genes with Pearson correlation coefficient threshold absolute values of >0.8 were selected as candidate genes.
Random Forest
Random forest is a generally acknowledged ensemble classifier for machine learning and can exploit large data repositories for the analysis of risk predictors and their intimate interactions and advancement risk prediction capability (Yang et al., 2021). It was used in evaluating cervical cancer candidate genes as the main objects of our research (Nunn et al., 2012; Wang et al., 2018; Deng et al., 2019b; Ru et al., 2019; Lv et al., 2020; Yang et al., 2020). In the random forest setting, we divided the expression profile data of 103 genes according to the properties of the samples and considered the samples before radiotherapy as “good” and the samples during radiotherapy “poor,” respectively. Then, 156 samples from GSE3578 were randomly divided into three groups. One group (52 samples) was used as the training set, whereas the other two groups (104 samples) were as the test sets. All samples were used for random forest scoring (Rigatti, 2017; Wang and Zhou, 2017).
Enrichment Analysis
Gene ontology (GO) enrichment was performed for the annotation of the three types of functions of genes, namely, biological process (BP), cytological component, and molecular function (Peng et al., 2017; Deng et al., 2019a). The gene list containing the top 30 scores in the random forest was regarded as an input file and utilized for GO function enrichment in the R package. On the basis of the reconstructed enrichment analysis, the results of pathways on cervical cancer radiosensitivity were further found.
Screening and Analysis of the Top 30 Genes
Medical literature retrieval is the basic and authoritative verification method in the verification process. After obtaining 30 candidate genes, we extensively reviewed the literature to determine whether the selected genes are indeed signature genes for cervical cancer radiotherapy resistance. We retrieved the name of the genes with the keywords “cervical cancer,” “radiotherapy resistance,” “irradiation,” and “radiosensitivity” in the databases to obtain relevant information about the top 30 genes. Ultimately, we obtained the signature genes for cervical cancer radiotherapy resistance through comprehensive analysis.
Results
Identification of Differentially Expressed Genes
First, to investigate the radiosensitivity in cervical cancer, we processed and converted the raw gene microarray expression dataset GSE3578 associated with the cancer. We extracted 17,097 genes from the obtained microarray dataset. Significance analysis with the limma package indicated that 316 DEGs were present in the samples, of which 237 were obviously down-regulated and 79 were up-regulated.
Pearson Correlation Analysis of DEGs
Pearson correlation analysis was subsequently performed between 316 differential genes for the assessment of the association (Cheng et al., 2019). Genes with a correlation coefficient threshold greater than 0.8 in absolute value were selected on the basis of the results. The analysis results were provided in Figure 2, and 103 genes met the threshold setting requirements. All the small nodes were the differential genes we calculated and each node represented a gene with a correlation coefficient greater than 0.8. According to the Pearson correlation coefficient, we obtained the correlations between genes, that is, the greater the connectivity, the greater the node size of the gene displayed.
Figure 2. Results of Pearson correlation analysis. According to the Pearson correlation coefficient, we obtained the correlations between genes, that is, the greater the connectivity, the greater the node size of the gene displayed.
Optimization of Candidate Genes With Random Forest
Random forest classification is a universally accepted machine learning procedure for regression tree classification and regression and model prediction (Cheng et al., 2018; Sapir-Pichhadze and Kaplan, 2020; Ubels et al., 2020; Zhao et al., 2020). It was first proposed by Breiman in 2001. It takes advantage of binary splits on predictor variables for value outcome forecasts (Steffens et al., 2006). To achieve cervical cancer-related candidate genes, we built a random forest model. We divided the 156 samples into three parts, of which one was the training set (52 samples) and the other was the test set (104 samples). In the test set, the “good group” contained 51 samples, meanwhile, the “poor group” contained 56 samples. Regarding the differential genes as features, the test set was classified by random forest method. When the number of “ntrees” was 1,000, we obtained the most ideal training result with a classification error rate of 9.62%, of which five cases were misclassified in the “good group” and five cases were misclassified in the “poor group.” As the results were obtained through the classification of the test set, we regarded this classified procedure noble ideal.
Furthermore, we comprehensively utilized random forest scoring on all samples to further verify the accuracy of our data. The results revealed that the minimum error rate of the classification results was only 1.92% and 78 samples were in the “good group” and one was misclassified. The “poor group” contained 78 samples, including two misclassified samples. After several times of analysis, the error rate of random forest classification results was less than 5%, indicating that the differential genes as features can classify samples well. Therefore, sample classification with the genes had a low error rate, indicating that our classification results were of great significance and had the potential as molecular markers of cervical cancer radiotherapy resistance.
Mean decrease gini (MDG), one of importance measures in random forest scoring, was employed to rank variables and for variable selection. In this study, MDG marked the degree of contributions of genes that could distinguish the characteristics of samples before and during radiotherapy. As a result, all genes were sorted in descending order according to MDG, and the top 30 candidate genes for radiotherapy resistance in cervical cancer were selected and shown in Figure 3.
Figure 3. Top 30 genes in random forest scoring results. All genes were sorted in descending order according to Mean decrease gini (MDG), and the top 30 candidate genes for radiotherapy resistance in cervical cancer were selected.
Enrichment Functional Analysis of Candidate Genes
To further uncover the molecular biology functions regulated by the genes, the top-ranked 30 genes according to the random forest scoring results were subsequently employed for the functional analysis of BP. As shown in Figure 4, the results indicated that biological functions were chiefly enriched in the corresponding tumor immune process pathways with the high fold overexpression, including cellular defense response, negative regulation of immune system process, T cell activation, neutrophil activation involved in immune response, regulation of antigen processing and presentation, and peptidyl-tyrosine autophosphorylation. These immune biological functions are significantly correlated with the examination of radiotherapy resistance. Researches on the radiotherapy resistance of cervical cancer from the perspective of the tumor immune microenvironment is extensive (Barker et al., 2015), showing that tumor radiotherapy has a certain impact on tumor microenvironment (TME), and changes in the TME, such as the generation of tumor-associated macrophages (TAMs) and cancer-associated fibroblasts can in turn affect the sensitivity of radiotherapy (Chu et al., 2014; Zhang et al., 2017; Choo et al., 2018). This shows that our enrichment results are very valuable. Given that the satisfying pathways were enriched, we used the results for the follow-up verification of the signature genes.
Figure 4. Enrichment analysis results of biological process (BP) was constructed by R package. The biological functions chiefly enriched in tumor immune responses.
Validation of Candidate Genes in Literature Retrieval
Considering the specific functions of the candidate genes in the sensitivity of cervical cancer radiotherapy, we comprehensively conducted a systematic medical literature search based on the top 30 genes. After verification in the literature, the candidate genes ranked 1, 5, 7, 9, 11, 13, 14, 18, 24, and 29 were confirmed with relevant experiments with cancers, demonstrating that these genes were our signature genes for cervical cancer radiotherapy sensitivity. In details, Killer cell lectin like receptor K1 (KLRK1), lymphocyte protein tyrosine kinase (LCK), kinesin family member 20A (KIF20A), CD247 molecule (CD247), Fas ligand gene (FASLG), CD163 molecule (CD163), zeta-chain associated protein 70 (ZAP70), CD8b molecule (CD8B), zinc finger protein 683 (ZNF683), and coagulation factor X (F10) were the top 10 genes. The scores of the 11th and 13th genes were closely associated to cervical cancer radiotherapy resistance, the 7th, 24th, and 29th genes were associated with cervical cancer, the first gene was correlated with cervical intraepithelial neoplasia, the 5th, 14th, and 18th genes were related to radiosensitivity for other cancers, and the 9th gene was associated to tumor immune microenvironment. In conclusion, this suggests that signature genes are expected to be biomarkers for cervical cancer radiotherapy resistance based on RNA sequencing data.
Discussion
Although radiotherapy for cervical cancer has been widely used, the emergence of radioresistance has deeply affected its treatment effect (Kaidar-Person et al., 2012; Han et al., 2021). Consequently, overcoming the radiotherapy resistance of cervical cancer and finding related accurate biomarkers are essential and urgently needed. Different gene expression characteristics may lead to different types of tumors with varying degrees of sensitivity to radiotherapy (Jiang et al., 2017; Prevo et al., 2017; Zhou et al., 2017). With the development of bioinformatics, gene expression microarrays have been widely used in studying cancer-related genes and considered promising tools for molecular prediction (Cheng et al., 2017; Liu et al., 2019a; Xing et al., 2020).
In the present study, the data from GSE3578 were used in identifying 316 DEGs, which we regarded as the origins. The relationships among these genes were considered, and Pearson correlation analysis was employed. A total of 103 genes with correlation coefficient thresholds >0.8 in absolute value were used for subsequent gene optimization. Identifying cancer-related candidate genes through the random forest method is important. The expression profile data of 103 genes were classified according to the properties of the samples. Then, the 156 samples in GSE3578 were randomly divided into a training set and two test sets, and all the samples were used for random forest scoring. The error rate of sample classification based on these genes was sufficiently small, which can help us to distinguish between normal and disease samples.
Gene ontology term enrichment analysis was conducted using the top 30 genes for the recognition of the roles of DEGs in radiotherapy initiation and progression in cervical cancer. The BP enrichment results were closely related to radiotherapy resistance in cervical carcinoma, suggesting that cervical cancer radiotherapy and immune system are mutually influential. That is, these results showed the specific molecular mechanisms of radiosensitivity in cervical cancer to a certain degree. According to the enriched GO terms, these genes were selected for subsequent verification and analysis. After thoroughly verifying the above genes in the literature, we observed that 10 genes were significant. These genes were selected as our signature genes for cervical cancer radiotherapy sensitivity.
Many studies have been carried out on the 10 candidate genes, and the ranking of the 11th and 13th genes were the most meaningful for our purpose. First, the genes related to cervical cancer or radiotherapy sensitivity were explored. FASLG, a type II membrane protein in the tumor necrosis factor superfamily of death receptors, can combine with and trimerize the death receptor Fas on cell surfaces and thereby activate authoritative extrinsic pathways for apoptosis (Nagata, 1997; Fulda and Debatin, 2006). Cervical cancer is sensitive to radiotherapy because this gene has different levels of expression. Di et al. stated that the proportion of the pro-apoptotic protein FASLG is up-regulated by the radiotherapy sensitizer diallyl disulfide, which induce the activation of the apoptotic pathway in HeLa cells; they also showed that FASLG plays a direct role in the study of human cervical cancer cell apoptosis induced by radiosensitizers (Fulda and Debatin, 2006). In addition, radiotherapy plays an important role in the adjuvant treatment of glioma cancer; Huang et al. (2019) showed that the radiosensitivity of human glioma cells increases through the FASLG expression-associated signaling pathway. FASLG may be a candidate signature for the mechanism of radiotherapy response in cervical cancer.
Changes in the TME are closely related to the biological effects of radiation and are important factors that affect the sensitivity of tumor cells to radiotherapy (Barker et al., 2015). CD163, one of the immune cell markers for macrophages in the TME, is often used in determining M1 and M2 macrophages (de Vos van Steenwijk et al., 2013). The biological behavior of FASLG is closely correlated with the radiosensitivity of cervical cancer. Lippens et al. (2020) found that patients with reducing CD163 scores showed more pathological complete response, besides, patients with high levels of CD163 cells before treatment that were decreased after therapy had better prognosis. Another study focused on impacts on TME subjected to radiotherapy against head and neck squamous cell carcinoma and pointed out the up-regulation of M2 macrophage-related genes (CD163 and CD206) in irradiated tumors (Moreira et al., 2021). M2 phenotype TAMs are associated with decrease in disease-free survival in advanced breast cancer, and the M2-shifted subgroup is unsuitable for partial breast irradiation (Schnellhardt et al., 2020). Radiation therapy, in addition to killing cancer cells, affects the TME, particularly contributing to tumor blood vessel damage and immunosuppressive pathway activation (Jarosz-Biej et al., 2019). The accumulation of cells can lead to the suppression of radioactive response specifically in TAMs with the M2 phenotype, which facilitates the survival of cancer cells (Fridman et al., 2017). Previous reports and our results have pointed out that CD163 may have a key role in exploring the effectiveness of radiotherapy in cervical cancer radiotherapy.
Four of the 10 candidate genes have been studied in cervical cancer, but their correlations with cervical cancer radiotherapy is unclear (Verhey and Hammond, 2009). KIF20A has a microtubule-dependent force, participating in different cellular processes, such as the formation of the biological behavior of spindle structures and chromosome composition. Overexpressed KIF20A increases the risk of cervical cancer. Through bioinformatics methods and meta-analysis, KIF20A expression was found to significantly vary among 363 normal tissues and cervical cancer tissues and thus considered correlated with cervical cancer development (Wu et al., 2018). He et al. (2020) detected the gene and protein expression levels of KIF20A and found that the up-regulation of KIF20A expression is positively associated with the proliferation and invasion of cervical cancer cells. KIF20A affects the proliferation, migration, and invasion of cervical cancer by participating in a new type of circ_0005576/miR-153-3p/KIF20A axis (Ma et al., 2019), and KIF20A expression is involved in HPV infection, clinical stage, tumor recurrence, pelvic lymph node metastasis, and poor outcomes in early-stage cervical squamous cell cancer (Zhang et al., 2016). Thus, aberrant KIF20A expression is a potential independent prognostic factor for cervical cancer and is correlated with other cancers; for instance, radiation resistance may be enhanced partly by inducing the expression of KIF20A (Xiu et al., 2018). Irradiation-induced IF20A overexpression is associated with liver fibrosis, and KIF20A is associated with advanced pancreatic cancer (Taniuchi et al., 2005). Therefore, KIF20A-targeted treatment may be of great significance and may serve as a potential method for cervical cancer or its radiotherapy.
ZNF683, as a transcription factor, disrupts transcriptional processes in tissue-resident memory T and natural killer T cells, supplying immune protection when the body is re-infected (Van Gisbergen et al., 2012; Mackay et al., 2016). Liu et al. (2018) conducted analysis to clarify the potential mechanism of biological behavior factors that cause cervical cancer. Their results elucidated that the ZNF683 mRNAs are linked to the GO term “adaptive immune response.” This result is consistent with ours. In our study, based on the gene list of the top 30 scores in the random forest model, GO function enrichment was constructed using the R package for the prediction of pathways on cervical cancer radiosensitivity and pathways related to immune response. ZNF683 probably mediated the regulation of immune response in malignant cervical cancer. F10 gene may lead to the malignant transformation of the hydatidiform mole and occurrence of gynecological malignancies (Zhang and Pang, 2017). In patients with F10 overexpression, poorly differentiated cervical cancer tissues are higher in proportion compared with well-differentiated cancer tissues, and the negative correlation may be related to the occurrence and development of cervical cancer. KLRK1 encodes NKG2D, serving as an activating receptor expressed by NK cells and T cell subsets (Raulet, 2003). KLRK1 was thought to promote cervical cancer and CIN lesion susceptibility in 195 patients from southern Brazil.
Although not related to cervical cancer, most of the other parts of the genes we obtained through bioinformatics analysis were found to be related to radiosensitivity for other cancers. LCK was verified as a member of the Src family of non-receptor tyrosine kinases (Veillette et al., 1991). In Chinese hamster V79 cells, the expression of LCK gene is a key response to ionizing radiation-induced apoptosis and increases with radiation dose and incubation period. The Src-like tyrosine kinase p56/Lck is involved in the regulation of apoptosis caused by ionizing radiation. Similarly, in radioresistant B-lymphoma cells, p56/Lck promotes radiation-associated cell apoptosis (Waddick et al., 1993; Belka et al., 1999). The B-cell receptor-associated Src-like kinase Btk exhibits this function in B cells (Uckun et al., 1996).
ZAP70 is involved in T cell receptor signal transduction. The ability of T lymphocytes to convert signals to the nucleus is essential to the initiation and preservation of immune response, and the stimulation of T-cell receptors requires ligands to stimulate events in cells organized by ZAP70 and other proteins (Mozaffari et al., 2007). According to our results, ZAP70 plays an important role in radiotherapy-related immune response pathway. Similarly, a longitudinal study focused on immunologic system effects of adjuvant chemotherapy and radiotherapy in breast cancer showed that regulatory T cells and ZAP70 remain in their low grade forms after adjuvant chemoradiotherapy. In radiotherapy for colorectal cancer, through RNA extraction and following bioinformatic analysis, overexpressed ZAP70 was found to influence prognostic factors, such as tumor size and lymph node metastasis and differentiation; this finding showed that ZAP70 gene may act as a sensitivity biomarker for radiotherapy in colorectal cancer (Huang et al., 2011). Most of our information about PD-1 signaling was obtained through the exploration of activated T cells. In addition, the inhibition of ZAP70 impedes the activation of transcription factors in PD-1 signaling transduction during the stimulation of T cells (Wu et al., 2019).
As for CD247 and CD8B, related reports have been published. Mori et al. (2021) performed univariate experiments and analysis and found that reduction in PD-L1 mRNA/CD8B mRNA in tumor tissues may be a specific prognostic factor for salvage surgery for local recurrence after definitive concurrent radiotherapy and chemotherapy. CD247, also called CD3z-chain or TCR-Z, is a 16 kDa molecule participating in the TCR complex (Germain and Stefanová, 1999). Abnormal CD247 causes the irregular stimulation of T cells upon appointment with the TCR. CD247 down-regulation in TILs and peripheral blood lymphocytes is related to several types of cancers, including gastric carcinoma, head and neck cancer, and ovarian cancer, suggesting that CD247 is a therapeutic target for these cancer types (Ishigami et al., 2002; Ye et al., 2019). However, the relationships of the four genes with cervical cancer are rarely investigated despite that they may play vital role in cervical cancer. Additional tests are needed to verify the mechanism of these genes in the radiosensitivity of cervical cancer.
We employed Pearson correlation analysis to assess the meaningful relationships among the DEGs. Random forest was used in revealing the error rate of classifying samples. The results are satisfying and suitable for distinguishing cervical cancer-related genes from the samples. The results of the top 30 gene enrichment pathways have been verified by medical literature search to confirm their correlation with the sensitivity of cervical cancer radiotherapy, indicating that our results are precise. Notably, the medical literature search further confirmed that the expression levels of the 10 risk genes are associated with cervical cancer radiotherapy response or different types of tumors to a certain degree. Hence, the results suggested that our analytical application is of great value in identifying radiotherapy resistant molecular biomarkers for cervical cancer.
However, our current study has limitations and drawbacks. In this research, we verified the top 30 genes through literature searching because this method is highly comprehensive and can be easily conducted. In subsequent studies, we will perform accurate experimental techniques in pre-irradiation and mid-irradiation cervical cancer cell lines to validate the mRNA and protein expression levels of the 10 signature genes in cervical cancer. According to the results of BP analyses, we concluded that our candidate genes were mainly abundant in the immune responses of tumors, which is a new perspective for our further study of the mechanisms. We will preserve on valuable researches of cervical cancer radiotherapy resistance to explore the potential therapeutic targets for it.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE3578.
Author Contributions
YF and YZ project design. NY and ZW project supervision. YF, SL, JY, and JS data generation and analysis. YF and YZ revisions of the manuscript. All authors reviewed the manuscript.
Funding
This study was supported by the National Natural Science Foundation (U20A20339).
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.
Footnotes
References
Barker, H. E., Paget, J. T., Khan, A. A., and Harrington, K. J. (2015). The tumour microenvironment after radiotherapy: mechanisms of resistance and recurrence. Nat. Rev. Cancer 15, 409–425. doi: 10.1038/nrc3958
Baskar, R., Lee, K. A., Yeo, R., and Yeoh, K. W. (2012). Cancer and radiation therapy: current advances and future directions. Int. J. Med. Sci. 9, 193–199. doi: 10.7150/ijms.3635
Bazan, J. G., Luxton, G., Kozak, M. M., Anderson, E. M., Hancock, S. L., Kapp, D. S., et al. (2013). Impact of chemotherapy on normal tissue complication probability models of acute hematologic toxicity in patients receiving pelvic intensity modulated radiation therapy. Int. J. Radiat. Oncol. Biol. Phys. 87, 983–991. doi: 10.1016/j.ijrobp.2013.09.017
Belka, C., Marini, P., Lepple-Wienhues, A., Budach, W., Jekle, A., Los, M., et al. (1999). The tyrosine kinase lck is required for CD95-independent caspase-8 activation and apoptosis in response to ionizing radiation. Oncogene 18, 4983–4992. doi: 10.1038/sj.onc.1202878
Cheng, J., Lu, X., Wang, J., Zhang, H., Duan, P., and Li, C. (2017). Interactome analysis of gene expression profiles of cervical cancer reveals dysregulated mitotic gene clusters. Am. J. Transl. Res. 9, 3048–3059.
Cheng, L., Hu, Y., Sun, J., Zhou, M., and Jiang, Q. (2018). DincRNA: a comprehensive web-based bioinformatics toolkit for exploring disease associations and ncRNA function. Bioinformatics 34, 1953–1956. doi: 10.1093/bioinformatics/bty002
Cheng, L., Yang, H., Zhao, H., Pei, X., Shi, H., Sun, J., et al. (2019). MetSigDis: a manually curated resource for the metabolic signatures of diseases. Brief. Bioinform. 20, 203–209. doi: 10.1093/bib/bbx103
Choo, Y. W., Kang, M., Kim, H. Y., Han, J., Kang, S., Lee, J. R., et al. (2018). M1 macrophage-derived nanovesicles potentiate the anticancer efficacy of immune checkpoint inhibitors. ACS Nano 12, 8977–8993. doi: 10.1021/acsnano.8b02446
Chopra, S., Gupta, M., Mathew, A., Mahantshetty, U., Engineer, R., Lavanya, G., et al. (2018). Locally advanced cervical cancer: a study of 5-year outcomes. Indian J. Cancer 55, 45–49. doi: 10.4103/ijc.IJC_428_17
Chu, T. Y., Yang, J. T., Huang, T. H., and Liu, H. W. (2014). Crosstalk with cancer-associated fibroblasts increases the growth and radiation survival of cervical cancer cells. Radiat. Res. 181, 540–547. doi: 10.1667/rr13583.1
de Vos van Steenwijk, P. J., Ramwadhdoebe, T. H., Goedemans, R., Doorduijn, E. M., Van Ham, J. J., Gorter, A., et al. (2013). Tumor-infiltrating CD14-positive myeloid cells and CD8-positive T-cells prolong survival in patients with cervical carcinoma. Int. J. Cancer 133, 2884–2894. doi: 10.1002/ijc.28309
Deng, L., Wang, J., and Zhang, J. (2019a). Predicting gene ontology function of human MicroRNAs by integrating multiple networks. Front. Genet. 10:3. doi: 10.3389/fgene.2019.00003
Deng, L., Yang, W., and Liu, H. (2019b). PredPRBA: prediction of protein-RNA binding affinity using gradient boosted regression trees. Front. Genet. 10:637. doi: 10.3389/fgene.2019.00637
Doja, M. N., Kaur, I., and Ahmad, T. (2020). Current state of the art for survival prediction in cancer using data mining techniques. Curr. Bioinform. 15, 174–186. doi: 10.2174/1574893614666190902152142
Dueñas-González, A., Zarbá, J. J., Patel, F., Alcedo, J. C., Beslija, S., Casanova, L., et al. (2011). Phase III, open-label, randomized study comparing concurrent gemcitabine plus cisplatin and radiation followed by adjuvant gemcitabine and cisplatin versus concurrent cisplatin and radiation in patients with stage IIB to IVA carcinoma of the cervix. J. Clin. Oncol. 29, 1678–1685. doi: 10.1200/jco.2009.25.9663
Fridman, W. H., Zitvogel, L., Sautès-Fridman, C., and Kroemer, G. (2017). The immune contexture in cancer prognosis and treatment. Nat. Rev. Clin. Oncol. 14, 717–734. doi: 10.1038/nrclinonc.2017.101
Fulda, S., and Debatin, K. M. (2006). Extrinsic versus intrinsic apoptosis pathways in anticancer chemotherapy. Oncogene 25, 4798–4811. doi: 10.1038/sj.onc.1209608
Gauthier, J., Vincent, A. T., Charette, S. J., and Derome, N. (2019). A brief history of bioinformatics. Brief. Bioinform. 20, 1981–1996. doi: 10.1093/bib/bby063
Germain, R. N., and Stefanová, I. (1999). The dynamics of T cell receptor signaling: complex orchestration and the key roles of tempo and cooperation. Annu. Rev. Immunol. 17, 467–522. doi: 10.1146/annurev.immunol.17.1.467
Green, J. A., Kirwan, J. M., Tierney, J. F., Symonds, P., Fresco, L., Collingwood, M., et al. (2001). Survival and recurrence after concomitant chemotherapy and radiotherapy for cancer of the uterine cervix: a systematic review and meta-analysis. Lancet 358, 781–786. doi: 10.1016/s0140-6736(01)05965-7
Han, X., Kong, Q., Liu, C., Cheng, L., and Han, J. (2021). SubtypeDrug: a software package for prioritization of candidate cancer subtype-specific drugs. Bioinformatics. doi: 10.1093/bioinformatics/btab011 Online ahead of print,
He, Q., Meng, J., Liu, S., Zeng, Q., Zhu, Q., Wei, Z., et al. (2020). Long non-coding RNA UCA1 upregulates KIF20A expression to promote cell proliferation and invasion via sponging miR-204 in cervical cancer. Cell Cycle 19, 2486–2495. doi: 10.1080/15384101.2020.1807666
Huang, C., Chen, D., Zhu, H., Lv, S., Li, Q., and Li, G. (2019). LITAF enhances radiosensitivity of human glioma cells via the FoxO1 pathway. Cell Mol. Neurobiol. 39, 871–882. doi: 10.1007/s10571-019-00686-4
Huang, M. Y., Wang, J. Y., Chang, H. J., Kuo, C. W., Tok, T. S., and Lin, S. R. (2011). CDC25A, VAV1, TP73, BRCA1 and ZAP70 gene overexpression correlates with radiation response in colorectal cancer. Oncol. Rep. 25, 1297–1306. doi: 10.3892/or.2011.1193
Ishigami, S., Natsugoe, S., Tokuda, K., Nakajo, A., Higashi, H., Iwashige, H., et al. (2002). CD3-zetachain expression of intratumoral lymphocytes is closely related to survival in gastric carcinoma patients. Cancer 94, 1437–1442. doi: 10.1002/cncr.10346
Jarosz-Biej, M., Smolarczyk, R., Cichoń, T., and Kułach, N. (2019). Tumor microenvironment as a “Game Changer” in cancer radiotherapy. Int. J. Mol. Sci. 20:3212. doi: 10.3390/ijms20133212
Jemal, A., Bray, F., Center, M. M., Ferlay, J., Ward, E., and Forman, D. (2011). Global cancer statistics. CA Cancer J. Clin. 61, 69–90. doi: 10.3322/caac.20107
Jiang, L. P., He, C. Y., and Zhu, Z. T. (2017). Role of microRNA-21 in radiosensitivity in non-small cell lung cancer cells by targeting PDCD4 gene. Oncotarget 8, 23675–23689. doi: 10.18632/oncotarget.15644
Kaidar-Person, O., Bortnyak-Abdah, R., Amit, A., Berniger, A., Ben-Yosef, R., and Kuten, A. (2012). Current principles for radiotherapy in cervical cancer. Med. Oncol. 29, 2919–2922. doi: 10.1007/s12032-012-0170-7
Kanehisa, M., and Bork, P. (2003). Bioinformatics in the post-sequence era. Nat. Genet. 33(Suppl.), 305–310. doi: 10.1038/ng1109
Lippens, L., Van Bockstal, M., De Jaeghere, E. A., Tummers, P., Makar, A., De Geyter, S., et al. (2020). Immunologic impact of chemoradiation in cervical cancer and how immune cell infiltration could lead toward personalized treatment. Int. J. Cancer 147, 554–564. doi: 10.1002/ijc.32893
Liu, H., Zhang, W., Zou, B., Wang, J., Deng, Y., and Deng, L. (2020). DrugCombDB: a comprehensive database of drug combinations toward the discovery of combinatorial therapy. Nucleic Acids Res. 48, D871–D881.
Liu, J., Nie, S., Gao, M., Jiang, Y., Wan, Y., Ma, X., et al. (2019a). Identification of EPHX2 and RMI2 as two novel key genes in cervical squamous cell carcinoma by an integrated bioinformatic analysis. J. Cell Physiol. 234, 21260–21273. doi: 10.1002/jcp.28731
Liu, J., Zhou, S., Li, S., Jiang, Y., Wan, Y., Ma, X., et al. (2019b). Eleven genes associated with progression and prognosis of endometrial cancer (EC) identified by comprehensive bioinformatics analysis. Cancer Cell Int. 19:136. doi: 10.1186/s12935-019-0859-1
Liu, Q., Lu, C., Dai, W., Li, K., Xu, J., Huang, Y., et al. (2018). Association of biobehavioral factors with non-coding RNAs in cervical cancer. Biosci. Trends 12, 24–31. doi: 10.5582/bst.2017.01325
Lv, Z. B., Zhang, J., Ding, H., and Zou, Q. (2020). RF-PseU: a random forest predictor for RNA pseudouridine sites. Front. Bioeng. Biotechnol. 8:10. doi: 10.3389/fbioe.2020.00134
Ma, H., Tian, T., Liu, X., Xia, M., Chen, C., Mai, L., et al. (2019). Upregulated circ_0005576 facilitates cervical cancer progression via the miR-153/KIF20A axis. Biomed. Pharmacother. 118:109311. doi: 10.1016/j.biopha.2019.109311
Mackay, L. K., Minnich, M., Kragten, N. A., Liao, Y., Nota, B., Seillet, C., et al. (2016). Hobit and Blimp1 instruct a universal transcriptional program of tissue residency in lymphocytes. Science 352, 459–463. doi: 10.1126/science.aad2035
Moreira, D., Sampath, S., Won, H., White, S. V., Su, Y. L., Alcantara, M., et al. (2021). Myeloid cell-targeted STAT3 inhibition sensitizes head and neck cancers to radiotherapy and T cell-mediated immunity. J. Clin. Invest. 131:e137001. doi: 10.1172/jci137001
Mori, T., Kumagai, K., Nasu, K., Yoshizawa, T., Kuwano, K., Hamada, Y., et al. (2021). Clonal expansion of tumor-infiltrating T cells and analysis of the tumor microenvironment within esophageal squamous cell carcinoma relapsed after definitive chemoradiation therapy. Int. J. Mol. Sci. 22:1098. doi: 10.3390/ijms22031098
Mozaffari, F., Lindemalm, C., Choudhury, A., Granstam-Björneklett, H., Helander, I., Lekander, M., et al. (2007). NK-cell and T-cell functions in patients with breast cancer: effects of surgery and adjuvant chemo- and radiotherapy. Br. J. Cancer 97, 105–111. doi: 10.1038/sj.bjc.6603840
Nunn, M. E., Fan, J., Su, X., Levine, R. A., Lee, H. J., and Mcguire, M. K. (2012). Development of prognostic indicators using classification and regression trees for survival. Periodontology 2000 58, 134–142. doi: 10.1111/j.1600-0757.2011.00421.x
Oliver, G. R., Hart, S. N., and Klee, E. W. (2015). Bioinformatics for clinical next generation sequencing. Clin. Chem. 61, 124–135. doi: 10.1373/clinchem.2014.224360
Peng, J., Wang, H., Lu, J., Hui, W., Wang, Y., and Shang, X. (2017). Identifying term relations cross different gene ontology categories. BMC Bioinform. 18:573. doi: 10.1186/s12859-017-1959-3
Prevo, R., Tiwana, G. S., Maughan, T. S., Buffa, F. M., Mckenna, W. G., and Higgins, G. S. (2017). Depletion of signal recognition particle 72kDa increases radiosensitivity. Cancer Biol. Ther. 18, 425–432. doi: 10.1080/15384047.2017.1323587
Raulet, D. H. (2003). Roles of the NKG2D immunoreceptor and its ligands. Nat. Rev. Immunol. 3, 781–790. doi: 10.1038/nri1199
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
Rose, P. G., Bundy, B. N., Watkins, E. B., Thigpen, J. T., Deppe, G., Maiman, M. A., et al. (1999). Concurrent cisplatin-based radiotherapy and chemotherapy for locally advanced cervical cancer. N. Engl. J. Med. 340, 1144–1153. doi: 10.1056/nejm199904153401502
Ru, X. Q., Li, L. H., and Zou, Q. (2019). Incorporating distance-based top-n-gram and random forest to identify electron transport proteins. J. Proteome Res. 18, 2931–2939. doi: 10.1021/acs.jproteome.9b00250
Sahu, R., and Pattanayak, S. P. (2020). Strategic developments & future perspective on gene therapy for breast cancer: role of mTOR and Brk/PTK6 as molecular targets. Curr. Gene Ther. 20, 237–258. doi: 10.2174/1566523220999200731002408
Sapir-Pichhadze, R., and Kaplan, B. (2020). Seeing the forest for the trees: random forest models for predicting survival in kidney transplant recipients. Transplantation 104, 905–906. doi: 10.1097/tp.0000000000002923
Schnellhardt, S., Erber, R., Büttner-Herold, M., Rosahl, M. C., Ott, O. J., Strnad, V., et al. (2020). Accelerated partial breast irradiation: macrophage polarisation shift classification identifies high-risk tumours in early hormone receptor-positive breast cancer. Cancers 12:446. doi: 10.3390/cancers12020446
Siegel, R. L., Miller, K. D., and Jemal, A. (2020). Cancer statistics, 2020. CA Cancer J. Clin. 70, 7–30. doi: 10.3322/caac.21590
Stapleton, S., Jaffray, D., and Milosevic, M. (2017). Radiation effects on the tumor microenvironment: implications for nanomedicine delivery. Adv. Drug Deliv. Rev. 109, 119–130. doi: 10.1016/j.addr.2016.05.021
Steffens, M., Lamina, C., Illig, T., Bettecken, T., Vogler, R., Entz, P., et al. (2006). SNP-based analysis of genetic substructure in the German population. Hum. Hered. 62, 20–29. doi: 10.1159/000095850
Tang, Q., Nie, F., Kang, J., and Chen, W. (2021). mRNALocater: Enhance the prediction accuracy of eukaryotic mRNA subcellular localization by using model fusion strategy. Mol. Ther. doi: 10.1016/j.ymthe.2021.04.004 Online ahead of print,
Tang, W., Wan, S., Yang, Z., Teschendorff, A. E., and Zou, Q. (2018). Tumor origin detection with tissue-specific miRNA and DNA methylation markers. Bioinformatics 34, 398–406. doi: 10.1093/bioinformatics/btx622
Taniuchi, K., Nakagawa, H., Nakamura, T., Eguchi, H., Ohigashi, H., Ishikawa, O., et al. (2005). Down-regulation of RAB6KIFL/KIF20A, a kinesin involved with membrane trafficking of discs large homologue 5, can attenuate growth of pancreatic cancer cell. Cancer Res. 65, 105–112.
Ubels, J., Schaefers, T., Punt, C., Guchelaar, H. J., and De Ridder, J. (2020). RAINFOREST: a random forest approach to predict treatment benefit in data from (failed) clinical drug trials. Bioinformatics 36, i601–i609. doi: 10.1093/bioinformatics/btaa799
Uckun, F. M., Waddick, K. G., Mahajan, S., Jun, X., Takata, M., Bolen, J., et al. (1996). BTK as a mediator of radiation-induced apoptosis in DT-40 lymphoma B cells. Science 273, 1096–1100. doi: 10.1126/science.273.5278.1096
Van Gisbergen, K. P., Kragten, N. A., Hertoghs, K. M., Wensveen, F. M., Jonjic, S., Hamann, J., et al. (2012). Mouse hobit is a homolog of the transcriptional repressor Blimp-1 that regulates NKT cell effector differentiation. Nat. Immunol. 13, 864–871. doi: 10.1038/ni.2393
Veillette, A., Abraham, N., Caron, L., and Davidson, D. (1991). The lymphocyte-specific tyrosine protein kinase p56lck. Semin. Immunol. 3, 143–152.
Verhey, K. J., and Hammond, J. W. (2009). Traffic control: regulation of kinesin motors. Nat. Rev. Mol. Cell Biol. 10, 765–777. doi: 10.1038/nrm2782
Waddick, K. G., Chae, H. P., Tuel-Ahlgren, L., Jarvis, L. J., Dibirdik, I., Myers, D. E., et al. (1993). Engagement of the CD19 receptor on human B-lineage leukemia cells activates LCK tyrosine kinase and facilitates radiation-induced apoptosis. Radiat. Res. 136, 313–319. doi: 10.2307/3578542
Wadlow, R., and Ramaswamy, S. (2005). DNA microarrays in clinical cancer research. Curr. Mol. Med. 5, 111–120. doi: 10.2174/1566524053152834
Wang, H., and Zhou, L. (2017). Random survival forest with space extensions for censored data. Artif. Intell. Med. 79, 52–61. doi: 10.1016/j.artmed.2017.06.005
Wang, H., Liu, C., and Deng, L. (2018). Enhanced prediction of hot spots at protein-protein interfaces using extreme gradient boosting. Sci. Rep. 8:14285.
Wu, D. M., Shi, J., Liu, T., Deng, S. H., Han, R., and Xu, Y. (2018). Integrated analysis reveals down-regulation of SPARCL1 is correlated with cervical cancer development and progression. Cancer Biomark. 21, 355–365. doi: 10.3233/cbm-170501
Wu, X., Gu, Z., Chen, Y., Chen, B., Chen, W., Weng, L., et al. (2019). Application of PD-1 Blockade in Cancer Immunotherapy. Comput. Struct. Biotechnol. J. 17, 661–674. doi: 10.1016/j.csbj.2019.03.006
Xing, L., Mi, W., Zhang, Y., Tian, S., Zhang, Y., Qi, R., et al. (2020). The identification of six risk genes for ovarian cancer platinum response based on global network algorithm and verification analysis. J. Cell Mol. Med. 24, 9839–9852. doi: 10.1111/jcmm.15567
Xiu, G., Sui, X., Wang, Y., and Zhang, Z. (2018). FOXM1 regulates radiosensitivity of lung cancer cell partly by upregulating KIF20A. Eur. J. Pharmacol. 833, 79–85. doi: 10.1016/j.ejphar.2018.04.021
Yang, H., Luo, Y., Ren, X., Wu, M., He, X., Peng, B., et al. (2021). Risk prediction of diabetes: big data mining with fusion of multifarious physical examination indicators. Inf. Fusion 75, 140–149. doi: 10.1016/j.inffus.2021.02.015
Yang, L., Wu, H., Jin, X., Zheng, P., Hu, S., Xu, X., et al. (2020). Study of cardiovascular disease prediction model based on random forest in eastern China. Sci. Rep. 10:5245. doi: 10.1038/s41598-020-62133-5
Yang, W., and Deng, L. (2020). PreDBA: a heterogeneous ensemble approach for predicting protein-DNA binding affinity. Sci. Rep. 10:1278.
Ye, W., Zhou, Y., Xu, B., Zhu, D., Rui, X., Xu, M., et al. (2019). CD247 expression is associated with differentiation and classification in ovarian cancer. Medicine 98:e18407. doi: 10.1097/md.0000000000018407
Zhang, H., Yue, J., Jiang, Z., Zhou, R., Xie, R., Xu, Y., et al. (2017). CAF-secreted CXCL1 conferred radioresistance by regulating DNA damage response in a ROS-dependent manner in esophageal squamous cell carcinoma. Cell Death Dis. 8:e2790. doi: 10.1038/cddis.2017.180
Zhang, L., Jiang, X., Li, Y., Fan, Q., Li, H., Jin, L., et al. (2020a). Clinical correlation of Wnt2 and COL8A1 with colon adenocarcinoma prognosis. Front. Oncol. 10:1504. doi: 10.3389/fonc.2020.01504
Zhang, M. Z., and Pang, Z. J. (2017). [F10 expressions in cervical cancer tissues]. Nan Fang Yi Ke Da Xue Xue Bao. 37, 792–796. doi: 10.3969/j.issn.1673-4254.2017.06.13
Zhang, S., He, X., Zhang, R., and Deng, W. (2020b). LncR2metasta: a manually curated database for experimentally supported lncRNAs during various cancer metastatic events. Brief. Bioinform. 20:bbaa178.
Zhang, W., He, W., Shi, Y., Gu, H., Li, M., Liu, Z., et al. (2016). High expression of KIF20A is associated with poor overall survival and tumor progression in early-stage cervical squamous cell carcinoma. PLoS One 11:e0167449. doi: 10.1371/journal.pone.0167449
Zhang, X., Feng, H., Li, Z., Li, D., Liu, S., Huang, H., et al. (2018). Application of weighted gene co-expression network analysis to identify key modules and hub genes in oral squamous cell carcinoma tumorigenesis. Oncol. Targets Ther. 11, 6001–6021. doi: 10.2147/ott.S171791
Zhang, Y. (2020). Artificial intelligence for bioinformatics and biomedicine. Curr. Bioinform. 15, 801–802. doi: 10.2174/157489361508201221092330
Zhang, Z. M., Tan, J. X., Wang, F., Dao, F. Y., Zhang, Z. Y., and Lin, H. (2020c). Early diagnosis of hepatocellular carcinoma using machine learning method. Front. Bioeng. Biotechnol. 8:254. doi: 10.3389/fbioe.2020.00254
Zhang, Z. M., Wang, J. S., Zulfiqar, H., Lv, H., Dao, F. Y., and Lin, H. (2020d). Early diagnosis of pancreatic ductal adenocarcinoma by combining relative expression orderings with machine-learning method. Front. Cell Dev. Biol. 8:582864. doi: 10.3389/fcell.2020.582864
Zhao, T., Hu, Y., Peng, J., and Cheng, L. (2020). DeepLGP: a novel deep learning method for prioritizing lncRNA target genes. Bioinformatics 36, 4466–4472. doi: 10.1093/bioinformatics/btaa428
Keywords: biomarkers, cervical cancer, radiotherapy resistance, bioinformatic, RNA sequencing data
Citation: Feng Y, Wang Z, Yang N, Liu S, Yan J, Song J, Yang S and Zhang Y (2021) Identification of Biomarkers for Cervical Cancer Radiotherapy Resistance Based on RNA Sequencing Data. Front. Cell Dev. Biol. 9:724172. doi: 10.3389/fcell.2021.724172
Received: 12 June 2021; Accepted: 14 July 2021;
Published: 03 August 2021.
Edited by:
Lei Deng, Central South University, ChinaCopyright © 2021 Feng, Wang, Yang, Liu, Yan, Song, Yang and Zhang. 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: Yunyan Zhang, zhangyunyan@hrbmu.edu.cn; Shanshan Yang, yangshanshan@hrbmu.edu.cn