Abstract

Introduction. Retinal structural and microvascular damages reflect damage to cerebral microvasculature and neurons. We aimed to investigate neovascular unit abnormalities among patients with large-artery atherosclerosis (LAA) or small-vessel occlusion (SAA) and control subjects. Methods. Twenty-eight LAA patients, forty-one SAA patients, and sixty-five age- and gender-matched controls were recruited. Based on optical coherence tomography angiography (OCTA), retinal capillary vessel density was assessed in the general and local sectors, and the thickness of individual retinal layer was extracted from retinal structural images. The differences between structural and microvascular were analyzed. Results. The superior peripapillary retinal nerve fiber layer (pRNFL) thickness was significantly different among the three groups, and the LAA group had the thinnest thickness. Compared to the control group, the deep retinal capillary vessel density in other two stroke subgroups were significantly reduced in all regions except in the inferior region ( < 0.05), and the fractal dimension in C2 and C4 regions of deep retina was significantly lower in the LAA group ( < 0.05). Discussion. Compared with superficial microvascular network, deep microvascular network is more sensitive to ischemic stroke. In addition, we have demonstrated quadrant-specific pRNFL abnormalities in LAA and SAA patients. Superior quadrant pRNFL thickness differences between stroke subgroups may suggest that changes in retinal nerve fiber layer are more sensitive to subtype identification than changes in retinal microvascular structure. All in all, the alteration in retinal structural and microvascular may further elucidate the role of the neovascular unit in ischemic stroke, suggesting that the combination of these two indicators could be used for subtype identification to guide prognosis and establish a risk prediction model.

1. Introduction

Stroke is the most common cause of serious disability in adults, and China bears the biggest burden globally [1]. As therapeutic options are limited, effective preventive strategies for early diagnosis are needed. The underlying subclinical pathologic process occurs much earlier before the onset of clinical stroke while current neuroimaging technologies may not be capable of directly observing subtle subclinical changes in stroke due to the resolution. Besides, current predictions of stroke are difficult to quantify. Therefore, there is an urgent need for additional surrogate techniques to detect the subtle changes in vivo. In addition to the importance of finding indicators to establish a model for stroke risk prediction and prognosis assessment, the identification of subtypes of ischemic stroke is vital for guiding clinical treatment and management. However, the current international classification Treatment of Acute Stroke Trial (TOAST) [2] requires a lot of auxiliary examinations which are expensive and time-consuming. The above shortcomings and possible examination contraindications make it challenging for patients to complete the examination in time when they are admitted to the hospital, affecting the early guiding role in clinical treatment. Therefore, finding an early, especially sensitive and effective, method is of great importance in disease prediction, treatment, and prognosis assessment. Because the retinal and cerebral vessels share similar anatomic, embryological, and physiological characteristics, the retina provides a unique “window” to assess the cerebral microvascular and neurons in vivo noninvasively.

Previous studies based on fundus photography have revealed an independent correlation between retinal vascular parameters and stroke [36]. Additionally, several studies also detected vascular changes vary according to stroke subtypes, suggesting the specific cerebral microvasculopathy of subtyping of stroke [7, 8]. Inconsistently, another study showed that vascular changes were similar between stroke subtypes [9]. The reason for the discordant results may be that the fundus photography is only a plane picture, which only reflects the large blood vessels of the retina without microvascular and quantitative retinal structural parameters. With the advancement of the technique, the optical coherence tomography angiography (OCTA) can reflect finer retinal capillary plexuses and choriocapillaris changes by generating three-dimensional images based on the comparisons of the motion of circulating blood cells. Consequently, we can observe retinal microvascular changes of stroke patients.

