Physiological changes in retinal layers thicknesses measured with swept source optical coherence tomography

Abstract

Purpose

To evaluate the physiological changes related with age of all retinal layers thickness measurements in macular and peripapillary areas in healthy eyes.

Methods

Wide protocol scan (with a field of view of 12x9 cm) from Triton SS-OCT instrument (Topcon Corporation, Japan) was performed 463 heathy eyes from 463 healthy controls. This protocol allows to measure the thickness of the following layers: Retina, Retinal nerve fiber layer (RNFL), Ganglion cell layer (GCL +), GCL++ and choroid. In those layers, mean thickness was compared in four groups of ages: Group 1 (71 healthy subjects aged between 20 and 34 years); Group 2 (65 individuals aged 35–49 years), Group 3 (230 healthy controls aged 50–64 years) and Group 4 (97 healthy subjects aged 65–79 years).

Results

The most significant thinning of all retinal layers occurs particularly in the transition from group 2 to group 3, especially in temporal superior quadrant at RNFL, GCL++ and retinal layers (p≤0.001), and temporal superior, temporal inferior, and temporal half in choroid layer (p<0.001). Curiously group 2 when compared with group 1 presents a significant thickening of RNFL in temporal superior quadrant (p = 0.001), inferior (p<0.001) and temporal (p = 0.001) halves, and also in nasal half in choroid layer (p = 0.001).

Conclusions

Excepting the RNFL, which shows a thickening until the third decade of life, the rest of the layers seem to have a physiological progressive thinning.

About the Authors:

Elisa Viladés

Roles Data curation, Investigation, Writing – original draft

Affiliations Department of Ophthalmology, Miguel Servet University Hospital, Zaragoza, Spain, Aragon Institute for Health Research (IIS Aragón), University of Zaragoza, Zaragoza, Spain

Amaya Pérez-del Palomar

Roles Formal analysis, Software

Affiliations Aragon Institute for Health Research (IIS Aragón), University of Zaragoza, Zaragoza, Spain, Mechanical Engineering Department, Aragon Institute of Engineering Research, University of Zaragoza, Zaragoza, Spain

José Cegoñino

Roles Formal analysis, Software

Affiliations Aragon Institute for Health Research (IIS Aragón), University of Zaragoza, Zaragoza, Spain, Mechanical Engineering Department, Aragon Institute of Engineering Research, University of Zaragoza, Zaragoza, Spain

Javier Obis

Roles Supervision, Writing – review & editing

Affiliation: Department of Ophthalmology, Miguel Servet University Hospital, Zaragoza, Spain

María Satue

Roles Supervision, Writing – review & editing

Affiliations Department of Ophthalmology, Miguel Servet University Hospital, Zaragoza, Spain, Aragon Institute for Health Research (IIS Aragón), University of Zaragoza, Zaragoza, Spain

Elvira Orduna

Roles Data curation, Resources, Writing – review & editing

Affiliations Department of Ophthalmology, Miguel Servet University Hospital, Zaragoza, Spain, Aragon Institute for Health Research (IIS Aragón), University of Zaragoza, Zaragoza, Spain

ORCID logo http://orcid.org/0000-0003-2710-1875

Luis E. Pablo

Roles Supervision, Writing – review & editing

Affiliations Department of Ophthalmology, Miguel Servet University Hospital, Zaragoza, Spain, Aragon Institute for Health Research (IIS Aragón), University of Zaragoza, Zaragoza, Spain

Marta Ciprés

Roles Writing – review & editing

Affiliations Department of Ophthalmology, Miguel Servet University Hospital, Zaragoza, Spain, Aragon Institute for Health Research (IIS Aragón), University of Zaragoza, Zaragoza, Spain

Elena Garcia-Martin

Roles Formal analysis, Writing – review & editing

* E-mail: [email protected]

Affiliations Department of Ophthalmology, Miguel Servet University Hospital, Zaragoza, Spain, Aragon Institute for Health Research (IIS Aragón), University of Zaragoza, Zaragoza, Spain

ORCID logo http://orcid.org/0000-0001-6258-2489

Introduction

