Open Access
23 September 2019 Distinguishing cognitive impairment by using singularity spectrum and lacunarity analysis of the retinal vascular network
Edmund Arthur, Gabor Mark Somfai, Maja Kostic, Susel Oropesa, Carlos Mendoza-Santiesteban, Delia Cabrera DeBuc
Author Affiliations +
Abstract

The development of effective therapies for cognitive impairment (CI), especially due to Alzheimer’s disease, demands diagnosing the condition during the prodromal phase. The diagnosis of CI involves expensive and invasive methods, such as positron emission tomography and cerebrospinal fluid assessment via spinal tap. Hence, a comparatively lower cost and noninvasive method of diagnosis is imperative. The human retina is an extension of the brain characterized by similarities in vascular and neural structures. The complications of CI are not only limited to the brain but also affect the retina for which the loss of retinal ganglion cells has been associated with neurodegeneration in the brain. The loss of retinal ganglion cells in individuals with CI may be related to reduced vascular demand and a potential remodeling of the retinal vascular branching complexity. Retinal imaging biomarkers may provide a low cost and noninvasive alternative for the diagnosis of CI. In this study, the retinal vascular branching complexity of patients with CI was characterized using the singularity spectrum multifractal dimension and lacunarity parameter. A reduced vascular branching complexity was observed in subjects with CI when compared to age- and sex-matched cognitively healthy controls. Significant associations were also found between retinal vascular and functional parameters.

1.

Introduction

There are approximately 46 million dementia patients worldwide according to the 2015 World Alzheimer Report.1 This number is likely to double every 20 years, and a projected increase to 131.5 million people by 2050 is expected.1 Alzheimer’s disease (AD) is the most common, progressive cause of dementia in the elderly, and a severe burden on the aging society worldwide.1 As one of the top ten leading causes of death in North America, AD has no proven preventive or curative interventions.2 Early- and cost-effective diagnosis is critical to the next stage of treatment and drug development.2 However, progress in the management and treatment of AD is limited by constraints surrounding the early diagnosis of the disease.38 Although there are antiamyloid drugs that can slow down the accumulation of amyloid plaques in the brain and hence the progression of the disease, a difficulty lies in the detection of the disease during the prodromal phase. The accumulation of amyloid plaques in the brain is detected by expensive and invasive methods, such as positron emission tomography (PET) and cerebrospinal fluid assessment via a spinal tap.2 Therefore, a comparatively lower cost and noninvasive method may aid in the diagnosis of the disease.

The human retina is an extension of the brain.9 The retinal ganglion cells are similar to that of the cerebral cortex neurons, and the cerebral small vessels are similar to that of the retinal vessels.9 Also, it is well established that complications of cognitive impairment (CI) are not only limited to the brain but also affect the retina, which is evident in the neurodegeneration in the brain concomitant with the loss of retinal ganglion cells10,11 and reduced bioelectrical activity of the retinal neurons.12,13 Compared with standard neuroimaging techniques, imaging the retina is noninvasive and can be low cost. Therefore, the human retina may serve as a noninvasive window to study and possibly diagnose CI early during the prodromal phase of the disease. The neurodegeneration that occurs in the retina of cognitively impaired individuals may indicate reduced oxygen and nutrient demand as well as a potential remodeling of the retinal vascular network or branching pattern.

The human retinal vascular network holds a self-similarity or scaling property with a characteristic branching pattern.1419 Therefore, the retinal vascular network produces similar patterns of different sizes at different magnification or scale.12,1419 This self-similar property of the human retinal vascular network can be characterized by using a mathematical computational method termed fractal dimension (FD) analysis. A retinal vascular network with a more complex branching pattern has a large FD and vice versa. The retinal vascular branching pattern in patients with CI,12,20,21 diabetes, 2226 and amblyopia27 has been investigated using FD analysis with the goal of introducing an automated diagnostic approach for these specific conditions.

The calculation of the FD of the retinal vascular network comprises the skeletonization of images of the vessels obtained with a fundus camera or a scanning laser ophthalmoscope (SLO). These images are fitted with rectangular grids consisting of boxes of varying sizes.12,1419 The slope of the regression line of the double log plot of the number of boxes versus the box side length is termed D0—the monofractal box-counting FD, which is a global measure that summarizes the whole branching pattern of the retinal vascular tree.14 The box-counting FD has been used to characterize the retinal vascular branching complexity in patients with CI.21,2830 However, the retinal vascular network has a multifractal characteristic rather than a monofractal pattern. Also, it is characterized by a hierarchy of scales rather than a single scale.17,23,25,27 Moreover, the multifractal behavior of the retinal vascular network is characterized by the generalized dimension spectrum (Dq versus q) and the singularity spectrum [f(α) versus α].17,23,25 In our previous study, the capacity or box-counting dimension (D0), the information dimension (D1), and the correlation dimension (D2) were reduced in the cognitively impaired participants compared to the cognitively healthy controls.12

The lacunarity (Λ) parameter is another metric used in the analysis of the retinal vascular network to measure the gap distribution or dispersion of the pixels within the retinal image.25,27,31 It characterizes the heterogeneity of the pixels within an image with a larger Λ, indicating a larger heterogeneity and vice versa.27 Also, it can distinguish between structures of similar FD.27,31 The multi-FD and Λ parameter can both be combined to assess the retinal vascular network of cognitively impaired participants versus age- and sex-matched cognitively healthy controls as initial steps toward an automated diagnostic approach for this condition.

