Impact Factor 2.889 | CiteScore 3.1

More on impact ›

Methods ARTICLE

Front. Neurol., 08 September 2020 | https://doi.org/10.3389/fneur.2020.00971

Eye Movement Evaluation in Multiple Sclerosis and Parkinson's Disease Using a Standardized Oculomotor and Neuro-Ophthalmic Disorder Assessment (SONDA)

Alessandro Grillini1*, Remco J. Renken2, Anne C. L. Vrijling3, Joost Heutink3,4 and Frans W. Cornelissen1
  • 1Laboratory for Experimental Ophthalmology, University Medical Center Groningen, University of Groningen, Groningen, Netherlands
  • 2Department of Ophthalmology, University Medical Center Groningen, University of Groningen, Groningen, Netherlands
  • 3Royal Dutch Visio, Center of Expertise for Blind and Partially Sighted People, Huizen, Netherlands
  • 4Department of Clinical and Developmental Neuropsychology, University of Groningen, Groningen, Netherlands

Evaluating the state of the oculomotor system of a patient is one of the fundamental tests done in neuro-ophthalmology. However, up to date, very few quantitative standardized tests of eye movements' quality exist, limiting this assessment to confrontational tests reliant on subjective interpretation. Furthermore, quantitative tests relying on eye movement properties, such as pursuit gain and saccade dynamics are often insufficient to capture the complexity of the underlying disorders and are often (too) long and tiring. In this study, we present SONDA (Standardized Oculomotor and Neurological Disorder Assessment): this test is based on analyzing eye tracking recorded during a short and intuitive continuous tracking task. We tested patients affected by Multiple Sclerosis (MS) and Parkinson's Disease (PD) and find that: (1) the saccadic dynamics of the main sequence alone are not sufficient to separate patients from healthy controls; (2) the combination of spatio-temporal and statistical properties of saccades and saccadic dynamics enables an identification of oculomotor abnormalities in both MS and PD patients. We conclude that SONDA constitutes a powerful screening tool that allows an in-depth evaluation of (deviant) oculomotor behavior in a few minutes of non-invasive testing.

Introduction

Eye movements are a fundamental component of vision. Their evaluation is an important aspect of assessments in neurology and ophthalmology. Two neurological conditions in which oculomotor assessment has great clinical relevance are Multiple Sclerosis (MS) and Parkinson's Disease (PD). In MS, the oculomotor assessment is relevant for diagnosis, for monitoring progress and for prognosis (13), while in Parkinson's Disease (PD), eye-movements can help in differentiating the “pure” form of PD from other forms of parkinsonism (46).

For both MS and PD, a confrontational assessment is the most commonly used technique for evaluating patients' eye motility. While the simplicity of a confrontational assessment is attractive, its qualitative and subjective nature may leave subtle abnormalities undetected. In cases like sub-clinical manifestations of internuclear ophthalmoplegia in MS, or in PD, which typically exhibit eye movement disorders that are also present in the healthy elderly population, quantitative approaches are crucially needed (7).

At present, quantitative assessments are relatively sparse and lack systematic, broadly accepted approaches. In pathological conditions, the oculomotor performance is usually parameterized using measures like the main-sequence [i.e., the relationship between saccadic amplitude and peak velocity (8, 9)], pursuit gain and onset latency, or generic positional errors, measured discretely by repeating many trials (10). These trial-based approaches are often complex, repetitive and time-consuming, and—although being very suitable for research [for a noticeable example see DEMoNS (11)]—they cannot be easily translated into daily clinical practice. Stressing the patients with such tests for longer than a few minutes can also exacerbate their symptoms: MS patients, for instance, are particularly affected by fatigue, which is reflected in changes in their saccadic peak velocity (12).

If we were to have a fast and simple yet powerful assessment of oculomotor quality, this could greatly benefit clinical practice. In this study we introduce SONDA, a Standardized Oculomotor and Neurological Disorders Assessment. This paradigm is partially based on the eye-movement cross-correlogram (13) applied to a continuous visual tracking task. This method measures several spatio-temporal properties (STP) of eye movements while patients track a moving dot with their eyes. In the analyses, the eye movements and their STP can be related to various aspects of visual functioning: e.g., the presence of visual field defects (14) or spatial (15). The main advantage of this approach is that it provides a broad range of oculomotor measurements in a time-efficient and patient-friendly manner. Additionally, being able to provide an in-depth characterization of multiple aspects of eye movements is critical, because the clinical manifestations of oculomotor abnormalities are complex and variegated.

In MS, inflammatory demyelinating lesions in different brain areas result in a wide range of oculomotor disorders: most commonly, static and dynamic ocular misalignment (16) coupled with dysmetric saccadic behavior (17). Furthermore, fixation (18), smooth pursuit (19), and vestibulo-ocular responses (20) are also often impaired. Most of these disorders are accentuated by the presence of internuclear ophthalmoplegia (INO). INO is a neuro-ophthalmic condition present in approximately one out of three MS patients (21) and is characterized by impaired adduction of conjugate lateral eye movements. The presence of INO is often determinant to confirm the diagnosis of MS, especially if it is bilateral (22), but its sub-clinical manifestations are difficult to identify without quantitative approaches (23).

Parkinson's Disease (PD) and other forms of parkinsonism share some common features, such as saccadic intrusions and square-wave jerks, but each clinical phenotype of parkinsonian disorders has a distinctive profile of oculomotor impairments (46). Eye-movement disorders typical of PD are impairment of self-paced saccades, saccadic hypometria, and saccadic fragmentation (24): all contribute to disrupting smooth pursuit as well, due to the resulting inadequacy of catch-up saccades. These disorders, with the exception of the impairment of self-paced saccades, are also commonly present in otherwise healthy elderly individuals (25). This implies that just finding an abnormality in eye movements will be insufficient to disentangle PD-related ones from those caused by normal aging, thus highlighting the need for multidimensional analysis of the oculomotor assessments.