Currently OCT is widely used in clinical practice and clinical trials accepting their measurements for the evaluation of the response to treatment and the progression of pathologies [1]. Retinal thickness or central macular thickness (CMT) measured with OCT is particularly used, which correlates with pathological changes and response to treatment for a variety of eye diseases [2].

Currently we still accept the thickness of the retina as the space between surfaces detected, but retinal image segmentation is challenging; structures such as vascular structures, macula, and microaneurysms have low contrast with their background. In contrast, other structures have high contrast with background tissues, but they are difficult to distinguish using classical segmentation techniques [1, 3].

Also, recent works have shown that the most commonly used algorithms in daily practice for retinal layer segmentation have a large number of segmentation errors, especially in the case of age-related macular degeneration, great disruptive pathology such as subretinal fluid, intraretinal cysts and retinal detachments that interrupt the structured logical organization of the retinal layers [4, 5].

Currently we can classify the existing segmentation algorithms into two clusters, mathematical modeling and machine learning approaches. Mathematical modeling is based on the previous anatomical, structural and clinical knowledge that is known about the retina. However pure machine learning algorithms for retinal layer segmentation classifies each pixel from an image on how they fall under a particular layer or boundary, that means that boundaries between layers are not linear [2].

Swept-source (SS)-OCT offers potential advantages due to a modified Spectral-domain (SD) and depth resolved technology which includes an improved imaging range, minor sensitivity roll-off with imaging depth, greater detection efficiencies, and an adjustability to longer imaging wavelengths of 1050nm, this allows a greater choroidal penetration and higher speed for image acquisition. The main difference of SS-OCT is that captures the interferences of the backscattered light from the retina thanks to a wavelength sweeping laser light source and a photodiode detector, in contrast to SD-OCT where a line scan camera and a spectrometer record the interferences between a broadband light source. Those improvements on SS-OCT enables higher density raster scan protocols and deeper image penetration, as a result on en-face reconstructions a better visualization of choroidal detail is possible [6].

Swept-source (SS)-OCT with a wavelength of 1,050 nm and 100,000 A-scans/sec has allowed in-depth visualization of the eye from the retina to the sclera even in patients with moderate to severe cataracts, as well as during eye blinking and/or ocular movement. SS-OCT systems have the potential for superior and simultaneous imaging of the retina and choroid because of the longer wavelength, potentially higher detection efficiency, and lower dispersion [7].

Choroidal thinning has been considered more and more important in the last months, because it has been associated with some ophthalmological pathologies such as age-related macular degeneration (AMD), and also in neurodegenerative diseases such as multiple sclerosis [8], Parkinson’s disease [9] or Fibromyalgia [10] and systemic conditions like diabetes mellitus [11] or pathologies with unclear physiopathology such as migraine [12, 13].

Also, age has been found to be negatively correlated with central choroidal thickness and with central choroidal volume [14], choroidal thickness and volume are also negatively statistically significant concerning the refractive error, and axial length measured with low-coherence reflectometry was also found to be negatively correlated with choroidal thickness and volume. On the other hand, sex has not been found to influence choroidal thickness significantly [15].

Imaging of the choroid was dramatically improved with the development of spectral domain optical coherence tomography (SD-OCT) and was further augmented with the advantage of enhanced depth imaging SD-OCT (EDI SD-OCT) by Spaide and colleagues [16].

Even though, standard cross-sectional (Bscan) OCT imaging are still limited, consequently the choroidal assessment is not as detailed as it could be. However, en-face OCT imaging, is able to provide a high-definition three-dimensional and depth-resolved reconstruction of the choroid, revealing choroidal vascular details not easily visible on cross-sectional OCT imaging. Despite this, SD-OCT is the gold standard for clinical assessment and management of chorioretinal disorders, nevertheless the limited depth of penetration (~850nm), could compromise the choroidal assessment although the selection of the enhanced depth imaging (EDI) method [7].

In the absence of automated segmentation software for SS OCT systems, previous investigators have used manual (mostly single-point) measurement techniques using in-built calibers or modification of retinal segmentation lines to evaluate choroidal thickness; given the high anatomic variability of the choroid, these are impractical for clinical use, are highly dependent on location of measurement, and may be subject to further operator effects [17]. Nevertheless, SS-OCT provides objective and automated measurements of the choroid.

