Skip to main content

ORIGINAL RESEARCH article

Front. Cardiovasc. Med., 12 July 2021
Sec. Cardiac Rhythmology
Volume 8 - 2021 | https://doi.org/10.3389/fcvm.2021.622483

Novel Clue to Locate Conduction Gaps in the Pulmonary Vein Isolation Ablation Line

Hai-yang Xie Xiao-gang Guo Jian-du Yang Yan-qiao Chen Zhong-jing Cao Qi Sun Jian Ma*
  • State Key Laboratory of Cardiovascular Disease, Arrhythmia Center, Fuwai Hospital, National Center for Cardiovascular Diseases, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China

Background: Several methods have been reported for locating the conduction gap (CG) in the pulmonary vein isolation (PVI) ablation line. However, the value of the interval between far-field atrial potential (FFP) and pulmonary vein potential (PVP) remains unknown.

Methods: Consecutive patients with a CG during observation on the table after PVI were included. The PVP, FFP, and the CG location were evaluated to develop a novel algorithm to identify the CG location in the left superior pulmonary vein. The performance of this novel algorithm was prospectively tested in a validation cohort of consecutive patients undergoing repeat PVI ablation.

Results: A total of 116 patients with atrial fibrillation (AF) were recruited, 56 of whom formed the validation cohort. The interval between FFP and PVP of the left superior pulmonary vein was associated with the CG location, and an interval <5 ms predicted the presence of CG in the upper portion of the ostium with a sensitivity of 92.9% and a specificity of 96.9%. In the prospective evaluation, the interval was able to correctly predict the site of CG in 89.6% of cases.

Conclusions: The interval between FFP and PVP is a novel and accurate index that can be used to predict the CG location in the left superior pulmonary vein. An far-field atrial potential and pulmonary vein potential (FFP–PVP) interval value of ≥5 ms could be used to exclude a CG in the upper portion of the ostium in the majority of patients undergoing AF ablation.

Background

Pulmonary vein isolation (PVI) is the cornerstone of catheter ablation in many patients with atrial fibrillation (AF) (1). Pulmonary veins (PVs) are an important source of ectopic beats, initiating frequent paroxysms of AF. Ectopic beats are found frequently, particularly in the left superior pulmonary veins (LSPVs) (47.7%) (2). However, non-isolation after initial PVI during observation on the table or in patients after discharge from the hospital is not rare (3). Reisolation of the recovered PVs has demonstrated improvement in eliminating AF and clinical outcomes (4, 5). Therefore, touch-up is frequently required to eliminate conduction gaps (CGs) (6). During a repeat PVI procedure, identifying CGs is sometimes critical and difficult. Predicting the CG location before ablation can significantly reduce atrial and PV lesions and decrease operation time.

The primary goal of this study was to systematically evaluate the role of the interval between the far-field atrial potential and pulmonary vein potential (FFP–PVP interval) in identifying CGs and to develop a novel algorithm based on these findings.

Methods

This study was conducted in two phases: (1) a development and hypothesis-generating phase whereby patients' electrophysiology study and Carto 3 systems (Carto, Biosense Webster, Diamond Bar, CA, USA) data were systemically recorded and analyzed during isolation of the LSPV, and then used to create a precise algorithm to predict the CG location; and (2) a prospective validation of this new algorithm. The study was approved by the institutional review board of Fuwai Hospital.

Patient Selection

Consecutive patients with a CG during observation on the table after PVI were included to create the new algorithm. Patients with LA-LSPV reconnection during repeat PVI were recruited to verify the results. Exclusion criteria were structural heart disease, left atrial diameter ≥50 mm, long-standing (≥1 year) persistent AF, AF that occurred on the table, and a common ostium of left/right PVs. Informed consent was obtained before the procedure.

Electrophysiology Study and Catheter Ablation