The concept of the neurovascular unit (NVU) was proposed in 2003 [10], which is composed of endothelial cells, neurons, astrocytes, and pericytes. The NVU provides new insights for the pathogenesis and diagnostic and treatment strategies of stroke [11, 12]. Spectral-domain optical coherence tomography (SD-OCT) with high-resolution retinal imaging can provide cross-sectional images of biologic structures and quantify the thickness of each retinal layer. One previous study observed that the transneuronal retrograde degeneration (TRD) of retinal ganglion cells (RGCs) assessed by SD-OCT is associated with cerebral infarction [13]. Most of the studies containing previous stroke patients with an increased risk of confounding factors concentrated on the association between large retinal vessels or neural structures changes and stroke separately [3, 13]. So far, there have been no in vivo studies on simultaneous microvascular and neural structures in stroke subjects.

In this study, we aimed to find the retinal microvascular and microstructural changes in the subtypes of initial acute stroke patients.

2. Methods

2.1. Study Population

In this study, a total of 85 patients with initial ischemic stroke within 14 days of an acute period were prospectively recruited from the neurology unit at the Second Affiliated Hospital & Yuying Children’s Hospital of Wenzhou Medical University from Jan 2017 to December 2018. One neurologist (Zhao Han) assessed the stroke severity with the National Institutes of Health Stroke Scale (NIHSS) [14] and classified large-artery atherosclerosis (LAA) and small-artery occlusion lacunar (SAA) according to a modified TOAST classification [2]. Besides, 65 age- and gender-matched controls with no self-reported history of stroke or transient ischemic attack or ophthalmic disease were enrolled consecutively from the relatives of patients or working staff at the Eye Hospital or the Second Affiliated Hospital & Yuying Children’s Hospital of Wenzhou Medical University between Jan 2017 to Aug 2019. Considering the effect of stroke site, parameters of ipsilateral eyes were selected for analysis. Additionally, random eyes were selected in nonunilateral stroke and control subjects. Written informed consent was obtained from patients or their next of kin, and the project was approved by the ethics committee of the Eye Hospital of Wenzhou Medical of Wenzhou Medical University.

2.2. Assessment of Cardiovascular Risk Factors

Patients finished detailed questionnaires for information on history of hypertension, diabetes mellitus, hypercholesterolemia, ischemic heart disease, cigarette smoking status, and medication use. All patients underwent usual examinations of stroke, including brain imaging, fasting blood samples for glycosylated hemoglobinA1C (HbA1C), total cholesterol (TC), total triglycerides (TG), homocysteine (HCY), creatinine (Cr), high-density lipoprotein (HDL-C), low-density lipoprotein (LDL-C), and body mass index (BMI). Besides, as a part of clinical care for stroke, blood pressure was measured three times after participants had been seated for at least 10 minutes at the same sitting. The mean blood pressure of three times was taken as the final result. The mean arterial pressure (MAP) is equal to one-third of systolic blood pressure (SBP) plus two-thirds of diastolic blood pressure (DBP).

Hypertension was diagnosed as SBP ≥140 mm Hg or DBP ≥90 mm Hg at examination or a self-reported history of physician-diagnosed hypertension or the use of antihypertensive medication. Diabetes mellitus was defined as fasting blood glucose ≥7.0 mmol/L and/or random blood glucose ≥11.1 mmol/L, hemoglobin A1C ≥7%, self-reported history of physician-diagnosed diabetes mellitus, or the use of antihyperglycemic medication. Hypercholesterolemia was defined as fasting total cholesterol ≥5.2 mmol/L, self-reported history of physician-diagnosed hypercholesterolemia, or the use of antilipemic medication. Current smokers were defined as people who smoke currently or quitted smoking less than one year before the examination.

2.3. Assessment of Ophthalmic Parameters and Microstructure of the Retina

All the subjects had detailed ophthalmologic examinations performed by two ophthalmologists (Ying Zhang and Ce Shi), including slit-lamp biomicroscopy, refraction diopter, best-corrected visual acuity (BCVA), and noncontact intraocular pressure (IOP). All patients were imaged by an OCT system (Optovue RTVue-XR Avanti; Optovue, Inc., Fremont, CA, USA) to obtain the OCTA images, with detailed steps described as follows: refraction data were converted to spherical equivalents (SEs) and calculated as the spherical dioptric power plus one-half of the cylindrical dioptric power. The exclusion criteria were presented as follows: patients with contraindications to magnetic resonance, hemorrhagic stroke, recurrent stroke, those who were unable to complete the eye examinations, and those with spherical equivalent (SE) under ± 5.00 D, IOP >21 mm Hg or previous ophthalmologic diseases (such as cataract, glaucoma, high myopia, and retinal diseases). Other exclusion criteria were systemic diseases that could affect the ocular structures, such as uncontrolled hypertension/diabetes and neurological diseases such as Parkinson’s disease and multiple sclerosis.