We investigated whether the STP quantified with continuous gaze tracking is a sufficiently sensitive measure of the oculomotor disorders usually present in MS and PD. Furthermore, we extended the continuous gaze tracking method in order to measure saccadic dynamics continuously over time and to extract their frequency distribution. This will provide both a direct comparison between the eye-movement cross-correlogram STP and the saccadic main-sequence, as well as allowing the possibility of integrating these measurements into a single assessment.

More specifically, we answer the following research question: how well can the continuous gaze-tracking task of SONDA identify oculomotor abnormalities associated with MS and PD?

To summarize our findings, in both MS and PD patients, SONDA identified both marked and subtle abnormalities in oculomotor behavior based on only a few minutes of non-invasive testing. Therefore, we conclude that SONDA allows an in-depth characterization of (deviant) oculomotor behavior and thus constitutes a powerful screening tool.

Methods

Observers

We assessed 71 participants divided into three groups. The first group comprised 50 healthy controls (age range 30–79, mean age 53 ± 13.7, 26 females). The second group comprised 12 MS patients (age range 40–76, mean age 54 ± 10.8, 8 females). The third group comprised 9 PD patients (age range 59–77, mean age 68 ± 6.8, 2 females). All observers but one control and one MS patient had normal or corrected to normal visual acuity (below or equal to 0.10 LogMAR). All control participants were tested and screened with questionnaires to assess their familiarity with neurological and ophthalmological disorders at Royal Dutch Visio, Haren, The Netherlands. All observers gave their written informed consent prior to participation. The study was approved by the Medical Ethical Committee of the University Medical Center Groningen and the Ethics Committee of Psychology of the University of Groningen. The study followed the tenets of the Declaration of Helsinki.

Apparatus

A “Tobii T60 XL” eye-tracker (Tobii Technology, Stockholm, Sweden) was used to record the eye movements of the observers. The recording was conducted binocularly, with a sampling frequency of 60 Hz. A custom-made 9-point calibration procedure was performed prior to each experimental session. The calibration was repeated until the average error was below 1° of visual angle and the maximum error below 2.5°. Stimuli were generated using custom-made scripts in MATLAB using the Psychophysics Toolbox (26, 27) and displayed on the integrated screen of the Tobii T60 XL at a refresh rate of 60 Hz and viewed from a distance of 60 cm. Normally, this relatively low refresh rate is less suitable for the study of saccade dynamics, but we addressed this issue using a custom-made saccadic detection algorithm (see “Eye-tracking data analysis” paragraph). Head movements were minimized using a chin-rest with forehead support. All data were analyzed with custom-made MATLAB scripts (R2017b, Mathworks). All display and analysis code is available upon request.

Continuous Tracking Task

The visual stimulus consists of a Gaussian blob of increased luminance moving along a random-walk trajectory on a uniform gray background (~140 cd/m2). The Gaussian blob is displayed at two contrast levels: at high contrast (50%) it has a peak luminance of ~385 cd/m2, while when presented at low contrast (10%), it has a peak luminance of ~160 cd/m2. Its full-width-at-half-maximum is 0.83°, roughly corresponding to the size III of a Goldman perimeter's stimulus.

There were two stimulus conditions: in the smooth pursuit condition, the stimulus moved continuously along the random-walk path, while in the saccadic pursuit condition, an additional positional displacement to a random location on the screen was added to the trajectory, that occurred every 2 s. Each observer performed 6 trials of 20 s for each one of the conditions for a total of 12 trials and 4 min of total data acquisition. Short breaks in between trials were allowed. Additional detail regarding the visual stimulation can be found in Supplementary Material.

Eye-Tracking Data Analysis

The extraction of the oculomotor spatio-temporal features is described in Supplementary Material.

Individual saccades were identified from the raw data using a customized MATLAB implementation based on the algorithm described by Behrens et al. (28). This algorithm allows saccade identification in eye movement time series using an adaptive threshold based on eye-movement acceleration. In short, the algorithm classifies any given time point at which the eye acceleration exceeds this adaptive threshold as being part of a saccadic movement. In their paper, the authors computed the threshold as a constant (K = 3.4) multiplied by the standard deviation of the acceleration distribution within a moving window of 200 frames, acquired at 1 kHz (i.e., a window of 0.2 s). Our adapted version uses a moving window of 60 frames acquired at 60 Hz (i.e., a window of 1 s). Empirically, we found our variant to be adequate for an accurate identification of saccades in our experiment. In contrast to Beherens et al., we did not apply a low-pass filter to the velocity signals.

Since we conducted the gaze recording binocularly, we applied the saccade detection algorithm to the time series of each eye. An example of the saccades detected for one eye during one tracking trial is shown in Figure 1. Once the individual saccades are identified, their amplitude, peak velocity, and direction are computed. Next, we used these parameters to compute the saccadic frequency and the saccadic dynamics distributions as a function of visual field eccentricity and direction. To determine whether the oculomotor behavior of patients differed from normal, we computed normative values based on the results of the control population. Because the age distribution of PD and MS patients differed, we computed separate normative values for each class of patients (see Results section, Figures 5B,C, 6B,C).

FIGURE 1
www.frontiersin.org

Figure 1. (A) Example of a time series of eye movements' positions recorded during a trial of 20 s of the tracking task. (B) Example of eye acceleration and the value of the adaptive saccadic threshold for one eye plotted as a function of time. The adaptive threshold (in red) at any given time point is computed as the standard deviation of acceleration distribution in the preceding 60 samples, multiplied by a constant (K = 3.4). Whenever acceleration exceeds the momentary threshold, the time point is marked as being part of a saccade. (C) Back-projection of the time points marked as being part of a saccade onto the original time series of eye positions (top figure showing the horizontal and the bottom figure showing the vertical component).

Statistical Analysis

The results of the patients' oculomotor behavior are expressed in modified z-scores computed as

zmodified=x-MED1.486·MAD    (1)

