Abstract

The lamina cribrosa (LC) is an active structure that responds to the strain by changing its morphology. Abnormal changes in LC morphology are usually associated with, and indicative of, certain pathologies such as glaucoma, intraocular hypertension, and myopia. Recent developments in optical coherence tomography (OCT) have enabled detailed in vivo studies about the architectural characteristics of the LC. Structural characteristics of the LC have been widely explored in glaucoma management. However, information about which LC biomarkers could be useful for the diagnosis, and follow-up, of other diseases besides glaucoma is scarce. Hence, this literature review aims to summarize the role of the LC in nonophthalmic and ophthalmic diseases other than glaucoma. PubMed was used to perform a systematic review on the LC features that can be extracted from OCT images. All imaging features are presented and discussed in terms of their importance and applicability in clinical practice. A total of 56 studies were included in this review. Overall, LC depth (LCD) and thickness (LCT) have been the most studied features, appearing in 75% and 45% of the included studies, respectively. These biomarkers were followed by the prelaminar tissue thickness (21%), LC curvature index (5.4%), LC global shape index (3.6%), LC defects (3.6%), and LC strains/deformations (1.8%). Overall, the disease groups showed a thinner LC (smaller LCT) and a deeper ONH cup (larger LCD), with some exceptions. A large variability between approaches used to compute LC biomarkers has been observed, highlighting the importance of having automated and standardized methodologies in LC analysis. Moreover, further studies are needed to identify the pathologies where LC features have a diagnostic and/or prognostic value.

1. Introduction

The lamina cribrosa (LC) is a mesh-like structure localized in the posterior scleral canal of the optic nerve head (ONH), allowing retinal ganglion cell (RGC) axons to pass through to the brain. It is a fenestrated complex that also accommodates vessels that nourish the retina. A large circumpapillary ring of collagen and elastin fibers, in the immediate peripapillary sclera, protects the LC against the mechanical strain, such as that induced by an imbalance between intraocular pressure (IOP) and intracranial pressure (ICP) [1, 2]. Due to its anatomical location, between two differently pressurized compartments, there is a pressure gradient along the LC, denominated translaminar pressure difference (TLPD), which can be calculated as the difference between IOP and ICP in the subarachnoid space (SAS) [3, 4]. Despite being an extremely relevant structure to the eye’s anatomy and function, little is known about the LC. LC morphology plays an important role in the development and progression of ophthalmic pathologies, notably on glaucomatous optic neuropathy, intraocular hypertension, and myopia [58]. The structural deformation and the correlated compression across the LC lead to blockade of axonal transport and eventually RGC death [9].

Recent advances in in vivo medical imaging techniques, such as optical coherence tomography (OCT), have allowed the visualization of deep connective tissues, including the LC, in greater detail (Figure 1) [10, 11].

Specific developments in OCT software, such as enhanced depth imaging (EDI), and light-attenuation correction software such as adaptive compensation (AC) significantly improved the visibility of the LC without compromising acquisition time. EDI-OCT was originally developed in order to improve the visualization of the choroid, although it has also been adopted to improve cross-sectional images of the LC. AC is a postprocessing technique developed to remove blood vessel shadows and enhance tissue contrast in order to facilitate posterior LC surface detection [12, 13]. In addition to these software developments, several studies have shown that swept-source OCT (SS-OCT) further improves the visualization of the LC [14, 15].

While an increasing number of works have studied the relevance of the LC (and its changes) in glaucoma [16], data on other diseases are still scarce. Hence, this review intends to provide a broader vision and a better comprehension of the measurable laminar structural features that have been identified as relevant for nonglaucomatous pathologies in the published literature.

2. Methods

2.1. Study Selection

A literature search was conducted in the MEDLINE (PubMed) bibliographic database on 15th May 2020. The search query was (optical coherence tomography NOT angiography) AND (lamina OR cribrosa). Only articles published in English were considered, and no publication date restriction was added. The exclusion criteria were (i) only included glaucomatous eyes in the experimental group; (ii) not conducted in humans; (iii) review articles or case reports; (iv) exclusive focus on imaging techniques and not presenting clinical data; (v) no evaluation of the lamina cribrosa; and (vi) no mention to LC structural parameters and how they were measured/extracted. This led to a total of 408 references, which were narrowed down to 56 after title/abstract screening, followed by a full-text screening (see Figure 2). The 56 included studies provided quantitative values for each of the analyzed features and described how the quantification was performed.

2.2. Data Collection

