Optical coherence tomography angiography in preclinical Alzheimer's disease.

BACKGROUND/AIMS: As a protrusion from the brain, the retina might reflect the status of the brain. Previous studies showed a decrease in vessel density and foveal avascular zone (FAZ) enlargement on optical coherence tomography angiography (OCTA) in individuals suffering from Alzheimer's disease (AD). This study aims to assess whether such changes are already present in preclinical stages of AD, in a population of monozygotic (MZ) twins.
METHODS: 124 cognitively healthy individuals (MZ twins, ages 60-93 years) underwent [18F]flutemetamol amyloid positron emission tomography (PET) scanning and OCTA. PET scans were visually rated for cortical amyloid-beta (Aβ) positivity. Parametric global cortical non-displaceable binding potential (BP ND) was used as a continuous measure for Aβ aggregation. FAZ size and vessel densities for the inner and outer ring of the macular ETDRS grid and in a 3-6 mm ring around the optic nerve head (ONH) were measured.OCTA measures were associated with visual Aβ score, BP ND and amyloid load estimated by twin concordance on visual Aβ score. Twin correlations were estimated as a measure of maximum heritability of OCTA measures.
RESULTS: 13 of 124 participants were Aβ+. Aβ+ individuals had significantly higher vessel density than Aβ- individuals in all regions but did not differ in FAZ size. Twin analyses showed a positive association between and vessel densities in all regions. BP ND tended to be associated with higher vessel density in the inner ring. Twin correlations were moderate/high for all OCTA parameters except vessel density around the ONH, which correlated weakly.
CONCLUSION: Retinal vessel density was higher in individuals with preclinical AD. © Author(s) (or their employer(s)) 2020. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

Introduction

Dementia, and in particular Alzheimer’s disease (AD), is the most common neurological disorder worldwide, and it is estimated that 1 in 3 of those born in developed countries today will develop dementia during their life.1 AD is a disease characterised by the accumulation of amyloid-beta (Aβ) and tau in the brain.2 At present, there is no curative treatment for AD, despite many trials targeting, for example, Aβ deposits in the brain.3 4 The most common explanation for the poor results is one of timing: once severe neurodegeneration has occurred, there is no way back. As most of the trials so far have focused on patients already suffering from AD, it is postulated that more success may be achieved by targeting those who still are cognitively intact.3 4

This, however, poses a new problem. How can we recognise those at risk for the development of AD when there are no clinical symptoms yet? PET scanning has provided a breakthrough in diagnosing these cases of ‘preclinical’ AD. The process of Aβ accumulation in the brain is a gradual one, which often has been ongoing for decades prior to onset of clinical symptoms.5 By using tracers sensitive to Aβ, PET enables visualisation of Aβ presence in vivo in cognitively healthy individuals.6 Presence of Aβ in the brain of healthy individuals constitutes a high risk of eventual developing AD, which is why these people are often labelled as suffering from preclinical AD.6 7

However, as reliable as this technique may be, it is currently not suitable for large-scale screening. It is a costly diagnostic procedure that is only available in larger hospitals.8 Additionally, PET scanning makes use of intravenously delivered radioactive tracers, which is unpleasant to the patient and subjects them to radioactivity.9

This illustrates the urgent need for an easy, non-invasive and reliable biomarker for preclinical AD. The eye, and especially the retina, may provide an answer to this need. The retina is closely related to the brain, as is reflected by the shared embryological origin and similarities such as a blood–retina barrier and glial cell connections.10 Many studies have already illustrated changes in the retina of individuals suffering from AD, such as retinal thinning and vascular changes.11–13

Conveniently, the retina is easily accessible through the optically clear media of the eye by using optical techniques. Optical coherence tomography (OCT) is such a technique, allowing imaging of the retina with a micrometre resolution.14 The recently developed functional extension of OCT termed OCT angiography (OCTA) enables a detailed angiographic view of the retinal vascular network. As AD is a disease that is also characterised by vascular changes, OCTA may reflect such vascular changes in the retina, making it of interest as a potential new and possibly early biomarker for AD.15 Some groups have already illustrated that vessel density is diminished and that the foveal avascular zone (FAZ) is enlarged in individuals with AD type dementia.16–18 Very few groups, however, have looked whether similar findings are present in preclinical stages. If this is the case, OCTA may prove to be the invaluable biomarker that is sorely needed. The purpose of this study was to elucidate the diagnostic capabilities of OCTA in recognising individuals with preclinical AD, as confirmed by the presence of cerebral Aβ on PET scanning. Furthermore, by including monozygotic (MZ) twins, secondary aims were to: (1) elucidate the contribution of genes to vessel density and FAZ size on OCTA and (2) gain more insight in the aetiology of possible relationships between cerebral Aβ and retinal vascular changes on OCTA.