2.4. MRI Analysis

All patients underwent the 3.0-T MRI (Signa HDxt GE Healthcare), which included T1-weighted and T2-weighted imaging, diffusion weighted imaging (DWI), and fluid attenuated inversion recovery (FLAIR). The slice thickness was 5 mm with an interslice gap of 1 mm. The high signal on DWI sequence of MRI indicates the presence of acute cerebral infarction. In addition, the size and the location of the lesion are conductive to the classification of stroke. In the control group, MRI was also used for the reason of homogeneous management. For ethical reasons, not all controls were willing to undergo the tests for fasting blood tests and magnetic resonance imaging. Therefore, the two indicators are not shown in the Table 1.

2.5. OCT and OCTA Acquisitions

All subjects remained seated under the same conditions, and examinations were performed by an expert examiner. The OCTA system, which employs the split-spectrum amplitude-decorrelation angiography (SSADA) algorithm, operated at a rate of 70000 A-scans per second, with the scan area of 3 × 3 mm2, and the results were obtained by orthogonal registration and merging of two consecutive B-scans. The size of the exported OCT images was 304 × 304 pixels. OCTA combines orthogonal fast-scan directions to correct motion artifacts based on the DualTrac Motion Correction technology [15]. A good set of scans with a signal strength index (SSI) over 40 was selected for further analysis.

2.5.1. Retinal Layer Thickness Analysis on Spectral-Domain Optical Coherence Tomography

Retinal thickness was imaged by the RTVue XR Avanti SD-OCT system (Optovue, Inc, Fremont, California, USA). Besides, the average, superior (S), temporal (T), inferior (I), and nasal (N) quadrants of peripapillary retinal nerve fiber layer (pRNFL) thickness were obtained. The ganglion cell complex (GCC) provides inner retinal thickness values from the internal limiting membrane (ILM) to the inner molecular layer (IPL), shown as average, superior, and inferior regions (Figure 1(f)).

2.5.2. Capillary Vessel Density and Fractal Dimension Analysis on Optical Coherence Tomography Angiography

Retinal microvasculature was evaluated by the RTVue XR with AngioVue (software version 2017.1.0.155; Optovue, Inc, Fremont, CA, USA). Patients underwent different types of scanning: 3 × 3 mm2 angioretina around the fovea, 4.5 × 4.5 mm2 angioretina scans around the optic nerve head (ONH), 3D retina scans, retinal map scans, and radial lines scans.

Vessel density (VD) is defined as the percentage of area occupied by OCTA detected vasculature. The software sets the superficial capillary plexuses (SCP) from 3 μm below the ILM to 15 μm below the IPL. The deep capillary plexuses (DCP) were set from 15 to 70 μm below the IPL (Figures 1(b) and 1(c)). In addition, the parafovea vessel density, defined as the area of annular circle with a diameter of 3 mm excluding the fovea zone (diameter = 1 mm), was divided automatically into whole and superior (S), temporal (T), inferior (I), and nasal (N) quadrants. Similarly, the 5 sectors of radial peripapillary capillary (RPC) vessel density were analyzed. The boundary of RPC ranges from ILM to the nerve fiber layer.