Materials and methods

Between 2015 and 2018, 480 healthy withe caucasian patients were recruited for this cross-sectional study. The inclusion criteria were age between 20 and 79 years, refractive error less than ± 5 diopters (D), axial length between 21 and 25 mm, intraocular pressure less than 21mmHg. The exclusion criteria were concomitant ocular disease (such as glaucoma or retinal pathology); systemic pathologies that could impair the visual system; ocular trauma; laser therapy; and pathology affecting the optic nerve and retina (such as glaucoma, optic neuritis, macular degeneration). We excluded eyes with morphometric parameters of optic disc suggestive of subclinical chronic glaucoma (cup to disc ration≥0.5).

All participants provided written informed consent to participate in this study. The written informed consent for participants and for the study protocol was approved by the Ethics Committee of Clinic Research in Aragon (CEICA) and by the Ethics Committee of Miguel Servet University Hospital, in Zaragoza, Spain, which specifically approved the study procedures. This study was conducted in accordance with the guidelines established by the principles of the Declaration of Helsinki.

Both eyes of each subject were evaluated, but only one of the eyes, randomly selected, was included in the statistical analysis to avoid potential bias by interrelation between eyes of the same subject. In the cases of subjects with exclusion criteria in only one eye, the other eye was selected for the analysis.

A total of 13 eyes were excluded because of not enough OCT quality or exclusion criteria detected during the exploration (morphometric parameters of optic disc compatibles with glaucoma, epiretinal membrane or macular hole). Finally, we included 467 eyes (233 males, 234 females), that were classified in four groups in steps of 15 years, as the difference between the youngest and the oldest is 60 years: Group 1 (composed by 71 eyes of 71 healthy subjects aged between 20 and 34 years); Group 2 (composed by 65 eyes of 65 healthy individual aged 35 to 49 years), Group 3 (with 230 eyes from 230 healthy controls aged between 50 and 64 years) and Group 4 (composed by 97 eyes from 97 healthy subjects aged 65–79 years).

OCT evaluation

Retinal measurements were obtained in all subjects using the DRI Triton SS-OCT device (Topcon, Tokyo, Japan). We performed the 3D+5LineCross protocol (12.0*9.0mm + 9.00mm overlap 8), which allows 100.000 A scans/sec. The DRI Triton SS-OCT provides a quality scale in the image to indicate the signal strength. The quality score ranges from 0 (poor quality) to 100 (excellent quality). Only images with a score >55 were analyzed in our study, and poor-quality images prior to data analysis were rejected.

We exported data of 5 layers using 3DH_DISC (Disc5.2x5.2-Superpixelgrid-200) and the thickness of the 26x26 grid were analyzed to get mean and standard deviation in the four age groups as shown in Fig 1.

[Figure omitted. See PDF.]

Fig 1.

A: Location of 3D+5line cross OCT scans on retina. B: 26*26 grid centered on optic disc, the center area is not shown because the OCT thickness is always zero. C: cross-sectional OCT image and segmentation boundaries.

https://doi.org/10.1371/journal.pone.0240441.g001

Using the Data Collector software of the Triton OCT, measurements of full layers and of 5 different layers were obtained (Fig 2): Retinal thickness (from the inner limiting membrane–ILM- to the retinal pigment epithelium boundaries), Retinal nerve fiber layer (RNFL) (between the ILM to the GCL boundaries), Ganglion cell layer (GCL) + (between RNFL to the inner nuclear layer boundaries) and GCL++ (between ILM to the inner nuclear layer boundaries), and choroid (from the Bruch membrane to the choroidal-scleral interface).

[Figure omitted. See PDF.]

Fig 2. Representation of the five layers measured by Triton optical coherence tomography.

Retina (from the inner limiting membrane–ILM- to the retinal pigment epithelium boundaries), Retinal nerve fiber layer (RNFL) (between the ILM to the GCL boundaries), Ganglion cell layer (GCL) + (between RNFL to the inner nuclear layer boundaries), GCL++ (between ILM to the inner nuclear layer boundaries), and choroid (from the Bruch membrane to the choroidal-scleral interface).

https://doi.org/10.1371/journal.pone.0240441.g002

