Abstract
Purpose The development of neuro-ophthalmic symptoms during long-duration spaceflight presents major risks to astronauts and their missions. Amid efforts to understand spaceflight associated neuro-ocular syndrome (SANS), uncovering the role of the choroid in its etiology is challenged by the accuracy of image segmentation. The present study extends deep learning-based choroid quantification from optical coherence to-mography (OCT) to the characterization of pulsatile and topological changes in the macular plane and investigates changes in response to prolonged microgravity exposure.
Methods We analyzed OCT macular videos and volumes acquired from astronauts before, during and after long-duration spaceflight. Deep learning models were fine-tuned for choroid segmentation and combined with further image processing towards vascularity quantification. Statistical analysis was performed to determine changes in time-dependent and spatially averaged variables from preflight baseline for OCT data.
Results For 13 astronauts with a mean age of 47 ±9 years, there were significant increases in choroid thickness (CT) and luminal area (LA) of 14 ±25 µm (P = 0.01) and 0.11 ±0.2 mm2 (P = 0.009) and no significant change in choroid vascularity index (CVI) averaged over OCT macular video segments. There was no significant change in the amplitude of pulsatile CT fluctuation (ΔCT), but significant increases in both pulsatile LA (ΔLA) and CVI (ΔCVI) of 0.077 ±0.14 mm2 (P = 0.009)±and 0.03 0.06 (P = 0.01), respectively. For a subgroup of 6 astronauts for which inflight imaging was available, choroid volume, luminal volume and CVI over the macular region all increased significantly during spaceflight with the greatest differences reached at 90 days following launch.±
Conclusions The findings show that localized choroid pulsatile changes occur following prolonged microgravity exposure. They suggest that the choroid vessels expand in a manner similar to the choroid layer across the macular region during spaceflight, accompanied by a relative increase in the space they occupy. The methods developed provide new tools and avenues for studying and establishing effective countermeasures to risks associated with long-duration spaceflight.
Introduction
Following prolonged exposure to microgravity, astronauts present with distinct neuro-ophthalmic findings collectively referred to as spaceflight associated neuro-ocular syndrome (SANS)21; 26. The symptoms defining SANS comprise unilateral and bilateral optic disc edema, posterior globe flattening, choroidal and retinal folds, hyperopic refractive error shifts, retina nerve fiber layer thickening and infarcts and increased cerebrospinal fluid volume in optic nerve sheaths12; 13; 28; 14; 21; 26. Depending on its severity, optic disc edema may result in permanent vision loss, and it has been identified has the chief risk associated with SANS28. Choroidal tissue undergoes acute variations in space over time25, and finite element modelling of the posterior layers of the eye suggests that increased pulsatile choroid volume fluctuation may elevate strains in prelaminar neural tissue. The increased strains could in turn lead to the development of edema4; 26.
Pulsatile choroid volume fluctuation has been used to compute ocular rigidity and evaluate its response to long-duration spaceflight34. The effect of microgravity exposure on pulsatile choroid volume fluctuation itself has not yet been scrutinized. However, there is evidence that the choroid expands as soon as it is exposed to microgravity33; 17; 26, and the peripapillary choroid remains thicker than the pre-spaceflight baseline for at least 30 days and up to 90 days after landing17; 26. The view that choroidal swelling may be in part responsible for the hyperopic shifts and choroidal folds observed in SANS19; 12; 13; 14, and the possible link between altered choroid pulsatility and edema warrant seeking a refined characterization of choroidal changes in response to microgravity.
The development of optical coherence tomography (OCT) with enhanced depth imaging has made possible the visualization of details of the choroid in most imaged subjects25 and enabled the performance of reproducible biometric measurements39; 25; 11. Several means of quantifying the morphology of the choroid from OCT representations have been devised31. Albeit routinely resorted to in clinical research contexts, manual segmentation of the choroid is labor-intensive35; 41; 42 and ill-suited to the considerable volume of data produced by OCT devices, entailing a pace vastly exceeded by recent automated approaches41. The automated choroid semantic segmentation landscape has evolved from the use of various classical image processing methods2; 22 to the application of deep neural networks42. Under the supervised deep learning paradigm, the transformer architecture has demon-strated superior capabilities compared to convolutional neural networks (CNNs) when applied to computer vision tasks including segmentation27; 36.
Within the context of choroid segmentation, a few methods enable the differentiation of vasculature lumen and walls. It is now possible to identify, within a single OCT B-scan, the luminal area (LA) and to calculate choroidal vascularity index (CVI), defined as the fraction of LA over the total choroid area. However, the successful application of such parameters is highly dependent upon the accurate delineation of the choroid1.
In the present work, we train deep learning models to reliably produce choroid semantic segmentations in macular images from OCT videos and OCT volumes acquired from astronauts. We complement our segmentation techniques with vascularity quantification, enhancing the ability to characterize the state of the choroid and its variations. We extend the capability of our approaches to OCT image timeseries and OCT volume reconstructions. For video acquisitions, we compare pre- and post-spaceflight quantifications and reveal significant changes in time-dependant variables. For volumes, we capture substantial mapping changes over several timepoints in a cohort of astronauts exposed to microgravity.
Methods
Datasets
The study was approved by the institutional review boards of the Maisonneuve-Rosemont hospital and NASA, and the clinical research ethics committee of Université de Montréal. The study was conducted according to the Declaration of Helsinki and its amendment. Two datasets were used (see table S1 for a summary). All acquisitions were performed using spectral-domain OCT (Spectralis OCT2; Heidelberg Engineering, Heidelberg, Germany). There was no pathology reported in any of the subjects participating in the study.
OCT macular videos were acquired pre- and post-spaceflight from both eyes, with up to 5 videos per eye at each timepoint. Each frame has dimensions 768 × 496 pixels, approximately corresponding to 8.4 × 1.9 mm, and each video is roughly 1 minute in length. A resting heart rate measurement was also obtained con-currently with video acquisition.
For the analysis of changes in 3-dimensional topology due to microgravity exposure, we used a set of OCT macular volumes composed of 97 or 193 B-scans. The macula block setting was used with placement over the fovea via the anatomic positioning system for each acquisition. For each timepoint, up to 2 volumes per eye were obtained. All frames in the set of OCT macular volumes measure 512 × 496 pixels, which physically approximates to 5.6 × 1.9 mm.
Segmentation performance evaluation
SegFormer models40 pretrained on the ImageNet database were retrained directly on both datasets to-wards choroid segmentation. For OCT videos, training data was gathered by randomly selecting 10 B-scans per movie where the choroid region was manually labelled by one trained grader and binarized. To implement 5-folds cross-validation, the training data was arranged 5 different ways, in each instance splitting the whole into training, validation and test sets. For each fold, 10 % of B-scans were assigned to the validation set, and 10 % to the test set. The training and validation sets were used for model training and the test sets were used for performance evaluation.
For OCT volumes, we devised a strategy to ensure adequate representation of the distribution of B-scans corresponding to each subject’s eye in the training data. The image set was split into 5 groups along the slow acquisition direction, and 4 B-scans from each group were randomly selected for a total of 20 per volume. For each volume, 4 B-scans out of 20 were manually labelled by one trained grader and binarized. The locations of the labelled B-scans were subsequently shifted as shown in figure S7. Training data 5-folds cross-validation was then implemented in the same manner as described for videos.
Post-processing was implemented to ensure continuous layer boundaries and minimize the effect of isolated false positive regions (details are provided in Supplementary Information B). The test set 5-fold cross-validation performance for both the videos and volumes datasets was evaluated based on similarity metrics, as well as choroid boundaries mean absolute error (MAE) and choroid thickness difference (TD) between the manually segmented and corresponding automatically segmented frames. The equations for MAE and TD are provided below where m is the B-scan width and y(i) and ŷ(i) are the inferred and manually labeled vertical boundaries of Bruch’s membrane (BM) or the choroid-sclera interface (CSI). Sørensen–Dice coefficient complement (DiceC) and Jaccard coefficient complement (JaccardC) were used as similarity metrics. We also defined mathematical expressions to identify inaccurate segmentation and segmentation failure associated with subtle and gross segmented boundary aberrations, respectively (equations are provided in Supplementary Information C).
Beyond choroid segmentation, we produced masks of LA to calculate CVI (details are provided in Supplementary Information D). CVI was calculated as LA divided by total choroid layer area. The LA processing performance was assessed on the test set of selected Seg-Former models for both videos and volumes using similar methods as for choroid segmentation performance evaluation. Additional test set performance evaluations relating to choroid layer and lumen segmentation are provided in Supplementary Information H.
Temporal analysis
Inference was only performed on B-scans which met noise specifications based on OCT device manufacturer defined metrics. Specifically, only B-scans with a quality score equal to or greater than 24 (0-40 score range; 0 indicating no signal and 40 indicating excellent quality) and with an automatic real time frame averaging number greater or equal to 2 were considered. The automatic real time averaging setting was 3 frames per B-scan.
B-scans with a signal-to-noise ratio beneath 0.55 were also excluded. The signal-to-noise ratio equation used is provided below where SNR is the signal-to-noise ratio, S is the mean and σ is the standard deviation of the intensities in a B-scan. B-scans for which the inferred choroid region was located within 10 % of the superior or inferior vertical edge of the image were excluded. Single B-scans with an inferred choroid area (CA) departing from the average of the manually segmented frames by greater or equal to 20 % were also excluded.
OCT videos were unevenly sampled and frequency analysis was performed with Lomb–Scargle periodograms and spectrograms generated over a 0.1-4 Hz frequency range (transforms implementation details and calibration test results are provided in Supplementary Information F). Movies with less than 50 good quality frames were discarded. To prevent large gaps in a video, we preserved only the longest video segment for which no zone of excluded frames amounting to greater than 1 s was present. Pulsatile CT, LA and CVI (ΔCT, ΔLA, ΔCVI) were computed as done by Beaton et al2. First, the locations of peaks and valleys across time series were obtained, where the minimum distance between them is given by the expression where dmin is the minimum distance between peaks and between valleys, Theart is the heart period and Mo, intervals is the mode of the 1D array of intervals between each element in a time series. Having obtained peaks and valleys locations, consecutive peaks between valleys and consecutive valleys between peaks were discarded. ΔCT, ΔLA and ΔCVI were then defined as the difference between the median of all remaining peaks and the median of all remaining valleys.
Linear mixed models were used to determine statistical significance of choroid quantification differences pre- and post-spaceflight. For each linear mixed model, timepoints were designated as the fixed effect. Subjects and eyes nested under each subject were specified as random effects3. Significance values were gathered using the Kenward-Roger approximation16.
Spatial analysis
For OCT volumes, the automatic real time averaging setting was 15 frames per B-scan. Following initial segmentation, smoothing of the surface of the CSI across each macular volume was implemented using a graph search algorithm similar to the one described by Mazzaferri et al22. The graph search algorithm was applied to each one of the volume’s cross-sectional image along the fast acquisition direction.
Having delineated the region between the BM and a smoothed CSI, choroid volume (CV), luminal volume (LV) and CVI were quantified (CT and LT are included in Supplementary Information K). In addition to a global level value, the result of each quantification was split into macular ETDRS subfields (center; superior inner; nasal inner; inferior inner; temporal inner; superior outer; nasal outer; inferior outer; temporal outer).
For better visualization, we created square maps of CT, LT and CVI of 512 × 512 pixels by interpolating the total number of B-scans of every OCT volume and we analysed variation across timepoints. Linear mixed models were used to determine statistical significance between timepoints. The same specifications were used to implement the linear mixed models as for the temporal choroid OCT data analysis.
Loss function
The training of the SegFormer models was achieved through the IoU loss. For binary segmentation, which was used exclusively in the present study, the IoU loss equation is given by where ℒIoU is the IoU loss, N is the number of pixels in an image, si is the predicted value of the ith pixel in an image, and gi is the targeted ground truth value of the ith pixel in an image.
Training process
The SegFormer models were trained using the Adam optimizer. A 0.0001 learning rate was used as well as early stopping based on validation loss for both Seg-Former models. No data augmentation strategy was carried as part of training the SegFormer models.
All training was implemented on a Windows 11 computer equipped with a NVIDIA GeForce RTX 2060 Graphics Processing Unit (GPU). Training of the Seg-Former transformer models was completed on Python (v3.10.0) using PyTorch Lightning (1.8.6).
Results
SegFormer-based choroid segmentation
We sought to quantify choroid morphology changes in response to microgravity exposure. For this, we retrained SegFormer models directly on astronaut OCT macular videos and volumes and implemented a simple post-processing strategy to analyze vascular changes. We manually labeled a subset of images representing 2 % of available images and assessed the test set performance of the pipeline using different metrics.
We leveraged a dataset of OCT macular videos to compare pre- and post-spaceflight time-dependant variables obtained from astronauts. We used images centered at the macula from 13 individuals, 9 male and 4 female of average age 47 ±9 yo. There were 8 astronauts who did not have prior spaceflight experience, 3 who had completed Space Shuttle short-duration missions, and 4 who had previous long-duration spaceflight experience (> 4 months15). All preflight videos were obtained within a 9 to 1 month before spaceflight window and all postflight videos within 2 to 30 days after spaceflight.±Figure 1a and b display illustrations of 3 individuals’ manual and automated segmentations showing excellent correspondence for both choroid and vascular lumen obtained from OCT videos. Mean DiceC and mean JaccardC for the OCT videos test set (figure 1c) were 0.044 ±0.02 and 0.083 ±0.03, respectively. Figure 1d provides segmentation mean absolute error (MAE), expressed as the deviation of the boundary (in pixel units) from the labeled traces. It also shows thickness difference (TD) which represents the deviation of the layer thickness (in pixels) from the labeled layer. Test set BM and CSI MAE were 0.93 ± 0.3 pixels and 6.2 ± 2.7 pixels, respectively, and TD was 6.1 ± 2.7 pixels. Test set LA mean DiceC and mean JaccardC were 0.32 ± 0.1 and 0.48 ± 0.09, respectively (figure 1e).
The performance was similar for the B-Scans from OCT volumes (manual and automated segmentations are shown in figure 1f and g). This second set contains OCT volumes centered at the macula from 6 of the 13 astronauts featured in on-Earth videos, 5 males and 1 female, which were of average age 48 ± 9 yo. There were 4 astronauts who did not have prior space-flight experience, 1 who had completed Space Shuttle short-duration missions, and 2 who had previous long-duration spaceflight experience. The volumes were obtained at up to 6 different timepoints before launch (launch - 21 to 18 months, launch - 9 to 6 months), during flight and before return (launch + 30 days, launch + 90 days, return - 30 days) and following spaceflight (return + 1 to 3 days). Test set choroid segmentation mean DiceC and mean JaccardC for OCT macular volumes were 0.057 ± 0.02 and 0.11 ± 0.04, respectively (figure 1h). Test set BM and CSI MAE for the same fold were 1.2 ± 0.3 pixels and 8 ± 4.2 pixels, respectively, and TD was 7.9 ± 3.9 pixels (figure 1i). Test set LA mean DiceC and mean JaccardC were 0.39 ± 0.07 and 0.55 ± 0.07, respectively (figure 1j). Additional choroid layer segmentation and lumen segmentation performance evaluations are provided in supplementary figures S8 and S9.
Quantification of choroid pulsatile changes
We investigated the presence of spaceflight-related changes in choroid pulsatility. For this, as we have done in the past2, we obtained videos of OCT images of about 1 minute, acquired at 14 Hz. The builtin eye tracker of the OCT device was used to assure that the same location in the macula was sampled throughout the video. The acquisition was automatically halted when a patient moved, only resuming after the laser was repositioned. Our choroid quantification approaches were implemented on the videos, yielding timeseries which displayed morphological changes at the frequency of the heart. Because images are not sampled at a constant rate, Lomb–Scargle peri-odograms and spectrograms of normalized power were expected to show peaks in the heart rate vicinity. Detection of those peaks, as illustrated in figure 2g-i, serves as quality control for segmentation. Narrow and prominent maxima can often be seen for CT, LA and CVI (also figures S11 and S12, second and third columns).
We found a global CT averaged over movie durations significant increase of 14 ±25 µm or 4.7 % (linear mixed model, P = 0.01) from pre-to postflight, and global LA significant increase of 0.11 ±0.2 mm2 or 6.8 % (linear mixed model, P = 0.009) (figure 2a and b), but no change in global CVI (figure 2c). For the analysis of the amplitude of pulsatile fluctuations, no significant differences in ΔCT were found globally (figure 2d). As shown in figure 2e and f, both global ΔLA and ΔCVI did show significant increases from pre-to postflight of 0.077 ±0.14 mm2 or 34 % (linear mixed model, P = 0.009) and 0.03 ±0.06 or 35 % (linear mixed model, P = 0.01), respectively. We provide baseline pre-spaceflight results of our choroid quantification methods for OCT videos from astronaut eyes, including their ETDRS subfield equivalents in supplementary table S2. The ET-DRS subfield equivalents of the measures displayed in figure 2a-f are provided in supplementary figure S10, and significance in table S4. Positive correlations were found between all global metrics and TRT250 but none were significant for the available data.
Spatial choroid changes associated with long-duration spaceflight
We extended our choroid quantification techniques to OCT volumes, for which large variations were expected in astronauts. Each volume is made up of 97 or 193 B-scans. In addition to obtaining cross-sectional image segmentation performance, we generated maps exposing choroid thickness topography and vascular patterns. Choroid vessel patterns were readily observable on LT and CVI maps in OCT volumes (figure 3d). For all eyes, there was a correspondence between CT and LT maps, with thicker choroid regions generally translating to thicker luminal regions. Different degrees of retinal vessel shadow artifacts were discernable in most scans.
When comparing individual timepoints with preflight, there were significant differences for global CV and LV (figure 3a and b) for all measurements. The same comparison for global CVI measurements indicated significant differences for inflight timepoints and normal values for return + 1 to 3 days (figure 3c). We provide baseline results of our choroid quantification methods for OCT volumes from astronauts, including their ET-DRS subfield equivalents in supplementary table S3.
The ETDRS subfield equivalents of measures displayed in figure 3a-c are provided in Supplementary Information K, and significance values in supplementary tables S5-S8.
Discussion
In addition to an increase in CT within the macular region during the month that follows long-duration spaceflight, we found an increase in the area corresponding to the lumen. Our results also indicate that while pulsatile fluctuation in CT did not change compared to preflight when measured within 30 days of return, the amplitude of pulsatile fluctuation in LA increases. During missions, we observed a global increase in the space occupied by lumen relative to the whole choroid layer across the macular plane.
Consistent with our results, a study which measured macular CT preflight and during spaceflight showed a 35 µm CT increase6. Previous reports demonstrated similar increases in peripapillary CT. Laurie et al found a mean peripapillary CT increase of 27 µm during space-flight compared to preflight10. Macias et al observed progressive peripapillary choroid thickening over the course of long-duration spaceflight, with a mean increase of 43 µm at 150 days17. In line with our findings for CT at the macula, they also showed still significantly increased peripapillary CT at 30 days following return to Earth. Unlike the present work, previous studies were limited to measuring CT from a single image. Here we complement previous reports by adding maps of thickness change, as well as an analysis of how lumen and CVI evolve and their response to cardiac pulsation.
Several mechanisms have been proposed to account for increased CT upon microgravity exposure. They include cephalad fluid shift resulting in venous congestion and accumulation of blood in the choroid18; 12; 14 as well as persistent buildup of choroidal interstitial fluid in the peripapillary area38; 26. Determining the source of optic disc edema observed in SANS has been described as challenging using available technology and approaches to data analysis17. Our results are in line with finite element modelling which suggested a link between choroid anatomy, increased pulsatile choroid volume fluctuation, prelaminar neural tissue strains and the development of edema4; 26.
A mathematical model of the cardiovascular response to long-duration spaceflight predicted increased upper body blood volume5. Multiple authors have demon-strated that the choroid expands in weightlessness, leaving open the question of choroidal pulsatility and its possible change in weightlessness. Our study directly shows increased amplitude of pulsatile fluctuation of the space occupied by choroid vessels after spaceflight, opening the possibility of investigating its relationship with peripapillary retinal thickness increase in SANS. Pulsatile deformation of the optic nerve head may be assessed using OCT32 and could also be studied in relation to choroid pulsatile changes, helping substantiate finite element modelling predictions.
Studies have linked the development of edema in SANS with single nucleotide polymorphisms (SNPs)43; 44. A lower enzyme functional activity in the 1-carbon path-way can increase B-vitamins requirements and the risk of vitamin insufficiency45. Folic acid and its active metabolite, 5-methyltetrahydrofolate, maintain endothelial nitric oxide synthesis and NADPH consumption coupling. Uncoupling decreases nitric oxide synthesis and increases the generation of reactive oxygen species8; 23. Oxidative stress in turn increases endothelial permeability, widening interendothelial junctions and resulting in protein and fluid leakage into the interstitial space7; 24. The present study did not differentiate between astronauts with specific genetic profiles, but the approaches we developed could be used to determine whether individuals with SANS-relevant SNPs are more likely to show choroid layer or lumen pulsatile fluctuation changes, and whether they are more vulnerable to them.
The application of lower body negative pressure has been explored as a means of mitigating the effects of cephalad fluid shift in SANS. While lower body negative pressure was shown to decrease intraocular pressure in astronauts during spaceflight, no change in macular CT was observed6. Numerical modelling suggests that under microgravity conditions, less rigid vessels are subjected to increased transmural pressure and greater than normal volume9. An artificial gravity experiment achieved partial alleviation of mice retinal tissue damage associated with weightlessness exposure20; 37. The effects of artificial gravity on the state of the choroid have not been explored, and the need for precise and reliable methods of assessing SANS-related eye alterations to validate its effectiveness as a countermeasure has been emphasized37. The approaches we developed for the measurement of pulsatile and volume choroid changes represent a step in this direction.
The significant postflight global and localized increases in ΔLA and ΔCVI we observed could be associated with increased choroidal vessel compliance following sustained greater than normal volume during long-duration spaceflight. Pulsatile choroid volume fluctuation correlates inversely with the measurement of ocular rigidity performed using an invasive procedure30. Combined with additional measurements, our choroid layer segmentation method has the potential to improve upon a key element in the reliable, non-invasive computation of ocular rigidity2; 29. An early exploration of ocular rigidity in a private astronaut suggested that it decreases following spaceflight33, and this was recently confirmed in professional astronauts34. Our work provides the means to build upon previous findings and further unravel the relationship between ocular rigidity and spaceflight, including its role a risk factor for developing of SANS.
Conclusions
Neuro-ophthalmic findings in spaceflight associated neuro-ocular syndrome (SANS) including optic disc edema constitute important risks for astronauts and their missions. Weightlessness results in rapid choroidal expansion which may partially account for concerning SANS findings. Macular region luminal area and the amplitude of its pulsatile fluctuation were both increased compared to preflight within 30 days after long-duration spaceflight. During long-duration space-flight, there is a global increase in both the choroid layer and the relative space occupied by lumen inside the choroid across the macular plane.
Data Availability
All data produced in the present study are available upon reasonable request to the authors.
Acknowledgements
Funding was provided by the Canadian Space Agency, the Canadian Institutes of Health Research, and the Fonds de Recherche en Ophtalmologie de l’Université de Montréal (FROUM). MMS received a scholarship from the Fonds de Recherche du Québec Santé. SC holds the Wolfe Professorship in Translational Research. We thank NASA and ESA for coordination and infrastructure.