Systematic Review ARTICLE
Prognostic and Clinic Pathological Value of Cx43 Expression in Glioma: A Meta-Analysis
- 1Department of Neurosurgery, Qilu Hospital, Shandong University, Jinan, China
- 2Brain Science Research Institute, Shandong University, Jinan, China
Gap junctional intercellular communication (GJIC) composed of connexin proteins is considered vital to cancer onset and progression since 50 years ago based on Lowenstein and Kano's works, however altered expression of connexins is still a lesser known “hallmark” of cancer. Although many studies support the hypothesis that connexins are tumor suppressors, recent evidence indicates that, in some tumor types including glioma, they may play contradictory role in some specific stages of tumor progression. We thus conduct a meta-analysis to evaluate the prognostic role of Cx43 in glioma for the unanswered questions that whether Cx43 is a beneficial or insalubrity factor for glioma. Eight studies with 1,706 patients were included for meta-analysis. The results showed that Cx43 expression was a clearly negative factor with tumor grades (I2 = 34%, P < 0.001) and beneficial for OS (n = 3, HR 2.62, 95%CI 1.47–4.68; P = 0.001). Subgroup analysis also found that Cx43 had different expression in Asian young patients vs. other groups. In conclusion, this article summarize the prognostic value of Cx43 and offer a clinical evidence for the notion that Cx43 is generally a tumor suppressor and beneficial for the patients' survival time.
Brain tumors account for 1.9% of all new cancer cases and 2.3% of cancer related deaths globally. Among them, glioblastoma (GBM) is the most common primary malignant brain tumor (more than 45% of all malignant brain tumors) (1). However, the 5-year survival individuals with traditional treatment is still <5% since past decades, ranking it the 6th most lethal of all types among tumors (2). Standard procedure of care for newly diagnosed GBM includes maximum surgery resection, followed by ionizing radiation and chemotherapy with temozolomide (TMZ). Evidences showed that surgery resection is still the most effective way to cure GBM and excision extension is closely tied to patients' overall survival, so the importance of early-diagnosing GBM and exact-safe resection appeal to more and more focus (3). Therefore, it is necessary to identify molecular markers to help to distinguish surgical range and give a predicted clinical outcomes of high risk patients.
Cellular communication is vital in numerous processes critical for tumor biological homeostasis, for cell survival, proliferation, differentiation, and invasion. And in turn, it facilitated disease through a passive intercellular transmit of small molecules, second messengers, ions, microRNAs, and electrical signals (4). The structure of Connexins are tetra-span integral plasma membrane, consisting of two extracellular loops and one intracellular loop along with cytoplasmic amino (NH2) and carboxyl-terminal (CT) tail domains. Up to now, numbers of Cx molecules have been explored in the CNS: Cx30, Cx32, Cx36, and Cx43 in neurons; Cx30, Cx40, Cx43, and Cx45 in astrocytes; Cx32, Cx36, and Cx43 in microglia; and Cx26, Cx32, Cx29, Cx36, and Cx47 in oligodendrocyte (5, 6).
Although the big family consisting of Connexins, Cx43 is one of the most important subtypes which highly expressed in GBM tissues in some patients and appeared to be a marker distinguishing glioma from other types of brain cancers such as oligodendrocytic (7, 8). Cx43 has also been reported to possess the tumor-suppressor-like activities, to cause cancer cell dedifferentiation (mesenchymal to epithelial transition), and to inhibit metastasis thus affecting the tumor progression (6). Recently, few cases reported that Cx43 expression was inversely correlated to tumor grade as a result Cx43 was considered to “normalize” the phenotype of rat and human glioma cells (4, 9, 10). However, these articles failed to find any association between Cx43 and OS in glioma in their study and no meta-analysis reported that. In order to give a convincing evidence for this issue, we thus gathered the most recent studies and performed a meta-analysis to pool the results from eligible studies to present a quantitative calculation.
Literature Search Strategy
A systematic review and meta-analysis was conducted under the recommends of PRISMA guidelines. Data extraction was systematically performed in PubMed, Web of Science and EMBASE/Medlin, with different combinations of the following key words [“Glioma,” “brain tumor” OR “glioblastoma”] AND [“Connexin43,” “gap junction channels 43” OR “Cx43”] and no language restrictions. Terms [“Glioma,” “brain tumor” OR “glioblastoma”] AND [“Connexin,” “gap junction channels” OR “Cx”] were also searched and after eliminating papers focusing on other Cx proteins, the included results were same as the combination of glioma and Cx43. This search was finally renewed till May 1, 2019. References consists of studies, clinical cases and review articles were also included to find additional patients data.
Our purpose of this study was to identify whether Cx43 expression decreased the risk of glioma. Therefore, the criteria for enrolled study was as follows: (1) study should concerned the relationship between Cx43 and clinical outcomes among patients with gliomas. (2) Gliomas should be diagnosed by standard criteria, histopathologic analysis. (3) Cx43 expression was examined in glioma tissue obtained from glioma patients not experimental cells. (4) Cx43 expression was examined by immunohistochemistry (IHC) methods, PCR or tissue microarray. (5) The data should provide detailed information to calculate overall survival (OS) or disease free survival (DFS). Studies that could not meet any one of the above inclusion criteria were excluded. Studies which did not meet the above requirements were excluded and animal studies, letters were also precluded.
Two authors independently evaluated all studies and selected eligible trials, and any discrepancies between the authors were resolved by discussion and consensus. Among the studies, data was extracted included following information: country of the population enrolled, first author's name, number of included patients, and year of publication, sample size, patients' age, histological grade, diagnosed technique, and positive standard, antibody used, overall survival rates and survival outcomes. Kaplan-Meier (K-M) curves of these studies were analyzed by Engauge Digitizer (http://sourceforge.net/projects/digitizer/) which was utilized to calculate the OR and DFS. When univariate analysis and multivariate analysis were all provided in studies, we prefer the former results to be analyzed in the study. If the extract information could not be retained directly from included studies, we would send an email to the correspond author to get the original relevant data, otherwise the item would be marked as “Not Documented (ND).”
The Revman (version 5.3) was used to conduct all the calculations during the whole process of meta-analysis. Odds risks (OR) with 95% confidence intervals (CIs) were calculated and conducted to evaluate the correlation between Cx43 and clinical pathological outcomes in forest plots. I2 test and Q test were used to estimate heterogeneity in these studies. If the I2 test (>50%) or the Q test (P < 0.05) were abnormal, which indicated a significant heterogeneity between the selected studies, random-effect model would be introduced to assess the results, and otherwise the fixed-effect model would be used. Sub-group analyses based on research techniques (IHC or PCR), ethnicity and sample size were conducted. Sensitivity analysis was conducted to evaluate heterogeneity and stability of enrolled data. Potential publication bias were assessed by the funnel plots and Egger's tests. Also the effect of Cx43 diagnostic sensitivity and specificity were also presented by forest plot and SROC curve.
Study Selection and Characteristics
The flow of study selection has been presented in Figure 1. Based on an extensive combination of keywords search and screened a total of 173 papers by article title as well as abstract, we picked up 8 (9, 11–17) published studies that fulfilled all inclusion criteria which required intact data, strict experiment design and minor publication bias in the present meta-analysis. The enrolled studies were well-controlled and accorded with selection criterions. Based on the expression level of Cx43, enrolled patients of all the studies were divided into different subgroups: the high Cx43 level patients confirming to every positive standard in papers were classified into the Cx43 High subgroup, and patients with low Cx43 levels were attributed into the Cx43 Low subgroup. Overall, 8 studies constituting 790 Cx43 High patients and 916 Cx43 Low patients were evaluated with tumor grades, ethnicity, research technique and overall survival (OS).
Study Characteristics and Quality Assessment
The enrolled studies and clinical characteristics of included articles are presented in Table 1. Eight studies were conducted within western countries, and two within Asia. Five studies including more than 100 patients while the other three studies had relatively smaller patient's numbers. Two studies examined Cx43 expression by RT-PCR and six studies used IHC methods. Three articles evaluated cancer survival and recurrence. The publication time of all papers ranged from 2003 to 2016. The number size of enrolled group ranged from 32 to 572, and the positive rates of Cx43 expression varied from 55.6 to 89.2%. To examine the quality of included studies, Newcastle-Ottawa Quality Assessment scores (NOS) were introduced and the data ranged from six to nine (detail listed in Table 2), which manifested that the quality of enrolled studies was high. Exacted clinical data could be browsed in Tables 1, 2.
Association Between Cx43 in Glioma and 3-Year Overall Survival
Methods described previously were introduced to evaluate the results, the association between Cx43 and tumor grades were estimated by forest plots. Cx43 showed clearly negative with tumor grades in glioma patients (I2 = 34%, P < 0.001) (Figure 2). Three studies enrolled 1,202 patients were examined for the relationship between Cx43 and 3-year OS. And data (Figure 3A) showed that less Cx43 level was related with poor prognosis of glioma patients (HR 2.62, 95%CI 1.47–4.68; P = 0.001). Otherwise, this meta-analysis indicated that Cx43 level was highly related with a higher OS rates. The P value and I2 of this results showed they were eligible. There was no notable heterogeneities among these studies likely existed [Heterogeneity: Chi2 = 10.53, df = 7 (P = 0.16); I2 = 34%] (Figure 3B), so the resulting was reliable and then random-effect model was introduced to assess the possible publication bias of clinic pathological features.
Figure 2. Forest plot of studies evaluating the relationship between Cx43 expression and tumor grades.
Figure 3. (A) Forest plot of the association between Cx43 and OS in glioma patients. (B) Funnel plot for publication bias test of Cx43 related studies.
Association Between Cx43 and Clinic Pathological Features
To further confirm the potential value of Cx43 in clinical practice, the relationship between Cx43 and patients' clinic pathological characters including tumor grades, gender and average age were explored precisely. As seen in Table 3 and Figure 4, forest plots of 8 eligible studies showed the downregulated Cx43 was associated with tumor grades (P < 0.001). No difference was found between Cx43 expression in gender group (P = 0.86) (Figure 4A), but in age group, it effected the positive rate of Cx43 (P = 0.002). The difference may come from the different morbidity in different age and also the criteria for young (<60) and old (>60) do have influence on the bias (Heterogeneity: Chi2 = 22.49, P = 0.002; I2 = 69%) (Figure 4B).
Figure 4. Association between Cx43 and clinical pathological features in glioma patients. (A) Gender. (B) Age.
Sensitivity and Subgroup Analyses
Subgroup analysis was performed to explore the potential sources of heterogeneity, which were divided by research techniques, ethnicity and sample sizes. Table 4 and Figure 5 present the results of subgroup analyzed elucidate the correlations between Cx43 level and tumor grades. These results also indicated that low Cx43 expression was related to a significantly poorer results as compared to high expression. While sample size (P = 0.21) and research technique (P = 0.20) (IHC vs. PCR) did not obviously effect the prognosis rate of Cx43, but there might be a difference between Asian group and Western country group (Chi2 = 4.31, P = 0.04, I2 = 76.8%) (Figure 5B).
Figure 5. Subgroup analysis of Cx43 with ethnicity, sample size and research technique: (A) sample size; (B) region; (C) research technique.
Potential publication bias was evaluated by using Begg's funnel plot and Egger's test. The results were presented in Figure 5. The results indicated that there was no significant publication bias for all studies, and no evidence of significant publication bias was found in this paper.
As seen in Figure 6 and Table 4, group size didn't I effect the correlation between Cx43 level and tumor grades (Test for subgroup differences: Chi2 = 1.57, P = 0.21, I2 = 36.5%). However, the heterogeneity of ethnicity existed in these subgroups (Heterogeneity: Chi2 = 29.29, P = 0.0001; I2 = 76%). Additionally, research technique didn't influence the conclusion. And there were little heterogeneity among these subgroups (IHC: I2 = 45%, P = 0.10; RT-PCR: I2 = 0%, P = 0.80). While in both subgroups divided by ethnicity, Cx43 expression was correlated to patients' ethnicity (Test for subgroup differences: Chi2 = 4.31, P = 0.04, I2 = 76.8%). Heterogeneity in these subgroups was obvious, so further studies were still needed (Heterogeneity: Chi2 = 29.29, P = 0.0001; I2 = 76%). With these results, the heterogeneity of tumor grades mainly resulted from the sample size and ethnicity.
Figure 6. Funnel plot for publication bias test of Cx43 related studies. (A) Sample size. (B) Region. (C) Research technique.
Cx43 Diagnostic Sensitivity and Specificity
Sensitivity and specificity analyses were conducted to evaluate the diagnostic effect of Cx43 in glioma patients'. As seen in Figure 7, forest plot and SROC curve of 8 eligible studies showed the Cx43 is highly related to low grades glioma patients and shows good potential of sensitivity and specificity to define the tumor grades.
Figure 7. Diagnosis effect of Cx43 in glioma patients. (A) The sensitivity and specificity of Cx43 in glioma patients. (B) The diagnostic curve of Cx43 among glioma patients.
A huge number of studies have confirmed the relationship between cancer and connexins which was hypothesized more than 40 years ago (8, 18–20). However, the pivotal role of Cx43 in influencing tumor progression used to be overlooked because downregulation of Cx43-mediated intercellular communication is normally leading to increased malignancy in tumor cells (8, 19).As a result, Cx43 is usually considered as a tumor suppressor by its effect on reducing glioma proliferation (21, 22), anti-metastasis (23), and pro-apoptosis (24) and anti-inflammation (25). And our study is corresponded to these existing research that Cx43 expression is negatively associated with tumor grades.
In some studies, restoring Cx43 to glioma cells in vitro inhibits their tumorigenicity but this tumor suppressor effect could be glossed by its' promotion on invasion, adhesion and migration in vivo (1). Over-expression of Cx43 has been reported to enhance glioma migration in a channel-dependent manner, especially within the help of astrocytes (26, 27). Some experiments showed that Co-culture of glioma cells with astrocytes enhanced the invasiveness of the glioma cells and silencing Cx43 could extenuate this effect (25). Further study demonstrated this result and confirmed that junctions between glioma-glioma suppressed its invasiveness, while GJs of glioma–astrocyte and astrocyte–astrocyte promoted invasion (28). This seemed to be contrary to our findings that high Cx43 expression was beneficial to overall survive rate (OS). But it could be explained by the clinical experience that Cx43 usually highly expressed in low grades patients who could get gross-total resection (GTR) or even supra-total resection (SupTR) (29).
Due to nature of Cx43 that locates at the cell surface, Cx43 is associated with the brain endothelial barrier formation and the loss of Cx43 in glioma may lead to the infiltration of inflammation cells as well as more concentration of chemotherapy drugs (30). And further study of the correlationship between Cx43 and the popular glioma biomarkers is undergoing a proof study which would be published soon by our team. Besides Cx43 also the induction of inflammatory signaling, such as the treatment of anti-cancer agents and exposure to cytokines. One report demonstrated that Cx43 promoted glioma cell resistance to temozolomide by a channel-dependent (31), and another study find that Cx43 interacted Bax signal pathway by modulating mitochondrial apoptosis (32). Therefore, it will be beneficial to examine Cx43 expression in glioma before physician apply temozolomide or develop a drug to target Cx43. Cx43 peptidomimetics provide a potential method to selectively modulate the activity of connexin GJs as well as the numerous factors that correlated with tumorigenesis, recently study shown that inhibiting Cx43 restored TMZ sensitivity in TMZ-resistant/Cx43-high GBM cells including GSCs (33).
Besides, Cx43 is the most abundant Cx isoform in adult astrocytes especially in reactive astrocytes (34, 35). It was also demonstrated that reactive astrocytes such as M2 cell (abbreviation of pro-tumorigenic macrophage) is vital for innate immune systems in brain tumor microenvironment (36) and might also be responsible for immune therapy results. As a result, a high Cx43 expression is often seen in low grade gliomas at the early stage (15).
Cx43 channels are well known for the direct passage of small ions and metabolites such as Ca2+, ATP, glutamate, glucose and peptides through (37). However, emerging evidence has shown that Cx43 is also permeable to oligonucleotides as long as 24 nucleotides in length. Other results also demonstrated that glioma–astrocyte GJs were capable to transfer miRNAs (28). These findings adding to recent studies confirmed that cancer cells might “reprogram” the normal stromal cells by miRNAs through Cx43 junctions and normal cell might be hijacked by tumor cells to form a favorable environment for tumor as a result.
One unique character of Cx43 is the relationship with metabolism. Under glucose uptake situation, GJs allow the transition of glucose and other metabolic substrates throughout the network (38). Recently, studies find the inhibition of gap junctional communication or the decreasing of Cx43 expression resulted in an increase in the rate of glucose uptake (39) and aerobic glycolysis, which may be an explanation of poor prognosis of OS and high negative relevance ratio of Cx43 in high grade gliomas.
The role of Cx43 in tumor hemostasis is actually a very complex relationship. The mechanisms underlying these effects suggest a complex balance of variety proteins in the whole process (18). Our meta-analysis directly indicated that lower Cx43 as associated with poorer patient prognosis and thus revealed the potential value of Cx43 in diagnosis and prognosis of gliomas, however need more studies are still needed to elaborate the exact mechanism of Cx43 to bring it into clinical practice.
CZ wrote the main part of the manuscript including introduction results, methods and materials, and abstract. CL collected the related papers. AC and ZY assessed the quality of included papers. WL, SX, and XM revised the paper. All authors reviewed the manuscript.
This paper was funded by Scientific and technological Innovation Plan for Clinical Medicine of Jinan City (201805037) and Natural Science Foundation of China (Grant 81702469).
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.
We appreciated to the clinical cases provided by CL, Dr. Yang, SX, and WL.
9. Caltabiano R, Torrisi A, Condorelli D, Albanese V, Lanzafame S. High levels of connexin 43 mRNA in high grade astrocytomas. Study of 32 cases with in situ hybridization. Acta Histochem. (2010) 112:529–35. doi: 10.1016/j.acthis.2009.05.008
10. Huang R, Liu YG, Lin Y, Fan Y, Boynton A, Yang D, et al. Enhanced apoptosis under low serum conditions in human glioblastoma cells by connexin 43 (Cx43). Mol Carcinog. (2001) 32:128–38. doi: 10.1002/mc.1072
12. Gielen PR, Aftab Q, Ma N, Chen VC, Hong X, Lozinsky S, et al. Connexin43 confers Temozolomide resistance in human glioma cells by modulating the mitochondrial apoptosis pathway. Neuropharmacology. (2013) 75:539–48. doi: 10.1016/j.neuropharm.2013.05.002
13. Reszec J, Szkudlarek M, Hermanowicz A, Bernaczyk PS, Mariak Z, Chyczewski L. N-cadherin, beta-catenin and connexin 43 expression in astrocytic tumours of various grades. Histol Histopathol. (2015) 30:361–71. doi: 10.14670/HH-30.361
14. Murphy SF, Varghese RT, Lamouille S, Guo S, Pridham KJ, Kanabur P, et al. Connexin 43 inhibition sensitizes chemoresistant glioblastoma cells to temozolomide. Cancer Res. (2016) 76:139–49. doi: 10.1158/0008-5472.CAN-15-1286
17. Crespin S, Fromont G, Wager M, Levillain P, Cronier L, Monvoisin A, et al. Expression of a gap junction protein, connexin43, in a large panel of human gliomas: new insights. Cancer Med. (2016) 5:1742–52. doi: 10.1002/cam4.730
18. Tabernero A, Gangoso E, Jaraiz-Rodriguez M, Medina JM. The role of connexin43-Src interaction in astrocytomas: a molecular puzzle. Neuroscience. (2016) 323:183–94. doi: 10.1016/j.neuroscience.2015.02.029
19. Mesnil M, Crespin S, Avanzo JL, Zaidan-Dagli ML. Defective gap junctional intercellular communication in the carcinogenic process. Biochim Biophys Acta. (2005) 1719:125–45. doi: 10.1016/j.bbamem.2005.11.004
21. Saito T, Sato H, Virgona N, Hagiwara H, Kashiwagi K, Suzuki K, et al. Negative growth control of osteosarcoma cell by Bowman-Birk protease inhibitor from soybean; involvement of connexin 43. Cancer Lett. (2007) 253:249–57. doi: 10.1016/j.canlet.2007.01.021
23. Kou Y, Ji L, Wang H, Wang W, Zheng H, Zou J, et al. Connexin 43 upregulation by dioscin inhibits melanoma progression via suppressing malignancy and inducing M1 polarization. Int J Cancer. (2017) 141:1690–703. doi: 10.1002/ijc.30872
24. Uzu M, Sato H, Yamada R, Kashiba T, Shibata Y, Yamaura K, et al. Effect of enhanced expression of connexin 43 on sunitinib-induced cytotoxicity in mesothelioma cells. J Pharmacol Sci. (2015) 128:17–26. doi: 10.1016/j.jphs.2015.04.002
26. Oliveira R, Christov C, Guillamo JS, de Bouard S, Palfi S, Venance L, et al. Contribution of gap junctional communication between tumor cells and astroglia to the invasion of the brain parenchyma by human glioblastomas. BMC Cell Biol. (2005) 6:7. doi: 10.1186/1471-2121-6-7
27. Lin JH, Takano T, Cotrina ML, Arcuino G, Kang J, Liu S, et al. Connexin 43 enhances the adhesivity and mediates the invasion of malignant glioma cells. J Neurosci. (2002) 22:4302–11. doi: 10.1523/JNEUROSCI.22-11-04302.2002
29. Roh TH, Kang SG, Moon JH, Sung KS, Park HH, Kim SH, et al. Survival benefit of lobectomy over gross-total resection without lobectomy in cases of glioblastoma in the noneloquent area: a retrospective study. J Neurosurg. (2019). doi: 10.3171/2018.12.JNS182558. [Epub ahead of print].
30. Kaneko Y, Tachikawa M, Akaogi R, Fujimoto K, Ishibashi M, Uchida Y, et al. Contribution of pannexin 1 and connexin 43 hemichannels to extracellular calcium-dependent transport dynamics in human blood-brain barrier endothelial cells. J Pharmacol Exp Ther. (2015) 353:192–200. doi: 10.1124/jpet.114.220210
31. Kim SW, Choi HJ, Lee HJ, He J, Wu Q, Langley RR, et al. Role of the endothelin axis in astrocyte- and endothelial cell-mediated chemoprotection of cancer cells. Neuro Oncol. (2014) 16:1585–98. doi: 10.1093/neuonc/nou128
32. Uzu M, Sato H, Shimizu A, Shibata Y, Ueno K, Hisaka A. Connexin 43 enhances Bax activation via JNK activation in sunitinib-induced apoptosis in mesothelioma cells. J Pharmacol Sci. (2017) 134:101–7. doi: 10.1016/j.jphs.2017.05.005
33. Rhett JM, Calder BW, Fann SA, Bainbridge H, Gourdie RG, Yost MJ. Mechanism of action of the anti-inflammatory connexin43 mimetic peptide JM2. Am J Physiol Cell Physiol. (2017) 313:C314–26. doi: 10.1152/ajpcell.00229.2016
34. Szilvasy-Szabo A, Varga E, Beliczai Z, Lechan RM, Fekete C. Localization of connexin 43 gap junctions and hemichannels in tanycytes of adult mice. Brain Res. (2017) 1673:64–71. doi: 10.1016/j.brainres.2017.08.010
35. Kozoriz MG, Bechberger JF, Bechberger GR, Suen MW, Moreno AP, Maass K, et al. The connexin43 C-terminal region mediates neuroprotection during stroke. J Neuropathol Exp Neurol. (2010) 69:196–206. doi: 10.1097/NEN.0b013e3181cd44df
36. Kolar K, Freitas-Andrade M, Bechberger JF, Krishnan H, Goldberg GS, Naus CC, et al. Podoplanin: a marker for reactive gliosis in gliomas and brain injury. J Neuropathol Exp Neurol. (2015) 74:64–74. doi: 10.1097/NEN.0000000000000150
38. Wolf A, Agnihotri S, Micallef J, Mukherjee J, Sabha N, Cairns R, et al. Hexokinase 2 is a key mediator of aerobic glycolysis and promotes tumor growth in human glioblastoma multiforme. J Exp Med. (2011) 208:313–26. doi: 10.1084/jem.20101470
Keywords: glioma, Cx43, gap junctional intercellular communication, glioma survival time, meta-analysis
Citation: Zhang C, Liu C, Chen A, Yao Z, Li W, Xu S and Ma X (2019) Prognostic and Clinic Pathological Value of Cx43 Expression in Glioma: A Meta-Analysis. Front. Oncol. 9:1209. doi: 10.3389/fonc.2019.01209
Received: 29 July 2019; Accepted: 23 October 2019;
Published: 12 November 2019.
Edited by:Liam Chen, Johns Hopkins University, United States
Reviewed by:Justin Lathia, Case Western Reserve University, United States
Maria Caffo, University of Messina, Italy
Copyright © 2019 Zhang, Liu, Chen, Yao, Li, Xu and Ma. 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: Xiang-yu Ma, firstname.lastname@example.org