Automated built-in calibration software, Topcon Advances Boundary Software (TABS) determined the distance between the delimiting lines in retina and choroidal plexus, establishing 7 boundaries and five layers.

Data analysis

Comparison between age groups was performed using analysis of variance (ANOVA) for each thickness measurement for four quadrants and four halves. Post-Hoc analysis was performed to obtain statistical differences in each comparison between the four age groups (Group 1 compared with Group 2; Group 1 compared with Group 3; Group 1 compared with Group 4; Group 2 compared with Group 3; Group 2 compared with Group 4; and Group 3 compared with Group 4). Correction for multiple comparisons was used in all analyses. In order to graphically see the evolution of the thickness of each segmented layer, a numerical method called the finite element method (FEM) was used. Numerical techniques such as the finite element method have been extensively implemented as effective and noninvasive methods to analyze biological tissues, and in particular in ophthalmology [18–20]. The finite element method is a computational tool which allows analysing the stress/strain behaviour of a structure subjected to different loads and boundary conditions. Here, this method has only been used to study the evolution of retinal layers thickness along time and to graphically see the evolution of each of them by a contour plot map. Thus, a finite element mesh using ABAQUS (Abaqus 6.14, Simulia, Rhode Island, USA) mimicking the Triton grid (Disc5.2x5.2-Superpixelgrid-200) was developed. The finite element mesh was constructed using membrane quadrilateral elements and the size of the mesh was 5.2x5.2mm. Afterwards, the thickness of each box of the grid for each group and each layer was introduced. In this way, the evolution of the thickness of each layer can be seen from a spatio-temporal point of view from Group 1 to Group 4.

Results

We analysed a total of 467 eyes from 467 healthy subjects, 71 eyes from 71 individuals between 20–34 years (group 1), 65 eyes from 65 subjects between 35–49 years (group 2), 230 eyes from 230 subjects between 50–64 years (group 3) and 97 eyes from 97 individuals between 65–79 years (group 4).

Mean and standard deviation for each layer were represented in Fig 3 for the four age groups.

[Figure omitted. See PDF.]

Fig 3. Representation of mean and standard deviation for each layer measured by Triton optical coherence tomography, for the four age groups.

Group 1 (composed by 71 healthy subjects with age between 20 and 34 years); Group 2 (composed by 65 individuals with 35–49 years), Group 3 (with 230 healthy controls with 50–64 years) and Group 4 (composed by 97 healthy subjects with 65–79 years). On the left column, the mean and standard deviation is shown for each layer and age group. On the right column, the quartiles are plotted to show data dispersion and the presence of outliers.

https://doi.org/10.1371/journal.pone.0240441.g003

Fig 4 represents the contour map obtained by finite element modeling of the evolution of the thickness of each layer for the different age groups, which shows progressive thinning of all layers with age, especially from 50 years onwards.

[Figure omitted. See PDF.]

Fig 4. Finite element contour map of the evolution of the thickness of each layer measured by Triton optical coherence tomography for the four age groups (20–34 years, 35–49 years, 50–64 years, and 65–80 years).

On the left, the contour bar represents the equivalence between colour and thickness value (in μm) for each layer. Blue zones correspond to thicker areas and red zones correspond to thinner ones. The optic nerve is represented in grey colour. Different scales have been used to represent the thickness value for each layer. It can be seen a progressive thinning of all layers with age, especially from 50 years onwards.

https://doi.org/10.1371/journal.pone.0240441.g004

A comparative analysis between age groups were performed for all layers to find thickness differences for four quadrants (temporal superior -TS-; temporal inferior -TI-; nasal superior -NS-; and nasal inferior -NI-) and four halves (superior, inferior, nasal and temporal).

Retinal full layer analysis

There are highly significant differences between groups in our analysis for each quadrant and halves. In the post-hoc analysis, a significant thickening that did not overcome the Bonferroni correction for multiple comparison was found for TI quadrant at group 2 compared with group 1 (p = 0.038) and also at temporal half (p = 0.043). Considerable differences that overcame the Bonferroni correction or multiple comparison were found comparing group 2 vs group 3, in which a thinning trend is visible for TS (p = 0.004), and TI (p = 0.007) quadrants and temporal (p = 0.005) halves. Comparing group 3 vs group 4, significant thinning was found for the four quadrants (p≤0.002), and the four halves (p<0.001). Comparing group 1 vs group 3 we did not find any statistical differences between groups although there is a presumed thinning in group 3. On the other hand, there is an obvious thinning in group 4 compared with group 1 and group 2 for all quadrants and halves (p<0.001).

