Combined registration and motion correction of longitudinal retinal OCT data

Longitudinal deep network for consistent OCT layer segmentation

Retinal layer thickness is an important bio-marker for people with multiple sclerosis (PwMS). In clinical practice, retinal layer thickness changes in optical coherence tomography (OCT) are widely used for monitoring multiple sclerosis (MS) progression. Recent developments in automated retinal layer segmentation algorithms allow cohort-level retina thinning to be observed in a large study of PwMS. However, variability in these results make it difficult to identify patient-level trends; this prevents patient specific disease monitoring and treatment planning using OCT. Deep learning based retinal layer segmentation algorithms have achieved state-of-the-art accuracy, but the segmentation is performed on each individual scan without utilizing longitudinal information, which can be important in reducing segmentation error and reveal subtle changes in retinal layers. In this paper, we propose a longitudinal OCT segmentation network which achieves more accurate and consistent layer thickness measurements for PwMS.

Trans-Synaptic Degeneration Following Acute Optic Neuritis in Multiple Sclerosis

OBJECTIVE

To explore longitudinal changes in brain volumetric measures and retinal layer thicknesses following acute optic neuritis (AON) in people with multiple sclerosis (PwMS), to investigate the process of trans-synaptic degeneration, and determine its clinical relevance.

METHODS

PwMS were recruited within 40 days of AON onset (n = 49), and underwent baseline retinal optical coherence tomography and brain magnetic resonance imaging followed by longitudinal tracking for up to 5 years. A comparator cohort of PwMS without a recent episode of AON were similarly tracked (n = 73). Mixed-effects linear regression models were used.

RESULTS

Accelerated atrophy of the occipital gray matter (GM), calcarine GM, and thalamus was seen in the AON cohort, as compared with the non-AON cohort (-0.76% vs -0.22% per year [p = 0.01] for occipital GM, -1.83% vs -0.32% per year [p = 0.008] for calcarine GM, -1.17% vs -0.67% per year [p = 0.02] for thalamus), whereas rates of whole-brain, cortical GM, non-occipital cortical GM atrophy, and T2 lesion accumulation did not differ significantly between the cohorts. In the AON cohort, greater AON-induced reduction in ganglion cell+inner plexiform layer thickness over the first year was associated with faster rates of whole-brain (r = 0.32, p = 0.04), white matter (r = 0.32, p = 0.04), and thalamic (r = 0.36, p = 0.02) atrophy over the study period. Significant relationships were identified between faster atrophy of the subcortical GM and thalamus, with worse visual function outcomes after AON.

INTERPRETATION

These results provide in-vivo evidence for anterograde trans-synaptic degeneration following AON in PwMS, and suggest that trans-synaptic degeneration may be related to clinically-relevant visual outcomes. ANN NEUROL 2023;93:76-87.

Association of Vitamin D Polygenic Risk Scores and Disease Outcome in People With Multiple Sclerosis

Background and Objectives

Observational studies suggest low levels of 25-hydroxyvitamin D (25[OH]D) may be associated with increased disease activity in people with multiple sclerosis (PwMS). Large-scale genome-wide association studies (GWAS) suggest 25(OH)D levels are partly genetically determined. The resultant polygenic scores (PGSs) could serve as a proxy for 25(OH)D levels, minimizing potential confounding and reverse causation in analyses with outcomes. Herein, we assess the association of genetically determined 25(OH)D and disease outcomes in MS.

Methods

We generated 25(OH)D PGS for 1,924 PwMS with available genotyping data pooled from 3 studies: the CombiRx trial (n = 575), Johns Hopkins MS Center (n = 1,152), and Immune-Mediated Inflammatory Diseases study (n = 197). 25(OH)D-PGS were derived using summary statistics (p < 5 × 10−8) from a large GWAS including 485,762 individuals with circulating 25(OH)D levels measured. We included clinical and imaging outcomes: Expanded disability status scale (EDSS), timed 25-foot walk (T25FW), nine-hole peg test (9HPT), radiologic activity, and optical coherence tomography-derived ganglion cell inner plexiform layer (GCIPL) thickness. A subset (n = 935) had measured circulating 25(OH)D levels. We fitted multivariable models based on the outcome of interest and pooled results across studies using random effects meta-analysis. Sensitivity analyses included a modifiedpvalue threshold for inclusion in the PGS (5 × 10−5) and applying Mendelian randomization (MR) rather than using PGS.

Results