To quantify the complexity of the branching pattern and density of the retinal capillary network in OCTA images, the automated fractal analysis system was employed to correct the image magnification based on the axial length [16, 17]. Briefly, the OCTA images in PNG format were imported to the custom automated algorithm software published previously [18]. Then, the grayscale of the two-dimensional OCTA images was first extended by bicubic interpolation to 1024 × 1024 pixels so as to improve the image details. The binary images of vessels were created by the algorithm. Subsequently, one binary image containing only large arteries and the other binary image containing both large and small vessels were subtracted to obtain the final binary image. Based on the final image of white-pixelated vasculature, a skeletonized image was created by detecting the central axis of each capillary. After the image processing, both the superficial and deep retinal capillary complexities were calculated based on the skeletonized images [19, 20]. The quantitative measured parameter of complexity, Dbox values, was obtained with the fractal analysis software (Benoit, Trusoft Benoit Fractal Analysis Toolbox; Trusoft International, Inc., St. Petersburg, FL). Both the general and local fractal dimensions were used to describe the complexity of capillary network. At first, after excluding the fovea avascular zone (FAZ) within the diameter of 0.6 mm, the fractal dimension (FD) was automatically calculated for the total annular zone (TAZ) within the 2.5 diameter and for the 4 parafoveal quadrant sectors (superior (S), temporal (T), inferior (I), and nasal (N)) and 6 concentric isometric annular rings (Figures 1(d) and 1(e)). The methods above were implemented using MATLAB v 7.10 (Mathworks, Inc., Natick, Massachusetts, USA).

2.6. Statistical Analysis

All statistical analyses were conducted using SPSS software (version 24.0; SPSS, Inc., Chicago, IL, USA). The data were expressed as the mean ± standard deviation (SD). One-way analysis of variance (ANOVA) was used to test the differences among patients with large-artery stroke and lacunar stroke and control subjects, and Bonferroni correction was used for pairwise comparisons. The differences in gender and medical history were determined by the χ2 test.

3. Results

A total of 85 patients were included in the present study. Among them, 16 patients were excluded due to poor image quality of OCTA scans. The remaining 28 patients with LAA, 41 patients with SAA, and 65 age-and gender-matched control subjects were further analyzed. The demographic and clinical characteristics are summarized in Table 1. Normally, distributed data are represented by mean plus or minus standard deviation, while nonnormally distributed data are represented by median and interquartile spacing. Differences in age, sex, BMI, SE, IOP, and DBP, together with the prevalence of hyperlipidemia, currently smoking, and previous ischemic heart disease were not significant among the three groups. Patients with SAA were more likely to have hypertension than the other two groups ( < 0.001). Besides, the values of SBP and MAP were significantly higher than those of the other two groups ( < 0.001). LAA patients were more likely to undergo diabetes ( = 0.007) than the controls. However, there was no statistical difference in demographic data between the two stroke subgroups. Not all subjects completed all tests. The number of eyes that completed each examination is shown in Table 2.

3.1. Intergroup Differences among the Three Groups
3.1.1. Retinal Microstructural Thicknesses

In terms of the quadrants, the superior pRNFL thickness was significantly thinner in the eyes of patients than that in the eyes of the control group ( = 0.01, Table 3, Figure 2(a)). In the eyes of LAA, the pRNFL thickness was significantly thinner in the superior quadrant compared to the eyes of SAA ( = 0.034) and control ( = 0.003). No significant superior pRNFL thinning was observed in SAA compared to control ( = 0.438). Additionally, the thicknesses of GCC were not significantly different among the three groups (all  > 0.05, Table 3, Figure 2(b)).

3.1.2. Vessel Density around ONH and Macula

The vessel density around the macula in the deep retinal capillary layer was significantly reduced in patients with LAA or SAA within all regions, except for the inferior region ( < 0.05, Table 4, Figure 3(b)). The significant differences mostly existed between the stroke group and control group ( < 0.05), and no significant difference was found in the two stroke subgroups, although the LAA group tended to have a lower vessel density. In addition to ONH capillary density, no significant difference was found in the superficial layer of all regions ( > 0.05, Table 4, Figure 3(a)).

3.1.3. Fractal Dimension around Macula

Differences in the fractal dimension were only statistically significant in C2 and C4 regions of the deep retina ( < 0.05, Table 5, Figure 4(b)) between patients with LAA and the controls. Compared with the control group, the fractal dimension of most regions in the stroke group showed a downward trend, although there was no statistical difference in most regions ( > 0.05, Table 5, Figures 4(a) and 4(b)).

