Evaluation of Corneal Tomographic, Densitometric and Aberrometric Features Following Accelerated Corneal Cross-Linking in progressive keratoconus

DOI: https://doi.org/10.21203/rs.3.rs-1692736/v1

Abstract

Background: To evaluate the 12 months changes in tomographic, densitometric, and aberrometric parameters in keratoconic eyes after accelerated corneal cross-linking (CCL) and classify densitometric course in different stages of the keratoconus separately (Mild, moderate and severe).

Methods: In a prospective observational study, 67 keratoconic eyes of 67 patients that underwent accelerated epithelium-off corneal cross-linking (9 mW/cm2 and 10 minutes) for treatment of progressive keratoconus were included. Corneal tomographic, densitometric, and aberrometric values obtained using the Pentacam HR were recorded at the baseline and 3, 6 and 12 months post-operatively.

Result: One year after treatment, corrected distance visual acuity (CDVA) was improved and maximum keratometry, thinnest pachymetry, higher order and total Root Mean Square (RMS) were significantly decreased. (p<0.001) Corneal densitometry values showed a significant elevation 3 months post-surgery compared to baseline and then decreases to baseline values at 1 year. Only the Anterior 0-2 mm zone densitometry at third month was different between three groups. RMS at 1 year had correlation with Anterior 0-2 mm, Anterior 2-6 mm, total corneal 0-2 mm and total corneal 2-6 mm densitometry values at third month. Final CDVA at 12th month follow up correlated with the Anterior 0-2 mm corneal densitometry on third month.

Conclusion: Anterior 0-2 mm zone densitometry at third month post accelerated CCL can be used to detect different staging of keratoconus. Due to correlation between final aberrometric and peak densitometric values in keratoconic eyes, peak densitometric values can be used as a prognostic factor for the final visual outcomes after accelerated CCL.

1. Background

Keratoconus is a progressive non-inflammatory disease of the cornea that defines by conical protrusion, progressive thinning and changes in biomechanical properties of the cornea (1). Although it has been reported in different ethnicities, published data shows more prevalence in Asian countries. (2) In Iran the prevalence has been estimated around 2.5 to 3.3 percent. (3)

Histopathology of the disease involves progressive degeneration of the corneal basal epithelial layer followed by bowman layer collagen fragmentation. The cornea progressively thins in the course of the disease and it suffers from collagen lamellar disarrangement. (4)

Treatment of the disease consists of refractive correction and halting the disease progression. Corneal collagen cross-linking (CCL) now approved by FDA (Food and Drug Administration), prevent the disease progression by means of ultraviolet light and riboflavin to cross link between corneal collagen fibers. (5) Previously it was thought that cross-linking was not that much effective for advanced stages of the disease, although cross-linking was shown to be safe and effective for advanced progressive keratoconus. (6, 7) Wollensak et al. developed a standard Corneal collagen cross-linking (CCL) protocol in 2003 (8), while accelerated protocols have been introduced since 2010 to achieve some advantages such as shortened operation time and reduced rate of complications (9). Accelerated protocols are carried out in 3, 5, and 10 minutes using 30, 18, and 9 mW/cm² irradiance, respectively without changing the cumulative irradiation dose of 5.4 J/cm². Kobashi et al, in the meta-analysis of randomized controlled trials showed accelerated CCL had comparable efficacy and safety profile to standard CCL, and both methods similarly stopped the KCN progression (10).

The Pentacam HR (Oculus Inc., Wetzlar, Germany) using a rotating Scheimpflug camera system can analyze high-resolution anterior segment images, and it is possible to perform a tomographic evaluation of the anterior and posterior surfaces of the cornea, change in elevation, anterior chamber angle width, anterior chamber depth and volume, corneal aberrations, and densitometric evaluation of the cornea (11).

There are a variety of studies evaluating changes in tomographic, densitometric, and aberrometric parameters in keratoconic eyes after accelerated corneal cross-linking (CCL) with the Pentacam HR imaging system (12, 13, 14, 15). Few studies in the literature investigated the outcomes of accelerated CCL using Pentacam (14, 15). The parameters most evaluated are maximum keratometry, central corneal thickness, visual acuity, spherical equivalent, and corneal biomechanical properties. Few studies have evaluated posterior corneal parameters and densitometric probertites, and they have mainly investigated the effects of standard CCL (16, 17). As far as the authors are aware, there is no study in the literature evaluating effects of accelerated CCL on topographic, densitometric, aberrometric and visual parameters in keratoconic eyes based on different stages of the disease.

In this prospective, observational study with 1 year follow up period, we study the natural course of tomographic, densitometric, aberrometric and visual parameters changes after accelerated CCL in keratoconic eyes and we classify densitometric course in different stages of the disease separately (Mild, moderate and severe).

2. Methods

In a prospective, observational study, 67 patients over 18 years of age with progressive keratoconus that underwent accelerated corneal cross linking, between December, 2020 and March, 2021 in the Farabi Eye Hospital were enrolled in this study. Approval for the study was obtained from the ethical committee of the Tehran University of Medical Sciences which is in compliance with the Helsinki Declaration. All patients received a thorough explanation of the study design and aims, and were provided with written informed consent.

All patients were selected according following inclusion and exclusion criteria

Inclusion criteria:

Age between 18 and 40 years, phakic patients with clear lenses, confirmed bilateral keratoconus (KCN) based on clinical and topography findings, bilateral minimum corneal thickness of 380 µm as measured with the Pentacam.

Exclusion criteria:

The following cases were excluded from the study: nonprogressive KC, pseudophakia, corneal thinnest point less than 380 µm, corneal scarring in either eye, ocular surface or tear problems, and coexistence of ocular pathologies other than KCN, systemic autoimmune diseases, pregnancy, history of corneal persistent epithelial defect, previous eye surgery, pterygium or other corneal surface irregularities.

Preoperative evaluation

It included the following clinical history (age, complaint, ocular trauma or disease, optical correction, any systemic medical diseases and etc.) and ophthalmic examination: (1) The BCVA: after refraction, BCVA with rigid gas permeable contact lens was measured as well. Visual acuity was then converted to logMAR scale for analysis. (2) Slit‑lamp biomicroscopy: the cornea was examined for evidence of corneal edema, corneal scars; the anterior chamber was examined for depth, regularity, aqueous flare, and cells; and goldometn applanation tonometery to record baseline intraocular pressure. (3) Fundus examination: auxiliary lenses (+ 90 D lenses) were used to examine central and mid-peripheral retina to exclude possible pathology, for example macular scars or retinal breaks.

