Article Text
Abstract
Introduction Although inflammation assessed by elevated C reactive protein (CRP) concentration is known to be associated with risk of cardiovascular disease, its association with microvascular and macrovascular dysfunction in diabetes and non-diabetes remains unclear. We examined the association between CRP and diabetes and associated microvascular and macrovascular dysfunction in sub-Saharan Africans with and without diabetes.
Research design and methods Cross-sectional analyses of baseline data from the multicenter RODAM study (Research on Obesity and Diabetes among African Migrants) including 5248 Ghanaians (583 with diabetes, 4665 without diabetes) aged 25–70 years were done. Logistic regression analyses were used to examine the associations between CRP Z-scores and diabetes and microvascular (nephropathy) and macrovascular (peripheral artery disease (PAD)) dysfunction, with adjustments for age, sex, site of residence, smoking, body mass index, systolic blood pressure, and low-density lipoprotein cholesterol.
Results In the fully adjusted models, higher CRP concentration was significantly associated with diabetes (adjusted OR 1.13; 95% CI 1.05 to 1.21, p=0.002). In participants with diabetes, higher CRP concentration was associated with PAD (1.19; 1.03 to 1.41, p=0.046) but not nephropathy (1.13; 0.97 to 1.31, p=0.120). Among participants without diabetes, higher CRP concentration was associated with higher odds of PAD (1.10; 1.01 to 1.21, p=0.029) and nephropathy (1.12; 1.04 to 1.22, p=0.004).
Conclusions In this study, higher CRP concentration was associated with higher odds of diabetes in sub-Saharan Africans. Also, higher CRP concentration was associated with higher odds of nephropathy and PAD in non-diabetes and higher odds of PAD in diabetes. CRP may be an important marker for assessment of risk of diabetes and risk for PAD and nephropathy in sub-Saharan Africans with and without diabetes.
- C-reactive protein
- inflammation
- microvascular and macrovascular complications
- adult diabetes
This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.
Statistics from Altmetric.com
Significance of this study
What is already known about this subject?
In European and Asian populations, inflammation measured by elevated C reactive protein (CRP) is known to be associated with diabetes.
Inflammation is a known factor in the development of atherosclerosis and subsequent atherosclerotic vascular events.
What are the new findings?
In sub-Saharan Africans, higher CRP concentration is associated with higher odds of diabetes, even after adjustments for the conventional cardiovascular risk factors.
Higher CRP concentration is significantly associated with higher odds of peripheral artery disease (PAD) and nephropathy in non-diabetes and higher odds of PAD in diabetes.
The conventional cardiometabolic risk factors did not explain the associations between CRP and PAD and nephropathy in diabetes and non-diabetes.
How might these results change the focus of research or clinical practice?
Our findings suggest that CRP could be explored as a potential marker to identify sub-Saharan Africans at risk of diabetes, PAD and nephropathy.
Introduction
Globally, microvascular and macrovascular diseases are important public health problems.1–4 In many regions of the world including sub-Saharan Africa, the rates of microvascular and macrovascular diseases are rising, contributing to the increasing rates of disability and death from cardiovascular disease (CVD).5 6 Specifically, macrovascular diseases including peripheral artery disease (PAD), coronary artery disease, and cerebrovascular disease may complicate critical limb ischemia, myocardial infarction, and stroke, respectively.7 Also, microvascular diseases such as retinopathy, nephropathy, and neuropathy may result in blindness, end-stage kidney disease, and lower limb amputation, respectively.8
C reactive protein (CRP), the most extensively studied biomarker of inflammation, is known to be significantly associated with CVD including diabetes in European and Asian populations; however, data on its role in CVD in sub-Saharan African populations are limited.9 10 Considering the substantial ethnic differences in the association between inflammation and diabetes, it may be valuable to investigate this association in sub-Saharan African populations, in the quest to integrate CRP to global risk scores for diabetes.11
Existing data on the association between inflammation and vascular dysfunction have focused on individuals with diabetes.12–14 In diabetes, endothelial injury from inflammation mediated by chronic hyperglycemia is known to play a key role in the development of vascular complications.4 15 However, the effect of glycemic control on macrovascular complication risk or progression in diabetes remains uncertain.13 16 Also, existing data suggest important associations between inflammation and CVD in non-diabetes.17 Given the above, it is plausible that inflammation triggered by causes other than hyperglycemia may be important in the pathogenesis of microvascular and macrovascular dysfunction; however, the biological basis for this association has not been clarified.17 We, therefore, assessed the associations between CRP and diabetes in Ghanaians. In addition, we assessed associations between CRP and microvascular (nephropathy) and macrovascular dysfunction (PAD) in Ghanaians with and without diabetes.
Materials and methods
Study design
The rationale, conceptual framework, design and methodology of the RODAM study (Research on Obesity and Diabetes among African Migrants) have been described in detail elsewhere.18 In brief, the study was conducted from 2012 to 2015 and comprised Ghanaians aged 25–70 years living in rural and urban Ghana as well as in three European cities, namely Amsterdam, Berlin and London. Data collection for the study was standardized across all sites. Informed written consent was also obtained from each participant prior to enrollment in the study. For the current analyses, only participants with complete data on CRP, microvascular disease (nephropathy) and macrovascular disease (PAD) were included. This comprised data from 583 participants with diabetes and 4665 participants without diabetes aged 25–70 years (figure 1).
Assessments
Questionnaires
A structured questionnaire18 was used to record the demographic, socioeconomic, and health-related behaviors of the study participants. Smoking was assessed as a positive reply to the question ‘Do you smoke at all?’ Alcohol intake in grams per day was estimated using standard portion sizes combined with frequencies of intake based on a standardized Food Propensity Questionnaire.19 Physical activity was derived for each participant using the International Physical Activity Questionnaire, and participants were classified into three categories of total physical activity, namely low, moderate and high level.20
Physical examination
According to standard operation procedures, weight was measured in light clothing and without shoes with SECA 877 scales to the nearest 0.1 kg. Height was measured without shoes with SECA 217 stadiometer to the nearest 0.1 cm. Body mass index (BMI) was calculated as weight (kg) divided by height squared (m2). Blood pressure (BP) was measured three times using the Microlife WatchBP home device, with appropriately sized cuffs after at least 5 min rest while seated. The mean of the last two BP measurements was used for the analyses. Hypertension was defined as systolic BP ≥140 mm Hg and/or diastolic BP ≥90 mm Hg, and/or being on antihypertensive medication treatment. All the anthropometrics were measured twice by the same assessor and the average of the two measurements was used for analyses.
Ankle-brachial pressure index (ABI) measurements were performed in supine position using a validated oscillometric device (WatchBP Office ABI, Microlife, Widnau) with appropriately sized cuffs after at least 10 min of supine rest.21 Systolic BP was measured twice in the right and left brachial artery and twice in the right and left posterior tibial arteries. ABI was calculated by taking the highest arm systolic BP as the denominator and the lowest ankle BP as the numerator.22 The lowest of the left and right ABI measurements was used for analyses. ABI obtained by the oscillometric method using the Microlife WatchBP Office ABI has been shown to correlate well with ABI acquired by Doppler ultrasound with 95% agreement between the two methods in diagnosing PAD.23
Biochemical analyses
According to standard operation procedures, fasting venous blood samples were processed and aliquoted into Sarstedt tubes after collection, registered and then temporarily stored at the research site at −20°C. Aliquoted blood samples and first early morning urine sample were transported to the local research centers, where they were checked, registered and stored at −80°C before being shipped to the central laboratory at Charité-University Medicine Berlin (Berlin, Germany) for determination of biochemical variables. Shipping of the samples from European sites was carried out using Styrofoam boxes filled with dry ice and from Ghana in dry shippers filled with liquid nitrogen. Extensive quality checks were done during the biochemical analysis, including blinded serial measurements. Fasting glucose, total cholesterol, low-density lipoprotein (LDL) cholesterol, high-density lipoprotein cholesterol, and triglycerides levels were determined using the ABX Pentra 400 Chemistry Analyzer (HORIBA ABX, Montpellier, France). Fasting plasma glucose concentration was measured using an enzymatic method (hexokinase). The concentration of total cholesterol was assessed using colorimetric test kits. Glycosylated hemoglobin (HbA1c) was measured by high-performance liquid chromatography (Tosoh G8 HPLC Analyzer). The concentration of urinary albumin (in µmol/L) was measured by an immunochemical turbidimetric method (Roche Diagnostics) and urinary creatinine concentration (in μmol/L) was measured by a kinetic spectrophotometric method (Roche Diagnostics). Urinary albumin to creatinine ratio (ACR; expressed in mg/mmol) was calculated by taking the ratio between urinary albumin and urinary creatinine. Serum creatinine concentration was determined by a kinetic colorimetric spectrophotometric isotope dilution mass spectrometry calibration method (Roche Diagnostics). The estimated glomerular filtration rate (eGFR) was calculated using the 2009 CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration) creatinine equation. High-sensitive CRP (hs-CRP) levels were measured in heparin plasma by a particle enhanced immunoturbidimetric assay. Human CRP agglutinates with latex particles were coated with monoclonal anti-CRP antibodies. The aggregates were determined turbidimetrically using Pentra 400 Chemistry Analyzer (HORIBA ABX). Z-scores were computed for each participant based on the participant’s CRP concentration (x), mean CRP concentration (μ) and SD (σ). We computed the Z-scores as: Z-score=(x−μ)/σ.
Definition of diabetes and microvascular and microvascular dysfunction
Diabetes was defined according to the WHO diagnostic criteria (self-reported diabetes, documented use of glucose-lowering medication(s), fasting plasma glucose ≥7.0 mmol/L) or HbA1c ≥6.5% or ≥48 mmol/mol.24 PAD was defined as ABI ≤0.90.22 ABI >1.3 was not considered normal as it could be suggestive of non-compressible vessels.22 Albuminuria and eGFR were classified according to the 2012 Kidney Disease: Improving Global Outcomes (KDIGO) guidelines.25 eGFR was classified as follows: G1, 90 mL/min/1.73 m2 (normal kidney function); G2, 60–89 mL/min/1.73 m2 (mildly decreased); G3a, 45–59 mL/min/1.73 m2 (mildly to moderately decreased); G3b, 30–44 mL/min/1.73 m2 (moderately to severely decreased); G4, 15–29 mL/min/1.73 m2 (severely decreased); and G5, <15 mL/min/1.73 m2 (kidney failure). Low eGFR was defined as eGFR <60 mL/min/1.73 m2 (category ≥G3).26 Albuminuria classifications were derived from ACR and were as follows: A1, <3 mg/mmol (normal to mildly increased); A2, 3–30 mg/mmol (moderately increased); and A3, >30 mg/mmol (severely increased). Because of the small number of participants in the severely increased albuminuria category (A3), we combined the moderately increased (A2) and severely increased (A3) categories.25 26 Nephropathy was defined as albuminuria and/or eGFR <60 mL/min/1.73 m2 based on the KDIGO guidelines.25 27
Statistical analysis
Data were analyzed using the IBM SPSS V.23 for Windows. Data with normal distributions were presented as mean±SD, whereas those not normally distributed presented as median (IQR). Categorical data were presented as frequencies (percentages). No significant interaction between site of residence and CRP was found; therefore, we did not stratify the analyses by site of residence. However, we used site as a covariate in our logistic regression models. Logistic regression analyses were used to examine the associations between CRP concentrations and diabetes, PAD and nephropathy with adjustments for covariates. ORs and their corresponding 95% CIs were estimated. The minimal sufficient adjustment sets for estimating the direct effect of CRP on PAD and nephropathy were determined by a directed acyclic graph (DAG) (online supplementary figure 1; DAG available at dagitty.net/mGETVa3). DAG was chosen because the traditional methods of adjusting for potential confounders may introduce conditional associations and bias rather than minimize it.28 Based on the DAG, the minimal sufficient adjustment sets for estimating the total effect of CRP on vascular dysfunction were age, smoking, obesity, hypertension, diabetes and dyslipidemia. Three models were used to examine the data. Model 1 was unadjusted, model 2 was adjusted for age and sex, and in model 3 we additionally adjusted for site of residence, smoking, BMI, systolic BP, diabetes and LDL-cholesterol. In the diabetes and non-diabetes subgroup analyses, we adjusted for age, sex, site of residence, smoking, BMI, systolic BP, and LDL-cholesterol in model 3.
Supplemental material
Results
Baseline characteristics
Table 1 shows the baseline characteristics of the study participants. Compared with participants without diabetes, participants with diabetes were less frequently female, less physically active, much older, consumed more alcohol and had a lower proportion of people who have never smoked. Also, the mean BMI, waist circumference, hip circumference, waist to hip ratio, systolic BP, diastolic BP, total cholesterol and LDL-cholesterol were higher in participants with diabetes than their peers without diabetes. The median CRP concentration and proportion of participants with CRP >3 mg/L were higher in participants with diabetes than those without diabetes. Participants with diabetes had higher rates of PAD and nephropathy than their peers without diabetes.
Association between CRP and diabetes
Higher CRP Z-score concentration was associated with higher odds of diabetes in the unadjusted model (figure 2). The association persisted in a fully adjusted model (adjusted OR (AOR) 1.13; 95% CI 1.05 to 1.21, p=0.002). In a sensitivity analysis, we excluded participants with CRP concentrations >10 mg/L; the SD of the mean CRP was 8.20 mg/L for the whole group and 2.07 mg/L for the subgroup with CRP concentrations ≤10 mg/L. In participants with CRP concentrations ≤10 mg/L (n=4962), higher CRP Z-score concentration was significantly associated with higher odds of diabetes in both the unadjusted (OR 1.31; 95% CI 1.22 to 1.42, p<0.001) and fully adjusted (AOR 1.24; 95% CI 1.14 to 1.36, p<0.001) models.
Association between CRP and microvascular and macrovascular dysfunction
In participants without diabetes, higher CRP Z-score concentration was associated with higher odds of PAD and nephropathy (table 2). The differences persisted in the fully adjusted model: PAD (AOR 1.10; 95% CI 1.01 to 1.21, p=0.029) and nephropathy (AOR 1.12; 95% CI 1.04 to 1.22, p=0.004). Among participants with diabetes, higher CRP Z-score concentration was also associated with PAD and nephropathy. The associations remained statistically significant in the fully adjusted models for PAD (AOR 1.19; 95% CI 1.03 to 1.41, p=0.046), but not for nephropathy (AOR 1.13; 95% CI 0.97 to 1.31, p=0.120).
In the analysis of the association between higher CRP Z-score concentration and the individual components of nephropathy, higher CRP Z-score concentration was associated with higher odds of albuminuria and low eGFR in the non-diabetes group. The differences persisted in the fully adjusted model for albuminuria (AOR 1.15; 95% CI 1.05 to 1.24, p=0.001), but not for low eGFR (AOR 1.07; 95% CI 0.94 to 1.21, p=0.316). Among participants with diabetes, higher CRP Z-score concentration was associated with higher odds of albuminuria, but not with low eGFR in the unadjusted model. In the fully adjusted model the association between higher CRP Z-score concentration and albuminuria was no longer statistically significant (AOR 1.16; 95% CI 0.99 to 1.36, p=0.066).
The sensitivity analysis showing the associations between CRP Z-score concentration and PAD and nephropathy in participants with CRP concentrations ≤10 mg/L is shown in table 3. In the non-diabetes group with CRP ≤10 mg/L, higher CRP Z-score concentrations were significantly associated with higher odds of PAD and nephropathy. The associations were no longer statistically significant in the fully adjusted model: PAD (AOR 1.06; 95% CI 0.93 to 1.21, p=0.372) and nephropathy (AOR 1.32; 95% CI 0.87 to 1.99, p=0.194). In participants with diabetes with CRP concentrations ≤10 mg/L, higher CRP Z-score concentration was not significantly associated with PAD or nephropathy: PAD (AOR 1.05; 95% CI 0.75 to 1.47, p=0.785) and nephropathy (AOR 1.69; 95% CI 0.77 to 3.71, p=0.189).
Discussion
Key findings
Using a sample of Ghanaians, we showed that among sub-Saharan Africans higher CRP concentration is associated with higher odds of diabetes, even after adjustments for the conventional cardiovascular risk factors. Additionally, higher CRP concentration is significantly associated with higher odds of nephropathy and PAD in non-diabetes and higher odds of PAD in diabetes. The conventional cardiometabolic risk factors did not explain the associations between CRP and PAD and/or nephropathy in diabetes and non-diabetes.
Discussion of key findings
This is the first study assessing the relationship between CRP and diabetes among a sub-Saharan African population. The results of our study in Ghanaians suggest an association between higher CRP and higher odds of diabetes in sub-Saharan Africans. A systematic review and meta-analysis by Wang et al10 had previously shown that elevated levels of inflammatory markers were significantly associated with increased risk of diabetes. The data on which this conclusion was based excluded sub-Saharan Africa origin populations. Considering the ethnic differences in the association between inflammation and diabetes, data on this association in sub-Saharan Africans are important in the quest to integrate CRP to global risk scores for diabetes.29 Our finding of a positive association between CRP and diabetes has confirmed earlier findings and expanded the evidence to sub-Saharan African populations. Considering the persistence of a statistically significant association between higher CRP concentration and higher odds of diabetes after excluding individuals with CRP concentrations >10 mg/L, our results could reflect the role of subclinical inflammation, instead of acute infection. Our cross-sectional analyses also suggest that a reciprocal association between inflammation and diabetes could exist as diabetes may lead to impaired immunity and increase the risk of chronic infections and/or infestations with subsequent rise in the levels of inflammatory markers.30
Although inflammation measured by CRP is a known risk factor for CVD, its associations with microvascular and macrovascular dysfunction in diabetes and in non-diabetes are unclear. In diabetes, most studies investigating this relationship have reported significant positive association between inflammation and vascular complications.12 31 32 Our results show that in sub-Saharan Africans with diabetes, higher CRP concentration is positively associated with PAD. The associations persisted after adjustments for the conventional cardiovascular risk factors. These findings resemble results previously reported in other populations.31 32 In this current study, we did not find a statistically significant association between CRP and nephropathy in participants with diabetes. Existing data in other ethnic groups suggest that higher CRP concentration is associated with nephropathy in diabetes.33 The direction of the association between CRP and nephropathy in previous studies is the same in the current study; therefore, it is possible that limited power restricted us from finding a significant association between CRP and nephropathy in diabetes.
There are limited data on the association between CRP and microvascular/macrovascular dysfunction in individuals without diabetes. Our study shows that in sub-Saharan Africans without diabetes, higher CRP concentration is significantly associated with PAD and nephropathy. This finding highlights the role of inflammation in the pathogenesis of microvascular/macrovascular dysfunction, independent of hyperglycemia.17 34 This has important implications in low-to-middle-income regions like sub-Saharan Africa, where in addition to the growing burden of obesity and diabetes there is a high burden of proinflammatory conditions, including chronic or recurrent infections and infestations.6 11 35 36 Chronic intravascular infections and infestations may trigger the inflammatory pathways, predisposing the macrovasculature and microvasculature to accelerated atherosclerosis.17 37 Also, chronic extravascular infections can enhance extravascular production of inflammatory cytokines that may accelerate the process of atherosclerosis.17 The conventional cardiometabolic risk factors did not significantly attenuate the association between CRP and PAD or nephropathy in non-diabetes. This further highlights the fact that the inflammatory pathways driving microvascular and macrovascular dysfunction may not be fully dependent on the conventional cardiovascular risk factors. Indeed, we have previously reported that the higher rates of microvascular and macrovascular dysfunction in sub-Saharan Africans living in Africa compared with their peers in Europe are not fully explained by glycemic control or the conventional cardiovascular risk factors.38 39
Overall, our results support the idea that some process related to persistent inflammation is associated with microvascular and macrovascular dysfunction, although the clinical utility of this remains unclear, especially in the setting of low-income to middle-income countries where the likelihood of low-grade inflammation from non-hyperglycemic causes including chronic infections and infestations is high.36 Based on the CANTOS trial (Canakinumab Anti-Inflammatory Thrombosis Outcomes Study), anti-inflammatory therapy targeting the interleukin-1β pathway significantly reduced macrovascular events.40 Although not previously investigated, the finding from the CANTOS trial may suggest that interventions aimed at controlling low-grade inflammation may be valuable in curbing the rates of microvascular and macrovascular dysfunction. This may include relatively less expensive interventions like control of chronic infections and infestations. Another clinical significance of our results is the potential of using CRP as a marker for assessing or monitoring the risk of developing vascular dysfunction in individuals with and without diabetes. This could aid early detection of vascular dysfunction, improve quality of life and reduce healthcare-related costs through preventive interventions.
Limitations
Our study has some limitations. First, we used hs-CRP as the marker of inflammation. hs-CRP does not necessarily reflect the degree of inflammation of individual tissues. Additionally, CRP is only one marker of inflammation and other markers such as fibrinogen and D-dimer were not measured in our study. Second, we did not have data on existing acute infections or infestations; however, we excluded participants with CRP concentrations >10 mg/L in a sensitivity analysis. Third, the sample size among individuals with diabetes was relatively small, limiting the power to detect associations between CRP and vascular complications in diabetes. Fourth, the microvascular dysfunction of retinopathy and neuropathy, as well as the macrovascular dysfunction of coronary artery disease and cerebrovascular disease, were not assessed. Fifth, we did not perform the 2-hour postload plasma glucose in identifying participants with diabetes; this could have underdiagnosed participants with diabetes.24 Finally, the duration of diabetes was not included as a covariate in the multivariable analysis for participants with diabetes because many study participants did not provide this information.
Conclusion
In our study, higher CRP concentration was associated with higher odds of diabetes. Higher CRP concentration was also associated with higher odds of PAD and nephropathy in non-diabetes and higher odds of PAD in diabetes. Our data from this first study on CRP in relation to diabetes and microvascular and macrovascular dysfunction in a sub-Saharan African population suggest that CRP could be explored as a potential marker to identify sub-Saharan African patients with diabetes at risk for macrovascular complications. Additionally, our finding of an association with higher odds of diabetes and both microvascular and macrovascular dysfunction in non-diabetes individuals implies that CRP could play a role in the assessment of risk before diabetes diagnosis, for example, in pre-diabetes individuals. In both individuals with and without diabetes, interventions aimed at controlling inflammation may help reduce the risk of diabetes and microvascular and macrovascular dysfunction.
Acknowledgments
We are grateful to the research assistants, interviewers, and other staff of the five research locations who have taken part in gathering the data, and most of all the Ghanaian volunteers participating in the RODAM (Research on Obesity and Diabetes Among African Migrants) study. We gratefully acknowledge the advisory board members for their valuable support in shaping the RODAM study methods, Jan van Straalen from the Amsterdam University Medical Centre with standardization of the laboratory procedures, and the Amsterdam University Medical Centre Biobank for their support in biobank management and high-quality storage of collected samples.
References
Footnotes
Contributors All authors have contributed substantially to this article and approved the submission. CFH-B, AHMZ, B-JB, AGBA and COA conceived the idea. CFH-B and COA were responsible for data acquisition and statistical analysis. CFH-B, AHMZ, B-JB, AGBA, KACM, KK-G, MBS, JS, ID, EJAJB, FM, and COA were responsible for data analysis/interpretation. Each author contributed important intellectual content during article drafting or revision and accepts accountability for the overall work by ensuring that questions pertaining to the accuracy or integrity of any portion of the work are appropriately investigated and resolved. CFH-B takes responsibility for the fact that this study has been reported honestly, accurately and transparently, that no important aspects of the study have been omitted, and that any discrepancies from the study as planned have been explained.
Funding This work was supported by the European Commission under the Framework Programme (grant number: 278901). KACM is supported by the Intramural Research Program of the National Institutes of Health at the Center for Research on Genomics and Global Health (CRGGH). The CRGGH is supported by the National Human Genome Research Institute, the National Institute of Diabetes and Digestive and Kidney Diseases, the Center for Information Technology, and the Office of the Director at the National Institutes of Health (1ZIAHG200362).
Competing interests None declared.
Patient consent for publication Not required.
Ethics approval Ethical approval of the study protocols was requested at all sites from the respective ethics committees in Ghana (School of Medical Sciences/Komfo Anokye Teaching Hospital Committee on Human Research, Publication and Ethical Review Board), the Netherlands (Institutional Review Board of the AMC, University of Amsterdam (reference number: W12_062 # 12.17.0086)), Germany (Ethics Committee of Charite-Universitätsmedizin Berlin), and the UK (London School of Hygiene and Tropical Medicine Research Ethics Committee) before data collection began in each country.
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement Data are available upon reasonable request. The data sets created and/or analyzed during the current study are available from the corresponding author on reasonable request.