Skip to main content

METHODS article

Front. Cardiovasc. Med., 06 June 2023
Sec. Cardiac Rhythmology
Volume 10 - 2023 | https://doi.org/10.3389/fcvm.2023.1205966

The effect of first step right atrial mapping (FRAM) on ablation duration and fluoroscopy exposure during cavotricuspid isthmus ablation of atrial flutter

Fu Guan1* Ardan M. Saguner1,2 Alexander Breitenstein1 Mia Wang3 Nadine Molitor1 Corinna Brunckhorst1 Thomas Wolber1,4 Firat Duru1,2,4
  • 1Cardiac Arrhythmia and Electrophysiology Division, Department of Cardiology, University Hospital Zurich, Zurich, Switzerland
  • 2Center for Translational and Experimental Cardiology (CTEC), University of Zurich, Schlieren, Switzerland
  • 3Department of Biostatistics, University of Michigan, Ann Arbor, MI, United States
  • 4Center for Integrative Human Physiology, University of Zurich, Zurich, Switzerland

Aim: To investigate the clinical significance of right atrial mapping prior to cavotricuspid isthmus (CTI) ablation in patients with typical atrial flutter (AFL).

Methods: Clinical and ablation parameters were retrospectively assessed and compared in patients undergoing CTI ablation with or without a first-step right atrial mapping (FRAM) by using the CARTO 3D mapping system.

Results: CTI block by radiofrequency ablation (RFA) was achieved in all 143 patients. In the FRAM group there was a shorter ablation duration and fluoroscopy exposure compared with the non-FRAM group. CHA2DS2-VASc score was associated with higher ablation durations, more ablation applications and increased fluoroscopy exposure. Body mass index (BMI) was associated with longer ablation duration and more ablation applications. Furthermore, patients with reduced left ventricular ejection fraction (LVEF) had longer ablation durations and more fluoroscopy exposure. One patient in the non-FRAM group developed cardiac effusion after ablation. None of the patients had recurrence after 6 months of follow-up.

Conclusions: Patients with high BMI, high CHA2DS2-VASc score and reduced LVEF may benefit from the FRAM approach by reducing ablation duration, number of ablation applications and fluoroscopy exposure.

Introduction

Cavotricuspid isthmus (CTI)—dependent atrial flutter (AFL) is common, carries a high thromboembolic risk and can cause tachycardia-mediated cardiomyopathy (1). Radiofrequency ablation (RFA) of the CTI is the first-line treatment recommended by current clinical guidelines (1) with a high success rate. However, as the indications for CTI ablation have become more widespread and the clinical profile of patients more diverse, interventional electrophysiologists may be challenged with difficult CTI ablations, as shown by the wide variability in procedure and fluoroscopy times in such patients in previous publications (2, 3). The current use of three-dimensional (3D) electroanatomical mapping system allows real-time visualization of cardiac anatomy during ablation, providing unprecedented high-quality images, and thus making the real-time visualization of abnormal cardiac structures possible (4).

Activation map of the entire right atrium (RA) by using a 3D electroanatomical mapping system together with entrainment maneuver can prove that the arrhythmia is CTI-dependent. The role of RA mapping performed prior to CTI ablation (first-step right atrial mapping, FRAM) in AFL has not been systematically studied to date. The objective of this study was to investigate whether FRAM might have a positive impact to facilitate CTI ablation procedures in patients with complex anatomies and clinical features.

Methods

Study population

This was a single-center, retrospective, observational study that included patients with typical or reverse-typical CTI-dependent AFL who underwent an RFA procedure at the University Heart Center Zurich between January 2018 and December 2022. The procedures were all performed by experienced electrophysiologists in our team. We reviewed the procedural characteristics and outcomes of two different groups: the group of patients who underwent RA electroanatomical mapping prior to CTI ablation (first-step right atrial mapping, FRAM group) and the group of patients who underwent direct CTI ablation without prior RA electroanatomical mapping (non-first-step right atrial mapping, non-FRAM group). All enrolled patients had ongoing AFL at the time of the RFA procedure. We did not include patients who underwent pulmonary vein isolation during the same procedure. Patients who were unavailable for clinical follow-up within six months after the procedure were not included.

