Sarcopenia, systemic immune-inflammation index and all-cause mortality in middle-aged and older people with COPD and asthma: a population-based study
- Elizabeth Benz1,2,
- Sara R.A. Wijnant1,3,4,
- Katerina Trajanoska1,2,
- Johnmary T. Arinze3,5,
- Emmely W. de Roos1,3,
- Maria de Ridder5,
- Ross Williams5,
- Frank van Rooij1,
- Katia M.C. Verhamme5,
- M. Arfan Ikram1,
- Bruno H. Stricker1,
- Fernando Rivadeneira2,7,
- Lies Lahousse1,4,7 and
- Guy G. Brusselle1,3,6,7⇑
- 1Dept of Epidemiology, Erasmus MC – University Medical Center Rotterdam, Rotterdam, The Netherlands
- 2Dept of Internal Medicine, Erasmus MC – University Medical Center Rotterdam, Rotterdam, The Netherlands
- 3Dept of Respiratory Medicine, Ghent University Hospital, Ghent, Belgium
- 4Dept of Bioanalysis, Faculty of Pharmaceutical Sciences, Ghent University, Ghent, Belgium
- 5Dept of Medical Informatics, Erasmus MC – University Medical Center Rotterdam, Rotterdam, The Netherlands
- 6Dept of Respiratory Medicine, Erasmus MC – University Medical Center Rotterdam, Rotterdam, The Netherlands
- 7These authors contributed equally
- Corresponding author: Guy G. Brusselle (guy.brusselle{at}uzgent.be)
Abstract
Background Increasing evidence suggests that sarcopenia and a higher systemic immune-inflammation index (SII) are linked with morbidity in patients with COPD. However, whether these two conditions contribute to all-cause mortality in middle-aged and older patients with COPD or asthma is unclear. Therefore, we investigated the association between sarcopenia, SII, COPD or asthma and all-cause mortality in a large-scale population-based setting.
Methods Between 2009 and 2014, 4482 participants (aged >55 years; 57.3% female) from the population-based Rotterdam Study were included. COPD and asthma patients were diagnosed clinically and based on spirometry. Six study groups were defined according to the presence or absence of COPD or asthma and sarcopenia. Cox regression models were used to assess all-cause mortality in the study groups, adjusted for sex, age, body mass index, SII, smoking, oral corticosteroid use and comorbidities. In addition, all participants were categorised into sex-specific quartiles of SII, and mortality in these groups was compared.
Results Over a median follow-up of 6.1 years (interquartile range 5.0–7.2 years), 466 (10.4%) persons died. Independent of the presence of sarcopenia, participants with COPD had a higher risk of all-cause mortality (hazard ratio (HR) 2.13, 95% CI 1.46–3.12 and HR 1.70, 95% CI 1.32–2.18 for those with and without sarcopenia, respectively). Compared to lower SII levels, higher SII levels increased mortality risk even in people without sarcopenia, COPD or asthma.
Conclusion Middle-aged and older people with COPD, higher SII levels or sarcopenia had an independently increased mortality risk. Our study suggests prognostic usefulness of routinely evaluating sarcopenia and SII in older people with COPD or asthma.
Abstract
Sarcopenia and COPD increased the risk of all-cause mortality in middle-aged and older populations. Additionally, compared to lower SII levels, higher SII levels increased mortality risk in people with and without sarcopenia and COPD. https://bit.ly/3d659m2
Introduction
COPD and asthma are heterogeneous chronic airway conditions in the older population and are highly prevalent worldwide [1]. COPD is the third leading cause of death [2], and older patients with asthma also experience a higher risk of morbidity and mortality [3]. A subgroup of patients with COPD and asthma shares a key feature, systemic inflammation, linked with lung function decline and the pathophysiology of musculoskeletal comorbidities [4–6]. Within the comorbidities framework of these chronic airway diseases, sarcopenia, characterised by the loss of muscle mass and decreased muscle function [7], is a common extrapulmonary condition, particularly in middle-aged and older patients [8]. Depending on diagnostic method, definition criteria and the population setting, a large range of prevalence of sarcopenia in subjects with COPD has been reported (between 8% and 63%) [9]; in older asthmatics from the general population, nearly 5% have reported sarcopenia [8].
In patients with COPD, systemic inflammatory biomarkers have been linked with sarcopenia [5, 10]; however, this effect has not been fully explored in asthmatic patients. Systemic inflammation encompasses a complex network of interactions, including immune-related cells such as neutrophils and lymphocytes [11]. A recent development is that the systemic immune-inflammation index (SII) appears to be a novel marker for chronic inflammation, reflecting increased absolute blood counts of neutrophils and platelets, and decreased counts of blood lymphocytes [12, 13]. Higher SII levels have been associated with worse prognosis and increased mortality in patients with cancer or cardiovascular diseases [14, 15]. Moreover, a higher neutrophil-to-lymphocyte ratio is a marker of systemic inflammation in COPD exacerbations [11, 16]. Subsequently, there is an increased interest in its usefulness in chronic conditions in the respiratory field [17].
Patients with sarcopenia, high SII levels, COPD and, to a lesser extent, asthma have been independently associated with increased risk of all-cause mortality compared to those without sarcopenia, low SII levels, COPD or asthma [11, 18–20]. However, our understanding of whether the combination of COPD or asthma and sarcopenia is related to the risk of all-cause mortality in older people remains limited. Such data might contribute to guiding early interventions as primary prevention.
We hypothesised that middle-aged and older people with COPD or asthma and sarcopenia would have increased risk of all-cause mortality, particularly those with higher SII levels. Therefore, our main objective was to investigate whether these combined conditions of COPD or asthma and sarcopenia increase mortality risk in a middle-aged and older population. Our secondary objective was to explore whether SII levels modulate the effect of sarcopenia and chronic airway diseases on all-cause mortality.
Methods
Study design and participants
This study was conducted within the Rotterdam Study, an ongoing population-based cohort that started in 1989, involving ∼15 000 participants aged ≥45 years, investigating determinants of several chronic diseases [21, 22]. Every 4–5 years, participants from the Ommoord district (in Rotterdam, the Netherlands) undergo follow-up examinations at the research centre (i.e. spirometry) [22].
Our study included all participants with valid data on spirometry, blood sample and sarcopenia, who visited the research centre between 2009 and 2014 (supplementary figure S1). Other relevant data were retrieved from general practitioners’ medical records and prescription data from pharmacies serving the Ommoord area. The Rotterdam Study was approved by the medical ethics committee of the Erasmus Medical Center, Rotterdam (registration number MEC 02.1015), and by the Dutch Ministry of Health, Welfare and Sport (Population Screening Act WBO, license number 1071272-159521-PG). All participants provided written informed consent to participate in the Rotterdam Study and obtain their information from treating physicians.
Exposure definition and measurements
COPD and asthma
Trained paramedic personnel performed pre-bronchodilator spirometry using a MasterScreen PFT Pro spirometer (CareFusion, Houten, the Netherlands) according to the American Thoracic Society and the European Respiratory Society guidelines [23]. Predicted values of forced expiratory volume in 1 s (FEV1) and forced vital capacity (FVC) were calculated employing the reference equations obtained from the Global Lung Function Initiative, which take into account age, sex, height and ethnicity [24]. COPD cases were defined by spirometry data (FEV1/FVC ratio <0.7) and medical records. Asthma was based on physician-diagnosis, as evidenced in the medical records [25]. Asthma–COPD overlap (ACO) cases were identified as asthmatic participants who met the obstructive spirometry definition of COPD. We excluded participants with other lung diseases, ACO, no data on lung function or uninterpretable spirometry. If there was no record for diagnosis of chronic airway diseases or normal spirometry (FEV1/FVC ≥0.7 and FEV1 ≥80%) at the time of sarcopenia measurement, participants were assumed to be without COPD, asthma, ACO or other lung diseases.
Sarcopenia
We defined sarcopenia according to the updated European Working Group of Sarcopenia in Older People (EWGSOP2) criteria (supplementary table S1).
Maximum handgrip strength was obtained from the nondominant hand after three trials using a hydraulic dynamometer (Fabrication Enterprises, White Plains, New York, NY, USA). Appendicular lean mass was measured by dual-energy X-ray absorptiometry (DXA) (iDXA total body-beam densitometer; GE Lunar Corp, Madison, WI, USA). Each scan was analysed using enCORE software V13.6. The sum of the lean mass from the upper and lower limbs is referred as appendicular lean mass (ALM). Appendicular lean mass index (ASMI) was defined as ALM divided by the height squared (m2). Probable sarcopenia was defined as having a low handgrip strength (<27 kg for males and <16 kg for females) [7], and was confirmed as having a low ASMI (<7.0 kg·m−2 for males and <5.5 kg·m−2 for females) [7]. We merged all participants with probable sarcopenia (low handgrip strength) and confirmed sarcopenia (low handgrip strength and low ASMI) as sarcopenia definition.
Study groups
All participants were categorised in six study groups according to the presence/absence of COPD or asthma and sarcopenia: 1) no COPD, asthma or sarcopenia (reference); 2) COPD without sarcopenia; 3) COPD with sarcopenia; 4) asthma without sarcopenia; 5) asthma with sarcopenia; and 6) sarcopenia only (supplementary table S2).
Systemic immune-inflammation index
Fasting blood samples were collected at the study centre. Full blood counts were performed on a Coulter Ac.T diff2 Hematology Analyzer (Beckman Coulter, San Diego, CA, USA). The SII was calculated from the platelets (×109 cells·L−1), granulocytes, as a proxy for neutrophils (×109 cells·L−1), and lymphocytes (×109 cells·L−1), using the formula SII=platelets×neutrophils/lymphocytes [12]. Sex-specific quartiles of SII were defined for all participants.
Outcome and covariates measurements
Information on vital status was obtained from general practitioners and municipal records. Follow-up time started at the date of DXA scan assessment and ended at the date of death or the end of the study (30 May 2018), whichever came first.
For each participant in this study, we retrieved information on the following covariates: age groups (<70 and ≥70 years); sex; smoking status (never, past and current); smoking pack-years (dichotomised at 20 pack-years: <20 pack-years and ≥20 pack-years); body mass index (BMI) (kg·m−2); oral corticosteroid use (Anatomical Therapeutic Chemical code: H02; never: no prescriptions in the complete medical history, past: at least one prescription >3 months before DXA date, current: at least one prescription within the 3 months preceding the DXA date); and physical activity using a validated adapted version of the Longitudinal Aging Study Amsterdam Physical Activity Questionnaire, expressed in metabolic equivalent of task per week.
Moreover, the comorbidity count was defined as none, one or two or more comorbidities, and a fourth category was used for participants with an unknown number of comorbidities (at least one missing). The comorbidities included were type 2 diabetes (fasting plasma glucose level ≥7 mmol·L−1, or a nonfasting plasma glucose level ≥11.1 mmol·L−1 or the use of blood glucose-lowering medication); hypertension (systolic blood pressure ≥140 mmHg or diastolic blood pressure ≥90 mmHg or use of antihypertensive medication); coronary heart disease (myocardial infarction, coronary artery bypass grafting or percutaneous coronary intervention); osteoporosis (T-score ≤−2.5, according to World Health Organization criteria); depression (using the validated Center for Epidemiologic Studies Depression Scale; scores ≥16 were considered as major depressive symptoms). Participants with diagnoses of cancer were excluded.
Statistical analysis
First, we described baseline characteristics for all participants according to the six study groups. Differences in continuous variables between the study groups were tested using one-way ANOVA, and only when the p-value was <0.05, post hoc pairwise comparisons were performed. For skewed continuous variables, the Dunn test was applied. Categorical variables were compared using the Chi-squared test (all p-values are based on nonmissing values).
Second, time to death (survival probability) was described using Kaplan–Meier curves and log-rank or Peto–Peto tests for participants: 1) with and without COPD or asthma; 2) with and without sarcopenia; 3) with different SII quartiles; and 4) in the six study groups.
Third, to evaluate the differences in all-cause mortality between study groups, we performed four Cox regression models, including an increasing number of covariates: model 0 was unadjusted. In model 1, we adjusted for age groups (<70 and ≥70 years) and sex. In model 2 we additionally adjusted for SII quartiles. In model 3, we accounted for previous confounders in model 2, adding BMI, smoking pack-years and oral corticosteroid use. In model 4, we additionally adjusted for the number of comorbidities. Likelihood ratio tests were performed to assess interaction terms between study groups and age groups, sex and SII quartiles into model 4. In addition, we analysed the association between lower versus higher SII quartiles and all-cause mortality in the six study groups. Cox regression models with progressive adjustment for the aforementioned confounders were fitted. We confirmed the assumptions of proportional hazards by statistical evaluation of Schoenfeld residual plots, inspecting for symmetry over time and the p-values.
Lastly, several sensitivity analyses were done. 1) As the vast majority of the sarcopenia cases were participants with low handgrip strength (probable sarcopenia), we repeated the analysis by excluding the cases of low handgrip strength and low ASMI (confirmed sarcopenia). 2) Since physical activity has been associated previously with mortality risk, we adjusted for physical activity. 3) We adjusted for lung function (FEV1 % predicted). 4) Finally, we explored the association between SII quartiles and all-cause mortality in all population and two separate subsets: 1) all groups without sarcopenia and 2) all groups with sarcopenia. Effect estimates were plotted as hazard ratios (HR) for time to death, alongside their corresponding 95% confidence intervals. For 11.3% (n=508) of missing data for physical activity, we used multiple imputation technique (20 imputations with 30 interactions for each imputation), with the mice package of the R software. Statistical analysis was performed using R version 3.5.1 (Foundation for Statistical Computing, Vienna, Austria). A p-value <0.05 was considered to indicate significance in all analyses.
Results
Baseline characteristics of the study population
Among 7162 participants of the Rotterdam Study, 5013 (70.0%) had complete data on spirometry, sarcopenia and SII between 2009 and 2014. Of these, we excluded 509 participants diagnosed with cancer and 22 with ACO, resulting in 4482 (63.0%) participants for the analysis. According to the six study groups, 3074 (68.6%) participants had no COPD nor asthma, nor sarcopenia (reference); 589 (13.1%) had COPD without sarcopenia; 92 (2.1%) had COPD with sarcopenia; 329 (7.3%) had asthma without sarcopenia; 62 (1.4%) had asthma with sarcopenia; and 336 (7.5%) had sarcopenia only (figure 1).
Compared to the reference group, participants of all sarcopenic categories were on average older and more likely to have three or more comorbidities. In addition, they had higher SII levels, except for the group with asthma and sarcopenia (table 1).
Mortality in the study population
The median (IQR) follow-up time was 6.1 (5.0–7.2) years. Over this period, 466 (10.4%) participants died. Of those, 211 (45.3%) had no COPD, asthma or sarcopenia; 99 (21.2%) had COPD without sarcopenia; 36 (7.7%) had COPD with sarcopenia; 21 (4.5%) had asthma without sarcopenia; 34 (7.3%) had asthma with sarcopenia; and 86 (18.5%) had sarcopenia only (supplementary table S2).
The survival probability was worse in subjects with COPD compared to those without COPD and without asthma (figure 2a; log-rank p<0.001) and in subjects with sarcopenia compared to those without sarcopenia (figure 2b; log-rank p<0.001). Moreover, there were significant differences in mortality between SII quartiles (figure 2c; log-rank p<0.001), and participants with higher values of SII (quartile 4) had worse survival probability than those with lower values of SII (quartile 1).
All-cause mortality between study groups
According to the study groups, Kaplan–Meier analyses for survival probability showed significant differences (figure 3; Peto–Peto p<0.001).
Table 2 describes the association of the study groups with all-cause mortality. Compared to the reference group, in the fully adjusted model (model 4), COPD with and without sarcopenia had a higher risk of all-cause mortality (HR 2.13, 95% CI 1.46–3.12 and HR 1.70, 95% CI 1.32–2.18), as well as participants with sarcopenia only (HR 2.09, 95% CI 1.61–2.70). Subjects with asthma and sarcopenia were at increased risk of all-cause mortality, but adjusting for the number of comorbidities attenuated this association (HR 1.75, 95% CI 0.99–3.10) (table 2, figure 4). We found no evidence for age, sex and SII interactions with study groups for all-cause mortality and no deviations from the proportional hazards assumption.
All-cause mortality according to levels of systemic immune-inflammation index
Regarding the association of SII levels and all-cause mortality, in the study groups, subjects with higher levels of SII had an increased risk of mortality than those with lower levels of SII, particularly subjects with sarcopenia only, COPD without sarcopenia and the reference group (model 2; table 3). This association is attenuated by the number of comorbidities (model 3; table 3).
Sensitivity analysis
When the analysis was repeated by excluding participants with confirmed sarcopenia (n=69) (low handgrip strength with normal ASMI), the estimator did not materially change in all study groups (supplementary figures S2, S3 and supplementary table S3). Additional adjustment for physical activity did not affect the effect estimate in all study groups (supplementary table S4), whereas FEV1 (%) attenuated the effect in COPD participants (supplementary table S4). When we focused on the association of SII and all-cause mortality in all participants (n=4482), subjects with higher levels of SII had a significant increase in mortality risk than those with lower levels of SII (HR 1.74, 95% CI 1.33–2.28) (supplementary figure S4a). Similarly, in both subsets of participants, i.e. participants without sarcopenia, COPD or asthma (n=3074), and participants with sarcopenia and COPD or asthma (n=1408), subjects with higher levels of SII had a significantly increased mortality risk compared to those with lower levels of SII (HR 1.50, 95% CI 1.02–2.21 and HR 1.87, 95% CI 1.27–2.75, respectively) (supplementary figure S4b,c). These associations were independent of age groups, sex, BMI, pack-years of smoking, oral corticosteroid use and the number of comorbidities.
Discussion
One of the main results of this population-based study was that middle-aged and older people with COPD (with and without sarcopenia) and sarcopenia only have a significantly higher risk of all-cause mortality than people without COPD, asthma or sarcopenia. In addition, we observed that subjects with higher levels of SII had an increased mortality risk compared to those with lower SII levels. Our results indicate that sarcopenia and high levels of SII are important risk factors for mortality risk in people with and without COPD.
In line with our results, previous studies have reported a significant effect of sarcopenia on mortality in COPD patients [26–30]. Borda et al. [31] reported that sarcopenic COPD patients had a higher risk of all-cause mortality after 3 years of follow-up than those without COPD and sarcopenia. Our results expand these findings from Borda et al. in a larger population (encompassing people with COPD or asthma) over a more extended follow-up period.
On one hand, our results shown that high SII levels increase mortality risk even in people without sarcopenia, COPD or asthma; however, asthma does not increase the risk of all-cause mortality, except for asthmatic participants with sarcopenia, by which mainly the number of comorbidities seems to attenuate the risk. On the other hand, two previous studies from the Rotterdam Study found that high SII levels were associated with risk of cancer and dementia [14, 32]. Importantly, two recent studies have reported that high SII levels were associated with increases all-cause mortality risk in general population [33, 34].
It is pivotal to examine how these biomarkers (i.e. SII) change during the pre-clinical phase of sarcopenia in older patients with chronic airway diseases. One aspect of immune senescence includes increased systemic inflammation, producing different cytokines (T-lymphocytes) linked with pro-inflammatory effects. For instance, elevated levels of tumour necrosis factor-α for a prolonged time are associated with muscle protein degradation at older ages [35, 36]. Moreover, neutrophils release free radicals during the immune system's activation, and muscle cells regulate this response by liberating nitric oxide [35, 37]. It may accumulate damaged protein during ageing due to lifelong exposure to free radicals, slowing muscle protein turnover [35, 38]. Clinical interest is increasing in the use of these haemogram indices, which might predict exacerbation or mortality in COPD. Although SII has been less studied [39], our study showed that levels of SII were increased in COPD and all-sarcopenic groups. Therefore, SII might be a useful index for further investigation and for clinical practice, because it is an additional inflammatory maker to complement the understanding of systemic inflammation pathways. Besides the clinical and population-based studies performed separately for sarcopenia and systemic immune-inflammatory index, further studies on chronic airway diseases are warranted to validate our results.
Strengths of this study include its population-based setting, robust clinical assessment of COPD and asthma (diagnosed by physicians with detailed medical records), high-quality criteria for the diagnosis of sarcopenia based on the EWGSOP2 statements (i.e. DXA scan and handgrip strength) and precise calculation regarding SII (blood sample). Nevertheless, this study has potential limitations that should be taken into account when interpreting these results. First, we had low numbers of events in the combined groups of sarcopenia with COPD or asthma; therefore, it was not always possible to conclude whether SII levels modulate their mortality risk. However, when we analysed two subsets of participants separately: one subset with more-sick people (n=1408, sarcopenia and COPD or asthma) and another subset with less-sick people (n=3074, no sarcopenia and no COPD or asthma), we found that high SII levels had more significant mortality risk than the low SII levels. Second, we did not consider the specific cause of death analyses; however, previous studies have reported that sarcopenia was associated with higher respiratory mortality risk [19]. Third, out of all participants with sarcopenia (n=490), most of them had probable sarcopenia (n=421), and few had confirmed sarcopenia (n=69). When we performed a sensitivity analysis by including only cases with probable sarcopenia (low handgrip strength), the estimator (HR) remained similar in all study groups, which might mean that having only a low handgrip strength increases mortality risk in older people. Nevertheless, more research, which includes only cases with confirmed sarcopenia (low handgrip strength and low ASMI), needs to be performed to quantify the effect on survival, particularly in the older COPD population. Although we believe these results are generalisable to Dutch middle-aged and older people, the applicability of our result to other populations is anticipated, but cannot be definitively addressed by our study. In addition, we observed a high percentage of obesity in all study groups; thus, further research might be focused on this phenotype: sarcopenic-obesity combined with SII levels [40].
In conclusion, our study has shown that COPD, sarcopenia and high SII levels are significantly associated with increased risk of all-cause mortality in middle-aged and older people. These results also reported that the presence of sarcopenia and high levels of SII are independent risk factors for all-cause mortality in middle-aged and older people with or without chronic airway diseases (COPD or asthma). In addition, our analysis suggests benefits in terms of prognostic information when systematically evaluating these independent factors such as sarcopenia and SII index in middle-aged and older patients, particularly those with COPD.
Supplementary material
Supplementary Material
Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.
Supplementary material 00628-2021.SUPPLEMENT
Acknowledgements
The authors are grateful for the dedication, commitment and contribution of the citizens, general practitioners and pharmacists of the Ommoord district who took part in the Rotterdam Study.
Footnotes
Provenance: Submitted article, peer reviewed.
This article has supplementary material available from openres.ersjournals.com
Conflict of interest: E. Benz has nothing to disclose.
Conflict of interest: S.R.A. Wijnant reports grants from GlaxoSmithKline outside the submitted work.
Conflict of interest: K. Trajanoska has nothing to disclose.
Conflict of interest: J.T. Arinze reports a doctoral research grant from Merck Sharp & Dohme, paid to their institution, outside the submitted work.
Conflict of interest: E.W. de Roos has nothing to disclose.
Conflict of interest: M. de Ridder has nothing to disclose.
Conflict of interest: R. Williams has nothing to disclose.
Conflict of interest: F. van Rooij has nothing to disclose.
Conflict of interest: K.M.C. Verhamme reports that they work for a research department that receives or has received unconditional grants from Yamanouchi, Pfizer/Boehringer Ingelheim, Novartis, GSK, UCB, Amgen and Chiesi, none of which are related to the content of this paper.
Conflict of interest: M.A. Ikram has nothing to disclose.
Conflict of interest: B.H. Stricker has nothing to disclose.
Conflict of interest: F. Rivadeneira has nothing to disclose.
Conflict of interest: L. Lahousse reports financial support for the current manuscript from the Fund for Scientific Research Flanders (grant 3G037618), and payment or honoraria for lectures, presentations, speakers’ bureaus, manuscript writing or educational events from Institute for Continuing Study of Pharmacists, outside the submitted work.
Conflict of interest: G.G. Brusselle reports receiving fees for advisory boards or lectures from AstraZeneca, Boehringer Ingelheim, Chiesi, GlaxoSmithKline, Novartis, Sanofi and Teva, outside the submitted work.
Support statement: The Rotterdam Study was supported by the Erasmus Medical Center, Erasmus University Rotterdam, the Netherlands Organisation for Scientific Research (NWO), the Netherlands Organisation for Health Research and Development (ZonMw), the Netherlands Heart Foundation, the Research Institute for Diseases in the Elderly (RIDE), the Netherlands Genomics Initiative (NGI), the Ministry of Education, Culture and Science, the Ministry of Health Welfare and Sports, the European Commission (DG XII), and the Municipality of Rotterdam; and financial support by the Fund for Scientific Research Flanders (grant 3G037618). The funding sources had no involvement in the collection, analysis, writing, interpretation or in the decision to submit the paper for publication. Funding information for this article has been deposited with the Crossref Funder Registry.
- Received November 10, 2021.
- Accepted November 21, 2021.
- Copyright ©The authors 2022
This version is distributed under the terms of the Creative Commons Attribution Non-Commercial Licence 4.0. For commercial reproduction rights and permissions contact permissions{at}ersnet.org