- 1Department of Neurosurgery, RWTH Aachen University, Aachen, Germany
- 2Institute for Neurophysiology, University of Cologne, Cologne, Germany
- 3Department of Medical Engineering and Technomathematics, FH Aachen University of Applied Sciences, Aachen, Germany
- 4Department of Ophthalmology, RWTH Aachen University, Aachen, Germany
- 5Department of Intensive Care and Intermediate Care, RWTH Aachen University, Aachen, Germany
- 6Department of Exercise and Health Sciences, University of Basel, Basel, Switzerland
- 7Department of Diagnostic and Interventional Neuroradiology, RWTH Aachen University, Aachen, Germany
Objective: Delayed cerebral ischemia (DCI) is a common complication after aneurysmal subarachnoid hemorrhage (aSAH) and can lead to infarction and poor clinical outcome. The underlying mechanisms are still incompletely understood, but animal models indicate that vasoactive metabolites and inflammatory cytokines produced within the subarachnoid space may progressively impair and partially invert neurovascular coupling (NVC) in the brain. Because cerebral and retinal microvasculature are governed by comparable regulatory mechanisms and may be connected by perivascular pathways, retinal vascular changes are increasingly recognized as a potential surrogate for altered NVC in the brain. Here, we used non-invasive retinal vessel analysis (RVA) to assess microvascular function in aSAH patients at different times after the ictus.
Methods: Static and dynamic RVA were performed using a Retinal Vessel Analyzer (IMEDOS Systems GmbH, Jena) in 70 aSAH patients during the early (d0−4), critical (d5−15), late (d16−23) phase, and at follow-up (f/u > 6 weeks) after the ictus. For comparison, an age-matched cohort of 42 healthy subjects was also included in the study. Vessel diameters were quantified in terms of the central retinal arterial and venous equivalent (CRAE, CRVE) and the retinal arterio-venous-ratio (AVR). Vessel responses to flicker light excitation (FLE) were quantified by recording the maximum arterial and venous dilation (MAD, MVD), the time to 30% and 100% of maximum dilation (tMAD30, tMVD30; tMAD, tMVD, resp.), and the arterial and venous area under the curve (AUCart, AUCven) during the FLE. For subgroup analyses, patients were stratified according to the development of DCI and clinical outcomes after 12 months.
Results: Vessel diameter (CRAE, CRVE) was significantly smaller in aSAH patients and showed little change throughout the whole observation period (p < 0.0001 vs. control for all time periods examined). In addition, aSAH patients exhibited impaired arterial but not venous responses to FLE, as reflected in a significantly lower MAD [2.2 (1.0–3.2)% vs. 3.6 (2.6–5.6)% in control subjects, p = 0.0016] and AUCart [21.5 (9.4–35.8)%*s vs. 51.4 (32.5–69.7)%*s in control subjects, p = 0.0001] on d0−4. However, gradual recovery was observed during the first 3 weeks, with close to normal levels at follow-up, when MAD and AUCart amounted to 3.0 [2.0–5.0]% (p = 0.141 vs. control, p = 0.0321 vs. d5−15) and 44.5 [23.2–61.1]%*s (p = 0.138 vs. control, p < 0.01 vs. d0−4 & d5−15). Finally, patients with clinical deterioration (DCI) showed opposite changes in the kinetics of arterial responses during early and late phase, as reflected in a significantly lower tMAD30 on d0−4 [4.0 (3.0–6.8) s vs. 7.0 (5.0–8.0) s in patients without DCI, p = 0.022) and a significantly higher tMAD on d16−23 (24.0 (21.0–29.3) s vs. 18.0 (14.0–21.0) s in patients without DCI, p = 0.017].
Conclusion: Our findings confirm and extend previous observations that aSAH results in sustained impairments of NVC in the retina. DCI may be associated with characteristic changes in the kinetics of retinal arterial responses. However, further studies will be required to determine their clinical implications and to assess if they can be used to identify patients at risk of developing DCI.
Trial Registration: ClinicalTrials.gov Identifier: NCT04094155.
Introduction
Aneurysmal subarachnoid hemorrhage (aSAH) due to sudden rupture of a cerebral aneurysm is responsible for only 5–7% of all stroke events (1) but is associated with a fatality rate of almost 50%, accounting for up to 20% of all cerebrovascular-related deaths (2). Owing to a complex pathophysiology that involves early and delayed forms of brain injury, a significant number of patients surviving aSAH remain disabled with or without compromise detected on detailed neuropsychological testing. Only 10% of patients recover completely. The most common complication responsible for long-term impairments is delayed cerebral ischemia (DCI), which is typically observed during the first two or at most 3 weeks after aSAH onset and is associated with the appearance of a new focal neurologic deficit or a persistent decline in the patient's Glasgow Coma Scale score (3). Although traditionally attributed to vasoconstriction of large caliber proximal vessels (i.e., angiographic vasospasm) secondary to degradation of subarachnoid blood and intrathecal immune activation, more recent findings also point to a critical role of microvascular dysfunction (4–7). In particular, several studies performed in animal models have demonstrated a progressive impairment and partial inversion of cerebral neurovascular coupling (NVC) (8–10), the process responsible for adjusting local cerebral blood flow (CBF) to regionally heterogeneous changes in metabolic demand due to neuronal activation. While there is evidence for similar impairment in human aSAH patients (11–13), clinical assessment of NVC remains a diagnostic challenge that is usually achieved through invasive CBF measurements or highly specialized and expensive functional imaging techniques. An alternative approach to evaluate NVC is to measure microvascular responses in the retina, an embryological part of the central nervous system that can be assessed non-invasively by retinal vessel analysis (RVA). Because cerebral and retinal microvasculature bear a close anatomical correlation and are governed by common regulatory mechanisms, retinal vascular changes are increasingly recognized as potential screening option for altered NVC in the brain (14, 15). For example, a potential role of RVA for predicting cerebrovascular changes has previously been shown in elderly adults (16, 17), after coarctation repair (18) and in patients with diabetes (19) or Alzheimer's disease (20). During aSAH, the sudden rise of intracranial pressure often forces subarachnoid blood into the pre-retinal space (21), which may expose retinal micro-vessels to the same blood degradation products and inflammatory cytokines thought to be responsible for microvascular dysfunction in the brain. In addition, recent discovery of a cerebral and retinal glymphatic system (22) argues for the existence of perivascular pathways through which vasoactive agents formed in the subarachnoid space may also directly reach the retinal microvasculature. We have previously demonstrated the feasibility of RVA in aSAH patients (23) and provided first evidence for retinal vasoconstriction and signs of impaired NVC in these patients compared to healthy control subjects (11). In this prospective study on a larger cohort of aSAH patients, we sought to validate our previous findings and evaluate the relationship between potential changes in NVC as detected by RVA and the occurrence of DCI and long-term clinical outcome.
Materials and Methods
Patients with aSAH admitted to the Department of Neurosurgery at the RWTH Aachen University Hospital between July 2015 and December 2019 were considered for enrollment into the present study. Patients were selected based on the following inclusion criteria: (1) age ≥ 18 years, (2) aneurysmal origin of SAH as confirmed by computed tomography (CT) angiography and/or digital subtraction angiography (DSA), (3) absence of intraocular comorbidities or contraindications that would have precluded application of mydriatic agents, and (4) absence of Terson syndrome unless it did not interfere with the analysis of selected retinal vessels. Before enrolment, written informed consent according to the Declaration of Helsinki was obtained from all patients or their legal representatives and the study protocol was approved by the local Research Ethics Committee of the RWTH Aachen University (EK 069/15).
A total of 70 patients with aSAH were prospectively recruited, in whom measurements were performed 0–4 days (early phase: d0−4), 5–15 days (critical phase: d5−15), 16–23 days (late phase: d16−23) and more than 6 weeks (follow-up: f/u) after the ictus. If more than one measurement was available for a given patient in a given phase, all measurements in this phase were averaged to calculate a mean valued. For comparison, an age-matched cohort of 42 healthy subjects who were measured by RVA at the Department of Exercise and Health Sciences at the University of Basel was also included in the study.
The clinical severity of aSAH on admission was assessed according to the Hunt and Hess (HH) grading scale and the degree of aneurysmal bleeding was scored according to the modified Fischer (mFS) scale based on CT scans. Patients were treated according to current multidisciplinary consensus guidelines (3, 24, 25). The aneurysm was secured by neurosurgical clipping or endovascular coiling based on an interdisciplinary consensus. Patients received prophylactic nimodipine therapy (180–360 mg per day) unless precluded by hemodynamic instability. All patients were continuously monitored for complications such as delayed cerebral ischemia (DCI), which was considered as a neurological worsening according to the definition by Vergouwen et al. (3) or territorial or watershed hypoperfusion on perfusion CT scans. In analgosedated patients, deterioration in cerebral microdialysis (71 high cut-off catheter, μdialysis, Stockholm, Sweden) or cerebral tissue oxygen saturation (Neurovent PTO, Raumedic AG, Helmbrechts, Germany) were used as surrogates to trigger CT perfusion. First-tier treatment in case of DCI consisted of induced hypertension (>180 mmHg systolic blood pressure) as per institutional protocol. In refractory cases without clinical or neuromonitoring improvement but confirmed misery perfusion on CTA or DSA, second-tier rescue therapy by balloon dilatation and/or continuous intra-arterial infusion of nimodipine was considered based on angiographic findings (26, 27). Clinical outcome after 12 months according to the extended Glasgow Outcome Scale (GOS-E) was assessed by an independent investigator based on a telephone interview, clinical investigation in the outpatient clinic or clinical status compiled from the medical reports. A summary of the demographic and clinical data obtained on admission or during in-hospital treatment is provided in Table 1.
Retinal Vessel Analysis
Static and dynamic retinal vessel analyses were performed using a maneuverable Retinal Vessel Analyzer (IMEDOS Systems GmbH, Jena, Germany) as described previously (11, 23). Briefly, after unilateral application of a mydriatic agent (Tropicamide, Mydriaticum Stulln UD, Pharma Stulln GmbH, Stulln, Germany), conscious patients were examined in a sitting position while analgo-sedated patients were examined in supine position with the head and torso partially turned. For assessment of NVC, retinal veins and arteries were focused and the vessel diameters recorded continuously in operator-selected regions of interest located equidistant to the papilla. Measurements were performed with a standard 350 s dynamic vessel analysis protocol comprised of 50 s of baseline recording and three cycles of flicker-light excitation (FLE). For extended analysis the recorded vessel diameters were expressed as percentage of the baseline diameter and plotted as a function of time as described in detail elsewhere (20, 28). To improve the signal-to-noise ratio, the individual response curves obtained for each subject during the three successive cycles of FLE were averaged. In healthy subjects, the typical arterial response curves thus obtained feature primary vasodilation upon initiation of FLE, which reaches a maximum after a characteristic latency and is followed, after termination of the stimulus, by a reflectory vasoconstriction (29). The general shape and amplitude of the response curve provide an index for the integrity of NVC, and were quantified by calculating maximum arterial and venous dilatation (MAD, MVD), the arterial and venous area-under-the-curve (AUCart, AUCven) and the time to reach 30% (tMAD30, tMVD30) and 100% (tMAD, tMVD) of MAD and MVD after flicker initiation (Figure 1A). Following the dynamic analysis, monochromatic fundus images for static retinal vessel analysis (Figure 1B) were obtained at a 30° camera angle with the optic nerve centered and processed using a dedicated software tool (VisualIS, IMEDOS Systems GmbH., Jena, Germany). Central retinal arterial and venous equivalents (CRAE and CRVE, respectively), and arterio-venous ratio (AVR) were calculated as described previously (30, 31). All dynamic and static examinations were routinely paralleled by multimodal monitoring of mean arterial blood pressure (MAP), heart rate, blood gases, oxygenation and intracranial pressure (ICP), ensuring physiological conditions at the time of image acquisition. Measurement was aborted in case of hemodynamic or ventilatory instability.
Figure 1. Static and dynamic retinal vessel analysis (RVA). (A) Average retinal vessel response to flicker light excitation (FLE) in healthy subjects with the different parameters determined for dynamic RVA. The stimulation period is indicated above and by light gray shading. Parameters used for quantification of the response comprised maximum arterial or venous dilation (MAD/MVD), time to maximum arterial or venous dilation (tMAD/tMVD) relative to flicker initiation, time to 30% of the maximum arterial or venous dilation (tMAD30/tMVD30) relative to flicker initiation, and arterial or venous area under the curve (AUCart/AUCven) during the FLE. (B) Example of a monochromatic fundus image used for RVA. Arterial (A) and venous (V) segments are indicated in red and blue, respectively.
Statistical Analysis
Normal distribution was tested using the Kolmogorow-Smirnow normality test. Quantitative parameter values are presented as mean ± standard deviation (SD) or as median [1 quartile−3 quartile], and as percentage in case of nominal data. Differences between two groups were analyzed using the two-sided student t-test, Mann–Whitney-test or Chi-Square test, respectively. For subgroup analyses with regard to DCI and outcome, all patients were stratified into groups according to the development of DCI (DCI vs. no DCI) and clinical outcome after 12 months (unfavorable: GOS-E1−4 vs. favorable: GOS-E5−8), respectively (see Supplementary Tables 1, 2). For analysis of the influence of nimodipine, all measurements (see Supplementary Table 3) or all patients (see Supplementary Table 4) were stratified into groups according to oral nimodipine treatment (measurements in nimodipine-treated patients vs. measurements in patients with no nimodipine treatment for at least 24 h). Statistical significance was set at p < 0.05; statistical results with p < 0.1 were accepted as a trend. All analyses were performed using Excel (Microsoft Office Excel 2010, Californian, USA), SPSS v. 21 (IBM Chicago, Illinois, USA), Numbers®, Apple Inc., Cupertino, USA, GraphPad Prism® and GraphPad Software, Inc., La Jolla, USA.
Results
The demographic data of patients included in this study are summarized in Table 1 In total, the aSAH cohort comprised 70 patients with a median age of 51 [48–59] years, of whom 51 (73%) were female and 19 (27%) were male. Twenty-three (32%) of the patients developed delayed cerebral ischemia (DCI), which progressed to infarction in one patient. The clinical outcome after 12 months was favorable (GOS-E5−8) in 64 (92%) and unfavorable (GOS-E1−4) in 6 (8%) of the patients. Retinal vessel analysis (RVA) could be performed in 32, 53, 23, and 40 patients during the early (d0−4), critical (d5−15), late (d16−23) phase, and follow-up, respectively. For comparison, a cohort of 42 healthy subjects with a median age of 50 [43–59] years, of whom 20 (48%) were female and 22 (52%) were male, was also included in the study (Table 1).
Vessel Diameters
When compared to control subjects, the arterial diameter, expressed as central retinal arterial equivalent (CRAE), was significantly reduced immediately after aSAH and remained so throughout the whole observation period (p < 0.0001 for all time periods examined), with only a minor, non-significant increase in vessel diameters observed in measurements performed at follow-up (f/u) more than 6 weeks after the ictus (Figure 2A, Table 2). Likewise, the central retinal venous equivalent (CRVE) was significantly reduced in aSAH patients compared to control subjects (p < 0.0001 for all time periods examined) and showed little change during the observation period (Figure 2B, Table 3). For relative assessment of vessel dimensions and to account for systematic variations in measurements, the retinal arterio-venous-ratio (AVR) was also calculated. As illustrated in Figure 2C, the AVR showed very similar changes, although the difference to the control subjects only reached statistical significance during the early phase (d0−4, p = 0.0396) and after more than 6 weeks (f/u, p = 0.0050).
Figure 2. Results from static RVA in aSAH patients and control subjects. Comparison of (A) central retinal arterial equivalent (CRAE) as a measure of arterial diameter, (B) central retinal venous equivalent (CRVE) as a measure of venous diameter and (C) retinal arterio-venous-ratio (AVR) in control subjects and aSAH patients during the early (aSAHd0−4), critical (aSAHd5−15) and late (aSAHd16−23) phase and at follow-up more than 6 weeks after the ictus (aSAHf/u). ***p < 0.001 vs. control; **p < 0.01 vs. control; *p < 0.05 vs. control.
Neurovascular Coupling
Dynamic retinal vessel analysis revealed dramatic and sustained impairments in the arterial response of aSAH patients to flicker light stimulation when compared to the control subjects (Figure 3A). In particular, aSAH patients exhibited an acute reduction in arterial vasodilation during the flicker, as reflected in a significantly lower MAD [2.2 (1.0–3.2)% vs. 3.8 (2.6–5.6)% in control subjects, p = 0.0016] and AUCart [21.5 (9.4–35.8)%*s vs. 51.4 (32.5–69.7)%*s in control subjects, p = 0.0001] during the early phase (d0−4, Figures 3B,C), while the timing of the response (i.e., tMAD and tMAD30) appeared to be largely unaffected (Table 2). Unlike the changes in arterial diameter and AVR, both MAD and AUCart showed a gradual recovery toward healthy controls, with a significant reduction only during the first 3 weeks after the ictus (Figures 3B,C). For example, AUCart increased to 22.4 [12.6–45.0]%*s (p < 0.0001 vs. Control) during the critical (d5−15) and 30.5 [3.2–53.9]%*s (p = 0.0033 vs. Control) during the late (d16−23) phase (Figure 3C). In measurements performed at follow-up more than 6 weeks after the ictus, AUCart already amounted to 44.5 [23.2–61.1]%*s, which is close to the value of 51.4 [32.5–69.7]%*s (p = 0.1380) observed in the control subjects and significantly higher than the values observed during the early [21.5 (9.4–35.8)%*s, p = 0.0066] and critical [22.4 (12.6–45.0)%*s, p = 0.0047] phase (Figure 3C).
Figure 3. Results from dynamic RVA in aSAH patients and control subjects. (A) Average relative retinal arterial responses to flicker light stimulation measured in control subjects and aSAH patients during the early (aSAHd0−4), critical (aSAHd5−15), and late (aSAHd16−23) phase and at follow-up more than 6 weeks after the ictus (aSAHf/u). Stimulation period is indicated by light gray shading. (B) Maximum arterial dilation (MAD) quantified from data shown in A. (C) Arterial area under the curve (AUCart ) quantified from the data shown in A. (D) Average retinal venous responses to flicker light stimulation measured in control subjects and aSAH patients at the same times as in A. (E) Maximum venous dilation (MVD) quantified from the data shown in D. (F) Venous area under the curve (AUCven) quantified from the data shown in A. ***p < 0.001 vs. control; **p < 0.01 vs. control; *p = 0.05 vs. control; ##p < 0.01 vs. aSAHf/u; #p < 0.05 vs. aSAHf/u.
Analysis of venous responses to flicker light revealed much less marked differences (Figure 3D, Table 3). In particular, neither MVD (Figure 3E) nor the timing of venous responses (Table 3) showed significant differences compared to control subjects. Moreover, there was only a tendency for AUCven during the early phase to be reduced when compared to control subjects [41.3 (16.4–46.8)%*s vs. 45.7 (31.2–62.2)%*s in control subjects, p = 0.0610], although a significant difference was detected when compared to the corresponding value at follow-up [41.3 (16.4–46.8)%*s vs. 52.4 (34.2–61.7)%*s at follow-up, p = 0.0205; Figure 3F].
Subgroup Analyses With Regard to Cardiovascular Risk Factors
Cardiovascular risk factors like diabetes mellitus have been shown to alter retinal vessel diameters and/or responses to stimulation (32) and could thus act as important confounding factors. There was only one patient with diabetes in our study, excluding further analysis of its potential as a confounding factor. However, to address the role of other cardiovascular risk factors like arterial hypertension, obesity and smoking, we examined their potential association with RVA parameters determined during the early phase. There was no difference in any of the parameters between patients stratified according to smoking status, but a tendency for AUCven to be lower in patients in hypertension [36.1 (12.8–43.2)%*s vs. 44.4 (27.6–56.5)%*s in normotensive patients, p = 0.0818], which may have contributed to the difference in AUCven between aSAH and control patients during the early phase described above. In addition, while AUCart tended to be lower in patients with BMI ≥ 30 [8.1 (−10.1–12.5)%*s vs. 26.9 (11.6–36.6)%*s in patients with BMI <30, p = 0.062], it actually tended to be higher in patients with arterial hypertension [27.8 (14.9–44.8)%*s vs. 11.9 (4.1–32.0)%*s in normotensive patients, p = 0.0502], although none of these differences reached statistical significance (for details see Supplementary Table 1).
Subgroup Analyses With Regard to DCI and Long-Term Clinical Outcome
To examine if differences in retinal vessel responses could facilitate early detection of patients at risk of complications and delayed brain injury, we also analyzed data stratified according to the occurrence of DCI and long-term clinical outcome. While most parameters were independent of the subgroup (Supplementary Tables 2, 3), there were some subtle but significant kinetic differences between arterial responses in patients with and without DCI. In particular, patients that developed DCI showed faster vasodilation on flicker initiation during the early phase (Figure 4A), which was most evident when the average responses in both groups were normalized by their maximum amplitude (inset and arrow in Figure 4A). Quantitatively, this difference was reflected in a significantly lower tMAD30 (Figure 4B), which amounted to 4.0 [3.0–6.8] s in patients with DCI (n = 10) vs. 7.0 [5.0–8.0] s (p = 0.022) in patients without DCI (n = 13). There was also a tendency for a lower tMAD in patients with DCI, but the difference was much less pronounced and did not reach statistical significance (p = 0.254, Figure 4C). Interestingly, these differences appeared to be restricted to the early phase and their direction actually reversed during the late phase (Figure 4D), where patients with DCI exhibited a higher tMAD30 (Figure 4E) and a significantly higher tMAD [24.0 (21.0–29.3) s, n = 6 vs. 18.0 (14.0–21.0) s in patients without DCI, n = 17, p = 0.017, Figure 4F]. In contrast, no significant differences were observed in venous retinal vessel parameters between patients with and without DCI (Supplementary Table 2). Comparison of results with regard to the clinical outcome after 12 months was hampered by the fact that RVA could only be performed in a total of six patients with a poor clinical outcome. However, patients with an unfavorable outcome were characterized by a tendency for lower arterial and venous diameters throughout the whole observation period (Supplementary Table 3), and the CRVE in these patients during the critical phase [d5−15: 171 (144–187) MU, n = 4] was significantly (p = 0.045) lower than the corresponding value in patients with a favorable outcome [199 (170–220) MU, n = 43]. Arterial and venous parameters from dynamic retinal vessel analysis on the other hand either showed no significant differences or the number of patients with a poor outcome was simply too small for statistical testing (Supplementary Table 3).
Figure 4. Comparison of retinal arterial responses in aSAH patients with and without DCI. (A) Average relative retinal arterial responses to flicker light stimulation measured during the early phase (days 0–4) in aSAH patients with or without delayed cerebral ischemia (DCI). Stimulation period is indicated by light gray shading. Inset shows the rising phase on an expanded time-scale and after scaling the responses to the same amplitude to highlight kinetic changes. (B) Time to 30% of the maximum arterial dilation (tMAD30) and (C) time to maximum arterial dilation (tMAD) determined during the early phase in patients with and without DCI. (D) Time-course of changes in median [q1–q3] tMAD30 (bottom) and tMAD (top) observed in patients with or without DCI. (E) tMAD30 and (F) tMAD determined during the late phase (days 16–23) in patients with and without DCI. For a summary of all parameters and p-values see Supplementary Table 2.
Effects of Nimodipine Treatment
Finally, to examine if retinal vessel diameters or their response to flicker stimulation were affected by treatment with the L-type Ca2+ channel antagonist nimodipine, we compared data obtained during nimodipine treatment and data obtained while no nimodipine was administered (Supplementary Table 4). However, neither static nor dynamic retinal vessel parameters appeared to be affected by nimodipine. For example, in 51 static measurements performed during nimodipine treatment, arterial and venous diameters (CRAE & CRVE) were 158 [136–174] MU and 199 [171–214] MU, while the corresponding values in 30 measurements performed while no nimodipine was administered were 162 [133–178] MU (p = 0.4868) and 200 [174–209] MU (p = 0.9692). Likewise, maximum arterial and venous dilation upon stimulation (MAD and MVD) amounted to 2.2 [1.4–3.6]% and 3.6 [2.6–5.0]% during nimodipine treatment as compared to 2.9 [1.2–4.6]% (p = 0.4343) and 3.7 [2.4–4.7]% (p = 0.8484) in measurements performed while no nimodipine was administered (Figures 5A,B, for remaining parameters see Supplementary Table 4). Similar results were obtained when the same analysis was performed after stratification of the patients according to the occurrence of DCI (see Supplementary Table 4). Finally, to ensure that this was not due to differences in the time-frame of measurements performed with and without nimodipine treatment, we also compared data obtained during early, critical and late phase respectively (Figures 5C–H, Supplementary Table 6). However, neither retinal vessel diameters (Figures 5C,D), nor the maximum amplitude (Figures 5E,F), AUC (Supplementary Table 6), or time-course (Figures 5G,H) of vasodilation in response to flicker light stimulation during any of the phases appeared to be affected by nimodipine treatment.
Figure 5. Dependence of retinal vessel properties in aSAH patients on nimodipine treatment. (A) Average relative arterial responses to flicker light stimulation measured with or without nimodipine (nimo) treatment. Inset shows the same data but scaled to the same amplitude for comparison of response kinetics. Stimulation period is indicated by light gray shading. (B) Average relative venous responses to flicker light stimulation measured with or without nimodipine (nimo) treatment. Inset shows the same data but scaled to the same amplitude for comparison of response kinetics. Stimulation period is indicated by light gray shading. (C–H) Comparison of (C) central retinal arterial equivalent (CRAE), (D) central retinal venous equivalent (CRVE), (E) maximum arterial dilation (MAD), (F) maximum venous dilation (MVD), (G) time to 30 and 100% of maximum arterial dilation (tMAD30 & tMAD), and (H) time to 30 and 100% of maximum venous dilation (tMVD30 & tMVD) in aSAH patients with or without nimodipine treatment during the early (aSAHd0−4), critical (aSAHd5−15), and late (aSAHd16−23) phase.
Discussion
The exact pathophysiological cascades underlying delayed cerebral ischemia (DCI) and other complications in aSAH patients surviving the initial bleed are still poorly understood. Evidence from animal models indicates that a mismatch between metabolic demand and local CBF due to impaired NVC could contribute to DCI, but the clinical relevance of these findings has remained elusive as there is still a lack of methods for non-invasive assessment of NVC in human subjects. A promising approach to evaluate NVC in aSAH patients is to measure microvascular responses in the retina, an embryological part of the central nervous system and possibly exposed to the same vasoactive hemoglobin metabolites and inflammatory cytokines as small cerebral vessels. Thus, recent discovery of a cerebral and retinal glymphatic system (22) indicates that there are perivascular pathways through which vasoactive agents formed in the subarachnoid space could directly reach the retinal microvasculature. Here, we used static and dynamic retinal vessel analysis (RVA) in aSAH patients to detect changes in basal diameter and responsiveness of retinal vessels to flicker light stimulation, both of which have been proposed as potential candidates for early detection of cerebrovascular dysfunction (14, 15).
Static RVA and Retinal Vessel Diameters After aSAH
Vasoconstriction of large caliber proximal vessels (i.e., angiographic vasospasm) is a common complication after aSAH and was long regarded as the main cause of DCI and poor clinical outcomes, but more recent studies have questioned this concept and proposed a more important role of microvascular changes (4, 5). Our present findings point to a sustained reduction in the diameter of retinal microvessels after aSAH compared to control subjects, which appeared to be more pronounced in the arterial compartment, as reflected in a parallel and significant decrease of the arterio-venous-ratio (AVR) during the early phase and at follow-up. The latter would be consistent with the general view that veins have only limited capability for active vasoconstriction and our previous proposal that the reduced venule diameter in the retina after aSAH is at least partly mediated by constriction of upstream arterioles and a resulting reduction in blood flow (11). In any case, the observed microvascular changes resemble aSAH-induced changes in large caliber vessels, namely early cerebral artery vasoconstriction and an associated reduction in cerebral perfusion which may contribute to DCI in some patients (33–35). Considering that increased intracranial pressure during the acute phase of aSAH forces subarachnoid blood into the preretinal space, a significant reduction in retinal vessel diameters observed during the early phase after aSAH is in line with previous in vitro and in vivo animal studies showing that blood-filled cerebrospinal fluid enhances constriction of both, the retinal (36) and cerebral (37) microvasculature. Interestingly and in contrast to angiographic vasospasm, which usually dissolves by week 2 after the ictus (34), however, there was almost no return of arterial and venous diameters toward the control level during the first 3 weeks after aSAH and only a partial, non-significant recovery at follow-up more than 6 weeks later. A possible explanation for these observations is that part of the reduced basal vessel diameter reflects pathophysiological mechanisms that are not readily reversible and outlast the impairments in vessel responsiveness to stimulation. This would be in line with the results from animal studies indicating that, at least in large caliber vessels, the initial arterial narrowing by vasoactive substances after aSAH leads to tissue damage and structural changes with an associated long-term alteration of the arterial tone (38, 39). In this context, it is also interesting to note that perivascular enlargement in the brain, a putative imaging biomarker for microvascular brain damage (40), has recently been demonstrated to be associated with narrower arterial calibers in the retina (41).
Dynamic RVA and Neurovascular Coupling After aSAH
Adequate NVC is critical for continuous adjustment of local CBF to regionally heterogeneous changes in metabolic demand but may be compromised by the pathophysiological cascades initiated during bleeding into the subarachnoid space. Our present findings demonstrate a sustained decrease of retinal arterial vasodilation with partial recovery in response to flicker light stimulation in aSAH patients, as reflected in a significantly reduced maximum amplitude and area under the response curve compared to control subjects. Unlike the changes in vessel diameters, arterial responsiveness in aSAH patients showed a gradual recovery during the first 3 weeks and reached almost normal levels at follow-up more than 6 weeks after the ictus. Given that complications like DCI and DCI-related infarctions, which have been proposed to be the result of microvascular dysfunction in the brain, are typically restricted to the first two or at most 3 weeks after the ictus as well, these findings support the assumption that retinal vessels could serve as a surrogate marker for changes in cerebral small vessels after aSAH. The observed time-course also matches the pathophysiological cascades traditionally thought to be responsible for vascular dysfunction after aSAH, namely hemolysis of subarachnoid blood and intrathecal immune activation during the early phase, secondary immune infiltration and formation of various vasoactive hemoglobin metabolites during the critical phase and their gradual clearance during the late phase (42–44). Further support for the idea comes from recent discovery of cerebral and retinal glymphatic systems (22), which could provide a perivascular pathway through which vasoactive agents formed during the critical phase reach the retinal microvasculature. If this turns out to be the case, dynamic RVA could not only provide a unique tool to monitor their impact on microvascular function in a relatively simple manner but also a read-out of disease- or treatment-related changes in cerebral small vessels that are otherwise difficult or impossible to assess directly and non-invasively. For example, our present findings suggest that nimodipine, the only FDA-approved pharmacological treatment for delayed cerebral vasospasm after aSAH (45), neither affects the basal diameter nor the responsiveness of retinal microvessels to stimulation, although this remains to be confirmed it future studies.
Likewise, further study will be required to evaluate the sensitivity of RVA with regard to the early detection of DCI.
Limitations and Confounding Factors
There are several confounding factors that have to be taken into account when interpreting our results. First and foremost, although adaptation of our current imaging technique allowed us to also recruit some more severely affected, analgo-sedated patients, cases with mild aSAH, good clinical grade (HH1−3), and/or good clinical outcome were still clearly overrepresented in our cohort, which may have biased the observed changes, likely toward milder alterations. To overcome this limitation, which is mainly brought about by the limited maneuverability of the required funduscope in supine patients, we are currently working on a miniaturized, wearable device for RVA that is comparable to commercially available “smartglasses” and should facilitate bedside testing of immobilized and sedated patients (46). That said, given the complexity of blood flow regulation and NVC, it is likely that numerous factors, some of which may still be unknown, can influence retinal vessel diameters and their responsiveness to stimulation. In order to reduce potential confounding factors and ensure standardized, physiological conditions at the time of image acquisition, our measurements were routinely paralleled by multimodal monitoring of mean arterial blood pressure (MAP), heart rate, blood gases, oxygenation, and intracranial pressure (ICP). In addition, our analysis indicates that there was no effect of nimodipine treatment on any of the parameters evaluated. With regard to static RVA, it is also important to keep in mind that direct comparison of CRAE and CRVE values between control and aSAH patients may have been influenced by the fact that measurements were performed using devices at two different facilities.
It is also worth noting that cardiovascular risk factors like diabetes or arterial hypertension could act as confounding factors and complicate the use of RVA in the context of aSAH. For example, in the present study, there was a tendency for AUCart to be higher in patients with arterial hypertension, which could partly mask the decrease of this parameter after aSAH. Consistent with previous studies (32), these differences were relatively subtle and they did not reach statistical significance, but further studies on the role of these and other cardiovascular risk factors will clearly be required before RVA can be reliably used for diagnostic applications in the context of aSAH.
Finally, it should be noted that in contrast to our previous study, we observed no significant effects of aSAH on the latency of arterial or venous responses, which could reflect the fact that these parameters are more sensitive to inherent technical limitations such as the finite signal-to-noise ratio. As such, our present findings that patients with DCI exhibit a reduced arterial latency (tMAD30) during the early but increased arterial latency (tMAD) during the late phase should be interpreted with care as well-until they can be confirmed in future studies.
In particular, given the subtle nature of these differences, it seems unlikely that they would allow for reliable prediction and/or confirmation of DCI in the clinical setting. However, based on the small number of patients we evaluated, future studies with larger patient populations may be able to identify more robust correlations between RVA parameters and DCI. Considering that DCI remains one of the main preventable causes of poor clinical outcomes after aSAH, such studies will almost certainly be performed in the near future and their results should help to conclusively rate the clinical value of RVA in this setting.
Conclusions
Our findings confirm and extend the results from previous studies that aSAH is associated with sustained vasoconstriction and impairments of NVC in retinal vessels. Retinal vessel responses may differ between patients with and without DCI. Although their clinical relevance remains to be firmly established, these findings suggest that RVA could be a unique tool for monitoring microvascular function in aSAH patients in a simple and non-invasive manner.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
All procedures involving human participants were performed in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Reference Number EK 069/15. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
WA, CC, MWe, KS, DK, and GS: conceived, designed, and performed the experiments. WA: first drafting of the manuscript. WA and FN: illustrations. WA, CC, MWe, KS, KK, and GS: data acquisition. WA, CC, KK, and GS: analysis and interpretation of data. WA, CC, TS, JH, AB, AS-T, FN, MWi, HC, and GS: critical review of the manuscript. All authors were critically revised and approved by the final manuscript.
Funding
WA was supported by funding from the German Society of Neurointensive and Emergency Medicine (Junior research award 2017). CC was supported by the START-Program of the Faculty of Medicine (RWTH Aachen University, Grant No. 691545) and by the Foundation of Neurosurgical Research (German Society of Neurosurgery, 2016).
Conflict of Interest
IMEDOS Systems GmbH provided the Retinal Vessel Analyzer for research purposes only. The company did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2021.690183/full#supplementary-material
References
1. Kaptain GJ, Lanzino G, Kassell NF. Subarachnoid haemorrhage: epidemiology, risk factors, and treatment options. Drugs Aging. (2000) 17:183–99. doi: 10.2165/00002512-200017030-00003
2. Weaver JP, Fisher M. Subarachnoid hemorrhage: an update of pathogenesis, diagnosis and management. J Neurol Sci. (1994) 125:119–31. doi: 10.1016/0022-510X(94)90024-8
3. Vergouwen MDI, Vermeulen M, van Gijn J, Rinkel GJE, Wijdicks EF, Muizelaar JP, et al. Definition of delayed cerebral ischemia after aneurysmal subarachnoid hemorrhage as an outcome event in clinical trials and observational studies: proposal of a multidisciplinary research group. Stroke. (2010) 41:2391–5. doi: 10.1161/STROKEAHA.110.589275
4. Meyers PM, Connolly ES. Disappointing results for clazosentan in CONSCIOUS-2. Nat Rev Neurol. (2011) 7:660–1. doi: 10.1038/nrneurol.2011.168
5. Naraoka M, Matsuda N, Shimamura N, Asano K, Ohkuma H. The role of arterioles and the microcirculation in the development of vasospasm after aneurysmal SAH. Biomed Res Int. (2014) 2014:253746. doi: 10.1155/2014/253746
6. Geraghty JR, Davis JL, Testai FD. Neuroinflammation and microvascular dysfunction after experimental subarachnoid hemorrhage: emerging components of early brain injury related to outcome. Neurocrit Care. (2019) 31:373–89. doi: 10.1007/s12028-019-00710-x
7. Suzuki H, Kanamaru H, Kawakita F, Asada R, Fujimoto M, Shiba M. Cerebrovascular pathophysiology of delayed cerebral ischemia after aneurysmal subarachnoid hemorrhage. Histol Histopathol. (2021) 36:143–58. doi: 10.14670/HH-18-253
8. Balbi M, Koide M, Wellman GC, Plesnila N. Inversion of neurovascular coupling after subarachnoid hemorrhage in vivo. J Cereb Blood Flow Metab. (2017) 37:3625–34. doi: 10.1177/0271678X16686595
9. Balbi M, Vega MJ, Lourbopoulos A, Terpolilli NA, Plesnila N. Long-term impairment of neurovascular coupling following experimental subarachnoid hemorrhage. J Cereb Blood Flow Metab. (2020) 40:1193–202. doi: 10.1177/0271678X19863021
10. Koide M, Bonev AD, Nelson MT, Wellman GC. Inversion of neurovascular coupling by subarachnoid blood depends on large-conductance Ca2+-activated K+ (BK) channels. Proc Natl Acad Sic USA. (2012) 109:E1387–95. doi: 10.1073/pnas.1121359109
11. Conzen C, Albanna W, Weiss M, Kürten D, Vilser W, Kotliar K, et al. Vasoconstriction and impairment of neurovascular coupling after subarachnoid hemorrhage: a descriptive analysis of retinal changes. Transl Stroke Res. (2018) 9:284–93. doi: 10.1007/s12975-017-0585-8
12. Winkler MKL, Chassidim Y, Lublinsky S, Revankar GS, Major S, Kang E-J, et al. Impaired neurovascular coupling to ictal epileptic activity and spreading depolarization in a patient with subarachnoid hemorrhage: possible link to blood-brain barrier dysfunction. Epilepsia. (2012) 53:22–30. doi: 10.1111/j.1528-1167.2012.03699.x
13. da Costa L, Fierstra J, Fisher JA, Mikulis DJ, Han JS, Tymianski M. BOLD MRI and early impairment of cerebrovascular reserve after aneurysmal subarachnoid hemorrhage. J Magn Reson Imaging. (2014) 40:972–9. doi: 10.1002/jmri.24474
14. Cabrera DeBuc D, Somfai GM, Koller A. Retinal microvascular network alterations: Potential biomarkers of cerebrovascular and neural diseases. Am J Physiol Circ Physiol. (2017) 312:H201–12. doi: 10.1152/ajpheart.00201.2016
15. Patton N, Aslam T, MacGillivray T, Pattie A, Deary IJ, Dhillon B. Retinal vascular image analysis as a potential screening tool for cerebrovascular disease: a rationale based on homology between cerebral and retinal microvasculatures. J Anat. (2005) 206:319–48. doi: 10.1111/j.1469-7580.2005.00395.x
16. Lipecz A, Csipo T, Tarantini S, Hand RA, Ngo BTN, Conley S, et al. Age-related impairment of neurovascular coupling responses: a dynamic vessel analysis (DVA)-based approach to measure decreased flicker light stimulus-induced retinal arteriolar dilation in healthy older adults. GeroScience. (2019) 41:341–9. doi: 10.1007/s11357-019-00078-y
17. Golzan SM, Goozee K, Georgevsky D, Avolio A, Chatterjee P, Shen K, et al. Retinal vascular and structural changes are associated with amyloid burden in the elderly: ophthalmic biomarkers of preclinical Alzheimer's disease. Alzheimers Res Ther. (2017) 9:13. doi: 10.1186/s13195-017-0239-9
18. Pressler A, Esefeld K, Scherr J, Ali M, Hanssen H, Kotliar K, et al. Structural alterations of retinal arterioles in adults late after repair of aortic isthmic coarctation. Am J Cardiol. (2010) 105:740–4. doi: 10.1016/j.amjcard.2009.10.070
19. Bettermann K, Slocomb J, Shivkumar V, Quillen D, Gardner TW, Lott ME. Impaired retinal vasoreactivity: an early marker of stroke risk in diabetes. J Neuroimaging. (2017) 27:78–84. doi: 10.1111/jon.12412
20. Kotliar K, Hauser C, Ortner M, Muggenthaler C, Diehl-Schmid J, Angermann S, et al. Altered neurovascular coupling as measured by optical imaging: a biomarker for Alzheimer's disease. Sci Rep. (2017) 7:12906. doi: 10.1038/s41598-017-13349-5
21. Fountas KN, Kapsalaki EZ, Lee GP, Machinis TG, Grigorian AA, Robinson JS, et al. Terson hemorrhage in patients suffering aneurysmal subarachnoid hemorrhage: predisposing factors and prognostic significance. J Neurosurg. (2008) 109:439–44. doi: 10.3171/JNS/2008/109/9/0439
22. Wostyn P, De Groot V, Van Dam D, Audenaert K, De Deyn PP, Killer HE. The glymphatic system: a new player in ocular diseases? Investig Opthalmology Vis Sci. (2016) 57:5426. doi: 10.1167/iovs.16-20262
23. Albanna W, Conzen C, Weiss M, Clusmann H, Fuest M, Mueller M, et al. Retinal vessel analysis (RVA) in the context of subarachnoid hemorrhage–a proof of concept study. PLoS ONE. (2016) 11:e0158781. doi: 10.1371/journal.pone.0158781
24. Le Roux P, Menon DK, Citerio G, Vespa P, Bader MK, Brophy GM, et al. Consensus summary statement of the international multidisciplinary consensus conference on multimodality monitoring in neurocritical care. Intensive Care Med. (2014) 40:1189–1209. doi: 10.1007/s00134-014-3369-6
25. Hutchinson P, O'Phelan K. International multidisciplinary consensus conference on multimodality monitoring: cerebral metabolism. Neurocrit Care. (2014) 21:148–58. doi: 10.1007/s12028-014-0035-3
26. Albanna W, Weiss M, Müller M, Brockmann MA, Rieg A, Conzen C, et al. Endovascular rescue therapies for refractory vasospasm after subarachnoid hemorrhage: a prospective evaluation study using multimodal, continuous event neuromonitoring. Neurosurgery. (2017) 80:942–9. doi: 10.1093/neuros/nyw132
27. Weiss M, Conzen C, Mueller M, Wiesmann M, Clusmann H, Albanna W, et al. Endovascular rescue treatment for delayed cerebral ischemia after subarachnoid hemorrhage is safe and effective. Front Neurol. (2019) 10:136. doi: 10.3389/fneur.2019.00136
28. Kotliar KE, Lanzl IM, Schmidt-Trucksäss A, Sitnikova D, Ali M, Blume K, et al. Dynamic retinal vessel response to flicker in obesity: a methodological approach. Microvasc Res. (2011) 81:123–28. doi: 10.1016/j.mvr.2010.11.007
29. Polak K. Evaluation of the zeiss retinal vessel analyser. Br J Ophthalmol. (2000) 84:1285–90. doi: 10.1136/bjo.84.11.1285
30. Nagel E, Vilser W, Fink A, Riemer T. Statische Gefäßanalyse in nonmydriatischen und mydriatischen Bildern. Klin Monbl Augenheilkd. (2007) 224:411–6. doi: 10.1055/s-2007-963093
31. Hubbard LD, Brothers RJ, King WN, Clegg LX, Klein R, Cooper LS, et al. Methods for evaluation of retinal microvascular abnormalities associated with hypertension/sclerosis in the atherosclerosis risk in communities study. Ophthalmology. (1999) 106:2269–80. doi: 10.1016/S0161-6420(99)90525-0
32. Heitmar R, Lip GYH, Ryder RE, Blann AD. Retinal vessel diameters and reactivity in diabetes mellitus and/or cardiovascular disease. Cardiovasc Diabetol. (2017) 16:56. doi: 10.1186/s12933-017-0534-6
33. Bederson JB, Levy AL, Hong Ding W, Kahn R, DiPerna CA, Jenkins AL, et al. Acute vasoconstriction after subarachnoid hemorrhage. Neurosurgery. (1998) 42:352–62. doi: 10.1097/00006123-199802000-00091
34. Weir B, Grace M, Hansen J, Rothberg C. Time course of vasospasm in man. J Neurosurg. (1978) 48:173–8. doi: 10.3171/jns.1978.48.2.0173
35. Ferguson S, Macdonald RL. Predictors of cerebral infarction in patients with aneurysmal subarachnoid hemorrhage. Neurosurgery. (2007) 60:658–67. doi: 10.1227/01.NEU.0000255396.23280.31
36. Liu Z, Li Q, Cui G, Zhu G, Tang W, Zhao H, et al. Blood-filled cerebrospinal fluid-enhanced pericyte microvasculature contraction in rat retina: a novel in vitro study of subarachnoid hemorrhage. Exp Ther Med. (2016) 12:2411–6. doi: 10.3892/etm.2016.3644
37. Wang KC, Tang SC, Lee JE, Tsai JC, Lai DM, Lin WC, et al. Impaired microcirculation after subarachnoid hemorrhage in an in vivo animal model. Sci Rep. (2018) 8:13315. doi: 10.1038/s41598-018-31709-7
38. Vorkapic P, Bevan JA, Bevan RD. Longitudinal in vivo and in vitro time-course study of chronic cerebrovasospasm in the rabbit basilar artery. Neurosurg Rev. (1991) 14:215–9. doi: 10.1007/BF00310660
39. Peeyush Kumar T, McBride DW, Dash PK, Matsumura K, Rubi A, Blackburn SL. Endothelial cell dysfunction and injury in subarachnoid hemorrhage. Mol Neurobiol. (2019) 56:1992–2006. doi: 10.1007/s12035-018-1213-7
40. Groeschel S, Chong WK, Surtees R, Hanefeld F. Virchow-Robin spaces on magnetic resonance images: normative data, their dilatation, and a review of the literature. Neuroradiology. (2006) 48:745–754. doi: 10.1007/s00234-006-0112-1
41. Mutlu U, Adams HHH, Hofman A, van der Lugt A, Klaver CCW, Vernooij MW, et al. Retinal microvascular calibers are associated with enlarged perivascular spaces in the brain. Stroke. (2016) 47:1374–6. doi: 10.1161/STROKEAHA.115.012438
42. Clark JF, Sharp FR. Bilirubin oxidation products (BOXes) and their role in cerebral vasospasm after subarachnoid hemorrhage. J Cereb Blood Flow Metab. (2006) 26:1223–33. doi: 10.1038/sj.jcbfm.9600280
43. Kikuchi T, Okuda Y, Kaito N, Abe T. Cytokine production in cerebrospinal fluid after subarachnoid haemorrhage. Neurol Res. (1995) 17:106–8. doi: 10.1080/01616412.1995.11740296
44. Roa JA, Sarkar D, Zanaty M, Ishii D, Lu Y, Karandikar NJ, et al. Preliminary results in the analysis of the immune response after aneurysmal subarachnoid hemorrhage. Sci Rep. (2020) 10:11809. doi: 10.1038/s41598-020-68861-y
45. Maruhashi T, Higashi Y. An overview of pharmacotherapy for cerebral vasospasm and delayed cerebral ischemia after subarachnoid hemorrhage. Expert Opin Pharmacother. (2021) 1–14. doi: 10.1080/14656566.2021.1912013
Keywords: aneurysmal subarachnoid hemorrhage, cerebral infarction, delayed cerebral ischemia, microvascular function, neurovascular coupling, retinal vessel analysis
Citation: Albanna W, Conzen C, Weiss M, Seyfried K, Kotliar K, Schmidt TP, Kuerten D, Hescheler J, Bruecken A, Schmidt-Trucksäss A, Neumaier F, Wiesmann M, Clusmann H and Schubert GA (2021) Non-invasive Assessment of Neurovascular Coupling After Aneurysmal Subarachnoid Hemorrhage: A Prospective Observational Trial Using Retinal Vessel Analysis. Front. Neurol. 12:690183. doi: 10.3389/fneur.2021.690183
Received: 02 April 2021; Accepted: 18 May 2021;
Published: 14 June 2021.
Edited by:
Chiara Robba, University of Genoa, ItalyReviewed by:
Masayo Koide, University of Vermont, United StatesAlessandra Splendiani, University of L'Aquila, Italy
Copyright © 2021 Albanna, Conzen, Weiss, Seyfried, Kotliar, Schmidt, Kuerten, Hescheler, Bruecken, Schmidt-Trucksäss, Neumaier, Wiesmann, Clusmann and Schubert. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Walid Albanna, WalidAlbanna@yahoo.de
†These authors have contributed equally to this work