Retinal nerve fiber layer analysis

Significant differences for all quadrants and halves were found between age groups, except for NS quadrant and nasal half. A significant RNFL thickening that overcame Bonferroni correction was observed at group 2 compared with group 1 in TS quadrant (p = 0.001), inferior half (p<0.001) and temporal half (p = 0.001). Group 3 presents a significant thickening at inferior half (p = 0.002) when compared with group 1, but a significant thinning at TS quadrant (p = 0.001) when compared with group 2. Group 4 presents just a significant thickening at NI quadrant (p = 0.002) compared with group 1, but compared with group 2, it is observed a significant thinning at TS (p<0.001) and TI (p<0.001) quadrants and at superior (p<0.001) and temporal (p<0.001) halves. Same conduct happens when this group is compared with group 3 at TI quadrant (p<0.001) (Table 1).

[Figure omitted. See PDF.]

Table 1. Mean ± standard deviation of four quadrants and halves for age groups in the retinal nerve fiber layer, GCL+ and GCL++, and comparison of thickness between age groups.

https://doi.org/10.1371/journal.pone.0240441.t001

RNFL to inner nuclear layer (GCL+)

GCL+ appears to be significantly thicker in group 2 when comparing with group 4 at TS quadrant (p = 0.001), inferior and temporal halves (p = 0.001), but not compared with groups 1 and 3. Group 3 present a significant thinner GCL+ layer at NI quadrant and inferior half (p<0.001) when compared with group 1, however GCL+ is significantly thicker in this group compared with older population (group 4) at TS quadrant and superior half (p = 0.001). Finally, group 4 exhibits a significant thinning of GCL+ layer for every quadrants and halves (p<0.001) (Table 1).

From inner limiting membrane to inner nuclear layer (GCL++)

Significant differences between groups in the GCL++ layer were found almost in every quadrant and half, except for NI quadrant and nasal half, that apparently are not predictive of changes in this layer (Table 1). Comparing group 3 vs group 1, significant differences were not found, but compared with group 2 significant thinning appears in TS quadrant (p<0.001). On the other hand, group 4 exhibits a thinner GCL++ layer compared with group 1 in TI quadrant (p = 0.001) and superior and temporal halves (p = 0.002). Compared with group 2 thinning gets more significant in TS and TI quadrants (p<0.001), and superior and temporal halves (p<0.001) (Table 1).

Retina (from RNFL to photoreceptors layer)

There are no significant differences when comparing group 1 with group 2, but when comparing with group 3 significant thinning is found at TS and NS quadrants (p = 0.003 and p = 0.001) and superior and nasal halves (p = 0.001 and p = 0.006). Furthermore, when comparing with group 4, significant thinning is found in all quadrants and halves (p<0.001) except for NI quadrant (p = 0.041) (Table 2).

[Figure omitted. See PDF.]

Table 2. Mean ± standard deviation of four quadrants and halves for age groups in the full layers, retina and choroid measurement, and comparison of thickness between age groups.

https://doi.org/10.1371/journal.pone.0240441.t002

By contrast comparing group 2 with group 3, there seems to be a significant thinning at TS (p<0.001), NS (p = 0.013) and NI (p = 0.009) quadrants and all halves (p<0.001, p = 0.024, p = 0.008, p = 0.003), but when comparing with group 4 significant thinning occurs in all quadrants (p<0.001, p<0.001, p<0.001 and p = 0.008) and all halves (p<0.001). Finally, when comparing group 3 with group 4 significant thinning is found at TS (p = 0.008) and TI (p<0.001) quadrants and superior (p = 0.016) and temporal (p<0.001) halves.

Choroid analysis