This study complied with the Declaration of Helsinki regarding investigation in humans. All patients gave written consent for data use.

Ablation procedure and follow-up

All RFA procedures were performed under local anesthesia. Analgesic sedation was used, if deemed clinically necessary. A steerable decapolar electrode catheter was used to record coronary sinus (CS) signals during AFL. A 3D electroanatomical mapping system (CARTO, Biosense Webster) was used for mapping and ablation in all patients. If necessary, fluoroscopic imaging was performed using a biplane fluoroscopy system (BICOR HS, Siemens Medical Systems). The system was uniformly set to pulse mode at 3 to 7.5 frames per second with a diameter of 23 cm image intensifier. Collimation and magnification were not utilized during the procedure.

In the FRAM group, electroanatomical mapping of RA during AFL was performed prior to CTI ablation to confirm a CTI-dependent AFL and identify the RA anatomy. A high-density mapping catheter (Pentaray, Biosense Webster) was used in this group. In the non-FRAM group, CTI-dependent atrial flutter was diagnosed on the basis of the clinical electrocardiogram in combination with entrainment maneuver. The local anatomy of the CTI in non-FRAM group was obtained retrospectively by reviewing point-by-point ablation line along the CTI by ablation catheter. The presence of a pouch in this study was defined as a depression (>2 mm in depth) within the CTI line in both groups. The pouch depth was measured both at right anterior oblique (RAO) and left anterior oblique (LAO) projections.

After confirmation of CTI-dependence, point-by-point ablation was performed from the tricuspid annulus (TA) to the inferior vena cava (IVC). An ablation catheter (Thermocool SmartTouch, Biosense Webster) was used in all patients. Mean contact force of 10 to 15 grams was targeted during ablation in both groups. During RF application, a uniform VISITAG module was used to mark the lesions. The settings were as follows: stabilization range of 25%, contact force (CF) of 3 g, label diameter of 3 mm, and point distance of 6 mm. The ablation power was set to 35–40 W. Saline flushing was set to 30 ml/min. Ablation lesion assessment was evaluated with target ablation index of 600 in combination with local potential amplitude reduction more than 50%, or elimination of local signals before ablation index reached 600. In the two groups, CTI ablation was performed using the same ablation catheter and with the same power and VISITAG settings.

All procedures were performed during AFL. After termination of AFL, RFA was continued during pacing from the proximal CS to achieve bidirectional CTI block. In the absence of first-round conduction block, the location of the conduction gap was determined by achieving bidirectional conduction block with additional RF application. Gaps were annotated to mark the location on the CARTO system. During a 20 minute waiting time, CTI block was repeatedly checked during CS proximal pacing. In the presence of reconnection, the gap was again marked and additional RF applications were applied.

After the procedure, all patients had bed rest for 6 h, and 12-lead ECG was recorded before discharge. A clinical follow-up was scheduled after 6 months with surface and Holter ECGs and for evaluation of symptoms.

Statistical analysis

Categorical variables were analyzed using the chi-square statistics and expressed as percentages. Quantitative variables were present as means with standard deviation, as appropriate for their distribution. Wilcoxon test was used for continuous variables, Chi-square for categorical variables, and Fisher's exact test for small cells. All statistical analyses were performed using the SPSS statistics (version 20.0, IBM Corp., Armonk, NY, USA). Correlation of various measurements were assessed using Spearman's correlation coefficient (R). A P value <0.05 was considered statistically significant for all tests.

Results

Patients characteristics

The baseline characteristics of the patients are summarized in Table 1. The median age was 65 ± 12 years, with CHA2DS2-VASc score was 2.1 ± 1.4. There were no differences in the baseline characteristics between the two groups.

TABLE 1
www.frontiersin.org

Table 1. Baseline patient characteristics.

Procedural outcomes

3D-mapping characterizations