In this review, our main aim was to identify potential biomarkers in the morphology of the LC that were associated with, and indicative of, certain pathologies. Therefore, we have opted to only report those that performed a statistical comparison between an experimental and a control group. The extracted data for each paper consist of the LC structural parameters, their mean and standard deviation (SD), and the values of the statistical analysis performed between experimental and control groups. Moreover, the image processing methodology applied to compute each feature was taken into account for posterior comparisons.

For all the included articles, the following characteristics are obtained and presented in Table 1: sample size, including the number of patients and eyes per group; age and statistical comparison ( value) between control and experimental groups; OCT device model, manufacturer, and light-source wavelength; cutoff value for the signal strength index (SSI) or similar qualitative image criteria used to exclude patients/eyes; and field of view of the OCT image [17]. Moreover, the procedures followed to measure the LC features, and their respective values, were also collected and compiled.

Data collection comprehended all structural components related to the LC and the surrounding ONH region that were included on the OCT B-scan images. Several locations, planes (superior, middle, and inferior), and sectors were considered for the measurements. The sectors were defined according to the Garway-Heath map [71]. The approach of data extraction by one investigator (ASP) with further verification by a senior author (JBB) has been used, as this has been demonstrated to be as accurate as double independent data extraction [72].

2.3. Data Analysis

The obtained data were used to calculate the frequency of each LC structural feature in the published literature and to determine the mean values of the most frequently reported features. Statistical relevance, given by the value, was also taken as a complement for study results and to comprehend the relation and differences found between study groups. To average the data for each LC parameter, pooled mean and pooled SD were determined according to equations (1) and (2), respectively, where represents the number of eyes included in the study, the mean value, the standard deviation, and the number of studies analyzed:

3. Results

All structural LC features were analyzed, and the studies were organized in three groups: healthy group (n = 19), nonophthalmic disease group (n = 6), and ophthalmic (nonglaucomatous) disease group (n = 31). Overall, LC depth (LCD) and LC thickness (LCT) have been the most studied features, appearing in 75% and 44.6% from the total articles. Other features, such as prelaminar tissue thickness (PTT) (21.4%) studied in ophthalmic and nonophthalmic diseases, and lamina cribrosa curvature index (LCCI) (5.4%), LC global shape index (3.6%), LC defects (3.6%), slope of the LC (3.6%), distance between the inner surface of the LC (3.6%), SAS (3.6%), and LC strains/deformations (1.8%) studied in ophthalmic diseases only, are also referenced but in fewer studies.

One bar chart, summarizing the most studied features in the ophthalmic and nonophthalmic disease groups, is presented in Figure 3. Hence, in the following sections, a detailed analysis on how these two biomarkers have been measured is provided. Moreover, a detailed explanation on how LCD and LCT measurements were carried out in each study is presented in Table 2. The normative values for the groups are also presented and discussed. In cases where more than one measurement was performed for the same feature (e.g., in different planes (superior, middle, or inferior), different scan directions (vertical and horizontal), or 2 eyes (left and right)), the pooled mean and SD were determined according to equations (1) and (2).

Only diseases for which at least 20 eyes were included in the studied group (cumulatively over all the evaluated papers) were considered for the average calculations (equations (1) and (2)).

3.1. LCT Measurements

LCT has been defined in the literature as the distance between the anterior and posterior borders of the highly reflective region visible below the optic disc cup in B-scan cross sections of the ONH (see the red arrow in Figure 4) [18]. However, a discrepancy between the LCT measurements, namely, between the locations used to calculate the LCT average, has been observed and reported in Table 2. For example, Lee et al. [74] considered three locations in each eye (midhorizontal, superior, and inferior midperipheral), with a separation of 100 m between the points. Bartolomé et al. [19] determined the points as close as possible to the vertical center of the ONH, which was identified as the point where the trunk of central retinal vessels extends from the ONH, as reported by Park et al. [12]. Other authors, such as Xiao et al. [20], considered LCT as the average of the central and paracentral points (150 m from the center point in the horizontal and vertical directions).

3.2. LCD Measurements