Methods

Participants

This study is part of the Amsterdam substudy of the European Medical Information Framework for Alzheimer’s Disease (EMIF-AD) PreclinAD cohort. From this cohort with 194 participants from the Netherlands Twin Registry,19 20 150 participants received a complete ophthalmological examination including OCTA. Written informed consent was obtained from all participants, and the study followed the Tenets of the Declaration of Helsinki.

Inclusion criteria for the EMIF-AD study were: age ≥60 years, monozygosity, cognitively healthy as defined by: Telephone Interview for Cognitive Status modified score >22, Geriatric Depression Scale score <11, Consortium to Establish a Registry for Alzheimer’s Disease 10-item word list immediate and delayed recall >−1.5 SD of age adjusted normative data and Clinical Dementia Rating (CDR) scale of 0 with a score on the memory sub domain of 0.

Exclusion criteria for the EMIF-AD study were: uncontrolled diabetes mellitus, alcohol consumption >35 units (1 unit=10 mL or 8 g of pure alcohol) per week, stroke with physical impairment, neurodegenerative disorders and/or cancer with a terminal life expectancy.20

Ophthalmological examination

The following ophthalmological data were obtained: best-corrected visual acuity, intraocular pressure, refraction data, slit lamp examination, indirect funduscopy, fundus photography and conventional OCT of the macula and optic disc. Pupil dilation was achieved using tropicamide 0.5%. OCT images and photographs were assessed by an experienced ophthalmologist (H-TN) for unexpected pathology. Participants suffering from ophthalmological conditions interfering with the retina or image quality were excluded prior to analyses (cataract, macular degeneration, glaucoma, diabetic retinopathy, vascular occlusions and refractive error <-8 dioptres or >+5 dioptres).

Optical coherence tomography angiography

Standard 6×6 mm angiographic scans consisting of 350 B-scans of the macular region and around the optic nerve head (ONH) were obtained using the Cirrus 5000 Angioplex (Zeiss Meditec, Inc, Germany). Vessel densities were measured in all scans, as well as the area of the FAZ in the macular scans. For macular scans, vessel densities in the ETDRS grid were averaged to obtain one value for the inner ring (1–3 mm around the fovea) and the outer ring (3–6 mm around the fovea), while for the ONH scans, vessel density was measured in a 3–6 mm ring around the ONH. Values were averaged over left and right eyes. If only one eye could be scanned reliably, only the values of that eye were used.

OCTA scans were acquired on average 22 months after the amyloid PET scans. This is due to OCTA being available in our hospital only recently, at the time of the 2-year follow-up.

Amyloid PET scanning

All PET scans were acquired using a TF PET-MRI scanner (Philips Medical Systems, Best, the Netherlands). After intravenous injection of 185 MBq (±10%) [18F]flutemetamol, all participants were scanned dynamically from 0 min to 30 min and then again from 90 min to 110 min.21 A T1-weighted gradient echo pulse MRI scan was obtained immediately prior to each part of the PET scan, which was used for attenuation correction. The first dynamic emission scan was reconstructed into 18 frames with progressive increase in frame length (6×5, 3×10, 4×60, 2×150, 2×300 and 1×600 s), the second part of the scan consisted of 4×5 min frames. Vinci viewing software 2.56 (Max Planck Institute for neurological research, Cologne, Germany) and in-house build tools were used to coregister and combine the two PET scans into a single multiframe image. Each individual’s T1 was coregistered to the dynamic PET image using the generic multimodality setting of Vinci with a linear rigid body scheme and normalised mutual information as the similarity measure. The entire image set was used to generate parametric binding potential relative to the non-displaceable compartment (non-displaceable binding potential or BP ND) images using the receptor parametric mapping22 23 implementation in PPET.24 Standard uptake value ratio (SUVr) images were constructed based on the 90–110 min acquired data. Cerebellar grey matter, defined by the Hammers atlas, was used as reference region for these analyses.25 Global values were computed based on the average of frontal (volume-weighted average of superior, middle and inferior frontal gyrus), parietal (volume weighted average of posterior cingulate, superior parietal gyrus, postcentral gyrus and inferolateral remainder of parietal lobe) and temporal (volume-weighted average of parahippocampal gyrus, hippocampus, medial temporal lobe, superior, middle and inferior temporal gyrus) regions.26