Where MED and MAD are the Median and Mean Absolute Deviation in the age-matched normative control population used as a reference. The individual saccades extracted with the algorithm described in the previous section, are binned in the function of their direction along the visual field, and used to display the saccadic frequency distribution (examples in Figures 5B, 6B) and the saccadic dynamics (amplitude and peak velocity, examples in Figures 5C, 6C). Specifically, we treated each saccade of each eye as a vector defined by the starting point and ending point of the saccades, expressed in polar coordinates. We computed the directions of the saccade vectors, and then binned them in 12 main directions, resulting in a bin size of 30°.

For MS patients, the entire group of controls is used as normative data, as their age distributions are very similar (see Observers section). Since as a group the PD patients were significantly older, for the age-matched normalization it was necessary to create a subgroup of controls excluding all participants under 50 years of age. To take into account the different sizes of the groups (Ncontrols = 50, NMS = 12, NPD = 9), the statistical comparison between patients and controls have been performed using one-tailed Welch's t-test for unequal variances with α = 0.05.

The analysis of oculomotor spatio-temporal properties (STP) resulted in 80 features (10 unique STP × 2 axis × 2 pursuit conditions × 2 eyes). See Supplementary Material for a description of the individual spatio-temporal properties. For the purposes of the clustering analysis, we did not normalize the patients' features with respect to the control group. Prior to performing the clustering analysis, we reduced the dimensionality of this dataset by employing a Principal Component Analysis (PCA) with the goal of removing redundant features. Based on the PCA, we identified the number of components sufficient to explain at least 95% of the variance. Following the PCA, we applied t-Distributed Stochastic Neighbor Embedding (tSNE) on the remaining components. This technique allows visualizing high-dimensional data by modeling each high-dimensional object as a two- (or three-) dimensional point in such a way that similar objects are modeled as nearby points and dissimilar objects are modeled as more distant points. Lastly, to identify clusters amongst the t-SNE points, we applied a k-means clustering algorithm. The appropriate number of clusters k was chosen using the Elbow Method (29) applied to the Within-clusters Sum of Squares (WSS). Based on the resulting k, the location and size of each of the clusters was computed.

Results

To summarize our findings, we find that the main-sequence alone is often insufficient to detect the presence of oculomotor abnormalities in either MS or PD. In contrast, we find that the STP obtained with our SONDA paradigm well separates PD patients from healthy controls. The STP of MS patients are insufficient to separate them from controls. Nevertheless, when we combine the STP with saccadic dynamics and frequency distributions, both PD and MS patients can be distinguished from controls.

Main-Sequences of Saccades in PD and MS Do Not Differ From Controls

First, we determined whether the oculomotor abnormalities of PD and MS patients are reflected in changes in the measures obtained with our continuous gaze-tracking task. Figure 2 shows the main-sequence plots of amplitude vs. peak velocity for each saccade made by every patient and control made during the tracking task (algorithm described in Methods, Figure 1). We found no differences at the group level: both the amplitude and peak velocity distributions of the two patient groups overlap with those of the control group. For all three groups tested, we find an exponential relationship between amplitude and peak velocity. This behavior is expected (8). The prediction intervals (95%) of the patients' group data overlap with those computed for the control groups' data.

FIGURE 2
www.frontiersin.org

Figure 2. Saccadic main-sequences in the three examined groups of observers. The individual saccades (scatter plot) display the classic exponential relationship between amplitude and peak velocity, with no statistically significant differences between Control, MS, and PD groups: all 95% prediction intervals overlap with each other. Also when taken separately (marginal histograms), the amplitude and peak velocity distributions do not differ between groups or between eyes.

Eye Movement Spatio-Temporal Properties of PD and MS Compared to Normative Data

Figure 3 shows the results of the analysis of the STP of eye movements made in the continuous tracking task.

FIGURE 3
www.frontiersin.org

Figure 3. Oculomotor spatio-temporal properties in MS and PD cohorts, compared to an age-matched normative baseline. For the spatio-temporal analysis of the continuous gaze tracking task, we applied the eye-movement crosscorrelogram method (13, 14). From all features computed, we selected the two temporal and two spatial features that were shown to be most strongly modulated by stimulus velocity (30). The selected features are: “lag” (F1), “temporal uncertainty” (F2), “error spread” (F3), and “dissimilarity” (F4). “Lag” represents the delay between eyes and stimulus, measured as the temporal offset of the peak of the cross-correlation between stimulus and eye velocities. “Temporal uncertainty” represents the consistency of the velocity of the eye in relation to the velocity of the stimulus. “Error spread” is the standard deviation of the distribution of eye-stimulus positional deviations. “Dissimilarity” is the inverse of the cosine similarity computed between the arrays of gaze coordinates and stimulus coordinates.

In MS patients (upper row), we find a substantial overlap of their features with the normative ones. We can observe a higher variability compared to controls and PD (lower row), especially in the Smooth condition. Statistically significant differences (indicated by an asterisk) are found for the temporal uncertainty and the movement dissimilarity features, but only in the Saccadic condition, where participants had to make saccades in order to track the rapid, random, displacements of the target.

In PD patients (lower row), nearly all features are significantly different from the normative ones. The magnitude and variance of the differences are larger in the Saccadic condition, similar to what we found in the MS patients.

For both MS and PD patients, there was neither significant difference for the left and right eyes, nor for the horizontal and vertical components of the eye movements.

Based on the amount of significantly deviating features, our results indicate that, at the group level, MS patients had abnormalities primarily during the saccadic pursuit, while PD patients had abnormalities in both saccadic and smooth pursuit conditions.

Classification of the Neurological Disorder Based on Oculomotor Abnormalities

We examined whether it is possible to identify the underlying neurological conditions of individual patients based exclusively on the eye-movements measurement. The results are shown in Figure 4. Of the 80 initial features (10 oculomotor properties × 2 axis × 2 pursuit conditions × 2 eyes), the PCA indicated that 11 components are sufficient to explain 95% of the variance in the data (Figure 4A). We then used the t-SNE algorithm to embed these 11 principal components into two final components for clustering and visual representation. We found that the optimal number of clusters to represent the data is five (Figure 4B). The results of the k-mean analysis are shown in Figure 4C. The leftmost cluster (1) contains nearly all PD patients and two of the MS patients. Being most separated from the other clusters, this cluster appears to group the cases with the most severe oculomotor impairments. The rightmost cluster (2) includes the vast majority of control subjects, thus grouping those without an oculomotor impairment. It partially overlaps with clusters (3) and (4) which contain participants from various groups and together cover nearly all MS patients. These clusters, therefore, appear to group patients with a suspected impairment of different degrees of severity. The remaining cluster (5) comprises outliers that, upon verification, were associated with very noisy data and/or exceptionally poor performance. Altogether, these results show that despite the relatively low number of subjects, a classification based exclusively on STP is reasonably effective for Parkinson's Disease but ineffective for MS.