Choroid howed significant differences between age groups for all quadrants and halves (Table 2). At second age group (35 to 49 years) it seems to exist a choroidal thickening at superior half (p = 0.001) when compared with group 1. By contrast, group 3 presents a thinner choroidal layer compared with previous groups, but this thinning is just significant when compared with group 2 for TS and TI quadrants, and for temporal half (p<0.001). The thinnest choroidal values can be seen for group 4, which is significant when comparing with group 1, 2 of 3 for every parameter analyzed (p<0.001).

Except the GCL+ with a linear thinning ratio of -0.050 μm/year, the tendency in the other layers differs. A thickening ratio of 0.476 μm/year for the full layers complex, 0.010 μm/year for the retina, 0.856 μm/year for the choroid, 0.125 μm/years for the RNFL and 0.087 μm/years for the GCL++ is observed until de third decade of life. Subsequently a thinning ratio of -1.406 μm/year for the full layers complex, -0.058 μm/year for the retina, -1.364 μm/year for the choroid, -0.084 μm/years for the RNFL and -0.131 μm/years for the GCL++ is observed from de third decade of life.

Discussion

Previous histopathological studies of the choroid have shown that choroidal thickness decreases with age; this finding has been confirmed by posterior in-vivo studies [7, 14, 17, 21]. Choroidal thinning on OCT has been associated with ophthalmological conditions such as age-related macular degeneration, age-related choroidal atrophy, axial length, or systemic conditions such as diabetes, Alzheimer’s disease, cognitive impairment, vascular diseases or obesity [11, 21, 22].

In addition, the loss of vascular perfusion related to age has been described in previous studies performed with OCT-Angiography [23], even this alteration of blood flow could be observed by magnetic resonance [24]. Furthermore, a significant reduction of retrobulbar circulation related to age has been found in Doppler studies, this reduction could be also related with the increased prevalence of cardiovascular risks related with age such as hypertension, diabetes, lipid disorders and sedentary lifestyle [25–28].

Traditionally, the age of 40 years has been established as the age of onset of the classic ophthalmological alterations, although other pathologies such as macular alterations are more common beyond 50 years, so these changes in the choroid can be a predictive factor of future clinical alterations [7, 29, 30].

In our study, we have seen that global retinal tendency is toward thickening in the third decade of life (group 2) and towards thinning after the fourth decade (group 3). Regarding the choroid, the same pattern was seen and also for the rest of layers, except at GCL+ layer, where there is a continuous thinning from group 1 to group 4. These results are consistent with those of Hanumunthadu et al. 2018, their age-stratified analysis suggested that choroidal thickness was smaller in children and younger adult population, and it appeared to decrease again in older adults [21].

In contrast with published studies, our analysis has focused on a division by quadrants and halves in the peripapillary area. We observed that superior half and quadrants are thicker than their inferior homologous, however, it seems that there is no difference when comparing nasal and temporal orientation. As we have seen previously, group 1 seems to have a very similar choroidal thickness to group 3 in practically all the quadrants and halves, except in the superior half, where apparently the thickness of the choroid is greater in this group.

Comparing the thinning of the total retina and the one produced specifically in the choroid, we can see that it is in group 4 where the thinning of this last layer is more pronounced, especially in the inferior half of the retina. Same thing happens for groups 1 and 2, but not in 3, where thinning is symmetric when comparing the superior quadrants and halves with the inferior ones and the nasal quadrants and halves with the temporal ones.

The potential clinical application of our findings is to know which thickness values in each layer should be considered pathological. This study may be the base to build a normative database that allows health personnel to improve the interpretation of retinal measurements, even with a color scale similar to that one used by most OCT devices (red color means pathological thinning, yellow color indicated slight thinning, green color means normal, blue indicates slight thickening and purple means abnormal thickening); this helps to determine which layers are out of normal limits in each range of age. This is an important limitation of the Triton current software because we can visualize numerous measurements of all retinal layers, but ophthalmologists do not know if these values are normal for the age of each subject or patient.

Our findings suggest that there is a progressive and physiological thinning of all retinal layers from the third decade of life. In any case, it would be useful to extend this study further with a larger population of ordinary individuals across a range of different ethnicities and from longer geographic areas, to clarify the physio-pathological mechanisms that affect retinal layers.

Supporting information

[Figure omitted. See PDF.]

S1 Data.