Special investigations

Scheimpflug Tomography imaging was obtained with pentacam device (Oculus Pentacam®; Oculus Optikgerate GmbH, Wetzlar, Germany).

Tomographic, densitometric, corneal wavefront aberrometric data was extracted from pentacam analysis printouts.

All examinations and measurements were repeated at 3, 6 and 12 months post-operative visits. The measurements were performed by an expert operator in a standard dim room that was equal for all patients. They were repeated three times for each patient.

Only examinations with ‘‘OK’’ quality check were considered for analysis. The Scheimpflug densitometry is Shown quantitatively in grayscale units (GSUs), which defines Corneal light backscatter on a scale of 0 (minimum scatter) to 100 (maximum scatter). Scheimpflug tomography can analyze the corneal densitometry at different concentric zones (0.0 to 2.0, 2.0 to 6.0, 6.0 to 10.0, 10.0 to 12.0 mm) and depths of the cornea; anterior stromal layer (120 micron), posterior stromal layer (60 micron) and middle layer between these anterior and posterior layers. (Fig. 1)

Surgical technique

Accelerated epithelium-off protocol was used for corneal cross linking.

Under topical anesthesia using Tetracine eye drops (Anestocine, SinaDaru pharmaceutical co, Tehran, Iran) standard preparation of eyelid skin with povidine iodine performed. After opening eyelids apart with lid speculum, for corneal debridement, 20% of alcohol was applied for 20 seconds on the cornea. Riboflavin 0.1% isotonic solutions in 20% dextran T-500 (SinaDaru pharmaceutical co, Tehran, Iran) was instilled every minutes on corneal surface for 20 minutes. Then the cornea was irradiated with UV light of 365nm wavelength for 10 minutes. Corneal collagen cross linking was performed with LightLink CCL device (LightMed, San clemente, CA, USA). The device irradiation power was adjusted to 9.0 mW/cm2. During the light-on period the corneal surface was moistened with riboflavin eye drops every 2 minutes.

In conclusion of the surgery antibiotic eye drops applied over ocular surface and Bandage contact lens (pure vision, Bausch and Lomb co.) was placed over the cornea.

The routine post-operative drug regimen was antibiotic eye drops (Ofloxacin ophthalmic eye drops 0.5%, SinaDaru pharmaceutical co, Tehran, Iran) every 6 hours and betamethasone eye drops (Bethasonate 0.1%,SinaDaru pharmaceutical co, Tehran, Iran) every 6 hours. Antibiotic was continued till the cornea was fully re-epithelialized and then contact lens was removed. Steroid drops was continued for 20 days after surgery.

Statistical analysis

Statistical analysis was performed using SPSS software version 24 (SPSS Inc., IBM). Kolmogorov-Smirnov test was used to assess the normality of the parameters. Wilcoxon test was used to evaluate the changes of parameters before and after surgery in each follow up visit. The potential correlation between parameters was assessed using spearman correlation coefficient. Kruskal Wallis H test and Mann-Whitney U test was used to evaluate the mean difference in parameters in different groups (mild, moderate and severe disease group). For all statistical analysis P value of less than 0.05 was considered statistically significant.

3. Results

From 92 patients underwent accelerated corneal cross linking for progressive keratoconus, 25 patients failed to complete their follow up visits and were excluded from the study.

67 Keratoconus eyes (43 male, 24 female) which completed the follow up were included for statistical analysis.

Baseline characteristics of patients in the study is summarized in Table 1.

 
 
Table 1

Baseline characteristic of kertoconus patients in the study

Parameter

Value

Sex, N (%)

Male

Female

43(64%)

24(36%)

Age (year)

Mean ± SD

Median

Range

21.40 ± 3.7

20

18–36

Number of patients with different Stages of keratoconus (%)

Mild (Amsler-krumeich stage1)

Moderate (Amsler-krumeich stage2)

Severe (Amsler-krumeich stage3,4)

37 (55.2)

19 (28.4)

11 (16.4)

CDVA (logMAR) [Mean ± SD]

0.28 ± 0.3

CDVA with contact lens (logMAR) [Mean ± SD]

0.07 ± 0.12

Mean refractive spherical equivalent [Mean ± SD]

-3.72 ± 3

maximum keratometry [Mean ± SD]

55.68 ± 7

Mean keraometry [Mean ± SD]

48.43 ± 4.7

The patients were classified to Mild, Moderate or severe keratoconus using Amsler-Krumeich classification (Table 2) and sage 1 were categorized as mild, stage 2 as moderate and stage 3 and 4 of the Amsler-Krumeich was categorized as severe in this study.

37 patients were classified as mild, 19 patients as moderate and 11 patients as severe keratoconus.

 
 
Table 2

The Amsler-Krumeich classification

Grades

Characteristics

Stage 1

Eccentric steeping

Myopia and astigmatism < 5.00 D

Mean central K readings < 48.00 D

Stage 2

Myopia and astigmatism from 5.00 to 8.00 D

Mean central K readings < 53.00 D

Absence of scarring

Minimum corneal thickness > 400 µm

Stage 3

Myopia and astigmatism from 8.00 to 10.00 D

Mean central K readings > 53.00 D

Absence of scarring

Minimum corneal thickness from 300 to 400 µm

Stage 4

Refraction not measurable

Mean central K readings > 55.00 D

Central corneal scarring

Minimum corneal thickness 200 µm

Stage is determined if one of the characteristics applies. D: Diopter; K: Keratometry.

Tomographic parameters change after accelerated CCL

The pre-operative and post-operative corneal tomographic parameters incliding Km (Mean Keratometry), Kmax (maximum keratometry of the cornea), thinnest pachymetry (Thinnest point pachymetric thickness) and corrected distance visual acuity (CDVA) and mean refraction spherical equivalent is summarized in Table 3.

There was an initial increase in mean CDVA (logMAR) in 3 month post-operative visit from baseline value from 0.285 to 0.287 (p = 0.03), then the value return to baseline on the sixth month post-operative visit (p = 0.374) and after 1 year the value decreases to less than baseline values. (p = 0.001; table-3; figure − 2) it shows statistically significant improved vision 1 year after treatment. 

There was a statistically significant decrease in maximum keratometry after 1 year of treatment. (p < 0.001) the mean keratometry was also decreased after 1 year and it was statistically significant. (p = 0.016; table-3; figure − 3)