The purpose of this study is twofold. First, we compared the singularity spectrum and Λ parameter of the retinal vascular network of patients with CI to age- and sex-matched cognitively healthy controls. We hypothesized that the singularity spectrum exponents will be reduced in cognitively impaired patients compared to cognitively healthy controls, indicating a less complex retinal branching pattern in the patients with CI. The Λ parameter, on the other hand, was expected to be more heterogeneous or larger in the participants with CI compared to cognitively healthy controls. Second, we investigated the associations between retinal vascular parameters (i.e., singularity spectrum exponents and Λ parameter) and functional parameters, such as the implicit time (IT) and amplitude of the electroretinogram (ERG), and the montreal cognitive assessment (MoCA) scores. These associations were hypothesized to be significant. The outcome of this study may serve as an initial step toward an automated diagnostic approach for the diagnosis of CI at low cost with the added value of an eye-screening.

2.

Materials and Methods

2.1.

Study Participants

Participants with CI were identified in a nonsystematic fashion as they appeared in the clinic or from a group of patients attending adult care centers and community clinics. Exclusion criteria included participants under the age of 55 years, history of any ophthalmic disease except for cataract surgery. The macular and optic disk regions of the fundus images of these subjects were examined for glaucoma, age-related macular degeneration, and other comorbid conditions and participants with such conditions were excluded from the study. Both hypertension and diabetes mellitus, as well as cardiovascular disease, were considered comorbid medical conditions related to retinal vascular alterations. In addition, current or history of study subject-reported smoking was considered in the study. All subjects underwent a cognitive function assessment with the MoCA test, ERG evaluation, as well as fundus imaging. Informed consent was obtained from all participants and for those whose CI affected their decision making on participation in the study, and an informed consent was obtained through an eligible proxy. The study adhered to the tenets of the declaration of Helsinki and was approved by the Human Research Ethics Committee of the University of Miami.

2.2.

Fundus Imaging and Skeletonization of Retinal Vessels

The analysis of the retinal fundus images and skeletonization of the retinal vessels of the subjects in this study have been described in detail in our previous study.12 Briefly, optic disk centered images (FOV of 45 deg) of size 1024×1024  pixels of both eyes were taken with a nonmydriatic SLO digital camera (EasyScan, iOptics, the Netherlands). To ensure the independence of samples, either eye of each subject was chosen based on good image scan quality index. The retinal vessels in the entire 45 deg images were skeletonized, as described in our previous study.12 This method of skeletonization provides more samples of vessels for further analysis unlike using only a select few vessels around the optic disk, as has been done previously.20,21,32

2.3.

ERG and Cognitive Function Assessment

The assessment of the bioelectrical activity of the retina was performed with a full-field ERG (RETeval™, LKC Technologies, Inc., Gaithersburg, Maryland, United States) based on the International Society for Clinical Electrophysiology of Vision (ISCEV) protocol.12,3335 The ISCEV 6 protocol was used, which first includes a light-adapted state, consisting of light-adapted ERG (stimulus strength, 3.0  cd.s/m2; frequency, 28.3 Hz flicker response) and dark-adapted ERG, including rod, maximal dark-adapted, and cone responses. ERG IT and amplitudes were elicited by 141 to 424 flashes separately for each eye.

CI is a condition that individuals face when they have trouble remembering, learning new things, concentrating, or making decisions that affect their everyday life. Therefore, assessment of the cognitive function of study participants was conducted using the MoCA test to identify CI.12 This test is a one-page 30-point assessment performed in 10 to 12 min with a score range of 0 to 30. Typical scores <26 indicate some form of CI.36,37 Detailed assessment of CI using MoCA as well as the MoCA scores of the participants in our current study have been provided previously.12

2.4.

Singularity Spectrum by Multifractal Analysis

The network of blood vessels of the human retina holds a fractal structure with a vascular branching process characterized by self-similarity and scaling. The FD analysis is a computational method that characterizes complexity in a natural phenomenon and has been extensively used to characterize the branching complexity of the human retinal vasculature.12,17,2327,38 The fractal analysis of the human retinal vasculature measures the retinal vascular branching pattern with a complex branching pattern indicating a larger FD and vice versa. This method has been extensively used in the characterization of the retinal vascular network in patients with CI,12,20,21 diabetes,2226 and amblyopia.27

Self-similarity over varying scales or magnification is an important characteristic of fractal structures. The concept is that at different magnification or scales, a fractal structure exhibits a similar pattern at different sizes.12,1419 This characteristic of fractal structures can be described by the equation:12

Eq. (1)

N(r)=constrD,
where N(r) refers to certain measurements applied on the complex structure of the fractal object at varying magnification or scaling r. The parameter D is the FD metric, which refers to the number of self-similar patterns generated from the fractal object as the magnification or scale r increases or decreases.12,1419 The FD calculated with the box-counting method is a common monofractal type of FD that used to characterize the human retinal vascular complexity.2022,26 In this method, the segmented or skeletonized retinal vessels from a retinal fundus or SLO camera are fitted with a rectangular grid consisting of large amount of boxes.12,26 The number of boxes in the rectangular grid in this scenario is N(r) with a box having side length r, as shown in Eq. (1). Once a double log plot of the number of boxes N(r) and the box side length r are obtained, the slope of the regression line in this plot indicates the monofractal FD (box counting-based) metric termed D0, as shown in Eq. (2):12,26

Eq. (2)

D0=limr0logNr/log(1r).

Because the retinal vessel network geometry has a finer texture or heterogeneity in its space-filling characteristics, its morphological property cannot be described sufficiently by a global parameter reflecting a simple FD obtained by monofractal analysis. Therefore, the human retinal vascular network is considered a geometric multifractal structure or a set of intertwined fractals characterized by a hierarchy of exponents rather than a single FD as in D0 (i.e, when the retinal vascular network is considered to have a monofractal pattern). Therefore, the multifractal property of the human retinal vascular network is characterized by the generalized dimension spectrum (Dq versus q, where Dq represents the FD at the q’th order or exponent) and the singularity spectrum [f(α) versus α].17,25,27 The singularity spectrum of the FD f(α) versus the singularity exponent α is defined as17,25

Eq. (3)