Initial analyses demonstrated a positive association between generated 25(OH)D-PGS and circulating 25(OH)D levels (per 1SD increase in 25[OH]D PGS: 3.08%, 95% CI: 1.77%, 4.42%;p= 4.33e-06; R2= 2.24%). In analyses with outcomes, we did not observe an association between 25(OH)D-PGS and relapse rate (per 1SD increase in 25[OH]D-PGS: 0.98; 95% CI: 0.87–1.10), EDSS worsening (per 1SD: 1.05; 95% CI: 0.87–1.28), change in T25FW (per 1SD: 0.07%; 95% CI: −0.34 to 0.49), or change in 9HPT (per 1SD: 0.09%; 95% CI: −0.15 to 0.33). 25(OH)D-PGS was not associated with new lesion accrual, lesion volume or other imaging-based outcomes (whole brain, gray, white matter volume loss or GCIPL thinning). The results were similarly null in analyses using otherpvalue thresholds or those applying MR.

Discussion

Genetically determined lower 25(OH)D levels were not associated with worse disease outcomes in PwMS and raises questions about the plausibility of a treatment effect of vitamin D in established MS.

Optical Coherence Tomography and Optical Coherence Tomography Angiography Findings After Optic Neuritis in Multiple Sclerosis

Background: In people with multiple sclerosis (MS), optic neuritis (ON) results in inner retinal layer thinning, and reduced density of the retinal microvasculature.

Objective: To compare inter-eye differences (IEDs) in macular optical coherence tomography (OCT) and OCT angiography (OCTA) measures in MS patients with a history of unilateral ON (MS ON) vs. MS patients with no history of ON (MS non-ON), and to assess how these measures correlate with visual function outcomes after ON.

Methods: In this cross-sectional study, people with MS underwent OCT and OCTA. Superficial vascular plexus (SVP) density of each eye was quantified using a deep neural network. IEDs were calculated with respect to the ON eye in MS ON patients, and with respect to the right eye in MS non-ON patients. Statistical analyses used mixed-effect regression models accounting for intra-subject correlations.

Results: We included 43 MS ON patients (with 92 discrete OCT/OCTA visits) and 14 MS non-ON patients (with 24 OCT/OCTA visits). Across the cohorts, mean IED in SVP density was −2.69% (SD 3.23) in MS ON patients, as compared to 0.17% (SD 2.39) in MS non-ON patients (p = 0.002). When the MS ON patients were further stratified according to time from ON and compared to MS non-ON patients with multiple cross-sectional analyses, we identified that IED in SVP density was significantly increased in MS ON patients at 1–3 years (p = < 0.001) and >3 years post-ON (p < 0.001), but not at <3 months (p = 0.21) or 3–12 months post-ON (p = 0.07), while IED in ganglion cell + inner plexiform layer (GCIPL) thickness was significantly increased in MS ON patients at all time points post-ON (p ≦ 0.01 for all). IED in SVP density and IED in GCIPL thickness demonstrated significant relationships with IEDs in 100% contrast, 2.5% contrast, and 1.25% contrast letter acuity in MS ON patients (p < 0.001 for all).

Conclusions: Our findings suggest that increased IED in SVP density can be detected after ON in MS using OCTA, and detectable changes in SVP density after ON may occur after changes in GCIPL thickness. IED in SVP density and IED in GCIPL thickness correlate well with visual function outcomes in MS ON patients.

Numerical method for axial motion artifact correction in retinal spectral-domain optical coherence tomography

A numerical method that compensates image distortions caused by random fluctuations of the distance to an object in spectral-domain optical coherence tomography (SD OCT) has been proposed and verified experimentally. The proposed method is based on the analysis of the phase shifts between adjacent scans that are caused by micrometer-scale displacements and the subsequent compensation for the displacements through phase-frequency correction in the spectral space. The efficiency of the method is demonstrated in model experiments with harmonic and random movements of a scattering object as well as during in vivo imaging of the retina of the human eye.

Alterations in the retinal vasculature occur in multiple sclerosis and exhibit novel correlations with disability and visual function measures

BACKGROUND

The retinal vasculature may be altered in multiple sclerosis (MS), potentially acting as a biomarker of disease processes.

OBJECTIVE

To compare retinal vascular plexus densities in people with MS (PwMS) and healthy controls (HCs), and examine correlations with visual function and global disability.

METHODS

In this cross-sectional study, 111 PwMS (201 eyes) and 50 HCs (97 eyes) underwent optical coherence tomography angiography (OCTA). Macular superficial vascular plexus (SVP) and deep vascular plexus (DVP) densities were quantified, and poor quality images were excluded according to an artifact-rating protocol.