The CTI length was measured as the distance between the TA and IVC in RAO view (Figure 1). The median length of CTI was 29.1 mm ± 5.9 mm. The presence of a pouch was revealed by CARTO mapping before ablation in the FRAM group or by ablation catheter mapping during CTI ablation in non-FRAM group.

FIGURE 1
www.frontiersin.org

Figure 1. CTI ablation guided by 3D mapping system with minimized fluoroscopic exposure. Cavotricuspid isthmus ablation with continuous lesions in AFL patients with/without prior RA high-density mapping [(A). non-FRAM vs. (B). FRAM)] (RAO, right anterior oblique view; LAO, left anterior oblique view). The starting point for CTI length was defined as 7 o'clock direction on LAO perspective, with the first local activation signal of annulus on TV side on RAO view (small A and big V, A1 narrow arrow). The end point for CTI length measurement on TV side was defined as 6-7 o'clock direction on LAO perspective with the last local activation signal of atrium (small A) on the Cava side (RAO) (A2, narrow arrow). The CTI was firstly identified a pouch by high-density map in FRAM group, as shown in B1 (white thick arrow). The CTI was later identified a pouch by ablation line in a patient from non-FRAM group, as shown in A1 (white thick arrow). The contiguous lesions are shown in a light red color on the screen, the blue dot was marked as the gap detected due to failure of first-pass block, or a reconnection during the waiting time after ablation, which was shown located in the pouch area in (A).

Ablation parameters