Data were recorded simultaneously by a digital multichannel system (LabSystem PRO, Bard Electrophysiology, Lowell, MA, USA). Bipolar signals were filtered at 30–500 Hz, and unipolar signals were filtered at 0.05–500 Hz. Electroanatomic mapping was performed using Carto 3 systems. A 3.5-mm-tip catheter (ThermoCool SmartTouch, Biosense Webster Inc.) was used for three-dimensional electroanatomic mapping and ablation. The PV ostium was identified by selective venography and focal potential, and tagged on a three-dimensional electroanatomic map. A Lasso catheter (Biosense-Webster, Inc.) was placed at the PV ostium to observe the change in PVP and confirm PVI during radiofrequency (RF) ablation.

The PVI technique has been previously described in detail (1). In all procedures performed by an experienced operator, irrigated RF ablation using a point-by-point technique was initially performed to create right-sided continuous circular lesions around PVs and subsequently on the left side. For the left side, ablation started at the inferior aspect and continued posteriorly, followed by the superior and anterior portions. Radiofrequency ablation was performed in the posterior wall ~1 cm and in the anterior wall ~5 mm from the PV ostium (7, 8). Radiofrequency energy was delivered in a power control mode and titrated from 20 to 35 W, 43°C maximum temperature with an irrigation rate of 17–25 ml/min. The target duration of RF delivery ranged between 30 and 60 s but was terminated in advance if local electrogram attenuation was >70% or double potentials were noted (9). The general consensus is that complete electrical PV isolation is the optimal endpoint when ablation is performed at the PV-LA junction. Once the LSPV was isolated, the PVP and the location of that ablation site were specifically marked for analysis.

For the validation group, the reconstruction of three-dimensional electroanatomic mapping in LA was achieved first, and then the previous electric reconduction of isolated PVs was assessed. Subsequently, PVP activation of the LSPV was recorded, and the CG location was mapped according to our findings. Finally, ablation was performed at which the earliest activation was detected at the Lasso catheter.

Definitions

The FFP–PVP interval was defined as the time from the FFP offset to the earliest PVP onset at the Lasso catheter (Figure 1). If those two potentials overlapped and their offset or onset could not be recognized, the interval was then defined as the value of 0 ms.

FIGURE 1
www.frontiersin.org

Figure 1. The interval between FFP and PVP at the Lasso catheter was measured from the FFP offset to the earliest PVP onset. (A) FFP and PVP overlapped, so the interval between them was defined as the value of 0 ms. (B) FFP and PVP were separated, so the interval between them was 46 ms. FFP, far-field atrial potential; PVP, pulmonary vein potential; CS, coronary sinus; L1.2, Lasso catheter distal electrode; L9.10, Lasso catheter proximal electrode.

According to the anatomical structure, the ostium of the LSPV was divided into two sections: upper and lower portions. In detail, the upper portion was defined as the top half of the ostium while the lower portion was the remaining section. Otherwise, from 10 to 2 o'clock of the ostium was defined as the PV roof.

Statistical Analysis

Continuous parameters are presented as the mean ±SD and were assessed using Student's t-test or the Mann–Whitney test. Categorical variables are presented as counts (percentage) and were analyzed using the chi-square test or Fisher's exact test. Sensitivity and specificity for each variable were determined using the receiver operator characteristics (ROC) curve. All clinically relevant variables were included in the multivariable logistic regression to identify the variables that were independently associated with the presence of a second loop. A value of P < 0.05 was considered significant. SPSS, version 19.0 (SPSS, Inc., Chicago, IL, USA) was used for statistical analysis.

Results

Patient Characteristics

The cohort consisted of 112 patients [96/112 male (85.7%)]. None had long-standing (≥1 year) persistent AF. Patients' characteristics are shown in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Baseline characteristics of patients in the development and prospective cohorts.

Development of the Novel Algorithm