RESULTS

Mean SVP density was 24.1% (SD = 5.5) in MS eyes (26.0% (SD = 4.7) in non-optic neuritis (ON) eyes vs. 21.7% (SD = 5.5) in ON eyes, p < 0.001), as compared to 29.2% (SD = 3.3) in HC eyes (p < 0.001 for all MS eyes and multiple sclerosis optic neuritis (MSON) eyes vs. HC eyes, p = 0.03 for MS non-ON eyes vs. HC eyes). DVP density did not differ between groups. In PwMS, lower SVP density was associated with higher levels of disability (expanded disability status scale (EDSS): R2 = 0.26, p = 0.004; multiple sclerosis functional composite (MSFC): R2 = 0.27, p = 0.03) and lower letter acuity scores (100% contrast: R2 = 0.29; 2.5% contrast: R2 = 0.40; 1.25% contrast: R2 = 0.31; p < 0.001 for all).

CONCLUSIONS

Retinal SVP density measured by OCTA is reduced across MS eyes, and correlates with visual function, EDSS, and MSFC scores.

Layer boundary evolution method for macular OCT layer segmentation

Optical coherence tomography (OCT) is used to produce high resolution depth images of the retina and is now the standard of care for in-vivo ophthalmological assessment. It is also increasingly being used for evaluation of neurological disorders such as multiple sclerosis (MS). Automatic segmentation methods identify the retinal layers of the macular cube providing consistent results without intra- and inter-rater variation and is faster than manual segmentation. In this paper, we propose a fast multi-layer macular OCT segmentation method based on a fast level set method. Our framework uses contours in an optimized approach specifically for OCT layer segmentation over the whole macular cube. Our algorithm takes boundary probability maps from a trained random forest and iteratively refines the prediction to subvoxel precision. Evaluation on both healthy and multiple sclerosis subjects shows that our method is statistically better than a state-of-the-art graph-based method.

The Relationship Between High-Order Aberration and Anterior Ocular Biometry During Accommodation in Young Healthy Adults

Purpose

This study investigated the anterior ocular anatomic origin of high-order aberration (HOA) components using optical coherence tomography and a Shack-Hartmann wavefront sensor.

Methods

A customized system was built to simultaneously capture images of ocular wavefront aberrations and anterior ocular biometry. Relaxed, 2-diopter (D) and 4-D accommodative states were repeatedly measured in 30 young subjects. Custom software was used to correct optical distortions and measure biometric parameters from the images.

Results