FIGURE 4
www.frontiersin.org

Figure 4. Results of STP-based clustering analysis for the entire tested population. (A) Principle Component Analysis (PCA): all the spatio-temporal features available from SONDA (see Supplementary Material) were used as an input for the PCA. The resulting components were processed with t-SNE to represent the high-dimensionality dataset into a lower-dimensionality space for clustering purposes. (B) Computation of the optimal number of clusters using the “elbow method” with the Within-clusters Sum of Squares (WSS) as a parameter. (C) Result of the k-means clustering algorithm (k = 5) applied to the outcome components of t-SNE.

Combination of STP With Statistical and Dynamic Properties of Saccades

Next, we explored whether taking into account the frequency distribution of saccades and their dynamic properties as a function of direction (i.e., visual field polar angle) in combination with STP obtained from the saccadic pursuit condition, would allow for a better characterization of the oculomotor disorders. The results for the left and right eye are displayed separately.

Figures 5, 6 show the results for MS and PD groups, respectively. At a group level, MS patients showed impaired STP (Figure 5A), no significant difference in the saccadic frequency distribution (Figure 5B) and significantly different amplitude and velocity (Figures 5C,D), primarily along the horizontal axis. Still at the group level, PD patients showed impaired STP (Figure 6A), significantly reduced saccadic frequency along the vertical axis (Figure 6B), and significantly different amplitude and velocity (Figures 6C,D). Overall the PD group showed more significant deviations from controls compared to the MS group.

FIGURE 5
www.frontiersin.org

Figure 5. (A) MS and PD group analysis of normalized spatio-temporal properties. (B) Normalized saccadic frequency distribution as a function of direction. (C) Normalized saccadic amplitude as a function of direction. (D) Normalized saccadic peak velocity as a function of direction. The results of the saccadic dynamics are stratified per quartile. An asterisk indicates a significant difference between normative values obtained from the control group.

FIGURE 6
www.frontiersin.org

Figure 6. Same as Figure 5.

We then performed the same analysis, but at the single subject level. Figures 7, 8 show the results for three MS and three PD patients, respectively, and two age-matched control observers. From each patient group, we selected three patients of interest that showed either mixed results regarding the STP, the saccadic frequency and saccadic dynamics, or that had an additional diagnosis (e.g., INO in addition to MS).

FIGURE 7
www.frontiersin.org

Figure 7. An in-depth overview of the oculomotor performance of three MS patients, three PD patients, and two age-matched controls. (A) Normalized spatio-temporal properties (STP). We consider an STP feature abnormal if its deviates 2 or more standard deviations from the normative value (indicated in green). (B) Normalized saccadic frequency distributions. (C) Normalized saccadic dynamics (amplitude and peak velocity), as a function of saccade direction. The results are stratified per quartile, separately for amplitude and peak velocity. Each quartile has its own normative boundaries.

FIGURE 8
www.frontiersin.org

Figure 8. Same as Figure 7.

Panel A shows part of the STP previously presented in Figure 3. We present the features for the horizontal and vertical components of the saccadic pursuit condition, which we already showed to be the condition with the most significant differences. Our results show that the STP and the saccadic dynamics are not always in agreement and thus convey different types of information that may be relevant for a clinical interpretation. Patient MS#1 had been previously diagnosed with internuclear ophthalmoplegia (INO) and has severely impaired STP (Figure 7A), a dysmetric saccadic behavior with reduced adduction frequency (Figure 7B), and an abnormal range of saccadic amplitudes and peak velocities of which both are hypermetric (Figure 7C). Patient MS#2 has normal STP (Figure 7A), an anisotropic saccadic frequency distribution skewed rightward and upward binocularly (Figure 7B). The saccades made by both eyes are in the normal velocity range, but their amplitude is mildly hypermetric, especially along a near vertical and in a horizontal rightward direction (Figure 7C). Patient MS#3 has moderately impaired STP (Figure 7A), a normal saccade frequency distribution (Figure 7B), and hypermetric amplitude of the vertical saccades, with a normal peak velocity for both eyes (Figure 7C).

Patient PD#1 has impaired STP, with a strong oculomotor lag in both horizontal and vertical components (Figure 8A), while both their saccadic frequency distribution (Figure 8B) and saccadic dynamics (Figure 8C) are within the normative range. Patient PD#2 shows severely impaired STP (Figure 8A), an unusual “diagonally skewed” saccadic frequency distribution for both eyes (with the upward component slanted toward the nasal side; Figure 8B) and the saccadic peak velocity within the normative range and nearly perfect normal saccadic amplitude (Figure 8C). Finally, patient PD#3 has normal STP (Figure 8A) and shows saccadic frequency distributions that deviate from normal (Figure 8B) and show slightly hypermetric vertical amplitude of the diagonal saccades. The peak velocities for both eyes were normal (Figure 8C).

Discussion

The main contribution of the present study is the introduction of a novel method to perform a time-efficient (<5 min), non-invasive, neuro-ophthalmic screening that is based on an in-depth characterization of oculomotor behavior in response to a simple visual tracking task.

First, our results indicate that the most commonly used conventional measures, such as saccadic amplitude and peak velocity by themselves, are often insufficient to fully characterize the oculomotor behavior of PD and MS patients. Moreover, we show that the computation of the spatio-temporal properties (STP) of eye-movements recorded during a continuous tracking task can help detect abnormalities in the oculomotor behavior. Therefore, we conclude that a combination of improved conventional (saccadic frequency distribution and dynamics) and our entirely new spatio-temporal analyses are an effective screening tool to comprehensively assess the oculomotor behavior of a patient. The clinical relevance of this will be discussed below for each class of patients.