https://doi.org/10.1371/journal.pone.0240441.s001

(SAV)

Citation: Viladés E, Pérez-del Palomar A, Cegoñino J, Obis J, Satue M, Orduna E, et al. (2020) Physiological changes in retinal layers thicknesses measured with swept source optical coherence tomography. PLoS ONE 15(10): e0240441. https://doi.org/10.1371/journal.pone.0240441

References

1. Venhuizen FG, van Ginneken B, Liefers B, van Grinsven MJJP, Fauser S Hoyng C, et al. Robust total retina thickness segmentation in optical coherence tomography images using convolutional neural networks. Biomed Opt Express. 2017; 8(7): 3292–3316. pmid:28717568

2. Keane PA, Liakopoulos S, Jivrajka RV, Chang KT, Alasil T, Walsh AC, et al. Evaluation of optical coherence tomography retinal thickness parameters for use in clinical trials for neovascular age-related macular degeneration. Invest. Ophthalmol. Vis. Sci. 2009; 50: 3378–3385. pmid:19264895

3. Almotiri J, Elleithy K, Elleithy A. A multi-anatomical retinal structure segmentation system for automatic eye screening using morphological adaptive fuzzy thresholding. IEEE J Transl Eng Health Med. 2018; 6: 3800123. pmid:29888146

4. Waldstein SM, Gerendas B S, Montuoro A, Simader C, Schmidt-Erfurth U. Quantitative comparison of macular segmentation performance using identical retinal regions across multiple spectral-domain optical coherence tomography instruments. Br. J. Ophthalmol. 2015; 99: 794–800. pmid:25563763

5. Hanumunthadu D, Wang JP, Chen W, Wong EN, Chen Y, Morgan WH, et al. Impact of retinal pigment epithelium pathology on spectral-domain optical coherence tomography-derived macular thickness and volume metrics and their intersession repeatability. Clin. Exp. Ophthalmol. 2017; 45: 270–279. pmid:28052542

6. Ahmad M, Kaszubski PA, Cobbs L, Reynolds H, Smith RT. Choroidal thickness in patients with coronary artery disease. PLoS One. 2017; 12(6): e0175691. pmid:28632734

7. Adhi M, Ferrara D, Mullins RF, Baumal CR, Mohler KJ, Kraus MF, et al. Characterization of choroidal layers in normal aging eyes using enface swept- source optical coherence tomography. PLoS One. 2015; 10(7): e0133080. pmid:26172550

8. Higashiyama T, Nishida Y, Ohji M. Optical coherence tomography angiography in eyes with good visual acuity recovery after treatment for optic neuritis. PLoS One. 2017; 12(2): e0172168. pmid:28192497

9. Garcia-Martin E, Pablo LE, Bambo MP, Alarcia R, Polo V, Larrosa JM, et al. Comparison of peripapillary choroidal thickness between healthy subjects and patients with Parkinson’s disease. PLoS One. 2017; 12(5): e0177163. pmid:28510576

10. Ulusoy MO, Kal A, Işik-Ulusoy S, Kal Ö. Choroidal thickness in patients with fibromyalgia and correlation with disease severity. Indian J Ophthalmol. 2018; 66(3): 428–432. pmid:29480257

11. Abadia B, Suñen I, Calvo P, Bartol F, Verdes G, Ferreras A. Choroidal thickness measured using swept-source optical coherence tomography is reduced in patients with type 2 diabetes. PLoS One. 2018; 13(2): e0191977. pmid:29394291

12. Ekinci M, Ceylan E, Cagatay HH, Keles S, Hüseyinoglu N, Tanyildiz B, et al. Retinal nerve fibre layer, ganglion cell layer and choroid thinning in migraine with aura. BMC Ophthalmol 2014; 14: 75. pmid:24885597

13. Colak HN, Kantarcı FA, Tatar MG, Eryilmaz M, Uslu H, Goker H, et al. Retinal nerve fiber layer, ganglion cell complex, and choroidal thicknesses in migraine. Arq Bras Oftalmol 2016; 79(2): 78–81. pmid:27224067

14. Akhtar Z, Rishi P, Srikanth R, Rishi E, Bhende M, Raman R. Choroidal thickness in normal Indian subjects using Swept source optical coherence tomography. PLoS One. 2018; 13(5): e0197457. pmid:29768485