In the development phase, 60 patients with AF who underwent PVI ablation in Fuwai Hospital from July 2015 to July 2017 were included in the development cohort. CGs were identified in 32 (53.3%) patients in the upper portion of the ostium, and 28 (46.7%) were identified in the lower portion. In the upper portion, most of the CGs (30/32, 93.75%) were located at the roof. Within the lower portion, all the CGs were located at the junction of PV. There were no significant differences in age, sex, left atrial diameter, or CHA2DS2-VASc between the two groups (Table 1). In the logistic regression analysis, FFP–PVP interval was identified as the only factor associated with the CG location (OR = 1.55, 95% CI 1.13–2.13, P = 0.006), which indicated that the interval in the lower portion of the ostium was greater than that in the upper portion. Receiver operator characteristics analysis of the FFP–PVP interval is shown in Figure 2. The AUC was 0.958, and an interval <5 ms predicted that there was a CG in the upper portion of the ostium with a sensitivity of 92.9% and a specificity of 96.9%. Based on the above findings, a preliminary algorithm was developed. Once nonisolation occurs after circumferential PV ablation, the FFP–PVP interval should be measured. If the interval is <5 ms, the CG is likely to be located in the upper portion; otherwise, it is located in the lower portion.

FIGURE 2
www.frontiersin.org

Figure 2. Receiver operating characteristic (ROC) analysis for different intervals between far-field atrial potential and pulmonary vein potential.

Validation of Novel Algorithm Incorporating Interval Index

From September 2017 to September 2019, 56 patients with recurrent AF and reconnected LA-LSPV were recruited in the validation cohort. Individual CG was found in the lower portion of the ostium in 20 patients and the upper portion in 28 patients. Eight patients had CGs in both portions. The scatter plot of cutoff values of the FFP–PVP interval in the validation cohort is shown in Figure 3. In 18 of 20 (90.0%) patients with the CG located in the lower portion of the ostium, the FFP–PVP interval was >5 ms, whereas in 25 of 28 (89.3%) patients with the CG located in the lower portion, the FFP–PVP interval was shorter than 5 ms. As such, this cutoff value was able to correctly predict the site of CG (upper portion versus lower portion of the ostium) in 89.6% (43 of 48) of cases, with a sensitivity of 92.6% and a specificity of 85.7%. An example of a CG that was correctly predicted by the interval in a re-do patient is shown in Figure 4.

FIGURE 3
www.frontiersin.org

Figure 3. Scatter plot of the cutoff value of the interval between far-field atrial potential and pulmonary vein potential in the validation cohort.

FIGURE 4
www.frontiersin.org

Figure 4. An example of a conduction gap that was correctly predicted by the interval in a re-do patient. (A) Lasso was placed at the ostium of the LSPV, and the PV and left atrium were reconnected. The interval between FFP and PVP was 0 ms. (B) Pacing at the distal CS catheter to distinguish and confirm the existence of PVP. (C) The successful ablation target was only located at the roof of the LSPV. The blue point represents the ablation target, the white point represents the pulmonary vein ostium, and the pink point represents the mitral valve. LSPV, left superior pulmonary vein; LIPV, left inferior pulmonary vein; RSPV, right superior pulmonary vein; RIPV, right inferior pulmonary vein; PV, pulmonary vein; LAA, left atrial appendage; CS, coronary sinus; FFP, far-field atrial potential; PVP, pulmonary vein potential; L1.2, Lasso catheter distal electrode; L9.10, Lasso catheter proximal electrode.

It is also interesting to note that in patients with CG located in both portions, the FFP–PVP intervals were all shorter than 5 ms. After closing the CG at the upper portion, a new CG was identified with a delay in the timing (22.3 ± 13.0 ms); all were eventually eliminated in the lower portion.

Discussion

Major Findings

The major finding of this study is that the interval between FFP and PVP of the LSPV provided a more accurate location of the CG on the PVI ablation line. The key findings in this series are as follows: (1) The upper portion of the LSPV is the dominant conduction area for the LSPV during sinus rhythm. (2) An FFP–PVP interval <5 ms predicted the presence of CG in the upper portion of the ostium with a sensitivity of 92.9% and a specificity of 96.9%.

Possible Anatomic Considerations

