Cell death due to apoptosis was identified via TUNEL assay using an in situ cell death detection kit (TMR Red; Roche, Nutley, NJ, USA). Corneas were harvested and fixed in 10% formalin overnight, and then in 30% sucrose for 3 days. Frozen corneal cross sections were cut at 12 μm using a cryostat, and mounted sections were stored at −80°C. For TUNEL staining, corneas were briefly washed in phosphate-buffered saline (PBS), permeabilized with PBS-Triton (T) for 15 minutes, and stained in accordance with the manufacturer's directions. Following TUNEL staining, corneas were processed for Ki67 protein (ab16667; Abcam, Cambridge, MA, USA), a cellular proliferation marker expressed only during active stages of the cell cycle. Corneal tissue was briefly washed in PBS, permeabilized with PBS-T for 15 minutes, and blocked with 5% normal donkey serum in PBS-T for 1 hour. Following blocking solution, 1:100 diluted primary Ki67 antibody was applied and incubated overnight. After brief washes in PBS, tissue was incubated with a FITC anti-rabbit secondary antibody (1:200; Thermo Fisher Scientific, Portsmouth, NH, USA) for 1 hour. Slides were then coverslipped with a Vectashield soft-set medium with 4′,6-diamidino-2-phenylindole (DAPI) (Vector Labs, Burlingame, CA, USA).
Actively proliferating and apoptotic cells were identified using a fluorescence microscope (Leica DM2500M) with LAS V4.6 software. For each cornea, four sagittal sections that crossed through the corneal apex were selected for analysis. In each section, two regions were defined based on their distance from the corneal apex. The central zone consisted of an area within 0.75 mm of the corneal apex, and the peripheral zone was located just outside of the measured corneal center. Harvested corneas did not include the limbus or conjunctiva. Two images were taken from each region by a researcher blinded to the animal's treatment and sex. From these images, TUNEL-positive (TUNEL+) and Ki67-positive (Ki67+) staining was counted by two blinded researchers. Overall counts from the two researchers differed by <5% and were averaged for the final analysis. Comparisons were made between treatment groups after calculating the mean of the counts from all eight regions.