The mean RA mapping time was 12.5 ± 3.8 min in the FRAM group. There was no difference in the total procedural time between two groups, which was recorded from the start of 3D mapping to the end of the CTI ablation. Fourteen patients (21.9%) revealed at least one pouch by CARTO mapping before ablation in the FRAM group and sixteen patients (23.0%) revealed pouches by mapping during CTI ablation in the non-FRAM group (Table 2). Median ablation duration was 602 ± 303 s. Ablation duration in the FRAM group was significantly shorter compared with in the non-FRAM group (P < 0.05). In addition, patients in this group also received less RF applications (P < 0.05). The fluoroscopy exposure during ablation was lower in the FRAM group than in the non-FRAM group (P < 0.05, Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. Comparison histogram of radiofrequency parameters in two groups. The distribution of ablation parameters in both groups including RF pulse applications (A), RF duration (B) and fluoroscopic exposure (C) Data were shown in Table 2.

TABLE 2
www.frontiersin.org

Table 2. Comparison of 3D mapping and ablation radiofrequency data.

Univariate and multivariate linear regression analysis were performed to identify the independent correlates of ablation duration, ablation pulse application and fluoroscopy exposure. CHA2DS2-VASc score was associated with higher ablation durations (r = 0.21, P < 0.05), more ablation applications (r = 0.24, P < 0.05) and increased fluoroscopy exposure (r = 0.20, P < 0.05). Body mass index (BMI) was associated with longer ablation durations (r = 0.80, P < 0.05) and more ablation applications (r = 0.22, P < 0.05). Furthermore, patients with reduced left ventricular ejection fraction (LVEF) had longer ablation durations (r = 0.55, P < 0.05) and more fluoroscopy exposure (r = 0.19, P < 0.05) (Figure 3). Meanwhile, we performed a comparative analysis between patients who were overweight and those with a normal body weight, and patients with higher CHA2DS2-VASc score and lower CHA2DS2-VASc score, as well as a trend test between patients with severe LVEF reduction, mild LVEF reduction and normal LVEF. Figure 4 showed the results, which were consistent with the above correlation analysis.

FIGURE 3
www.frontiersin.org

Figure 3. Correlation presentation by spearman coefficient (A) Association between patients’ clinical features and ablation parameters determined by coefficient with regression line. (A) The relationship between BMI and number of RF application. (B) The relationship between CHA2DS2-VASc score and number of RF application. (C) The relationship between CHA2DS2-VASc score and fluoroscopy exposure. (D) The relationship LVEF and fluoroscopy exposure.

FIGURE 4
www.frontiersin.org

Figure 4. Comparison presentation of dichotomous correlation factors. Comparative analysis between patients with different clinical features. (A) Obese patients with BMI ≥30 received more RF applications than with normal body weight (21.5 ± 8.5 and 18.1 ± 7.9, P = 0.0094). (B) Patients of CHA2DS2-VASc score ≥ 2 received more RF applications compared with patients of CHA2DS2-VASc score<2 (20.3 ± 8.7 and 16.6 ± 6.7, P = 0.0054). (C): Patients of CHA2DS2-VASc score ≥ 2 had more fluoroscopic exposure compared with patients of CHA2DS2-VASc score<2 (69.5 ± 33.0 µGy.m2 and 22.2 ± 10.7 µGy.m2, P = 0.0206). (D): Patients with normal, mild reduced or severely reduced LVEF manifested a trend of increased fluoroscopic exposure (40.0 ± 14.7 µGy.m2, 48.0 ± 32.6 µGy.m2, 104.0 ± 62. 2 µGy.m2, respectively, P = 0.0004).

Reconnection during the waiting time occurred in 6 of 64 (9%) patients in the FRAM group and in 9 of 79 (11%) patients in the non-FRAM group (P = 0.55). Reconnected gaps were located at the pouch area in 11 of 15 (73%) patients (Figure 1B). Right AFL aside from CTI-dependent AFL was induced during electrophysiological study after CTI ablation in 4 of 64 (6%) patients in the FRAM group and in 7 of 79 (9%) patients in the non-FRAM group. Among those patients with another right-sided AFL, 5 of 7 (71%) were remapped with the mapping catheter.

None of the patients developed cardiac tamponade during ablation, and only one patient in non-FRAM group developed asymptomatic medium cardiac effusion after ablation and recovered with conservative management. Bidirectional block was achieved in all patients at the end of the procedure. None of the patients had clinical or ECG recurrence of AFL after 6 months of follow-up in both groups.

Discussion

Our study demonstrated that in patients with CTI-dependent AFL, real-time guided linear ablation with RA mapping (FRAM) prior to RFA was associated with reduced ablation durations and applications as well as fluoroscopy exposure during ablation. For patients with high BMI, high CHA2DS2-VASc scores and reduced LVEF, these benefits were more prominent.

A number of studies have recently reported CTI anatomical abnormalities, such as pouches detected by 3D transesophageal echocardiography (TEE), intracardiac ultrasound (ICE) and cardiac magnetic resonance imaging, which were associated with prolonged ablation times and energy delivery during CTI ablation (57). The authors pointed out that anatomical abnormalities of the CTI led to prolonged ablation time. RA pouches were usually associated with a prominent Eustachian valve or Eustachian ridge (5, 6) and were previously found to be more prevalent in patients with enlarged LA and reduced left ventricular systolic function (8). However, there are no reports on how to simply identify this particular anatomical variant of pouch in this population without the above-mentioned imaging modalities.

Our results suggest that RA mapping prior to CTI ablation was associated with a significantly reduced ablation time and fluoroscopy exposure, which may facilitate the ablation procedure and hence interventional success rate. In patients with high BMI and heart failure, there was a prolonged time to reach CTI block, and a high number of ablation applications. Those patients were also exposed to increased amount of fluoroscopy. This may be at least partially related to poor catheter stability due to poor respiratory compensation during catheter ablation aside from CTI abnormalities. A higher incidence of AFL has also been reported in patients with chronic obstructive pulmonary disease (COPD) (9). Furthermore, it has been shown that patients with heart failure have a more remodeled CTI, manifesting as prolonged or concave deformation of this region (3). A study using modified long-axis 2D TTE view also reported that the presence of a long CTI, extensive Eustachian ridge, RA pouch and the presence of tricuspid regurgitation were associated with ablation difficulties (10). In our study, lower fluoroscopic exposure and shorter ablation duration were documented in the population with the FRAM protocol, suggesting that FRAM might guide efficient ablation by identifying abnormal CTI anatomy intraoperatively. Although 3D TEE was reported to detect specific pouch structures, this detailed examination is not routine in most clinical centers due to its time-consuming nature (5). ICE is helpful for identification of CTI structures during the intervention (6), but it adds further cost to the procedure. For these reasons, FRAM protocol with CARTO mapping may be considered a convenient approach facilitating the ablation procedure.

Notably, our study revealed that the CHA2DS2-VASc score showed correlations with increased difficulty of CTI ablation (11). Risk factors present in the CHA2DS2-VASc score have been associated with atypical AFL including hypertension, cardiomyopathy, prior atrial fibrillation, prior surgery, and COPD (12, 13). Hence, there is an increased likelihood that these patients develop both typical and atypical AFL, either with concurrent atypical AFL or presenting a clinical ECG that is misinterpreted as typical AFL (14, 15). In these cases, a single CTI ablation with bidirectional block was not able to solve the problem. The FRAM protocol in our study not only showed advantages for patients with high CHA2DS2-VASc score in terms of ablation parameters, but also enabled detection and ablation of concurrent atypical AFL.

This study is a small, single-center, observational cohort study. The correlation coefficient sizes are due to the size of our study population. Our follow-up was rather short and consisted of only ECG/Holter recording and symptom evaluation. The recurrence identification was only based on readmission due to typical AFL. Thus, our follow-up may overestimate the success rate of CTI ablation.

In conclusion, in patients with CTI-dependent AFL, the FRAM approach for CTI ablation, which was originally intended for practical training of our electrophysiology fellows in our center, was associated with reduced ablation duration and fluoroscopy exposures. These advantages may be especially important in patients with high BMI, high CHA2DS2-VASc scores and reduced left ventricular systolic function.

Data availability statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving human participants were reviewed and approved by University Hospital Zurich. The patients/participants provided their written informed consent to participate in this study.

Author contributions

FG and FD contributed to conception and design of the study. FG organized the database and performed the statistical analysis. FG wrote the first draft of the manuscript. AS, AB, NM, CB and TW wrote sections of the manuscript. MW contributed to the statistical writing of this part. All authors contributed to the article and approved the submitted version.

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.

References

1. Hindricks G, Potpara T, Dagres N, Arbelo E, Bax JJ, Lundqvist CB, et al. 2020 ESC guidelines for the diagnosis and management of atrial fibrillation developed in collaboration with the European association for cardio-thoracic surgery (EACTS): the task force for the diagnosis and management of atrial fibrillation of the European society of cardiology (ESC) developed with the special contribution of the European heart rhythm association (EHRA) of the ESC. Eur Heart J. (2021) 42(5):373–498. doi: 10.1093/eurheartj/ehaa612

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Santoro A, Baiocchi C, Sisti N, Zaca V, Pondrelli CR, Falciani F, et al. Catheter inversion during cavotricuspid isthmus catheter ablation: the new shaft visualization catheter reduces fluoroscopy use. J Arrhythm. (2021) 37(5):1311–7. doi: 10.1002/joa3.12596

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Asirvatham SJ. Correlative anatomy and electrophysiology for the interventional electrophysiologist. J Cardiovasc Electrophysiol. (2009) 20(1):113–22. doi: 10.1111/j.1540-8167.2008.01344.x

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Marazzato J, Marazzi R, Doni AL, Angeli F, Bagliani G, Leonelli MF, et al. Mapping and ablation of atypical atrial flutters. Card Electrophysiol Clin. (2022) 14(3):471–81. doi: 10.1016/j.ccep.2022.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Regoli R, Faletra F, Marcon S, Leo LA, Dequarti MC, Caputo ML, et al. Anatomic characterization of cavotricuspid isthmus by 3D transesophageal echocardiography in patients undergoing radiofrequency ablation of typical atrial flutter. Eur Heart J Cardiovasc Imaging. (2018) 19(1):84–91. doi: 10.1093/ehjci/jew336

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Shimizu Y, Yoshitani K, Murotani K, Kujira K, Kurozumi Y, Fukuhara R, et al. The deeper the pouch is, the longer the radiofrequency duration and higher the radiofrequency energy needed-cavotricuspid isthmus ablation using intracardiac echocardiography. J Arrhythm. (2018) 34(4):410–7. doi: 10.1002/joa3.12075

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Lim KT, Murray C, Liu H, Weerasooriya R. Pre-ablation magnetic resonance imaging of the cavotricuspid isthmus. Europace. (2007) 9(3):149–53. doi: 10.1093/europace/eul187

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Sorgente A, Moccetti T. Abnormal right atrial pouch in a patient with heart failure and cavotricuspid isthmus-dependent atrial flutter. Heart Rhythm. (2011) 8(4):639–40. doi: 10.1016/j.hrthm.2010.02.030

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Hayashi T, Fukamizu S, Hojo R, Komiyama K, Tanabe Y, Tejima T, et al. Prevalence and electrophysiological characteristics of typical atrial flutter in patients with atrial fibrillation and chronic obstructive pulmonary disease. Europace. (2013) 15(12):1777–83. doi: 10.1093/europace/eut158

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Pedro MA, Damian SQ, Cabrera JA, Cabrera JA, Farre J, Rubio MJ, de Agustin JA, et al. Two-dimensional echocardiographic features of the inferior right atrial isthmus: the role of vestibular thickness in catheter ablation of atrial flutter. Eur Heart J Cardiovasc Imaging. (2014) 15(1):32–40. doi: 10.1093/ehjci/jet112

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Cetinkal G, Kocas BB, Ser OS, Kilci H, Keskin K, Ozcan SN, et al. Assessment of the modified CHA2DS2VASc risk score in predicting mortality in patients hospitalized with COVID-19. Am J Cardiol. (2020) 15(135):143–9. doi: 10.1016/j.amjcard.2020.08.040

CrossRef Full Text | Google Scholar

12. Ipek EG, Marine J, Yang E, Habibi M, Chrispin J, Spragg D, et al. Predictors and incidence of atrial flutter after catheter ablation of atrial fibrillation. Am J Cardiol. (2019) 124(11):1690–6. doi: 10.1016/j.amjcard.2019.08.026

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Roca-Luque I, Rivas-Gándara N, Dos-Subirà L, Francisco-Pascual J, Pijuan-Domenech A, Perez-Rodon J, et al. Predictors of acute failure ablation of intra-atrial re-entrant tachycardia in patients with congenital heart disease: cardiac disease, atypical flutter, and previous atrial fibrillation. J Am Heart Assoc. (2018) 7(7):e008063. doi: 10.1161/JAHA.117.008063

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Huang HK, Chang SL, Lin YJ, Lo LW, Hu YF, Chung FP, et al. Atrial fibrillation originating from superior vena cava with atrial flutter-electrocardiogram pattern. Pacing Clin Electrophysiol. (2017) 40(7):754–61. doi: 10.1111/pace.13092

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Chugh A, Latchamsetty R, Oral H, Elmouchi D, Tschopp D, Reich S, et al. Characteristics of cavotricuspid isthmus-dependent atrial flutter after left atrial ablation of atrial fibrillation. Circulation. (2006) 113(5):609–15. doi: 10.1161/CIRCULATIONAHA.105.580936

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: arrhythmia, atrial flutter, catheter ablation, electroanatomical mapping, outcome

Citation: Guan F, Saguner AM, Breitenstein A, Wang M, Molitor N, Brunckhorst C, Wolber T and Duru F (2023) The effect of first step right atrial mapping (FRAM) on ablation duration and fluoroscopy exposure during cavotricuspid isthmus ablation of atrial flutter. Front. Cardiovasc. Med. 10:1205966. doi: 10.3389/fcvm.2023.1205966

Received: 14 April 2023; Accepted: 22 May 2023;
Published: 6 June 2023.

Edited by:

Roberto Vazquez-Padron, University of Miami, United States

Reviewed by:

Martin Manninger, Medical University of Graz, Austria
Philip L. Mar, Saint Louis University, United States

© 2023 Guan, Saguner, Breitenstein, Wang, Molitor, Brunckhorst, Wolber and Duru. 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: Fu Guan fu.guan@usz.ch

Download