The anterior ocular biometry changed during 2-D accommodation, in which central lens thickness, ciliary muscle thicknesses at 1 mm posterior to the scleral spur (CMT1), and the maximum value of ciliary muscle thickness increased significantly, whereas anterior chamber depth, CMT3, radius of anterior lens surface curvature (RAL), and radius of posterior lens surface curvature (RPL) decreased significantly. The changes in the anterior ocular parameters during 4-D accommodation were similar to those for the 2-D accommodation. \begin{document}\newcommand{\bialpha}{\boldsymbol{\alpha}}\newcommand{\bibeta}{\boldsymbol{\beta}}\newcommand{\bigamma}{\boldsymbol{\gamma}}\newcommand{\bidelta}{\boldsymbol{\delta}}\newcommand{\bivarepsilon}{\boldsymbol{\varepsilon}}\newcommand{\bizeta}{\boldsymbol{\zeta}}\newcommand{\bieta}{\boldsymbol{\eta}}\newcommand{\bitheta}{\boldsymbol{\theta}}\newcommand{\biiota}{\boldsymbol{\iota}}\newcommand{\bikappa}{\boldsymbol{\kappa}}\newcommand{\bilambda}{\boldsymbol{\lambda}}\newcommand{\bimu}{\boldsymbol{\mu}}\newcommand{\binu}{\boldsymbol{\nu}}\newcommand{\bixi}{\boldsymbol{\xi}}\newcommand{\biomicron}{\boldsymbol{\micron}}\newcommand{\bipi}{\boldsymbol{\pi}}\newcommand{\birho}{\boldsymbol{\rho}}\newcommand{\bisigma}{\boldsymbol{\sigma}}\newcommand{\bitau}{\boldsymbol{\tau}}\newcommand{\biupsilon}{\boldsymbol{\upsilon}}\newcommand{\biphi}{\boldsymbol{\phi}}\newcommand{\bichi}{\boldsymbol{\chi}}\newcommand{\bipsi}{\boldsymbol{\psi}}\newcommand{\biomega}{\boldsymbol{\omega}}\({\rm Z}_4^0\)\end{document} decreased significantly during 2-D accommodation, and \begin{document}\newcommand{\bialpha}{\boldsymbol{\alpha}}\newcommand{\bibeta}{\boldsymbol{\beta}}\newcommand{\bigamma}{\boldsymbol{\gamma}}\newcommand{\bidelta}{\boldsymbol{\delta}}\newcommand{\bivarepsilon}{\boldsymbol{\varepsilon}}\newcommand{\bizeta}{\boldsymbol{\zeta}}\newcommand{\bieta}{\boldsymbol{\eta}}\newcommand{\bitheta}{\boldsymbol{\theta}}\newcommand{\biiota}{\boldsymbol{\iota}}\newcommand{\bikappa}{\boldsymbol{\kappa}}\newcommand{\bilambda}{\boldsymbol{\lambda}}\newcommand{\bimu}{\boldsymbol{\mu}}\newcommand{\binu}{\boldsymbol{\nu}}\newcommand{\bixi}{\boldsymbol{\xi}}\newcommand{\biomicron}{\boldsymbol{\micron}}\newcommand{\bipi}{\boldsymbol{\pi}}\newcommand{\birho}{\boldsymbol{\rho}}\newcommand{\bisigma}{\boldsymbol{\sigma}}\newcommand{\bitau}{\boldsymbol{\tau}}\newcommand{\biupsilon}{\boldsymbol{\upsilon}}\newcommand{\biphi}{\boldsymbol{\phi}}\newcommand{\bichi}{\boldsymbol{\chi}}\newcommand{\bipsi}{\boldsymbol{\psi}}\newcommand{\biomega}{\boldsymbol{\omega}}\({\rm{Z}}_3^{ - 1}\)\end{document}, \begin{document}\newcommand{\bialpha}{\boldsymbol{\alpha}}\newcommand{\bibeta}{\boldsymbol{\beta}}\newcommand{\bigamma}{\boldsymbol{\gamma}}\newcommand{\bidelta}{\boldsymbol{\delta}}\newcommand{\bivarepsilon}{\boldsymbol{\varepsilon}}\newcommand{\bizeta}{\boldsymbol{\zeta}}\newcommand{\bieta}{\boldsymbol{\eta}}\newcommand{\bitheta}{\boldsymbol{\theta}}\newcommand{\biiota}{\boldsymbol{\iota}}\newcommand{\bikappa}{\boldsymbol{\kappa}}\newcommand{\bilambda}{\boldsymbol{\lambda}}\newcommand{\bimu}{\boldsymbol{\mu}}\newcommand{\binu}{\boldsymbol{\nu}}\newcommand{\bixi}{\boldsymbol{\xi}}\newcommand{\biomicron}{\boldsymbol{\micron}}\newcommand{\bipi}{\boldsymbol{\pi}}\newcommand{\birho}{\boldsymbol{\rho}}\newcommand{\bisigma}{\boldsymbol{\sigma}}\newcommand{\bitau}{\boldsymbol{\tau}}\newcommand{\biupsilon}{\boldsymbol{\upsilon}}\newcommand{\biphi}{\boldsymbol{\phi}}\newcommand{\bichi}{\boldsymbol{\chi}}\newcommand{\bipsi}{\boldsymbol{\psi}}\newcommand{\biomega}{\boldsymbol{\omega}}\({\rm{Z}}_3^1\)\end{document}, \begin{document}\newcommand{\bialpha}{\boldsymbol{\alpha}}\newcommand{\bibeta}{\boldsymbol{\beta}}\newcommand{\bigamma}{\boldsymbol{\gamma}}\newcommand{\bidelta}{\boldsymbol{\delta}}\newcommand{\bivarepsilon}{\boldsymbol{\varepsilon}}\newcommand{\bizeta}{\boldsymbol{\zeta}}\newcommand{\bieta}{\boldsymbol{\eta}}\newcommand{\bitheta}{\boldsymbol{\theta}}\newcommand{\biiota}{\boldsymbol{\iota}}\newcommand{\bikappa}{\boldsymbol{\kappa}}\newcommand{\bilambda}{\boldsymbol{\lambda}}\newcommand{\bimu}{\boldsymbol{\mu}}\newcommand{\binu}{\boldsymbol{\nu}}\newcommand{\bixi}{\boldsymbol{\xi}}\newcommand{\biomicron}{\boldsymbol{\micron}}\newcommand{\bipi}{\boldsymbol{\pi}}\newcommand{\birho}{\boldsymbol{\rho}}\newcommand{\bisigma}{\boldsymbol{\sigma}}\newcommand{\bitau}{\boldsymbol{\tau}}\newcommand{\biupsilon}{\boldsymbol{\upsilon}}\newcommand{\biphi}{\boldsymbol{\phi}}\newcommand{\bichi}{\boldsymbol{\chi}}\newcommand{\bipsi}{\boldsymbol{\psi}}\newcommand{\biomega}{\boldsymbol{\omega}}\({\rm{Z}}_4^0\)\end{document}, and \begin{document}\newcommand{\bialpha}{\boldsymbol{\alpha}}\newcommand{\bibeta}{\boldsymbol{\beta}}\newcommand{\bigamma}{\boldsymbol{\gamma}}\newcommand{\bidelta}{\boldsymbol{\delta}}\newcommand{\bivarepsilon}{\boldsymbol{\varepsilon}}\newcommand{\bizeta}{\boldsymbol{\zeta}}\newcommand{\bieta}{\boldsymbol{\eta}}\newcommand{\bitheta}{\boldsymbol{\theta}}\newcommand{\biiota}{\boldsymbol{\iota}}\newcommand{\bikappa}{\boldsymbol{\kappa}}\newcommand{\bilambda}{\boldsymbol{\lambda}}\newcommand{\bimu}{\boldsymbol{\mu}}\newcommand{\binu}{\boldsymbol{\nu}}\newcommand{\bixi}{\boldsymbol{\xi}}\newcommand{\biomicron}{\boldsymbol{\micron}}\newcommand{\bipi}{\boldsymbol{\pi}}\newcommand{\birho}{\boldsymbol{\rho}}\newcommand{\bisigma}{\boldsymbol{\sigma}}\newcommand{\bitau}{\boldsymbol{\tau}}\newcommand{\biupsilon}{\boldsymbol{\upsilon}}\newcommand{\biphi}{\boldsymbol{\phi}}\newcommand{\bichi}{\boldsymbol{\chi}}\newcommand{\bipsi}{\boldsymbol{\psi}}\newcommand{\biomega}{\boldsymbol{\omega}}\({\rm{Z}}_6^0\)\end{document} shifted to negative values during 4-D accommodation. The change in \begin{document}\newcommand{\bialpha}{\boldsymbol{\alpha}}\newcommand{\bibeta}{\boldsymbol{\beta}}\newcommand{\bigamma}{\boldsymbol{\gamma}}\newcommand{\bidelta}{\boldsymbol{\delta}}\newcommand{\bivarepsilon}{\boldsymbol{\varepsilon}}\newcommand{\bizeta}{\boldsymbol{\zeta}}\newcommand{\bieta}{\boldsymbol{\eta}}\newcommand{\bitheta}{\boldsymbol{\theta}}\newcommand{\biiota}{\boldsymbol{\iota}}\newcommand{\bikappa}{\boldsymbol{\kappa}}\newcommand{\bilambda}{\boldsymbol{\lambda}}\newcommand{\bimu}{\boldsymbol{\mu}}\newcommand{\binu}{\boldsymbol{\nu}}\newcommand{\bixi}{\boldsymbol{\xi}}\newcommand{\biomicron}{\boldsymbol{\micron}}\newcommand{\bipi}{\boldsymbol{\pi}}\newcommand{\birho}{\boldsymbol{\rho}}\newcommand{\bisigma}{\boldsymbol{\sigma}}\newcommand{\bitau}{\boldsymbol{\tau}}\newcommand{\biupsilon}{\boldsymbol{\upsilon}}\newcommand{\biphi}{\boldsymbol{\phi}}\newcommand{\bichi}{\boldsymbol{\chi}}\newcommand{\bipsi}{\boldsymbol{\psi}}\newcommand{\biomega}{\boldsymbol{\omega}}\({\rm{Z}}_4^0\)\end{document} negatively correlated with those in CMT1, and the negative change in \begin{document}\newcommand{\bialpha}{\boldsymbol{\alpha}}\newcommand{\bibeta}{\boldsymbol{\beta}}\newcommand{\bigamma}{\boldsymbol{\gamma}}\newcommand{\bidelta}{\boldsymbol{\delta}}\newcommand{\bivarepsilon}{\boldsymbol{\varepsilon}}\newcommand{\bizeta}{\boldsymbol{\zeta}}\newcommand{\bieta}{\boldsymbol{\eta}}\newcommand{\bitheta}{\boldsymbol{\theta}}\newcommand{\biiota}{\boldsymbol{\iota}}\newcommand{\bikappa}{\boldsymbol{\kappa}}\newcommand{\bilambda}{\boldsymbol{\lambda}}\newcommand{\bimu}{\boldsymbol{\mu}}\newcommand{\binu}{\boldsymbol{\nu}}\newcommand{\bixi}{\boldsymbol{\xi}}\newcommand{\biomicron}{\boldsymbol{\micron}}\newcommand{\bipi}{\boldsymbol{\pi}}\newcommand{\birho}{\boldsymbol{\rho}}\newcommand{\bisigma}{\boldsymbol{\sigma}}\newcommand{\bitau}{\boldsymbol{\tau}}\newcommand{\biupsilon}{\boldsymbol{\upsilon}}\newcommand{\biphi}{\boldsymbol{\phi}}\newcommand{\bichi}{\boldsymbol{\chi}}\newcommand{\bipsi}{\boldsymbol{\psi}}\newcommand{\biomega}{\boldsymbol{\omega}}\({\rm{Z}}_3^1\)\end{document} correlated with changes in RAL and CMT1.