SONDA Identifies Oculomotor Abnormalities in Multiple Sclerosis

The parameters of the main-sequences—saccadic amplitude and saccadic peak velocity—of MS patients showed an almost identical distribution to those of the controls (see Figure 2). This is consistent with the literature, which indicates that most CNS lesions of this disease have not been associated with unique identifiable clinical symptoms in eye-movements (2). Thus far, the majority of the studies on oculomotor abnormalities in MS patients focused on the detection of internuclear ophthalmoplegia (INO). This disorder may be present (often subclinically) in ~30% of MS patients and is usually identified by measuring the versional dysconjugacy index (VDI), which is the ratio between abducting and adducting eye saccadic dynamics (21, 31, 32). However, despite that relying on the main-sequence parameters alone is insufficient, standardized screening tests capable of detecting other manifestations of oculomotor abnormalities in MS do not exist.

With our novel “SONDA” approach introduced in this study, in MS we found preserved smooth pursuit responses and two abnormal features in the saccadic pursuit condition (see Figure 3). This is consistent with the notion that the most common saccadic issue in MS is saccadic dysmetria (17), present especially in patients with demyelination affecting the cerebellar peduncles (33). However, given the relatively small sample (NMS = 12) the significant difference seems to be driven mostly by a few severely impaired patients, rather than a general decrease in performance.

An interesting result emerged from the clustering analysis shown in Figure 4. In this case, we performed a multivariate analysis based on the whole feature-space of STPs obtained with SONDA, rather than selecting a few parameters a priori. We found that MS patients, although not well separable from controls, nevertheless tended to cluster [clusters (3) and (4)] based solely on their STPs. This suggests the presence of subtle oculomotor abnormalities that require more detailed quantification of the eye movements.

To further explore this hypothesis we also quantified per eye the saccadic frequency and the saccadic dynamics distributions in the function of direction (Figures 58).

These extended analyses revealed the presence of abnormalities even in patients falling within the normative STP ranges, or in the saccadic dynamic when taking into account their directionality. While at a group level, few significant deviations are detected (see Figures 5B–D), at the single subject level taking into account the direction of the saccades can aid in identifying the underlying oculomotor disorder. For instance, patient MS#2 has a normal STP (Figure 7A), only slightly abnormal saccadic dynamics (Figure 7C), but a very peculiar saccadic frequency distribution (Figure 7B). This latter analysis revealed an abnormal prevalence of horizontal rightward saccades for both eyes. This is analogous to what has been observed in pathological nystagmus: a slow-drift of eccentric gaze in one direction followed by a resetting saccade in the other, which is caused by a failure of the common neural integrator of eye movements in the brainstem (2). While the other two MS patients (MS#1 and MS#3) have an abnormal STP (Figure 7A). A further inspection of their saccadic frequency distributions does differentiate between them: MS#3 has a normal saccadic frequency distribution (Figure 7B) and only slightly hypermetric vertical saccades (Figure 7C). The oculomotor abnormalities of patient MS#3 might have gone unnoticed if the assessment had been limited to the saccadic dynamics. On the other hand, besides an abnormal STP, patient MS#1 also displays overall hypermetric saccades, as well as the suppression of adducting saccades (Figure 7B) that is typical of INO.

SONDA Identifies Oculomotor Abnormalities in Parkinson's Disease

In PD, slower saccades should be expected in advanced cases (25, 34, 35). Yet, at the group level, the prediction intervals for the saccadic main-sequence of the PD patients overlap with those of the control group (see Figure 2). In contrast a significant difference is found in all examined STP, except for the spread of the positional error (Figure 3). Noticeably, we found marked differences along the vertical axis of the smooth pursuit condition for almost the entire group, although mild cases of PD are expected to be similar to age-matched control observers (25, 36). This can be due to the fact that while unpredictable pursuit should be poor in PD, the predictive component of pursuit is generally preserved (37). However, our stimulus trajectory is completely unpredictable (by design it does not contain any autocorrelation; see Supplementary Material). This explains why our approach was sensitive to PD related abnormalities.

Regarding saccades, one of the characteristic oculomotor disorders in PD is the impairment of self-generated saccades with the relative sparing of visually-driven saccades, which only deteriorate as the disease progresses (10). In our tracking task we tested only visually-driven saccades but still found significant differences in the STP of most patients. This is important, because at the single-subject level (Figure 8C) the saccadic dynamics of the PD in some patients were remarkably normal. It suggests that SONDA is a sensitive paradigm.

For the PD patients, the clustering analysis based on SONDA's STP resulted in the separation of most PD individuals from the controls (Figure 4). This is quite remarkable considering that, for pure Parkinson's Disease (which all patients in our group were) the eye movements usually only show minor disturbances that are common in healthy elderly subjects as well (10). Note that this contrasts with other forms of parkinsonism. Furthermore, one of the main oculomotor abnormalities in PD is the difficulty in generating internally driven saccades. Most errors are made when PD patients are required to switch between instructions [e.g., switching between pro- and antisaccades (3840)]. Such switching is not demanded by our paradigm. Nonetheless, abnormal oculomotor behavior could be clearly identified in all PD patients.