4. Discussion

Assuming that the retinal vasculature mirrors the cerebral vasculature, OCTA enables noninvasive imaging of retinal capillaries in multiple layers invisible on fundus images. Based on fundus photography, retinal abnormalities, including arteriovenous nicking and generalized and localized arteriolar thinning, a lower arteriolar/venular diameter ratio and geometric parameters have been demonstrated to be significantly related to the incidence of stroke [5, 9, 21, 22]. The previous studies using fundus photos could not qualitatively detect the subtle changes at the capillary level while OCTA provides the opportunity to investigate the retinal capillary microcirculation at the micrometer resolution [23].

In terms of vessel density around the macula, our findings demonstrate that the changes of vessel density are more obvious in the deep layers than the superficial layer between stroke patients and control subjects. Our results point to a preferential involvement of the deep layer in patients with stroke which may be attributed to the fact that the deep network consists of a dense and complex system of smaller vessels [24]. It can be speculated that deep retinal capillaries are more susceptible to ischemic and hypoxia. As per our expectation, patients with SAA displayed much smaller changes than patients with LAA. However, there exists no significant difference between the LAA and SAA. The finding suggested that the retinal vasculopathy may result from downstream effects of large-artery pathology in the cerebral circulation. In addition, we found that fractal dimension was not helpful in identifying stroke subtypes and was significantly lower in the LAA group than in controls only in individual regions of the deep retina. Our results were in consistence with those of others [9], which found that decreased FD was correlated with stroke, suggesting a loss of complexity. However, previous results of fractal dimension based on the fundus photograph remain controversial. Some considered lacunar stroke subtype was associated with decreased retinal FD [8] while others demonstrated that the lacunar stroke was positively associated with higher FD [7]. We found no significant difference between stroke subgroups, which may indicate that the FD is not applicable to differentiate the disease subgroups who have already had stroke.

Regarding peripapillary vessel density, there was no significance in RPC among the groups. The finding could be related to the anatomical differences between the different areas. The parafoveal superficial capillary plexus originates largely from the retinal circulation, whereas the RPC receives additional blood supply from the choroid [25]. Larger vascular channels around the optic disc may have masked subtle changes in the capillary network. Additionally, RPC contains multiple layers of capillaries that overlap on en-face OCTA images, lacking the ability to detect tiny vascular losses.

In addition to retinal capillary changes, our study observed that the pRNFL thickness was statistically reduced in the superior quadrant of stroke patients, and there was a statistical difference between the stroke subgroups which may indicate different patterns of nerve damage in the two stroke subtypes. Additionally, some research recently also observed that both acute and previous stroke were significantly associated with retinal nerve fiber layer defects (RNFLDs) [13]. These findings were also in accordance those of with others [26], which reported the associations between severity and laterality of RNFLD and laterality of hemispheric damage as well as arterial territory of infarct. They found RNFLDs were significant in the temporal sector of the ipsilateral side and in the nasal sector of the contralateral side of the stroke. Furthermore, they also confirmed that the degree of the transneuronal retrograde degeneration (TRD) was time-dependent. However, we find that the significant RNFLDs only exist in the superior sector in our study although the ipsilateral sides were included.

As is known to all, over 30 morphological types of RGCs compose the structure of the retina. The midget RGCs (80%), with wide retinal dendritic fields located in the peripheral retina, mainly project to the magnocellular layers of the lateral geniculate body. The parasol RGCs (5–15%) predominantly present in the papillomacular bundle (macula) and project to suprachiasmatic nucleus of the hypothalamus. Axons from midget RGCs enter the superior, inferior, and nasal sectors of the optic nerve, and parasol RGCs gather their axons and enter the temporal sector of the optic nerve [27, 28]. The different degenerative patterns of the ganglion cells and nerve fibers have been demonstrated in several neurodegenerative diseases [29, 30]. Recent studies have described that the parasol RGCs are more involved in the pathogenesis of Alzheimer’s disease [3133]. However, as for the Parkinson syndrome and mitochondrial optic neuropathy, the midget RGCs are mainly involved [30, 34, 35]. Both the above and our study suggest that mean thickness measurement may not reflect the disease well and reduce the diagnostic efficacy of the ocular biomarker. We speculate that detailed analysis of the focal nerve structure alteration may be developed as an ocular imaging biomarker for monitoring disease progression and evaluating prognosis of these diseases.

