- 1Tianjin Key Laboratory of Ionic-Molecular Function of Cardiovascular Disease, Department of Cardiology, Tianjin Institute of Cardiology, Second Hospital of Tianjin Medical University, Tianjin, China
- 2Laboratory of Cardiac Electrophysiology, Second Department of Cardiology, “Evangelismos” General Hospital of Athens, Athens, Greece
- 3Divisions of Cardiology, Queen's University, Kingston General Hospital, Kingston, ON, Canada
- 4Department of Medicine and Therapeutics, Chinese University of Hong Kong, Hong Kong, China
- 5Li Ka Shing Institute of Health Sciences, Chinese University of Hong Kong, Hong Kong, China
- 6Department of Cardiology and Cardiovascular Research Institute, Renmin Hospital of Wuhan University, Wuhan, China
- 7Masonic Medical Research Laboratory, Utica, NY, United States
Fragmented QRS (fQRS) is an electrocardiographic marker related to ventricular fibrillation (VF) and sudden cardiac death (SCD) in various clinical settings. Current data regarding the prognostic significance of fQRS in Brugada syndrome (BrS) are contradictory. This meta-analysis aimed to evaluate the presence of fQRS as a risk stratification tool in BrS. Electronic databases (PubMed, EMBASE, and Cochrane Library) were searched until May 2016. Eight observational studies accumulating data on 1,637 BrS patients (mean age: 47 ± 11 years) were included in this meta-analysis. The mean follow-up duration ranged from 21 to 96 months. fQRS was found to be an independent predictor of future arrhythmic events in BrS (RR:3.88, 95% CI 2.26 to 6.65, p < 0.00001) with moderate heterogeneity (I2 = 54%, P = 0.03). When analyzing VF as independent end-point, the RR for VF was 3.61, and its 95% CI was 2.11 to 6.18, p < 0.00001. This meta-analysis showed that BrS patients with fQRS are at high risk for future arrhythmic events. The presence of fQRS warrants prospective evaluation as valid arrhythmogenic risk marker in BrS.
Introduction
Brugada syndrome (BrS) is a primary arrhythmic syndrome characterized by ST-segment elevation in the right precordial leads on the surface electrocardiogram (ECG) in the absence of overt structural heart disease (Brugada and Brugada, 1992; Mizusawa and Wilde, 2012; Antzelevitch and Patocskai, 2016; Tse et al., 2016). It is associated with a higher risk of ventricular fibrillation (VF) and sudden cardiac death (SCD) compared to the general population. BrS is a genetically heterogeneous ion channelopathy. Up to now, mutations in 19 genes have been identified in subjects with BrS (Mizusawa and Wilde, 2012). This disease typically manifests as cardiac arrest or syncope, occurring in the third and fourth decades of life (Mizusawa and Wilde, 2012; Antzelevitch and Patocskai, 2016). The majority of BrS patients are asymptomatic and diagnoses are often made incidentally. Risk stratification of BrS patients, and particularly those who are asymptomatic, remains challenging (Raju et al., 2011). Fragmented QRS (fQRS) is a simple non-invasive ECG depolarization marker used to identify individuals at high risk of ventricular arrhythmias and SCD in various clinical settings, including coronary artery disease, BrS, long QT syndrome, arrhythmogenic right ventricular cardiomyopathy, and cardiac sarcoidosis (Pietrasik and Zareba, 2012; Jain et al., 2014). It is defined as the presence of an additional R wave (R') or notching in the nadir of the S wave or the presence of >1 R' in two contiguous leads or the presence of more than two notches in the R or S waves in two consecutive leads in the presence of bundle branch block (Jain et al., 2014). Previous studies addressing the prognostic significance of fQRS in subjects with BrS have demonstrated conflicting results (Priori et al., 2012; Take et al., 2012; Tokioka et al., 2014; Calo et al., 2016). Therefore, we performed a comprehensive systematic review and meta-analysis of the current evidence regarding the prognostic significance of fQRS.
Methods
Two reviewers (L. M. and Q. Z.) independently and systematically performed a literature search on the PubMed, EMBASE, and Cochrane Library databases, to identify relevant studies. We searched the related studies published from May 2006 to 2016 using the following keywords: “Brugada” and “syndrome” or “Brugada syndrome” and “electrocardiography” or “ECG” or “implantable cardioverter defibrillator therapy.” In order to identify all potentially relevant studies, the titles, abstracts, and reference lists of all articles were carefully reviewed. The definition of fQRS included in this meta-analysis, according to previous studies what we synthesized (Das et al., 2007, 2010; Priori et al., 2012), was the presence of an abnormal fragmentation within QRS complex as ≥ two spikes in leads V1 to V3, or the presence of an additional wave or notching in the leads.
Studies were included if they met the following criteria: (a) the study design was a prospective or retrospective observational study; (b) patients with either a spontaneous or a drug-induced type 1 ECG BrS pattern; (c) fQRS was assessed with ECG or Holter monitoring; (d) the follow-up duration was sufficiently long that the arrhythmic events would be detected (the duration of follow-up ≥1 years); (e) endpoint events [appropriate implantable cardioverter-defibrillator therapy (ICD), VF, and SCD] were clearly defined; (f) studies with full-text; and (g) risk ratio (RR), hazard ratio (HR), odds ratio (OR), and the corresponding 95% confidence interval (CI) was reported, or could be calculated. We included both published and unpublished studies without language restriction. In the case of numerous reports by the same group of authors, only the study with the largest number of patients was included. The potential relevant studies were retrieved as full text and assessed for compliance with the inclusion criteria by two investigators. Any uncertainties or discrepancies were resolved through consensus after rechecking the source data and consultation with a third reviewer (T. L.).
Using a standard data extraction form, two blinded reviewers (L. M. and Q. Z.) independently performed data extraction to determine eligibility for inclusion. The extracted data elements of this meta-analysis consisted of: (a) publication details: surname of first author, publication year and location; (b) type of study: multicenter or single center study; (c) study design; (d) follow-up duration; (e) definition of fragmented QRS; (f) endpoint events (arrhythmic events were defined as ventricular fibrillation or flutter, SCD, and the combination of those two during the follow-up period; (g) the quality score; and (h) the characteristics of the population including sample size, gender, age, number of subjects with spontaneous or drug-induced BrS ECG pattern, number of subjects with ICD, number of subjects with family history of SCD or syncope, positive SCN5A gene mutation, detailed information in relation to programmed ventricular stimulation (PVS), positive number of inducible VF, and the presence of fQRS.
Quality assessment of these studies included in our meta-analysis was performed using the Newcastle–Ottawa Quality Assessment Scale (NOS; Gussak et al., 2000). The point score system evaluated the categories of study participant selection, comparability of the results, and quality of the outcomes. The following characteristics were assessed: (a) representativeness of the exposed cohort; (b) selection of the non-exposed cohort; (c) ascertainment of exposure; (d) demonstration that outcome of interest was not present at start of study; (e) comparability of cohorts on the basis of the design or analysis; (f) assessment of outcome; (g) follow-up period was sufficiently long for outcomes to be detected; and (h) adequacy of follow up of cohorts. This scale varied from zero to nine stars, which indicated that studies were graded as poor quality if they met <5 criteria, fair if they met 5 to 7 criteria, and good if they met >8 criteria. Studies with a score equal to or higher than six were considered to be concluded.
We extracted and analyzed all the multivariate adjusted RRs, ORs, or HRs with 95% confidence interval (CI) for each study. Pooled effect sizes were presented as the RR with 95% CI for each trial, using the random effects model. Since the related data were occasionally absent, raw data were used to calculate unadjusted risk estimates (Morita et al., 2008; Maury et al., 2013; Tokioka et al., 2014; Calo et al., 2016; Rivard et al., 2016). The HR values in multivariate Cox proportional hazards models in each primary study were directly considered as RR values (Greenland, 1987). To evaluate the heterogeneity across studies, I2 derived from the standard chi-square test, which described the percentage of the variability in effect estimates resulting from heterogeneity, was used. I2 > 50% is an indicator of significant statistical heterogeneity. In this case, the random-effects model using the inverse variance heterogeneity method was used. Subgroup analyses were performed based on the different primary endpoints and whether pooled effect sizes were adjusted in the patients with BrS. In addition, we performed the sensitivity analysis in a random predefined manner. Publication bias was evaluated using the funnel plot. Statistical significance was defined at P-values ≤ 0.05. All analyses were performed using Review Manager, version 5.0.12 (Revman; The Cochrane Collaboration, Oxford, United Kingdom).
Results
A flow diagram of the data search and study selection is presented in Figure 1. A total of 5,843 records were identified using our search criteria. 3,026 duplicate studies were discarded. After screening the titles and abstracts, 2,796 studies were excluded as they were guidelines, editorials, case reports, laboratory studies, animal studies, review articles, or irrelevant to the present study. Therefore, 21 potentially relevant studies were retrieved for detailed evaluation. Of these, 13 were further excluded from the analysis for the following reasons: in 10 studies, RRs, ORs, or HRs were not provided or couldn't be calculated, or 95% confidence intervals were not included; one study did not clearly define the type of the abnormal QRS complex; one failed to clearly define the endpoints, and one was an abstract without full-text.
A total of 1,637 patients (mean age: 47 ± 11 years) with BrS from eight studies (6 prospective cohort studies and 2 retrospective cohort studies) were included in our meta-analysis (Morita et al., 2008; Priori et al., 2012; Take et al., 2012; Maury et al., 2013; Apiyasawat et al., 2014; Tokioka et al., 2014; Calo et al., 2016; Rivard et al., 2016). The baseline characteristics of these studies are listed in Table 1. The mean follow-up duration ranged from 21 to 96 months. The definition of fQRS was similar in all studies. Table 2 shows the patients' characteristics included in each study.
Five out of eight studies demonstrated that the presence of fQRS predicted future arrhythmic events in patients with a spontaneous or drug-induced type 1 ECG pattern of BrS (Morita et al., 2008; Priori et al., 2012; Take et al., 2012; Apiyasawat et al., 2014; Tokioka et al., 2014). On the contrary, three studied failed to show any prognostic significance of fQRS (Maury et al., 2013; Calo et al., 2016; Rivard et al., 2016). As shown in Figure 2, the pooled meta-analysis of eight included studies demonstrated that the presence of fQRS is an independent predictor of future arrhythmic events in BrS (RR: 3.88, 95% CI 2.26 to 6.65, p < 0.00001) with moderate heterogeneity (I2 = 54%, P = 0.03).
Figure 2. Forest plot demonstrating the association between fQRS and the fatal arrhythmias events and sudden cardiac death in the patients with Brugada syndrome.
Subgroup analysis was performed based on different arrhythmic endpoints, When analyzing studies using VF/SCD as a combined endpoint, pooled analysis of four studies showed that fQRS was an independent predictor of future arrhythmic events (RR 3.36, 95% CI 2.03 to 5.56, p < 0.00001) with moderate heterogeneity (p = 0.005 and I2 = 76%; Morita et al., 2008; Maury et al., 2013; Tokioka et al., 2014; Calo et al., 2016). When analyzing VF and SCD as different end-points, the RRs were the following: VF from three studies (RR 3.61, 95% CI 2.11 to 6.18, p < 0.00001; Priori et al., 2012; Take et al., 2012; Apiyasawat et al., 2014) and SCD in one study (RR 2.36, 95% CI 0.69 to 8.07, p = 0.9502; Rivard et al., 2016).
We also performed subgroup analysis based on whether pooled effect sizes was adjusted from each study. Regardless of whether adjustment of potential confounding variables were made, the outcomes were similar: unadjusted RR (4.23, 95% CI 1.68 to 10.61, p = 0.002; Morita et al., 2008; Maury et al., 2013; Tokioka et al., 2014; Calo et al., 2016; Rivard et al., 2016) and adjusted HR (3.61, 95% CI 2.11 to 6.18, p < 0.00001; Priori et al., 2012; Take et al., 2012; Apiyasawat et al., 2014).
Sensitivity analysis was then performed to identify possible causes of the significant heterogeneity in our meta-analysis. After excluding the studies which the RR value was calculated according to the raw data, no significant heterogeneity was found in the remaining studies (I2 = 3%). As a consequence, a potential source of the heterogeneity observed in our meta-analysis may be related to the different measures used for risk assessment. As shown in the Figure 3, the results of the funnel plot suggested that little publication bias was present.
Discussion
The main findings of the present meta-analysis are that: (1) overall, the presence of fQRS is associated with a 3-fold increased risk of future arrhythmic events in patients with BrS, and (2) the presence of fQRS confers a 3- and 2-fold increase for VF and SCD events, respectively. We further performed the subgroup and sensitivity analyses to identify possible causes of the significant heterogeneity in the analysed studies. The results showed that the calculation method of risk ratios is a possible origin of heterogeneity observed in our meta-analysis.
Asymptomatic BrS patients display an annual event rate of arrhythmic events between 0.5 and 1% (Priori et al., 2012). Such arrhythmic events occur in about 50% of cases as VF without prior warning symptoms (Jain et al., 2014). Risk stratification of BrS patients is therefore of paramount importance. Several ECG markers have been proposed for risk stratification of BrS patients (Probst et al., 2010; Tse, 2016; Tse and Yan, 2016a,b; Tse et al., 2017), but the majority have not been tested in a prospective manner (Mizusawa and Wilde, 2012).
fQRS is a relatively new arrhythmogenic ECG marker validated in different clinical settings (Raju et al., 2011; Pietrasik and Zareba, 2012). In a recent meta-analysis of 5,009 patients with coronary artery disease and non-ischemic cardiomyopathy, the presence of fQRS was associated with a relative risk for SCD of 2.2 (Rosengarten et al., 2015). Data regarding the prognostic significance of fQRS in BrS are limited. Morita et al. have demonstrated that fQRS is more commonly seen in BrS patients with VF (85%) and syncope (50%) compared to asymptomatic ones (34%; Morita et al., 2008). The PRELUDE study confirmed these findings and showed that fQRS is an independent predictor of future arrhythmic events (HR: 4.94; Priori et al., 2012). Tokioka et al. recently demonstrated that the presence of QRS-fragmentation lead to a 5-fold increase in the incidence of arrhythmic events (VF, SCD; Tokioka et al., 2014). ICD is the golden therapeutic method for cardiac arrest. The study by Apiyasawat et al. had confirmed that fQRS was directly associated with appropriate ICD therapy (Apiyasawat et al., 2014). However, the choice of an appropriate ICD intervention as surrogate endpoint for SCD event may lead to overestimation of the potential risk.
The presence of fQRS on surface ECG has been attributed to scar-related inhomogeneous conduction of action potentials through the ventricles. The different morphologies of fQRS are caused by shifting of the QRS vector in and around the fibrotic areas during depolarization, depending on their extent and location in the ventricles (Jain et al., 2014). The pathophysiology of BrS is only partially resolved. So far, repolarization, depolarization, and the current mismatch hypotheses are thought to underlie the development of VF in BrS (Nishii et al., 2010). There is increasing evidence suggesting that mild structural abnormalities observed in right ventricular outflow tract provide the arrhythmia substrate in BrS. Thus, BrS has been recently associated with increased collagen and fibrosis, and reduced gap junction expression in the right ventricular outflow tract. It is therefore quite possible that fQRS in BrS reflects conduction delay in the right ventricle. Indeed, using an isolated canine right ventricular tissue model of BrS, Morita et al. demonstrated that activation delay in the epicardium reproduces fQRS in the transmural ECG (Morita et al., 2008).
Study Limitations
Our study has several potential limitations. Firstly, the small size of the study population may have an important impact in our findings. Second, a potential publication bias was identified in the funnel plot. Thirdly, some studies included in our studies are retrospective studies, which may more recall bias. Fourthly, a potential overlap exists in the studies by Take et al. (2012) and Tokioka et al. (2014). Fifthly, the control variables used for adjusted and unadjusted RR/HR in the different studies may influence the final results of our meta-analysis. However, the subgroup analysis based on whether risk ratio was adjusted from the potential confounding variables had similar final outcomes. Finally, the different measures used for risk assessment may have contributed to the heterogeneity in this meta-analysis.
Conclusion
In conclusion, the presence of fQRS is predictive of future arrhythmic events in patients with BrS. The presence of fQRS warrants prospective evaluation and validation as a risk marker for clinical use in BrS.
Author Contributions
LM and KL participated in study design, searched databases, extracted and assessed data, and drafted the manuscript. AB, GT, DH, and GL revised the manuscript, QS, NZ, and ZZ participated in extracted data, TL designed the study and revised the manuscript. All authors approved the final version of the manuscript.
Conflict of Interest Statement
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.
References
Antzelevitch, C., and Patocskai, B. (2016). Brugada syndrome: clinical, genetic, molecular, cellular, and ionic aspects. Curr. Probl. Cardiol. 41, 7–57. doi: 10.1016/j.cpcardiol.2015.06.002
Apiyasawat, S., Sahasthas, D., Ngarmukos, T., Chandanamattha, P., and Likittanasombat, K. (2014). Fragmented QRS as a predictor of appropriate implantable cardioverter-defibrillator therapy. Indian Pacing Electrophysiol. J. 14, 4–11. doi: 10.1016/S0972-6292(16)30710-0
Brugada, P., and Brugada, J. (1992). Right bundle branch block, persistent ST segment elevation and sudden cardiac death: a distinct clinical and electrocardiographic syndrome. A multicenter report. J. Am. Coll. Cardiol. 20, 1391–1396.
Calo, L., Giustetto, C., Martino, A., Sciarra, L., Cerrato, N., Marziali, M., et al. (2016). A new electrocardiographic marker of sudden death in brugada syndrome: the S-Wave in lead I. J. Am. Coll. Cardiol. 67, 1427–1440. doi: 10.1016/j.jacc.2016.01.024
Das, M. K., Maskoun, W., Shen, C., Michael, M. A., Suradi, H., Desai, M., et al. (2010). Fragmented QRS on twelve-lead electrocardiogram predicts arrhythmic events in patients with ischemic and nonischemic cardiomyopathy. Heart Rhythm 7, 74–80. doi: 10.1016/j.hrthm.2009.09.065
Das, M. K., Saha, C., El Masry, H., Peng, J., Dandamudi, G., Mahenthiran, J., et al. (2007). Fragmented QRS on a 12-lead ECG: a predictor of mortality and cardiac events in patients with coronary artery disease. Heart Rhythm 4, 1385–1392. doi: 10.1016/j.hrthm.2007.06.024
Greenland, S. (1987). Quantitative methods in the review of epidemiologic literature. Epidemiol. Rev. 9, 1–30. doi: 10.1093/oxfordjournals.epirev.a036298
Gussak, I., Chaitman, B. R., Kopecky, S. L., and Nerbonne, J. M. (2000). Rapid ventricular repolarization in rodents: electrocardiographic manifestations, molecular mechanisms, and clinical insights. J. Electrocardiol. 33, 159–170. doi: 10.1016/S0022-0736(00)80072-2
Jain, R., Singh, R., Yamini, S., and Das, M. K. (2014). Fragmented ECG as a risk marker in cardiovascular diseases. Curr. Cardiol. Rev. 10, 277–286. doi: 10.2174/1573403X10666140514103451
Maury, P., Rollin, A., Sacher, F., Gourraud, J. B., Raczka, F., Pasquie, J. L., et al. (2013). Prevalence and prognostic role of various conduction disturbances in patients with the Brugada syndrome. Am. J. Cardiol. 112, 1384–1389. doi: 10.1016/j.amjcard.2013.06.033
Mizusawa, Y., and Wilde, A. A. (2012). Brugada syndrome. Circ. Arrhythm. Electrophysiol. 5, 606–616. doi: 10.1161/circep.111.964577
Morita, H., Kusano, K. F., Miura, D., Nagase, S., Nakamura, K., Morita, S. T., et al. (2008). Fragmented QRS as a marker of conduction abnormality and a predictor of prognosis of Brugada syndrome. Circulation 118, 1697–1704. doi: 10.1161/CIRCULATIONAHA.108.770917
Nishii, N., Nagase, S., Morita, H., Kusano, K. F., Namba, T., Miura, D., et al. (2010). Abnormal restitution property of action potential duration and conduction delay in Brugada syndrome: both repolarization and depolarization abnormalities. Europace 12, 544–552. doi: 10.1093/europace/eup432
Pietrasik, G., and Zareba, W. (2012). QRS fragmentation: diagnostic and prognostic significance. Cardiol. J. 19, 114–121. doi: 10.5603/CJ.2012.0022
Priori, S. G., Gasparini, M., Napolitano, C., Della Bella, P., Ottonelli, A. G., Sassone, B., et al. (2012). Risk stratification in Brugada syndrome: results of the PRELUDE (PRogrammed ELectrical stimUlation preDictive valuE) registry. J. Am. Coll. Cardiol. 59, 37–45. doi: 10.1016/j.jacc.2011.08.064
Probst, V., Veltmann, C., Eckardt, L., Meregalli, P. G., Gaita, F., Tan, H. L., et al. (2010). Long-term prognosis of patients diagnosed with Brugada syndrome: results from the FINGER Brugada Syndrome Registry. Circulation 121, 635–643. doi: 10.1161/CIRCULATIONAHA.109.887026
Raju, H., Papadakis, M., Govindan, M., Bastiaenen, R., Chandra, N., O'Sullivan, A., et al. (2011). Low prevalence of risk markers in cases of sudden death due to Brugada syndrome relevance to risk stratification in Brugada syndrome. J. Am. Coll. Cardiol. 57, 2340–2345. doi: 10.1016/j.jacc.2010.11.067
Rivard, L., Roux, A., Nault, I., Champagne, J., Roux, J. F., Tadros, R., et al. (2016). Predictors of ventricular arrhythmias and sudden death in a Québec cohort with Brugada syndrome. Can. J. Cardiol. 32, 1355.e1–1355.e7. doi: 10.1016/j.cjca.2016.03.012
Rosengarten, J. A., Scott, P. A., and Morgan, J. M. (2015). Fragmented QRS for the prediction of sudden cardiac death: a meta-analysis. Europace 17, 969–977. doi: 10.1093/europace/euu279
Take, Y., Morita, H., Toh, N., Nishii, N., Nagase, S., Nakamura, K., et al. (2012). Identification of high-risk syncope related to ventricular fibrillation in patients with Brugada syndrome. Heart Rhythm 9, 752–759. doi: 10.1016/j.hrthm.2011.11.045
Tokioka, K., Kusano, K. F., Morita, H., Miura, D., Nishii, N., Nagase, S., et al. (2014). Electrocardiographic parameters and fatal arrhythmic events in patients with Brugada syndrome: combination of depolarization and repolarization abnormalities. J. Am. Coll. Cardiol. 63, 2131–2138. doi: 10.1016/j.jacc.2014.01.072
Tse, G. (2016). Novel conduction-repolarization indices for the stratification of arrhythmic risk. J. Geriatr. Cardiol. 9, 811–812. doi: 10.11909/j.issn.1671-5411.2016.09.008
Tse, G., Gong, M., Wong, W. T., Georgopoulos, S., Letsas, K. P., Vassiliou, V. S., et al. (2017). The Tpeak – Tend interval as an electrocardiographic risk marker of arrhythmic and mortality outcomes: a systematic review and meta-analysis. Heart Rhythm 14, 1131–1137. doi: 10.1016/j.hrthm.2017.05.031
Tse, G., Liu, T., Li, K. H., Laxton, V., Chan, Y. W., Keung, W., et al. (2016). Electrophysiological mechanisms of brugada syndrome: insights from pre-clinical and clinical studies. Front. Physiol. 7:467. doi: 10.3389/fphys.2016.00467
Tse, G., and Yan, B. P. (2016a). Novel arrhythmic risk markers incorporating QRS dispersion: QRSd x (Tpeak - Tend)/QRS and QRSd x (Tpeak - Tend)/(QT x QRS). Ann. Noninvasive Electrocardiol. 13, 811–812. doi: 10.1111/anec.12397
Keywords: fragmented QRS, Brugada syndrome, predictor, arrhythmic events, meta-analysis
Citation: Meng L, Letsas KP, Baranchuk A, Shao Q, Tse G, Zhang N, Zhang Z, Hu D, Li G and Liu T (2017) Meta-analysis of Fragmented QRS as an Electrocardiographic Predictor for Arrhythmic Events in Patients with Brugada Syndrome. Front. Physiol. 8:678. doi: 10.3389/fphys.2017.00678
Received: 23 March 2017; Accepted: 24 August 2017;
Published: 12 September 2017.
Edited by:
Ruben Coronel, University of Amsterdam, NetherlandsReviewed by:
Flavia Ravelli, University of Trento, ItalyAndrew F. James, University of Bristol, United Kingdom
Copyright © 2017 Meng, Letsas, Baranchuk, Shao, Tse, Zhang, Zhang, Hu, Li and Liu. 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) or licensor 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: Tong Liu, bGl1dG9uZ2RvY0AxMjYuY29t