In the literature, LCD (also named as anterior lamina cribrosa depth in several studies) is defined as the perpendicular distance from the BMO plane to the maximum depth point of the anterior LC surface. All articles included in this review provide the measurements relative to Bruch’s membrane opening (BMO). Only the measurements relative to BMO were considered for the calculations since all articles provide the measurement relative to this plane. Two studies, Rhodes et al. [21] and Luo et al. [22], also considered the scleral plane and the ASCO as the reference for depth measurements. These differences have been shown to lead to measurement bias, as reported by Luo et al. [22], who obtained 402  91 m, 309  88 m, and 332  96 m for BMO, ASCO, and scleral reference planes, respectively. The number of selected points and B-scan planes to average the measure also influence the precision of the results. Other authors, such as Park et al. [23], obtained measurements as the average from 11 equidistant planes that divided the optic disc diameter into 12 equal parts vertically in each eye. A line was drawn from each of the two LC insertion points perpendicularly to the line connecting the two Bruch’s membrane edges (see line D in Figure 4). The area surrounded by these two lines was measured (see area S in Figure 4). The mean LCD was approximated by dividing S by the length of D for each of the 11 horizontal OCT scans. Finally, Lee et al. [24] defined LCD as the mean of 3 values obtained from the 3 upper B-scans ( to the scan), the 3 central B-scans ( to the scan), and at the 3 lower B-scans (from the to the scan) passing through the ONH. Commonly, temporal adjacent points were selected because the maximally depressed point was often close to the central vessel trunk, and its shadow obscured the LC [25].

3.3. Features’ Applicability and Measurements

This section details the mean values for the two dominating LC structural features (LCT and LCD) in the three groups (healthy controls, ophthalmic, and nonophthalmic diseases). The values were calculated based on the articles presented in Table 2 for each group and disease, and the mean and SD values for each group are summarized in Figure 5.

3.3.1. Healthy Group Measurements

Analysis of healthy subjects is very important to establish normative values for the healthy population, and hence facilitate the diagnosis and follow-up of the pathology. The studies that included only healthy subjects, as well as those comparing patients to a healthy control group, were selected, and the LCT and LCD average values were determined. The observed averages were 261  39 m (range: 211–323 m) and 386  91 m (range: 293–441 m) for the LCT and LCD, respectively. Figures 5(a) and 5(b) show a comparison between the three groups for both features. These parameters seem to be influenced by several factors, such as age and racial ancestry. Rhodes et al. [21] conducted a study in healthy eyes and concluded that the LC was significantly anteriorly displaced with increasing age in those with European ancestry in contrast to those with African ancestry.

3.3.2. Ophthalmic Disease Group Measurements

The ophthalmic disease group represented the largest group (n = 31) and included a large number of conditions, the most common being myopia, retinal vein occlusion (RVO), nonarteritic anterior ischaemic optic neuropathy (NAION), pseudoexfoliation syndrome (PXS), superior segmental optic nerve hypoplasia (SSOH), compressive optic neuropathy (CON), age-related macular degeneration (AMD), autosomal dominant optic atrophy (ADOA), and diabetic macular edema (DME). For ophthalmic patients, mean LCT and LCD were 211  33 m and 403  90 m, respectively (Figure 5). The graphics in Figure 6 presents the mean values for different ophthalmic diseases in comparison with the healthy population (horizontal dashed green line). Regarding the LCT, its mean was lower for every pathology in this group except for SSOH (Figure 6(a)). Overall, the studies that reported LCD in nonglaucomatous ophthalmic diseases showed a slightly higher mean LCD compared with healthy controls (Figure 5(b)). Nonetheless, this trend is not significant because the standard deviations of all the diseases cross the standard deviation of the healthy LCD (Figure 6(b)). For example, Rebolleda et al. [26] found a lower average LCD in the superior, middle, and inferior planes in eyes with NAION compared to healthy eyes (), which was also true between unaffected fellow eyes when compared to healthy eyes (). The maximum value for LCD was found in cases of Graves’ orbitopathy with proptosis and/or compressive optic neuropathy. Seo et al. [27] conducted a study in these patients and reported 462.79  95.96 m and 621.39  78.39 m values, at baseline, for the muscle-dominant and fat-dominant group, respectively.

3.3.3. Nonophthalmic Disease Group Measurements

The number of studies in this group was smaller than in the other groups. The registered diseases were diabetes mellitus [28], Parkinson’s disease (PD) [29], obstructive sleep apnea syndrome (OSAS) [30], Alzheimer’s disease (AD) [31, 32], mild cognitive impairment (MCI) [31], and migraine [33]. The mean LCT and LCD were 234  36 m and 390  68 m, respectively, as shown in Figure 5. The graphics in Figure 7 presents LCT and LCD for each pathology in comparison to the healthy population (horizontal dashed green line).

