key: cord-0972288-j93fc0ap authors: Asgarpur, Golschan; Treskatsch, Sascha; Angermair, Stefan; Danassis, Michaela; Nothnagel, Anna Maria; Toepper, Christoph; Trauzeddel, Ralf Felix; Nordine, Michael; Heeschen, Julia; Al-Chehadeh, Alaa; Landmesser, Ulf; Sander, Leif Erik; Kurth, Florian; Berger, Christian title: Echocardiographic Evaluation of Right Ventricular (RV) Performance over Time in COVID-19-Associated ARDS—A Prospective Observational Study date: 2021-05-01 journal: J Clin Med DOI: 10.3390/jcm10091944 sha: 150caa6c423b30cc50ab32bbc9d8cc05476a6edf doc_id: 972288 cord_uid: j93fc0ap (1) Background: To evaluate time-dependent right ventricular (RV) performance in patients with COVID-19-associated acute respiratory distress syndrome (ARDS) undergoing intensive care (ICU) treatment. (2) Methods: This prospective observational study included 21 ICU patients with COVID-19-associated ARDS in a university hospital in 2020 (first wave). Patients were evaluated by transthoracic echocardiography at an early (EE) and late (LE) stage of disease. Echocardiographic parameters describing RV size and function as well as RV size in correlation to P(a)O(2)/F(i)O(2) ratio were assessed in survivors and nonsurvivors. (3) Results: Echocardiographic RV parameters were within normal range and not significantly different between EE and LE. Comparing survivors and nonsurvivors revealed no differences in RV performance at EE. Linear regression analysis did not show a correlation between RV size and P(a)O(2)/F(i)O(2) ratio over all measurements. Analysing EE and LE separately showed a significant increase in RV size correlated to a lower P(a)O(2)/F(i)O(2) ratio at a later stage of COVID-19 ARDS. (4) Conclusion: The present study reveals neither a severe RV dilatation nor an impairment of systolic RV function during the initial course of COVID-19-associated ARDS. A trend towards an increase in RV size in correlation with ARDS severity in the second week after ICU admission was observed. The 2019 Coronavirus disease caused by the novel SARS Coronavirus 2 (SARS-CoV-2) is responsible for a global pandemic, infecting millions of people worldwide. In certain cases, a fulminant SARS-CoV-2 infection led to severe acute respiratory distress syndrome (ARDS), as well affecting other organ systems, such as the nervous and the cardiovascular system [1, 2] . The incidence of hospitalized patients with COVID-19 developing ARDS is approximately 33% [3] and is associated with poor clinical outcome and a high mortality rate, reaching up to 45% according to recent research [4] . 2 The pathophysiology of non-COVID-associated ARDS has been described as an increased permeability to liquid, protein and cellular compounds across the lung endothelium, leading to interstitial edema, which further translocate across the alveolar barrier into the alveolar space, impairing oxygenation [5] . The resulting hypoxemia caused by a ventilation-to-perfusion mismatch as well as right-to-left intrapulmonary shunting provokes compensatory hypoxic pulmonary vasoconstriction (HPV) with consecutive pulmonary hypertension and increased afterload, leading to right ventricular (RV) impairment or failure [6, 7] . Non-COVID ARDSassociated RV failure is associated with a mortality rate of up to 60-70% [2, 6, 8] . Virus-associated RV impairment has been documented during the outbreak of H1N1 in 2009 [9] . However, in contrast to the pathophysiological course of non-COVID ARDS, early experiences with COVID-19-associated ARDS raised the suspicion for a loss of pulmonary vascular tone. Therapeutic strategies aiming to dilate alveolar pulmonary vessels in the area of alveoli still participating in gas exchange in order to improve gas exchange, with for example, inhaled nitric oxide, were found to be not as clinically effective as deemed [10] . Thus, we hypothesized that hypoxemia in COVID-19-associated ARDS may present as a different hemodynamic phenotype without increases in afterload, thus impairing RV performance. To investigate this hypothesis, we evaluated time-dependent RV performance in COVID-19-associated ARDS intensive care (ICU) patients using transthoracic echocardiography (TTE) as a noninvasive diagnostic and readily available bedside tool [7, 8] . For evaluating the impact of COVID-19 ARDS on RV performance, we conducted a monocenter prospective observational study at the department of anesthesiology and intensive care medicine, Campus Benjamin Franklin, Charité-Universitätsmedizin Berlin. This study was constituted as a substudy of the PA-COVID-19 trial [11] , approved by Charité's Ethics committee (EA2/066/20). Informed consent was obtained from all participants or their authorized representative. The trial was registered with the following number DRKS00021688 on 13 May 2020 (WHO International Clinical Trials Registry Platform). Adult patients with positive SARS-CoV-2 polymerase chain reaction (PCR) testing admitted to our intensive care unit (ICU) requiring treatment due to COVID-19 ARDS in accordance with the Berlin Definition [4] were included in this study during the first wave in Germany in 2020. Patients with a palliative approach to therapy were excluded. All participants were scheduled for an early TTE evaluation (EE) of RV performance within the first week of ICU treatment. For evaluating a possible time-dependent impact on RV performance, an additional "late" transthoracic echocardiography (LE) was conducted in the second week after ICU admission. All TTE examinations were performed by echocardiography-trained ICU practitioners according to national standards [12] and actual guidelines [13] . For evaluation of RV performance, the 2-dimensional apical-4-chamber (A4C) or subcostal 4-chamber view (SC4C) were obtained using a VIVID S60 ultrasound system (GE Healthcare, Chicago, IL, USA). The following parameters were analyzed using EchoPac (GE Healthcare, Chicago, USA) according to recent guidelines: RV medial diameter (RVMD), RV end diastolic area index (RVEDAi), RV/LV medial diameter ratio (RLDR), RV/LV area ratio (RLAR), tricuspid annular plain systolic excursion (TAPSE) and RV fractional area change (RFAC) [13] . Demographic, morphometric, laboratory, respiratory/ventilatory and hemodynamic data were obtained from two patient data management systems (COPRA System GmbH, Sasbachwalden, Germany and SAP AG, Walldorf, Germany) at two time points: EE and LE. All data are available on demand. Descriptive analyses and statistical testing were performed using IBM SPSS Statistics (version 25; IBM, Armonk, NY, USA) with a p-value below 0.05 regarded as significant. Unless otherwise stated, all data are presented as median and interquartile range (IQR). Statistical significance among groups was analyzed by the exact nonparametric Mann-Whitney U test or Wilcoxon single rank test. Exact chi-square tests were used for qualitative data. Linear regression (r 2 ) was calculated to detect correlations between parameters describing RV performance and oxygenation index as surrogate for ARDS severity. All tests should be understood as constituting explorative analysis, and no adjustment for multiple testing was performed. Between March and May 2020, 28 patients were assessed for eligibility. Twenty-one patients were included in this study, while six did not meet inclusion criteria (no ARDS according to Berlin definition), and one patient was transferred to a different ICU after initial evaluation (Figure 1 ). Statistical significance among groups was analyzed by the exact nonparametric Mann-Whitney U test or Wilcoxon single rank test. Exact chi-square tests were used for qualitative data. Linear regression (r 2 ) was calculated to detect correlations between parameters describing RV performance and oxygenation index as surrogate for ARDS severity. All tests should be understood as constituting explorative analysis, and no adjustment for multiple testing was performed. Between March and May 2020, 28 patients were assessed for eligibility. Twenty-one patients were included in this study, while six did not meet inclusion criteria (no ARDS according to Berlin definition), and one patient was transferred to a different ICU after initial evaluation (Figure 1 ). Table 1 . Except for hyperlipoproteinemia (HLP), which found was more amongst nonsurvivors (survivors 7%, nonsurvivors 43%; p = 0.007), both groups exhibited comparable age, gender, BMI and comorbidities (for detailed and patient individual information for comorbidities, see Supplementary Table S1 ). Laboratory parameters at ICU admission showed no differences between survivors and nonsurvivors (Supplemental Table S2 ). Median (IQR) age of all included patients was 68 (59/76) years with a body mass index (BMI) of 27.8 (24.0/31.8) kg/m 2 and a gender distribution of 9/12 (female/male). Seven out of 21 patients died, accounting for a mortality rate of 33%. Median time from ICU admission to death was 20 (12/30) days. Detailed morphometric and demographic data of survivors and nonsurvivors at ICU admission are presented in Table 1 . Except for hyperlipoproteinemia (HLP), which found was more amongst nonsurvivors (survivors 7%, nonsurvivors 43%; p = 0.007), both groups exhibited comparable age, gender, BMI and comorbidities (for detailed and patient individual information for comorbidities, see Supplementary Table S1 ). Laboratory parameters at ICU admission showed no differences between survivors and nonsurvivors (Supplemental Table S2 ). BMI: body mass index; CAD: coronary artery disease; AHT: arterial hypertension; HF: heart failure; COPD: chronic obstructive pulmonary disease; IDDM: insulin-dependent diabetes mellitus; NIDDM: non-insulin-dependent diabetes mellitus; CKD: chronic kidney disease; CLD: chronic liver disease; HLP: hyperlipoproteinemia; PVD: peripheral vascular disease. Increased APACHE II and SAPS II scores at ICU admission depict an increased disease severity in nonsurvivors despite comparable ratios of arterial partial pressure of oxygen to fraction of inspired oxygen (P a O 2 /F i O 2 ) ( Table 2 ). Seventeen out of 21 patients (81%) required invasive mechanical ventilation at ICU admission (survivors: n = 11; nonsurvivors n = 6) ( Table 2) , and nearly all patients had bacterial superinfections during ICU treatment (Supplementry Table S2 ). Time from ICU admission to EE was 1.7 (0.4/3.7) days and 11.0 (7.9/12.8) days to LE. No patient died before the LE examination. Clinical conditions were comparable between EE and LE (Table 3) as well as between survivors and nonsurvivors at EE (Table 4 ). Computer tomography (CT) confirmed segmental pulmonary embolism (PE) in 5 of 21 (24%) patients during the observation period without significant differences between survivors (n = 3; 21%) and nonsurvivors (n = 2; 29%) (p = 0.717). Echocardiographic parameters evaluating RV size (RVMD, RVEDAi, RLDR, RLAR) and systolic function (TAPSE, RFAC) were not significantly different between both time points (Table 5 ). Comparing survivors with nonsurvivors also revealed no differences in RV performance parameters at EE (Table 6 ). Paradoxical movement of the interventricular septum (IVS) was observed in only one patient at EE. Linear regression analysis including all echocardiographic measurements did not detect a correlation between RV size and P a O 2 /F i O 2 ratio (RLDR: r 2 = 0.000, p = 0.907; RLAR: r 2 = 0.02, p = 0.466) (Figure 2a) . Separate evaluation of time-dependent RV performance also found no correlation between RV size and P a O 2 /F i O 2 ratio at EE (RLDR: r 2 = 0.149, p = 0.084; RLAR: r 2 = 0.002, p = 0.835) (Figure 2b ). Nevertheless, a trend towards a slight, though significant increase in RV size associated with a decreased P a O 2 /F i O 2 ratio became noticeable in the second week of ICU treatment at LE (RLDR: r 2 = 0.27, p = 0.047; RLAR: r 2 = 0.168, p = 0.129) (Figure 2c ). Linear regression analysis including all echocardiographic measurem tect a correlation between RV size and P a O 2 /F i O 2 ratio (RLDR: r 2 = 0.000, p r 2 = 0.02, p = 0.466) (Figure 2a) . Separate evaluation of time-dependent R also found no correlation between RV size and P a O 2 /F i O 2 ratio at EE (RLD = 0.084; RLAR: r 2 = 0.002, p = 0.835) (Figure 2b ). Nevertheless, a trend to though significant increase in RV size associated with a decreased P a O 2 /F i O noticeable in the second week of ICU treatment at LE (RLDR: r 2 = 0.27, p = = 0.168, p = 0.129) (Figure 2c ). Linear regression analysis including all echocardiographic measurements did not detect a correlation between RV size and PaO2/FiO2 ratio (RLDR: r 2 = 0.000, p = 0.907; RLAR: r 2 = 0.02, p = 0.466) (Figure 2a) . Separate evaluation of time-dependent RV performance also found no correlation between RV size and PaO2/FiO2 ratio at EE (RLDR: r 2 = 0.149, p = 0.084; RLAR: r 2 = 0.002, p = 0.835) (Figure 2b ). Nevertheless, a trend towards a slight, though significant increase in RV size associated with a decreased PaO2/FiO2 ratio became noticeable in the second week of ICU treatment at LE (RLDR: r 2 = 0.27, p = 0.047; RLAR: r 2 = 0.168, p = 0.129) (Figure 2c) . This observational study is the first to evaluate the time-dependent impact of COVID-19-associated ARDS on RV performance and demonstrated a virtually unaffected RV with absence of severe RV dilatation within the first two weeks after ICU admission using two-dimensional echocardiography. We were also unable to detect any discernible RV differences between survivors and nonsurvivors during the observation period via This observational study is the first to evaluate the time-dependent impact of COVID-19-associated ARDS on RV performance and demonstrated a virtually unaffected RV with absence of severe RV dilatation within the first two weeks after ICU admission using two-dimensional echocardiography. We were also unable to detect any discernible RV differences between survivors and nonsurvivors during the observation period via echocardiographic bedside examinations. Although RV values remained within normal or acceptable ranges, according to proposed values in COVID-19 ARDS, non-COVID ARDS and/or ASE/ESC guidelines [13] [14] [15] , a trend towards RV enlargement correlating with a decrease in P a O 2 /F i O 2 ratio-as surrogate marker for ARDS severity-was observed during the second week of evaluation. These results may support previous assumptions that COVID-19 ARDS may present as an altered hemodynamic phenotype without severe impairment of RV performance. RV dilatation or dysfunction in non-COVID ARDS is understood as a result of physiological changes in the pulmonary circulation [16] . The mechanism of acute cor pulmonale in non-COVID ARDS has been established as refractory pulmonary edema due to endothelial cell swelling and hypoxemia, leading to HPV in precapillary arterioles and consecutively increasing RV afterload. Additional microvascular thrombosis deriving from endothelial cell activation can further promote RV deterioration, and later remodeling can lead to persistent pulmonary hypertension [17, 18] . According to this pathophysiological pathway, we assumed that RV function during COVID-19-associated ARDS would follow suit, however, we were not able to detect any meaningful RV impairment in COVID-19 ARDS patients in our results. For higher accuracy, we determined RV size not only in terms of absolute diameter, but also in correlation to LV size using two different methods: end-diastolic RV/LV diameter (RLDR) and enddiastolic RV/LV area (RLAR). With a calculated RLDR and RLAR between 0.68 and 0.72 in conjunction with a preserved systolic RV function, our COVID-19 ARDS cohort revealed, at most, only a slight increase in RV size. Whether or not this increase implies a clinically relevant RV impairment remains controversial, as former studies regarding disease impact on RV set a much higher limit for a relevant RV/LV ratio increase [19] [20] [21] . Analyzing timeddependent RV performance in our ICU cohort, an absence of severe RV impairment was also observed at a later stage of disease. Further, to investigate if a more severe COVID-19 ARDS may correlate with a pronounced RV impairment, we analyzed linear regression of P a O 2 /F i O 2 ratios with RV/LV ratios. Based on experience from non-COVID ARDS pathophysiology, a more severe ARDS should lead to a higher RV/LV ratio. Again, this assumption could not be confirmed in our COVID-19 ARDS cohort as we observed only a slight increase in RV in relation to LV size correlating with a more impaired oxygenation at a later stage of disease. Interestingly, critically ill patients suffering from H1N1 infection frequently exhibit RV dilatation and failure [9] , whereas a recent study showed a reduced incidence of RV failure in ventilated COVID-19 patients [22] . These findings may reflect the varying clinical manifestation due to different pathophysiological changes in non-COVID and COVID-19 ARDS. The precise pathological mechanism in COVID-19 ARDS is still unknown, but comparisons to other viral infections support the hypothesis of a different clinical phenotype with altered pulmonary vascular reactivity [22, 23] promoting increased vascular permeability [24] finally leading to vasoplegia [10] . Impaired endothelial function without HPV may thus provide an explanation for the sustained absence of severe RV dysfunction among the COVID-19 ARDS population over time, as demonstrated by our findings. In this context, the course of COVID-19 as well as non-COVID ARDS may be assumed as being possibly bimodal in nature, or as some authors have suggested, with different phenotypes [10, 25] . This theory may be supported by our observed shift to a slightly higher correlation of increased RV/LV ratio with a decreased P a O 2 /F i O 2 ratio in the second week after ICU admission. Furthermore, it should be taken into consideration that during the course of virally induced ARDS, bacterial superinfection is a common complication [22, 26] . This was also present in our study population, where 80.9% of patients had a proven bacterial superinfection during their ICU treatment (Supplementary Table S3 ). Hence, it cannot be ruled out that the ARDS phenotype may change over time. It is thus reasonable that the later course of COVID-19 ARDS may present as a mixture of viral and bacterial ARDS, which may be reflected in the observed change in linear regression analysis at LE. The destruction of pulmonary tissue, which we frequently observed in our COVID-19 ARDS (Supplemental Figure S1 ), may provide another explanation for our findings, especially in the late course. Such destruction may lead to a rarefication of the pulmonary vascular bed resulting in an increased RV afterload. Another possible explanation for a change in RV performance over time may be pulmonary embolism with consecutive RV impairment. Segmental PE was detected in five (21%) patients during the observation period, which is within the reported incidence of PE in COVID-19 [27] . No central PE occurred among our population and comparison of RV parameters in PE and non-PE patients revealed no significant differences in RV performance (Supplementary Table S4 ). Therefore, additional PE seems to be unlikely to have a relevant impact on the findings in this study. Only limited data for RV performance in COVID-19 ARDS exist to date and the studies vary in their finding. Other than the above-mentioned pathophysiological theories, some methodical explanations may account for these differences. Recently presented results for RV function in COVID-19 were obtained from single, inconsistent time points among COVID-19 infected patients but without differentiation between presence or absence of ARDS or sepsis [28] [29] [30] [31] . Other working groups have presented data with partly mechanically ventilated COVID-19 patients (30%) [29] or exclusion of ventilated COVID-19 patients [30] , while echocardiographic examinations were performed. Another work from Bagate et al. reported at a single time point-with a focus on filling pressures-a high incidence of cor pulmonale, but RV functional parameters like TAPSE were found to be comparable to our findings [28] . These different approaches may lead to a heterogeneity of studied populations and may impact RV evaluation findings in COVID-19 patients. Therefore, we assessed RV performance in every patient, irrespective of ventilation status at an early state after onset of COVID-19-associated ARDS followed by a second evaluation within the second week of treatment. Nevertheless, we did not find any differences in parameters describing sole RV performance as well as in correlation with ARDS severity between survivors and nonsurvivors. The observed slight but not significant decrease of TAPSE in nonsurvivors is somewhat in contrast to a recent study from D'Alto and colleagues. They investigated TAPSE in relation to echocardiographic estimations of systolic pulmonary arterial pressure (sPAP) among a comparable cohort of mostly mechanically ventilated COVID-19 ARDS patients [32] . A significant reduction in TAPSE and TAPSE/sPAP ratio between survivors and nonsurvivors was observed, but despite statistical significance, TAPSE in nonsurvivors was also still within normal ranges confirming our results. In addition, a significant increase in C-reactive protein and procalcitonin were observed only among nonsurvivors in their study, which could be interpreted as a higher incidence of bacterial superinfections at the time of RV assessment. This finding cannot be confirmed by our results as most patients suffered bacterial superinfections irrespective of survival. This study has some limitations. Due to the course of ARDS, with the need of invasive ventilation, including high PEEP and periodic change to prone position, obtaining adequate TTE views at exactly the same time point was not possible. Further, TTE evaluation in such patients is challenging, therefore, our results are somewhat limited due to the absence of sPAP evaluation. The assessment of RV function via TTE amongst ARDS patients has been observed to be inferior to transesophageal echocardiography (TEE), as TTE is prone to interference due to anatomical barriers such as the chest wall [33] . Additionally, Evrard et al. found that cardiac assessment via TTE was suboptimal compared with TEE assessment amongst COVID-19 patients, and the primary reason for this is the ability to obtain short axis cardiac images with TEE [22] . Because of the increased complexity to obtain proper Doppler-derived sPAP in a positive pressure ventilated ARDS population, we decided to utilize a simple and reproducible protocol for TTE evaluation over time. Furthermore, all participants required (non)invasive ventilatory support at the time of RV evaluation. It may be reasonable to suggest that this intervention improved P a O 2 /F i O 2 ratios, thereby preserving RV performance. In this context a baseline TTE evaluation before initiation of ARDS therapy would have been beneficial, but ethically not feasible. Because of the nature of this study, a previous sample size calculation was not possible. Additionally, due to the small sample size in this study, drawing broader conclusions concerning RV function in the majority of COVID-19 patients may be limited. A higher number of cases might have led to more statistically significant differences between groups, as well as increasing the general power of the study. Our work reveals neither severe RV dilatation nor impairment of systolic RV function during the initial course of COVID-19-associated ARDS. In contrast to the initial evaluation, a trend towards an increase in RV size in correlation with ARDS severity in the second week after ICU admission was detectable. These findings support the assumptions of impaired regulation of pulmonary vascular and/or endothelial dysfunction due to SARS-CoV-2 and hint towards a change of the COVID-19-associated ARDS phenotype over time. However, final evidence for these assumptions is still missing, and further research into the hemodynamic changes during the course of COVID-19 ARDS needs to be performed. The following are available online at https://www.mdpi.com/article/10 .3390/jcm10091944/s1, Table S1 : Known medical history, Table S2 : Laboratory data on ICU admission from survivors vs. nonsurvivors, Table S3 : Presumed superinfections and proven superinfections, pulmonary and extrapulmonary during ICU stay, Table S4 : TTE evaluated parameters for right ventricular function at EE in patients without PE compared to patients with PE, Figure S1 : Pat. No 4, CT Scan of the lung at LE, after first week of ICU stay due to COVID-19 ARDS (mechanically ventilated). Endothelial dysfunction in COVID-19: A position paper of the ESC Working Group for Atherosclerosis and Vascular Biology, and the ESC Council of Basic Cardiovascular Science The Right Ventricle in ARDS Incidence of ARDS and outcomes in hospitalized patients with COVID-19: A global literature survey Acute respiratory distress syndrome: The Berlin Definition Acute respiratory distress syndrome Acute cor pulmonale in ARDS: Rationale for protecting the right ventricle Echocardiographic assessment of the right ventricle and its correlation with patient outcome in acute respiratory distress syndrome COVID-19, and echocardiography Right and left heart failure in severe H1N1 influenza A infection COVID-19 pneumonia: ARDS or not? Studying the pathophysiology of coronavirus disease 2019: A protocol for the Berlin prospective COVID-19 patient cohort (Pa-COVID-19) Perioperative fokussierte Echokardiographie: Modul 1-5. A&I 2017 Guidelines for the echocardiographic assessment of the right heart in adults: A report from the American Society of Echocardiography endorsed by the European Association of Echocardiography, a registered branch of the European Society of Cardiology, and the Canadian Society of Echocardiography Echocardiographic assessment of the right ventricle in COVID-related acute respiratory syndrome Prevalence and prognosis of cor pulmonale during protective ventilation for acute respiratory distress syndrome Right ventricular performance in adult respiratory distress syndrome Incidence of thrombotic complications in critically ill ICU patients with COVID-19 Pathophysiology of pulmonary hypertension in acute lung injury Right-to-left ventricular end diastolic diameter ratio in severe sepsis and septic shock Prognostic value of echocardiographic right/left ventricular end-diastolic diameter ratio in patients with acute pulmonary embolism: Results from a monocenter registry of 1416 patients Right-toleft ventricular end-diastolic diameter ratio and prediction of right ventricular failure with continuous-flow left ventricular assist devices Cardiovascular phenotypes in ventilated patients with COVID-19 acute respiratory distress syndrome Acute respiratory failure in COVID-19: Is it "typical" ARDS? The Enigma of Endothelium in COVID-19 Compliance Phenotypes in Early Acute Respiratory Distress Syndrome before the COVID-19 Pandemic Clinical features of patients infected with 2019 novel coronavirus in Acute Pulmonary Embolism and COVID-19 Advanced echocardiographic phenotyping of critically ill patients with coronavirus-19 sepsis: A prospective cohort study Right Ventricular Dilation in Hospitalized Patients with COVID-19 Infection Echocardiographic features of patients with COVID-19 infection: A cross-sectional study Spectrum of Cardiac Manifestations in COVID-19: A Systematic Echocardiographic Study Right ventricular-arterial uncoupling independently predicts survival in COVID-19 ARDS Prevalence and prognostic value of acute cor pulmonale and patent foramen ovale in ventilated patients with early acute respiratory distress syndrome: A multicenter study Informed consent was obtained from all subjects involved in the study or their authorized representatives. Written informed consent has been obtained from the patients to publish this paper. The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ethical restrictions.