N(α)=rf(α),
where N(α) is the number of boxes, such that the probability Pi (r) of finding a pixel within a given region i scales, as observed in Eq. (4):17,25

Eq. (4)

Pi=rαi.

f(α) is hence the FD of all the regions with singularity strengths between α and α+dα, where the singularity exponent α takes on values within the interval and +.

The relationship between the D(q) spectrum and the f(α) spectrum is established through the Legendre transformation, as shown in Eq. (5):1719

Eq. (5)

f(α(q))=qα(q)τ(q),
where α(q) represents the singularity exponent at the q’th order moment expressed1719

Eq. (6)

α(q)=dτ(q)/dq,
and τ(q) represents the mass correlation exponent of the q’th order related to Dq, as observed in Eq. (7):1719

Eq. (7)

τ(q)=(q1)Dq,
where for q=1 and τ(1)  =0.

The singularity FD [f(α(q))] and exponent α(q) of the q’th order can be computed using the software Image J (Wayne Rasband, National Institutes of Health in Bethesda, Maryland, United States) together with the FracLac plugin (A. Karperien—Charles Sturt University, Australia)39 with the settings, as previously described.12 Hence, q represents values from 10 to +10 with an increment of 1. These values are computed and then averaged with their standard deviations at each value of q. Then, a plot of f[α(q)] versus α(q) with standard deviation error bars is obtained, representing the singularity spectrum.17,25 For a multifractal object, the singularity spectrum is typically a parabola with concavity facing down.17,23,25

Typically, the height (Δf), width (Δα), and asymmetry (A) of the parabola are the measures used to describe the singularity spectrum, as shown in Eqs. (8)–(10) below.23,25,40 The α values (α0, α1, and α2) of the singularity spectrum at q=0, 1, 2 can also be compared between groups with greater α values indicating greater singularities or maxima and vice versa:17

Eq. (8)

Δf(α)=f(α)minf(α)max,
where f(α)max and f(α)min represent the maximum and minimum f(α) values, respectively,23,40

Eq. (9)

Δα=αmaxαmin.

Also, αmax and αmin represent the maximum and minimum α values, respectively.25 The higher the Δα, the stronger is the multifractality, and the more complex is the pixel distribution within the image:25,41,42

Eq. (10)

A=(α0αmin)/(αmaxα0).

The parabola of the singularity or F spectrum is symmetric when A=1, left skewed when A>1, and right skewed when A<1.23,25 A left-skewed F spectrum means that there is a stronger presence of high fractal exponents and a significant fluctuation, while a right-skewed F spectrum, on the other hand, posits low fractal exponents and a slight fluctuation.25

2.5.

Lacunarity Analysis

The Λ parameter measures the gap dispersion or the heterogeneity within an object.25,27,31 It can distinguish between two objects with similar fractality and it is a parameter that describes the coarseness or texture of an image.27,31 The Λ parameter of the skeletonized images of the participants was computed with the ImageJ software together with the FracLac plugin with the settings used expatiated, as previously described.12 The ImageJ software together with the FracLac plugin computes the Λ parameter based on the variation in pixel density at different box sizes in fixed and sliding scans. A lower Λ parameter indicates more homogeneity of the pixel distribution in the image and vice versa.27 The mean Λ parameter computed from the Image J software together with the FracLac plugin is then summarized as shown in Eq. (11):25,27,31

Eq. (11)

Λ=1n(j=1gi=1n[1+(σμ)^2]),
where σ is the standard deviation of the number of pixels that were in a box of size ε; μ is the mean for pixels per box at this size ε, in a box count at an orientation g; and then, n is the number of box sizes.25,27

2.6.

Statistical Analysis

All statistical analyses were performed using IBM SPSS Statistics for Windows, Version 24.0 (IBM Corporation, Armonk, New York, USA). All values are presented as mean±SD. The Shapiro–Wilk test of normality was used to test the Gaussian distribution of α0, α1, α2, Δα, A, Δfα, and Λ. Independent sample t-test was used to compare these parameters between the cognitively impaired and cognitively healthy participants assuming Levene’s test for equality of variance was not statistically significant. In cases where the Shapiro–Wilk test was significant for any of the parameters, the Mann–Whitney U-test was used to compare that parameter between the two groups. Cohen’s d was used as an effect size measure when a parameter was significant when compared between the two groups, calculated using GPower calculator 3.1.43 A Cohen’s d value of 0.2 was considered small, 0.5 was considered medium, and 0.8 was considered large.44 Pearson product–moment correlation was used to assess the relationship between the retinal vascular parameters (α0, α1, α2, and Λ) and functional parameters (ERG IT and amplitude, and MoCA) for the cognitively impaired participants with the value 0.7 considered as high. A p value <0.05 was considered statistically significant.

3.

Results

A total of 20 individuals with CI and 19 cognitive healthy participants were involved in this study. The average age of the CI participants (81±6  years) did not significantly differ from that of the cognitively healthy participants (80±7  years), t(37)=0.38, p=0.70. The proportion of males (20%) and females (80%) in the participants with CI did not significantly differ from the proportion of males (16%) and females (84%) in the cognitive healthy participants (Fisher’s exact test, p=1.00). The singularity spectrum of both the cognitively impaired and cognitive healthy participants followed a parabola with concavity facing down demonstrating the multifractal property of the retinal vascular network of the two groups (Fig. 1). The singularity spectrum of the participants with CI also showed a trend or a shift toward lower singularity exponents or α range and had lower maxima compared to the singularity spectrum of the cognitively healthy participants even though there was quite an overlap for Δf and Δα between the two groups [Fig. 1(c)].

Fig. 1

