Article Text
Abstract
Objective This study aims to compare aortic morphology between repaired coarctation patients and controls, and to identify aortic morphological risk factors for hypertension and cardiovascular events (CVEs) in coarctation patients.
Methods Repaired coarctation patients with computed tomography angiography (CTA) or magnetic resonance angiography (MRA) were included, followed-up and compared with sex-matched and age-matched controls. Three-dimensional aortic shape was reconstructed using patients’ CTA or MRA, or four-dimensional flow cardiovascular magnetic resonance in controls, and advanced geometrical characteristics were calculated and visualised using statistical shape modelling. In patients, we examined the association of geometrical characteristics with (1) baseline hypertension, using multivariable logistic regression; and (2) cardiovascular events (CVE, composite of aortic complications, coronary artery disease, ventricular arrhythmias, heart failure hospitalisation, stroke, transient ischaemic attacks and cardiovascular death), using multivariable Cox regression. The least absolute shrinkage and selection operator (LASSO) method selected the most informative multivariable model.
Results Sixty-five repaired coarctation patients (23 years (IQR 19–38)) were included, of which 44 (68%) patients were hypertensive at baseline. After a median follow-up of 8.7 years (IQR 4.8–15.4), 27 CVEs occurred in 20 patients. Aortic arch dimensions were smaller in patients compared with controls (diameter p<0.001, wall surface area p=0.026, volume p=0.007). Patients had more aortic arch torsion (p<0.001) and a higher curvature (p<0.001). No geometrical characteristics were associated with hypertension. LASSO selected left ventricular mass, male sex, tortuosity and age for the multivariable model. Left ventricular mass (p=0.014) was independently associated with CVE, and aortic tortuosity showed a trend towards significance (p=0.070).
Conclusion Repaired coarctation patients have a smaller aortic arch and a more tortuous course of the aorta compared with controls. Besides left ventricular mass index, geometrical features might be of importance in long-term risk assessment in coarctation patients.
- Aortic Coarctation
- Cardiac Imaging Techniques
- Aneurysm, Dissecting
- Aortic Aneurysm
- Computed Tomography Angiography
Data availability statement
Data are available upon reasonable request.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
- Aortic Coarctation
- Cardiac Imaging Techniques
- Aneurysm, Dissecting
- Aortic Aneurysm
- Computed Tomography Angiography
WHAT IS ALREADY KNOWN ON THIS TOPIC
Patients with repaired coarctation of the aorta have a lifetime increased risk of developing cardiovascular events and hypertension, even without residual stenosis or re-coarctation.
Previous studies have shown that the distinct aortic geometry of coarctation patients may be associated with the increased risk of cardiovascular events; however, prospective studies evaluating the prognostic consequences of aortic geometry in this population are missing.
WHAT THIS STUDY ADDS
This study describes that repaired coarctation patients have a smaller aortic arch and a more tortuous course of the aorta compared with age-matched and sex-matched controls, stressing that in spite of successful repair, their aortic geometry is not similar compared with a normal aorta.
Aortic tortuosity was identified as a possible additional risk factor besides left ventricular mass and age.
HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY
Findings of this study suggest that aortic geometry deserves more attention in the follow-up of repaired coarctation patients, and especially patients with more aortic tortuosity might require closer surveillance.
Therefore, a large multicentre study is necessary to validate prognostic implications of aortic geometry in repaired coarctation patients.
Introduction
Surgical and percutaneous treatment of aortic coarctations is generally anatomically successfully relieved.1 Despite effective repair, patients’ life expectancy is diminished compared with the general population,1 2 and patients have an increased risk of developing coronary artery disease, aortic complications, arrhythmias, heart failure and cerebrovascular events.3
Up to 57% of the patients have persistent or develop arterial hypertension, even without overt residual stenosis or re-coarctation.1 4 5 Repaired aortas have abnormal vascular function, such as endothelial dysfunction, increased carotid intima-media thickness, reduced arterial compliance, aortic stiffness, reduced baroreceptor sensitivity and dysregulation of the renin–angiotensin system, and this has been associated with surrogate endpoints such as hypertension and increased left ventricular (LV) mass.6 7 Aortic shape has been related to both vascular dysfunction and LV ejection fraction, LV end-diastolic volume, LV mass and hypertension.8–10
We hypothesise that coarctation patients are at increased risk of cardiovascular morbidity and mortality because of their distinct aortic geometry and accompanied vascular dysfunction. This study aims (1) to describe a comprehensive method to analyse aortic geometry using current standard image protocols and state-of-the-art analytical techniques; (2) to compare aortic geometrical characteristics of repaired coarctation patients with healthy controls matched on sex and age; and (3) to investigate associations of geometrical features in repaired coarctation patients with hypertension and cardiovascular events (CVEs).
Methods
Patient population
Patients with repaired coarctation (aged ≥18 years) from the Erasmus Medical Center registered in the prospective follow-up CONCOR registry were included when CT angiography (CTA) or magnetic resonance angiography (MRA) of the complete thoracic aorta within 5 years from inclusion was available (online supplemental figure 1). A post-hoc analysis was performed in patients with available imaging data of their thoracic aorta. Patients were included if they had either surgical or transcatheter repair of the aortic coarctation before their scan. Baseline characteristics including cardiovascular risk factors were collected at the time of inclusion in the registry. Hypertension was defined as a right brachial systolic blood pressure (SBP) >140 mm Hg, diastolic blood pressure (DBP) >90 mm Hg or prescription of antihypertensive medication, and was assessed at baseline and at follow-up.11 LV mass was calculated using the Devereux formula based on echocardiographic reports and was indexed for body surface area. The degree of aortic valve regurgitation or stenosis was based on current guideline recommendations.12 13 Atrial fibrillation was defined as described by treating specialist according to the European Society of Cardiology guidelines.14 Patients were followed up for the occurrence of hypertension, reintervention of the coarctation and CVE from the time of the CTA/MRA until the last patient visit. CVE included aneurysms, dissections, coronary artery disease, ventricular arrhythmia, heart failure hospitalisation, cerebrovascular events, transient ischaemic attack and cardiovascular death. Definitions of CVE are described in the online supplemental appendix. Reinterventions included surgery, stenting, balloon angioplasty or a combination of the latter, aimed to relieve the residual or re-coarctation. CVEs occurring within 30 days from an invasive procedure were considered postprocedural and were excluded. CVEs were scored when reported by the treating physician. Healthy controls, without history of cardiovascular disease, were identified to describe aortic geometry differences between coarctation and the healthy population and were sex and age matched on a group level.
Supplemental material
Image acquisition and segmentation
Scans were acquired between 2003 and 2018, on either Siemens CT scanner or GE MRI 1.5T scanner (figure 1). Follow-up scans were made to monitor aneurysm formation or ventricular function, evaluate residual or re-coarctation or for research purposes. Patients underwent the standard clinical aortic imaging protocol, including either CTA or MRA. Contrast-enhanced four-dimensional (4D) flow MRI was acquired in healthy controls following a standardised study protocol as reported previously.15 Typical spatial resolution characteristics were 0.7×0.7×1.0 mm for CTA, 1.3×1.7×1.4 mm for MRA and 1.8×2.1×1.4 mm for 4D flow. Previously, excellent agreement was found between CTA and MRA-based aortic diameters in our institution.16
Three-dimensional geometry of the thoracic aorta and its branches was extracted from CTA or MRA using the semiautomated segmentation software ITK-snap (ITK-SNAP Medical Image Segmentation Tool) between the aortic annulus and the level of the diaphragm (figure 1). In healthy controls, a similar methodology was followed in CAAS MR Solutions V.5.1 (Pie Medical Imaging, Maastricht, the Netherlands) (figure 1). In Meshmixer (Autodesk, San Rafael, USA), the same aortic volume of interest was extracted, between the anatomical landmarks (annulus aortic valve and diaphragm), in all subjects (figure 1).
Geometry analysis
The global aortic arch geometry was visually categorised by two researchers (SCSM and RvM) into Romanesque, Crenel or Gothic (figure 2).5 The type of coarctation was recorded based on surgical reports of the initial repair. The segmentations were processed using the Vascular Modeling Toolkit (VMTK, Orobix, Bergamo, Italy) (figure 1). The thoracic aorta was divided into the ascending aorta between the aortic valve annulus and the innominate artery, into the aortic arch between the innominate artery and the left subclavian artery and into the descending aorta between the left subclavian artery and the diaphragm. Mean aortic diameter, minimal luminal area, aortic wall surface area and volume of these segments were calculated using an in-house-developed Python script using VMTK functions. Segmentation of the descending aorta was incomplete in the healthy controls; consequently, all measurements of the descending aorta in this group were disregarded. The ratio between ascending and descending aortic diameter and percentage residual stenosis were derived from mean diameter measurements (figure 1A,B). The percentage residual stenosis was defined as 100 × (1−diameterstenosis/diameterreference) in per cent. The average aortic diameter of the descending aorta was taken as the reference diameter. Residual stenosis was considered significant if larger than 20%. Aortic arch was defined as hypoplastic when the ratio of the aortic arch diameter to the reference area in the descending aorta was smaller than 0.9.17
The centreline of each aortic volume was determined and based on this centreline, the curvature, torsion and tortuosity of the centreline were calculated using VMTK functions (figure 1). The centreline is a curve and the deviation of a curve from linearity and planarity is described by the curvature and torsion. Curvature is the degree of bending of the centreline, the deviation from a straight line in one plane and was defined as the inverse of the radius of the osculating circle (figure 1C). The torsion is the twist of the centreline in space, the degree of deviation of the curve from the osculating plane (figure 1D). Curvature and torsion were calculated with the standard VMTK functions (vmtkcenterlinegeometry). Tortuosity was defined as the ratio of the centreline length and the direct distance between the centreline endpoints (figure 1E).
Average shape analysis via statistical shape modelling
In order to visualise the numerical results of our geometrical analysis and fully exploit three-dimensional shape information provided by the reconstructed aortic arches, average aortic shapes of subgroups of interest were calculated via statistical shape modelling (SSM).9 The SSM method is a non-parametric method that analyses how a representative prototype shape (the template) deforms into each of the target shapes that are present in the population.
In short, of the healthy controls, the subject closest to the average volume, aortic wall surface area and centreline was chosen as the initial template. All aortas included for a subgroup were rigidly registered (i.e., aligned) to this initial template using the iterative closest point algorithm provided in VMTK. Average aortic shape (template) per group and its deformations towards each included subject shape were then calculated in the SSM framework Deformetrica (www.deformetrica.org), which used an alternate two-step minimisation algorithm.18 First, the initial template is independently deformed towards each target shape in order to minimise the distance between deformed template and the respective target shapes with respect to the deformations. In the second step, a new, updated template is computed based on these deformations and the initial template in order to minimise the distance between the deformed template shapes (matched shapes) and the original target shapes with respect to the template.19
Once the template for a subgroup was computed, these steps were repeated with the newly calculated template aorta as initial template. This process was repeated until the final template, i.e. the final average aorta’s shape did not change anymore following the principles of Generalised Procrustes Analysis.20 The final template was assumed to be the average (‘prototype’) aortic shape, representative for the considered subgroup. As results are provided as three-dimensional computational surface meshes, this processing allows visualising a subgroup’s average shape in three dimension comprising all relevant shape features instead of merely reporting average morphometric parameters.
Statistical analysis
Normality was visually checked on histograms. For continuous variables, mean and standard deviation (SD) in case of normal distribution, or median and the interquartile range (IQR (25th–75th)) in case of a non-normal distribution were calculated. Categorical variables were summarised with frequencies and percentages. Differences in baseline characteristics and geometrical measurements between patients and healthy controls were tested for significance using an unpaired t-test if normally distributed or a Mann-Whitney test if not normally distributed. To analyse differences of the original morphology, characteristics of the initial repair and geometrical characteristics between subgroups (descriptive morphology subtypes, tricuspid aortic valves (TAVs) vs bicuspid aortic valves (BAVs), surgical technique, no vs aortic hypoplasia), differences were tested with one-way analysis of variance (ANOVA) (when normally distributed) or Kruskal-Wallis test (when not normally distributed) and subsequently pairwise comparisons were performed between separate groups with an unpaired t-test or a Mann-Whitney test. For categorical variables, Χ2 test or Fisher’s exact test was performed, as appropriate. To determine the optimal cut-off value for the risk of hypertension and CVE, the Youden index was estimated at each percentage of residual stenosis and minimal luminal area of descending aorta and maximised.
To assess repeatability of our measurements, intraobserver variability of the geometrical measures was assessed in 15 randomly selected aortas (five CTA, five MRA and five 4D flow MRI). Repeatability data were reported using intraclass correlation coefficients (ICCs). The ICC was calculated using a two-way ANOVA model to determine observer variability based on a single measure and absolute agreement.
Univariable logistic regressions were used to investigate the association between the type of coarctation and presence of hypoplasia, characteristics of the initial repair (type of surgery and age at initial surgery) and baseline and geometrical characteristics with hypertension at the baseline moment. Univariable Cox proportional hazard regressions were used to investigate determinants of CVE and reintervention. Patients were followed from the time of the scan until the occurrence of a CVE and reintervention. If no CVE or reintervention occurred, they were censored at the time of the last clinic visit. For these analyses, residual stenosis was entered as a continuous variable. To obtain the most informative multivariable logistic regression and Cox proportional hazards models, while also considering the limited number of events per variable, we applied least absolute shrinkage and selection operator (LASSO) penalised regression with the penalty parameter selected by 10-fold cross-validation. Besides the geometrical risk factors, other risk factors were included based on their clinical relevance and previous reports; these included sex, age at inclusion, age at repair, aortic valve type, type of surgery and LV mass indexed.2
Statistical analyses were conducted using R Statistical Software V.4.1.0 (R Foundation for Statistical Computing, Vienna, Austria). Two-tailed p values of <0.05 were considered statistically significant.
Results
Study population
In total, 65 patients were eligible for inclusion. Baseline variables of patients and healthy controls are shown in table 1. Fourteen patients had aortic arch hypoplasia, of which the majority (10 out of 14) had proximal aortic arch hypoplasia. Other than SBP, no differences were statistically significant. Time between inclusion and scan was 1.6 years (IQR 0.5–3.0). During the follow-up between CTA/MRA until censoring (median 8.7 years (IQR 4.8–15.4)), 27 CVEs occurred in 20 patients. Four patients died because of a cardiovascular reason, and aortic aneurysm was the most frequently occurring event (22%, 6 of 27 times) (online supplemental table 1). Five of the six aneurysms occurred in the ascending aortic in patients with a BAV and not at the site of the surgical repair. Four (6%) patients went for reintervention during follow-up at a median follow-up period of 1.2 years (IQR 0.6–1.7). Given the low number of reinterventions, no further univariable and multivariable analyses were performed on risk factors of reinterventions. At the last clinical visit, 45 (69%) patients had hypertension or antihypertensive medication (of which 40 used antihypertensive medication) and the median blood pressure was 132/79 mm Hg. Comparing baseline versus follow-up, SBP decreased (p=0.045), DBP was similar (p=0.357), antihypertensive medication increased (p=0.036) and overall patients with hypertension was similar (p=0.853). Patients were matched with 23 healthy controls.
Aortic geometry
An overview of the aortic geometry parameters in healthy controls and coarctation patients is presented in table 2 and is visually reflected in figure 2A. Ascending aortic diameters, wall surface area and volume were not different in coarctation patients with a TAV compared with healthy controls, but were higher in coarctation patients with a BAV compared with healthy controls (all p<0.001) (table 2). Aortic arch diameters, wall surface area and volumes were smaller in coarctation patients compared with healthy controls (diameter p<0.001, wall surface area p=0.026 and volume p=0.007). Patients had more torsion of their thoracic aorta (entire thoracic aorta patients 30 dm−1 (IQR 27–33), healthy controls 22 dm−1 (IQR 20–25; p<0.001); aortic arch patients 17 dm−1 (IQR 11–21), healthy controls 10 dm−1 (IQR 8–13; p<0.001)) (table 2 and figure 2). The aortic arch curvature was significantly higher in patients (4.2 dm−1 (IQR 3.8–4.9)) compared with healthy controls (3.7 dm−1 (IQR 3.4–4.0); p<0.001).
Differences stratified per descriptive aortic morphology are presented in online supplemental table 2 and figure 2B. Patients with a Crenel-shaped aorta were significantly older at their initial intervention (p=0.046) and had more commonly a hypoplastic aortic arch (p=0.007) compared with patients with a Romanesque-shaped aorta. LV mass indexed was higher in patients with a Gothic-shaped aorta compared with Romanesque-shaped aortas (p=0.023). Aortic torsion was higher in patients with a Gothic-shaped aorta compared with Crenel (p=0.008) or Romanesque-shaped aortas (p=0.022), and aortic arch torsion was lower in patients with a Romanesque-shaped aorta (vs Crenel p=0.027, vs Gothic p=0.010). Aortic arch hypoplasia was not significantly different in geometry groups. No other clinical or geometrical parameters were associated with a certain aortic arch morphology.
Besides larger aortic arch wall surface and volume in patients with BAV, there were no geometrical differences between patients with a BAV and a TAV (online supplemental table 3). Online supplemental table 4 shows the analysis of aortic shape features of the three most common surgical techniques in our cohort. Patients with a patch angioplasty were significantly older, had a higher LV mass, more residual stenosis and more torsion of their aortic arch compared with end-to-end anastomosis patients. Finally, the aortic shape was compared of patients with and without hypoplastic aortic arches (online supplemental table 5). Patients with hypoplastic aortic arches were more often hypertensive and had significantly more tortuosity and a lower mean aortic arch curvature compared with those without aortic arch hypoplasia.
Intraobserver agreement of the aortic shape metrics was in general good to excellent (online supplemental table 6), except for the aortic arch volume (ICC 0.789 (95% CI 0.493 to 0.923)) and aortic torsion (ICC 0.530 (95% 0.057 to 0.812)) (online supplemental table 6).
Risk factors for hypertension
In univariable logistic regression, male sex, age, age at repair, LV mass indexed and aortic arch hypoplasia were positively associated with hypertension (male sex OR 1.33 (95% CI 1.07 to 1.66), p=0.013; age per 5-year increase OR 1.05 (95% CI 1.01 to 1.10), p=0.025; age at repair per 5 years OR 1.11 (95% CI 1.05 to 1.17), p<0.001; LV mass indexed per 25 g/m2 (LV mass) OR 1.08 (95% CI 1.01 to 1.16); and hypoplastic aortic arch 1.29 (95% CI 1.00 to 1.66), p=0.023) (online supplemental table 7). After LASSO regression selection, in the most informative multivariable model, male sex, age per 5 years and tortuosity were associated with hypertension. Based on the highest Youden index, the optimal cut-off value for a residual stenosis was 25% to distinguish patients with and without hypertension (Youden index 0.20). For the minimal luminal area, 134 mm2 was the optimal cut-off value.
Risk factors of CVEs
Figure 3 shows the average aortic shape of patients with and without CVE. Patients without CVE had a Romanesque-shaped aorta more similar to the aortic shape of healthy controls, while patients with CVE had a more angulated aortic arch and a more tortuous course of the descending aorta. Based on the highest Youden index, the optimal cut-off value for a residual stenosis was 15% to distinguish patients with and without CVE (Youden index 0.11). Online supplemental table 8 shows regression analysis of risk factors for CVE. In the univariable analysis of the geometrical features, aortic tortuosity and minimal luminal area of the descending aorta were associated with increased risk of CVE (aortic tortuosity HR 2.20 (95% CI 1.18 to 4.10), p=0.013; and minimal luminal area of the descending aorta HR 1.00 (95% CI 1.00 to 1.01), p=0.035). Residual stenosis was not associated with CVE. Using LASSO, age, LV mass and aortic tortuosity were selected for the most informative multivariable model. In this model, LV mass was independently associated with CVE (HR 1.38 (95% CI 1.07 to 1.78), p=0.014). There was a trend towards an association between aortic tortuosity and CVE; however, this was not significant (HR for aortic tortuosity 2.20 (95% CI 0.94 to 5.17), p=0.070) (online supplemental table 8).
Discussion
The current study shows that comprehensive analysis of three-dimensional geometrical features of the aorta is feasible using standard imaging protocols and state-of-the-art analysis techniques such as SSM. Geometrical aortic characteristics of coarctation patients (with and without CVE) and healthy age-matched controls were quantitatively compared and visually summarised via computed three-dimensional average shapes. Residual stenosis was not associated with the occurrence of CVE, while there was a trend towards an association between aortic tortuosity and CVE.
The presented methodology provides a complete overview of the aortic geometry in coarctation patients by looking beyond simple aortic dimensions, and investigate for example, cross-sectional vessel area, wall surface area and volumes. With three-dimensional metrics, such as the curvature and torsion, the course of the aorta through the thorax is quantified. This study demonstrates that the aorta of coarctation patients has a distinct shape and a more complex course compared with the healthy population with more torsion, a higher curvature and a smaller aortic arch. These differences underline that coarctation of the aorta is not just a ‘narrowing’ at the aortic isthmus, but an aortopathy affecting a larger region of the aorta.
No studies have reported the relationship between complex aortic shape features and CVE. In other patient groups, i.e., patients with Marfan syndrome, ascending aortic curvature was associated with thoracic aneurysm growth.21 Previous studies in coarctation patients associated a more curved, Gothic-shaped aortic arch with surrogate endpoints such as cardiac parameters (e.g., lower LV ejection fraction) and vascular parameters (e.g., abnormal intima thickness).8 9 In this study, although aortic arch curvature appeared to be higher in patients with CVE, there was no significant association between curvature and/or Gothic-shaped aorta and CVE. As in a previous study in patients with a type B aortic dissection, tortuosity was associated with CVE in this study which might be a stronger predictor than curvature and/or Gothic-shaped aorta.22
Interplay between aortic shape, flow and pathophysiology has long been recognised both in embryology, vessel homeostasis and in the onset of disease.23 24 In coronary arteries, vessel torsion specifically has shown to promote helical flow.25 Helical flow is associated with higher wall shear stress values and is considered atheroprotective, but also with plaque rupture and aneurysm formation.15 26–29 Pathological wave reflections were previously associated with aortic radius and not with curvature.30 In our population, potentially altered haemodynamics, because of the altered vessel shape, impact the aortic wall by inducing vessel wall remodelling and increasing LV afterload. These two processes are associated with the occurrence of CVE. However, the assessment of the direct relationship between altered flow haemodynamics and the risk on CVE was beyond the scope of this research.
For the risk of hypertension, a cut-off for a residual stenosis of 25% narrowing of the aortic diameter seems to be appropriate and the minimal luminal area in the descending aorta was not better to predict hypertension than the degree of residual stenosis. The low Youden index of both parameters supports the theory that other factors (e.g., vascular remodelling) play an important role in the development of hypertension.8 Clinical parameters such as male sex and age at repair appeared to have stronger associations with hypertension than residual stenosis and arch curvature. The degree of residual stenosis is currently central in the treatment of repaired coarctation patients, is associated with LV afterload and aortic wall shear stresses proximal to the coarctation and was expected to predict the risk of CVE.31 32 However, the present study did not show an association between the degree of residual stenosis and the risk on CVE. The low Youden index of residual stenosis also stresses that the predictive value of residual stenosis for CVE is not strong. Future larger studies are required to confirm this finding.
Study limitations
For aortic shape reconstructions, different imaging modalities were used for coarctation patients and controls which might have impacted the aortic dimension measurements and hindered the evaluation of features involving the distal descending aorta in controls. Only patients with a CT or MRI available acquired in a certain time frame were selected and this might have resulted in inclusion of a selected population; the results might not be applicable to all postoperative coarctation patients. Despite LASSO penalised regression, there is a possibility that our multivariable models are influenced by characteristics unique to our small sample size. Also, our study population was rather small and the total number of events was small. Therefore, our findings should be confirmed in larger studies. Finally, ambulatory blood pressure measurements were not available in all patients; therefore, office blood pressure measurements were used.
Conclusion
After successful coarctation repair, patients have a distinct aortic arch geometry compared with matched controls with a more tortuous course and smaller aortic arch. The degree of residual stenosis can only partly explain the presence of hypertension in patients and was not associated with CVE. However, to predict CVE, other geometrical features (e.g., tortuosity) might be helpful. Gaining insight into which factors attribute to the altered aortic morphology in coarctation patients may be beneficial in reducing lifelong cardiovascular burden in these patients and improving their quality of life.
Data availability statement
Data are available upon reasonable request.
Ethics statements
Patient consent for publication
Ethics approval
This study involves human participants and was approved by the Erasmus Medical Center Ethical Committee (MEC-2019-0664). All participating patients provided written informed consent.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
SCSM and RvM contributed equally.
Contributors Conceptualisation—SCSM, RvM, TAM, MV and AEvdB. Data curation—SCSM, RvM and TAM. Data processing—S-AK. Formal analysis—SCSM and RvM. Supervision—AEvdB. Writing (original draft)—SCSM, RvM and AEvdB. Writing (review and editing)—all authors. Guarantor—SCSM.
Funding This research was funded by the Thorax Foundation.
Competing interests None declared.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.