LCT measurements seem to be lower relative to the healthy group, with the exception of diabetes mellitus. Akkaya et al. [28] described a significantly higher mean LCT in diabetic patients when compared to a healthy group, 271.61  33.96 m vs. 248.50  5.40 m, respectively (). Regarding LCD, diabetes mellitus showed significantly lower mean values in comparison to healthy controls in one study [28] (351  59 m vs. 420  90 m; ). The maximum mean LCD absolute value (deepest ONH cup) was described in patients with migraine (see Figure 7(b)). Sirakaya et al. [33] reported significantly higher mean LCD values for both migraine groups (412.15  58.80 m with aura and 405.57  55.39 m without aura) when compared to the healthy group (355.34  65.53 m; ).

4. Discussion

The present study highlights which LC structural parameters have been analyzed in the literature with a focus on nonglaucomatous diseases. Overall, the most commonly studied parameters were LCD and LCT. The disease groups (ophthalmic and nonophthalmic) presented lower values for mean LCT, relative to the healthy population (Figure 5(a)). In parallel, mean LCD values were higher (deeper ONH cup) for these groups (Figure 5(b)). An exception in the nonophthalmic disease group was DM, which presented a shallower cup and thicker LC, when compared to healthy subjects. Akkaya et al. proposed that this evidence supports the “neuroprotective effect of DM on glaucomatous optic neuropathy and suggests that LCT and lamina cribrosa position mediate this protective effect.” [28]

This study shows that LC structural features are significantly different between healthy patients and some (nonglaucomatous) ocular and systemic pathologies. As such, there is a potential to add them as additional clinical features for clinical diagnosis. Nonetheless, being patient-specific features, LC features might hold an even better role for patient follow-up, signalling disease status’ change. Unfortunately, we did not find any longitudinal studies focusing on this matter in this review. This fact highlights the need for longitudinal studies linking LC parameters and diseases, similarly to what is now common in glaucoma-related studies [75].

LC features are influenced by factors such as age, race, and also by the way measurements are carried out. However, these factors were not used as segmentation criteria in this study due to the lack of this information in several studies. For future works in this field, it is important to take into account these factors when analysing and comparing results between studies since they are a potential source of bias. Moreover, the current methods are heterogeneous (see Table 2), which may lead to imprecise comparisons between studies. For depth measurements, the consensus is to use the BMO plane as a reference, but the way the feature is measured is not consistent among research groups. One of the causes for this heterogeneity is the fact that the analysis of LC features still requires a considerable amount of manual input. This causes measurement bias due to the inherent difficulty of the manual delineation of the structure. This lack of automation increases the likelihood that each research group adopts their own reference points and methods. Besides, studies usually report averages of a limited number of B-scans without capturing the whole LC. The distance between these B-scan slices, as well as their number and position, may also be a source of discrepancies when comparing studies. Finally, some authors have pointed out the fact that the BMO reference place might be biased due to choroidal thickness changes and that perhaps the anterior sclera reference plane would be better suited for these calculations [76, 77]. Ideally, similar measurement methods should be adopted across all research groups. Currently, the measurement of the LC features is laborious and time consuming. As such, automation might hold the key to reduce bias in LC feature measurement. There is a need for easy-to-use software that can automatically measure LC features (possibly starting with LCT and LCD), ideally capturing all the information from OCT volumes, instead of selecting some of the B-scans. Providing such a tool with a fast and repeatable computation would contribute to making LC features a part of everyday clinical practice.

The main limitation of this review is the reduced number of studies, mainly in the nonophthalmic disease group, which precludes definite conclusions. Moreover, due to the lack of individual study data, it was not possible to perform statistical comparisons between groups and pathologies. As such, our results point towards differences that need to be better clarified. Nonetheless, LC features’ ability to discriminate between these groups is supported by results presented by several individual studies, as reported in Table 2. Lastly, it is noteworthy to mention that the statistical analysis performed on groups (ophthalmic and nonophthalmic diseases) may be biased by different pathologies comprised in each group.

5. Conclusion

There is a growing interest in LC features outside the glaucoma field. The results of this meta-analysis show several promising features (mainly, LCT and LCD) that may be relevant for clinical practice. Nevertheless, further studies are needed to validate these findings, and longitudinal data are needed to clarify the potential for use in patient follow-up. Moreover, efforts should be employed to develop automated tools that can capture LC features from OCT data in a standardized manner, thus allowing more accurate comparisons between studies. These efforts should enable to further explore the potential of LC parameters for use in daily clinical practice.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Horizon 2020 Research and Innovation Programme (grant agreement no. 780989: Multi-modal, multi-scale retinal imaging project) and was funded by Portuguese National Funds through the FCT, Fundação Para a Ciência e a Tecnologia, I.P., in the scope of the project UIDB/04559/2020.