The mean thinnest pachymetry after 1 year decreased approximately 9 microns and it was statistically significant. (p < 0.001)

The absolute value of mean Spherical equivalent decreased from baseline through 3, 6 and 12 month post-operatively.

  
 
Table 3

Refractive, Tomographic and aberrometric parameters change after accelerated CCL

Parameter

(Mean ± SD)

Baseline

3month

6month

12month

CDVA (logMAR)

0.285 ± 0.3

0.287 ± 0.3

0.285 ± 0.3

0.232 ± 0.22

Spherical equivalent (D)

-3.72 ± 3

-3.2 ± 2.6

-3.1 ± 2.8

-2.9 ± 2.9

Kmax (D)

55.6 ± 7

55.7 ± 7.1

55.8 ± 7

54.6 ± 6.9

Kmean (D)

48.4 ± 4.6

48.5 ± 4.7

48.5 ± 4.6

48.2 ± 4.7

Thinnest pachymetry

445.9 ± 42

433.5 ± 41

446.2 ± 42

436.8 ± 40

Total Root Mean Square

9.48 ± 4.99

12.52 ± 6.5

9.61 ± 5.1

9.1 ± 4.7

high order aberration Root mean Square

2.35 ± 1.3

2.6 ± 1.5

2.3 ± 1.3

2.14 ± 1.4

 

Corneal aberrometric values change after accelerated CCL

Table 3 shows changes of mean Root Mean Square (RMS) changes for both total and high order corneal aberrations at baseline, on the third month, on the sixth month and 1 year post operatively. The total RMS values increase on the third month from baseline value. (p < 0.001) it decreases after month 3 through month 12 and at 1 year follow up visit the total RMS value decreases significantly from baseline. (p = 0.009; Table 3)