An in-depth inspection of the oculomotor behavior that takes into account the direction of saccades, revealed that at a group level some significant deviations in amplitude and peak velocity can be detected (Figures 6C,D). Further more the single-subject analysis (Figure 8) showed that, analogously to MS, an analysis of the saccadic frequency distribution can reveal subtle abnormalities that other measures (STP and saccadic dynamics) cannot (see Figure 8, PD#3 for an example). Nevertheless, in the majority of the PD patients, the STP by itself was sufficient to identify eye movement abnormality (as already revealed by the clustering analysis).

SONDA Is Clinically Relevant

The SONDA paradigm allows for a comprehensive examination of eye-movements that provides detailed information on the dynamics and statistical properties of saccades, as well as a novel set of spatio-temporal properties that quantify both smooth and saccadic pursuit behavior. Given the large number of parameters, our technique also opens up possibilities for machine-learning-aided diagnosis of neuro-ophthalmic and neurological impairment based on eye movements. We and others have previously demonstrated the feasibility of such an approach for identifying (simulated) ophthalmic disorders (14, 41), and for neurological conditions [the latter based on participants viewing natural scenes (42)]. Compared to other high-throughput eye-tracking tests, SONDA has two main advantages. First, excluding calibration (that is common for all eye-tracking tests), the whole assessment takes only 4 min [compared to ~20 min for the approach of Tseng et al. (42) and ~25 min for DEMoNS (11)]. Second, our task and stimulus were designed to be as intuitive and simple as possible to avoid confounding factors that may be caused by natural or more complex stimuli, or complicated instructions resulting in a higher cognitive load.

Future Applications

Objective, quantitative eye-movement-based tests are increasingly becoming an important aid in diagnosing and monitoring neuro-ophthalmic conditions of which the oculomotor manifestations cannot always be clinically detected with confrontational assessments (43). For instance, in some patients oscillopsia occurs in conjunction with eye movements that are too small to be detected using an ophthalmoscope, yet that are large enough to cause visual complaints (44). Another example of the use of eye-movement-based metrics is the assessment of surgical or pharmacological interventions for nystagmus (45, 46). Such disorders are rather common in MS when lesions involve the brainstem and the cerebellum (2) and can occur in PD as well, sometimes causing oscillopsia (47).

Besides MS and PD, the new parameters provided by SONDA can be relevant in the evaluation of other neurological and psychiatric disorders, such as mild traumatic brain injury (48), schizophrenia (49), or dementia with Lewy bodies (50). In these disorders, the parameters of the saccadic main-sequence of patients often overlap with those of controls. For this reason, outcomes cannot be considered definitively diagnostic (17). A multi-dimensional characterization of oculomotor properties, such as provided by SONDA, could thus also be beneficial in these conditions.

Limitations

Our present study still has some limitations that should be addressed in future follow-up. First, the number of observers tested in comparison to the number of features analyzed was only just enough to use t-SNE and k-mean clustering (N = 71 observers and 11 PCA components; see Figure 4). Combined with the imbalance within groups (Ncontrol = 50, NMS = 12, NPD = 9) this may have prevented an optimal separation of the MS patients from the controls. Another aspect that requires further investigation is the sensitivity/specificity of the measures that were expressed as a function of the visual field polar angle (Figures 7B,C, 8B,C). The number of bins used to determine the polar histograms was chosen arbitrarily and might affect the false negative/positive rate depending on the data available in each bin.

Conclusions

We presented SONDA, a novel paradigm to assess oculomotor abnormalities that can be linked to neurological disorders. We tested patients affected by Multiple Sclerosis and Parkinson's Disease and showed that: (1) the parameters of the saccadic main-sequence alone are insufficient to separate patients from healthy controls; (2) spatio-temporal and statistical properties of the eye movements obtained with SONDA, when combined with detailed analysis of saccadic dynamics allow for a correct identification of the oculomotor abnormalities and a possible identification of the underlying neurological disorder. Our test constitutes an attractive approach for quantitative neuro-ophthalmic screening, being non-invasive, intuitive and fast to administer.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by Medical Ethical Committee of the University Medical Center Groningen and the Ethics Committee of Psychology of the University of Groningen. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

AG, AV, and FC designed the experiment. AG and RR analyzed the data. AG and FC wrote the manuscript. AV, JH, and FC wrote the grant. All authors reviewed the manuscript.

Funding

This project has received funding from the Programmaraad Visuele Sector (grant agreement STIP-MS), the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie (grant agreement No. 641805 NextGenVis), and the Graduate School Medical Sciences (GSMS) of the University Medical Center Groningen.

Conflict of Interest

AG and RR are listed as inventors on the European patent application (51) which is partially based on the content of the Supplementary Material of this manuscript. AG is currently a majority shareholder of Reperio B.V. (Dutch chamber of commerce registration: 76175596), but was not at the time the research was conducted.

The remaining 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.

Acknowledgments

We are thankful to Nadine Naumann and Teatske Hoekstra from Royal Dutch Visio for their support in assessing patients. We are thankful to Rijul S. Soans for his help in setting up the customized calibration procedure of the Tobii eye-tracker.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2020.00971/full#supplementary-material

References

1. Derwenskus J, Rucker JC, Serra A, Stahl JS, Downey DL, Adams NL, et al. Abnormal eye movements predict disability in MS: two-year follow-up. Ann N Y Acad Sci. (2005) 1039:521–3. doi: 10.1196/annals.1325.058

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Frohman EM, Frohman TC, Zee DS, McColl R, Galetta S. The neuro-ophthalmology of multiple sclerosis. Lancet Neurol. (2005) 4:111–21. doi: 10.1016/S1474-4422(05)00992-0

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Graves J, Balcer LJ. Eye disorders in patients with multiple sclerosis: natural history and management. Clin Ophthalmol. (2010) 4:1409–22. doi: 10.2147/OPTH.S6383

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Abadi RV, Gowen E. Characteristics of saccadic intrusions. Vision Res. (2004) 44:2675–90. doi: 10.1016/j.visres.2004.05.009

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Pretegiani E, Optican LM. Eye movements in Parkinson's disease and inherited Parkinsonian syndromes. Front Neurol. (2017) 8:592. doi: 10.3389/fneur.2017.00592

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Rascol O, Sabatini U, Simonetta-Moreau M, Montastruc JL, Rascol A, Clanet M. Square wave jerks in Parkinsonian syndromes. J Neurol Neurosurg Psychiatry. (1991) 54:599–602. doi: 10.1136/jnnp.54.7.599

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Marandi RZ, Gazerani P. Aging and eye tracking: in the quest for objective biomarkers. Fut Neurol. (2019) 14:FNL33. doi: 10.2217/fnl-2019-0012

CrossRef Full Text | Google Scholar

8. Bahill AT, Terry Bahill A, Clark MR, Stark L. The main sequence, a tool for studying human eye movements. Math Biosci. (1975) 24:191–204. doi: 10.1016/0025-5564(75)90075-9

CrossRef Full Text | Google Scholar

9. Collewijn H, Erkelens CJ, Steinman RM. Binocular co-ordination of human horizontal saccadic eye movements. J Physiol. (1988) 404:157–82. doi: 10.1113/jphysiol.1988.sp017284

CrossRef Full Text | Google Scholar

10. Leigh RJ, Zee DS. The Neurology of Eye Movements. Oxford, UK: Oxford University Press (2015).

Google Scholar

11. Nij Bijvank JA, Petzold A, Balk LJ, Tan HS, Uitdehaag BMJ, Theodorou M, et al. A standardized protocol for quantification of saccadic eye movements: DEMoNS. PLoS ONE. (2018) 13:e0200695. doi: 10.1371/journal.pone.0200695

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Finke C, Pech LM, Sömmer C, Schlichting J, Stricker S, Endres M, et al. Dynamics of saccade parameters in multiple sclerosis patients with fatigue. J Neurol. (2012) 259:2656–63. doi: 10.1007/s00415-012-6565-8

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Mulligan JB, Stevenson SB, Cormack LK. Reflexive and voluntary control of smooth eye movements. In: Rogowitz BE, Pappas, TN, de Ridder H, editors. Human Vision and Electronic Imaging XVIII. Bellingham, WA: SPIE (2013). doi: 10.1117/12.2010333

CrossRef Full Text | Google Scholar

14. Grillini A, Ombelet D, Soans RS, Cornelissen FW. Towards using the spatio-temporal properties of eye movements to classify visual field defects. In: Proceedings of the 2018 ACM Symposium on Eye Tracking Research & Applications–ETRA '18. New York, NY (2018). doi: 10.1145/3204493.3204590

CrossRef Full Text | Google Scholar

15. Bonnen K, Burge J, Yates J, Pillow J, Cormack LK. Continuous psychophysics: target-tracking to measure visual sensitivity. J Vision. (2015) 15:14. doi: 10.1167/15.3.14

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Serra A, Chisari CG, Matta M. Eye movement abnormalities in multiple sclerosis: pathogenesis, modeling, and treatment. Front Neurol. (2018) 9:31. doi: 10.3389/fneur.2018.00031

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Serra A, Derwenskus J, Downey DL, Leigh RJ. Role of eye movement examination and subjective visual vertical in clinical evaluation of multiple sclerosis. J Neurol. (2003) 250:569–75. doi: 10.1007/s00415-003-1038-8

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Mallery RM, Poolman P, Thurtell MJ, Full JM, Ledolter J, Kimbrough D, et al. Visual fixation instability in multiple sclerosis measured using SLO-OCT. Invest Ophthalmol Visual Sci. (2018) 59:196–201. doi: 10.1167/iovs.17-22391

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Lizak N, Clough M, Millist L, Kalincik T, White OB, Fielding J. Impairment of smooth pursuit as a marker of early multiple sclerosis. Front Neurol. (2016) 7:206. doi: 10.3389/fneur.2016.00206

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Huygen PLM, Theunissen EJJ, Hommes OR. The vestibulo-ocular reflex in multiple sclerosis. In: Hommes OR, editor. Multiple Sclerosis Research in Europe. Dordrecht: Springer (1986). p. 363–72. doi: 10.1007/978-94-009-4143-4_52

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Jozefowicz-Korczynska M, Łukomski M, Pajor A. Identification of internuclear ophthalmoplegia signs in multiple sclerosis patients. J Neurol. (2008) 255:1006–11. doi: 10.1007/s00415-008-0819-5

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Bolanos I, Lozano D, Cantu C. Internuclear ophthalmoplegia: causes and long-term follow-up in 65 patients. Acta Neurol Scand. (2004) 110:161–5. doi: 10.1111/j.1600-0404.2004.00278.x

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Frohman TC, Frohman EM, O'Suilleabhain P, Salter A, Dewey RB Jr, Hogan N, et al. Accuracy of clinical detection of INO in MS: corroboration with quantitative infrared oculography. Neurology. (2003) 61:848–50. doi: 10.1212/01.WNL.0000085863.54218.72

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Winograd-Gurvich C, Georgiou-Karistianis N, Fitzgerald PB, Millist L, White OB. Self-paced saccades and saccades to oddball targets in Parkinson's disease. Brain Res. (2006) 1106:134–41. doi: 10.1016/j.brainres.2006.05.103

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Rottach KG, Riley DE, DiScenna AO, Zivotofsky AZ, Leigh RJ. Dynamic properties of horizontal and vertical eye movements in parkinsonian syndromes. Ann Neurol. (1996) 39:368–77. doi: 10.1002/ana.410390314

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Brainard DH. The psychophysics toolbox. Spatial Vision. (1997) 10:433–6. doi: 10.1163/156856897X00357

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Pelli DG. The VideoToolbox software for visual psychophysics: transforming numbers into movies. Spatial Vision. (1997) 10:437–42. doi: 10.1163/156856897X00366

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Behrens F, Mackeben M, Schröder-Preikschat W. An improved algorithm for automatic detection of saccades in eye movement data and for calculating saccade parameters. Behav Res Methods. (2010) 42:701–8. doi: 10.3758/BRM.42.3.701

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Thorndike RL. Who belongs in the family? Psychometrika. (1953) 18:267–76. doi: 10.1007/BF02289263

CrossRef Full Text | Google Scholar

30. Grillini A, Kromm M, Renken RJ, Cornelissen FW. Motion sensitivity assessment based on an analysis of the spatio-temporal features of eye movements. PsyArXiv [Preprint]. (2020). Available online at: psyarxiv.com/wqzms

Google Scholar

31. Frohman EM, Frohman TC, O'Suilleabhain P, Zhang H, Hawker K, Racke MK, et al. Quantitative oculographic characterisation of internuclear ophthalmoparesis in multiple sclerosis: the versional dysconjugacy index Z score. J Neurol Neurosurg Psychiatry. (2002) 73:51–5. doi: 10.1136/jnnp.73.1.51

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Nij Bijvank JA, van Rijn LJ, Balk LJ, Tan HS, Uitdehaag BMJ, Petzold A. Diagnosing and quantifying a common deficit in multiple sclerosis: Internuclear ophthalmoplegia. Neurology. (2019) 92:e2299–308. doi: 10.1212/WNL.0000000000007499

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Frohman EM, Frohman TC, Fleckenstein J, Racke MK, Hawker K, Kramer PD. Ocular contrapulsion in multiple sclerosis: clinical features and pathophysiological mechanisms. J Neurol Neurosurg Psychiatry. (2001) 70:688–92. doi: 10.1136/jnnp.70.5.688

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Vidailhet M, Rivaud S, Gouider-Khouja N, Pillon B, Bonnet A-M, Gaymard B, et al. Eye movements in parkinsonian syndromes. Ann Neurol. (1994) 35:420–6. doi: 10.1002/ana.410350408

CrossRef Full Text | Google Scholar

35. White OB, Saint-Cyr JA, Sharpe JA. Ocular motor deficits in Parkinson's disease. Brain. (1983) 106:555–70. doi: 10.1093/brain/106.3.555

CrossRef Full Text | Google Scholar

36. Waterston JA, Barnes GR, Grealy MA, Collins S. Abnormalities of smooth eye and head movement control in Parkinson's disease. Ann Neurol. (1996) 39:749–60. doi: 10.1002/ana.410390611

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Gilaie-Dotan S, Saygin AP, Lorenzi LJ, Egan R, Rees G, Behrmann M. The role of human ventral visual cortex in motion perception. Brain. (2013) 136:2784–98. doi: 10.1093/brain/awt214

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Cameron IGM, Pari G, Alahyane N, Brien DC, Coe BC, Stroman PW, et al. Impaired executive function signals in motor brain regions in Parkinson's disease. Neuroimage. (2012) 60:1156–70. doi: 10.1016/j.neuroimage.2012.01.057

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Rivaud-Pechoux S, Vidailhet M, Brandel JP, Gaymard B. Mixing pro- and antisaccades in patients with parkinsonian syndromes. Brain. (2006) 130:256–64. doi: 10.1093/brain/awl315

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Terao Y, Fukuda H, Ugawa Y, Hikosaka O. New perspectives on the pathophysiology of Parkinson's disease as assessed by saccade performance: a clinical review. Clin Neurophysiol. (2013) 124:1491–506. doi: 10.1016/j.clinph.2013.01.021

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Crabb DP, Smith ND, Zhu H. What's on TV? Detecting age-related neurodegenerative eye disease using eye movement scanpaths. Front Aging Neurosci. (2014) 6:312. doi: 10.3389/fnagi.2014.00312

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Tseng P-H, Cameron IGM, Pari G, Reynolds JN, Munoz DP, Itti L. High-throughput classification of clinical populations from natural viewing eye movements. J Neurol. (2013) 260:275–84. doi: 10.1007/s00415-012-6631-2

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Bedell HE, Stevenson SB. Eye movement testing in clinical examination. Vision Res. (2013) 90:32–7. doi: 10.1016/j.visres.2013.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Smith JL. Neuro-Ophthalmology Now! Berlin: Springer (1986).

45. Hertle RW, Dell'Osso LF, FitzGibbon EJ, Thompson D, Yang D, Mellow SD. Horizontal rectus tenotomy in patients with congenital nystagmus: results in 10 adults. Ophthalmology. (2003) 110:2097–105. doi: 10.1016/S0161-6420(03)00802-9

PubMed Abstract | CrossRef Full Text | Google Scholar

46. McLean R, Proudlock F, Thomas S, Degg C, Gottlob I. Congenital nystagmus: randomized, controlled, double-masked trial of memantine/gabapentin. Ann Neurol. (2007) 61:130–8. doi: 10.1002/ana.21065

CrossRef Full Text | Google Scholar

47. Kaski D, Bronstein AM. Ocular tremor in Parkinson's disease: discussion, debate, and controversy. Front Neurol. (2017) 8:134. doi: 10.3389/fneur.2017.00134

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Thiagarajan P, Ciuffreda KJ, Ludlam DP. Vergence dysfunction in mild traumatic brain injury (mTBI): a review. Ophthal Physiol Opt J Br Coll Ophthal Optic. (2011) 31:456–68. doi: 10.1111/j.1475-1313.2011.00831.x

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Holzman PS, Proctor LR, Hughes DW. Eye-tracking patterns in schizophrenia. Science. (1973) 181:179–81. doi: 10.1126/science.181.4095.179

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Mosimann UP, Müri RM, Burn DJ, Felblinger J, O'Brien JT, McKeith IG. Saccadic eye movement changes in Parkinson's disease dementia and dementia with Lewy bodies. Brain. (2005) 128:1267–76. doi: 10.1093/brain/awh484

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Grillini A, Hernández-García A, Renken JR. Method, System and Computer Program Product for Mapping a Visual Field. EP19209204.7 (2019).

Keywords: eye movement, assessment, Parkinson's disease, multiple sclerosis, standardized, clinical, screening

Citation: Grillini A, Renken RJ, Vrijling ACL, Heutink J and Cornelissen FW (2020) Eye Movement Evaluation in Multiple Sclerosis and Parkinson's Disease Using a Standardized Oculomotor and Neuro-Ophthalmic Disorder Assessment (SONDA). Front. Neurol. 11:971. doi: 10.3389/fneur.2020.00971

Received: 07 April 2020; Accepted: 24 July 2020;
Published: 08 September 2020.

Edited by:

Tim Anderson, University of Otago, Christchurch, New Zealand

Reviewed by:

Yasuo Terao, Kyorin University, Japan
Aasef G. Shaikh, Case Western Reserve University, United States

Copyright © 2020 Grillini, Renken, Vrijling, Heutink and Cornelissen. 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: Alessandro Grillini, a.grillini@rug.nl