Conclusions

HOA components altered during step-controlled accommodative stimuli. Ciliary muscle first contracted during stepwise accommodation, which may directly contribute to the reduction of spherical aberration (SA). The lens morphology was then altered, and the change in anterior lens surface curvature was related to the variation of coma.

Image registration reveals central lens thickness minimally increases during accommodation

Purpose

To evaluate anterior chamber depth, central crystalline lens thickness and lens curvature during accommodation.

Setting

California Retina Associates, El Centro, CA, USA.

Design

Healthy volunteer, prospective, clinical research swept-source optical coherence biometric image registration study of accommodation.

Methods

Ten subjects (4 females and 6 males) with an average age of 22.5 years (range: 20–26 years) participated in the study. A 45° beam splitter attached to a Zeiss IOLMaster 700 (Carl Zeiss Meditec Inc., Jena, Germany) biometer enabled simultaneous imaging of the cornea, anterior chamber, entire central crystalline lens and fovea in the dilated right eyes of subjects before, and during focus on a target 11 cm from the cornea. Images with superimposable foveal images, obtained before and during accommodation, that met all of the predetermined alignment criteria were selected for comparison. This registration requirement assured that changes in anterior chamber depth and central lens thickness could be accurately and reliably measured. The lens radii of curvatures were measured with a pixel stick circle.