All SUVr images were visually read by an experienced nuclear physician (BvB), and a positive or negative rating was awarded according to criteria defined by the manufacturer (GE Healthcare).27

Statistical analysis

The values for BP ND were log-transformed to normalise their distribution. Therefore, analyses with these as dependent variables are reported in ratios instead of regression coefficients. A ratio can take any value from 0 onward, where 1 means there is no relationship, <1 means there is a negative relationship and >1 means there is a positive relationship. Associations between amyloid positivity/BP ND and vessel density/FAZ were analysed using generalised estimating equations (GEE) to correct for clustering in the data from twin pairs as well as confounders such as age, gender, scan quality and a diagnosis of diabetes. As multiple yet related variables (vessel density in several regions and FAZ) were analysed, data were corrected for multiple testing using the Holm-Bonferroni method.28 Receiver operating characteristic (ROC) curves were used to test overall predictive accuracy. To estimate the maximal genetic contribution to OCTA variables, the resemblance of MZ twins was quantified using Pearson correlations. For OCTA measures that were associated with amyloid status, post hoc analyses were performed, which used the status of the twin pair as a measure of ‘Aβ load’. We defined four groups, depending on concordance or discordance for Aβ status within pairs.29 Twins were grouped with increasing Aβ load starting with twins from pairs who were both Aβ– (group 1=concordant negative, ie, the ‘healthy group’), followed by twins from discordant pairs who were Aβ– (group 2=discordant negative) or Aβ+ (group 3=discordant positive) and ending with twins from pairs who were both Aβ+ (group 4=concordant positive, ie, the group with the highest load). This classification assumed that discordant Aβ– participants were in a very early stage of amyloid accumulation. The resulting variable with four different loads for Aβ pathology was analysed with GEE to look for associations with OCTA markers. All statistical analyses were performed in SPSS (V.22).

Results

Of the 150 participants scanned, 26 were excluded due to ophthalmological disease (N=12, mostly glaucoma and age-related macular degeneration) or low-quality scans of both eyes (N=14). table 1 shows demographic information of the 124 participants included in the analyses.

Table 1

Demographics of the study population

Number of participants (N)	124
Age (years)	68.6 (±6.3)
Sex, female N (%)	66 (53.2)
BCVA (logarithm of the minimum angle of resolution)	0.01 (±0.09)
IOP (mm Hg)	14.2 (±2.6)
Spherical equivalent	0.40 (±1.82)
MMSE (median, IQR)	29.0 (29.0–30.0)
Brain amyloid status, positive N (%)	13 (10.5)
Global BP ND of Aβ (median, IQR)	0.122 (0.095–0.177)
Retinal vessel density
Inner ring macula	17.3 (±1.3)
Outer ring macula	17.6 (±1.2)
Around optic nerve head	18.1 (±1.3)
Foveal avascular zone (mm2)	0.22 (±0.09)

Data are means unless otherwise specified.

Aβ, amyloid-beta; BCVA, best-corrected visual acuity (both eyes averaged); BPND, non-displaceable binding potential; IOP, intraocular pressure (both eyes averaged); MMSE, Mini-Mental State Exam.

Vessel density in all regions was significantly higher for Aβ+ individuals versus Aβ– individuals. There was no significant difference in FAZ size between Aβ+ and Aβ– individuals. table 2 shows all vessel density and FAZ differences for both groups.

Table 2

Difference in vessel densities and foveal avascular zone area between amyloid-positive and amyloid-negative participants

	Aβ+ versus Aβ–
	Difference	Raw p value	95% CI	Adjusted p value
Retinal vessel density
Inner ring macula	0.81	<0.001	0.36 to 1.26	0.002
Outer ring macula	0.50	0.012	0.11 to 0.90	0.024
Around optic nerve head	0.83	0.005	0.25 to 1.42	0.015
Foveal avascular zone (mm2)	−0.02	0.387	−0.05 to 0.02	0.387

Adjusted p values were calculated using the Holm-Bonferroni method.

GEE, corrected for age, gender, scan quality and diagnosis of diabetes.

Aβ, amyloid-beta; GEE, generalised estimating equation.