Singularity or F spectrum of cognitively impaired and cognitively healthy participants showing a parabola with concavity facing down demonstrating the multifractality of the retinal vessels of both groups. (a) The singularity spectrum of the cognitively healthy participants. The singularity spectrum shows a multifractal behavior, as demonstrated by the parabola with concavity facing down. (b) Singularity spectrum of the participants with CI showing a multifractal behavior, as demonstrated by the parabola with concavity facing down. (c) Singularity spectrum of both groups showing a trend or shift of the singularity spectrum in participants with CI toward a lower α range or maxima but quite an overlap of Δf and Δα between the two groups. Error bars represent standard deviation.

NPH_6_4_041109_f001.png

The values of α0, α1, and α2 were significantly greater in the cognitively healthy participants than in the cognitively impaired participants, with greater than moderate Cohen’s d or effect size for α0 and α1 (d=0.63, 0.63, respectively) and a large Cohen’s d for α2 (d=0.97) (Table 1). Thus, the participants with CI had lower singularity or α values compared to the cognitively healthy participants and this difference had a high effect size for α2 but not α0 and α1. The values of Δf, Δα, and A of the F spectrum did not significantly differ between the cognitively healthy participants and the participants with CI (Table 1). The singularity spectrum of both the cognitively impaired and cognitively healthy participants was right-skewed (A<1), indicating the presence of low fractal exponents and a slight fluctuation. Moreover, the Λ parameter was not significantly different between the participants with CI and the cognitively healthy participants (Table 1).

Table 1

Multifractal and Lacunarity parameters (mean±SD) obtained for the cognitively impaired and cognitive healthy participants.

Multifractal and lacunarity parametersCognitively impaired groupCognitively healthy groupp valueCohen’s d
α01.60±0.061.63±0.030.030.63
α11.56±0.061.59±0.030.030.63
α21.54±0.051.58±0.030.020.97
Δf0.87±0.100.87±0.070.93NA
Δα0.74±0.050.71±0.050.14NA
Aa0.35±0.060.35±0.050.99NA
Λ0.35±0.050.34±0.030.48NA
Note: α0, α1, and α2 represent the singularity exponents at q=0, 1, 2, respectively. The Δf, Δα, and A represent the height, width, and asymmetry of the singularity spectrum, respectively. The Λ parameter represents lacunarity, indicating the gap dispersion within the image. NA: not applicable.

aMann–Whitney U-test was performed, otherwise an independent sample t-test was performed.

There was a significant positive association between retinal vascular singularity exponents (α0, α1, and α2) and the IT, which was one of the functional measures obtained with the ERG assessments (Table 2). In particular, there was a high Pearson correlation coefficient between α2 versus ERG IT. However, a negative association was observed between the Λ parameter versus the ERG IT (Table 2). There were no such associations between α0, α1, and α2, and Λ versus the ERG amplitude and MoCA, p>0.05, Table 2.

Table 2

Association between retinal vascular measures (i.e., multifractal and lacunarity parameters) and functional measures (ERG IT, ERG amplitude, and MoCA).

ParametersPearson’s correlation (r)p value
α0 versus ERG IT0.610.004
α1 versus ERG IT0.670.001
α2 versus ERG IT0.71<0.001
Λ versus ERG IT0.510.022
α0 versus ERG amplitude0.410.07
α1 versus ERG amplitude0.350.13
α2 versus ERG amplitude0.290.22
Λ versus ERG amplitude0.230.33
α0 versus MoCA0.430.06
α1 versus MoCA0.430.06
α2 versus MoCA0.420.07
Λ versus MoCA0.180.44

4.

Discussion

There is an emerging interest in the identification of retinal vascular and neural biomarkers for the early diagnosis of CI, especially due to AD, during the prodromal phase of the disease because the neurodegeneration in the brain of AD patients occurs years before clinical symptoms appear.38,12,21,4550 Although there is active neurodegeneration occurring in the brain, patients in the prodromal phase of the disease typically can go about their daily activities without showing any signs of CI. Therefore, the target of treatment and management to slow down disease progression is to diagnose the disease during the prodromal phase.51 However, AD and other related dementia diseases are diagnosed through expensive and invasive methods, such as PET and cerebrospinal fluid assessment via a spinal tap.2 Therefore, there is a desire for a comparatively lower cost and noninvasive method that may aid in the early diagnosis of the disease. As the human retina is an extension of the brain and retinal imaging can be attained noninvasively at low cost, the human retina offers a noninvasive window to identify retinal vascular and neural biomarkers of brain diseases.

In this study, we identified changes in retinal vascular parameters (FD and Λ) in patients with CI that were significantly associated with a particular functional measure obtained with the ERG assessment (IT). Also, changes in some fractal exponents and their association with the bioelectrical activity of the retina (ERG IT) were found to have a large effect size. These associations interestingly pointed to a clear perturbation of the neurovascular component as a result of abnormal conditions mediated by the individual’s disease status affecting both the brain and eye structures. We also found significantly reduced singularity spectrum exponents (α0, α1, and α2) in the participants with CI compared to the cognitively healthy participants with a large effect size difference for α2 (Table 1), which may indicate that the reduced retinal vascular branching complexity for patients with CI, possibly due to reduced retinal neural demand, could be incorporated as a clinical tool for the diagnosis of CI. The multifractal behavior of the retinal vasculature of both groups is consistent with that found in previous studies.17,23,25 However, there was quite an overlap or no significant difference between the height (Δf), width (Δα), and asymmetry (A) of the singularity spectrum of both groups (Fig. 1 and Table 1). As in a previous study comparing patients with diabetic retinopathy and healthy controls, the lack of a significant difference in the width and asymmetry of the singularity spectrum as well as in the gap dispersion of the pixels (Λ) could be attributed to the small sample size.25