Nevertheless, as this is the first study on ischemic stroke subtypes based on OCTA, further large samples are required to confirm the universality in the disease. The discrepancies existing between different studies could be ascribed to differences in course and severity of disease as well as various OCT devices and study designs.

Besides, a previous animal study also revealed that the TRD of retinal ganglion cells occurred after shrinkage of the optic tract, and degeneration of the RGC was progressing slowly in the next few years [36]. Thus, our negative results in GCC thickness might be related to the fact that patients were tested within two weeks after stroke.

To conclude, our study provides a side view that the neurovascular unit is affected in ischemic stroke and severely in LAA patients. The components of neurovascular units are interrelated in the microenvironment. Studies have shown that the signal transmission between neurons, astrocytes, and microvascular endothelial cells regulates the brain microenvironment [37] [29]. Briefly, NVU as a structural and functional whole, the relationship between members changes in the state of illness. The mechanism of cerebral ischemia is complex and involves multiple cascading reactions. Therefore, monitoring the NVU as a whole and improving its function help maintain brain cell function and make stroke treatment more ideal. The retina provides a visualization of the neurovascular units simultaneously to reflect the changes of the brain, which has important implications for disease surveillance.

In addition, we also acknowledge the limits in this study. The cross-sectional study with a small sample limits our ability to identify the different pathogenesis of the capillaries and microstructures of the stroke. The outer structure of the retina which is mainly supplied by the choroid has not been analyzed. Due to the practical difficulty of recruiting and examining the patients, milder patients may be recruited, which may be more likely with milder vascular lesions to differentiate the two subtypes. Despite the weakness, there are several strengths in our study. We recruited different ischemic stroke subtypes with strict inclusion and exclusion criteria. Otherwise, to minimize the confounding of stroke, the patients with a history of previous stroke were excluded. We completed the detailed ophthalmic and clinical examinations within two weeks after the onset of the stroke and maintained the blinding of retinal and brain images to each other. Finally, this is the first attempt to simultaneously observe retinal microvascular and neurological changes in different stroke subtypes in vivo.

Moreover, further longitudinal studies with a larger sample size are required to show differences in the characteristics of retinal microstructure and capillaries. It remains to be seen whether the retinal signs are indicative of cerebrovascular risk beyond the conventional risk indicators and whether the retinal imaging technique will be a surrogate or accessory examination system in clinical settings and ultimately become a part of the routine stroke risk or evaluation of treatment assessment.

Data Availability

All data generated or analyzed during this study are included within published article.

Ethical Approval

This study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee Board of the Wenzhou Medical University (2019-027-K-26).

All subjects, recruited voluntarily, were informed about the purposes, methods, and potential risks of the study. A signed consent form was obtained from each patient.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Ying Zhang, Fan Lu, and Meixiao Shen designed the study. Ying Zhang, Fan Lu, and Ce Shi conducted the study. Ying Zhang, Xianda Lin, Weicheng Wang, and Shenghai Huang collected the data. Shenghai Huang wrote the image processing algorithm for data analysis. Ying Zhang and Ce Shi analyzed and interpreted the data. Fan Lu and Meixiao Shen were the main contributors to manuscript discussion. Meixiao Shen, Weicheng Wang, and Zhao Han revised the manuscript. All authors read and approved the final manuscript.

Acknowledgments

This study was supported by the research grants from the National Key Project of Research and Development Program of Zhejiang Province (Grant no. 2019C03045), Wenzhou Municipal Science and Technology Bureau (Grant no. 2018ZY016), and Zhejiang Provincial Natural Science Foundation of China (Grant no. LQ17H180004).