The Area Under the Curve (AUC) values corresponding to the ROC curves for vessel densities in inner and outer ring of the macula and around the ONH were 0.651, 0.640 and 0.764, respectively (online supplementary figure 1). Sensitivity and specificity were both around 0.8 at a cut-off of 18.9 for the vessel density around the ONH.

Table 3 shows the linear twin post hoc analyses using the Aβ load score for vessel densities in the three regions; they were significant for vessel densities all regions. Online supplementary table 1 shows these analyses when comparing group means, and online supplementary figure 2 shows the boxplots.

Table 3

Twin post hoc Aβ load analyses for retinal vessel densities in three regions

Retinal vessel density	Twin Aβ load score
	Beta	Raw P value	95% CI	Adjusted Pp value
Inner ring macula	0.27	0.009	0.07 to 0.47	0.018
Outer ring macula	0.18	0.037	0.01 to 0.35	0.037
Around optic nerve head	0.30	0.002	0.11 to 0.49	0.006

Obtained using GEEee, corrected for age, gender, scan quality and diagnosis of diabetes.

Aβ, amyloid-beta; GEE, generalised estimating equations.

The continuous measure of BP ND of [18F]flutemetamol to Aβ was significantly associated with vessel density in the inner ring, but this association did not retain statistical significance after correction for multiple testing (table 4).

Table 4

Associations between global cortical non-displaceable binding potential of Aβ and vessel densities/foveal avascular zone, given in ratios due to a log transformation applied to the dependent variable

	Global cortical non-displaceable binding potential of Aβ
	Ratio	Raw p value	95% CI	Adjusted p value
Retinal vessel density
Inner ring macula*	2.28	0.016	1.21 to 6.41	0.064
Outer ring macula*	3.74	0.052	0.99 to 14.1	0.156
Around optic nerve head*	0.86	0.654	0.43 to 1.69	0.654
Foveal avascular zone (mm2)	0.59	0.356	0.19 to 1.80	0.712

GEE, corrected for age, gender, scan quality and diagnosis of diabetes.

*Reported in steps of 10. Adjusted p values were calculated using the Holm-Bonferroni method.

Aβ, amyloid-beta; GEE, generalised estimating equations.

All OCTA measures correlated moderately to highly within twin pairs, except for vessel density around the ONH (table 5). Correction for APOE genotype (APOE4 carrier vs non-carrier) did not change these correlations.

Table 5

Intratwin pair correlation coefficients

	Correlation coefficient	P value
Retinal vessel density
Inner ring macula	0.50	<0.001
Outer ring macula	0.62	<0.001
Around optic nerve head	0.34	0.014
Foveal avascular zone	0.85	<0.001

Discussion

The main finding of this study was that Aβ+ participants had a significantly higher vessel density on OCTA than Aβ– participants in both the inner and outer ring of the macula and around the ONH. No significant differences in FAZ size were found.

Earlier studies demonstrated a lower retinal vessel density and larger FAZ in patients suffering from AD in the dementia stage.16–18 Based on these findings, we anticipated similar differences in the preclinical phase of AD. Yet, we found the opposite: a higher vessel density in all retinal areas (inner and outer ring of the macula, around the ONH), remaining statistically significant after correction for multiple testing. An explanation for this unexpected finding may lie in an inflammatory state of the retina in the early stages of amyloid accumulation, as many studies have found that the cerebral events occurring during the development of AD and accumulation of Aβ are often of an inflammatory nature.30 Supposing such events occur in synchrony within the retina, it may be that in the preclinical stages this inflammatory reaction with hypoxia leads to increased retinal blood flow. Microvessels normally not detected on OCTA due to a blood flow, which is below the detection level, will become visible. This increase in the number of microvessels detected on OCTA will result in a higher vessel density.31 After this initial phase, continued inflammation and Aβ build-up may cause further damage resulting in a loss of (micro)vasculature in turn causing a decrease in vessel density, explaining why individuals suffering from full-blown AD may have lower vessel density.16 17 32

In our cohort, no differences were found in FAZ size between Aβ+ and Aβ− individuals. This is in contrast with the findings of a previous study, where a FAZ enlargement was found in Aβ+ participants.33 Two other studies looking at AD also found FAZ enlargement.16 17 It is possible that our participants were at an earlier stage in the disease process, where FAZ size is still unchanged but that does not explain why the other group already found this difference in preclinical AD cases. The FAZ size being unchanged does support our theory that in an early inflammatory process, blood flow in microvasculature may be increased, as one would not expect an enlargement of FAZ if this were the case. Further studies need to be performed to elucidate possible changes in FAZ size in preclinical AD.