Our findings might be explained by anatomical structure and interatrial conduction. Usually, Bachmann's bundle (Figure 5), known as the interatrial band, is composed of a nearly parallel alignment of myocardial strands, which accounts for its role as the prevalent interatrial conduction pathway for propagation of the sinus impulse to the anterior left atrial wall. To the right, Bachmann's bundle connects the right atrial appendage (AA), and the superior arm can extend toward the sinus node. To the left, Bachmann's bundle branches to pass around the neck of the left AA, reuniting to continue into the musculature of the lateral and posteroinferior atrial walls (10). Therefore, as for the LSPV, the earliest breakthrough usually occurs at the roof and then the impulse travels down anteriorly and posteriorly to the inferior portion. The left AA is activated simultaneously when LSPV is activated from the upper portion. Therefore, PVP and FFP in general overlap without any intervention. Once the upper portion was blocked, especially the roof, LSPV could be activated only from the lower portion where a CG existed. As a result, the activation pathway of the LSPV was prolonged, and PVP separated from FFP. Conversely, if the upper portion was not blocked, the activation pathway remained unchanged, and the FFP–PVP interval would not change. This agrees with the results of the present study on the prediction of CG in the upper portion of the ostium, using the cutoff value of an FFP–PVP interval <5 ms.

FIGURE 5
www.frontiersin.org

Figure 5. Bachmann's bundle branches to the left atrium and passes around the neck of the left atrial appendage, reuniting to continue into the musculature of the lateral and posteroinferior atrial walls. The red point represents the pulmonary vein ostium. LSPV, left superior pulmonary vein; RSPV, right superior pulmonary vein; LAA, left atrial appendage.

Strategy for Locating CGs

A previous study reported that conduction over a single gap with a longer conduction time was frequently masked by another gap with a shorter conduction time, when there were more than two gaps with different conduction times (11). In our study, we found that the FFP–PVP intervals were all shorter than 5 ms in eight patients with ≥2 CGs. After eliminating the upper portion gaps, the FFP–PVP interval was suddenly prolonged substantially and the residual gaps finally closed at the lower portion. It seems that the upper portion is the dominant conduction area of the LSPV. Thus, based on the above findings, a novel algorithm was developed for locating CGs (Figure 6). According to this algorithm, once the LSPV is not isolated or reconnected after circumferential PV ablation, the FFP–PVP interval should be measured first. If the interval is longer than 5 ms, a CG located superiorly can be excluded. If the interval is shorter than 5 ms, the CG may be located in the upper portion. When PVP with a delay was detected after eliminating the upper CG, the remaining ones located in the lower portion should be considered.

FIGURE 6
www.frontiersin.org

Figure 6. Stepwise diagnostic for locating conduction gap(s).

A Simple Step for Locating the CG

Although experienced physicians can achieve 100% acute PVI with standard mapping techniques, reconnection between the PV and left atrium is not rare in patients during observation on the table, or after discharge from the hospital (4). Previous studies (4, 12) have underlined the importance of permanent PVI in treating AF. Individuals who have received permanent PV isolation are unlikely to experience an arrhythmia recurrence. In addition, most of the atrial triggers were found to originate in the PVs, particularly in the LSPV, which accounted for 47.7% (2). Hence, CG at LSPV should be focused (1315). Guided by Nav X (11), Rhythmia (16), and the CARTO mapping system (17), several methods have been reported to allow for acquisition of multiple 3D anatomical points along with the local activation time to visually identify of the CG. However, Miyamoto et al. (11) demonstrated that this visual method could sometimes be compromised by low amplitude and/or fractionated electrograms, which requires additional effort to accurately distinguish between atrial and PV electrograms. As Benito et al. (18) and Masuda et al. (16) reported, the voltage map failed to identify gap localization in two-thirds of the cases. Furthermore, Furnkranz et al. (15) used the interval between P-wave onset and the earliest PVP at Lasso to demonstrate that the CG at the LA roof resulted in no delay or minimal delay of PV activation. However, the timing of P-wave onset to the earliest PVP varied between individuals and was greatly affected by intra-atrial conduction. In this study, we found that there was a short or no delay of PV activation in the presence of a CG at the superior LSPV. We also provided a simple and intuitive method to determine the CG(s) by comparing the FFP–PVP interval. A separation between FFP and PVP suggested that no CG was likely located superiorly. This method is especially valuable for cryoballoon ablation. The combined use of the interval could facilitate the recognition of the CG location without an additional mapping tool, which could decrease the need for touch-ups, and reduce the number of cryoablation applications.