Previous studies have investigated the fractal properties of the retinal vasculature of patients with CI.12,20,21 However, most studies that investigated the FD of cognitively impaired patients used the monofractal box-counting dimension approach rather than a multifractal spectrum methodology.20,21 Also, these studies have limited their analyses by sampling few vessels around the optic disk rather than using the whole skeletonized vessel tree in the entire fundus image of the participants.20,21 We have previously characterized the multifractal behavior of the retinal vasculature of participants with CI versus cognitively healthy participants using only the generalized dimension spectrum.12 The results obtained for α0, α1, and α2 demonstrated that these singularity exponents were significantly reduced in the cognitively impaired individuals versus controls, which is consistent with the results of our previous study, where D0, D1, and D2 were also significantly reduced in the cognitively imapired group.12 The significantly reduced FD found in the group with CI is also consistent with that found in other studies that investigated the retinal vascular branching pattern of cognitively impaired participants although a monofractal box-counting approach was used in these studies.20,21

The neurovascular coupling concept of the retina calls for increased vascular demand with associated increased neural activity.52,53 It is well established that the neurodegeneration that occurs in the brain of patients with CI is associated with loss of retinal ganglion cells.10,11 This trend implies that there will be reduced vascular demand in terms of nutrients and oxygen from the retinal neurons in these subjects and a potential remodeling of the branching pattern complexity of the retinal vessels in these subjects. This tendency may explain the significantly reduced singularity spectrum exponents, as well as a trend toward lower maxima in the subjects with CI versus the age- and sex-matched cognitively healthy controls. Interestingly, this trend toward lower singularity spectrum exponents or maxima for the cognitively impaired participants is consistent with that found in another study that investigated the singularity spectrum for various retinal pathological conditions.17 The significant difference in the singularity exponents between the two groups had a high effect size for α2 despite the small sample size in our study. This particular result implies that the retinal vessel branching pattern of the participants with CI was less complex compared with the cognitively healthy participants and this difference had a large effect size beyond statistical significance. Therefore, the singularity exponents of the F spectrum multifractal dimension could be incorporated into an automated diagnostic approach for participants with CI using the information of the retinal vascular branching complexity.

The significant associations between retinal vascular parameters and the ERG IT followed a positive trend for the retinal vascular FD parameters (α0, α1, and α2) but a negative trend for Λ (Table 2). Specifically, the negative association between the Λ parameter and the ERG IT is expected as FD and Λ are hypothesized to be negatively associated.54 The significant positive associations between the retinal vascular FD parameters (α0, α1, and α2) and ERG IT had a large effect size or Pearson correlation coefficient for α2 (Table 2), which is consistent with the concept of neurovascular coupling.52,53 These associations are also consistent with the results of our previous study, which found significant positive associations between the generalized dimension spectrum and ERG IT.12

The results of this study should be interpreted in the light of some limitations. First, the sample size of our study was small, and this may have accounted for the overlap or lack of significant difference in the Λ parameter, height, width, and asymmetry of singularity spectrum between the two groups as well as the lack of significant associations between the retinal vascular FD and the Λ parameter versus the ERG amplitude and the MoCA scores for the cognitively impaired group. However, the sample size did not affect the significant difference and large effect size found for the retinal vascular FD parameters between the two groups as well as the associations between the retinal vascular and functional parameters. Second, the cause of CI for our patients is unknown—this could be due to AD, vascular dementia, frontotemporal lesion, and other dementia-related diseases. The causation inquiry may require a more comprehensive approach, including the use of PET, magnetic resonance imaging, and cerebrospinal fluid analysis.

The process of finding biomarkers is complicated, and despite the small sample size of our study’s population, we found significant associations between retinal vascular and functional parameters with large effect sizes that can be incorporated into longitudinal studies with a larger sample size for an automated diagnostic approach that could aid the detection of CI using retinal biomarkers. These multimodal parameters could be implemented in an automated screening setting using low-cost fundus imaging devices that could even be coupled with low-cost electrophysiological screening tools. All this together, used in an appropriate setting and targeting populations at high risk could, in turn, help not only obtain an early diagnosis, and thus a reduced burden on society, but also serve in the better understanding of the development of CI.

Disclosures

The University of Miami and Delia Cabrera DeBuc have filed a provisional patent application 62/776,693. The other authors report no conflict of interest and have no proprietary interest in any of the materials published in this article.

Acknowledgments

This study was supported in part by the Finker Frenkel Legacy Foundation, U.S. National Institutes of Health Center Grant No. P30-EY014801 to the University of Miami, an unrestricted grant to the University of Miami from Research to Prevent Blindness, Inc., and the Alzheimer’s Association (Grant No. AARGD-17-531255, DCD).

References

1. 

M. J. Prince, “World Alzheimer Report 2015: the global impact of dementia: an analysis of prevalence, incidence, cost and trends,” (2019). http://www.alz.co.uk/research/WorldAlzheimerReport2011.pdf Google Scholar

2. 

E. M. Lad et al., “Evaluation of inner retinal layers as biomarkers in mild cognitive impairment to moderate Alzheimer’s disease,” PLoS One, 13 (2), e0192646 (2018). https://doi.org/10.1371/journal.pone.0192646 POLNCL 1932-6203 Google Scholar

3. 

A. La Rue and L. F. Jarvik, “Cognitive function and prediction of dementia in old age,” Int. J. Aging Hum. Dev., 25 (2), 79 –89 (1987). https://doi.org/10.2190/DV3R-PBJQ-E0FT-7W2B Google Scholar

4. 

R. T. Linn et al., “The ‘preclinical phase’ of probable Alzheimer’s disease. A 13-year prospective study of the Framingham cohort,” Arch. Neurol., 52 (5), 485 –490 (1995). https://doi.org/10.1001/archneur.1995.00540290075020 Google Scholar

5. 

D. A. Snowdon et al., “Linguistic ability in early life and cognitive function and Alzheimer’s disease in late life. Findings from the Nun Study,” JAMA, 275 (7), 528 –532 (1996). https://doi.org/10.1001/jama.1996.03530310034029 JAMAAP 0098-7484 Google Scholar