Considering the rather pronounced differences in vessel density between Aβ+ and Aβ− individuals, ROC curves for these measures were generated (online supplementary figure 1). Due to the low N in the Aβ+ group, caution is warranted in their interpretation. Vessel density around the ONH showed the best ROC curve (AUC=0.764, online supplementary figure 1), and with a sensitivity and specificity of around 0.8 suggests diagnostic potential, although findings need to be replicated.

When staging concordant and discordant twins in a presumed model of increasing amyloid accumulation (twin post hoc Aβ load analysis29) vessel density increased further along the Aβ load groups in a more or less linear fashion (online supplementary table 1, online supplementary figure 2). A statistically significant positive association between this Aβ load score and vessel density was found for all regions (table 3), supporting our theory that with more Aβ pathology, retinal vessel density increases.

We found an association between BP ND and vessel density in the inner ring of the macula, and a trend with vessel density in the outer ring of the macula, but both lost statistical significance after correction for multiple testing. Considering the differences between the two groups (Aβ+ vs Aβ−) in vessel density, we had expected a stronger association between retinal vessel density and BP ND (as a continuous measure for cortical Aβ). The explanation for the lack of this (strong) association may be found in the distribution of BP ND; many people have little to no cortical Aβ, whereas just a few individuals have high cortical Aβ. This results in little variation in BP ND, which might explain why no associations were found.

Intratwin pair correlations were moderate to high for all OCTA measures except vessel density around the ONH (table 5). In earlier studies, we found higher intratwin pair correlations for other retinal markers such as retinal vascular parameters on fundus photography.34 The reason for these lower correlations of OCTA markers may lie in lower reproducibility of vessel density measures, as we noticed that vessel density is very much affected by scan quality (lower quality=lower vessel density) as well as clouding of the media (floaters and cataract), making this a less reliable measure. This also explains the high number of participants who were excluded due to insufficient quality scans of both eyes (N=14). Although the intratwin pair correlations were lower for OCTA measures, they still suggest a moderate to high contribution of either genes or shared environment to these traits. As no dizygotic twin pairs were included in this study, no heritability estimates could be made. However, the estimated contribution of shared environmental effects is lower in older twins, as these twins have spent most of their life apart.29 35 Therefore, due to the older age of this cohort, the twin correlations in OCTA measures are likely of a genetic origin. Nevertheless, a substantial amount of variability (39%–66%, table 5) is related to unique environmental factors, which remain to be identified. To the best of our knowledge, no other studies on the subject of heritability of OCTA measures have been performed.

Strengths of the study include the thorough characterisation of participants from both a neurological and ophthalmological perspective. Extensive screening for cognitive health guaranteed that only cases of preclinical AD (ie, cognitively healthy) were included. We excluded (glaucoma, vascular occlusions and age-related macular degeneration) or controlled (diabetes) for confounders known to affect vessel density.

The main limitation of this study lies in its low inclusion and resulting (statistical) power, especially for Aβ+ individuals. This is due to the nature of the study: only around 20% of individuals around ages 70 years are positive for Aβ, making this a common problem in studies assessing preclinical AD.5 Yet, several OCTA measures showed significant differences between the groups, which remained after correction for multiple testing. However, the low number of participants warrants caution in interpreting these results, as the results found in this (small) subset of participants might not be applicable to the population as a whole. Further studies with a larger number of subjects are needed to confirm these findings. Additionally, OCTA is a relatively new technology, and the parameters that are currently included in the imaging software (vessel density and FAZ size) may end up not being the most useful metrics for biomarker assessment. Future studies should also include newly developed, and perhaps more suitable, OCTA measures to examine their feasibility as potential biomarkers for preclinical AD.

Another limitation was that OCTA scans were made 22 months after the PET scans. At that time, four individuals had a CDR score of 0.5, indicating mild cognitive impairment. Exclusion of these individuals yielded similar results. It is possible that some individuals may have become amyloid positive at follow-up, but a previous study in amyloid negative, cognitively normal individuals showed that the risk for this in 2 years is very low and therefore unlikely to have an impact on our findings.36

In conclusion, we found a significantly higher retinal vessel density in all measured regions in Aβ+ versus Aβ− participants, suggesting retinal vessel density to be of interest as a possible biomarker for preclinical AD. No differences were found in FAZ size between the two groups.