High order RMS values increases from baseline at month 3. (p < 0.001) and it decreases from month 3 to month 12 and at 1 year follow up visit the high order RMS decreased from pre-operative value that was statistically significant. (p = 0.008; Fig. 4

Corneal densitometric value changes after accelerated CCL

Table-4 shows mean corneal densitometric values for different concentric zones and different depths.

Anterior corneal densitometry 0–2 mm increases from baseline to its peak on third month (from 22.48 to 30.1) (p < 0.001), then it decreases from month 3 to 6 but it’s still increased compare to pre-operative value (p < 0.001). After 1 year the densitometry value return to baseline and there is no statistically significant difference from pre-operative value (p = 0.9; Table 4; Fig. 5).

Anterior corneal densitometry 2–6 mm zone increases from baseline to its peak at month 3 (from 18.82 to 27.58) (p < 0.001) then it decreases from month 3 to 6 but it’s still increased compare to pre-operative value. (p = 0.039) after 1 year the densitometry value return to baseline and there is no statistically significant difference from pre-operative value (p = 0.075).

Middle corneal densitometry in 0–2 mm zone increases from baseline to its peak at month 3 (from 14.9 to 18.34) (p < 0.001) then it decreases from month 3 to 6 but it’s still increased compare to pre-operative value. (p < 0.001) after 1 year the densitometry value return to baseline and there is no statistically significant difference from pre-operative value. (p = 0.1; Table 4; Fig. 6)

Middle corneal densitometry in 2–6 mm zone increases from baseline to its peak at month 3 (from 12.07 to 16.3; p < 0.001), then it decreases from month 3 to 6 but it’s still increased compare to pre-operative value. (p < 0.001) after 1 year the densitometry value return to baseline and there is no statistically significant difference from pre-operative value (p = 0.21).

The densitometric value changes in other zones are summarized in Table 4.

 
 
Table 4

Densitometric parameters change after accelerated CCL

Corneal densitometry

Mean (gray scale unit) ± SD

Baseline

3month

6month

12month

Anterior

0–2 mm zone

2–6 mm zone

Total

22.48 ± 2.6

18.82 ± 1.5

19.80 ± 2.1

30.1 ± 3.9

27.58 ± 2.8

23.75 ± 2.7

24.82 ± 3.1

18.98 ± 1.7

21.84 ± 2.2

22.45 ± 2.6

18.92 ± 1.7

19.71 ± 2.3

Middle layer

0–2 mm zone

2–6 mm zone

Total

14.09 ± 2.0

12.07 ± 0.82

13.25 ± 1.8

18.34 ± 3.0

16.30 ± 1.3

17.37 ± 2.7

14.88 ± 2.3

13.20 ± 1.0

14.51 ± 2.1

14.18 ± 2.1

11.96 ± 1.1

13.33 ± 1.8

Total corneal 0–2 zone

16.41 ± 1.8

20.35 ± 2.4

18.15 ± 2.1

16.35 ± 1.8

Total corneal 2–6 mm zone

14.08 ± 0.98

18.62 ± 1.87

15.57 ± 1.1

14.09 ± 1.1

Corneal densitometric changes in different stages of keratoconus (mild, moderate and severe) is summarized in Table 5, and the difference between these stages in different follow-ups is summarized in Tables 6, 7 and 8. The results of Kruskal Wallis H test show that among all the analyzes performed, in the follow-up on the third month, the Anterior 0–2 mm zone densitometry of all three groups was different and in this follow-up, using Mann Whitney U test the densitometric results of mild and severe group were shown to be statistically different.(Table 7)

 
 
Table 5

Corneal densitometric changes in different stages of keratoconus after accelerated CCL

Corneal densitometry

Mean (gray scale unit) ± SD

Baseline

3month

6month

12month

Anterior 0–2 mm zone

Mild

Moderate

Severe

21.82 ± 2.0

23.01 ± 2.5

23.82 ± 3.8

29.07 ± 2.8

30.58 ± 4.2

32.74 ± 5.4

24.01 ± 2.1

25.47 ± 2.9

26.44 ± 5.0

21.9 ± 2.1

22.71 ± 2.5

23.75 ± 3.5

Anterior 2–6 mm zone

Mild

Moderate

Severe

18.72 ± 1.4

19.01 ± 1.6

18.86 ± 1.8

27.24 ± 2.5

27.69 ± 2.9

28.55 ± 3.6

18.78 ± 1.6

19.29 ± 1.9

19.13 ± 2.1

18.85 ± 1.5

19.13 ± 1.8

18.80 ± 2.1

Total anterior

Mild

Moderate

Severe

19.87 ± 2.0

19.84 ± 2.1

19.5 ± 2.4

23.67 ± 2.6

24.06 ± 2.7

23.44 ± 3.4

21.94 ± 2.1

21.9 ± 2.3

21.43 ± 2.6

19.86 ± 1.9

19.68 ± 2.5

19.26 ± 2.9

Middle layer 0–2 mm

Mild

Moderate

severe

14.00 ± 1.1

13.87 ± 1.0

14.84 ± 4.5

18.07 ± 1.7

18.00 ± 1.2

19.82 ± 6.4

14.92 ± 1.3

14.43 ± 1.3

15.53 ± 4.9

14.12 ± 1.3

13.93 ± 1.1

14.8 ± 4.5

Middle layer 2–6 mm

Mild

Moderate

Severe

12.11 ± 2.8

12.05 ± 0.8

11.96 ± 1.2

16.13 ± 1.1

16.32 ± 1.1

16.8 ± 2.0

13.20 ± 0.96

13.18 ± 0.95

13.21 ± 1.5

11.98 ± 0.96

11.84 ± 1.2

12.1 ± 1.1

Total middle layer

mild

moderate

severe

13.14 ± 2.1

13.16 ± 1.2

13.79 ± 3.6

17.29 ± 1.9

17.68 ± 2.1

17.41 ± 5.2

14.31 ± 1.5

14.55 ± 1.5

15.16 ± 3.9

13.22 ± 1.2

13.33 ± 1.1

13.72 ± 3.6

Total corneal 0–2 zone

Mild

Moderate

Severe

16.21 ± 1.2

16.32 ± 1.3

17.15 ± 3.5

20.14 ± 1.8

20.24 ± 2.0

21.23 ± 4.2

17.84 ± 1.4

15.67 ± 1.0

19.42 ± 4

16.23 ± 1.3

16.2 ± 1.0

16.99 ± 3.8

Total corneal 2–6 zone

Mild

Moderate

Severe

14.1 ± 0.9

14.15 ± 0.9

13.86 ± 1.1

18.40 ± 1.7

18.6 ± 1.8

19.41 ± 2.1

15.58 ± 1.1

15.67 ± 1.0

15.36 ± 1.3

14.02 ± 1.0

14.01 ± 1.1

14.44 ± 1.5

 
 
Table 6

Corneal densitometry parameters at baseline visit for each group

Corneal densitometry at baseline

mild

moderate

severe

P

Anterior 0–2 mm zone

21.82 ± 2.0

23.01 ± 2.5

23.82 ± 3.8

0.117

Anterior 2–6 mm zone

18.72 ± 1.4

19.01 ± 1.6

18.86 ± 1.8

0.825

Total anterior

19.87 ± 2.0

19.84 ± 2.1

19.5 ± 2.4

0.791

middle 0–2 mm zone

14.00 ± 1.1

13.87 ± 1.0

14.84 ± 4.5

0.917

middle 2–6 mm zone

12.11 ± 2.8

12.05 ± 0.8

11.96 ± 1.2

0.752

Total middle

13.14 ± 2.1

13.16 ± 1.2

13.79 ± 3.6

0.851

Total corneal 0–2 mm

16.21 ± 1.2

16.32 ± 1.3

17.15 ± 3.5

0.756

Total corneal 2–6 mm

14.1 ± 0.9

14.15 ± 0.9

13.86 ± 1.1

0.772

*Significant differences at P < 0.05. Kruskal Wallis H test P: Difference between 3 groups
 
 
 
Table 7

Corneal densitometry parameters at 3 month follow up visit for each group

Corneal densitometry (GSU) at 3 month follows up

mild

moderate

severe

p

P1

P2

P3

Anterior 0–2 mm zone

29.07 ± 2.8

30.58 ± 4.2

32.74 ± 5.4

0.035*

0.118

0.013*

0.245

Anterior 2–6 mm zone

27.24 ± 2.5

27.69 ± 2.9

28.55 ± 3.6

0.396

0.446

0.206

0.491

Total anterior

23.67 ± 2.6

24.06 ± 2.7

23.44 ± 3.4

0.925

0.924

0.685

0.779

middle 0–2 mm zone

18.07 ± 1.7

18.00 ± 1.2

19.82 ± 6.4

0.645

0.899

0.308

0.559

middle 2–6 mm zone

16.13 ± 1.1

16.32 ± 1.1

16.8 ± 2.0

0.247

0.222

0.171

0.358

Total middle

17.29 ± 1.9

17.68 ± 2.1

17.41 ± 5.2

0.439

0.234

0.623

0.387

Total corneal 0–2 mm

20.14 ± 1.8

20.24 ± 2.0

21.23 ± 4.2

0.741

0.742

0.451

0.636

Total corneal 2–6 mm

18.40 ± 1.7

18.6 ± 1.8

19.41 ± 2.1

0.376

0.748

0.175

0.270

*Significant differences at P < 0.05. Kruskal Wallis H test (P: Difference between 3 groups), P1: Difference between Group mild and Group moderate (Mann-Whitney U test), P2: Difference between Group mild and Group severe (Mann-Whitney U test), P3: Difference between Group moderate and Group severe. (Mann-Whitney U test)
 
 
Table 8

Corneal densitometry parameters at 12 month follow up visit for each group

Corneal densitometry at 12 months

mild

moderate

severe

p

Anterior 0–2 mm zone

21.9 ± 2.1

22.71 ± 2.5

23.75 ± 3.5

0.239

Anterior 2–6 mm zone

18.85 ± 1.5

19.13 ± 1.8

18.80 ± 2.1

0.643

Total anterior

19.86 ± 1.9

19.68 ± 2.5

19.26 ± 2.9

0.491

middle 0–2 mm zone

14.12 ± 1.3

13.93 ± 1.1

14.8 ± 4.5

0.716

middle 2–6 mm zone

11.98 ± 0.96

11.84 ± 1.2

12.1 ± 1.1

0.753

Total middle

13.22 ± 1.2

13.33 ± 1.1

13.72 ± 3.6

0.667

Total corneal 0–2 mm

16.23 ± 1.3

16.2 ± 1.0

16.99 ± 3.8

0.959

Total corneal 2–6 mm

14.02 ± 1.0

14.01 ± 1.1

14.44 ± 1.5

0.677

*Significant differences at P < 0.05. Kruskal Wallis H test P: Difference between 3 groups

In order to assess the prognostic value of peak densitometric numbers on final visual outcomes, Correlation between peak densitometric values at 3 month and final visual acuity, total RMS and higher order RMS at 12th month were evaluated and the Table 9 and Table 10 shows this correlation.

Best corrected visual acuity at month 12 had correlation with month 3 middle 0–2 mm zone corneal densitometry (p = 0.013), but there was no correlation between final visual acuity and other densitometric parameters. (Table 9)

Higher order RMS at 1 year had correlation with Anterior 0–2 mm, Anterior 2–6 mm, total corneal 0–2 mm and total corneal 2–6 mm densitometry values at 3th month visit. (Table 10)

 
 
Table 9

Correlation between peak densitometric values at 3 month and final visual acuity

Densitometric values at 3 month

Best corrected visual acuity (logMAR) at 12 month (Correlation coefficient)

P value

Anterior 0–2 mm zone

0.113

0.362

Anterior 2–6 mm zone

0.102

0.409

Total anterior

0.001

0.910

Middle 0–2 mm zone

0.303

0.013*

Middle 2–6 mm zone

0.076

0.540

Total middle

0.033

0.789

Total corneal 0–2 mm

0.171

0.167

Total corneal 2–6 mm

0.146

0.238

*significant correlation with p < 0.05. Spearman test used for correlation evaluation.
 
  
Table 10

Correlation between peak densitometric values at month 3 and Higher order RMS at month 12

Densitometric values at 3 month

Higher order RMS at 12 month (Correlation coefficient)

P value

Anterior 0–2 mm zone

0.485

< 0.001*

Anterior 2–6 mm zone

0.409

0.001*

Total anterior

0.122

0.325

middle 0–2 mm zone

0.222

0.071

middle 2–6 mm zone

0.210

0.088

Total middle

0.157

0.205

Total corneal 0–2 mm

0.339

0.005*

Total corneal 2–6 mm

0.249

0.042*

*Significant correlation with p < 0.05. Spearman test used for correlation evaluation.

Correlation between high order RMS and total RMS at 12th month and densitometric parameters at 3 month was analyzed separately for each disease stages (mild, moderate and severe).

In severe keratoconus there was no correlation between densitometric values and RMS. In moderate KCN subgroup analysis (N = 19) there was correlation between high order RMS at 12th month and third month densitometric values including anterior 0–2 mm (p = 0.008), anterior 2–6 mm densitometry (p = 0.002) and total middle layer densitometry (p = 0.030). In mild keratoconus (N = 37) there was correlation between high order RMS and anterior 2–6 mm densitometry (p = 0.28) and total corneal 0–2 mm (p = 0.035).

These results show that, we should use different densitometric zones at third month to predict the final high order RMS at 1 year in different stage of the disease.

4. Discussion

Corneal cross-linking (CCL) was approved in 2016 by the Food and Drug Administration (FDA) as a drug and device combination for the treatment of progressive keratoconus and corneal ectasia after refractive surgery. CCL procedure strengthens the cornea with the purpose of preventing the progression of keratoconus or post-refractive corneal ectasia. In addition, some studies have shown that CCL may also lead to improvement of topographic and visual outcomes in this patient (18).

Most studies in keratoconus have compared the results of different CCL methods such as standard corneal cross-linking or epithelium-off CCL (epi-off CCL) and trans-epithelial corneal cross-linking (TE-CCL) (19, 20, 21). In recent years, several clinical trials have compared the therapeutic effects between the two methods (19, 22, 23), but studies on accelerated CCL are very limited.

Our study aim was to evaluate tomographic, densitometric and aberrometric parameters in patients with progressive keratoconus after accelerated CCL and, due to the fact that appearance and natural course of postoperative haze after corneal cross-linking will assist physicians in knowing what to anticipate postoperatively and in counseling patients regarding postoperative expectations, firstly, we classify densitometric outcomes of accelerated CCL in different stages of the KCN separately (Mild, moderate and severe), and we have answered the question whether the severity of keratoconus affects the corneal haze after CCL?

Accelerated CCL resulted in stabilization of the CDVA during the 12-month follow-up, which is similar to the results obtained in a number of other studies (24, 25), although most of the comparative studies have reported improvement in CDVA. Although our results demonstrated minimal changes in refractive and keratometric values, CDVA were significantly improved, respectively. This visual improvement after accelerated CCL could be explained by the improvement in the Root mean Square (high order aberration).

The two most significant indicators of visual acuity improvement after corneal cross-linking were the low CDVA (≤ 20/40) preoperatively and high maximum keratometry values (≥ 55 D). In the current study, final CDVA at month 12 had correlation with month 3 middle 0–2 mm zone corneal densitometry.

In our study, we found statistically significant flattening in maximum keratometry value after a one year follow-up. In addition, the flattening effect on maximum and mean keratometry were in agreement with the results of previous well-known studies (13, 24, 25). Current study results are consistent with those of recent studies that have decreased Thinnest pachymetry values after accelerated CCL (21, 27). The physiology of this corneal thinning after CCL still not identified: however, structural changes occur in corneal collagen fibrils (28), such as changes in corneal hydration (29), compression of collagen fibrils (30), edema (31), and keratocyte apoptosis (31) are discussed in the literature. It should be borne in mind that previous studies have shown that due to the difference between the Pentacam and ultrasound pachymetry, pachymetric values obtained from Pentacam printout should be interpreted with extreme discretion (32).

Reported in our study, as in the epidemiological study (33), progressive KCN was more common in men than in women.

Root mean Square (high order aberration) is part of the refractive errors which are not correctable with sphere and cylinder corrections, and it is among errors of the optical system of the eye which can exacerbate the quality of the retinal image (34, 35). Since RMS (HOA) can has a considerable effect on the visual function and contrast sensitivity; It’s considered an important value in the context of vision quality (35, 36). In the present study High order RMS value was significantly increased at 3 months after CCL and then decreased at one year, relative to the baseline. Naderan et al(38), and Ahmed et al(39), reported the efficacy of corneal cross-linking in improving HOA parameters in eyes with progressive KCN, and their results showed that there is a statistically significant correlations between pre-operative values of HOAs parameters with corrected distance visual acuity.

In order to determine the prognostic strength of peak densitometric values on final visual outcomes, Correlation between peak densitometric numbers at 3 month post-op and final visual acuity, total RMS and higher order RMS at 12th month were evaluated, and the results showed that RMS at 1 year had correlation with densitometric values of Anterior 0–2 mm area, Anterior 2–6 mm area and total corneal 0–2 mm and 2–6 mm areas at third month. This relationship means that the aberrometric effect in final visual acuity is due to its effect on densitometric factors.

Mathews et al (40) showed for the first time a correlation between densitometric value changes at 6 month and higher order aberrations after CCL with Dresden protocol. We found a similar correlation in eyes with keratoconus after accelerated CCL.

Also, for the first time we studied densitometric factors in different stages of keratoconus patients, our results showed that relationship between Higher order RMS at 12th months and peak densitometric values has different nature in different stages of progressive KCN.

For example, densitometric values on the severe stage had no correlation with RMS at 12th months. But we showed correlation in multiple zones in mild and moderate stage. These results may bring in to mind that in addition to the possible effect of densitometric factor on aberrometric values, there must be other factors that affect aberrometric factors and consequently the patient's CDVA in more advanced stages of KCN and the prognostic value of peak densitometrc parameters is weekend in more advanced stage of the disease, while this is maybe due to high aberration and distorted retinal images in advanced cases.

CCL often leads to a preliminary decrease of visual acuity postoperatively. This reduction of visual acuity can be attributed to loss of transparency of the corneal stroma (41, 42). During the first months after CCL, an increased corneal densitometry numbers; which is clinically detectable as haze, can often be observed. Pircher et al. showed that standard CCL causes various changes in the corneal stroma which leads to an increased densitometric value, particularly in the anterior (120 µm) and central corneal zone (from 0 to 2 mm) (26).

This study evaluates corneal densitometry values obtained from Pentacam and evaluates their effect on parameters, such as visual acuity and aberrometric features after accelerated CCL in keratoconic eyes at different stages. In the study, corneal densitometry was measured by the Pentacam imagining system which offers the possibility of analyzing densitometric values at different corneal depths and concentric zones.

Our results reveal a significant increase in corneal densitometry numbers in anterior and middle stromal layers 3 months after accelerated CCL. The highest change in densitometry was captured in the anterior layer of the stroma, especially in two central concentric zones (0–2 mm and 2–6 mm zone). This is in agreement with a previous report by Böhm et al who also reported a significant increase in densitometric values especially 3 months after accelerated corneal cross-linking, mostly in the anterior layer of the two central zones in eyes with progressive keratoconus, and their results showed the greatest corneal densitometry changes in the anterior stromal layer within the central 0 to 2 mm zone of the cornea (43).

Our result consistent with the previous studies of corneal densitometry, showed that the significant increase in densitometric values at third month postoperative visit after accelerated CCL is more pronounced in the anterior layer. It can be assumed that stray light in the anterior stromal layer of the cornea is more prominent because of the UV-A light intensity, oxygen concentration decline toward the middle and deeper stromal layers and riboflavin concentration.

Our study first examined and classify densitometric course after accelerated CCL in different stage of KCN. Our result showed no correlation between disease staging and densitometric parameters in different concentric corneal zones and different corneal depths in keratotonic eye before CCL, but after accelerated CCL, during 12 month follow up, only the Anterior 0–2 mm zone densitometry of all three groups was different, and patients in higher stage of keratoconus had higher densitometry values, therefore, it seems that Anterior 0–2 mm zone densitometry at third month post accelerated CCL can be used to detect different staging of KCN. The most likely reason for this finding seems to be a complex corneal healing process after CCL comprising of a transformation of keratocytes into myofibroblasts, which are associated with stromal remodeling, that is followed by haziness in the cornea. These factors are different in nature from the Amsler-Krumeich variables that use for staging keratoconus patients.

Although recent studies did not demonstrate a correlation between corneal densitometry values and visual acuity outcomes (29, 33, 43), in our study, final CDVA at 12th month follow up correlated with the changes in corneal densitometry of the Anterior 0–2 mm zone densitometry on third month (Correlation coefficient: 0.303; P: 0.013). However, due to the lack of a control group, the results of our study could not be confirmed nor rejected; thus, it is necessary to conduct future, further research with a larger sample size, a control group and a longer follow-up period studies in the future.

5. Conclusion

In conclusion, corneal changes induced by accelerated CCL seemed to be greatest in the anterior stromal layer in Anterior 0–2 mm zone, and given the Anterior 0–2 mm zone densitometry difference at third months post CCL; In different stages of keratoconus patients, Anterior 0–2 mm zone densitometry at third month post accelerated CCL can be used to detect different staging of KCN. The increase in Anterior 0–2 mm zone densitometry at third month post CCL and, thus, decrease in corneal transparency, however, influence final visual acuity, since correlation between the changes of corneal densitometry of the anterior stromal layer (0–2 mm), and CDVA was found.

Abbreviations

CCL

Corneal collagen cross-linking

FDA

Food and Drug Administration

KCN

Keratoconus

GSU

GrayScale Units

RMS

Root Mean Square

HOA

Higher order Aberration

CDVA

Corrected Distance Visual Acuity

Kmax

Maximum Keratometry

Kmean

Mean Keratometry

Declarations

Approval for the study was obtained from the ethical committee of the Tehran University of Medical Sciences which is in compliance with the Helsinki Declaration. All patients received a thorough explanation of the study design and aims, and were provided with written informed consent.

References

  1. McGhee CN. 2008 Sir Norman McAlister Gregg lecture: 150 yearsof practical observations on the conical cornea – What have welearned? Clin Exp Ophthalmol 2009;37:16076.
  2. Pearson AR, Soneji B, Sarvananthan N, SandfordSmith JH. Does ethnic origin influence the incidence or severity of keratoconus?Eye (Lond) 2000;14(Pt 4):6258.
  3. Hashemi H, Khabazkhoob M, Fotouhi A. Topographic keratoconus is not rare in an Iranian population: The Tehran eye study. Ophthalmic Epidemiol 2013;20:38591.
  4. Torricelli AA, Singh V, Santhiago MR, Wilson SE. The corneal epithelial basement membrane: Structure, function, and disease. Invest Ophthalmol Vis Sci 2013;54:6390400.
  5. Spoerl E, Huhle M, Seiler T. Induction of crosslinks in corneal tissue. Exp Eye Res 1998;66:97103.
  6. Moghadam RS, Akbari M, Alizadeh Y, Medghalchi A, Dalvandi R. The Outcome of Corneal Collagen Cross-Linking in Patients with Advanced Progressive Keratoconus: A 2-Year Follow-up Study. Middle East Afr J Ophthalmol. 2019 Jan-Mar;26(1):11–16. doi: 10.4103/meajo.MEAJO_101_18. PMID: 31114118; PMCID: PMC6507378.
  7. Giacomin NT, Netto MV, Torricelli AA, Marino GK, Bechara SJ, Espíndola RF, Santhiago MR. Corneal Collagen Cross-linking in Advanced Keratoconus: A 4-Year Follow-up Study. J Refract Surg. 2016 Jul 1;32(7):459 – 65. doi: 10.3928/1081597X-20160429-01. PMID: 27400077
  8. Wollensak G, Spoerl E, Seiler T. Riboflavin/ultraviolet-a-induced collagen crosslinking for the treatment of keratoconus. Am J Ophthalmol. 2003 May;135(5):620-7. doi: 10.1016/s0002-9394(02)02220-1. PMID: 12719068.
  9. Wen D, Li Q, Song B, Tu R, Wang Q, O'Brart DPS, McAlinden C, Huang J. Comparison of Standard Versus Accelerated Corneal Collagen Cross-Linking for Keratoconus: A Meta-Analysis. Invest Ophthalmol Vis Sci. 2018 Aug 1;59(10):3920–3931. doi: 10.1167/iovs.18-24656. PMID: 30073363.
  10. Kobashi H, Tsubota K. Accelerated Versus Standard Corneal Cross-Linking for Progressive Keratoconus: A Meta-Analysis of Randomized Controlled Trials. Cornea. 2020 Feb;39(2):172–180. doi: 10.1097/ICO.0000000000002092. PMID: 31369463.
  11. Belin MW, Ambrósio R. Scheimpflug imaging for keratoconus and ectatic disease. Indian J Ophthalmol. 2013 Aug;61(8):401–6. doi: 10.4103/0301-4738.116059. PMID: 23925323; PMCID: PMC3775073.
  12. Koc M, Tekin K, Tekin MI, Uzel MM, Kosekahya P, Ozulken K, Yilmazbas P. An Early Finding of Keratoconus: Increase in Corneal Densitometry. Cornea. 2018 May;37(5):580–586. doi: 10.1097/ICO.0000000000001537. PMID: 29384807.
  13. Hashemi H, Miraftab M, Seyedian MA, Hafezi F, Bahrmandy H, Heidarian S, Amanzadeh K, Nikbin H, Fotouhi A, Asgari S. Long-term Results of an Accelerated Corneal Cross-linking Protocol (18 mW/cm2) for the Treatment of Progressive Keratoconus. Am J Ophthalmol. 2015 Dec;160(6):1164–1170.e1. doi: 10.1016/j.ajo.2015.08.027. Epub 2015 Aug 24. PMID: 26314662.
  14. Sadoughi MM, Einollahi B, Baradaran-Rafii A, Roshandel D, Hasani H, Nazeri M. Accelerated versus conventional corneal collagen cross-linking in patients with keratoconus: an intrapatient comparative study. Int Ophthalmol. 2018 Feb;38(1):67–74. doi: 10.1007/s10792-016-0423-0. Epub 2016 Dec 29. PMID: 28035498.
  15. KOSEKAHYA, PINAR, et al. "Long Term Evaluation of Posterior Corneal Surface Parameters after Accelerated Corneal Crosslinking with a Comparison with Uncross-Linked Keratoconic Eyes." (2021).
  16. J Lai M, Greenstein SA, Gelles JD, Hersh PS. Corneal Haze After Transepithelial Collagen Cross-linking for Keratoconus: A Scheimpflug Densitometry Analysis. Cornea. 2020 Sep;39(9):1117–1121. doi: 10.1097/ICO.0000000000002334. PMID: 32282356.
  17. Akkaya Turhan S, Toker E. Changes in Corneal Density After Accelerated Corneal Collagen Cross-linking With Different Irradiation Intensities and Energy Exposures: 1-Year Follow-up. Cornea. 2017 Nov;36(11):1331–1335. doi: 10.1097/ICO.0000000000001362. PMID: 28872519.
  18. Alenezi B, Kazaili A, Akhtar R, Radhakrishnan H. Corneal biomechanical properties following corneal cross-linking: Does age have an effect? Exp Eye Res. 2022 Jan;214:108839. doi: 10.1016/j.exer.2021.108839. Epub 2021 Nov 14. PMID: 34785203.
  19. Cunha AM, Sardinha T, Torrão L, Moreira R, Falcão-Reis F, Pinheiro-Costa J. Transepithelial Accelerated Corneal Collagen Cross-Linking: Two-Year Results. Clin Ophthalmol. 2020 Aug 14;14:2329–2337. doi: 10.2147/OPTH.S252940. PMID: 32884233; PMCID: PMC7434572.
  20. Madeira C, Vasques A, Beato J, Godinho G, Torrão L, Falcão M, Falcão-Reis F, Pinheiro-Costa J. Transepithelial accelerated versus conventional corneal collagen crosslinking in patients with keratoconus: a comparative study. Clin Ophthalmol. 2019 Mar 1;13:445–452. doi: 10.2147/OPTH.S189183. PMID: 30880905; PMCID: PMC6402612.
  21. Li W, Wang B. Efficacy and safety of transepithelial corneal collagen crosslinking surgery versus standard corneal collagen crosslinking surgery for keratoconus: a meta-analysis of randomized controlled trials. BMC Ophthalmol. 2017 Dec 28;17(1):262. doi: 10.1186/s12886-017-0657-2. PMID: 29282020; PMCID: PMC5745766.
  22. Nawaz S, Gupta S, Gogia V, Sasikala NK, Panda A. Trans-epithelial versus conventional corneal collagen crosslinking: A randomized trial in keratoconus. Oman J Ophthalmol. 2015 Jan-Apr;8(1):9–13. doi: 10.4103/0974-620X.149855. PMID: 25709267; PMCID: PMC4333560.
  23. Akbar B, Intisar-Ul-Haq R, Ishaq M, Fawad A, Arzoo S, Siddique K. Comparison of transepithelial corneal crosslinking with epithelium-off crosslinking (epithelium-off CXL) in adult Pakistani population with progressive keratoconus. Taiwan J Ophthalmol. 2017 Oct-Dec;7(4):185–190. doi: 10.4103/tjo.tjo_38_17. PMID: 29296550; PMCID: PMC5747228.
  24. Kırgız A, Atalay K, Çabuk KŞ, Kaldırım H, Taşkapılı M. Factors affecting visual acuity after accelerated crosslinking in patients with progressive keratoconus. Arq Bras Oftalmol. 2016 May-Jun;79(3):151-4. doi: 10.5935/0004-2749.20160046. PMID: 27463624.
  25. Omar IAN, Zein HA. Accelerated Epithelium-Off Corneal Collagen Cross-Linking For Keratoconus: 12-Month Results. Clin Ophthalmol. 2019 Dec 4;13:2385–2394. doi: 10.2147/OPTH.S232118. PMID: 31824132; PMCID: PMC6900281.
  26. Pircher N, Pachala M, Prager F, Pieh S, Schmidinger G. Changes in straylight and densitometry values after corneal collagen crosslinking. J Cataract Refract Surg. 2015 May;41(5):1038–43. doi: 10.1016/j.jcrs.2014.07.043. Epub 2015 May 5. PMID: 25953473.
  27. Greenstein SA, Shah VP, Fry KL, Hersh PS. Corneal thickness changes after corneal collagen crosslinking for keratoconus and corneal ectasia: one-year results. J Cataract Refract Surg. 2011 Apr;37(4):691–700. doi: 10.1016/j.jcrs.2010.10.052. PMID: 21420594.
  28. Seiler T, Hafezi F. Corneal cross-linking-induced stromal demarcation line. Cornea. 2006 Oct;25(9):1057-9. doi: 10.1097/01.ico.0000225720.38748.58. PMID: 17133053.
  29. Bottós KM, Dreyfuss JL, Regatieri CV, Lima-Filho AA, Schor P, Nader HB, Chamon W. Immunofluorescence confocal microscopy of porcine corneas following collagen cross-linking treatment with riboflavin and ultraviolet A. J Refract Surg. 2008 Sep;24(7):S715-9. doi: 10.3928/1081597X-20080901-14. PMID: 18811117.
  30. Wollensak G, Aurich H, Pham DT, Wirbelauer C. Hydration behavior of porcine cornea crosslinked with riboflavin and ultraviolet A. J Cataract Refract Surg. 2007 Mar;33(3):516 – 21. doi: 10.1016/j.jcrs.2006.11.015. PMID: 17321404.
  31. DOHLMAN CH, HEDBYS BO, MISHIMA S. The swelling pressure of the corneal stroma. Invest Ophthalmol. 1962 Apr;1:158–62. PMID: 13886960.
  32. Kosekahya P, Koc M, Yalcinsoy KO, Kocabas DO, Toker MI. Comparative Evaluation of Central Corneal Thickness in Cross-Linked Keratoconic Eyes. Cornea. 2020 Sep;39(9):1080–1085. doi: 10.1097/ICO.0000000000002339. PMID: 32304430.
  33. Rafati S, Hashemi H, Nabovati P, Doostdar A, Yekta A, Aghamirsalim M, Khabazkhoob M. Demographic profile, clinical, and topographic characteristics of keratoconus patients attending at a tertiary eye center. J Curr Ophthalmol. 2019 Feb 23;31(3):268–274. doi: 10.1016/j.joco.2019.01.013. PMID: 31528760; PMCID: PMC6742613.
  34. Hofer H, Chen L, Yoon GY, Singer B, Yamauchi Y, Williams DR. Improvement in retinal image quality with dynamic correction of the eye's aberrations. Opt Express. 2001 May 21;8(11):631 – 43. doi: 10.1364/oe.8.000631. PMID: 19421252.
  35. Thibos LN. Retinal image quality for virtual eyes generated by a statistical model of ocular wavefront aberrations. Ophthalmic Physiol Opt. 2009 May;29(3):288 – 91. doi: 10.1111/j.1475-1313.2009.00662.x. PMID: 19422560.
  36. Yamaguchi T, Negishi K, Ohnuma K, Tsubota K. Correlation between contrast sensitivity and higher-order aberration based on pupil diameter after cataract surgery. Clin Ophthalmol. 2011;5:1701–7. doi: 10.2147/OPTH.S21819. Epub 2011 Dec 7. PMID: 22205829; PMCID: PMC3245194.
  37. Oshika T, Okamoto C, Samejima T, Tokunaga T, Miyata K. Contrast sensitivity function and ocular higher-order wavefront aberrations in normal human eyes. Ophthalmology. 2006 Oct;113(10):1807–12. doi: 10.1016/j.ophtha.2006.03.061. Epub 2006 Jul 31. PMID: 16876865.
  38. Naderan M, Jahanrad A. Higher-order aberration 4 years after corneal collagen cross-linking. Indian J Ophthalmol. 2017 Sep;65(9):808–812. doi: 10.4103/ijo.IJO_21_17. PMID: 28905822; PMCID: PMC5621261.
  39. MOHAMED, M. K., AHMED, M. M., KHALED, A., & HISHAM, A. (2021). Changes in Ocular Higher-Order Aberrations after Corneal Collagen Cross-Linking for Keratoconus. The Medical Journal of Cairo University, 89(September), 1367–1371.
  40. Mathews PM, De Rojas JO, Rapuano PB, Zemsky CJ, Florakis GJ, Trokel SL, Suh LH. Correlation of Scheimpflug densitometry changes with clinical outcomes after corneal crosslinking. J Cataract Refract Surg. 2018 Aug;44(8):993–1002. doi: 10.1016/j.jcrs.2018.05.016. PMID: 30115299.
  41. Asri D, Touboul D, Fournié P, Malet F, Garra C, Gallois A, Malecaze F, Colin J. Corneal collagen crosslinking in progressive keratoconus: multicenter results from the French National Reference Center for Keratoconus. J Cataract Refract Surg. 2011 Dec;37(12):2137-43. doi: 10.1016/j.jcrs.2011.08.026. PMID: 22108109.
  42. Seiler TG, Schmidinger G, Fischinger I, Koller T, Seiler T. Komplikationen der Vernetzung der Hornhaut [Complications of corneal cross-linking]. Ophthalmologe. 2013 Jul;110(7):639 – 44. German. doi: 10.1007/s00347-012-2682-0. PMID: 23760422.
  43. Böhm M, Shajari M, Remy M, Kohnen T. Corneal densitometry after accelerated corneal collagen cross-linking in progressive keratoconus. Int Ophthalmol. 2019 Apr;39(4):765–775. doi: 10.1007/s10792-018-0876-4. Epub 2018 Mar 26. PMID: 29582259.