6. 

H. Braak and E. Braak, “Frequency of stages of Alzheimer-related lesions in different age categories,” Neurobiol. Aging, 18 (4), 351 –357 (1997). https://doi.org/10.1016/S0197-4580(97)00056-0 NEAGDO 0197-4580 Google Scholar

7. 

M. F. Elias, “The preclinical phase of alzheimer disease: a 22-year prospective study of the Framingham cohort,” Arch. Neurol., 57 (6), 808 –813 (2000). https://doi.org/10.1001/archneur.57.6.808 Google Scholar

8. 

C. H. Kawas et al., “Visual memory predicts Alzheimer’s disease more than a decade before diagnosis,” Neurology, 60 (7), 1089 –1093 (2003). https://doi.org/10.1212/01.WNL.0000055813.36504.BF NEURAI 0028-3878 Google Scholar

9. 

D. C. DeBuc, G. M. Somfai and A. Koller, “Retinal microvascular network alterations: potential biomarkers of cerebrovascular and neural diseases,” Am. J. Physiol. Heart. Circ. Physiol., 312 (2), H201 –H212 (2017). https://doi.org/10.1152/ajpheart.00201.2016 Google Scholar

10. 

B. Katz and S. Rimmer, “Ophthalmologic manifestations of Alzheimer’s disease,” Surv. Ophthalmol., 34 (1), 31 –43 (1989). https://doi.org/10.1016/0039-6257(89)90127-6 SUOPAD 0039-6257 Google Scholar

11. 

H. Hampel et al., “Revolution of Alzheimer precision neurology. Passageway of systems biology and neurophysiology,” J. Alzheimer Dis., 64 (s1), S47 –S105 (2018). https://doi.org/10.3233/JAD-179932 Google Scholar

12. 

D. C. DeBuc et al., “Investigating multimodal diagnostic eye biomarkers of cognitive impairment by measuring vascular and neurogenic changes in the retina,” Front. Physiol., 9 1721 (2018). https://doi.org/10.3389/fphys.2018.01721 FROPBK 0301-536X Google Scholar

13. 

M. M. Moschos et al., “Structural and functional impairment of the retina and optic nerve in Alzheimer’s disease,” Curr. Alzheimer Res., 9 (7), 782 –788 (2012). https://doi.org/10.2174/156720512802455340 Google Scholar

14. 

G. Liew et al., “The retinal vasculature as a fractal: methodology, reliability, and relationship to blood pressure,” Ophthalmology, 115 (11), 1951 –1956.e1 (2008). https://doi.org/10.1016/j.ophtha.2008.05.029 OPANEW 0743-751X Google Scholar

15. 

V. F. Cosatto et al., “Retinal vascular fractal dimension measurement and its influence from imaging variation: results of two segmentation methods,” Curr. Eye Res., 35 (9), 850 –856 (2010). https://doi.org/10.3109/02713683.2010.490628 CEYRDM 0271-3683 Google Scholar

16. 

G. N. Thomas et al., “Measurement of macular fractal dimension using a computer-assisted program,” Invest. Ophthalmol. Visual Sci., 55 2237 –2243 (2014). https://doi.org/10.1167/iovs.13-13315 IOVSDA 0146-0404 Google Scholar

17. 

T. Stosic and B. D. Stosic, “Multifractal analysis of human retinal vessels,” IEEE Trans. Med. Imaging, 25 (8), 1101 –1107 (2006). https://doi.org/10.1109/TMI.2006.879316 ITMID4 0278-0062 Google Scholar

18. 

S. Talu, “Characterization of retinal vessel networks in human retinal imagery using quantitative descriptors,” Hum. Vet. Med., 5 (2), 52 –57 (2013). Google Scholar

19. 

S. Talu, “Multifractal geometry in analysis and processing of digital retinal photographs for early diagnosis of human diabetic macular edema,” Curr. Eye Res., 38 (7), 781 –792 (2013). https://doi.org/10.3109/02713683.2013.779722 CEYRDM 0271-3683 Google Scholar

20. 

C. Y. Cheung et al., “Retinal vascular fractal dimension is associated with cognitive dysfunction,” J. Stroke Cerebrovasc. Dis., 23 (1), 43 –50 (2014). https://doi.org/10.1016/j.jstrokecerebrovasdis.2012.09.002 Google Scholar

21. 

C. Y. Cheung et al., “Microvascular network alterations in the retina of patients with Alzheimer’s disease,” Alzheimer Dement., 10 (2), 135 –142 (2014). https://doi.org/10.1016/j.jalz.2013.06.009 Google Scholar

22. 

M. Kostic et al., “Investigating the fractal dimension of the foveal microvasculature in relation to the morphology of the foveal avascular zone and to the macular circulation in patients with Type 2 diabetes mellitus,” Front. Physiol., 9 1233 (2018). https://doi.org/10.3389/fphys.2018.01233 FROPBK 0301-536X Google Scholar

23. 

T. Zhu et al., “Multifractal and lacunarity analyses of microvascular morphology in eyes with diabetic retinopathy: a projection artifact resolved optical coherence tomography angiography study,” Microcirculation, 26 (3), e12519 (2019). https://doi.org/10.1111/micc.2019.26.issue-3 Google Scholar

24. 

A. Avakian et al., “Fractal analysis of region-based vascular change in the normal and non-proliferative diabetic retina,” Curr. Eye Res., 24 (4), 274 –280 (2002). https://doi.org/10.1076/ceyr.24.4.274.8411 CEYRDM 0271-3683 Google Scholar

25. 

E. V. Costa and R. A. Nogueira, “Fractal, multifractal and lacunarity analysis applied in retinal regions of diabetic patients with and without nonproliferative diabetic retinopathy,” Fractal Geom. Nonlinear Anal. Med. Biol., 1 (3), 112 –119 (2015). https://doi.org/10.15761/FGNAMB.1000118 Google Scholar

