The Impact of Comorbidities on Pulmonary Function Measured by Spirometry in Patients After Percutaneous Cryoballoon Pulmonary Vein Isolation Due to Atrial Fibrillation
Monika Różycka-Kosmalska, Marcin Kosmalski, Michał Panek, Alicja Majos, Izabela Szymczak-Pajor, Agnieszka Śliwińska, Jacek Kasznicki, Jerzy Krzysztof Wranicz, Krzysztof Kaczmarek

TL;DR
This study examines how comorbidities affect lung function after a heart procedure for atrial fibrillation.
Contribution
It identifies subgroups of patients where the procedure impacts respiratory function differently.
Findings
Overall spirometry parameters did not change significantly after the procedure.
Patients with CAD and high BMI showed a decline in MEF75, while those with HF and low EF showed improved FEV1/FVCex.
Abstract
Background/Objectives: Pulmonary vein isolation (PVI) via cryoballoon ablation (CBA) is a recommended therapeutic strategy for patients with symptomatic paroxysmal and persistent atrial fibrillation (AF) who are refractory to antiarrhythmic drugs. Although PVI has demonstrated efficacy in reducing AF recurrence and improving patients’ quality of life, its impact on respiratory function is not well understood, particularly in patients with comorbid conditions. The aim of the study was to search for functional predictors of the respiratory system in the process of evaluating the efficiency of clinical assessment of CBA in patients with AF. Methods: We conducted a prospective study on 42 patients with symptomatic AF who underwent CBA, assessing their respiratory function through spirometry before and 30 days after the procedure. Exclusion criteria included pre-existing lung disease and…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8- —Medical University of Lodz
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsAtrial Fibrillation Management and Outcomes · Cardiac Arrhythmias and Treatments · Venous Thromboembolism Diagnosis and Management
1. Introduction
Atrial fibrillation (AF) is the most common cardiac arrhythmia, characterized by rapid, irregular atrial excitation that leads to dyssynchronous atrial contraction and uneven ventricular activation [1]. It affects over 59 million people worldwide, with prevalence projected to rise significantly by 2050 [2]. AF increases the risk of stroke, heart failure, and death and impairs quality of life [3]. AF is classified as duration into paroxysmal (PAF), persistent (AF-PE), and permanent with treatment goals focusing on symptom control, reducing long-term complications, and improving quality of life. Catheter ablation (CA) prevents AF recurrences, reduces AF burden, and improves the quality of life in symptomatic PAF and AF-PE patients who do not tolerate or do not respond to antiarrhythmic drugs [4]. One of the methods of CA is pulmonary vein isolation (PVI) using a balloon-mounted cryoablation system. Complete pulmonary vein occlusion during balloon ablation has been demonstrated to predict successful electrical isolation [5]. The European Society of Cardiology and Heart Rhythm Society Guideline identify PVI as the main strategy for treating symptomatic, drug-refractory AF. They also note that cryoballoon CA (CBA) can be an effective alternative to the traditional radiofrequency ablation. PVI is usually performed via the inferior vena cava through the access of the femoral veins [6,7].
The modern transvenous cryocatheter system includes a deflectable hollow shaft with a cooling tip, ring electrodes, and a thermocouple. Ablation is performed by delivering pressurized cryorefrigerant from an external console through a thin injection tube. The refrigerant absorbs heat from the heart tissue via convective cooling, causing cellular injury through ice crystal formation and microcirculatory failure [8,9].
PVI with high-energy ablation is effective for AF, with success rates of 50–95%. However, it is invasive and carries risks such as cardiac tamponade, thromboembolism, pulmonary vein stenosis, atrioesophageal fistula, and other complications [10,11,12]. Mugnai et al. reported a serious adverse event rate of 2.9%, including vascular complications (1.1%), tamponade (1.0%), and thromboembolism (0.3%). Rare complications like atrioesophageal fistula and nerve palsy occurred in 0.1% of cases, with no procedure-related deaths. The overall complication rate was similar for both radiofrequency and cryoballoon techniques [13]. A multicenter study in Poland found a low adverse event rate of 6.4%, including mainly local complications, and a low incidence of post-procedure atrial flutter or tachycardia (5%). Most procedures used radiofrequency ablation (89%) [14]. There is a lack of data on the effect of CBA on lung function. Furthermore, it is not known whether or to what extent this relationship is influenced by comorbidities and cardiac dysfunction.
The aim of this study was to search for functional predictors of the respiratory system in the process of evaluating the efficiency of clinical assessment of CBA in patients with AF.
2. Materials and Methods
2.1. Study Population
This study included a total of 42 consecutive patients (29 males and 13 females) aged between 36 and 77 years (mean: 51.25; range: 36–77) who were qualified for PVI as the first interventional antiarrhythmic treatment for symptomatic atrial fibrillation at the Department of Electrocardiology, Central University Hospital of the Medical University of Lodz, between 2020 and 2021. The inclusion criteria covered the following: (1) diagnosed FA (PAF/AF-PE) before PVI; (2) age over 18 years; (3) indications for the PVI procedure (lack of effectiveness of antiarrhythmic therapy, side effects caused by antiarrhythmic therapy, a serious threat to the patient’s life posed by arrhythmia); (4) presence of an active inflammatory condition; (5) cryoballoon as the primary planned PVI technique; (6) eligibility and willingness for spirometry testing; (7) signed informed consent to participate in this study. The exclusion criteria included the following: (1) lung diseases such as pneumonia, bronchial asthma, chronic obstructive pulmonary disease (COPD), and lung cancer; (2) phrenic nerve injury during CB-PVI; (3) previous PVI; (4) clinical signs and symptoms of respiratory insufficiency. PAF was defined according contemporaneous guidelines [4]. All patients provided written informed consent. This study was approved by the Medical University of Lodz Committee on the Ethics of Research in Human Experimentation (approval number RNN/182/17/KE) and was conducted in accordance with the Helsinki Declaration for human research. This study was not registered as a clinical trial because it was a non-interventional study that only added a respiratory function test to routine medical examinations. Spirometry was performed in each qualified patient before PVI and repeated 30 days after the procedure.
2.2. Spirometry
Respiratory function tests were performed according to European Respiratory Society (ERS) and the American Thoracic Society (ATS) standards [15]. All spirometry measurements were performed using the Lungtest 1000 device (producent: MES Sp. z o.o. of Krakow, Poland). Prior to testing, patients were instructed to breathe normally, and each measurement was repeated at least three times to ensure reproducibility. The best score of three spirometry readings was selected for analysis, following the guidelines of the Polish Society of Lung Diseases [16]. During spirometry, the following parameters were assessed: FEV_1_/FVC_ex_ (forced expiratory volume during the first second of expiration to forced vital capacity), FEV_1_ (forced expiratory volume during the first second of expiration), FVC_ex_ (expiratory forced vital capacity), PEF (peak expiratory flow), and MEF75, MEF50, and MEF25 (maximal expiratory flow at 75%, 50%, and 25% of FVC) to assess small and large airway function. For each parameter, the percentage of predicted values was calculated based on reference values adjusted for age, sex, height, and ethnicity, according to standard reference equations. The predicted values served for comparison to assess the degree of airway obstruction or restriction in individual patients.
2.3. Pulmonary Vein Isolation
PVI was performed percutaneously using the cryoballoon technique [17], according to the standard of care at our institution. The procedure began with the insertion of two vascular sheaths into the right femoral vein under sterile conditions. Prior to the transseptal puncture procedure, transesophageal echocardiography (TEE) was performed to exclude the presence of a thrombin in the left atrial appendage, thereby reducing the risk of embolic complications. The transseptal puncture was then carried out under combined guidance with fluoroscopy and TEE to ensure precise access to the left atrium. A dedicated transseptal sheath for the cryoballoon procedure was inserted into the left atrium and continuously flushed with heparinized saline to prevent clot formation. Following standard preparation, the balloon catheter was advanced into the pulmonary veins. The catheter was positioned to occlude each pulmonary vein individually, confirmed by contrast injection and imaging. Once proper occlusion was achieved, cryoablation was performed for a duration ranging from 180 to 300 s, with the application time adjusted based on the minimum temperature reached during each freeze cycle—with colder temperatures typically requiring shorter application times to optimize safety and efficacy. Throughout the entire procedure—from the transseptal puncture to the removal of catheters from the left atrium—the activated clotting time (ACT) was maintained above 300 s to minimize thromboembolic risk. The ACT was monitored regularly and adjusted with heparin infusion as needed.
2.4. Statistical Analysis
The measured respiratory parameters were expressed as medians (quartile 1; quartile 3). The null hypothesis that the data in the tested population would be normally distributed was rejected using the Shapiro–Wilk test. The relationship between respiratory parameters measured before and 30 days after ablation was evaluated using Wilcoxon signed-rank two-tailed test. The differences in the spirometry parameters, defined as the “Δ = (value on day 30)—(value on day 0)” for each spirometry parameter, analyzed in relation to selected coexisting factors, were assessed using the Mann–Whitney U test or Kruskal–Wallis test in the case of multiple comparisons. All analyses were carried out using the Statistica 13PL software (StatSoft Polska Sp. z o.o. Cracow, Poland). The significance level was assumed to be 0.05 in all the statistical tests.
3. Results
Out of a total of 42 consecutive patients included in this study, only 27 patients completed the entire study, of whom 14 patients did not return for follow-up visits and 1 patient died during the study without any relation to the PVI procedure. Median age in the whole group was 58 years (IQR, 48.0–69.0; min.-max., 38.0–77.0) with a median body mass index (BMI) of 28.7 kg/m^2^ (IQR 25.7–32.5, min.-max. 21.5–48.4) and median EF of 56.0% (IQR, 52.0–62.0%; min.-max., 22.0–69.0%).
The characteristics of patients who completed the entire study are presented in Table 1.
The comparison of spirometry parameters pre-PVI and 30 days after PVI for the whole study group is presented in Table 2; no statistically significant differences were observed in this approach to this study.
A different approach (Δ of each spirometry parameter) was applied for analysis according to the occurrence of comorbidities. The analysis of Δ of the assessed spirometry parameters before and 30 days after PVI in relation to BMI, coronary artery disease (CAD), heart failure (HF), and ejection fraction (EF) value was performed (data presented in Supplementary Materials Tables S1–S4). As shown in Figure 1 and Figure 2, among the tested parameters, Δ MEF75 (l/s) (p = 0.02) was significantly elevated in the group with BMI < 30 kg/m^2^ as compared to the group with BMI ≥ 30 kg/m^2^. A similar statistically significant difference was found in Δ MEF75 (% of predicted) (p = 0.02). As shown in Figure 3 and Figure 4, in the group with CAD, Δ MEF75 (l/s) (p = 0.04) and Δ MEF75 (% of predicted) (p = 0.03) were markedly lower compared to those in the group without CAD. Regarding the comparison of the group with HF versus the group without HF, presented in Figure 5 and Figure 6, Δ FEV_1_/FVCex (p = 0.05) and Δ FEV_1_/FVCex (% of predicted) (p = 0.03) showed a significant increase in the group with HF. Figure 7 compares the parameters between the group with EF < 50% and the group with EF ≥ 50%. The results reveal a significant increase in Δ FEV_1_/FVCex (p = 0.02) in the group with EF < 50%, suggesting that the change in this parameter was notably greater in patients with reduced ejection fraction. Similarly, Figure 8 depicts the comparison of parameters between the two groups. The findings demonstrate a significant increase in Δ FEV_1_/FVCex (% of predicted) (p = 0.05) in the EF < 50% group, indicating that the change in this parameter was also significantly higher in patients with lower EF.
We did not observe statistically significant differences for other comorbidities such as type 2 diabetes mellitus (T2DM), dyslipidemia (DLP), and hypertension (HA) as well as EHRA (European Hearth Rhythm Association) score and AF type (Tables S5–S9 in Supplementary Materials).
4. Discussion
Studies have shown a correlation between impaired lung ventilation, as confirmed by spirometry, and the risk of AF. This applies to individuals both with and without diagnosed lung diseases like COPD or bronchial asthma. Furthermore, it has been shown that this relationship is independent of other cardiovascular risk factors and concerns reduced values of FEV_1_, FVC, and the FEV_1_/FVC ratio, as well as the percentage of predicted FEV_1_ and FVC [18,19,20,21,22,23]. The exact mechanisms linking decreased lung function and AF are unclear, but potential factors include hypoxia, systemic inflammation, increased sympathetic activity, and treatments like beta-2 agonists or steroids [24]. Recent research suggests that ectopic beats triggering AF often originate in the pulmonary veins, possibly due to changes in gas composition or pulmonary hypertension. Even mild to moderate reductions in lung function have been associated with AF, and reduced lung capacity is also an independent predictor of coronary artery disease and stroke. The connection may involve shared pathways like atherosclerosis, influenced by early-life exposures or vascular and airway diseases. The close anatomical and physiological relationship between the lungs, the heart, and blood vessels suggests that damage to one component can affect overall cardiovascular health [25,26,27,28,29,30].
However, there is no data assessing the effect of CA on lung function in patients with AF. In a study conducted by Oguri et al., it was found that CA performed by various techniques, including radiofrequency catheter ablation (RFCA), CBA, hot balloon ablation, and laser balloon-mediated ablation, significantly improved VC, FVC, and PEF in patients with PAF and non-PAF defined as persistent AF (defined as AF-PE or long-standing AF-PE during 6-month observation period). In the patients with PAF, a significant increase in FEV_1_ was observed. CBA significantly improved VC and FVC. A significant improvement in these pulmonary parameters was specifically observed in patients with PAF treated with a CBA protocol but not with RFCA or hot balloon ablation. A significant decrease in FEV_1_ was observed with hot balloon ablation. Significant improvements in pulmonary parameters were seen specifically in patients with PAF who underwent CBA [31]. Additionally, we conducted an analysis of the impact of the CBA procedure on spirometry parameters to determine whether it might cause ventilation disorders in short-term observation and, if present, the possibility of applying appropriate pharmacotherapy. However, it should be emphasized that we did not find a significant effect of selected ablation techniques on the spirometry results in the entire cohort of patients with PAF and AF-PE. It is worth noting that, to minimize confounding factors, we excluded patients with previously diagnosed lung diseases and smokers from our study.
It has been proven that, apart from lung diseases such as asthma or COPD, factors such as age, sex, race, HA, congestive HF, CAD, obesity, T2DM, DLP, physical inactivity, obstructive sleep apnea, and smoking are associated with increased odds of developing AF [32,33,34,35,36,37]. Moreover, the Prospective Global GLORIA-AF Registry revealed that these comorbidities influence the management and long-term prognosis of patients with AF [38]. Despite the established association between cardiovascular risk factors such as obesity, dyslipidemia, T2DM, HA, HF, CAD, and the presence of AF, data assessing the impact of these diseases on the efficacy and safety of CA, particularly PVI, in both PAF and AF-PE groups are still lacking.
An annual observation conducted by Sargent et al. showed that the PVI procedure, performed using either CBA or radiofrequency ablation, was more successful in patients with lower weight and BMI, particularly those with a BMI < 27 kg/m^2^ and in those without obstructive sleep apnea. The only independent predictor of the primary outcome (successful isolation of the pulmonary veins and freedom from AF without repeated ablation or ongoing antiarrhythmic therapy at 12 months) was BMI. Furthermore, over 75% of complications occurred in patients with a BMI ≥ 27 kg/m^2^ [39]. It should be emphasized that, despite the proven association between obesity and long-term success of pulmonary vein isolation performed using transcatheter radiofrequency for drug-refractory AF [40], it has been shown that for obese patients, CBA is as safe and effective as RFCA. The significantly shorter procedure time for CBA may minimize potential obesity-related complications. However, the lower contrast medium quantity and fluoroscopy dose in RFCA should be considered. AF recurrence rates are comparable between CBA and RFCA [41]. However, some data indicate that CBA in overweight and obese patients is safe, with similar levels of complications and AF recurrence rates compared to those in patients of healthy weight. However, obese patients and operators are exposed to higher radiation doses [42,43,44]. Our results indicate that a BMI ≥ 30 kg/m^2^ significantly impairs respiratory ventilation in patients undergoing the CBA procedure, due to significant decrease in the MEF_75_ and % of predicted MEF_75_ values. It is important to note that the MEF_75_ value reflects the functional state of the medium and small airways, and a decrease in this indicator is common even in the early stages of cardiovascular diseases [45].
Similar observations were also made for patients with T2DM who underwent CBA [46]. It has been demonstrated that T2DM is an independent predictor of AF recurrence after catheter ablation, particularly in patients with persistent AF [47]. Recent data indicate that baseline HbA1c level is an independent predictor of AF recurrence following cryoablation, both in patients with and without diabetes [48]. Moreover, both existing T2DM and newly diagnosed T2DM were identified as independent risk factors for AF recurrence following ablation. This study underscores the urgent need for careful management of diabetes mellitus in individuals undergoing AF ablation to reduce the likelihood of arrhythmic relapse [49]. Considering the significant number of patients with T2DM and obesity, it is important to note the possibility of undiagnosed obstructive sleep apnea (OSA) in the analyzed group [50,51]. It is worth emphasizing that OSA may significantly impact spirometry parameters, and due to factors such as regions of interest (ROIs) and hypoxemia, it could lead to AF [52,53]. Nevertheless, it is worth emphasizing that our study did not demonstrate any influence of T2DM on spirometry parameters.
Despite the proven beneficial effect of CBA for AF in patients with heart failure with mid-range ejection fraction (HFmrEF) and heart failure with preserved ejection fraction (HFpEF) [54,55], no studies have assessed the influence of heart function on short- and long-term prognosis after this procedure. In our study, we found a beneficial effect of the CBA procedure on spirometry parameters such as FEV_1_/FVC_ex_ value in the group of patients with HF and % of predicted FEV_1_/FVC_ex_. A similar positive effect was noted in the group of patients with EF <50%, where a significant improvement in the index was observed. It is highlighted that this observation is clinically significant because, in patients with HF, pulmonary function assessment, particularly spirometry, can serve as a non-invasive and sensitive method for detecting and monitoring early pulmonary congestion. This approach could complement the measurement of natriuretic peptides (BNP and NT-proBNP). Such diagnostic procedures provide an opportunity for early pharmacological and non-pharmacological interventions. The respiratory system is an integral part of the cardiorespiratory interactions in HF, which is why pulmonary function tests, including spirometry and other methods (e.g., diffusing capacity of the lungs for carbon monoxide—DLCO measurement), should be considered complementary to clinical assessment, echocardiography, cardiac biomarkers such as natriuretic peptides, and chest imaging. Furthermore, pulmonary function parameters provide prognostic information in patients with HF and can assist in treatment decisions. Therefore, spirometry and DLCO should be recommended as part of the initial diagnosis and monitoring of patients with HF. Longer follow-up periods, such as one year, using pulmonary function tests may help identify risk markers for the CBA procedure [56].
Data on the impact of CAD on the effectiveness and safety of CBA suggest that CBA is effective for treating AF in patients with stable CAD. The presence of CAD does not appear to affect AF recurrence after CA [57]. In our study, we demonstrated for the first time that the presence of CAD before PVI is a risk factor for deterioration of lung ventilation during a 30-day follow-up after PVI according to a significant reduction in MEF_75_ and the % of the predicted FEV_1_/FVC_ex_ value. We did not find any significant effect of HA and DLP, as well as the type of AF and EHRA score on spirometry parameters within 30 days after PVI.
Study Limitations
A significant limitation of this study is the small sample size, which is largely due to the limited access to various PVI techniques and the coexistence of lung diseases in patients undergoing ablation procedures. Another limitation is the relatively short observation period, but this is related to the focus on assessing the indications for the rapid diagnosis of ventilation disorders in patients after PVI, as well as the ability to implement appropriate pharmacotherapy as soon as possible. Therefore, we plan a long-term (one-year) follow-up of patients participating in this study to assess PVI complications associated with possible pulmonary remodeling, taking into account comorbidities and cardiometabolic risk markers.
It should also be mentioned that this study excluded individuals with known lung disease and smokers, which may be appropriate for internal validity but reduces the generalizability to the typical AF population. Additionally, the presence of OBS and its correlation with T2DM and obesity represent a major unmeasured variable that could substantially bias the respiratory results. Moreover, time to ablation is a crucial factor in assessing the consequences of PVI. Available data suggest that a shorter time to ablation (≤1 year) is associated not only with a lower risk of repeat ablation and new cardioversion but also with a reduced risk of cardiovascular hospitalization compared to late ablation (>1 year). Additionally, shorter time to ablation correlates with a better CHA_2_DS_2_-VASc score, lower incidence of HF, and higher ejection fraction (EF). Given the observed associations with PVI, long-term follow-up appears warranted [58].
5. Conclusions
In our study, we observed that CBA for AF was not associated with significant changes in spirometry parameters during a one-month observation period. However, we observed that in obese patients, as well as in those with CAD, there were indications that post-procedural spirometry parameters might deteriorate significantly. This preliminary observation suggests the potential need for close monitoring of respiratory function in these patient groups following the procedure to enable timely implementation of appropriate pharmacotherapy. On the other hand, there is a hypothesis that PVI may have the potential to improve lung function in patients with HF and reduced EF. Given the limited sample size and the preliminary nature of these findings, further investigation in larger cohorts with extended follow-up is necessary to better understand these associations and evaluate possible preventive and therapeutic strategies for at-risk populations.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Bizhanov K.A. Abzaliyev K.B. Baimbetov A.K. Sarsenbayeva A.B. Lyan E. Atrial fibrillation: Epidemiology, pathophysiology, and clinical complications (literature review)J. Cardiovasc. Electrophysiol.20233415316510.1111/jce.1575936434795 · doi ↗ · pubmed ↗
- 2Linz D. Gawalko M. Betz K. Hendriks J.H. Lip G.Y.H. Vinter N. Guo Y. Johnsen S. Atrial fibrillation: Epidemiology, screening and digital health Lancet Reg. Health Eur.20243710078610.1016/j.lanepe.2023.10078638362546 PMC 10866942 · doi ↗ · pubmed ↗
- 3Ko D. Chung M.K. Evans P.T. Benjamin E.J. Helm R.H. Atrial Fibrillation: A Review JAMA 202533332934210.1001/jama.2024.2245139680399 PMC 11774664 · doi ↗ · pubmed ↗
- 4Joglar J.A. Chung M.K. Armbruster A.L. Benjamin E.J. Chyou J.Y. Cronin E.M. Deswal A. Eckhardt L.L. Goldberger Z.D. Gopinathannair R. 2023 ACC/AHA/ACCP/HRS Guideline for the Diagnosis and Management of Atrial Fibrillation: A Report of the American College of Cardiology/American Heart Association Joint Committee on Clinical Practice Guidelines Circulation 2024149 e 1e 15610.1161/CIR.000000000000119338033089 PMC 11095842 · doi ↗ · pubmed ↗
- 5Siklódy C.H. Minners J. Allgeier M. Allgeier H.J. Jander N. Keyl C. Weber R. Schiebeling-Römer J. Kalusche D. Arentz T. Pressure-guided cryoballoon isolation of the pulmonary veins for the treatment of paroxysmal atrial fibrillation J. Cardiovasc. Electrophysiol.20102112012510.1111/j.1540-8167.2009.01600.x 19793145 · doi ↗ · pubmed ↗
- 6Packer D.L. Kowal R.C. Wheelan K.R. Irwin J.M. Champagne J. Guerra P.G. Dubuc M. Reddy V. Nelson L. Holcomb R.G. Cryoballoon ablation of pulmonary veins for paroxysmal atrial fibrillation: First results of the North American Arctic Front (STOP AF) pivotal trial J. Am. Coll. Cardiol.2013611713172310.1016/j.jacc.2012.11.06423500312 · doi ↗ · pubmed ↗
- 7Kirchhof P. Benussi S. Kotecha D. Ahlsson A. Atar D. Casadei B. CastelláM. Diener H.C. Heidbuchel H. Hendriks J. 2016 ESC Guidelines for the Management of Atrial Fibrillation Developed in Collaboration with EACTS Rev. Esp. Cardiol. (Engl. Ed.)2017705010.1016/j.rec.2016.11.03328038729 · doi ↗ · pubmed ↗
- 8Andrade J.G. Cryoballoon ablation for pulmonary vein isolation J. Cardiovasc. Electrophysiol.2020312128213510.1111/jce.1445932239557 · doi ↗ · pubmed ↗