Results

Images from only 3 of 10 subjects met the predetermined criteria for registration. Mean anterior chamber depth decreased, −67 μm (range: −0.40 to −110 μm), and mean central lens thickness increased, 117 μm (range: 100–130 μm). The lens surfaces steepened, anterior greater than posterior, while the lens, itself, did not move or shift its position as appeared from the lack of movement of the lens nucleus, during 7.8 diopters of accommodation, (range: 6.6–9.7 diopters).

Conclusion

Image registration, with stable invariant references for image correspondence, reveals that during accommodation a large increase in lens surface curvatures is associated with only a small increase in central lens thickness and no change in lens position.

Longitudinal Analysis of Mouse SDOCT Volumes

Spectral-domain optical coherence tomography (SDOCT), in addition to its routine clinical use in the diagnosis of ocular diseases, has begun to find increasing use in animal studies. Animal models are frequently used to study disease mechanisms as well as to test drug efficacy. In particular, SDOCT provides the ability to study animals longitudinally and non-invasively over long periods of time. However, the lack of anatomical landmarks makes the longitudinal scan acquisition prone to inconsistencies in orientation. Here, we propose a method for the automated registration of mouse SDOCT volumes. The method begins by accurately segmenting the blood vessels and the optic nerve head region in the scans using a pixel classification approach. The segmented vessel maps from follow-up scans were registered using an iterative closest point (ICP) algorithm to the baseline scan to allow for the accurate longitudinal tracking of thickness changes. Eighteen SDOCT volumes from a light damage model study were used to train a random forest utilized in the pixel classification step. The area under the curve (AUC) in a leave-one-out study for the retinal blood vessels and the optic nerve head (ONH) was found to be 0.93 and 0.98, respectively. The complete proposed framework, the retinal vasculature segmentation and the ICP registration, was applied to a secondary set of scans obtained from a light damage model. A qualitative assessment of the registration showed no registration failures.