26. 

N. Popovic et al., “Fractal dimension and lacunarity analysis of retinal microvascular morphology in hypertension and diabetes,” Microvasc. Res., 118 36 –43 (2018). https://doi.org/10.1016/j.mvr.2018.02.006 MIVRA6 0026-2862 Google Scholar

27. 

S. Tălu, C. Vlăduţiu and C. A. Lupaşcu, “Characterization of human retinal vessel arborisation in normal and amblyopic eyes using multifractal analysis,” Int. J. Ophthalmol., 8 (5), 996 –1002 (2015). https://doi.org/10.3980/j.issn.2222-3959.2015.05.26 Google Scholar

28. 

S. Frost et al., “Retinal vascular biomarkers for early detection and monitoring of Alzheimer’s disease,” Transl. Psychiatry, 3 e233 (2013). https://doi.org/10.1038/tp.2012.150 Google Scholar

29. 

Y. T. Ong et al., “Retinal vascular fractals and cognitive impairment,” Dement. Geriatr. Cognit. Disord., 4 (2), 305 –313 (2014). https://doi.org/10.1159/000363286 Google Scholar

30. 

M. A. Williams et al., “Retinal microvascular network attenuation in Alzheimer’s disease,” Alzheimers Dement., 1 (2), 229 –235 (2015). https://doi.org/10.1016/j.dadm.2015.04.001 Google Scholar

31. 

Ş. Ţălu et al., “Fractal and lacunarity analysis of human retinal vessel arborisation in normal and amblyopic eyes,” Hum. Vet. Med., 5 (2), 45 –51 (2013). Google Scholar

32. 

C. Y. Cheung et al., “Quantitative and qualitative retinal microvascular characteristics and blood pressure,” J. Hypertens., 29 (7), 1380 –1391 (2011). https://doi.org/10.1097/HJH.0b013e328347266c JOHYD3 0263-6352 Google Scholar

33. 

M. F. Marmor et al., “Standard for clinical electroretinography (2004 update),” Doc. Ophthalmol., 108 (2), 107 –114 (2004). https://doi.org/10.1023/B:DOOP.0000036793.44912.45 Google Scholar

34. 

G. E. Holder et al., “ISCEV standard for clinical pattern electroretinography—2007 update,” Doc. Ophthalmol., 114 (3), 111 –116 (2007). https://doi.org/10.1007/s10633-007-9053-1 Google Scholar

35. 

D. C. Hood et al., “ISCEV guidelines for clinical multifocal electroretinography (2007 edition),” Doc. Ophthalmol., 116 (1), 1 –11 (2008). https://doi.org/10.1007/s10633-007-9089-2 Google Scholar

36. 

T. Smith, N. Gildeh and C. Holmes, “The Montreal cognitive assessment: validity and utility in a memory clinic setting,” Can. J. Psychiatry, 52 (5), 329 –332 (2007). https://doi.org/10.1177/070674370705200508 CPAJAK 0706-7437 Google Scholar

37. 

Z. S. Nasreddine et al., “The Montreal cognitive assessment, MoCA: a brief screening tool for mild cognitive impairment,” J. Am. Geriatr. Soc., 53 (4), 695 –699 (2005). https://doi.org/10.1111/jgs.2005.53.issue-4 JAGSAF 0002-8614 Google Scholar

38. 

Ş. Ţălu and S. Giovanzana, “Fractal and multifractal analysis of human retinal vascular network: a review,” Hum. Vet. Med., 3 (3), 205 –212 (2011). Google Scholar

39. 

A. Karperien, “FracLac for ImageJ—FracLac advanced user’s manual,” (1993) http://rsb.info.nih.gov/ij/plugins/fraclac/fraclac-manual.pdf April 2019). Google Scholar

40. 

L. G. França et al., “Fractal and multifractal properties of electrographic recordings of human brain activity: toward its use as a signal feature for machine learning in clinical applications,” Front. Physiol., 9 1767 (2018). https://doi.org/10.3389/fphys.2018.01767 FROPBK 0301-536X Google Scholar

41. 

K. Shi, C. Q. Liu and N. S. Ai, “Monofractal and multifractal approaches in investigating temporal variation of air pollution indexes,” Fractals, 17 (4), 513 –521 (2009). https://doi.org/10.1142/S0218348X09004454 FRACEG 0218-348X Google Scholar

42. 

M. G. Hu, J. F. Wang and Y. Ge, “Super-resolution reconstruction of remote sensing images using multifractal analysis,” Sensors, 9 (11), 8669 –8683 (2009). https://doi.org/10.3390/s91108669 SNSRES 0746-9462 Google Scholar

43. 

F. Faul et al., “Statistical power analyses using G∗Power 3.1: tests for correlation and regression analyses,” Behav. Res. Methods, 41 (4), 1149 –1160 (2009). https://doi.org/10.3758/BRM.41.4.1149 Google Scholar

44. 

J. Cohen, “A power primer,” Psychol. Bull., 112 (1), 155 –159 (1992). https://doi.org/10.1037/0033-2909.112.1.155 PSBUAI 0033-2909 Google Scholar

45. 

B. E. O’Bryhim et al., “Association of preclinical Alzheimer disease with optical coherence tomographic angiography findings,” JAMA Ophthalmol., 136 (11), 1242 –1248 (2018). https://doi.org/10.1001/jamaophthalmol.2018.3556 Google Scholar

46. 

N. J. Hart et al., “Ocular indicators of Alzheimer’s: exploring disease in the retina,” Acta Neuropathol., 132 (6), 767 –787 (2016). https://doi.org/10.1007/s00401-016-1613-6 ANPTAL 1432-0533 Google Scholar

47. 