15. Michalewski J, Michalewska Z, Nawrocka Z, Bednarski M, Nawrocki J. Correlation of choroidal thickness and volume measurements with axial length and age using swept source optical coherence tomography and optical low-coherence reflectometry. Biomed Res Int. 2014; 2014: 639160. pmid:25013793

16. Mazzaferri J, Beaton L, Hounye G, Sayah DN, Costantino S. Open-source algorithm for automatic choroid segmentation of OCT volume reconstructions. Sci Rep. 2017; 7: 42112. pmid:28181546

17. Wakatsuki Y, Shinojima A, Kawamura A, Yuzawa M. Correlation of aging and segmental choroidal thickness measurement using swept source optical coherence tomography in healthy eyes. PLoS One. 2015; 10(12): e0144156. pmid:26632821

18. Pandolfi A, Holzapfel GA. Three-dimensional modeling and computational analysis of the human cornea considering distributed collagen fibril orientations. J Biomech Eng. 2008; 130(6): 061006. pmid:19045535

19. Shin A, Yoo L, Park J, Demer JL. Finite element biomechanics of optic nerve sheath traction in adduction. J Biomech Eng. 2017; 139(10).

20. Schutte S, van den Bedem SP, van Keulen F, van der Helm FC, Simonsz H J. A finite-element analysis model of orbital biomechanics. Vision Res. 2006; 46(11): 1724–31. pmid:16413594

21. Hanumunthadu D, Ruiz-Medrano J, Dumpala S, Jabeen A, Jabeen A, Goud A, et al. Comparison of choroidal vessel thickness in children and adult eyes by enhanced-depth imaging optical coherence tomography imaging. Int J Ophthalmol. 2018; 11(4): 681–686. pmid:29675391

22. Zhao J, Wang YX, Zhang Q, Wei WB, Xu L, Jonas JB. Macular choroidal small-vessel layer, Sattler's layer and Haller's layer thicknesses: The Beijing Eye Study. Sci Rep. 2018; 8(1): 4411. pmid:29535365

23. Falavarjani KG, Shenazandi H, Naseri D, Anvari P, Kazemi P, et al. Foveal avascular zone and vessel density in healthy subjects: an optical coherence tomography angiography study.J Ophthalmic Vis Res. 2018; 13(3): 260–265. pmid:30090182

24. Emeterio Nateras OS, Harrison JM, Muir ER, Zhang Y, Peng Q, Chalfin S, et al. Choroidal blood flow decreases with age: an MRI study. Curr Eye Res. 2014; 39(10): 1059–1067. pmid:24655028

25. Modrzejewska M, Siesky B, Amireskandari A, Holland S, Grzesiak W, Zaborski D et al. Parameters characterizing age-dependent retrobulbar circulation in healthy subjects measured by color Doppler ultrasonography. Curr Eye Res 2015; 40:729–736. pmid:25197962

26. Steigerwalt RD Jr., Belcaro GV, Laurora G, Cesarone MR, De Sanctis MT, Incandela L. Ocular and orbital blood flow in patients with essential hypertension treated with trandolapril. Retina Phila Pa 1998; 18:539–545.

27. Foroozan R, Savino PJ, Sergott RC. Embolic central retinal artery occlusion detected by orbital color Doppler imaging. Ophthalmology 202; 109:744–747. pmid:11927432

28. Costa VP, Kuzniec S, Molnar LJ, Cerri GG, Puech-Leão P, Carvalho CA. Clinical findings and hemodynamic changes associated with severe occlusive carotid artery disease. Ophthalmology 1997; 104:1994–2002. pmid:9400757

29. Nikolopoulou E, Lorusso M, Micelli Ferrari L, Cicinelli MV, Bandello F, Querques G, et al. Optical coherence tomography angiography versus dye angiography in age-related macular degeneration: sensitivity and specificity analysis. Biomed Res Int. 2018; 2018: 6724818. pmid:29707575

30. Chalam KV, Sambhav K. Optical coherence tomography angiography in retinal diseases. J Ophthalmic Vis Res. 2016; 11(1): 84–92. pmid:27195091