Study Limitations

This study has several limitations. First, the sample size in the prospective cohort was relatively small and the findings should be further validated in larger populations with multiple centers. Second, although the new algorithm can accurately distinguish between the upper and the lower portion of the ostium for locating CG, the precise location could not be determined. Furthermore, this new algorithm is only applicable to the LSPV, which tends to have a CG. The applicability in other PVs should be further validated.

Conclusions

The interval between FFP and PVP is a novel and accurate index for predicting the location of a CG in LSPV, and an interval ≥5 ms could exclude a CG in the upper portion of the ostium in a majority of patients.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional Review Board of Fuwai Hospital, National Center for Cardiovascular Diseases, Chinese Academy of Medical Sciences and Peking Union Medical College. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

H-yX and X-gG conceived the study. J-dY and Y-qC analyzed the data. Z-jC collected the data. QS critically revised the article. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by a grant from the National Natural Science Foundation of China (#81670309).

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.

References

1. Ouyang F, Bänsch D, Ernst S, Schaumann A, Hachiya H, Chen M, et al. Complete isolation of left atrium surrounding the pulmonary veins: new insights from the double-Lasso technique in paroxysmal atrial fibrillation. Circulation. (2004) 110:2090–6. doi: 10.1161/01.CIR.0000144459.37455.EE

CrossRef Full Text | Google Scholar

2. Haïssaguerre M, Jaïs P, Shah DC, Takahashi A, Hocini M, Quiniou G, et al. Spontaneous initiation of atrial fibrillation by ectopic beats originating in the pulmonary veins. N Engl J Med. (1998) 339:659–66. doi: 10.1056/NEJM199809033391003

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Kuck KH, Hoffmann BA, Ernst S, Wegscheider K, Treszl A, Metzner A, et al. Impact of complete versus incomplete circumferential lines around the pulmonary veins during catheter ablation of paroxysmal atrial fibrillation: results from the gap-atrial fibrillation-German Atrial fibrillation competence network 1 trial. Circ Arrhythm Electrophysiol. (2016) 9:e003337. doi: 10.1161/CIRCEP.115.003337

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Ouyang F, Antz M, Ernst S, Hachiya H, Mavrakis H, Deger FT, et al. Recovered pulmonary vein conduction as a dominant factor for recurrent atrial tachyarrhythmias after complete circular isolation of the pulmonary veins: lessons from double Lasso technique. Circulation. (2005) 111:127–35. doi: 10.1161/01.CIR.0000151289.73085.36

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Katritsis D, Wood MA, Giazitzoglou E, Shepard RK, Kourlaba G, Ellenbogen KA. Long-term follow-up after radiofrequency catheter ablation for atrial fibrillation. Europace. (2008) 10:419–24. doi: 10.1093/europace/eun018

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Sandorfi G, Rodriguez-Manero M, Saenen J, Baluja A, Bories W, Huybrechts W, et al. Less pulmonary vein reconnection at redo procedures following radiofrequency point-by-point antral pulmonary vein isolation with the use of contemporary catheter ablation technologies. JACC Clin Electrophysiol. (2018) 4:1556–65. doi: 10.1016/j.jacep.2018.09.020

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Pappone C, Oreto G, Rosanio S, Vicedomini G, Tocchi M, Gugliotta F, et al. Atrial electroanatomic remodeling after circumferential radiofrequency pulmonary vein ablation: efficacy of an anatomic approach in a large cohort of patients with atrial fibrillation. Circulation. (2001) 104:2539–44. doi: 10.1161/hc4601.098517

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Ernst S, Ouyang F, Löber F, Antz M, Kuck KH. Catheter-induced linear lesions in the left atrium in patients with atrial fibrillation: an electroanatomic study. J Am Coll Cardiol. (2003) 42:1271–82. doi: 10.1016/S0735-1097(03)00940-9

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Pappone C, Santinelli V, Manguso F, Vicedomini G, Gugliotta F, Augello G, et al. Pulmonary vein denervation enhances long-term benefit after circumferential ablation for paroxysmal atrial fibrillation. Circulation. (2004) 109:327–34. doi: 10.1161/01.CIR.0000112641.16340.C7

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Ho SY, Ernst S. Anatomy for Cardiac Electrophysiologists: A Practical Handbook. Minneapolis, MN: Cardiotext Publishing (2012).

Google Scholar

11. Miyamoto K, Tsuchiya T, Yamaguchi T, Nagamoto Y, Ando S, Sadamatsu K, et al. A new method of a pulmonary vein map to identify a conduction gap on the pulmonary vein antrum ablation line. Circ J. (2011) 75:2363–71. doi: 10.1253/circj.CJ-11-0198

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Willems S, Steven D, Servatius H, Hoffmann BA, Drewitz I, Müllerleile K, et al. Persistence of pulmonary vein isolation after robotic remote-navigated ablation for atrial fibrillation and its relation to clinical outcome. J Cardiovasc Electrophysiol. (2010) 21:1079–84. doi: 10.1111/j.1540-8167.2010.01773.x

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Ciconte G, Mugnai G, Sieira J, Velagić V, Saitoh Y, Irfan G, et al. On the quest for the best freeze: predictors of late pulmonary vein reconnections after second-generation cryoballoon ablation. Circ Arrhythm Electrophysiol. (2015) 8:1359–65. doi: 10.1161/CIRCEP.115.002966

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Park CI, Lehrmann H, Keyl C, Weber R, Schiebeling J, Allgeier J, et al. Mechanisms of pulmonary vein reconnection after radiofrequency ablation of atrial fibrillation: the deterministic role of contact force and interlesion distance. J Cardiovasc Electrophysiol. (2014) 25:701–8. doi: 10.1111/jce.12396

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Furnkranz A, Julian JK, Schmidt B, Wohlmuth P, Tilz R, Kuck KH, et al. Ipsilateral pulmonary vein isolation performed by a single continuous circular lesion: role of pulmonary vein mapping during ablation. Europace. (2011) 13:935–41. doi: 10.1093/europace/eur067

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Masuda M, Fujita M, Iida O, Okamoto S, Ishihara T, Nanto K, et al. The identification of conduction gaps after pulmonary vein isolation using a new electroanatomic mapping system. Heart Rhythm. (2017) 14:1606–14. doi: 10.1016/j.hrthm.2017.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Dong J, Liu X, Long D, Yu R, Tang R, Lu F, et al. Single-catheter technique for pulmonary vein antrum isolation: is it sufficient to identify and close the residual gaps without a circular mapping catheter? J Cardiovasc Electrophysiol. (2009) 20:273–9. doi: 10.1111/j.1540-8167.2008.01324.x

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Benito EM, Andreu D, Mont L, Berruezo A. Correlation between functional electrical gaps identified by ultrahigh-density mapping and by late gadolinium enhancement cardiac magnetic resonance in repeat atrial fibrillation procedure. Heart Rhythm Case Rep. (2017) 3:282. doi: 10.1016/j.hrcr.2017.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: atrial fibrillation, pulmonary vein ablation/isolation, catheter ablation–atrial fibrillation, arrhythima, conduction gap

Citation: Xie H-y, Guo X-g, Yang J-d, Chen Y-q, Cao Z-j, Sun Q and Ma J (2021) Novel Clue to Locate Conduction Gaps in the Pulmonary Vein Isolation Ablation Line. Front. Cardiovasc. Med. 8:622483. doi: 10.3389/fcvm.2021.622483

Received: 28 October 2020; Accepted: 04 June 2021;
Published: 12 July 2021.

Edited by:

Shimon Rosenheck, Hebrew University of Jerusalem, Israel

Reviewed by:

Osmar Antonio Centurion, National University of Asunción, Paraguay
François Regoli, University of Zurich, Switzerland

Copyright © 2021 Xie, Guo, Yang, Chen, Cao, Sun 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: Jian Ma, majian_fuwai@126.com

Download