Y. Koronyo et al., “Retinal amyloid pathology and proof-of-concept imaging trial in Alzheimer’s disease,” JCI Insight, 2 (16), e93621 (2017). https://doi.org/10.1172/jci.insight.93621 Google Scholar

48. 

D. Liu et al., “Thinner changes of the retinal nerve fiber layer in patients with mild cognitive impairment and Alzheimer’s disease,” BMC Neurol., 15 (1), 14 (2015). https://doi.org/10.1186/s12883-015-0268-6 Google Scholar

49. 

C. Y. Cheung et al., “Retinal ganglion cell analysis using high-definition optical coherence tomography in patients with mild cognitive impairment and Alzheimer’s disease,” J Alzheimer Dis., 45 (1), 45 –56 (2015). https://doi.org/10.3233/JAD-141659 Google Scholar

50. 

M. Koronyo-Hamaoui et al., “Identification of amyloid plaques in retinas from Alzheimer’s patients and noninvasive in vivo optical imaging of retinal plaques in a mouse model,” Neuroimage, 54 S204 –S217 (2011). https://doi.org/10.1016/j.neuroimage.2010.06.020 NEIMEF 1053-8119 Google Scholar

51. 

M. Prince, R. Bryce and C. Ferri, “World Alzheimer report 2011: the benefits of early diagnosis and intervention,” (2011). Google Scholar

52. 

H. Radhakrishnan and V. J. Srinivasan, “Multiparametric optical coherence tomography imaging of the inner retinal hemodynamic response to visual stimulation,” J. Biomed. Opt., 18 (8), 086010 (2013). https://doi.org/10.1117/1.JBO.18.8.086010 JBOPFO 1083-3668 Google Scholar

53. 

M. Hammer et al., “Retinal venous oxygen saturation increases by flicker light stimulation,” Invest. Ophthalmol. Visual Sci., 52 (1), 274 –277 (2011). https://doi.org/10.1167/iovs.10-5537 IOVSDA 0146-0404 Google Scholar

54. 

Jr. T. G. Smith, G. D. Lange and W. B. Marks, “Fractal methods and results in cellular morphology—dimensions, lacunarity and multifractals,” J. Neurosci. Methods, 69 (2), 123 –136 (1996). https://doi.org/10.1016/S0165-0270(96)00080-5 JNMEDT 0165-0270 Google Scholar

Biography

Edmund Arthur received his doctor of optometry (OD) degree from Kwame Nkrumah University of Science and Technology-Ghana (2014) and his PhD from Indiana University Bloomington (2018). His PhD dissertation focused on detecting preclinical retinal changes in diabetic patients using advanced retinal imaging. He is currently a postdoctoral research associate at the Quantitative Imaging and Eye-Brain Research Lab at Bascom Palmer Eye Institute, University of Miami, researching retinal biomarkers for early diagnosis of Alzheimer’s disease.

Gabor Mark Somfai received his MD degree (2001) and his PhD degree in ophthalmology (2007) at Semmelweis University in Budapest, Hungary. As an ophthalmologist, his main clinical interest is the medical and surgical therapy of retinal diseases, mostly diabetic retinopathy, while as a researcher, he is focusing on retinal imaging analysis for the understanding of retinal pathophysiology. He is currently working as a consultant at the Retinology Unit of the Pallas Kliniken, in Olten, Switzerland.

Maja Kostic received her doctor of medicine (MD) degree, masters, and PhD from the University of Novi Sad, Serbia. Her PhD dissertation focused on morphometric fractal analysis of retinal blood vessels development in prematurely born babies. Currently, she is a pediatric ophthalmology fellow at the Bascom Palmer Eye Institute, University of Miami. Her research mainly focuses on cognitive functions and retinal biomarkers for the early diagnosis of Alzheimer’s disease.

Susel Oropesa received her MD (2018) from Columbia University College of Physicians and Surgeons and her BS degree in molecular biophysics and biochemistry (2014) from Yale University. She is currently an ophthalmology resident at Mount Sinai Hospital. She has broad interests in the field of ophthalmology including retinal and oculoplastic diseases. Her most recent project focuses on discovering patterns in retinal vascular changes that are associated with cognitive impairment.

Carlos Mendoza-Santiesteban received his MD (1994) from Havana Medical University, Cuba. He is currently an assistant professor of clinical ophthalmology at Bascom Palmer Eye Institute, University of Miami. He has more than 20 years of experience in clinical neuro-ophthalmology and medical retina with a special interest in diagnostic tools, including visual electrophysiology and retinal imaging. Also, he has extensive experience in conducting comprehensive evaluations in patients with optic neuropathies, familial dysautonomia, and autosomal dominant optic atrophy.

Delia Cabrera DeBuc received her PhD in applied physics (2002) from the University of Michigan. She is currently a research associate professor leading the Quantitative Imaging and Eye-Brain Research Laboratory at the Bascom Palmer Eye Institute, University of Miami. She is a biophysicist with specific training in imaging processing and data science. She is also leading research on ocular imaging biomarkers in CNS neurodegenerative diseases, such as multiple sclerosis, Parkinson’s disease, and Alzheimer’s disease.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Edmund Arthur, Gabor Mark Somfai, Maja Kostic, Susel Oropesa, Carlos Mendoza-Santiesteban, and Delia Cabrera DeBuc "Distinguishing cognitive impairment by using singularity spectrum and lacunarity analysis of the retinal vascular network," Neurophotonics 6(4), 041109 (23 September 2019). https://doi.org/10.1117/1.NPh.6.4.041109
Received: 24 May 2019; Accepted: 27 August 2019; Published: 23 September 2019
Lens.org Logo
CITATIONS
Cited by 6 scholarly publications.
Advertisement
Advertisement
KEYWORDS
Cognitive informatics

Information technology

Brain

Retina

Fractal analysis

Control systems

Diagnostics

Back to Top